Navier Stokes Equation In Cylindrical Coordinates

Navier Stokes Equation in Cylindrical Coordinates

The Navier‑Stokes equations describe how the velocity field of a fluid evolves under the influence of internal stresses, external body forces, and inertia. While the Cartesian form is familiar, many engineering problems—such as flow in pipes, rotating machinery, or vortex dynamics—exhibit natural symmetry around an axis. In these cases expressing the governing equations in cylindrical coordinates ((r,\theta,z)) simplifies boundary‑condition implementation and highlights the physical mechanisms that dominate the motion. Below we derive the full set of Navier‑Stokes equations in cylindrical coordinates, discuss the accompanying continuity equation, explore common simplifications, and outline practical solution strategies.

1. Derivation Overview

Starting from the vector form of the momentum balance

[ \rho\left(\frac{\partial \mathbf{v}}{\partial t}+ \mathbf{v}!\cdot!\nabla \mathbf{v}\right)= -\nabla p + \mu \nabla^{2}\mathbf{v} + \mathbf{f}, ]

where (\rho) is density, (\mu) dynamic viscosity, (p) pressure, and (\mathbf{f}) body‑force per unit volume, we substitute the cylindrical‑coordinate expressions for the gradient, divergence, and Laplacian operators. The velocity vector is written as

[ \mathbf{v}=v_r,\mathbf{e}r+v\theta,\mathbf{e}_\theta+v_z,\mathbf{e}_z, ]

with unit vectors (\mathbf{e}r,\mathbf{e}\theta,\mathbf{e}_z) that vary with (\theta). Because these unit vectors are not constant, extra terms appear when differentiating the velocity components—a key feature that distinguishes the cylindrical form from its Cartesian counterpart.

2. Navier‑Stokes Equations in Cylindrical Coordinates

After carrying out the differentiation and collecting terms, the momentum balance yields three scalar equations, one for each direction.

2.1 Radial ((r)) Component

[\boxed{ \begin{aligned} \rho\Bigg(&\frac{\partial v_r}{\partial t}

• v_r\frac{\partial v_r}{\partial r}
• \frac{v_\theta}{r}\frac{\partial v_r}{\partial \theta}
• v_z\frac{\partial v_r}{\partial z}

• \frac{v_\theta^{2}}{r}\Bigg) \ =&-\frac{\partial p}{\partial r} +\mu\Bigg[ \frac{1}{r}\frac{\partial}{\partial r}!\left(r\frac{\partial v_r}{\partial r}\right) +\frac{1}{r^{2}}\frac{\partial^{2} v_r}{\partial \theta^{2}} +\frac{\partial^{2} v_r}{\partial z^{2}} -\frac{v_r}{r^{2}} -\frac{2}{r^{2}}\frac{\partial v_\theta}{\partial \theta} \Bigg]

• f_r . \end{aligned}} ]

Key points

• The term (-v_\theta^{2}/r) is the centrifugal acceleration arising from curvature of the (\theta)‑direction.
• The viscous term contains extra geometric contributions (-v_r/r^{2}) and (-2/r^{2},\partial v_\theta/\partial \theta) that vanish in Cartesian coordinates.

2.2 Azimuthal ((\theta)) Component

[ \boxed{ \begin{aligned} \rho\Bigg(&\frac{\partial v_\theta}{\partial t}

• v_r\frac{\partial v_\theta}{\partial r}
• \frac{v_\theta}{r}\frac{\partial v_\theta}{\partial \theta}
• v_z\frac{\partial v_\theta}{\partial z}
• \frac{v_r v_\theta}{r}\Bigg) \ =&-\frac{1}{r}\frac{\partial p}{\partial \theta} +\mu\Bigg[ \frac{1}{r}\frac{\partial}{\partial r}!\left(r\frac{\partial v_\theta}{\partial r}\right) +\frac{1}{r^{2}}\frac{\partial^{2} v_\theta}{\partial \theta^{2}} +\frac{\partial^{2} v_\theta}{\partial z^{2}} -\frac{v_\theta}{r^{2}} +\frac{2}{r^{2}}\frac{\partial v_r}{\partial \theta} \Bigg]
• f_\theta . \end{aligned}} ]

Key points

• The Coriolis‑like term (+v_r v_\theta/r) couples radial and azimuthal motions.
• Viscous diffusion includes the (+2/r^{2},\partial v_r/\partial \theta) term, reflecting shear generated by radial gradients in the (\theta) direction.

2.3 Axial ((z)) Component

[ \boxed{ \begin{aligned} \rho\Bigg(&\frac{\partial v_z}{\partial t}

• v_r\frac{\partial v_z}{\partial r}
• \frac{v_\theta}{r}\frac{\partial v_z}{\partial \theta}
• v_z\frac{\partial v_z}{\partial z}\Bigg) \ =&-\frac{\partial p}{\partial z} +\mu\Bigg[ \frac{1}{r}\frac{\partial}{\partial r}!\left(r\frac{\partial v_z}{\partial r}\right) +\frac{1}{r^{2}}\frac{\partial^{2} v_z}{\partial \theta^{2}} +\frac{\partial^{2} v_z}{\partial z^{2}} \Bigg]
• f_z . \end{aligned}} ]

