# Temperature Control Valve (G)

Temperature control valve in a gas network

**Libraries:**

Simscape /
Fluids /
Gas /
Valves & Orifices /
Flow Control Valves

## Description

The Temperature Control Valve (G) block represents an orifice with a thermostat as a flow control mechanism. The thermostat contains a temperature sensor and an opening mechanism. The sensor is at the inlet and responds with a slight delay, captured by a first-order time lag, to variations in temperature.

When the sensor reads a temperature in excess of a preset activation value, the
opening mechanism actuates and the valve begins to open or close, depending on the
operation mode specified by the **Valve operation** parameter. The
change in opening area continues up to the limit of the temperature range of the valve,
beyond which the opening area is a constant. Within the temperature range, the opening
area is a linear function of temperature.

The flow can be laminar or turbulent, and it can reach up to sonic speeds. The maximum velocity happens at the throat of the valve where the flow is narrowest and fastest. The flow chokes and the velocity saturates when a drop in downstream pressure can no longer increase the velocity. Choking occurs when the back-pressure ratio reaches the critical value characteristic of the valve. The block does not capture supersonic flow.

### Control Temperature

The temperature reading at the inlet is the control signal for the valve. The more
the temperature reading rises over the activation temperature, the more the opening
area diverges from maximally closed when **Valve operation** is
`Opens above activation temperature`

, or from fully
open when **Valve operation** is ```
Closes above
activation temperature
```

.

The difference between the sensor temperature reading and the activation temperature is the temperature overshoot. The block normalizes this variable against the temperature regulation range of the valve. The fraction of valve opening is

$$\widehat{T}=\frac{{T}_{S}-{T}_{Act}}{\Delta T},$$

where:

*T*is the value of the_{Act}**Activation temperature**parameter.*T*is the sensor temperature reading._{S}When

**Temperature sensing**is`Valve inlet temperature`

,*T*is the upstream temperature of the valve._{S}When

**Temperature sensing**is`Gas sensing port`

,*T*is the temperature of the gas network where it connects to port_{S}**T**.When

**Temperature sensing**is`Thermal sensing port`

,*T*is the temperature of the thermal network where it connects to port_{S}**T**.

*ΔT*is the value of the**Temperature regulation range**parameter.

**Numerical Smoothing**

When the **Smoothing factor** parameter is nonzero, the block
applies numerical smoothing to the fraction of valve opening, $$\widehat{T}$$. Enabling smoothing helps maintain numerical robustness in
your simulation.

For more information, see Numerical Smoothing.

**Sensor Dynamics**

To emulate a real temperature sensor, which can only register a shift in
temperature gradually, the block adds a first-order time lag to the temperature
reading, *T _{S}*. The lag gives the sensor
a transient response to variations in temperature. This expression for

*T*is

_{S}$$\frac{d}{dt}{T}_{\text{S}}=\frac{{T}_{\text{In}}-{T}_{\text{S}}}{\tau},$$

where *T _{In}* is the
actual inlet temperature at the current time step of the simulation and

*τ*is the value of the

**Sensor time constant**parameter. The smaller this parameter is, the faster the sensor responds.

### Valve Parameterizations

The block behavior depends on the **Valve parametrization**
parameter:

`Cv flow coefficient`

— The flow coefficient*C*_{v}determines the block parameterization. The flow coefficient measures the ease with which a gas can flow when driven by a certain pressure differential.`Kv flow coefficient`

— The flow coefficient*K*_{v}, where $${K}_{v}=0.865{C}_{v}$$, determines the block parameterization. The flow coefficient measures the ease with which a gas can flow when driven by a certain pressure differential.`Sonic conductance`

— The sonic conductance of the resistive element at steady state determines the block parameterization. The sonic conductance measures the ease with which a gas can flow when*choked*, which is a condition in which the flow velocity is at the local speed of sound. Choking occurs when the ratio between downstream and upstream pressures reaches a critical value known as the*critical pressure ratio*.`Orifice area`

— The size of the flow restriction determines the block parametrization.

The block scales the specified flow capacity by the fraction of valve opening. As
the fraction of valve opening rises from `0`

to
`1`

, the measure of flow capacity scales from its specified
minimum to its specified maximum.

### Momentum Balance

The block equations depend on the **Orifice parametrization**
parameter. When you set **Orifice parametrization** to
`Cv flow coefficient parameterization`

, the mass
flow rate, $$\dot{m}$$, is

