[sahilbhapkar30]

Hello,
Thank-you for your reply. I did not get you completely.
There is only wall shear stress functor implemented. I dont have an idea about how to implement my own functor. Could you give some hints?

Thank-you

[sahilbhapkar30]

Hello,

Thank-you for your answer. I understood it completely.

[sahilbhapkar30]

Hello,

Thank-you for your prompt reply. I did not understand the meaning of the difference is a constant.

Thank-you

[sahilbhapkar30]

Hello,
Thank-you for your prompt response.

In the following forum discussion,
https://www.openlb.net/forum/topic/poiseuille2d-interpolated-and-local-valuesvelocity-and-pressure/

they have discussed about the density and pressure being related by the equation of state and the offset of 1 in the equation.

The formula for lattice pressure in the code is given as: p= cs2 * (rho-1). As per one of your colleagues, he mentioned that ” Instead, it simply seems to be since one wants to correspond density=1 to pressure=0 as a default. Consequently, one has to modify the offset when a different char. phys. pressure is defined. The UnitConverter does this.”

I did not get this completely, should I remove this offset of 1 in this equation as I am having a different pressure values?

Thank you for your earlier suggestion but I need to have velocity constraints on the input and output as well.

I would appreciate your help,

Thank-you.

[sahilbhapkar30]

Dear Stephan,

Thank you for your response. I have another doubt.I want to calculate the wall shear stress at the inlet and outlet of my 3D geometry.

I implemented this thing in the main:

SuperLatticePhysWallShearStress3D<T,DESCRIPTOR> wss(sLattice, superGeometry, 3, converter, inflow1);
SuperLatticePhysWallShearStress3D<T,DESCRIPTOR> wss1(sLattice, superGeometry, 3, converter, inflow2);
SuperLatticePhysWallShearStress3D<T,DESCRIPTOR> wss2(sLattice, superGeometry, 4, converter, outflow_wall);
SuperLatticePhysWallShearStress3D<T,DESCRIPTOR> wss3(sLattice, superGeometry, 4, converter, outflow_wall_2);

where inflow1, inflow2, outflow_wall and outflow_wall_2 are the indicator cuboids which I have used for the inlet and outflow.

In addition, I have implemented the functor (vtmWriter.addFunctor( wss)) for each of them in the get results.

But after the simulation, I am getting the wall shear stress value for the inlet flow only. If I comment the inflow region area, I do get the outflow wall shear stress value if I rerun the simulation.

Could you tell me, what I am doing wrong?

Thank you.

[Author]

Posts