modelWetCoilEffectivenessNTU

Heat exchanger with effectiveness - NTU relation and with moisture condensation

Extends from Buildings.Fluid.Interfaces.PartialFourPortInterface, Buildings.Fluid.Interfaces.FourPortFlowResistanceParameters.

Information

This model describes a cooling coil applicable for fully-dry, partially-wet, and fully-wet regimes. The model is developed for counter flow heat exchangers but is also applicable for the cross-flow configuration, although in the latter case it is recommended to have more than four tube rows (Elmahdy and Mitalas, 1977 and Braun, 1988). The model can also be used for a heat exchanger which acts as both heating coil (for some period of time) and cooling coil (for the others). However, it is not recommended to use this model for heating coil only or for cooling coil with no water condensation because for these situations, Buildings.Fluid.HeatExchangers.DryCoilEffectivenessNTU computes faster.

Main equations

The coil model consists of two-equation sets, one for the fully-dry mode and the other for the fully-wet mode. For the fully-dry mode, the ε-NTU approach (Elmahdy and Mitalas, 1977) is used. For the fully-wet mode, equations from Braun (1988) and Mitchell and Braun (2012a and b), which are essentially the extension of the ε-NTU approach to simultaneous sensible and latent heat transfer, are utilized. The equation sets are switched depending on the switching criteria described below that determines the right mode based on a coil surface temperature and dew-point temperature for the air at the inlet of the coil. The transition regime between the two modes, which represents the partially-wet and partially-dry coil, is approximated by employing a fuzzy modeling approach, so-called Takagi-Sugeno fuzzy modeling (Takagi and Sugeno, 1985), which provides a continuously differentiable model that can cover all fully-dry, partially-wet, and fully-wet regimes.

The switching rules are:

  • R1: If the coil surface temperature at the air inlet is lower than the dew-point temperature of air at inlet, then the cooling coil surface is fully-wet.
  • R2: If the coil surface temperature at the air outlet is higher than the dew-point temperature of air at inlet, then the cooling coil surface is fully-dry.
  • R3: If any of the conditions in R1 or R2 is not satisfied, then the cooling coil surface is partially wet.

For more detailed descriptions of the fully-wet coil model and the fuzzy modeling approach, see Buildings.Fluid.HeatExchangers.BaseClasses.WetCoilWetRegime. and Buildings.Fluid.HeatExchangers.BaseClasses.WetCoilDryWetRegime.

Assumptions and limitations

This model contains the following assumptions and limitations:

Medium 2 must be air due to the use of various psychrometric functions.

When parameterizing this model with rated conditions (with the parameter use_UA_nominal set to false), those should correspond to a fully-dry or a fully-wet coil regime, because the model uncertainty yielded by partially-wet rated conditions has not been assessed yet.

The model uses steady-state physics. That is, no dynamics associated with water and coil materials are considered.

The Lewis number, which relates the mass transfer coefficient to the heat transfer coefficient, is assumed to be 1.

The model is not suitable for a cross-flow heat exchanger of which the number of passes is less than four.

By default, the flow regime, such as counter flow or parallel flow, is kept constant based on the parameter value configuration. If a flow reverses direction, it is not changed, e.g., a heat exchanger does not change from counter flow to parallel flow if one flow changes direction. To dynamically change the flow regime, set the constant use_dynamicFlowRegime to true. However, use_dynamicFlowRegime=true can cause slower simulation due to events.

Validation

Validation results can be found in Buildings.Fluid.HeatExchangers.Validation.WetCoilEffectivenessNTU.

References

Braun, James E. 1988. "Methodologies for the Design and Control of Central Cooling Plants". PhD Thesis. University of Wisconsin - Madison. Available online.

Mitchell, John W., and James E. Braun. 2012a. Principles of heating, ventilation, and air conditioning in buildings. Hoboken, N.J.: Wiley.

Mitchell, John W., and James E. Braun. 2012b. "Supplementary Material Chapter 2: Heat Exchangers for Cooling Applications". Excerpt from Principles of heating, ventilation, and air conditioning in buildings. Hoboken, N.J.: Wiley. Available online.

