# Ejector (G)

**Libraries:**

Simscape /
Fluids /
Gas /
Turbomachinery

## Description

The Ejector (G) block models an ejector in a gas network. The ejector forces the high-pressure primary flow through a nozzle with a specific throat area, which creates a low-pressure jet that draws in the secondary flow. The expansion of the primary jet creates a throat for the secondary flow. The block assumes that the two flows mix under uniform pressure and the combined flow exits through a diffuser to recover pressure. The primary and secondary flows must be the same gas.

You can use an ejector as a pump or compressor to drive a secondary fluid. Unlike a pump or compressor, there is no external source of mechanical energy, and rather the high pressure primary fluid is the energy source that drives the secondary fluid.

### Ejector Geometry

The high-pressure primary flow enters the ejector at location p and flows through a
convergent-divergent nozzle with throat area
*S _{pt}*, which is the value of the

**Primary nozzle throat area**parameter. The ejector chokes the primary flow at the throat (location pt) and accelerates the flow to a low-pressure supersonic jet at the primary nozzle exit (location pn). The low pressure of the supersonic primary flow jet draws in the secondary flow at location s around the nozzle.

The supersonic primary flow jet continues to expand and decrease in pressure beyond the nozzle exit until it reaches its largest area at the expanded primary flow jet (location pe). The expansion of the supersonic primary flow jet creates a fictitious aerodynamic throat for the secondary flow at location st, which is at the same location as pe. The secondary flow also chokes at the aerodynamic throat. After passing the aerodynamic throat, the primary and secondary flow mix under uniform pressure (location m). If the mixed flow is still supersonic, the flow causes a shock in the mixing chamber at location s and the diffuser (location d) decelerates the mixed flow to raise the outlet pressure.

This schematic describes the ejector where:

Location p is the primary inlet.

Location s is the secondary inlet.

Location pt is the primary nozzle throat.

Location pn is the primary nozzle exit.

Location pe is the expanded primary flow jet.

Location st is the secondary aerodynamic throat.

Location m is the mixed flow.

Location ms is the shock in the mixing chamber.

Location d is the diffuser outlet.

### Mass Flow Rate

The ejector accelerates the primary gas flow through the nozzle. The mass flow rate of the primary flow is

$${\dot{m}}_{p}={S}_{pt}\sqrt{\left(\frac{2\gamma}{\gamma -1}\right){\eta}_{p}{p}_{p}{\rho}_{p}\left[{\left(\frac{{p}_{m}}{{p}_{p}}\right)}^{\frac{2}{\gamma}}-{\left(\frac{{p}_{m}}{{p}_{p}}\right)}^{\frac{\gamma +1}{\gamma}}\right]},$$

where:

*S*is the value of the_{pt}**Primary nozzle throat area**parameter.$$\gamma =\frac{{c}_{p}}{{c}_{v}}$$ is the ratio of heat capacities where

*c*is the specific heat at constant pressure and_{p}*c*is the specific heat at constant volume._{v}*η*is the value of the_{p}**Efficiency for primary flow through nozzle**parameter.*p*is the pressure at the primary flow inlet._{p}*p*is the pressure after the primary and secondary flow mix._{m}*ρ*is the density at the primary flow inlet._{p}

The primary flow draws the secondary gas flow into the suction chamber and accelerates it towards the fictitious aerodynamic throat created by the expanding primary jet. The mass flow rate of the secondary flow is

$${\dot{m}}_{s}={S}_{st}\sqrt{\left(\frac{2\gamma}{\gamma -1}\right){\eta}_{s}{p}_{s}{\rho}_{s}\left[{\left(\frac{{p}_{m}}{{p}_{s}}\right)}^{\frac{2}{\gamma}}-{\left(\frac{{p}_{m}}{{p}_{s}}\right)}^{\frac{\gamma +1}{\gamma}}\right]},$$

where:

*S*is the area of the fictitious aerodynamic throat. The block calculates the value based on the expansion of the primary flow jet._{st}*η*is the value of the_{s}**Efficiency for secondary suction flow**parameter.*p*is the pressure at the secondary flow inlet._{s}*ρ*is the density at the secondary flow inlet._{s}

When $$\frac{{p}_{m}}{{p}_{p}}<{B}_{crit},$$ where $${B}_{crit}={\left(\frac{2}{\gamma +1}\right)}^{\frac{\gamma}{\gamma -1}},$$ the primary flow becomes choked and the mass flow rate reaches the maximum value

$${\dot{m}}_{p}^{*}={S}_{pt}\sqrt{{\eta}_{p}{p}_{p}{\rho}_{p}\gamma {\left(\frac{2}{\gamma +1}\right)}^{\frac{\gamma +1}{\gamma -1}}.}$$