Key points

• No extra geometric terms appear because the unit vector (\mathbf{e}_z) is constant; the axial component resembles the Cartesian form apart from the cylindrical Laplacian.

3. Continuity Equation in Cylindrical Coordinates

Mass conservation for a compressible fluid reads

[ \frac{\partial \rho}{\partial t} +\nabla!\cdot!(\rho\mathbf{v})=0 . ]

In cylindrical coordinates the divergence expands to

[ \boxed{ \frac{\partial \rho}{\partial t} +\frac{1}{r}\frac{\partial}{\partial r}!\left(r\rho v_r\right) +\frac{1}{r}\frac{\partial}{\partial \theta}!\left(\rho v_\theta\right) +\frac{\partial}{\partial z}!\left(\rho v_z\right)=0 . } ]

For incompressible flow ((\rho=) constant) this simplifies to

[ \frac{1}{r}\frac{\partial}{\partial r}!\left(r v_r\right) +\frac{1}{r}\frac{\partial v_\theta}{\partial \theta} +\frac{\partial v_z}{\partial z}=0 . ]

4. Common Simplifications and Applications

4.1 Axisymmetric Flow ((\partial/\partial\theta =0), (v_\theta=0))

Many engineering problems—such as laminar flow in a straight pipe or a steady jet impinging on a plate—possess axisymmetry. Setting (\partial/\partial\theta=0) and (v_\theta=0) eliminates all (\theta)‑derivatives and the centrifugal/Coriolis terms, leaving:

3. Simplified Equationsfor Axisymmetric Flow

Many practical fluid flow problems exhibit symmetry around the z-axis, known as axisymmetric flow. This symmetry allows significant simplification of the governing equations. The defining characteristics are:

1. No θ-dependence: (\frac{\partial}{\partial \theta} = 0)
2. No azimuthal velocity: (v_\theta = 0)

Applying these conditions to the momentum equations:

3.1 Radial Component (r)

The original radial momentum equation contained terms involving (\theta):

• The Coriolis-like term: (+\frac{v_r v_\theta}{r})
• The viscous diffusion term: (+\frac{2}{r^{2}}\frac{\partial v_r}{\partial \theta})

Setting (v_\theta = 0) and (\frac{\partial}{\partial \theta} = 0) eliminates these terms entirely. The simplified radial momentum equation becomes:

[ \boxed{ \rho\Bigg(\frac{\partial v_r}{\partial t} + v_r \frac{\partial v_r}{\partial r} + v_z \frac{\partial v_r}{\partial z}\Bigg) = -\frac{\partial p}{\partial r} + \mu \left[ \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial v_r}{\partial r} \right) + \frac{\partial^{2} v_r}{\partial z^{2}} \right] + f_r } ]

3.2 Axial Component (z)

The axial momentum equation originally had no (\theta)-dependent terms. Setting (v_\theta = 0) and (\frac{\partial}{\partial \theta} = 0) removes all (\theta)-related terms. The simplified axial momentum equation is:

[ \boxed{ \rho\Bigg(\frac{\partial v_z}{\partial t} + v_r \frac{\partial v_z}{\partial r} + v_z \frac{\partial v_z}{\partial z}\Bigg) = -\frac{\partial p}{\partial z} + \mu \left[ \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial v_z}{\partial r} \right) + \frac{\partial^{2} v_z}{\partial z^{2}} \right] + f_z } ]

3.3 Continuity Equation

The continuity equation for incompressible flow ((\rho = \text{constant})) simplifies further under axisymmetry:

[ \frac{1}{r} \frac{\partial}{\partial r} (r v_r) + \frac{\partial v_z}{\partial z} = 0 ]

For compressible flow, the general continuity equation remains:

[ \frac{\partial \rho}{\partial t} + \frac{1}{r} \frac{\partial}{\partial r} (r \rho v_r) + \frac{\partial}{\partial z} (\rho v_z) = 0 ]

4. Applications and Significance

Axisymmetric flow analysis is fundamental in numerous engineering disciplines:

• Internal Flows: Laminar flow in straight pipes, annular ducts, and axisymmetric nozzles.
• External Flows: Steady jet impingement on a flat plate, flow around axisymmetric bodies (e.g., spheres, cylinders).
• Rotating Flows: Centrifugal pumps, turbines, and other rotating machinery where the flow develops axisymmetric patterns.
• Geophysical Flows: Atmospheric and oceanic circulation models often use axisymmetric approximations for large-scale phenomena.

The reduction from the full 3D cylindrical equations to the simplified

The reduction from the full three‑dimensional cylindrical equations to the axisymmetric set not only lowers the algebraic complexity but also reveals a clear physical hierarchy: the dominant transport mechanisms are now confined to two independent variables, (r) and (z). This simplification enables analytical progress in regimes where additional constraints can be imposed without sacrificing essential dynamics.

4.1 Steady, Fully‑Developed Axial Flow

