mrt_d2q9 Function

public pure function mrt_d2q9(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

Called by

proc~~mrt_d2q9~~CalledByGraph proc~mrt_d2q9 mrt_d2q9 proc~mrt_advrel_d2q9 mrt_advRel_d2q9 proc~mrt_advrel_d2q9->proc~mrt_d2q9 proc~mrt_advrel_d2q9_incomp mrt_advRel_d2q9_incomp proc~mrt_advrel_d2q9_incomp->proc~mrt_d2q9