> Diff(rho, t) = -(Diff(rho * u, x) + Diff(rho * v, y)); d /d \ /d \ -- rho = -|-- rho u| - |-- rho v| dt \dx / \dy / > Diff(u, t) = -u * Diff(u, x) - v * Diff(u, y) + (1 / rho) * (-Diff(p, x) + Diff(tau[xx], x) + Diff(tau[yx], y)); /d \ /d \ /d \ -|-- p| + |-- tau[xx]| + |-- tau[yx]| d /d \ /d \ \dx / \dx / \dy / -- u = -u |-- u| - v |-- u| + ------------------------------------- dt \dx / \dy / rho > Diff(v, y) = -u * Diff(v, x) - v * Diff(v, y) + (1 / rho) * (-Diff(p, y) + Diff(tau[xy], x) + Diff(tau[yy], y)); /d \ /d \ /d \ -|-- p| + |-- tau[xy]| + |-- tau[yy]| d /d \ /d \ \dy / \dx / \dy / -- v = -u |-- v| - v |-- v| + ------------------------------------- dy \dx / \dy / rho > Diff(e + V^2 / 2, t) = -u * Diff(e + V^2 / 2, x) - v * Diff(e + V^2 / 2, y) + Diff(q, t) + (1 / rho) * (Diff(k * Diff(T, x), x) + Diff(k * Diff(T, y), y) - Diff(p * u, x) - Diff(p * v, y) + Diff(u * tau[xx], x) + Diff(u * tau[yx], y) + Diff(v * tau[xy], x) + Diff(v * tau[yy], y)); d 2 /d 2 \ /d 2 \ -- (e + 1/2 V ) = -u |-- (e + 1/2 V )| - v |-- (e + 1/2 V )| dt \dx / \dy / /d \ //d /d \\ /d /d \\ /d \ + |-- q| + ||-- k |-- T|| + |-- k |-- T|| - |-- p u| \dt / \\dx \dx // \dy \dy // \dx / /d \ /d \ /d \ - |-- p v| + |-- u tau[xx]| + |-- u tau[yx]| \dy / \dx / \dy / /d \ /d \\ + |-- v tau[xy]| + |-- v tau[yy]||/rho \dx / \dy // >