$$\dot{m}={C}_{v}{N}_{6}Y\sqrt{({p}_{in}-{p}_{out}){\rho}_{in}},$$

where:

*C*is the flow coefficient._{v}*N*is a constant equal to 27.3 for mass flow rate in kg/hr, pressure in bar, and density in kg/m_{6}^{3}.*Y*is the expansion factor.*p*is the inlet pressure._{in}*p*is the outlet pressure._{out}*ρ*is the inlet density._{in}

The expansion factor is

$$Y=1-\frac{{p}_{in}-{p}_{out}}{3{p}_{in}{F}_{\gamma}{x}_{T}},$$

where:

*F*is the ratio of the isentropic exponent to 1.4._{γ}*x*is the value of the_{T}**xT pressure differential ratio factor at choked flow**parameter.

The block smoothly transitions to a linearized form of the equation when the
pressure ratio, $${p}_{out}/{p}_{in}$$, rises above the value of the **Laminar flow pressure
ratio** parameter,
*B _{lam}*,

$$\dot{m}={C}_{v}{N}_{6}{Y}_{lam}\sqrt{\frac{{\rho}_{avg}}{{p}_{avg}(1-{B}_{lam})}}({p}_{in}-{p}_{out}),$$

where:

$${Y}_{lam}=1-\frac{1-{B}_{lam}}{3{F}_{\gamma}{x}_{T}}.$$

When the pressure ratio, $${p}_{out}/{p}_{in}$$, falls below $$1-{F}_{\gamma}{x}_{T}$$, the orifice becomes choked and the block switches to the equation

$$\dot{m}=\frac{2}{3}{C}_{v}{N}_{6}\sqrt{{F}_{\gamma}{x}_{T}{p}_{in}{\rho}_{in}}.$$

When you set **Orifice parametrization** to ```
Kv
flow coefficient parameterization
```

, the block uses these same
equations, but replaces *C _{v}* with

*K*by using the relation $${K}_{v}=0.865{C}_{v}$$. For more information on the mass flow equations when the

_{v}**Orifice parametrization**parameter is

```
Kv
flow coefficient parameterization
```

or ```
Cv flow
coefficient parameterization
```

, see [2][3].When you set **Orifice parametrization** to
`Sonic conductance parameterization`

, the mass flow
rate, $$\dot{m}$$, is

$$\dot{m}=C{\rho}_{ref}{p}_{in}\sqrt{\frac{{T}_{ref}}{{T}_{in}}}{\left[1-{\left(\frac{\frac{{p}_{out}}{{p}_{in}}-{B}_{crit}}{1-{B}_{crit}}\right)}^{2}\right]}^{m},$$

where:

*C*is the sonic conductance.*B*is the critical pressure ratio._{crit}*m*is the value of the**Subsonic index**parameter.*T*is the value of the_{ref}**ISO reference temperature**parameter.*ρ*is the value of the_{ref}**ISO reference density**parameter.*T*is the inlet temperature._{in}

The block smoothly transitions to a linearized form of the equation when the
pressure ratio, $${p}_{out}/{p}_{in}$$, rises above the value of the **Laminar flow pressure
ratio** parameter
*B _{lam}*,

$$\dot{m}=C{\rho}_{ref}\sqrt{\frac{{T}_{ref}}{{T}_{avg}}}{\left[1-{\left(\frac{{B}_{lam}-{B}_{crit}}{1-{B}_{crit}}\right)}^{2}\right]}^{m}\left(\frac{{p}_{in}-{p}_{out}}{1-{B}_{lam}}\right).$$

When the pressure ratio, $${p}_{out}/{p}_{in}$$, falls below the critical pressure ratio,
*B _{crit}*, the orifice becomes
choked and the block switches to the equation

$$\dot{m}=C{\rho}_{ref}{p}_{in}\sqrt{\frac{{T}_{ref}}{{T}_{in}}}.$$

For more information on the mass flow equations when the **Orifice
parametrization** parameter is ```
Sonic conductance
parameterization
```

, see [1].

When you set **Orifice parametrization** to
`Orifice area parameterization`

, the mass flow
rate, $$\dot{m}$$, is

$$\dot{m}={C}_{d}{S}_{r}\sqrt{\frac{2\gamma}{\gamma -1}{p}_{in}{\rho}_{in}{\left(\frac{{p}_{out}}{{p}_{in}}\right)}^{\frac{2}{\gamma}}\left[\frac{1-{\left(\frac{{p}_{out}}{{p}_{in}}\right)}^{\frac{\gamma -1}{\gamma}}}{1-{\left(\frac{{S}_{R}}{S}\right)}^{2}{\left(\frac{{p}_{out}}{{p}_{in}}\right)}^{\frac{2}{\gamma}}}\right]},$$

