# System-Level Condenser Evaporator (2P-TL)

Heat exchanger between two-phase fluid and thermal liquid networks, with model based on performance data

**Library:**Simscape / Fluids / Fluid Network Interfaces / Heat Exchangers

## Description

The System-Level Condenser Evaporator (2P-TL) block models a heat exchanger between a two-phase fluid network and a thermal liquid network.

The block can act as a condenser or as an evaporator in a refrigeration system, depending on the direction of heat transfer. The block model is based on performance data from the heat exchanger datasheet, rather than on the detailed geometry of the exchanger, and therefore lets you easily adjust the size and performance of the heat exchanger during design iterations, or model heat exchangers with uncommon geometries. You can also use this block to model heat exchangers with a certain level of performance at an early design stage, when detailed geometry data is not yet available.

You parameterize the block by the nominal operating condition. The heat exchanger is sized to match the specified performance at the nominal operating condition at steady state.

The Two-Phase Fluid 1 side approximates the liquid zone, mixture zone, and vapor zone based on the change in enthalpy along the flow path.

This block is similar to the Condenser Evaporator (TL-2P) block but uses a different parameterization model. The table provides a comparison of the two blocks, to help you choose the right block for your application.

Condenser Evaporator (TL-2P) | System-Level Condenser Evaporator (2P-TL) |
---|---|

Block parameters are based on the heat exchanger geometry | Block parameters are based on performance and operating conditions |

Heat exchanger geometry may be limited by the available geometry parameter options | Model is independent of the specific heat exchanger geometry |

You can adjust the block for different performance requirements by tuning geometry parameters, such as fin sizes and tube lengths | You can adjust the block for different performance requirements by directly specifying the desired heat and mass flow rates |

Lets you select between parallel, counter, or cross flow configurations | Lets you select between parallel, counter, or cross flow arrangement at nominal operating conditions, to help with sizing |

Predictively accurate results over a wide range of operating conditions, subject to the applicability of the E-NTU equations and the heat transfer coefficient correlations | Very accurate results around the specified operating condition; accuracy may decrease far away from the specified operating conditions |

Heat transfer calculations account for the variation of temperature along the flow path by using the E-NTU model | Heat transfer calculations approximate the variation of temperature along the flow path by dividing it into three segments |

Accounts for different fluid properties and heat transfer coefficients for subcooled liquid, liquid-vapor mixture, and superheated vapor | Accounts for different fluid properties and heat transfer coefficients for subcooled liquid, liquid-vapor mixture, and superheated vapor |

Keeps track of variable zone length fractions for subcooled liquid, liquid-vapor mixture, and superheated vapor regions based on the geometry | Approximates the effect of subcooled liquid, liquid-vapor mixture, and superheated vapor regions using weighting factors based on the difference in enthalpy between inlet and outlet |

Does not model the wall thermal mass; you can approximate the effect by connecting a pipe block with a thermal mass downstream | Includes an option to model the wall thermal mass |

### Heat Transfer

The two-phase fluid flow and the thermal liquid flow are each divided into three segments of equal size. Heat transfer between the fluids is calculated in each segment. For simplicity, the equation for one segment is shown here.

If the wall thermal mass is off, then the heat balance in the heat exchanger is

$${Q}_{seg,2P}+{Q}_{seg,TL}=0,$$

where:

*Q*_{seg,2P}is the heat flow rate from the wall (that is, the heat transfer surface) to the two-phase fluid in the segment.*Q*_{seg,TL}is the heat flow rate from the wall to the thermal liquid in the segment.

If the wall thermal mass is on, then the heat balance in the heat exchanger is

$${Q}_{seg,2P}+{Q}_{seg,TL}=-\frac{{M}_{wall}{c}_{{p}_{wall}}}{N}\frac{d{T}_{seg,wall}}{dt},$$

where:

*M*_{wall}is the mass of the wall.*c*_{pwall}is the specific heat of the wall.*N*= 3 is the number of segments.*T*_{seg,wall}is the average wall temperature in the segment.*t*is time.