When the flow is steady ((\partial/\partial t =0)) and fully developed in the axial direction ((\partial/\partial z =0) for all velocity components), the governing equations collapse to ordinary differential equations in the radial coordinate alone. For a pressure‑driven pipe flow, the axial momentum balance reduces to

[ 0 = -\frac{dp}{dz} + \mu\left[\frac{1}{r}\frac{d}{dr}!\left(r\frac{dv_z}{dr}\right)\right] + f_z . ]

Integrating twice with respect to (r) and applying the no‑slip condition at the wall ((v_z=0) at (r=R)) together with the regularity condition at the axis ((dv_z/dr|_{r=0}=0)) yields the classic Hagen–Poiseuille parabolic profile. The same procedure, when coupled with the simplified continuity relation, provides the volumetric flow rate (Q = \pi R^{4},(-\Delta p)/(8\mu L)), a result that is recovered instantly from the reduced system without revisiting the full three‑dimensional derivation.

4.2 Swirl‑Free Jet Impingement

In many process‑engineering applications a jet of fluid issues from a circular nozzle and strikes a planar surface. The velocity field can be approximated as axisymmetric and steady, with the radial component dominating near the nozzle exit while the axial component decays rapidly after impact. By assuming that the radial diffusion term balances the pressure gradient in the radial momentum equation and neglecting inertial convection, the governing balance becomes

[ 0 = -\frac{\partial p}{\partial r} + \mu\left[\frac{1}{r}\frac{\partial}{\partial r}!\left(r\frac{\partial v_r}{\partial r}\right)\right] . ]

Integrating once yields a shear stress distribution that is directly proportional to the imposed mass flux, and a second integration furnishes the radial velocity field. The resulting expression for the wall‑shear stress (\tau_w) can be expressed in closed form as a function of the jet radius (R_j) and the discharge coefficient, offering a quick design tool for nozzle shaping and impact‑load prediction.

4.3 Dimensionless Form and Similarity Solutions

Introducing the characteristic length (L_c) (e.g., pipe radius or nozzle exit diameter) and velocity scale (U_c) (e.g., mean axial speed), the reduced equations can be rendered dimensionless through the substitutions

[ \hat{r}= \frac{r}{L_c},\qquad \hat{z}= \frac{z}{L_c},\qquad \hat{v}_r = \frac{v_r}{U_c},\qquad \hat{v}_z = \frac{v_z}{U_c},\qquad \hat{p}= \frac{p}{\rho U_c^{2}} . ]

The resulting nondimensional groups are the Reynolds number (Re = \rho U_c L_c/\mu) and the ratio of body‑force to inertial forces (Fr = U_c^{2}/(gL_c)) (if gravity is present). In many canonical problems the governing system admits similarity variables that reduce the partial differential equations to ordinary differential equations. For instance, in unsteady boundary‑layer growth over a flat plate aligned with the axis, the similarity variable (\eta = z\sqrt{\nu t}/r) transforms the axisymmetric momentum and continuity equations into the well‑known Falkner–Skan equation, whose solution provides the time‑dependent velocity profile without resorting to full transient simulations.

4.4 Numerical Strategies

Because the axisymmetric system retains only two spatial dimensions, it is ideally suited for finite‑volume or finite‑element discretisations on structured or unstructured grids that exploit rotational symmetry. A common practice is to store flow variables at cell centres and employ axisymmetric fluxes, where the radial coordinate appears as a geometric weighting factor in the continuity and momentum balances. Implicit time‑marching schemes (e.g., Crank–Nicolson) combined with pressure‑correction algorithms (SIMPLE, PISO) converge rapidly, especially when the grid is refined near the axis to resolve the singular behaviour of the (1/r) terms. Adaptive mesh refinement (AMR) can be deployed to capture sharp gradients during transient events such as vortex shedding or shock impingement, while still maintaining the reduced computational cost relative to a full 3‑D simulation.

Conclusion

Axisymmetric flow analysis provides a powerful bridge between the idealised, analytically tractable cases of inviscid potential flow and the computationally intensive, fully three‑dimensional Navier–Stokes simulations of contemporary engineering problems. By exploiting the inherent symmetry of the flow field, the governing equations are reduced

...to a manageable two-dimensional problem, significantly reducing computational demands without sacrificing the ability to capture essential flow features. This approach is particularly valuable for preliminary design iterations, rapid prototyping, and understanding the fundamental behavior of nozzle geometries. The development of sophisticated numerical techniques, including adaptive mesh refinement and implicit time-stepping schemes, further enhances the accuracy and efficiency of axisymmetric simulations.

In essence, axisymmetric analysis offers a pragmatic and efficient pathway to understanding and optimizing nozzle designs. It allows engineers to explore a wide range of flow conditions and geometries with a balance of accuracy and computational cost, accelerating the design process and enabling the development of innovative and high-performance nozzles for various applications. As computational power continues to increase, the role of axisymmetric analysis will likely expand, serving as a crucial tool for addressing complex flow challenges in aerospace, automotive, and other engineering disciplines.