Elmahdy, A.H. and Mitalas, G.P. 1977. "A Simple Model for Cooling and Dehumidifying Coils for Use In Calculating Energy Requirements for Buildings". ASHRAE Transactions. Vol.83. Part 2. pp. 103-117.

Takagi, T. and Sugeno, M., 1985. Fuzzy identification of systems and its applications to modeling and control.  IEEE transactions on systems, man, and cybernetics, (1), pp.116-132.

Parameters

TypeNameDefaultDescription
Booleanuse_dynamicFlowRegimefalseIf true, flow regime is determined using actual flow rates
Buildings.Fluid.Types.HeatExchangerConfigurationconfigurationBuildings.Fluid.Types.HeatExchangerConfiguration.CounterFlowHeat exchanger configuration
Realr_nominal2/3Ratio between air-side and water-side convective heat transfer coefficient
Nominal thermal performance
Booleanuse_Q_flow_nominalfalseSet to true to specify Q_flow_nominal and inlet conditions, or to false to specify UA_nominal
Modelica.Units.SI.HeatFlowRateQ_flow_nominalNominal heat flow rate (positive for heat transfer from 1 to 2)
Modelica.Units.SI.TemperatureT_a1_nominalWater inlet temperature at a rated condition
Modelica.Units.SI.TemperatureT_a2_nominalAir inlet temperature at a rated condition
Modelica.Units.SI.MassFractionw_a2_nominalHumidity ratio of inlet air at a rated condition (in kg/kg dry air)
Modelica.Units.SI.ThermalConductanceUA_nominalThermal conductance at nominal flow, used to compute heat capacity
Dynamics › Conservation equations
Modelica.Fluid.Types.DynamicsenergyDynamicsModelica.Fluid.Types.Dynamics.SteadyStateType of energy balance: dynamic (3 initialization options) or steady state

Components

TypeNameDefaultDescription
Modelica.Units.SI.HeatFlowRateQ1_flow-dryWetCalcs.QTot_flowHeat input into water stream (positive if air is cooled)
Modelica.Units.SI.HeatFlowRateQ2_flowdryWetCalcs.QTot_flowTotal heat input into air stream (negative if air is cooled)
Modelica.Units.SI.HeatFlowRateQSen2_flowdryWetCalcs.QSen_flowSensible heat input into air stream (negative if air is cooled)
Modelica.Units.SI.HeatFlowRateQLat2_flowBuildings.Utilities.Psychrometrics.Constants.h_fg*mWat_flowLatent heat input into air (negative if air is dehumidified)
RealSHRQSen2_flow/noEvent(if (Q2_flow > 1E-6 or Q2_flow < -1E-6) then Q2_flow else 1)Sensible to total heat ratio
Modelica.Units.SI.MassFlowRatemWat_flowdryWetCalcs.mCon_flowWater flow rate of condensate removed from the air stream
RealdryFradryWetCalcs.dryFraDry fraction, 0.3 means condensation occurs at 30% heat exchange length from air inlet

Revisions

  • March 31, 2026, by Michael Wetter:
    Corrected unit propagation error that causes Dymola 2026x to not show certain units.
    See #2100.
  • February 7, 2025, by Jelger Jansen:
    Removed import statement. This is for IBPSA, #1961.
  • February 3, 2023, by Jianjun Hu:
    Added noEvent() in the assertion function to avoid Optimica to not converge.
    This is for issue 1690.
  • January 24, 2023, by Hongxiang Fu:
    Set flowRegime to be equal to flowRegime_nominal by default. Added an assertion warning to inform the user about how to change this behaviour if the flow direction does need to change.
    This is for issue 1682.
  • March 3, 2022, by Michael Wetter:
    Removed massDynamics.
    This is for issue 1542.
  • November 2, 2021, by Michael Wetter:
    Corrected unit assignment during the model instantiation.
    This is for issue 2710.
  • Jan 21, 2021, by Donghun Kim:
    First implementation.