The heat flow rate from the wall to the two-phase fluid in the segment is

$${Q}_{seg,2P}=U{A}_{seg,2P}\left({T}_{seg,wall}-{T}_{seg,2P}\right),$$

where:

*UA*_{seg,2P}is the weighted-average heat transfer conductance for the two-phase fluid in the segment.*T*_{seg,2P}is the weighted-average fluid temperature for the two-phase fluid in the segment.

The heat flow rate from the wall to the thermal liquid in the segment is

$${Q}_{seg,TL}=U{A}_{seg,TL}\left({T}_{seg,wall}-{T}_{seg,TL}\right),$$

where:

*UA*_{seg,TL}is the heat transfer conductance for the thermal liquid in the segment.*T*_{seg,TL}is the average liquid temperature in the segment.

### Two-Phase Fluid Heat Transfer Correlation

If the segment is subcooled liquid, then the heat transfer conductance is

$$U{A}_{seg,L,2P}={a}_{L,2P}{\left({\mathrm{Re}}_{seg,L,2P}\right)}^{{b}_{2P}}{\left({\mathrm{Pr}}_{seg,L,2P}\right)}^{{c}_{2P}}{k}_{seg,L,2P}\frac{{G}_{2P}}{N},$$

where:

*a*_{L,2P},*b*_{2P}, and*c*_{2P}are the coefficients of the Nusselt number correlation. These coefficients appear as block parameters in the**Correlation Coefficients**section.*Re*_{seg,L,2P}is the average liquid Reynolds number for the segment.*Pr*_{seg,L,2P}is the average liquid Prandtl number for the segment.*k*_{seg,L,2P}is the average liquid thermal conductivity for the segment.*G*_{2P}is the geometry scale factor for the two-phase fluid side of the heat exchanger. The block calculates the geometry scale factor so that the total heat transfer over all segments matches the specified performance at the nominal operating conditions.

The average liquid Reynolds number is

$${\mathrm{Re}}_{seg,L,2P}=\frac{{\dot{m}}_{seg,2P}{D}_{ref,2P}}{{\mu}_{seg,L,2P}{S}_{ref,2P}},$$

where:

$${\dot{m}}_{seg,2P}$$ is the mass flow rate through the segment.

*μ*_{seg,L,2P}is the average liquid dynamic viscosity for the segment.*D*_{ref,2P}is an arbitrary reference diameter.*S*_{ref,2P}is an arbitrary reference flow area.

**Note**

The *D*_{ref,2P} and
*S*_{ref,2P} terms are included in this equation
for unit calculation purposes only, to make
*Re*_{seg,L,2P} nondimensional. The values of
*D*_{ref,2P} and
*S*_{ref,2P} are arbitrary because the
*G*_{2P} calculation overrides these values.

Similarly, if the segment is superheated vapor, then the heat transfer conductance is

$$U{A}_{seg,V,2P}={a}_{V,2P}{\left({\mathrm{Re}}_{seg,V,2P}\right)}^{{b}_{2P}}{\left({\mathrm{Pr}}_{seg,V,2P}\right)}^{{c}_{2P}}{k}_{seg,V,2P}\frac{{G}_{2P}}{N},$$

where:

*a*_{V,2P},*b*_{2P}, and*c*_{2P}are the coefficients of the Nusselt number correlation. These coefficients appear as block parameters in the**Correlation Coefficients**section.*Re*_{seg,V,2P}is the average vapor Reynolds number for the segment.*Pr*_{seg,V,2P}is the average vapor Prandtl number for the segment.*k*_{seg,V,2P}is the average vapor thermal conductivity for the segment.

The average vapor Reynolds number is

$${\mathrm{Re}}_{seg,V,2P}=\frac{{\dot{m}}_{seg,2P}{D}_{ref,2P}}{{\mu}_{seg,V,2P}{S}_{ref,2P}},$$

where *μ*_{seg,V,2P} is the average vapor dynamic
viscosity for the segment.

If the segment is liquid-vapor mixture, then the heat transfer conductance is

