\(\frac{\partial\left(\rho_{L}\varepsilon_{L}u_{L}\right)}{\partial t}+\nabla.\left(\rho_{L}\varepsilon_{L}u{{}_{L}u}_{L}-\mu_{L}\varepsilon_{L}\left(\nabla u_{L}+\left(\nabla u_{L}\right)^{T}\right)\right)=-\varepsilon_{L}\nabla p-M_{G,L}+\rho_{L}\varepsilon_{L}g\), and (3)