Re: mrtCollision operator in OpenLB

[mathias]

Dear Eric,

1/ Actually 0,3, and 5 are ZERO:

template<typename T>
const T MRTD2Q9DescriptorBase<T>::S[MRTD2Q9DescriptorBase<T>::q_] =
{ T(), (T)1.1, (T)1.1, T(), (T)1.1, T(), (T)1.1, T(), T() };
// s7=s8 to have a shear viscosity nu
// and the bulk viscosity depends on s2.

2/ In openLB we save f_OLB = f – f^eq. So \sum_q f_OLB + 1 = rho.

Best
Mathias