$$U{A}_{seg,M,2P}={a}_{M,2P}{\left({\mathrm{Re}}_{seg,SL,2P}\right)}^{{b}_{2P}}CZ{\left({\mathrm{Pr}}_{seg,SL,2P}\right)}^{{c}_{2P}}{k}_{seg,SL,2P}\frac{{G}_{2P}}{N},$$

where:

*a*_{M,2P},*b*_{2P}, and*c*_{2P}are the coefficients of the Nusselt number correlation. These coefficients appear as block parameters in the**Correlation Coefficients**section.*Re*_{seg,SL,2P}is the saturated liquid Reynolds number for the segment.*Pr*_{seg,SL,2P}is the saturated liquid Prandtl number for the segment.*k*_{seg,SL,2P}is the saturated liquid thermal conductivity for the segment.*CZ*is the Cavallini and Zecchin term.

The saturated liquid Reynolds number is

$${\mathrm{Re}}_{seg,SL,2P}=\frac{{\dot{m}}_{seg,2P}{D}_{ref,2P}}{{\mu}_{seg,SL,2P}{S}_{ref,2P}},$$

where *μ*_{seg,SL,2P} is the saturated liquid
dynamic viscosity for the segment.

The Cavallini and Zecchin term is

$$CZ=\frac{{\left(\left(\sqrt{\frac{{\nu}_{seg,SV,2P}}{{\nu}_{seg,SL,2P}}}-1\right)\left({x}_{seg,out,2P}+1\right)\right)}^{1+{b}_{2P}}-{\left(\left(\sqrt{\frac{{\nu}_{seg,SV,2P}}{{\nu}_{seg,SL,2P}}}-1\right)\left({x}_{seg,in,2P}+1\right)\right)}^{1+{b}_{2P}}}{\left(1+{b}_{2P}\right)\left(\sqrt{\frac{{\nu}_{seg,SV,2P}}{{\nu}_{seg,SL,2P}}}-1\right)\left({x}_{seg,out,2P}-{x}_{seg,in,2P}\right)},$$

where:

*ν*_{seg,SL,2P}is the saturated liquid specific volume for the segment.*ν*_{seg,SV,2P}is the saturated vapor specific volume for the segment.*x*_{seg,in,2P}is the vapor quality at the segment inlet.*x*_{seg,out,2P}is the vapor quality at the segment outlet.

The expression is based on the work of Cavallini and Zecchin [5], which derives a heat
transfer coefficient correlation at a local vapor quality *x*. Equations
for the liquid-vapor mixture are obtained by averaging Cavallini and Zecchin’s correlation
over the segment from *x*_{seg,in,2P} to
*x*_{seg,out,2P}.

### Two-Phase Fluid Weighted Average

The two-phase fluid flow through a segment may not be entirely represented as either
subcooled liquid, superheated vapor, or liquid-vapor mixture. Instead, each segment may
consist of a combination of these. The block approximates this condition by computing
weighting factors (*w*_{L},
*w*_{V}, and
*w*_{M}) based on the change in specific enthalpy
across the segment and the saturated liquid and vapor specific enthalpies. The block assumes
that the specific enthalpy across the segment varies piecewise linearly from inlet to
outlet, with the breakpoints corresponding to the saturation boundaries for liquid and
vapor. The zone with a larger heat transfer coefficient has a steeper slope than the zone
with a lower heat transfer coefficient.

$$\begin{array}{l}{w}_{L}=\frac{{\Delta}_{L}}{{\Delta}_{L}+{\Delta}_{M}+{\Delta}_{V}}\\ {w}_{V}=\frac{{\Delta}_{V}}{{\Delta}_{L}+{\Delta}_{M}+{\Delta}_{V}}\\ {w}_{M}=1-{w}_{L}-{w}_{V}\end{array}$$

