\(\frac{\partial\left(\varepsilon_{L}\rho_{L}\right)}{\partial t}+\nabla.\left(\varepsilon_{L}\rho_{L}u_{L}\right)=-S_{\text{GL}}\) (1)