# Local Restriction (2P)

Fixed flow resistance

## Library

Two-Phase Fluid/Elements

## Description

The Local Restriction (2P) block models the pressure drop due to a fixed flow resistance such as an orifice. Ports A and B represent the restriction inlet and outlet. The restriction area, specified in the block dialog box, remains constant during simulation.

The restriction consists of a contraction followed by a sudden expansion in flow area. The contraction causes the fluid to accelerate and its pressure to drop. The expansion recovers the lost pressure though only in part, as the flow separates from the wall, losing momentum in the process.

**Local Restriction Schematic**

### Mass Balance

The mass balance equation is

$${\dot{m}}_{A}+{\dot{m}}_{B}=0,$$

where:

$${\dot{m}}_{A}$$ and $${\dot{m}}_{B}$$ are the mass flow rates into the restriction through port A and port B.

### Energy Balance

The energy balance equation is

$${\varphi}_{A}+{\varphi}_{B}=0,$$

where:

*ϕ*_{A}and*ϕ*_{B}are the energy flow rates into the restriction through port A and port B.

The local restriction is assumed to be adiabatic and the change in specific total enthalpy is therefore zero. At port A,

$${u}_{A}+{p}_{A}{\nu}_{A}+\frac{{w}_{A}^{2}}{2}={u}_{R}+{p}_{R}{\nu}_{R}+\frac{{w}_{R}^{2}}{2},$$

while at port B,

$${u}_{B}+{p}_{B}{\nu}_{B}+\frac{{w}_{B}^{2}}{2}={u}_{R}+{p}_{R}{\nu}_{R}+\frac{{w}_{R}^{2}}{2},$$

where:

*u*_{A},*u*_{B}, and*u*_{R}are the specific internal energies at port A, at port B, and the restriction aperture.*p*_{A},*p*_{B}, and*p*_{R}are the pressures at port A, port B, and the restriction aperture.*ν*_{A},*ν*_{B}, and*ν*_{R}are the specific volumes at port A, port B, and the restriction aperture.*w*_{A},*w*_{B}, and*w*_{R}are the ideal flow velocities at port A, port B, and the restriction aperture.

The ideal flow velocity is computed as

$${w}_{A}=\frac{{\dot{m}}_{ideal}{\nu}_{A}}{S}$$

at port A, as

$${w}_{B}=\frac{{\dot{m}}_{ideal}{\nu}_{B}}{S}$$

at port B, and as

$${w}_{R}=\frac{{\dot{m}}_{ideal}{\nu}_{R}}{{S}_{R}},$$

inside the restriction, where:

$${\dot{m}}_{ideal}$$ is the ideal mass flow rate through the restriction.

*S*is the flow area at port A and port B.*S*_{R}is the flow area of the restriction aperture.

The ideal mass flow rate through the restriction is computed as:

$${\dot{m}}_{ideal}=\frac{{\dot{m}}_{A}}{{C}_{D}},$$

where:

*C*_{D}is the flow discharge coefficient for the local restriction.

**Local Restriction Variables**

### Momentum Balance

The change in momentum between the ports reflects in the pressure loss across the restriction. That loss depends on the mass flow rate through the restriction, though the exact dependence varies with flow regime. When the flow is turbulent:

$$\dot{m}={S}_{\text{R}}\left({p}_{\text{A}}-{p}_{\text{B}}\right)\sqrt{\frac{2}{\left|{p}_{\text{A}}-{p}_{\text{B}}\right|{\nu}_{\text{R}}{K}_{\text{T}}}},$$

where *K*_{T} is defined as:

$${K}_{\text{T}}=\left(1+\frac{{S}_{\text{R}}}{S}\right)\left(1-\frac{{\nu}_{\text{in}}}{{\nu}_{\text{out}}}\frac{{S}_{\text{R}}}{S}\right)-2\frac{{S}_{\text{R}}}{S}\left(1-\frac{{\nu}_{\text{out}}}{{\nu}_{\text{R}}}\frac{{S}_{\text{R}}}{S}\right),$$

in which the subscript `in`

denotes the inlet
port and the subscript `out`

the outlet port. Which port serves as
the inlet and which serves as the outlet depends on the pressure differential across
the restriction. If pressure is greater at port **A** than at port
**B**, then port **A** is the inlet; if
pressure is greater at port **B**, then port **B**
is the inlet.

When the flow is laminar:

$$\dot{m}={S}_{\text{R}}\left({p}_{\text{A}}-{p}_{\text{B}}\right)\sqrt{\frac{2}{\Delta {p}_{\text{Th}}{\nu}_{\text{R}}{\left(1-\frac{{S}_{\text{R}}}{S}\right)}^{2}},}$$

where *Δp*_{Th} denotes the
threshold pressure drop at which the flow begins to smoothly transition between
laminar and turbulent:

$$\Delta {p}_{\text{Th}}=\left(\frac{{p}_{\text{A}}+{p}_{\text{B}}}{2}\right)\left(1-{B}_{\text{L}}\right),$$

in which *B*_{Lam} is the
**Laminar flow pressure ratio** block parameter. The flow is
laminar if the pressure drop from port **A** to port
**B** is below the threshold value; otherwise, the flow is
turbulent.

The pressure at the restriction area, *p*_{R}
likewise depends on the flow regime. When the flow is turbulent:

$${p}_{\text{R,L}}={p}_{\text{in}}-\frac{{\nu}_{\text{R}}}{2}{\left(\frac{\dot{m}}{{S}_{\text{R}}}\right)}^{2}\left(1+\frac{{S}_{\text{R}}}{S}\right)\left(1-\frac{{\nu}_{\text{in}}}{{\nu}_{\text{R}}}\frac{{S}_{\text{R}}}{S}\right).$$

When the flow is laminar:

$${p}_{\text{R,L}}=\frac{{p}_{\text{A}}+{p}_{\text{B}}}{2}.$$

### Variables

To set the priority and initial target values for the block variables prior to simulation, use
the **Initial Targets** section in the block dialog box or Property
Inspector. For more information, see Set Priority and Initial Target for Block Variables.

### Assumptions and Limitations

The restriction is adiabatic. It does not exchange heat with its surroundings.

## Parameters

**Restriction area**Area normal to the flow path at the restriction aperture—the narrow orifice located between the ports. The default value,

`0.01`

m^2, is the same as the port areas.**Cross-sectional area at ports A and B**Area normal to the flow path at the restriction ports. The ports are assumed to be identical in cross-section. The default value,

`0.01`

m^2, is the same as the restriction aperture area.**Flow discharge coefficient**Ratio of the actual to the theoretical mass flow rate through the restriction. The discharge coefficient is an empirical parameter used to account for non-ideal effects such as those due to restriction geometry. The default value is

`0.64`

.**Laminar flow pressure ratio**Ratio of the outlet to the inlet port pressure at which the flow regime is assumed to switch from laminar to turbulent. The prevailing flow regime determines the equations used in simulation. The pressure drop across the restriction is linear with respect to the mass flow rate if the flow is laminar and quadratic (with respect to the mass flow rate) if the flow is turbulent. The default value is

`0.999`

.

## Ports

A pair of two-phase fluid conserving ports labeled A and B represent the restriction inlet and outlet.

## Extended Capabilities

## Version History

**Introduced in R2015b**