EBSILON®Professional Online Documentation
EBSILON Professional Components / Components - General and Categories / Turbines / Component 6: Steam Turbine / General Expander
In This Topic
    Component 6: Steam Turbine / General Expander
    In This Topic

    Component 6: Steam Turbine / General Expander


    Specifications

    Line connections

    1

    Steam Inlet

    2

    Steam Outlet

    3

    Extraction 1                   

    4

    Extraction 2  Steam /"NONE"

    5 (8)

    Shaft inlet  Shaft / "NONE"

    If shaft inlet 5 should be active, then that is
    Shaft outlet 8 = "NONE" to be set.

    6

    Shaft outlet

    7

    Control inlet for efficiency (as H)   Control inlet / NONE

    8

    Shaft outlet   Shaft / "NONE"

    If shaft outlet 8 should be active, then that is
    Shaft outlet 5 = "NONE" to be set.

     

    General       Shaft connection       Calculation       Total vs Static Isentropic Efficiency       User Input Values       Characteristic Lines       Physics Used       Displays      Results      Displays       Example

    General

    Component 6 serves the purpose of converting thermal/potential energy of a process stream into mechanical energy to a shaft. It can be applied to water (incompressible flow, hydraulic turbine), to steam or any other gases as defined by stream type Gas, Flue Gas, Universal Fluid, and 2-phase vapor/liquid ( compressible flow, turbo machinery). As such it is a most versatile component in EBSILON to model energy conversion by means of expansion of a process stream.

    Component 6 represents a single expansion stage, a stage group, or a complete expansion section of the modelled equipment. Any extraction or admission along the expansion process must be modelled using multiple instances of component 6.   Admission must be added by means of a mixer into the connecting process stream between two adjacent components 6. Extraction streams can be modelled, using any or both of ports 3 and 4, which have been added for convenience to spare addition of a splitter.

     

    Shaft connection

    Previously, the second shaft pin of Component 6 (Steam turbine) was a shaft inlet. Thus in the case of Components 6, it was to connect several turbine disks in series, so the shaft power added up..

    The reverse case (power distribution for Component 6) could not be graphically represented so far. However, there was a flag FQ (power flow) for the calculation that allowed to switch it over although it had the imperfection that the graphic representation did not match the calculation then.

    It is now possible, to also graphically represent the reverse power flow has now been created for these three components. An additional shaft pin has been implemented for this purpose:

    To enable a corresponding display, the previously existing pin can now be hidden and the new pin has been positioned in the same place. Usually, either the inlet or the outlet will be used and the unused pin should be hidden then. As a rule, however, the software also allows to simultaneously use both pins.

    As in the case of the previous second shaft pin, the power must be specified on the new pin too. The turbine can only calculate the power on the main shaft outlet.

    The new pin has made the flag FQ obsolete. For reasons of compatibility, however, it is still available but has been marked as “deprecated“. Moreover, a comment may be output that points out the options for using the new pin. For the new shaft pin, the flag produces a reversal of the calculation direction as well.

    A new QSHAFT result value has been implemented for component 6 (also for component 23, 58), which outputs the shaft power generated in the component, regardless of which connections it is distributed over or which shaft power is added.

    Calculation

    The calculation of the steam turbine has two objectives:  the determination of (1) the flow characteristics which describe the correlation of throughput and inlet pressure, and (2) the power output which is determined as shaft power by means of an efficiency model.

    Pressure-Flow Characteristic

    The flow characteristic (inlet pressure as function of mass flow) is determined according to Stodola's law, which in essence means that the flow coefficient at the inlet, defined by

    is constant for all operating modes. With m1 = inlet flow, p1 = inlet pressure and v1 = inlet specific volume.

    In design mode, CF will be established using the inputs for flow pressure and specific volume. The user has two choices to specify the inlet pressure:

    The outlet pressure P2 is always defined by external components. An external component can be the next expansion stage, a condenser, or a direct input be means of component 33 of 46.

    P1 - Calculation in part load

    At partial load, component 6 calculates the inlet pressure p1 as a function of the mass flow, outlet pressure and its specific volume from the Stodola law: 

    See also and for the formulation of the Stodola law: Turbines - OffDesign - Stodola 

    In the chapter "Part load calculation of the steam turbine", M1N, P1N, P2N and V1N designate the nominal values in the design case or M1, P1, P2 and V1 the corresponding quantities under the current conditions. As in the design case, the outlet pressure P2 is always determined by external components.

    Power Calculation

    Component 6 uses an isentropic efficiency model to calculate the enthalpy drop during the expansion.

    Figure 1: Enthalpy-entropy diagram of a simple expansion stage

    The mechanical energy output at the shaft evaluates as follows:

    Using isentropic efficiency this converts to:

    With m1 = inlet flow, h1 and h2 = inlet and outlet enthalpy, ηis = isentropic efficiency and mechanical losses as defined by QlossM in absolute terms or ηmech in relative terms. In design mode the user specifies exit enthalpy either directly by setting it outside component 6 on the exit stream by means of components 33 or 46 or by direct input of the isentropic efficiency with the parameter ETAIN. The off-design mode estimates the off-design efficiency relative to the nominal efficiency by means of correction curves (characteristic lines), which correlate the change of efficiency with the ratio of mass flow rate, or the ratio of volumetric flow rate, or the of change expansion pressure ratio.

     

    Identification modes

    Previously, identification modes of the steam turbine could be activated by a negative value of the flag for the characteristic line type (FCHR). This, however, is not just difficult to find for the user but also has the disadvantage that when determining the result value ETACL is is unknown which characteristic line type should be used. Previously ETACL was calculated for the characteristic line type FCHR=0 (mass flow-dependent characteristic line) for historical reasons.

    Analogous to other components, a flag FIDENT has now been implemented also for the steam turbine.  It has the following settings:

    To prevent the behaviour of existing models from changing, the identification modes can still be activated by negative FCHR settings. In this case, the FIDENT setting is ignored. You will be notified by a comment.

    When using FIDENT, the characteristic line type set in FCHR is now used for calculating ETACL.

     

    Calculation of the mass flow from the power output

    It is possible to calculate the required steam mass flow for an externally given shaft power. This allows to directly model a feed water pump with power turbine. A controller had to be used to harmonize the required power of the pump with the turbine output.

    This mode of calculation is activated by setting the new flag FSPECQ to 1. However, it can only be used for a single turbine disk. If several turbines and extractions exist, a controller has to be used for setting a required power output.

    Total vs Static Isentropic Efficiency

    It is important to note that in EBSILON property calls are conducted without any reference to the flow velocity of the process streams, consequently all stream properties are based upon the total or stagnation enthalpies. The underlying convention is that EBSILON works on basis of total enthalpies.

    The energy balance of turbo machinery can only be closed precisely and physically correct on basis of the total enthalpy values, unless exact knowledge of the flow velocities and inlet and exit swirls exists. Total and static enthalpies of a process stream correlate as follows:

    Figure 2 shows the corresponding expansion line in the enthalpy-entropy diagram, exactly relating static and total properties by the portion of kinetic energy.

    Figure 2: Total and static properties in the enthalpy-Entropy diagram for component 6

    Isentropic efficiencies must be build consistently using either static or total properties.

    Application to Plant Measurements

    In EBSILON applications for performance monitoring, where calculations are correlated with actual plant measurements, it is important to recognize that pressure and temperature measurements are rarely available in form of total pressure and total temperature measurements but rather in terms of static pressure and static temperature measurements. To be precise, these static measurements would need to be converted into stagnation properties, using above relationship. Figure 3 shows a balance, using the exact relationships, with total (stagnation) properties on the streams.

     

    Figure 3: Energy balance around component 6, using correct definition of static vs total properties.

    Yet, exact information about flow velocities at the point of measurement is rarely available. For practical purposes, these static measurements are frequently used directly to evaluate enthalpies for input into EBSILON. Figure 4 shows the same balance, using static measurements in lieu of total properties, feeding into component 6.

     

    Figure 4: Energy balance around component 6, using static measurements.

     

    The impact is small - in particular in relation to the absolute number - in regular applications where Component 6 is used to evaluate stage group efficiencies rather than at stage by stage expansion calculations. Typical design values for flow velocities in steam pipes and steam turbine intake and outlet structures are in the range of 50 to 70 m/s. This is roughly the same for inlet and outlet structures as indicated by significantly larger outlet flanges compared to the diameter of the inlet flanges.

    The impact, whether performance is evaluated based upon static or total values is quite small in numeric terms. The uncertainty of pressure and temperature measurement is most of the time bigger than this ambiguity. Figure 5 shows an Enthalpy-Entropy diagram with expansion lines for both, one based on total properties and one based on static properties in true proportions. This clearly shows that the order of magnitude of the effects is negligibly small.

    Figure 5: Enthalpy-Entropy diagram showing expansion lines based upon static and total properties in true proportions.

    The enthalpy drop in above example, using static properties, evaluates to 350.442 kJ/kg, while the enthalpy drop, using total properties, evaluates to 350.135 kJ/kg a difference of less than 1/10 of a percent.

    Total Isentropic Efficiency in Component 6

    Nevertheless, EBSILON offers the user to account for effect of kinetic energy content by using the specification FSPEC= Total Isentropic Efficiency ,
    CKIN1 and CKIN2' (FSPEC = 1). Figure 6 shows an energy balance around component 6, using the “Total Isentropic Efficiency” method.

    Figure 6: Energy balance around component 6, using static properties and total isentropic efficiency method.

    Particular care must be taken if component 6 was to be used as the last stage before a condenser. At this position steam velocities leaving the expansion stage are in the order of 100 to 260 m/s, in some part load cases even up to sonic velocity, which is in the order of 410 m/s for steam, with associated kinetic energies ranging from 5 to 80 kJ/kg. In such case the kinematic component of the enthalpy cannot be neglected, particularly since its energy cannot be recovered by the turbine. It has to be considered as exhaust loss, which is dissipated to the environment through the condenser.

    If enough information was available from DCS data, the user might want to reverse-engineer the efficiency characteristics of a particular expansion turbine. In such case, the user is advised to switch to component 122 for the task of modeling the condensing section. Component 122 uses industry standards for properly correcting expansion efficiencies for expansion into the wet steam region and to account for exhaust and ventilation losses.  In addition, component 122 would also limit last stage expansion to Mach number equal 1, which actually constitutes a physical limit of steam turbine expansion in the last stage.

     

    Polytropic efficiency

    The FETA switch can be used to switch between isentropic (ETAIN) and polytropic (ETAPN) efficiency. Both (ETAI and ETAP) are calculated as the result value.

    Switching to external presetting of efficiency on logic port 7 with switch FVALETAI was possible only when using isentropic efficiency (FETA=0). Now the external presetting is also applied to the polytropic efficiency (with FETA=1).

    Since the calculation of the polytropic efficiency ETAP is very time-consuming (in some cases with complex material data, this multiplies the total calculation time of the model), the result value ETAP is no longer calculated automatically for each simulation. As was already the case with the compressor (component 24), the FOUTETAP switch has now also been introduced for component 6, with which it is possible to set whether ETAP should be calculated or not.

     

    Logic inlet (Connection point 7) for controlling component properties

    (see also : Editing components  --> Ports)

    To make component properties like efficiencies or heat transfer coefficients (variation quantity) accessible from the outside (for control or reconciliation) it is possible to place the respective value on an auxiliary line as an indexed measured value (specification value FIND). In the component, the same index must then be entered as specification value IPS.

    It is also possible to place this value on a logic line that is directly connected to the component (please see FVALETAI=2, Variation variable: ETAIN, Dimension: Enthalpy). The advantage is that the allocation is graphically visible, and errors (e.g. when copying) are thus avoided.

    The activation of this logic line can also be made conditional on the mode of calculation. This way, this feature can also be used for designs without having to switch manually all the time.
    For this, the flag FVALETAI features the settings

    This option is available for Components 2, 6, 8, 13, 18, 19, 23, 24, and 94.

    In the case of the heat exchangers these setting options are unnecessary as using the logic line for specifying KAN is only possible in off-design calculations anyway. In design mode, KAN is calculated by Ebsilon, of course.

    If connection point 7 is not used for efficiency specification but is available, the efficiency is output there as enthalpy. Depending on FETA, the isentropic efficiency ETAIN or the polytropic efficiency ETAPN is output.

     

    User Input Values

    FP1N

    Flag for the type of pressure definition

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: P1 calculated from P1NSET (by Stodola equation)

    =1: P1-given from outside

    P1NSET

    Inlet pressure (nominal)

    FIDENT

    Component identification
    Like in Parent Profile (Sub profile option only)

    Expression

    =0: No identification
    =1: Identification of the efficiency by specifying the Power output
    =2: Identification of the efficiency by specifying the outlet enthalpy (H2 given)

    FETA

    Efficiency type

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: Isentropic efficiency

    =1: Polytropic efficiency

    FVALETAI

    Validation of the isentropic efficiency

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: ETAIN / ETAPN used without validation

    =1: (Deprecated) IPS used instead of ETAIN / ETAPN(can be validated)

    =2: ETAIN / ETAPN given by enthalpy on control inlet 7 (enthalpy)

    =4: Enthalpy on control inlet 7 used in design, specification value ETAIN / ETAPN in off-design

    =5: Specification value ETAIN /ETAPN used in design, enthalpy on control inlet 7 in off-design

    ETAIN

    Isentropic efficiency (nominal)

    ETAPN

    Polytropic efficiency (nominal)

    IPS

    Index for pseudo measuring point

    ETAMN

    Mechanical efficiency (nominal)

    QLOSSM

    Mechanical loss (constant fraction), restricted to 5% (warning, if QLOSSM > 5%: "QLOSSM too high")  

    DH1LN

    Specific kinetic energy at the inlet (nominal)

    DH2LN

    Specific kinetic energy at the outlet (nominal)

    = outlet loss (nominal)

    FQ

    Flag for power flow at point 5

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: Add power to power of shaft on HP side           Q6=DQ65+Q5

    =1: Distribute power to both shafts          Q6=DQ65-Q5

    FSPECQ

     Flag for type of efficiency calculation

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: Flow given externally, power calculated

    =1: Power given externally, flow calculated

    FSPEC

    Flag for type of efficiency calculation

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: EBSILON standard formalism (CKIN1 and CKIN2 are not used), (VM1N is used)

    =1: Special mode using total isentropic efficiency (CKIN1 and CKIN2 are used)

    Further note for FSPEC=1:  VM1N is determined from M1N and the "nominal total specific Volume", which in turn is calculated from the total pressure and the total enthalpy.

    FCHR

    Flag for the type of characteristic

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: ETAI/ETAIN = f(M1/M1N)

    =1: ETAI/ETAIN = f((P1/P2)/(P1N/P2N))

    =2: ETA/ETAN = f (VM1/VM1N)

    = Deprecated: -1 : Power definition (Q6), identification mode

    = Deprecated: -2 : Enthalpy given (H2), identification mode

    FMODE

    Flag for calculation mode

    =0: GLOBAL

    =1: local off-design

    =-1: local design

    FADAPT

    Flag for adaption polynomial / adaption function

    Like in Parent Profile (Sub profile option only)

    Expression

    =0: not used and not evaluated

    =1: correction factor [ETAI = ETAIN * char line factor *polynomial]

    =2: Replace [ETAI = ETAIN * polynomial]

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

    =-1: correction factor [ETAI = ETAIN * char line factor *function]

    =-2: Replace [ETAI = ETAIN * function]

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

    EADAPT

    Adaption function
    function evalexpr: REAL;
    begin
           evalexpr:= 1.0;
    end;

    A2

    Outlet cross section (used for the calculation of the outlet velocity, not for outlet losses)

    FOUTETAP Calculate result value ETAP only if this value is 1

    P2N                

    Outlet pressure (nominal)

    H1N                

    Total enthalpy at inlet (nominal)

    H2N               

    Total enthalpy at outlet (nominal)

    T1N                

    Inlet temperature (nominal)

    M1N               

    Inlet mass flow (nominal)

    VM1N            

    Inlet volume flow (nominal)

    VM2N            

    Outlet volume flow (nominal)

    The identification values marked in blue are reference values for off-design mode. These values are referred to for the actual off-design values used in the equations respectively. If these identification values are stream data, then these values are often taken from attached pipes or calculated values.

    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


    Characteristics Lines

    Characteristic line 1  efficiency characteristic line

    FCHR=0:       ETAI/ETAIN  =          f(M1/M1N)                          =f(value)

    FCHR=1:       ETAI/ETAIN  =          f((P1/P2)/(P1N/P2N))          =f(value)

    FCHR=2:       ETAI/ETAIN   =         f (VM1/VM1N)                    =f(value)

    Characteristic line 1:  Efficiency

         X-Axis         1          Value                           1st point
                            2          Value                           2nd point
                            .
                            N         value                           last point


         Y-Axis          1          ETAI/ETAIN                1st point
                            2          ETAI/ETAIN                 2nd point

                            .
                            N         ETAI/ETAIN                 last point
     

    Characteristic line 2  Kinetic portions at inlet (characteristic line relevant only in the mode "total isentropic efficiency")

    Characteristic line 2:  kinetic portions at the inlet

         X-Axis         1          VM1/VM1N               1st point
                            2          VM1/VM1N               2nd point
                            .
                            N         VM1/VM1N               last point


         Y-Axis         1          DH1L/DH1LN            1st point
                            2          DH1L/DH1LN            2nd point

                            .
                            N         DH1L/DH1LN            last point
     

    Characteristic line 3  Kinetic portions at the outlet (characteristic line relevant only in the mode "total isentropic efficiency")

    Characteristic line 3:  kinetic portions at the outlet

         X-Axis         1          VM2/VM2N               1st point
                            2          VM2/VM2N               2nd point
                            .
                            N         VM2/VM2N               last point


         Y-Axis          1          DH2L/DH2LN           1st point
                            2          DH2L/DH2LN            2nd point

                            .
                            N         DH2L/DH2LN            last point  
     

    Characteristic line 4: Correction characteristic line for the Stodola-pressure

    Characteristic line 4:  Correction characteristic line for the Stodola-pressure

         X-Axis         1          M1/M1N                 1st point
                            2          M1/M1N                 2nd point
                            .
                            N         M1/M1N                 last point


         Y-Axis         1        P1/P1_Stodola            1st point
                            2         P1/P1_Stodola           2nd point

                            .
                            N         P1/P1_Stodola          last point
     

    This characteristic line defines a correction factor, with which the pressure calculated according to the Stodola's steam turbine expansion path model is multiplied. This makes it simple to adapt the pressure course to the real plant states.


    Physics Used

    Equations

    Design case

    Simulation flag: GLOBAL = design case and FMODE = GLOBAL

     when FP1N = "P1 given from outside"

     

    P1 = P from outside (P inlet)                           (1)           

     
     

     when FP1N = "P1N=P1NSET"

     

    P1 = P1NSET                                                   (1)           

     
       

     

     

    Off-design case

    Simulation flag: GLOBAL = off-design mode or FMODE = local off-design mode

     when FP1N = "P1 given from outside"

     

    P1 = P from outside (P inlet)                          (1)                     

     
     

    when FP1N = "P1N=P1NSET"

     

    according to setting

    P1 calculation according to  setting Ideal Gas or Real Gas ( Modell Settings -->Simulation Calculation )

    and according to the equation according to Stodola. See also: Part-load - Turbine  

     

     
         

     

    All cases

     

    X1  = f (P1,H1)

    S1  = f (P1,H1)

    V1  = f (P1,H1)

    VM1 = M1 * V1

    S2S = S1

    H2S = f (P2,S2S)

    DHS = H1 - H2S

    if GLOBAL= design mode and FMODE = GLOBAL, then

           ETAI = ETAIN

           else {

           if FCHR = 0

                    FAK = M1/M1N

                    ETAI = ETAIN * f (FAK)   

                    from characteristic line

           if FCHR = 1

                    FAK = (P1/P2)/(P1N/P2N)

                    ETAI = ETAIN * f (FAK)   

                    from characteristic line

           if FCHR = 2

                    FAK = VM1/VM1N

                    ETAI = ETAIN * f (FAK)   

                    from characteristic line

                     }

    Outlet losses

    If GLOBAL= design and FMODE = GLOBAL, then

           DH2L = DH2LN

           else {

           FAK = VM2/VM2N

           DH2L = DH2LN*FAK*FAK

                    }

    DH  = DHS * ETAI

    H2  = H1 - DH + DH2L                             

    X2  = f (P2,H2)

    T2  = f (P2,H2)

    M1  = M2 + M3 + M4

    P3  = P2                                                       

    T3  = T2

    H3  = H2                                                     

    Q3  = M3 * H3

    P4  = P2                                                      

    T4  = T2

    H4  = H2                                                    

    Q4  = M4 * H4

    M2  = M1 - M3 - M4                                   

    Q2  = M2 * H2

    If FQ1 = 0,   FAKT =  1

    If FQ1 = 1,   FAKT = -1

    H6 = (M1*(H1-H2)*ETAM + M5*H5*FAKT)/M6   

     

    Total isentropic efficiency

     
     
    Design case

    given:      P1static,         H1total,              (c1N)2/2

                P2static,         ETAN,                   (c2N)2/2

    with:

                            c1        =          velocity at inlet

                            c2        =          velocity at outlet

                            (c1N)2/2           =          DH1LN

                            (c2N)2/2           =          DH2LN

    H1static =H1total - (c1N)2/2

    S1=f(P1static, H1static)

    P1total=f(H1total ,S1)

    H2isentropic_total=f(S1,P2total)    xxx)

              Total isentropic efficiency :

               ETA =             H1total - H2total_____

                                H1total - H2isentrop_total

    ETAN given in design mode

    H2total = H1total - ETAN* (H1total-H2isentrop_total)

    S2=f(P2static, H2total -(c2N)2/2)

    P2total=f(H2total,S2) 

    Iteration: go to xxx)

     
     
    Off-design

    given:      H1total,  (c1N)2/2 - characteristic line,  M1

                            P2static,  ht - characteristic line,

                            (c2N)2/2 - characteristic line

    P1total=f(Stodola)      xx)

    V1total=f(P1total,H1total)

    Iteration: go to xx)

    S1=f(P1total,H1total)

    H2isentrop_total =f(S1,P2total)    xxx)

    ht from characteristic line

    H2total = H1total - ht* (H1total-H2isentrop_total)

    H2static =H2total -(c2N)2/2

    S2=f(P2static, H2static)

    P2total=f(H2total,S2) 

    Iteration: go to xxx)

    H1static =H1total -(c1N)2/2

    P1static=f(S1, H1static)

     

    Identification mode:

    For FCHR = -1 or FCHR = -2, the efficiency is calculated backward. If, for instance, the power is given (enthalpy), then ETAI is calculated. In this case ETAIN and the specified characteristic lines are ignored.


     

    Results

     

     Isentropic efficiency

    ETAI -
     Polytropic efficiency ETAP -
     Mechanical efficiency (including QLOSSM) ETAM -
     Inlet steam content X1 -
     Outlet steam content X2 -
     Kinetic Energy at inlet DH1L kJ/kg
     Outlet loss DH2L kJ/kg
     Total pressure at inlet P1T bar
     Total pressure at outlet P2T bar
     Entropy at inlet S1 kJ/kgK
     Entropy at outlet S2 kJ/kgK
     Inlet volume flow VM1 m3/s
     Outlet volume flow VM2 m3/s
     Outlet speed WAX m/s
     Relative mass flow M1M1N -
     Relative pressure ratio P1P2N -
     Relative volume flow VM1VM1N -
     Used value for nominal isentropic efficiency ETAINR -
     Used value for nominal polytropic efficiency ETAPNR -
     Relative efficiency ETAETAN -
     Efficiency due to characteristic ETACL -
     Result of ADAPT / EADAPT RADAPT -
     For FSPEC=0; virtual P1N for adaption calculations P1NST bar
     Generated mechanical power DQ65 kW
     Generated mechanical power DQ65 kW

     

     

     

     

    Component Displays

    Display Option 1

    Display Option 2

    Display Option 3

    Display Option 4

    Example

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

    See Also