[PanPeter]

Hello,

I am investigating the buoyantBoussinesqPimpleFoam solver and believe the gravity term is not calculated:

The pressure and gravity terms in the Navier-Stokes equations are redefined as:

,where


Now for the pimple loop:
1: Solving for the velocity field, where NSE with above pressure and gravity is solved.
2: The heat equation is solved, where is updated.
3: Updating the pressure field in the nested pressure-correction loop, where is updated.

Now for my problem:
I initialize the pressure field to be zero everywhere and disable the the pressure correction loop (nCorrectors=0 in fvSolution.PIMPLE). So the field is zero and will always be zero. That kills the pressure gradient term in the NSE. The dynamics in the systems should then be induced though the gravity term alone but no motion develops! Does anybody have an idea why?

I attached a lock-exchange case example, in which the pressure loop is deactivated and the Prandtl number is set to 10^9 to disable heat diffusion.

[MSF]

Hi,

Nothing unexpected happens. Check buoyantBoussinesqPisoFoam: A detailed explanation

Best Moritz