$$\begin{array}{l}{\Delta}_{L}=\left|\mathrm{min}\left({h}_{seg,out,2P},{h}_{seg,SL,2P}\right)-\mathrm{min}\left({h}_{seg,in,2P},{h}_{seg,SL,2P}\right)\right|\cdot U{A}_{seg,M,2P}\cdot U{A}_{seg,V,2P}\\ {\Delta}_{M}=\left|\mathrm{min}\left(\mathrm{max}\left({h}_{seg,out,2P},{h}_{seg,SL,2P}\right),{h}_{seg,SV,2P}\right)-\mathrm{min}\left(\mathrm{max}\left({h}_{seg,in,2P},{h}_{seg,SL,2P}\right),{h}_{seg,SV,2P}\right)\right|\cdot U{A}_{seg,L,2P}\cdot U{A}_{seg,V,2P}\\ {\Delta}_{V}=\left|\mathrm{max}\left({h}_{seg,out,2P},{h}_{seg,SV,2P}\right)-\mathrm{max}\left({h}_{seg,in,2P},{h}_{seg,SV,2P}\right)\right|\cdot U{A}_{seg,L,2P}\cdot U{A}_{seg,M,2P}\end{array}$$

where:

*h*_{seg,in,2P}is the specific enthalpy at the segment inlet.*h*_{seg,out,2P}is the specific enthalpy at the segment outlet.*h*_{seg,SL,2P}is the saturated liquid specific enthalpy for the segment.*h*_{seg,SV,2P}is the saturated vapor specific enthalpy for the segment.

The weighted-average two-phase fluid heat transfer conductance for the segment is therefore

$$U{A}_{seg,2P}={w}_{L}\left(U{A}_{seg,L,2P}\right)+{w}_{V}\left(U{A}_{seg,V,2P}\right)+{w}_{M}\left(U{A}_{seg,M,2P}\right).$$

The weighted-average fluid temperature for the segment is

$${T}_{seg,2P}=\frac{{w}_{L}\left(U{A}_{seg,L,2P}\right){T}_{seg,L,2P}+{w}_{V}\left(U{A}_{seg,V,2P}\right){T}_{seg,V,2P}+{w}_{M}\left(U{A}_{seg,M,2P}\right){T}_{seg,M,2P}}{U{A}_{seg,2P}},$$

where:

*T*_{seg,L,2P}is the average liquid temperature for the segment.*T*_{seg,V,2P}is the average vapor temperature for the segment.*T*_{seg,M,2P}is the average mixture temperature for the segment, which is the saturated liquid temperature.

### Thermal Liquid Heat Transfer Correlation

The heat transfer conductance is

$$U{A}_{seg,TL}={a}_{TL}{\left({\mathrm{Re}}_{seg,TL}\right)}^{{b}_{TL}}{\left({\mathrm{Pr}}_{seg,TL}\right)}^{{c}_{TL}}{k}_{seg,TL}\frac{{G}_{TL}}{N},$$

where:

*a*_{TL},*b*_{TL}, and*c*_{TL}are the coefficients of the Nusselt number correlation. These coefficients appear as block parameters in the**Correlation Coefficients**section.*Re*_{seg,TL}is the average Reynolds number for the segment.*Pr*_{seg,TL}is the average Prandtl number for the segment.*k*_{seg,TL}is the average thermal conductivity for the segment.*G*_{TL}is the geometry scale factor for the thermal liquid side of the heat exchanger. The block calculates the geometry scale factor so that the total heat transfer over all segments matches the specified performance at the nominal operating conditions.

The average Reynolds number is

$${\mathrm{Re}}_{seg,\text{TL}}=\frac{{\dot{m}}_{seg,\text{TL}}{D}_{ref,\text{TL}}}{{\mu}_{seg,\text{TL}}{S}_{ref,\text{TL}}},$$

where:

$${\dot{m}}_{seg,\text{TL}}$$ is the mass flow rate through the segment.

*μ*_{seg,TL}is the average dynamic viscosity for the segment.*D*_{ref,TL}is an arbitrary reference diameter.*S*_{ref,TL}is an arbitrary reference flow area.

**Note**

The *D*_{ref,TL} and
*S*_{ref,TL} terms are included in this equation
for unit calculation purposes only, to make
*Re*_{seg,TL} nondimensional. The values of
*D*_{ref,TL} and
*S*_{ref,TL} are arbitrary because the
*G*_{TL} calculation overrides these values.