When $$\frac{{p}_{m}}{{p}_{s}}<{B}_{crit},$$ the secondary flow becomes choked and the mass flow rate reaches the maximum value

$${\dot{m}}_{s}{}^{*}={S}_{st}\sqrt{{\eta}_{s}{p}_{s}{\rho}_{s}\gamma {\left(\frac{2}{\gamma +1}\right)}^{\frac{\gamma +1}{\gamma -1}}}.$$

### Primary and Secondary Flow Relations

The flow area expansion of the supersonic primary flow jet comes from the relation

$$\frac{{S}_{pe}}{{S}_{pt}}={\Psi}_{pe}\sqrt{\frac{\left(\frac{\gamma -1}{2}\right){\left(\frac{2}{\gamma +1}\right)}^{\frac{\gamma +1}{\gamma -1}}}{{\left(\frac{{p}_{m}}{{p}_{p}}\right)}^{\frac{2}{\gamma}}-{\left(\frac{{p}_{m}}{{p}_{p}}\right)}^{\frac{\gamma +1}{\gamma}}}},$$

where *ψ _{pe}* is the value of the

**Efficiency for primary flow expansion**parameter.

The block uses the area of the expanded supersonic primary flow jet,
*S _{pe}*, to calculate the flow area of
the fictitious secondary flow aerodynamic throat,

*S*

_{st}$$\frac{{S}_{st}}{{S}_{pt}}=\frac{{S}_{m}}{{S}_{pt}}-\mathrm{max}\left\{\frac{{S}_{pe}}{{S}_{pt}},\frac{{S}_{pm}}{{S}_{pt}}\right\},$$

where:

$$\frac{{S}_{m}}{{S}_{pt}}$$ is the value of the

**Area ratio of mixing chamber to throat**parameter.$$\frac{{S}_{pn}}{{S}_{pt}}$$ is the value of the

**Area ratio of nozzle exit to throat**parameter.

The primary flow velocity of the expanded jet is

$${\text{v}}_{pe}=\sqrt{\frac{2\gamma}{\gamma -1}\left(\frac{{p}_{p}}{{\rho}_{p}}\right)\left[1-{\left(\frac{{p}_{m}}{{p}_{p}}\right)}^{\frac{\gamma -1}{\gamma}}\right]}.$$

The secondary flow velocity at the aerodynamic throat is

$${\text{v}}_{st}=\sqrt{\frac{2\gamma}{\gamma -1}\left(\frac{{p}_{s}}{{\rho}_{s}}\right)\left[1-{\left(\frac{{p}_{m}}{{p}_{s}}\right)}^{\frac{\gamma -1}{\gamma}}\right]}.$$

The secondary flow velocity is limited to the choked value. When $$\frac{{p}_{m}}{{p}_{s}}<{B}_{crit},$$, the secondary flow reaches the maximum value of

$${\text{v}}_{st}^{*}=\sqrt{\frac{2\gamma}{\gamma +1}\left(\frac{{p}_{s}}{{\rho}_{s}}\right)}.$$

### Mixing Model

After the fictitious aerodynamic throat, the primary and secondary flows mix under uniform pressure. The mixed flow temperature is

$${\text{T}}_{m}=\frac{{\dot{m}}_{p}{T}_{p}+{\dot{m}}_{s}{T}_{s}}{{\dot{m}}_{m}}-\frac{\gamma -1}{2\gamma R}{v}_{m}^{2},$$

where:

*T*is the primary flow inlet temperature. The block assumes that this value is equal to the primary flow stagnation temperature because kinetic energy at the inlet is negligible relative to the inside of the ejector._{p}*T*is the mixed flow temperature._{m}*T*is the secondary flow inlet temperature. The block assumes that this value is equal to the secondary flow stagnation temperature because kinetic energy at the inlet is negligible relative to the inside of the ejector._{s}*R*is the specific gas constant.*v*is the flow velocity after the primary and secondary flow mix._{m}$$\dot{m}$$

_{m}is the mass flow rate after the primary and secondary flow mix, meaning $${\dot{m}}_{m}={\dot{m}}_{p}+{\dot{m}}_{s}.$$

The momentum balance in the mixing region is

$${p}_{m}{S}_{m}+{\dot{m}}_{m}{v}_{m}={p}_{m}{S}_{pe}+{\dot{m}}_{p}{v}_{pe}+{p}_{m}{S}_{st}+{\dot{m}}_{s}{v}_{st},$$

where *S _{m}* is the
cross-sectional area after the primary and secondary flow mix.

The pressure terms drop out due to the uniform pressure mixing assumption, and the mixed flow velocity is

$${\text{v}}_{m}={\psi}_{m}\frac{{\dot{m}}_{p}{v}_{pe}+{\dot{m}}_{s}{v}_{st}}{{\dot{m}}_{m}},$$

