Skip to main content
Have a personal or library account? Click to login
CFD Study of Ejector Nozzle Models Cover
Open Access
|Sep 2026

Figures & Tables

Fig.1.

Ejector nozzle thrust augmentation ratio vs. α*: 1 – first solution, ideal; 2 – second solution, ideal; 3 – second solution, with nozzle underexpansion losses from 1D theory; 4 – second solution, with nozzle underexpansion losses from 2D theory.

Fig. 2.

Ejector nozzle thrust augmentation ratio vs. flight Mach number (without taking into account losses).

Fig. 3.

Computational domain (configurations 1–5).

Fig. 4.

The five meshes considered: a – 20 mm; b – 15 mm; c – 10 mm; d – 7.5 mm; e – 5 mm.

Fig. 5.

Configurations considered (α*= 5).

Fig. 6.

Configurations considered (α*= 4).

Fig. 7.

Model boundaries.

Table 1.

Values of corresponding boundary parameters.

Boundary/ParameterTotal/static pressure, PaTurbulence intensity, %Mach numberFlow rate, kg/sTotal/static temperature, K
Left (Pressure-far-field)–/1013251.00.6––/288.15
Top (Pressure-far-field)–/1013251.00.6––/288.15
Right (Pressure outlet)–/1013251.0––308.89/–
AirOutlet (Mass-flow-outlet)–––0.99–
InletPrimary (Pressure-inlet)209787/–10.0––993.05/–
Table 2.

Parameter values on the corresponding meshes.

Mesh,Integral parameter, NRepresentative mesh cell size, m
1 (20 × 0.40 mm)395.651.24×10−2
2 (15 × 0.30 mm)396.211.01×10−2
3 (10 × 0.20 mm)396.317.46×10−3
4 (7.5 × 0.15 mm)394.665.94×10−3
5 (5 × 0.10 mm)395.604.23×10−3
Table 3.

Comparative estimation of the computational results.

ParameterMeshes 1–2–3Meshes 2–3–4Meshes 3–4–5
h1, m7.46×10−35.94×10−34.23×10−3
h2, m1.01×10−27.46×10−35.94×10−3
h3, m1.24×10−21.01×10−27.46×10−3
P1, N396.31394.66395.60
P2, N396.21396.31394.66
P3, N395.65396.21396.31
r321.2291.3501.255
r211.3501.2551.403
s1−1−1
p8.52517.7882.119
CR0.192−15.573−0.568
[0; 1) – monotonic convergence(−∞; −1) – oscillating divergence[−1; 0) – oscillating convergence
fex, N396.32394.63396.49
e32, %0.1400.0270.419
e21, %0.0270.4190.238
GCI32, %0.0360.0000.848
GCI21, %0.0030.0090.283
< 1% – for detailed calculations< 1% – for detailed calculations< 1% – for detailed calculations
AR1.0000.00031.463
[0.95; 1.05] – optimum(−∞; 0.95) – mesh should be improved(1.05; ∞) – mesh should be improved
hreq3.74×10−34.57×10−32.33×10−6
Table 4.

Differences of integral parameter calculations.

ParameterMesh 20 mmMesh 15 mmMesh 10 mmMesh 7.5 mmMesh 5 mm
Δm, %−4.438−4.173−4.076−4.605−4.329
Δv, %−17.564−16.510−16.232−17.036−16.542
Δp, %25.71724.28923.87624.89424.288
ΔNozzle Th, %−2.668−2.488−2.442−2.888−2.640
ΔTh, %−3.494−3.259−3.333−3.737−3.507
Fig. 8.

Pathline patterns for configurations 3 and 6 (α* = 5).

Fig. 9.

Fields of Mach number for the configurations considered (α* = 5).

Fig. 10.

Fields of Mach number for the configurations considered (α* = 4).

Fig. 11.

Fields of static pressure for the configurations considered (α* = 5).

Fig. 12.

Fields of static pressure for the configurations considered (α* = 4).

Fig. 13.

Fields of static temperature for some of the configurations considered (α* = 5).

Table 5.

Results of integral parameter calculations.

