[mchurchf]

To whomever can answer,

I have two questions about buoyantBoussinesqPisoFoam.

1). Why is pEqn.setReference(pRefCell,pRefValue) not called after setting up the pressure equation? If you take a look at pisoFoam, the pressure referencing is done, but why isn't it done in buoyantBoussinesqPisoFoam?

2). Why is the temperature equation solved before the PISO loop. Currently, the T equation is solved right after the velocity predictor step. The predicted velocity is not divergence-free, and the T equation uses this incorrect velocity in its convection term. Wouldn't it be more appropriate to solve for T after the PISO loop with a good divergence-free velocity? I guess that would mean that velocity and pressure are solved with a time lagged T.

Thanks,

Matt