Viewing 12 posts - 1 through 12 (of 12 total)
• Author
Posts
• #7210
H.Yu
Participant

Hello
I am trying to calculate aerodynamic performance using OpenLB.
I studied some examples. In the 3D cylindrical laminar flow example, I am trying to proceed with calculation by defining it as turbulent flow, but some problems occur.

In the example of a 3D laminar flow cylinder, the result values ​​such as pressure 1, pressure 2, pressure drop, drag force, and lift are not calculated as -nan.

I try to calculate by adjusting the physDeltaX and time step, but I can’t overcome -nan.

What methods can I try?

The code of the part defining the turbulence is also attached.

CODE : bulk dynamics = SmagorinskyBGKdynamics
// Material=1 –>bulk dynamics
sLattice.defineDynamics<BulkDynamics>(superGeometry, 1);

// Material=2 –>bounce back
sLattice.defineDynamics<BounceBack>(superGeometry, 2);

// Setting of the boundary conditions

//setLocalVelocityBoundary(sLattice, omega, superGeometry, 3);
//setLocalPressureBoundary(sLattice, omega, superGeometry, 4);

//inlet
sLattice.defineDynamics<BulkDynamics>(superGeometry, 3);
setInterpolatedVelocityBoundary<T,DESCRIPTOR>(sLattice, omega, superGeometry, 3);
//outlet
sLattice.defineDynamics<BulkDynamics>(superGeometry.getMaterialIndicator({4}));
setInterpolatedPressureBoundary<T,DESCRIPTOR>(sLattice, omega, superGeometry.getMaterialIndicator({4}));
//setInterpolatedPressureBoundary(sLattice, omega, superGeometry, 4);

// Material=5 –>bouzidi / bounce back
//airfoil
#ifdef BOUZIDI
sLattice.defineDynamics<NoDynamics>(superGeometry, 5);
#else
sLattice.defineDynamics<BounceBack>(superGeometry, 5);
#endif

help!

#7211
FBukreev
Keymaster

Hello,

I think you have not defined the Smagorinsky Constant there
`sLattice.setParameter<collision::LES::Smagorinsky>(T(0.15));`

Greetings
Fedor

#7213
H.Yu
Participant

None of the above, but I already defined it as the Smagorinsky constant 0.1

#7214
H.Yu
Participant

Is the meaning of ‘-nan’ uncalculated or diverged?

#7215
H.Yu
Participant

What is the relationship between voxel size and stl size if I use stlreader? voxel size =<stl size ??

#7216
H.Yu
Participant

I attach a photo.
https://ifh.cc/g/n5rT5x.jpg

This picture is a speed contour. What are the reasons why it looks like that picture?

I am testing various things such as voxel size, stl size, maxVel, resolution, and time discretization refinement, but it is still calculated at first, and -nan occurs as time passes. There is also the problem that yPlusMax keeps increasing.

Anything else to consider?

#7217
FBukreev
Keymaster

Hello,

you can look how the stl part is implemented in particles/bifurcation/eulerEuler example with IndicatorLayer function. It must ease the stl inclusion for you.
You can use the UnitConverterFromResolutionAndLatticeVelocity and set the lattice velocity to 0,1 there or with some smaller value. It can make the simulation more stable.

Greetings
Fedor

#7220
H.Yu
Participant

thank you!
I have one more question. Are there any usage examples to learn how to use IndicatorLayer?

#7221
FBukreev
Keymaster

You can also look into examples/laminar/cylinder3d.
This indicator just enlarges the stl geometry to lattice points depending on cell size.

#7222
H.Yu
Participant

thank you.
I set it up like the example above. Also, when I set ‘UnitConverterFromResolutionAndLatticeVelocity’, ‘-nan’ does not occur anymore.

#7223
H.Yu
Participant

I have one more question.
I want to plot various things using paraview for post-processing, but even if I load the calculated ‘.vtm’ file into paraview, only cuboid is visible. What can I think about to solve this problem?

• This reply was modified 1 year, 6 months ago by H.Yu.
#7227
FBukreev
Keymaster

You need to use the pvd file to see the results.

Viewing 12 posts - 1 through 12 (of 12 total)
• You must be logged in to reply to this topic.