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