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.

Arguments

Type	Intent	Optional	Attributes		Name
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

Contents

None