modelAnalyticWetExchanger

Information

WetExchanger

This components allows to model heat and moisture laden transfers between a moist air media and a cooling fluid (commonly water). The mathematical model under this component takes inspiration from the paper: Methodologies d’identification d’économies d’energie : application aux systèmes de climatisation à eau glacee where the theoretical approach under the wet exchanger modelling is summarized here. The module embedes a modelling of a fully dry exchanger (variable with suffix '_1' ex:TA_out_1) and a model of a partially or fully wet (variable with suffix '_2' ex:TA_out_2). Both the fully dry exchange and the wet exchange are running simultaneously during the simulation. The output values are taken from one of the two configurations based on the value of the wet surface. A null or negative value means a fully dry exchange. To avoid chattering between the configurations, the outputs are computed by a polynomial interpolation close to zero value for the wet surface.
For information regarding the dry exchange, please refer to the dry exchanger module. The following description focuses only on the wet exchange.

The physical model under this module uses the following assumptions:

  • The lewis number is constant and equal to one
  • The sensible term for the water vapour within the expression of the mixing specific enthalpy is neglected
  • The conduction resistance within the tube between the moist air and the cooling fluid is neglected regarding the great thermal conductivity of used material. The heat regime is considered as quasistatic
  • The exchanger is adiabatic regarding the environment
  • The convective heat exchange coefficients between the tube and the cooling fluid and the moist air and the condensation film are constant
  • The exchanger is linkened to an counter flow exchanger

To sum up the approach of a wet echange modlling described by the paper and the article above, The link between the mass and heat transfer is made via a Lewis number (Le) equal to one. It is discussed in the paper the relevancy of the assumptions of a Le = 1.

image Lewis number

The founding principles of the representation method were established according to the work of [THRELKELD, 1970]. This method aggregates the phenomena of heat and mass transfer in an enthalpy exchange between air and water represented by a fictitious enthalpy.

Exchange between the moist air and the condensation film

It is assumed that the air immediately against the surface of the condensing film is saturated at the surface temperature of the condensing film. The water vapor and the condensation water being in equilibrium, the temperature of the condensation film is the temperature of the saturated air Tsat. The exchange between the air and the condensing film is expressed as follows:

In the temperature range [15 - 30°C] on which we are working, the term [cp T] represents only about 1% of the total enthalpy and can therefore be neglected. The enthalpy of humid air can then be written as follows:

By introducing the expression of the number of LEWIS, we obtain:

If we consider a number of LEWIS of 1 in the domain considered, the expression becomes:

It can be seen that, for the air / condensation film exchange, the two exchange potentials for temperature gradient and humidity gradient have been replaced by a single enthalpy potential.

Exchange between the condensation film and the cooling fluid

Regarding the exchange between the condensation film and the cooling fluid, the principle of representation is to be reduced to a homogeneous expression with that of the exchange between the air and the condensing film. We then assume that in a small temperature scale, we can represent the enthalpy of saturated air as a linear function of temperature:

If the conduction in the tube and the condensing film is neglected due to the high conductivity of the materials used and of the water regarding the convective phenomena, the heat exchange between the condensing film and the water current expresses as follows:

Using the assumption that cpsat does not vary over the small temperature interval taken into account, water is then represented by a fictitious enthalpy corresponding to the enthalpy that saturated air would have at water temperature:

Considering the expression of the heat flow in the different exchanges, we then conceive of a direct exchange between air and temperature by introducing the expression of two resistors in series, one comprising the convective exchange between the air and the condensation film (Uext), the other comprising the conductive exchange (neglected) and the convective exchange between the tube and the water (Uint).
The global conductance is therefore written, associating transfer conductance and exchange surface:

The heat and mass exchange within the cold battery is then described by the following expression of the power:

The whole of this approach therefore amounts to substituting for the two thermodynamic forces generating the heat flow, a single force derived from the enthalpy of humid air. It then becomes possible to apply the calculation techniques developed for heat exchangers using this unified expression of heat and mass transfer.

Methodology of the heat exchange coefficient calculation

As presented above, the value of the convective heat exchange coefficient hv_A and hv_B are required. However, their value are often not available in the manufacturer's data sheets. The only available data is sometimes the global heat exchange coefficient (K). The choice has been to have as parameter K and hv_A to computed hcv_B from K and hcv_A as follows. If the configuration flow is not counter flow, the coefficient K here has to be the coefficient that would have a counter flow exchanger with the same output properties.

