Hi,
I’m working on BoussinesqThermalProcessor3D and have studied the BoussineqThermal3D case.
I have some doubts.

For the parameter “gravity” for BoussinesqThermalProcessor3D, the example gave RayleighBenardFlowParam.getLatticeGravity(). But RayleighBenardFlowParam.getLatticeGravity()= getDeltaT() * getDeltaT() / getDeltaX(), it seems that it is just the nondimension parameter for the gravity. So, for the real gravity, should I use 9.8 * RayleighBenardFlowParam.getLatticeGravity() instead?

For the parameter “dir” for BoussinesqThermalProcessor3D, in my opinion, it represents the direction vector for the gravity constant, for the real three dimensional should it be (0,0,1) or (0,0,1)? I saw the example employed (0,0,1). And when I tried (0,0,1), it did not work.

We know that Boussinesq approximation is F=  beta * g * (TT0), where beta is the coefficient of thermal expansion and usually beta=1/T0 [K]. But in the source code, I found F= gravity * dir / deltaTemp * (TT0), where deltaTemp=HighTempLowTemp. Why do we use deltaTemp=HighTempLowTemp here, while not beta=1/T0?

For AdvectionDiffusionBGKDynamics, I think its omega is obtained from kappa= nu/Pr. But when I’m using LES for the flow field, we also have kappa_sgs=nu_sgs/Pr_sgs. Is there any AdvectionDiffusion Dynamics that can read omega_total from (kappa + kappa_sgs), just like the omega_total in the flow field?
with best wishes,
steed188