$\begingroup$

After doing quite some research I think I have the correct equations. Unfortunately all of the equations above (except for the continuity) are wrong.

Linear momentum:

The linear momentum conservation equation in the most compact and general form can be written as:

$$\rho \frac{D \check{v}}{D t}=\check{

abla} \check{\sigma} \tag{5}$$

Where the $\check{\sigma}$ is the Cauchy stress tensor $\boldsymbol{\sigma}$ in matrix form. The Cauchy stress tensor can be divided into a volumetric/dilatational part and a deviatoric part:

$$\check{\sigma}= \check{\tau} -P\check{I} \tag{6}$$

In a perfect fluid P is the hydrostatic pressure and $\boldsymbol{\tau}$ is the viscous stress tensor.

Constitutive equation:

In an ideal Newtonian fluid viscous stress tensor can be calculated from:

$$ \check{\tau}=\eta\left( \check{

abla}^T \check{v}+ \left( \check{

abla}^T \check{v} \right)^T \right) +\lambda \left(\check{

abla}\check{v}^T\right) \check{I} \tag{7}$$

Where $\eta$ is the absolute/dynamic viscosity and $\lambda$ is the bulk viscosity defined as:

$$\lambda=\kappa-\frac{2}{3}\eta \tag{8}$$

Where $\kappa$ is the dilation/expansion viscosity, which from Stokes’ hypothesis for monatomic gases at low density it is negligible impyling:

$$ \check{\tau}=\eta\left( \check{

abla}^T \check{v}+ \left( \check{

abla}^T \check{v} \right)^T -\frac{2}{3} \left(\check{

abla}.\check{v}^T\right) \check{I} \right) \tag{9}$$

$\eta$ is not necessarily constant. There are some empirical models for example Sutherland:

$$\eta \approx C_S \frac{T^{\left(3/2\right)}}{T+T_S} \tag{10}$$

Where $C_S$ and $T_S$ are constants.

Combining the equations 5, 6 and 9 and also considering the axisymmetry, ideal gas and steady-state assumptions equations of linear momentum in axial and radial directions can be written as:

$$\begin{gather*} \rho \left( v_r \frac{\partial v_r}{\partial r} +v_x \frac{\partial v_r}{\partial x} \right) = \frac{\partial}{\partial r} \left( \eta \left( -\frac{2}{3}\left( \frac{1}{r}\frac{\partial}{\partial r} \left( r v_r \right) +\frac{\partial v_x}{\partial x} \right) +2 \frac{\partial v_r}{\partial r} \right) \right) \\ +\frac{\partial}{\partial x}\left( \eta\left( \frac{\partial v_r}{\partial x}+\frac{\partial v_x}{\partial r} \right) \right)+ \frac{2 \eta}{r}\left( \frac{\partial v_r}{\partial r}- \frac{v_r}{r} \right) -\frac{\partial }{\partial r}\left(\mathring{R}\rho T \right) \end{gather*} \tag{11}$$

In radial direction and

$$\begin{gather*} \rho \left( v_r \frac{\partial v_x}{\partial r} +v_x \frac{\partial v_x}{\partial x} \right) = \frac{\partial}{\partial x} \left( \eta \left( -\frac{2}{3}\left( \frac{1}{r}\frac{\partial}{\partial r} \left( r v_r \right) +\frac{\partial v_x}{\partial x} \right) +2 \frac{\partial v_x}{\partial x} \right) \right) \\ +\frac{\partial}{\partial r}\left( \eta\left( \frac{\partial v_r}{\partial x}+\frac{\partial v_x}{\partial r} \right) \right)+ \frac{ \eta}{r}\left( \frac{\partial v_r}{\partial x}+ \frac{\partial v_x}{\partial r} \right) -\frac{\partial }{\partial x}\left(\mathring{R}\rho T \right) \end{gather*} \tag{12}$$

In axial direction. (ref1, ref2)

Energy equation:

I'm still not completely sure about the correct form of the energy equation as I have explained here. But assuming my post here is valid then the energy equation in compact form can be written as:

$$ \rho \frac{D e}{D t}=\check{\sigma} : \check{

abla}^T \check{v} \tag{13}$$

Considering that for an ideal gas $e=c_v T$ (more info here and here) and axisymmetry and steady-state assumptions it expands to:

$$ \begin{gather*} \rho c_v\left( v_x \frac{\partial T}{\partial x} + v_r \frac{\partial T}{\partial r} \right)+ \mathring{R} \rho T \left( \frac{1}{r}\frac{\partial}{\partial r} \left( r v_r \right)+ \frac{\partial v_x}{\partial x} \right)= \\ \eta \left( 2 \left( \frac{\partial v_x}{\partial x} \right)^2+ 2 \left( \frac{\partial v_r}{\partial r} \right)^2+ \left( \frac{\partial v_r}{\partial x}+ \frac{\partial v_x}{\partial r} \right)^2 -\frac{2}{3}\left( \frac{1}{r} \frac{\partial}{\partial r}\left( r v_r \right) + \frac{\partial v_x}{\partial x} \right)^2 \right) \end{gather*} \tag{14}$$