Taylor Green Vortex Velocity Initialization

[liu]

Dear Community,

I would like to reproduce the Taylor-Green vortex case. In this case, the velocity initialization is given by:
ux​(x,y)=−cos(2πx)sin(2πy)
uy(x,y)=sin(2πx)cos(2πy)
Which case should I refer to, and how should I initialize the velocity field? thank you!

[stephan]

Dear liu,

thank you for your post.
The best example to start would be examples/turbulence/tgv3d/.
I suggest to change the 3D functors to 2D, and remove all the lines you do not require.
The intialization should also be simplified, such that only defineU and defineRho is left.
The velocity field for u_0 is constructed as functor and then passed to defineU.
The pressure for p_0 (given for example in the book of Krüger et al. (2017)) is constructed as a functor, then transformed to density and passed to defineRho.

If you want to have specific support on creating your own applications with OpenLB, please consider taking part in our upcoming Spring School: https://www.openlb.net/spring-school-2025/

BR
Stephan

[liu]

thank you

[Author]

Posts