In exchanger dealing with air, it is common to used fins to increase the surface of exchange and thus the exchanged heat. The surface of exchange between the cooling fluid and the tube (piping) and between the pipe and the air is not the same. In that case, the supplier mainly refers the value of K to the largest surface (surface tube / air).
With the above formula, the convective exchange coefficient hcv_B is not related to the exchange surface between the cooling fluid and the tube but to the exchange surface between the tube and the air (the largest surface ). In such way, the knowledge of the exchange surface between the cooling fluid and the tube is not required:

Parameters

TypeNameDefaultDescription
Realphi_out0.8
SI.SpecificEnergyLl2501e3Latent heat of liquifaction at 0°C
SI.CoefficientOfHeatTransferhcv_B(1/K_global - 1/hcv_A)^(-1)
Geometrical parameters
SI.AreaCrossSectionA4.26Cross section of the pipe for the fluid A
SI.AreaCrossSectionB4.26Cross section of the pipe for the fluid B
SI.AreaExchangeSurface0.0Largest exchange surface area
Flow parameters
Realksi_fixedA1.0
Realksi_fixedB1.0
Exchange parameters
SI.CoefficientOfHeatTransferK_globalGlobal heat transfer coeffcient
SI.CoefficientOfHeatTransferhcv_AHeat transfer coeffcient Wall<->Fluid A

Connectors

TypeNameDefaultDescription
Modelica.Fluid.Interfaces.FluidPort_aport_in_A
Modelica.Fluid.Interfaces.FluidPort_bport_out_A
Modelica.Fluid.Interfaces.FluidPort_aport_in_B
Modelica.Fluid.Interfaces.FluidPort_bport_out_B

Components

TypeNameDefaultDescription
MediumA.ThermodynamicStatestateA_inState of fluid A a inlet
MediumB.ThermodynamicStatestateB_inState of fluid B a inlet
SI.SpecificHeatCapacitycpA
SI.SpecificHeatCapacitycpB
SI.MassFlowRatem_flowAMass flow rate of fluid A
SI.MassFlowRatem_flowBMass flow rate of fluid B
SI.ThermalConductanceQcAThermal flow rate unit
SI.ThermalConductanceQcBThermal flow rate unit
SI.TemperatureTA_in
SI.TemperatureTA_out
SI.TemperatureTB_in
SI.TemperatureTB_out
SI.PowerPexexchanged power
RealNTU_1Number of transfer unit
RealCr_1Ration of thermal condutance
SI.EfficiencyEff_1Exchanger effectiveness
SI.TemperatureTA_out_1
SI.TemperatureTB_out_1
SI.PowerPex_1
SI.TemperatureTA_mid_2
SI.TemperatureTA_mid_buffer
SI.TemperatureTB_mid_2
SI.TemperatureTB_out_2
SI.TemperatureTA_out_2
SI.TemperatureTdew
SI.TemperatureTsat_out
SI.Pressurep_water
RealNTU_2
RealNTU_wetNumber of transfer unit
RealCr_wet
SI.EfficiencyEff_2
SI.EfficiencyEff_wetExchanger effectiveness
SI.AreaS_sensibleSensible surface to achieved saturation on moist air
SI.AreaS_wetSensible surface to achieved saturation on moist air
SI.MassFractionwsat_eq_in
SI.MassFractionwA_in
SI.MassFractionwsat_out
SI.MassFractionwA_out_2
SI.MassFractionXA_out_2
SI.MassFractionwA_outMoisture content peer kg of dry air
SI.SpecificEnthalpyhsat_eq_in
SI.SpecificEnthalpyhA_mid_2
SI.SpecificEnthalpyhA_out_2
SI.SpecificEnthalpyhcond_out
SI.SpecificEnthalpyhA_sat_inEnthalpies peer kg of dry air
SI.SpecificHeatCapacitycp_eqSpecific heat capacity of the fictive fluid
SI.MassFlowRatem_flow_eqMass flow rate of the fictive fluid
SI.PowerQ_flow_wet
SI.PowerQ_flow_dry
SI.PowerQ_flow_wet_2
SI.PowerQ_flow_dry_2
SI.PowerP_sensibleSensible exchanged power
SI.PowerP_latentLatent exchanged power

Contents

NameDescription
MediumA
MediumB