Skip to content

Natural convection in square cavity : trouble with the discretization

OpenLB – Open Source Lattice Boltzmann Code Forums on OpenLB General Topics Natural convection in square cavity : trouble with the discretization

Viewing 2 posts - 1 through 2 (of 2 total)
  • Author
  • #8763


    I’m currently a beginner in OPENLB and I started recently natural convection simulations using the “squareCavity2d” example provided with the last release of the code (/example/thermal/squareCavity2d directory). I set the value of the Rayleigh number to 10e3 and I tried to test the influence of the spatial discretization by changing the number N in the .cpp file. I varied the number N between N= 128 and N = 140 and I encountered some trouble with the allocation of the material number to the right vertical wall of the 2D square. The allocation is well done except for the values N = 131 – 132 – 135 and 139. In that cases, there is only one count associated with the right wall (material number = 3) whereas the count number is correct (= N) for the left wall (material number = 2). I don’t understand why the allocation is ok for certain values and not for others (I also tested N values between N = 5 and 10 and I have alse a trouble for the values N = 5 – 8 – 9 and 10). I precise that I run the calculations in parallel mode (mpi option) on 8 cores.
    Could anyone provide guidance or resources that could help me to address this issue?
    Thank you in advance for your help!

    Best regards


    Dear Mutelle,

    sorry for the late response.
    I looked into the implementation of the case and it seems that that example uses a deprecated way of assigning material numbers. I would refer to the laminar/cavity2d example regarding this issue. We will work on updating the deprecated case in the next release.
    If you need more specific assist, please let us know.

    Best wishes,
    Shota Ito

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