### Pressure Loss

The pressure losses on the two-phase fluid side are

$$\begin{array}{l}{p}_{A,2P}-{p}_{2P}=\frac{{K}_{2P}}{2}\frac{{\dot{m}}_{A,2P}\sqrt{{\dot{m}}^{2}{}_{A,2P}+{\dot{m}}^{2}{}_{thres,2P}}}{2{\rho}_{avg,2P}}\\ {p}_{B,2P}-{p}_{2P}=\frac{{K}_{2P}}{2}\frac{{\dot{m}}_{B,2P}\sqrt{{\dot{m}}^{2}{}_{B,2P}+{\dot{m}}^{2}{}_{thres,2P}}}{2{\rho}_{avg,2P}}\end{array}$$

where:

*p*_{A,2P}and*p*_{B,2P}are the pressures at ports**A1**and**B1**, respectively.*p*_{2P}is internal two-phase fluid pressure at which the heat transfer is calculated.$${\dot{m}}_{A,2P}$$ and $${\dot{m}}_{B,2P}$$ are the mass flow rates into ports

**A1**and**B1**, respectively.*ρ*_{avg,2P}is the average two-phase fluid density over all segments.$${\dot{m}}_{thres,2P}$$ is the laminar threshold for pressure loss, approximated as 1e-4 of the nominal mass flow rate. The block calculates the pressure loss coefficient,

*K*_{2P}, so that*p*_{A,2P}–*p*_{B,2P}matches the nominal pressure loss at the nominal mass flow rate.

The pressure losses on the thermal liquid side are

$$\begin{array}{l}{p}_{A,\text{TL}}-{p}_{\text{TL}}=\frac{{K}_{\text{TL}}}{2}\frac{{\dot{m}}_{A,\text{TL}}\sqrt{{\dot{m}}^{2}{}_{A,\text{TL}}+{\dot{m}}^{2}{}_{thres,\text{TL}}}}{2{\rho}_{avg,2P}}\\ {p}_{B,\text{TL}}-{p}_{\text{TL}}=\frac{{K}_{\text{TL}}}{2}\frac{{\dot{m}}_{B,\text{TL}}\sqrt{{\dot{m}}^{2}{}_{B,\text{TL}}+{\dot{m}}^{2}{}_{thres,\text{TL}}}}{2{\rho}_{avg,\text{TL}}}\end{array}$$

where:

*p*_{A,TL}and*p*_{B,TL}are the pressures at ports**A2**and**B2**, respectively.*p*_{TL}is internal thermal liquid pressure at which the heat transfer is calculated.$${\dot{m}}_{A,TL}$$ and $${\dot{m}}_{B,TL}$$ are the mass flow rates into ports

**A2**and**B2**, respectively.*ρ*_{avg,TL}is the average thermal liquid density over all segments.$${\dot{m}}_{thres,TL}$$ is the laminar threshold for pressure loss, approximated as 1e-4 of the nominal mass flow rate. The block calculates the pressure loss coefficient,

*K*_{TL}, so that*p*_{A,TL}–*p*_{B,TL}matches the nominal pressure loss at the nominal mass flow rate.

### Two-Phase Fluid Mass and Energy Conservation

The mass conservation equation for the overall two-phase fluid flow is

$$\left(\frac{d{p}_{2P}}{dt}{\displaystyle \sum _{segments}\left(\frac{\partial {\rho}_{seg,2P}}{\partial p}\right)}+{\displaystyle \sum _{segments}\left(\frac{d{u}_{seg,2P}}{dt}\frac{\partial {\rho}_{seg,2P}}{\partial u}\right)}\right)\frac{{V}_{2P}}{N}={\dot{m}}_{A,2P}+{\dot{m}}_{B,2P},$$

where:

$$\frac{\partial {\rho}_{seg,2P}}{\partial p}$$ is the partial derivative of density with respect to pressure for the segment.

$$\frac{\partial {\rho}_{seg,2P}}{\partial u}$$ is the partial derivative of density with respect to specific internal energy for the segment.