where:

*S*is the orifice or valve area._{r}*S*is the value of the**Cross-sectional area at ports A and B**parameter.*C*is the value of the_{d}**Discharge coefficient**parameter.*γ*is the isentropic exponent.

The block smoothly transitions to a linearized form of the equation when the
pressure ratio, $${p}_{out}/{p}_{in}$$, rises above the value of the **Laminar flow pressure
ratio** parameter,
*B _{lam}*,

$$\dot{m}={C}_{d}{S}_{r}\sqrt{\frac{2\gamma}{\gamma -1}{p}_{avg}^{\frac{2-\gamma}{\gamma}}{\rho}_{avg}{B}_{lam}^{\frac{2}{\gamma}}\left[\frac{1-\text{\hspace{0.17em}}{B}_{lam}^{\frac{\gamma -1}{\gamma}}}{1-{\left(\frac{{S}_{R}}{S}\right)}^{2}{B}_{lam}^{\frac{2}{\gamma}}}\right]}\left(\frac{{p}_{in}^{\frac{\gamma -1}{\gamma}}-{p}_{out}^{\frac{\gamma -1}{\gamma}}}{1-{B}_{lam}^{\frac{\gamma -1}{\gamma}}}\right).$$

When the pressure ratio, $${p}_{out}/{p}_{in}$$, falls below$${\left(\frac{2}{\gamma +1}\right)}^{\frac{\gamma}{\gamma -1}}$$ , the orifice becomes choked and the block switches to the equation

$$\dot{m}={C}_{d}{S}_{R}\sqrt{\frac{2\gamma}{\gamma +1}{p}_{in}{\rho}_{in}\frac{1}{{\left(\frac{\gamma +1}{2}\right)}^{\frac{2}{\gamma -1}}-{\left(\frac{{S}_{R}}{S}\right)}^{2}}}.$$

For more information on the mass flow equations when the **Orifice
parametrization** parameter is ```
Orifice area
parameterization
```

, see [4].

### Mass Balance

The block assumes the volume and mass of fluid inside the valve is very small and ignores these values. As a result, no amount of fluid can accumulate in the valve. By the principle of conservation of mass, the mass flow rate into the valve through one port equals that out of the valve through the other port

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

where $$\dot{m}$$ is defined as the mass flow rate into the valve through the port
indicated by the **A** or **B** subscript.

### Energy Balance

The resistive element of the block is an adiabatic component. No heat exchange can
occur between the fluid and the wall that surrounds it. No work is done on or by the
fluid as it traverses from inlet to outlet. Energy can flow only by advection,
through ports **A** and **B**. By the principle of conservation of energy, the sum of the port
energy flows is always equal to zero

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

where *ϕ* is the energy flow rate into the valve through ports
**A** or **B**.

### 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.

Nominal values provide a way to specify the expected magnitude of a variable in a model.
Using system scaling based on nominal values increases the simulation robustness. Nominal
values can come from different sources, one of which is the **Nominal
Values** section in the block dialog box or Property Inspector. For more
information, see Modify Nominal Values for a Block Variable.

### Assumptions and Limitations

The

`Sonic conductance`

setting of the**Valve parameterization**parameter is for pneumatic applications. If you use this setting for gases other than air, you may need to scale the sonic conductance by the square root of the specific gravity.The equation for the

`Orifice area`

parameterization is less accurate for gases that are far from ideal.This block does not model supersonic flow.

## Ports

### Conserving

## Parameters

## References

[1] ISO 6358-3. "Pneumatic fluid power – Determination of flow-rate characteristics of components using compressible fluids – Part 3: Method for calculating steady-state flow rate characteristics of systems". 2014.

[2] IEC 60534-2-3. "Industrial-process control valves – Part 2-3: Flow capacity – Test procedures". 2015.

[3] ANSI/ISA-75.01.01. "Industrial-Process Control Valves – Part 2-1: Flow capacity – Sizing equations for fluid flow underinstalled conditions". 2012.

[4] P. Beater. *Pneumatic
Drives*. Springer-Verlag Berlin Heidelberg. 2007.

## Extended Capabilities

## Version History

**Introduced in R2018b**