where *ψ _{m}* is the value
of the

**Efficiency for mixing**parameter.

The mixed flow Mach number is

$${\text{M}}_{m}=\frac{{v}_{m}}{\sqrt{\gamma R{T}_{m}}}.$$

### Mixed Diffuser Flow Relations

If the flow is supersonic after mixing, meaning
*M _{M}* > 1, a shock wave occurs. The
pressure rise after the shock is

$$\frac{{p}_{ms}}{{p}_{m}}=1+\frac{2\gamma}{\gamma +1}\left({M}_{m}^{2}-1\right),$$

where *M _{m}* is the Mach
number after the primary and secondary flow mix.

The Mach number after the shock is

$${\text{M}}_{ms}^{2}=\frac{\frac{\gamma -1}{\gamma +1}\left({M}_{m}^{2}-1\right)+1}{\frac{2\gamma}{\gamma +1}\left({M}_{m}^{2}-1\right)+1}.$$

After the shock, the diffuser decelerates the flow to raise the pressure. The block assumes
that the kinetic energy at the diffuser outlet is negligible and decelerates the
flow to stagnation pressure at *p _{d}*:

$$\frac{{p}_{d}}{{p}_{ms}}={\left(1+\frac{\gamma -1}{2}{M}_{ms}^{2}\right)}^{\frac{\gamma}{\gamma -1}}.$$

If *M _{M}* ≤ 1, there is no shock and

$$\begin{array}{l}\frac{{p}_{ms}}{{p}_{m}}=1\\ {M}_{ms}^{2}={M}_{m}^{2}.\end{array}$$

The relationship between the mixing pressure, *p _{m}*,
and the diffuser outlet pressure,

*p*, is

_{d}$$\frac{{p}_{d}}{{p}_{m}}=\{\begin{array}{cc}\frac{{\left[\frac{\gamma +1}{2}{M}_{m}^{2}\right]}^{\frac{\gamma}{\gamma -1}}}{{\left[\frac{2\gamma}{\gamma +1}\left({M}_{m}^{2}-1\right)+1\right]}^{\frac{1}{\gamma -1}}},& {M}_{m}>1\\ {\left(1+\frac{\gamma -1}{2}{M}_{m}^{2}\right)}^{\frac{\gamma}{\gamma -1}},& {M}_{m}\le 1\end{array}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}$$

### Critical Mode Operation

In normal operation, the primary flow is high pressure and is choked during operation. The block can also operate if the secondary flow is choked. If both the primary and secondary flows are choked, the block operates in critical mode.

The critical mixing pressure,
*p _{m}^{*}*, is
a low pressure threshold that causes choking in both primary and secondary flow:

$${\text{p}}_{m}^{*}=\mathrm{min}\left\{{p}_{p},{p}_{s}\right\}{\left(\frac{2}{\gamma +1}\right)}^{\frac{\gamma}{\gamma -1}}.$$

The block uses the same equations as above, with
*p _{m}* replaced by

*p*, to calculate the critical diffuser outlet pressure,

_{m}^{*}*p*. If the actual outlet pressure is greater than the critical diffuser outlet pressure, then the ejector is subcritical and the block relates the mixing pressure to diffuser outlet pressure. If the actual outlet pressure is less than the critical diffuser outlet pressure, then the ejector is critical and the block limits the mixing pressure to the critical mixing pressure. This means that:

_{d}^{*}$${\text{p}}_{m}=\{\begin{array}{cc}\frac{{p}_{d}}{\left(\frac{{p}_{d}}{{p}_{m}}\right)}& {p}_{d}\ge {p}_{d}^{*}\\ {p}_{m}^{*}& {p}_{d}<{p}_{d}^{*}\end{array}$$

### Assumptions and Limitations

The block uses isentropic ideal gas relations to derive the model equations.

The block accounts for losses due to friction, mixing, and expansion waves and shocks at the nozzle exit using empirical loss coefficients.

After the primary flow exits the nozzle, the flows do not mix until the primary flow has fully expanded.

After expansion, the flows mix under uniform pressure.

Kinetic energy at the primary flow inlet, secondary flow inlet, and diffuser outlet are negligible compared to the kinetic energy of the flow inside the ejector.

The flow is steady and one-dimensional.

The flow is adiabatic.

Model results with reversed flow may not be accurate.

## Ports

### Conserving

## Parameters

## References

[1] Huang, B. J., et al. "A 1-D
analysis of ejector performance." *International journal of
refrigeration* 22.5 (1999): 354-364.

[2] Chen, WeiXiong, et al. "A 1D
model to predict ejector performance at critical and sub-critical operational regimes."
*International journal of refrigeration* 36.6
(2013): 1750-1761.

## Extended Capabilities

## Version History

**Introduced in R2023a**