mrt_d3q15 Function

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

This function returns mrt relaxation diagonal matrix for d3q15 layout Parameters are taken from: D’Humières, D., Ginzburg, I., Krafczyk, M., Lallemand, P., & Luo, L.-S. (2002). Multiple-relaxation-time lattice Boltzmann models in three dimensions. Philosophical Transactions. Series A, Mathematical, Physical, and Engineering Sciences, 360(1792), 437–51.


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