[Alexmine]

Thank you for your response.

[Alexmine]

Dear Adrian,

Thank you for your reply. I would like to know the CUDA-ARCH value of RTX4070, and whether RTX4070 can be used for parallel computing in OpenLB.

Sincerely,
Alexmine.

[Alexmine]

Dear Adrian,

The characteristic length ( the radius of sphere) D= 0.02m, resolution N=20, so the deltaX=0.001m,the size of Computational domain is 15D*8D*8D. The basic unit of my STL is mm,so I set stlsize=0.001.I think there is no problem .
Then I checked indicatorF3D.hh,from line 621 to line 624:
621IndicatorCuboid3D<S>::IndicatorCuboid3D(Vector<S,3> extend, Vector<S,3> origin)
622 : _center(origin+.5*extend),_xLength(extend[0]), _yLength(extend[1]), _zLength(extend[2])
623{
624 assert(_xLength>0 && _yLength>0 && _zLength>0);
As _xLength= extend[0], _yLength= extend[1], _zLength= extend[2],the assertion error may caused by the component coordinates of the vector extend, at least one of them is zero or negative value.
So I replaced the original cylinder3d.stl file in the example cylinder3d with the exact same model created by the modeling software CROE, and the same assertion error occurred.
It looks like that the source of the error seems to be caused by offset of coordinate systems between OpenLB and CROE, so what should I do?

Sincerely,
Alex

[Alexmine]

Dear Adrian,

Thankyou for your reply.

Sincerely,
Alex

[Alexmine]

Dear Mathias,
Thank you for your reply. In fact, I have written my own program through the case of laminar flow, thanks to the structure of OpenLB is clear and easy to understand. What bothers me now is how to achieve a stable layered fluid flow around the obstacles in the channel. I don’t know how to correctly implement the layered fluid in OpenLB. Given that density is a very sensitive physical quantity for OpenLB, obviously I should not set layered fluid density directly, so what should I do?

Sincerely,
Alex

[Alexmine]

Thanks, FBukreev. This is indeed very helpful.

[Alexmine]

Thank you very much for your reply!
I have another question in the case laminar cylinder3d,the code:
maxVelocity[0] = 2.25*frac[0]*converter.getCharLatticeVelocity();
what is the specific value of frac[0] ?

[Author]

Posts