EBSILON®Professional Online Documentation
EBSILON Professional Components / Components - General and Categories / Heat Exchanger / Component 112: Cooling Tower (Mechanical Draft)
In This Topic
    Component 112: Cooling Tower (Mechanical Draft)
    In This Topic

    Component 112: Forced Draft Cooling Tower (Merkel)


     

    Specifications

    Line connections

    or (if FTYPEFAN = 1)

      

    1

    Air Inlet

    2

    Air Outlet

    3

    Cooling Water Inlet

    4

    Cooling Water Outlet

    5

    Makeup Inlet

    6

    Blow-down

    7

    Fan Electrical Power

    8

    Logic-inlet, Control inlet for performance factor (PACKPERF)

    9

    Logic-inlet, Control inlet for reduction of fan speed (T4MIN)

     

    General       Transient modeling       User Input Values       Physics Used       Displays       Example

     

    General

    Component 112 simulates a mechanical draft cooling tower with a counter current wet cooling zone. The wet cooling zone model is based on the Merkel equation (e.g. VDI Wärmeatlas,  section Mj).

     

    There is as of Release 11 an identification mode for these cooling towers.

    The flag FIDENT allows to control whether
    • the cooling water outlet temperature T4 is to be calculated (FIDENT=0)
    • the cooling water outlet temperature T4 is to be specified, and the air flow is to be calculated (FIDENT=1)
    • the cooling water outlet temperature T4 is to be specified, and the Merkel number is to be calculated (FIDENT=2)
    It is also possible to specify a wind correction.

    Also, a kernel expression can be specified for the packing characteristic factor C, the performance factor K, and the wind correction.

    Now the Merkel cooling towers can also be used in the INTMAT mode (integration of the material data into the equation system). Thus for models containing these components a material data reconciliation is possible too.

     

    Annotation

    There is a performance factor PACKPERF by means of which a change of the cooling tower performance can be simulated. This parameter has been made accessible via a logic line, so that an adjustment is possible.


    Transient modeling

    Component 112 enables the modeling of transient regimes (time-dependent computations). This type of calculation is activated with the switch FINST.

    A cooling tower features transient behavior mostly in the water basin which acts as a big water storage. The wet zone with fill or packing material also features some significant mass and acts as indirect storage delaying the cooled water temperature changes e.g. at load changes. Finally, the cooling tower shell wall also acts as an indirect storage influencing however only the air outlet state.

    Cooling tower transient computations are possible for drain mode (FCIRC=1) and without hybrid mode (FHYBRID=0) only. In cooling tower transient computations the air outlet humidity is always computed (FHUM=1).

    The user input values for the transient modeling are classified in the following groups.

    Water basin definition - here the water basin geometric details are specified. The water basin transient mass and energy balance is computed assuming homogeneous mixing of water inside the basin. Depending on the flag value FSPIN one can either specify the water basin level (FSPIN=0) or compute the level specifying both the inlet and the outlet water mass flow (FSPIN=1).

    Fill or packing material definition - here the geometric and material details of the fill or packing material (wet zone) are specified. The fill or packing material part is defined by the volume of the wet zone WZVOL, the void volume fraction PHI (i,e, the volume portion of the wet zone not occupied by the fill or packing material) as well as the fill or packing material specific surface area FV.

    Heat transfer coefficients - in this group the heat and mass transfer coefficients, in particular for the simulation of the evaporation and the heat transfer in the wet zone, are specified. The heat transfer between water and air is controlled by the coefficient ALPHIWAN at Design conditions (defined by the air to water ratio) and is scaled by the corresponding parameters PACKFACT, PACKEXP, PACKPERF at off-design conditions. The mass transfer (water evaporation) is controlled by the coefficient BETAWAN at Design conditions and is scaled by the corresponding parameters PACKFACT, PACKEXP, PACKPERF at off-design conditions. There are two alternatives to specify the mass transfer coefficient. One can either specify the value of BETAWAN directly or by specifying the value of ALPHIWAN and the Lewis factor LEF defined as

    \[ Le_f = \frac{\alpha}{\rho \beta c_p} \]

    Here \(\alpha\) and \(\beta\) are the heat and mass transfer coefficients respectively. \(\rho\) and \(c_p\) are the air density and heat capacity.

    Shell wall material definition - here the geometric and material details of the cooling tower shell wall are specified.

    After spefifying all geometric details of the cooling tower one can match the results of the steady state Merkel equation model implemented in the component 112 for the wet zone by adjusting the values of ALPHIWAN, BETAWAN (or LEF)  at steady state conditions.  

    In case of identification calculation - FIDENT > 0 - the values of the heat and the mass transfer coefficient water to air (\(\alpha\) and \(\beta\)) are adjusted to match the specified temperature T4.

    User Input Values

     

    FFU

    Switch (on/off)

    Like in Parent Profile (Sub Profile option only)
    Expression

    =0: Off (No Air calculations, set DT34N = T3-T4)

    =1: On

    DT34N

    Cooling range (Only for FFU = 0)

    FMODE

    Calculation mode (design / off-design)

    Like in Parent Profile (Sub Profile option only)
    Expression 

    =0: Global

    = 1: local Off-Design

    =-1: local Design

    FIDENT

    Identification mode

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: None

    = 1: T4 given, identify Airflow

    = 2: T4 given, identify Merkel Number

    FWETZONE

    Switch for wet zone mode

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Set cooling water outlet temperature T4

    =1: wet-zone temperature approach

    =2: Merkel number (Me)

    =3: NTU (NTU = Me * M_Water_in / M_Air_dry_in

    PWETZONE

    Parameter for wet-zone

    Cooling water outlet temperature T4 (for FWETZONE=0)

    wet-zone temperature approach (for FWETZONE=1)

    Merkel number (for FWETZONE=2)

    NTU (for FWETZONE=3)

    AWR

    Dry Air / Water Ratio, Design

    FPACKFACT

    Flag for Packing characteristic factor C (off design)

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Internally by the specification value

    =1: PACKFACT= Port_8.P/(1 bar)

    =2: Expression EPACKFACT

    PACKFACT

    Packing-Characteristic Factor (C), off design

    EPACKFACT

    Expression for Packing-Characteristic Factor (C)

    PACKEXP

    Packing-Characteristic Exponent (M), off design

    FPACKPERF

    Specification of performance factor , off-design

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Internally by the specification value

    =1: PACKPERF= Port_8.H/(1 kJ/kg)

    =2: Expression EPACKPERF

    PACKPERF

    Packing-Characteristic Performance Factor (K), off design

    EPACKPERF

    Expression for PACKPERF

    FMERKEL

    Switch for Merkel equation

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Standard Merkel Equation

    =1: Enhanced Merkel Equation

    FHUM

    Handling of wet-zone outlet humidity

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: wet-zone outlet relative humidity is user input

    =1: wet-zone outlet relative humidity will be calculated

    PHI2

    wet-zone outlet relative humidity (FHUM=0)

    DRIFT

    Drift loss fraction

    FCIRC

    Switch for water circulation type

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Closed circulation. Determine makeup and purge flow rate via the parameter COC

    =1: Open circulation. No makeup and purge flow

    COC

    Cycles of concentration in a closed circulation loop

    Hint:

    According to VGB R 455 P is

        COC = makeup water / purge water = Z/A

    This is a simplification that is not valid for components 111/112, because here the calculation is based on the real increase in concentration of salts under consideration of drift losses (DRIFT).

    NBAYS

    Number of bays active (Design)

    FSPECFAN

    Off-design calculation mode for fan speed

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: manual input

    =1: reduce fan speed to have Wet zone water temperature outlet > T4MIN (in off-design)

    =2: reduce fan speed to have Wet zone water temperature outlet > H9 (in off-design)

    T4MIN

    Minimum wet zone water exit temperature (FSPECFAN=1 or 2)

    DTBAND

    Switch hysteresis (FSPECFAN=1 or 2)
    (not used yet)

    NBAYSH

    Number of bays with reduced volume flow

    NBAYSOFF

    Number of bays shut off (water-side bypass)

    FTYPEFAN

    Fan type

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Induced draft

    =1: Forced draft

    DP12

    Air side pressure drop

    FRVMH

    Reduced volume flow fraction

    EFFFANN

    Fan efficiency at nominal volume flow

    EFFFANH

    Fan efficiency at reduced volume flow

    LGEARN

    Gear losses (fraction) at nominal volume flow

    LGEARH

    Gear losses (fraction) at reduced volume flow

    EFFMOT

    Motor efficiency

    FCFWIND

    Specification of wind correction factor

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: None

    =1: value CFWIND

    =2: value CWTDT

    =3: Airflow factor = Port_8.M/(1kg/s)

    =4:Cooling water outlet offset = Port_8.M/(1 kg/s)*(1K)

    =5: Char Line lookup CWINDAIRFLOW, wind speed x = Port_8.M/(1 kg/s)*(1 m/s)

    =6: Charlie lookup CWINDCWTDT, wind speed x = Port_8.M/(1 kg/s)*(1 m/s)

    =7: Expression ECFWIND corrects airflow (factor)

    =8: Expression ECFWIND corrects cooling water temperature (offset)

    CFWIND

    Wind correction factor

    CWTDT

    Wind correction temperature offset  for cooling water outlet

    ECFWIND

    Expression wind correction factor

    FHYBRID

    Switch for hybrid mode

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Off

    =1: Manual input

    =2: Set margin to plume formation (temperature approach)

    MGNPLUME

    margin to plume formation (FHYBRID=2)

    FHX

    Switch for heat exchanger mode

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Off

    =1: Effectiveness

    =2: Temperature change

    =3: Area

    PHX

    Heat exchanger parameter

    Effectiveness (if FHX=1)

    Temperature change (if FHX=2)

    Area (if FHX=3)

    WATRAT

    Dry zone/total water - flow ratio

    AIRRAT

    Dry zone/total air - flow ratio

    KAIR

    Air side heat transfer coefficient

    KWAT

    Water side heat transfer coefficient

    FOUL

    Fouling factor

    PLGASEXP

    Off-design exponent gas side for heat transfer coefficient

    PLWATEXP

    Off-design exponent water side for heat transfer coefficient

    HXFRAC

    Active fraction of heat exchanger (Off-Design)

    FINIT

    Flag: Initializing state

    =0: Global, which is controlled via global variable "Transient mode" under Model Options
          "Extras" ->"Model Options" -> "Simulation" -> "Transient" -> Combo Box "Transient mode"

            (See -> Used Physics / equations -> Global Initialization of Transient Components )

    =1: First run -> Initializing while calculating steady state values
    =2: Continuation run -> Values from previous time step are input for the present ones

    SHEIG Shell height
    SVOL Shell volume
    THWALL Shell wall thickness
    RHOS Density of the shell wall material
    LAMS Heat conductivity of the shell wall material
    CS Heat capacity of the shell wall material
    NFLS Number of points in flow direction: dry zone
    WZVOL Wet zone volume
    PHI Free cross section fraction
    FV Outer surface to volume ratio
    RHOFM Density of fill or packing material
    LAMFM Heat conductivity of fill or packing material
    CFM Heat capacity of fill or packing material
    NFLWZ Number of points in flow direction: wet zone
    FSPIN

    Transient balance calculation mode

    0: Liquid level given, mass flows computed

    1: Mass flows given, liquid level computed

    VF Liquid volume fraction (liquid level) at the end of the time step
    VMIN Volume at liquid volume fraction value 0
    VMAX Volume at liquid volume fraction value 1
    FLVCALC

    Liquid volume calculation mode

    0: linear between VMIN and VMAX

    1: Using ELV

    ELV Function for the liquid volume computation
    TWBEG Temperature in water basin at time step begin (relevant for transient time steps only)
    ALPHIS Inner heat transfer coefficient between air and shell wall
    ALPHIFP Inner heat transfer coefficient between water and fill or packing material
    ALPHIWAN Heat transfer coefficient water to air (nominal value)
    FCALCMT

    Mass transfer coefficient water to air calculation

    0: compute mass transfer coefficient value using Lewis factor input value (LEF)

    1: use directly BETAWAN value and Off-design factors and exponents

    LEF Lewis factor input value
    BETAWAN Mass transfer coefficient water to air (nominal value)
    ALPHO Outer heat transfer coefficient (to ambient)

    AWRN             

    Air to water ratio (nominal)

    MERKELN     

    Merkel number (nominal)

    DP12N            

    Air side pressure drop (nominal)

    RHO1N           

    Ambient air density

    PHI1RAT        

    Ambient air humidity ratio

    PHI2RAT        

    Air outlet humidity ratio

    RHO2              

    Air outlet density

    M2N                

    Outlet air flow (nominal)

    M1N                

    Inlet air flow (nominal)

    MDRYWZ       

    Wet zone dry air mass flow (nominal)

    MAIRHXN      

    Heat Exchanger air mass flow (nominal)

    MWATHXN    

    Heat Exchanger water mass flow (nominal)

    KNAIR            

    Air side heat transfer coefficient (nominal)

    KNWAT          

    Water side heat transfer coefficient (nominal)

    AHX                

    Heat exchanger area (nominal)

    NBAYSN        

    Number of Bays (nominal)

    VMSTACK     

    Stack volume flow (nominal)

    The identification value marked in blue is a reference value for off-design calculations. The actual off-design values refer to the values used in the equations.  
     

    Generally, all inputs that are visible are required. But, often default values are provided.

    For more information on colour of the input fields and their descriptions see Edit Component\Specification values

    For more information on design vs. off-design and nominal values see General\Accept Nominal values

     

     

    Characteristic Lines

     

    Characteristic Line 1, CWINDAIRFLOW: CF(air mass flow) = f (wind speed)

    CF = Correction factor

     

        X-Axis          1         wind speed                         1st point
                            2          wind speed                        2nd point
     
                           .
     
                          N         wind speed                         last point
     
        Y-Axis          1          CF (air mass flow)              1st point
                             2          CF (air mass flow)             2nd point   
                             .                                                                .
                                    CF (air mass flow)              last point

     

    Characteristic Line  2, CWINDCWTDT: dT (cooling water temperature) = f (wind speed)
                                       

     

        X-Axis                   wind speed                                       1st point
                            2          wind speed                                       2nd point

     
                            .
     
                          N          wind speed                                       last point

        Y-Axis          1          dT (cooling water temperature)        1st point
                             2          dT (cooling water temperature)       2nd point
                                                                       .
                                    dT (cooling water temperature)        last point

                                                             



    Physics Used

    The mechanical draft cooling tower consists of multiple bays with a common water basin. Each bay consists of three zones. These zones are the wet cooling zone, the (optional) dry cooling zone and the fan.

    The core to the model of the wet cooling zone is the so-called Merkel Equation [VDI Wärmeatlas, 9th ed., Chapter Mj] which describes the process of combined heat and mass transfer between the water and the air flow through the packing:

    The left-hand part of the equation (the so-called Merkel-Number) describes the ratio of thermal duties to reach the desired cooling duty (with the enthalpy difference between saturation and the respective condition of the air as basis) while the right-hand side reflects the properties of the packing and the operating mode of the cooling zone.

    In the design mode, the user can specify through the parameters for the calculation mode FWETZONE and the respective parameter value PWETZONE if the characteristics of the cooling zone shall be pre-defined (either as Merkel-Number or NTU Number which is equal Me/AWR), or if it shall be calculated from inputs of either the water exit temperature or the approach to the wet bulb temperature of the air.

    In order to model also hybrid cooling towers, which are used to avoid plume formation, the dry cooling zone can be activated by selecting the hybrid mode in the parameter FHYBRID.

    In design mode, the distribution of the cooling duty between the dry and the wet zone of the hybrid cooling tower will be set either by definition of the air and water flow ratios or by the desired margin to plume formation (specified in degrees difference to the dew point).

    In off-design mode, the air distribution will be calculated based on the definition of the fan operation (number of fans at half-speed or off, respectively) and the pressure loss characteristics.

    The physics of the dry zone are modelled according to the rules for a normal air-water cross-flow heat exchanger, the physics of the wet cooling zone according to the rules of Merkel's equation.

    The part-load performance of the wet zone is a characteristic of the Merkel number as a function of the dry air to water ratio.

    The energy balance of the common water basin including the make-up water determines the final exit temperature of the cooled water stream.

     

     

     Wet zone

    Design:

    Water outlet temperature given:

    The Merkel Number will be calculated, by integration of the Merkel equation.

    Merkel number given:

    Find the corresponding water outlet temperature to get the desired Merkel number.

    Off-Design:

    Via the relationship: Me = MeDesign * K*C*(AWR/AWRDesign)m the actual Merkel number will be calculated. Then a solver finds the corresponding water outlet for this Merkel number. K, C and m are defined as PACKPERF, PACKFACT and PACKEXP.

    Integration of the Merkel equation:

    If the wet zone outlet humidity is set by the user (FHUM=0) the Merkel number can be determined by numerical quadrature. An Implementation of algorithm 699 from TOMS; TRANSACTIONS ON MATHEMATICAL SOFTWARE, VOL. 17, NO. 4, DECEMBER, 1991, PP. 457-461; is used.

    If the wet zone outlet humidity is calculated (FHUM=1) a System of ordinary differential equations needs to be solved. Here a Runge-Kutta Method of order 5 is used (Dormand-Prince).

     Heat exchanger

    The heat exchanger is modelled as single pass with the  NTU-Effectiveness Method. (VDI Wärmeatlas, Section Ca; „Compact Heat Exchangers by W.M. Kays and A.L. London)

    The Off-Design correlation for the heat transfer coefficients is

    for the air-side:  

    k Gas Off-Design = k Gas Design *(v Gas Off-Design / v Gas Design)m1 (Default m1=0.8,  m1 corresponds to PLGASEXP)

    for the water side:

    k Water Off-Design = k Water Design *(v Water Off-Design / v Water Design)m2 (Default m2=0.8, m2 corresponds to PLWATEXP)

    Water density is considered constant

     

    Fan and air flow:

    Design:

    In design mode, air flow and pressure drop are user input (AWR and DP12).
    The volumetric flow through the fan will be calculated.

    Off-Design:

    The model assumption is a constant volumetric flow rate at nominal fan speed.
    At reduced fan speed the volumetric flow rate will be scaled by the parameter  FRVM.

    The relationship for the pressure drop is described by the following equations:

     

    with rho = density, and A = cross-sectional area:

    and elimination of leads to the final off-design correlation for the pressure drop:

     


    Component Displays

    Display Option 1

    Example

    Click here >> Component 112 Demo << to load an example.