Viewing 2 posts - 1 through 2 (of 2 total)
• Author
Posts
• #5929
achodankar
Participant

Hello Developers,
I have a question in regard to computing the mass flux. I get the mass flux and area(line in 2D) values as zero(i.e the first and second index of output array). Also, I would like to know the significance of the last three values in the output array. In the handbook, it says that “flow vector (ie. vector of summed quantities)”. What would these values represent in velocity, pressure, and mass flux. Lastly, we specify the origin and direction vector for computing the flux, so it means we compute the flux along the line of the origin depending on the direction of the vector, Correct?

Here is the code snippet for computing the mass flux:

{
SuperLatticeVelocity2D<T, DESCRIPTOR> velocityF( sLattice);
SuperLatticeDensity2D<T,DESCRIPTOR> densityF(sLattice);

T conversationFactorMass = converter.getConversionFactorMass();
T conversationFactorTime = converter.getConversionFactorTime();

int Lx = converter.getLatticeLength( lx );
int Ly = converter.getLatticeLength( ly );

T point;
point = Lx/2;
point = Ly/2;

ofstream mflux1;
mflux1.open(“mflux_num_lat_”+ std::to_string(iT) + “.dat”);

SuperPlaneIntegralFluxMass2D<T> massFlux(velocityF, densityF, superGeometry, conversationFactorMass, conversationFactorTime, {Lx/2, Ly/2}, {0., 1.}, BlockDataReductionMode::Analytical );

AnalyticalFfromSuperF2D<T> intpolatemass_lat( massFlux, true );

T numerical_mflux_lat;

intpolatemass_lat( numerical_mflux_lat,point);

mflux1<<“————————————————————“<<endl;
mflux1<<“The value of x is “<<point<<endl;
mflux1<<“The value of y is “<<point<<endl;
mflux1<<“The mass flow rate value is “<<numerical_mflux_lat<<endl;
mflux1<<“The size of the area value is “<<numerical_mflux_lat<<endl;
mflux1<<“The flow vector (ie. vector of summed quantities) value is “<<numerical_mflux_lat<<endl;
mflux1<<“The flow vector (ie. vector of summed quantities) value is “<<numerical_mflux_lat<<endl;
mflux1<<“The flow vector (ie. vector of summed quantities) value is “<<numerical_mflux_lat<<endl;
mflux1<<“————————————————————“<<endl;

mflux1.close();

}

Thank you.

Yours sincerely,

Abhijeet C.

#5951
Specifying the origin and direction is one possibility of parametrizing the hyperplane used by the plane integral functors. You can see other possibilities that are overloaded by constructors of `SuperLatticeIntegralFluxMass2D` in e.g. Doxygen. So yes, the flux is computed through the plane that is orthogonal to the direction vector and intersects the origin.
As for the output vector: `SuperPlaneIntegralFluxMass2D` scales the integral obtained by `SuperPlaneIntegral2D` by the given factors. The other documented components are returned unchanged.