“…where ω i is the weight associated with e i and c s is the speed of sound, c s = x/ (3) t, and the local fluid density, ρ, and velocity, u, at the lattice node are given by ρ = i f i and ρu = i f i e i +τ * ρg, where g is the strength of an external force [37] and τ * = 0.5+3ν * t/ x 2 . The kinematic viscosity of the non-Newtonian fluid is described as ν * = 2 n−1 |Π D | n−1 2 ξ [34,38], in which ξ is the consistency, Π D is the second invariant of the rate of strain tensor, n is the fluid-type identifier, n < 1, n = 1, and n > 1 correspond to pseudoplastic (shear-thinning), Newtonian, and dilatant (shear-thickening) fluids, respectively. Π D is computed as…”