Skip to content

Reply To: Thermal Poiseuille 2D.

#6348
ramirofreile
Participant

I have also noticed that in AdvectionDiffusionTRTDynamics the collision function is:

cell[iPop] -= _omega2 * (fPlus[iPop] – fEqPlus[iPop]) + this->_omega * (fMinus[iPop] – fEqMinus[iPop]);

whereas in TRTDynamics:

cell[iPop] -= _omega * (fPlus[iPop] – fEqPlus[iPop]) + _omega2 * (fMinus[iPop] – fEqMinus[iPop]);

The omega related to the viscosity and the free parameter omega are switched in the collision.