*u*_{seg,2P}is the specific internal energy for the segment.*V*_{2P}is the total two-phase fluid volume.

The summation is over all segments.

**Note**

Although the two-phase fluid flow is divided into *N*=3 segments for
heat transfer calculations, all segments are assumed to be at the same internal pressure,
*p*_{2P}. That is why
*p*_{2P} is outside of the summation.

The energy conservation equation for each segment is

$$\frac{d{u}_{seg,2P}}{dt}\frac{{M}_{2P}}{N}+{u}_{seg,2P}\left({\dot{m}}_{seg,in,2P}-{\dot{m}}_{seg,out,2P}\right)={\Phi}_{seg,in,2P}-{\Phi}_{seg,out,2P}+{Q}_{seg,2P},$$

where:

*M*_{2P}is the total two-phase fluid mass.$${\dot{m}}_{seg,in,2P}$$ and $${\dot{m}}_{seg,out,2P}$$ are the mass flow rates into and out of the segment.

*Φ*_{seg,in,2p}and*Φ*_{seg,out,2p}are the energy flow rates into and out of the segment.

The mass flow rates between segments are assumed to be linearly distributed between the values of$${\dot{m}}_{A,2P}$$ and $${\dot{m}}_{B,2P}$$.

### Thermal Liquid Mass and Energy Conservation

The mass conservation for the overall thermal liquid flow is

$$\left(\frac{d{p}_{TL}}{dt}{\displaystyle \sum _{segments}\left(\frac{\partial {\rho}_{seg,TL}}{\partial p}\right)}+{\displaystyle \sum _{segments}\left(\frac{d{T}_{seg,TL}}{dt}\frac{\partial {\rho}_{seg,TL}}{\partial T}\right)}\right)\frac{{V}_{TL}}{N}={\dot{m}}_{A,TL}+{\dot{m}}_{B,TL},$$

where:

$$\frac{\partial {\rho}_{seg,TL}}{\partial p}$$ is the partial derivative of density with respect to pressure for the segment.

$$\frac{\partial {\rho}_{seg,TL}}{\partial T}$$ is the partial derivative of density with respect to temperature for the segment.

*T*_{seg,TL}is the temperature for the segment.*V*_{TL}is the total thermal liquid volume.

The summation is over all segments.

**Note**

Although the thermal liquid flow is divided into *N*=3 segments for
heat transfer calculations, all segments are assumed to be at the same internal pressure,
*p*_{TL}. That is why
*p*_{TL} is outside of the summation.

The energy conservation equation for each segment is

$$\begin{array}{l}\left(\frac{d{p}_{TL}}{dt}\frac{\partial {u}_{seg,TL}}{\partial p}+\frac{d{T}_{seg,TL}}{dt}\frac{\partial {u}_{seg,TL}}{\partial T}\right)\frac{{M}_{TL}}{N}+{u}_{seg,TL}\left({\dot{m}}_{seg,in,TL}-{\dot{m}}_{seg,out,TL}\right)=\\ {\Phi}_{seg,in,TL}-{\Phi}_{seg,out,TL}+{Q}_{seg,TL},\end{array}$$

where:

$$\frac{\partial {u}_{seg,TL}}{\partial p}$$ is the partial derivative of specific internal energy with respect to pressure for the segment.

$$\frac{\partial {u}_{seg,TL}}{\partial T}$$ is the partial derivative of specific internal energy with respect to temperature for the segment.

*M*_{TL}is the total thermal liquid mass.$${\dot{m}}_{seg,in,TL}$$ and $${\dot{m}}_{seg,out,TL}$$ are the mass flow rates into and out of the segment.

*Φ*_{seg,in,TL}and*Φ*_{seg,out,TL}are the energy flow rates into and out of the segment.

The mass flow rates between segments are assumed to be linearly distributed between the values of$${\dot{m}}_{A,TL}$$ and $${\dot{m}}_{B,TL}$$.

## Ports

### Output

### Conserving

## Parameters

## Extended Capabilities

## Version History

**Introduced in R2022a**