where U is the control volume. r is the density of fluid. m is the dynamic viscosity of fluid. cp is the specific heat capacity of fluid. C0, A and B are Lagrange multipliers. V$U ¼ 0 is the constraint condition of mass conservation. lV2T rcpU$VT ¼ 0 is the constraint condition of energy conservation