mrt_d2q9_incomp Function

private pure function mrt_d2q9_incomp(omegaKine, omegaBulk, QQ) result(s_mrt)

This function returns mrt relaxation diagonal matrix for d2q9 layout Parameters are taken from: Lallemand, P., & Luo, L. (2000). Theory of the lattice boltzmann method: dispersion, dissipation, isotropy, galilean invariance, and stability. Physical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics, 61(6 Pt A), 6546–62.

Another paper for d2q9 with different set of parameters Chen, S., Peng, C., Teng, Y., Wang, L. P., & Zhang, K. (2016). Improving lattice Boltzmann simulation of moving particles in a viscosity flow using local grid refinement. Computers and Fluids, 136, 228–246.


real(kind=rk), intent(in) :: omegaKine

omega related to kinematic viscosity

real(kind=rk), intent(in) :: omegaBulk

omega related to bulk viscosity

integer, intent(in) :: QQ

number of directions

Return Value real(kind=rk)(QQ)

output mrt diagonal matrix