Parameter/ConfigurationSecondary flow Mach number at the entrance of mixing chamberPrimary flow rate, kg/sEntrainment ratioThrust, NAugmentation ratio, φ
Alperin0.4119—2.868—1.72
Base engine—1.0102—409.978—
1 (α* = 5)—0.96537—395.65—
2 (α* = 5)—0.96804—396.21—
3 (α* = 5)—0.96903—396.31—
4 (α* = 5)—0.96368—394.66—
5 (α* = 5)—0.96647—395.60—
6 (α* = 5)—0.97104—394.35—
7 (α* = 5)—0.97484—395.92—
8 (α* = 5)0.320820.968392.31028299.260.72994
9 (α* = 5)—0.61498—251.28—
10 (α* = 5)0.327500.973102.33550301.400.73516
11 (α* = 5)0.344700.966672.42548305.870.74607
12 (α* = 5)0.339970.964082.41586305.310.74470
13 (α* = 5)0.467430.2381713.57310227.010.55370
14 (α* = 5)0.342940.965642.40173304.530.74279
15 (α* = 5)0.354250.967382.42944298.820.72887
16 (α* = 5)0.345830.946222.83003365.900.89247
17 (α* = 5)0.347220.940402.97538382.910.93398
18 (α* = 5)0.344720.931972.98313383.830.93622
19 (α* = 5)0.354230.934402.95817385.230.93964
20 (α* = 5)0.316340.926132.84300381.470.93047
21 (α* = 5)0.339550.931602.83899386.710.94323
22 (α* = 5)0.344510.926282.85432386.430.94256
23 (α* = 5)0.373770.977952.66855407.850.99482
24 (α* = 5)0.402190.981812.62750405.870.98997
25 (α* = 5)0.399480.981022.61869401.610.97958
26 (α* = 5)0.371020.979632.64032410.451.00115
27 (α* = 5)0.370750.980062.62426404.570.98681
28 (α* = 5)0.367160.980292.60928402.620.98206
29 (α* = 5)0.374230.977422.60485406.130.99062
30 (α* = 5)0.368180.979932.61096402.310.98130
31 (α* = 5)0.369760.980772.61755403.160.98337
32 (α* = 5)0.363900.979552.54888378.820.92399
Alperin0.379—2.008—1.62
Base engine—1.0102—409.978—
1 (α* = 4)0.356920.978622.04513402.430.98158
2 (α* = 4)0.359620.979302.05995398.650.97238
3 (α* = 4)0.358930.978052.06253399.830.97525
4 (α* = 4)0.356520.977932.04692404.070.98559
5 (α* = 4)0.357500.978542.04527397.830.97036
6 (α* = 4)0.358140.978482.05250399.170.97363
7 (α* = 4)0.366320.979722.03578405.580.98928
8 (α* = 4)0.365160.979182.04474406.690.99198
9 (α* = 4)0.350500.976292.06375400.550.97699
10 (α* = 4)0.350070.978772.05843400.610.97715
11 (α* = 4)0.344510.974902.03757402.300.98127
12 (α* = 4)0.350060.976482.06037403.690.98467
13 (α* = 4)0.350340.980142.04143402.450.98165
14 (α* = 4)0.356510.978612.04413398.080.97097
15 (α* = 4)0.355970.980162.03369395.000.96346
16 (α* = 4)0.355250.979892.04758396.510.96716
17 (α* = 4)0.336910.979772.04349395.180.96389
18 (α* = 4)0.383280.964842.33288464.501.13298
19 (α* = 4)0.409150.948742.63502512.451.24994
20 (α* = 4)0.095720.962090.61172375.470.91584
21 (α* = 4)0.221710.968761.35933400.380.97658
22 (α* = 4)0.580580.980282.89229385.550.94040
23 (α* = 4)0.716140.982273.14938315.270.76900
24 (α* = 4)0.708020.976913.19422167.340.40818
25 (α* = 4)0.736420.982363.1746713.620.03323
26 (α* = 4)0.352330.999941.94880390.920.95351
27 (α* = 4)0.348800.994591.98429399.940.97553
28 (α* = 4)0.376240.982982.32974521.441.00290
29 (α* = 4)0.391630.970162.59136612.900.98818
Fig. 14.

Pathline patterns for the configurations considered (α* = 5).

Fig. 15.

Pathline patterns for the configurations considered (α* = 4).

Fig. 16.

Thrust augmentation ratio (φ) vs. the ratio of the common exhaust unit exit area to the mixing chamber area (αD) (configurations 20, 21, 2, 22, 23, 24, 25) (α* = 4).

Fig. 17.

Thrust augmentation ratio (φ) vs. the ratio of the primary convergent-divergent nozzle exit area to its throat area (αDp) (configurations 2, 27, 26) (α* = 4).

Fig. 18.

Fields of Mach number for configurations 28, 28w, 29 and 29w (α* = 4).

Fig. 19.

Thrust augmentation ratio (φ) vs. primary nozzle pressure ratio πN (the value at the point πN = 4 was interpolated).

Language: English
Page range: 41 - 68
Submitted on: May 10, 2026
Accepted on: Jul 21, 2026
Published on: Sep 28, 2026
Published by: ŁUKASIEWICZ RESEARCH NETWORK – INSTITUTE OF AVIATION
In partnership with: Paradigm Publishing Services

© 2026 Ruslan Tsukanov, Sergiy Yepifanov, published by ŁUKASIEWICZ RESEARCH NETWORK – INSTITUTE OF AVIATION
This work is licensed under the Creative Commons Attribution 4.0 License.