EBSILON®Professional Online Documentation
EBSILON Professional Components / Components - General and Categories / Steam Generator / Component 89: Steam Generator - Main Heating Surface
In This Topic
    Component 89: Steam Generator - Main Heating Surface
    In This Topic

    Component 89: Steam Generator - Main Heating Surface


    Specifications

    Line connections

    1

    Cold side, inlet

    2

    Cold side, outlet

    3

    Heat flow from flue gas

    Connection to the attached flue gas zone, component 88 or component 90

    4

    Radiation exchange with component upstream in flue gas path

    5

    Radiation exchange with component downstream in flue gas path

    6

    Connection to auxiliary heating surface 1

    7

    Connection to auxiliary heating surface 2

     

    General       User Input Values       Characteristic Lines       Physics Used       Displays       Example

     

     

    General

    Component 89, together with components 88 (steam generator, flue gas zone) or 90 (steam generator, reaction zone) allows the modeling of the heat transfer in a heating surface of a steam generator.

    It can be used as economizer, evaporator or super heater. It is used to model membrane walls in radiation areas or combustion chambers. I can also be used to model tube bundles. Water cooled walls or supporting tubes in the area of tube bundles can be represented using component 91 (steam generator, auxiliary heating surface).

    Component 89 maps only the cold side (water / steam) side of the heating surface. Usage without components 88 or 90 to model the secondary (flue gas) side will therefore result in errors. Besides water/steam, also thermal oil, two-phase fluid, and universal fluid are supported as primary fluid.

     

    Extensions for the Pressure Drop Calculation of the Flue Gas Side (for Components 88 and 89)

    In Component 89, the pressure drop calculation of the flue gas side can be specified as Kernel expression, EDPO. As geometrical details of the heating surface are specified in Component 89, it is possible to access the corresponding specification values (like e.g. pipe diameter, partition, details on fins, etc.) and to implement an own pressure drop calculation (e.g. according to manufacturer information) in the Kernel expression EDPO.

    The pressure drop on the flue gas side is calculated in Component 89 and transferred to Component 88. In Component 88 there is a flag for this, FDP, which switches between the previous pressure drop calculation according to Bernoulli’s principle and EDPO of Component 89.

     

    Extensions for the Calculation of the Heat Transfer Coefficients of the Heating Surfaces (for Component 89)

    As an alternative to the previous implementation according to VDI formalisms, the calculation of the individual heat transfer coefficients of the heating surface (alpha numbers) can be implemented by means of a corresponding Kernel expression (e.g. according to manufacturer information) in Component 89.

    The following pair of specification values exists for this purpose in Component 89:

     

    For the bundle heating surface the mean fluid flow velocity in the pipe is output as result value VELM.

    MODEL:

    Like other heat exchangers, the calculation model is based on the fundamental equations for the mass flow, pressure and enthalpy

    -  (1)      M1 =  M2

    -  (2)      P1 = P2 + DP12

    -  (3)     (k*A)*LMTD+QR_4+QR_5 = M2*H2 - M1*H1

    with the indices

    - 1 for the primary inlet (cold side)

    - 2 for the primary outlet (cold side)

    - 4 radiation flow from predecessor (with respect to flue gas arrangement)

    - 5 radiation flow from successor (with respect to flue gas arrangement)

    For calculating the mean logarithmic temperature difference LMTD, the flue gas inlet and outlet temperatures of that component 88 or 90 are used, which is connected using logic line 3. Therefore component 89 will work properly only when (PIN 3) is connected to (PIN 3) of component 88 or 90.

    In case of mass flow and pressure, the component assumes that M1 / M2 or P1 / P2 respectively are defined on the connected lines. In each case, the missing value is calculated.

     

    Pressure drop limitations in off-design (Extras --> Model Options--> Calculation -->Relative pressure-drop maximum) :
    As the pressure drop rises quadratically with the mass flow, pressure drops that are significantly too high can quickly arise in the event of a transgression of the nominal mass flow. These will then cause phase transitions and convergence problems. For this reason, pressure drop limitations have been installed.

     

    Pinch point violation :
    In case of a pinch point violation KA is reduced automatically (like for other heat exchangers) to avoid the pinch point violation. A warning is issued in this case.

     

    CALCULATION RULES

    The heat transfer coefficient  k*A is calculated from geometric data of the heating surface as well as process data (mass flows, temperatures, pressures,  flue gas composition) following published calculations rules. In this aspect, component 89 differs from the other heat exchanger components, where k*A is calculated in design mode from given enthalpies or terminal temperature differences (see HeatExchanger,General Notes ). Depending on the parameter FIDENT, the behaviour of component 89 is as follows:

    For the identification mode there are new variations FIDENT=“T2 given validable“, both for calculating the fouling and for calculating the rating number. Previously, for the validation a temperature specified from outside was treated as a constant quantity. In the new modes, an equation is generated into which this temperature is integrated, so that a statistical reconciliation can take place.

    Thus, the calculation mode DESIGN or OFF-DESIGN has no influence on the behaviour of component 89 concerning the heat transfer calculation.

    The k-value results from the inner heat transfer ALI, the outer heat transfer by radiation ALR and by convection ALC as well as from heat conductivity LAMDBA, efficiency factor EFF and fouling factor KFOUL:

    AA, AI, AM are the specific surface areas of the heat-exchanger tubes with respect to outer, inner and average surface. S is the wall thickness.

     

    LAMBDA is stored taking into account the dependency on material and temperature. The heat transfer coefficients are calculated from the VDI Heat Atlas (7th edition, 1994), where a differentiation is made between the tube bundles, radiation areas and combustion areas:

    • According to Gnielinski (VDI-Gb) for single phase flow,                 
    • Bubble boiling in vertical tubes (VDI-Hbb) for evaporation zones

    The calculation rule for the inner heat transfer can be chosen by the switch FPHASE. Alternatively, there is the option to switch automatically according to the enthalpies at inlet and outlet.      
    If FPHASE = AUTOMATIC, then bubble boiling is assumed, if (H1+H2)/2 is in the wet steam area.

                          The convective heat transfer is set to 0 for radiation zones.

    • gas radiation and effective layer thickness according to VDI-Kc,

                          Radiation fractions of coke, ash, soot according to VDI-Ke. Particle radiation is taken into account only in the combustion zone.

    Besides the heat transfer by radiation within a flue gas zone, component 89 can also take into account the radiation exchange between a heating surface and the flue gas volume of components upstream or downstream, in the direction of  flue gas flow.  For this purpose, PIN 5 must be attached to PIN 7 of  the components 88 or 90, which represents the flue gas part of the successor and, if necessary, PIN 4 with PIN 6 of the predecessor. The calculation of the radiation exchange is based on the following assumptions:

    • for tube bundle or radiation zone, the average of flue gas inlet and outlet temperature
    • for combustion zone, it can be adjusted between adiabatic temperature (ALPHA_T=1.0) and flue gas outlet temperature (ALPHA_T=0.0) using the model parameter ALPHA_T
    • according to VDI-Kb for radiation and combustion zones,  
    • for bundles of plain tubes da/tl*n/(n+1) with da = outer diameter, tl = transversal pitch, n = number of tube layers,
    • for finned tubes 1.

     

    COMPONENT IDENTIFICATION

    Identification of the component means calculating EFF or KFOUL using e.g. equations (3) and (5), if one of the flue gas temperatures and both water / steam in- and outlet temperature are given. Currently, the component allows only for calculating the efficiency factor EFF with a given fixed fouling KFOUL. The calculated efficiency is displayed as the result REFF.

    In a steam generator there are sometimes heating surfaces, which form a unit for the process (e.g. reheater 1 or evaporator), but which consist of several constructional groups, installed in separate flue gas zones. Most of the time there is no temperature / enthalpy information available from measurements among the constructional groups. Thus, an identification of each group using e.g. equations (3) and (5) is not possible. However, a common efficiency factor for the unit as a whole can be calculated  from the temperatures at inlet and outlet.

    For this purpose, a group index (IGROUP) for the identification can be assigned in component 89. All heating surfaces, which should have a common efficiency factor, have the same group index, and up to 10 groups are allowed, 1<=IGROUP<=10.  One of the heating surfaces within the group is then set in the identification mode (FIDENT = T2 given, efficiency factor calculated), all other ones in FIDENT = no identification. All heating surfaces within a group defined in this way will have the same efficiency factor as a result.

    In case the model is to be used in the simulation mode without an identification of the groups, then for one heating surface in the identification mode, FIDENT should be set to "no identification" and for ALL heating surfaces of the group, IGROUP = 0 is to be set.

     

     AUXILIARY HEATING SURFACES

    In the flue gas zone, there are frequently smaller heating surfaces along with the main heating surfaces in the form of construction parts (supporting tubes, cooled walls), through which water and steam flow. These auxiliary heating surfaces are often arranged in a different part of the water / steam circuit than the main surface, but share the same flue gas inlet and outlet conditions with the main heating surface. Component 91 is designed to represent these auxiliary heating surfaces.

    The heat transfer in a flue gas zone is dominated by the geometry of the main heating surface. Thus, the heat transfer coefficient calculated in component 89 is transferred to component 91, in case it is connected. To support this, an optional auxiliary heating surface (with its PIN 4) can be attached to the main heating surface using the logic lines PIN 6 or 7 of component 89. A maximum of two auxiliary surfaces for each component 89 is allowed.

    Using an auxiliary heating surface in a flue gas zone in addition to a main heating surface makes sense only if the main heating surface is configured as a bundle tube by setting FHEAT = bundle. If FHEAT = radiation / combustion zone, the main heating surface normally represents the membrane wall.

     

    Note the result value REFF

    In the case of Component 89 the possible heat exchange is explicitly calculated at a specified geometry. REFF represents here the ratio between the real (fouled) and ideal (clean) heat exchanger. Thus here REFF is a measure for the condition (performance factor) of the heat exchanger.

     

    Further result values and plausibility checks

    The pressure drop is displayed as result value DP, independently of whether it was calculated in the component or inlet and outlet pressure were given from the outside.

    Note: 

    The formulae used in the component are only suitable for Reynolds numbers <10^6. The Reynolds number is limited to 2 * 10 ^ 6 and a warning will be output.

     

    Further results are:

    are now output as result values,

    For a better comprehension of the heat transfer by radiation on the flue gas side, the following result values are calculated and displayed besides

    the following result values are calculated and displayed besides

    Minimum terminal temperature difference

    The minimum value for the terminal temperature difference has been reduced from 3 to 1 K.

     

    CHARACTERISTIC LINES AND ADAPTATION POLYNOMIAL

    The part-load behaviour of component 89 results from the implemented model equations. For an adaptation to the real plant behaviour, characteristic lines and adaptation polynomials are provided. The switches FCHR and FADAPT allow you to choose between a characteristic line and an adoption polynomial.

    For the emissivity and the absorption ratio it is possible for the user to implement their own calculation algorithms by means of kernel expressions. The previous implementation is stored as default.


     The list of available types of tube and fin material has been expanded by 58 new types of material derived from a research paper published by Mannesmann.


     

    New Result Values

    To make it easier to determine the characteristic line, result values M1M1N and EFFEFFN have been added in this component; they show the ratio of the prevailing mass flow and the prevailing effectivity respectively to the corresponding nominal value.

     

    User Input Values

     

     

    FFLOW

    Flag for the heat-exchanger type (see HeatExchanger,General Notes ).

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: Counter Current flow

    =1: Concurrent flow

    FMODE

    Flag for calculation mode Design/Off-design

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: set as global
    =1: local off-design (always off-design, even if design is set globally)
    =-1: Local design

    FDP

    Flag for the calculation of the pressure loss

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: calculated in the component from DP12N
    =-1: inlet and outlet pressure given externally

    DP12N

    Pressure loss between ports 1 and 2 in design mode. If design mode is set, the prescribed value is used. If off-design is set, pressure will be corrected with mass flow and specific volume.

    FIDENT

    Flag for component identification

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: no identification or identification in a group with another component 89 with a common efficiency factor. In case of a group identification, the related components shall have the same IGROUP number. In the case of no identification, IGROUP shall be 0

    =1: T2 given (constant), efficiency is calculated. In case identification of several components 89 is to be done with a common efficiency factor, exactly one component within the group shall have FIDENT = 1, the others shall have FIDENT=0.

    =2: T2 given (constant), fouling is calculated. This option is currently not implemented.

    =3: T2 given (validable), efficiency is calculated.

    =4: T2 given (validable), fouling is calculated. This option is currently not implemented.

    =5: T2 and efficiency factor given, heating surface area calculated (The calculated surface is output in the result value RA.)

    IGROUP

    Group Index: Component 89 allows for a group identification of several components with a common efficiency factor. The membership of a component 89 in such a group is marked by a common group index, 0<=IGROUP<=10. IGROUP = 0 does not mark any group, but must be set, if no identification is to be done. A maximum of 10 identification groups are possible.

    DIAM

    Inner diameter of the heat exchanger tubes

    THWALL

    Wall thickness of the heat exchanger tubes

    PITCHL

    Longitudinal pitch of the tube bundle (distance of successive tubes in the flow direction of the flue gas. PITCHL will be evaluated only if FHEAT = tube bundle. It is not relevant if FHEAT = radiation or combustion zone, because component 89 represents a membrane wall in these cases)

    PITCHT

    Transversal pitch of the tube bundle (distance of successive tubes perpendicular to the  flow direction of the flue gas. PITCHT will be evaluated only if FHEAT = tube bundle. It is not relevant if FHEAT = radiation or combustion zone, because component 89 represents a membrane wall in these cases)

    FTUBE

    Tube arrangement

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: inline

    =1: staggered

    FTUBE will be evaluated only if FHEAT = tube bundle. It is not relevant if FHEAT = radiation or combustion zone, because component 89 represents a membrane wall in these cases

    KFOULN

    Additional thermal resistance to be taken into account on calculating the heat transfer coefficient. Nominal value that can be adjusted using characteristic lines and/or adaption polynomial depending on the specification values FADAPT and FCHR

    FHEAT

    Code for heating surface type:

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: tube bundle, represents tube coils that are arranged perpendicular to the direction of flue gas flow

    =1: radiation zone, represents cooled water walls (membrane wall); radiation zones do not house tube bundles across the direction of flow of the flue gas. Thus, specification values referring to the geometry of tube bundles are meaningless; convective heat transfer is ignored.

    =2: combustion zone, similar to the radiation zone; the radiation zone takes into account radiation emission of gas, ash and soot, the combustion zone in addition the emission of fuel particles (coke)

    FMATT

    Code for tube material :   Material Properties of Steel

    =-1 : Lambda calculated by Kernel expression EMATT
             
    (Own algorithm for calculating the thermal conductivity of the material of the fin. ) 

    EMATT

    Function for heat conductivity for tube material

    function evalexpr:REAL;
                var
                t_wall:real;
                begin
                t_wall:=keGetInternal("T_WALL");
                t_wall:=max(20.0,min(t_wall,600.0));
                evalexpr := -0.0000286963*t_wall*t_wall - 0.0186106198*t_wall + 58.2256566351;
                end;

    FPHASE

    Code for calculating the inner heat transfer type

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: Automatic; depending on the average enthalpy 0.5*(H1+H2), one of the other options is selected

    =1: Single phase, water; calculation according to Gnielinski under the assumption of a single phase flow in the heating surface. In case the average enthalpy is in the wet steam area, the heat transfer is calculated with water.

    =2: Single phase steam;  calculation according to Gnielinski supposing single phase flow in the heating surface. In case the average enthalpy is in the wet steam area, the heat transfer is calculated with water.

    =3: Nucleate boiling; Calculation of a heat transfer coefficient that depends only on the heat flux density by assuming bubble boiling within the complete heating surface.

    A

    Heat exchange area of the tube bundle or the membrane wall

    FFIN

    Code for fin type ( Relevant only if FHEAT = 0): 


    Like in Parent Profile (Sub profile option only)

    Expression

    =0: Plain pipe

    =1: Circular fins

    =2: Rectangular fins

    =3: Continuous fins

    CFIN

    C- value for the calculating the outer heat transfer coefficient of finned tubes according to Schmidt: NU =CFIN * RE^0.625 *PR^0.333.Typical value is CFIN=0.28  for in-line, CFIN=0.3333 for staggered arrangement.

    HFIN

    Fin height, relevant only if FFIN = CIRCULAR FIN

    LFIN

    Fin length

    - for rectangular fins LFIN>=WFIN;
    - Connected fins (Length and width are determined by the arrangement of the pipes) :

    •   In the case of an aligned arrangement of the pipes, the fin length equals the longitudinal pitch of the pipes
    •   In the case of a staggered arrangement, each pipe is allocated a hexagonal section of the fin plate.
        The fin length thus results as follows:
      • LFIN = (PITCHL^2 +0.25*PITCHT^2)^0.5

    (see VDI Wärmeatlas, Kapitel M1)

    WFIN

    Fin width (Rectangular fins)

    - for Rectangular fins WFIN = Specification value

    - Connected fins (Length and width are determined by the arrangement of the pipes) :

    •   In the case of an aligned arrangement of the pipes is the fin width equals the transverse pitch of the pipes.
    •   In the case of a staggered arrangement, the fin width results as follows:
      •  WFIN = PITCHT, falls PITCHL>=0.5*PITCHT
      •  WFIN = 2*PITCHL, falls PITCHL<0.5*PITCHT

    (see VDI Heat Atlas, Chapter M1)

     

    PFIN

    Fin pitch (number of fin / tube length)(distance between neighbouring fins), relevant only if FFIN = CIRCULAR FIN

    THFIN

    Fin thickness, relevant only if FFIN = CIRCULAR FIN

    FMATF

    Code for fin material :   Material Properties of Steel
     

    =-1 : Lambda calculated by Kernel expression EMATF
             
    (Own algorithm for calculating the thermal conductivity of the material of the fin. ) 

    Relevant only when FFIN=1.

    EMATF

    Function for heat conductivity for fin material

    PHIINC

    Angle between the flow direction of the flue gas and the tube coils; valid only if FHEAT = 0; the standard value 90° corresponds to perpendicular flow.

    NTUBES

    Number of tubes, through which flow is parallel on the water / steam-side; does  not necessarily always correspond to the number of tubes in one layer (seen from the flue gas side), but it is usually an integer multiple of the same (flow = 1/2, 1, 2, 3, 4, ...)

    EFFN

    Efficiency factor for correcting the inaccuracies in the model or the contamination states in accordance with: EFF/K = 1/ALO + 1/ALI*AA/AI + S/LAMBDA*AA/AM + KFOUL; depending upon the specification values FADAPT and FCHR, it can be corrected with characteristic curves and/or adaptation polynomials.

    CSASH

    Relative cross-section of the ash for absorptions, typical value 0.25

     

    The parameters CSASH, DIAASH and DISTASH configure the fraction of ash radiation in the total emission of the flue gas.

    DIAASH

    Mean particle diameter of the ash, typical value 16*10**-6m

     

    The parameters CSASH, DIAASH and DISTASH configure the fraction of ash radiation in the total emission of the flue gas.

    DISTASH

    Control parameter of the Rosin-Rammler distribution of the ash particles, typical value 1.5

     

    The parameters CSASH, DIAASH and DISTASH configure the fraction of ash radiation in the total emission of the flue gas.

    SOOTCONT

    Soot content of the flue gas; used for calculating the fraction of soot emission in solid form in the total emission of the flue gas

    FADAPT

    code for using the adaptation polynomial ADAPT/ adaptation function EADAPT

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: not used and not evaluated

    =1: EFF=EFFN*char. line* polynomial

    =2: EFF=EFFN*polynomial; if FCHR=0 is used at the same time, the characteristic line is ignored.

    =3: KFOUL=KFOULN*char. line*polynomial

    =4: KFOUL=KFOULN*polynomial;  if FCHR=1 is used at the same time, the characteristic line is ignored.

    =5: KA= polynomial*KAN

    =1000: not used but ADAPT evaluated as RADAPT (Reduction of the computing time)

     

    = -1: EFF=EFFN*char. line* adaptation function

    = -2: EFF=EFFN*adaptation function; if FCHR=0 is used at the same time, the characteristic line is ignored.

    = -3: KFOUL=KFOULN*char. line*adaptation function

    = -4: KFOUL=KFOULN*adaptation function;  if FCHR=1 is used at the same time, the characteristic line is ignored.

    = -5: KA=adaptation function*KAN

    = -1000: not used but EADAPT evaluated as RADAPT (Reduction of the computing time)

    EADAPT

    Adaptation function (input)

    FCHR

    Code for using the characteristic lines

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: EFF/EFFN=f(M1/M1N); characteristic line is ignored, if FADAPT=2 is set.

    =1: KFOUL/KFOULN=f(M1/M1N); characteristic line is ignored, if FADAPT=4 is set.

    FVALKA

    Validation of k*A (only in off-design)

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: KAN is used without validation
    =1: Pseudo-measurement point is used, indicated by IPS

    IPS

    Index of the pseudo-measurement point for the validation of k*A

    FEM

    Flag for Calculation of emissivity and absorption ration

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: By internal formulas

    =1: By internal formulas with correction factors EEM and EABS

    =2: By Kernel expressions EEM and EABS

    EEM

    Function for emissivity

    EABS

    Function for absorption

    EMWALL

    Wall emissivity

    FALI

    Flag for switching the calculation of the internal heat transfer coefficient according to EALI on and off respectively

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: By internal formulas

    =1: By Kernel expressions EALI

    EALI

    Kernel expression for the internal heat transfer coefficient (result value ALI)

    FALC

    Flag for switching the calculation of the external heat transfer coefficient by convection according to EALC on and off respectively

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: By internal formulas

    =1: By Kernel expressions  EALC

    EALC

    Kernel expression for the external heat transfer coefficient by convection (result value ALC).

    FALR

    Flag for switching the calculation of the external heat transfer coefficient by radiation according to EALR on and off respectively

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: By internal formulas

    =1: By Kernel expressions  EALR

    EALR

    Kernel expression for the external heat transfer coefficient by radiation (result value ALR)

    EDPO

    Function for external (flue gas) pressure drop (see under "General")

    M1N           

    mass flow at inlet in the design case

    V1N           

    Specific volume at the inlet in the design case

    KAN           

    Heat transfer coefficient (nominal)

     

    The identification value marked in blue is the reference value for the off-design mode. 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

    It can be seen in the following table, whether T2, efficiency factor or heating surface area is to be specified:
     

     

    Design

    Off-design

     

    FIDENT = no identification

    FIDENT = T2 specification, efficiency factor calculation

    FIDENT = T2 and efficiency factor specification, heating surface area calculation

    FIDENT = no identification

    FIDENT = T2 specification, efficiency factor calculation

    FIDENT = T2 and efficiency factor specification, heating surface area calculation

    efficiency factor

    EFFN to be specified by user, REFF=EFFN

    ignored, REFF calculated

    EFFN to be specified by user, REFF=EFFN

    EFFN and char. line or adaption polynomial to be specified by user, REFF = f(EFFN, char. line or polynomial)

    ignored, REFF calculated

    EFFN and char. line or adaption polynomial to be specified by user, REFF = f(EFFN, char. line or polynomial)

    T2

    calculated

    to be specified by user

    to be specified by user

    calculated

    to be specified by user

    to be specified by user

    heating surface area

    to be specified by user

    to be specified by user

    calculated

    to be specified by user

    to be specified by user

    calculated

     

     

    The following medium combinations are possible:

     

    Primary

    Secondary

    Water or steam

     

    in the connected component 88 or 90:

    Flue gas

     


    Characteristic Lines

    1st characteristic line     FAC/FAC0 = f (M1/M1N)

     (EFF/EFFN = f(M1/M1N) (FCHR=0) or KFOUL/KFOULN=f(M1/M1N) (FCHR=1)

     

    Characteristic line 1: (k*A)-characteristic line:  (k*A)1/(k*A)N = f (M1/M1N)

     

         X-axis          1         M1/M1N                      1st point
                            2         M1/M1N                     2nd point
                            .
                            N        M1/M1N                      last point
     
         Y-axis          1         (k*A)1/(k*A)N              1st point
                            2        (k*A)1/(k*A)N               2nd point
                            .
                            N       (k*A)1/(k*A)N                last point
     

      


    Physics Used

    Equations

    Design

    (Simulation flag:

    GLOBAL = Design

    and

    FMODE = Design)

     

    P2 = P1-DP12N                   

     

      

    Off-design

    (Simulation flag:

    GLOBAL = Off-design

    or

    FMODE = off-design)

     

    DP12 = DP12N * V1/V1N * (M1/M1N)**2

    P2 = P1 - DP12                     

     

      

    All cases

     

     

    M1 = M2                                                                

     

    -M1*H1+M1*H2-M3*H3-M4*H4-M5*H5= 0                    

     

    Calculation of average pipe values

    t_rg = (T3+T4)/2.0

    p_rg = (P3+P4)/2.0

    h_med = (H1+H2)/2.0

    p_med = (P1+P2)/2.0

    t_med = (T1+T2)/2.0

     

    Estimation of wall temperature

    t_wall = t_rg - 0.9*(t_rg-t_med)

     

    Calculation of layer thickness

    according to VDI-Kc

     

    Calculation of form factors for predecessor and successor

    according to VDI-Kb

     Fixed point iteration 1:

    (efficiency factor)

     Fix point iteration 2:

    (Wall temperature with average pipe values)

    calculate emissivity

    according to VDI- Kc (gas radiation) or according to VDI-Ke (combustion zones)

    heat transfer through convection

    according to VDI-Gf (plain tube) or VDI-Mb (finned tube)

    heat transfer through radiation

    ALR = BOLTZMANN*EMMIWAND/(1-(1-EMMIWAND)*(1-absorptionsvh))

    *(emissivitaet*t_rg+TABS)**4-absorptionsvh*(t_wand+TABS)**4)

    /(t_rg-t_wall);

    inner heat transfer

    according to VDI-Gb (single phase) or VDI-Hbb (bubble boiling)

    k-value

    EFF/K = 1/ALO + 1/ALI*AA/AI + S/LAMBDA*AA/AM + KFOUL

    Wall temperature

    t_wall = t_rg - kvalue/(BEWERTUNG*alpha)*(t_rg-t_med)

    END Fix point iteration 2

     

    IF (no identification)

     Calculation of H2 from  (k*A)*LMTD +QR_4+QR_5= M2*H2 - M1*H1

    ELSE

    Calculation of (k*A)_ist from (k*A)*LMTD +QR_4+QR_5 = M2*H2 - M1*H1

    Calculation of EFF from k_ist and EFF/K = 1/ALO + 1/ALI*AA/AI + S/LAMBDA*AA/AM + KFOUL

     

    END Fix point iteration 1

     

    H3 = Q2-Q1-QR_4-QR_5     

     

     

     

     

     

        

     

     


    Component Displays

    Display Option 1

    Example

    Click here >> Component_89 Demo << to load an example.