Validation Case: Compressible Flow In a de Laval Nozzle

This validation case belongs to fluid dynamics. The aim of this test is to validate the following parameters for a compressible, steady-state turbulent flow through a de Laval Nozzle:

• Mach number across the nozzle.
• Pressure across the nozzle.

The simulation results from SimScale were compared to analytical results obtained from methods presented in [1].

Geometry

A typical configuration of the de Laval nozzle with a non-smooth throat that was chosen as the geometry can be seen below:

The nozzle is axisymmetric and an angular slice of the complete geometry with a wedge angle of 18\(°\) was used. The following table contains the coordinates of the reference points:

Analysis Type and Mesh

Tool Type : OpenFOAM®

Analysis type : Compressible with k-omega SST turbulence model

Time dependency: Steady state

Mesh and Element types:

The hexahedral mesh was created externally, and then was imported on to the SimScale platform, with an inflation of \(y^+\) = 30 near the walls, as the ‘wall functions’ model was chosen. In order to keep the flow two-dimensional, the mesh was designed to have only one layer in the \(y\) direction.

Simulation Setup

Fluid:

• Custom He psi thermo type fluid:
  • Amount of substance: 1 mol
  • Molar mass: \(M_m\) = 28.9 \(kg \over \ kmol \)
  • Transport model: Sutherland
  • Reference viscosity: \(\nu_0\) = 1.716e-5 \(kg \over \ {s \times \ m}\)
  • Reference temperature: \(T_0\) = 273.15 \(K\)
  • \(T_s\) = 116 \(K\)
  • Thermo model: hconst
  • Specific heat: \(C_p\) = 1005 \(J \over\ {kg \times \ K}\)
  • Heat of formation: \(H_f\) = 0 \(J \over\ kg\)

Boundary Conditions:

The wedge boundary condition (BC) is used in computational fluid dynamics to define an axisymmetric situation. You can learn more about it here.

Reference Solution

The results are calculated from Bernoulli’s equation and the ideal gas equation. The speed of sound can be calculated as:

$$c = \sqrt{\frac{\gamma RT}{m}} \tag{1}$$

Where:

• \(\gamma\) is the specific heat ratio
• \(T\) is the temperature
• \(m\) is the molecular weight of the fluid
• \(R\) is the universal gas constant

The Mach number of a cutting plane be calculated as:
$$M = \frac{ U}{c} \tag{2}$$

Where:

• \(U\) is the average velocity magnitude of the plane

Finally, the static pressure can be computed as:

$$P_{stat} = P_{stag}\ – \frac{1}{2}\rho \times U^2 \ \tag{3}$$

Where:

• \(P_{stag}\) : the stagnation pressure, assumed to be 0.1301 \(MPa\)
• \(\rho\) : the density of the fluid

Result Comparison

A comparison of the Mach number and pressure variation in the nozzle obtained with SimScale with analytical results is given below:

A source of deviation between the calculated and analytical results exists because the latter is calculated with a one-dimensional hypothesis. This means that all parameters are assumed to be uniform in the radial direction, while in the following figures, it seems that there exists radial variation too in the flow variables.

One can observe the radial variation causing results to slightly deviate:

Last updated: January 7th, 2021