EBSILON®Professional Online Documentation
EBSILON Professional Components / Components - General and Categories / Renewable / Component 113: Linear Focusing Solar Collector
In This Topic
    Component 113: Linear Focusing Solar Collector
    In This Topic

    Component 113: Line focusing solar collector


    Specifications

    Line connections

    1

    Fluid inlet

    2

    Fluid outlet

    3

    Limit input

    4

    Logic output

     

    General       User Input Values       Physics Used       Solar input: QSOLAR       Characteristic Lines       Displays       Example

     

    General

    This component represents a single line-focusing solar collector which can be of parabolic trough or linear Fresnel type. For the most common systems (Eurothrough, LS-2, IST in the case of parabolic troughs, Industrial Solar (formerly Mirroxx) in the case of Fresnel), the data are already stored in the database ( Component properties window "Basic Properties"--> "Load default-values").
    The underlying models calculate the energy balance from direct solar irradiation to usable heat in the heat transfer fluid. Solar and ambient data are provided by a component 117 (Sun). The component includes models with different degrees of detail which can be switched on or off by several flags. For efficiency data, the user has the possibility to define the coefficients in standard formulations, to use an adaptation function or to define data tables for interpolation.

    In parallel to the thermodynamic behaviour, the component includes a model-based approach for pressure loss calculation. Empirical correlations for single and two-phase flow are available.

    For controlling the focusing, a logic inlet (as previously in component 116) has been added in order to avoid a temperature transgression at the outlet. For this, the flag FFOCUS is to be set to 1 and the desired FOCUS value is to be specified as enthalpy at (logic line) pin 3. If required, this can be changed via a controller in such a way that a requested outlet state is achieved.

    In addition, a logic outlet 4 for QEFF has been added.

    Note on the consideration of the wind impact:

    As the impact of wind speed and wind direction individually depends on the respective plant and the ambient landscape, there are no general formulae for calculating this impact. Therefore the tuning factor CORWIND for the wind impact has the default value 1, i.e. no wind impact is considered. If it is known how great the wind impact is, the corresponding tuning factor for the wind impact can be entered manually in CORWIND.

    If it is even known how the wind impact depends on wind speed and wind direction, Ebsilon allows to store this dependency in a Kernel expression EWIND. For EWIND to be used, the flag FWIND has to be switched to 1 (use of EWIND).

    For creating the Kernel expression, the component provides the following variables as KE_Internals:

    The flag FSWIND determines whether the local values from the component or the superordinate values from the sun component are used for VWIND and AWIND.

    The result of the Kernel expression EWIND must be between 0 and 1. This value is then multiplied by the specification value CORWIND and subtracted from 1 in order to obtain the efficiency for the wind impact:

    ETAWIND = 1 – CORWIND * EWIND

    Thus in the case of EWIND=0, there is no wind impact and
    in the case of EWIND=1, the maximum wind impact exists.
    The efficiency then decreases by a fraction CORWIND.

    Please note that for historical reasons CORWIND is interpreted differently in the case of FWIND=0 and FWIND=1.

    In the case of FWIND = 0, CORWIND is the tuning factor for the efficiency (i.e. CORWIND=1 means no correction); in the case of FWIND = 1, CORWIND is the fraction by which the efficiency is maximally decreased (i.e. CORWIND=0 means no correction).

     


     

    User Input Values

    Geometric Parameters

    COLSET

    Name of collector set load

    FTYPE

    Collector type
    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Parabolic Trough

    =1: Linear Fresnel

    LENGTH

    Gross length of the collector module

    AWIDTH

    Gross aperture width of the collector module

    NRATIO

    Optical active portion of aperture: Ratio of active reflective area to gross collector area given by LENGTH*AWIDTH
    (This factors considers spaces and shadowing by the absorber, e.g.)

    LFOCAL

    Focal length of collector (Parabolic Trough) / Height of absorber tube over mirror plane (Linear Fresnel) (used for endless calculation)

    DINNER

    Inner diameter of absorber tube

    ROWDIST

    Distance of the axis of two parallel collector rows (used for shading calculation)

    CDIST

    Distance between two collectors in series
    (used for calculation of end gains if FELOSS=2..4)

    CAZIM

    Collector azimuth angle: Direction of the positive collector axis. Towards north=0°, positive in eastern direction
    (used for calculation of incident and transversal angle if FSPHI=2)

    CSLOP

    Collector axis slope: Angle between collector axis and horizontal plane (used for calculation of incident and transversal angle if FSPHI=2)

    FMODE

    Flag for calculation mode (design / off-design)

    Like in Parent Profile (Sub Profile option only)

    Expression

    =  0: Global
    =  1:   local off-design (i.e. always off-design mode, even when a design calculation has been done globally)
    = -1: local design

    Pressure Loss Parameters

    FDP12N

    Method for calculation of nominal pressure loss

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Given by parameter DP12N

    =1: Model based calculation

    DP12N

    Nominal pressure loss (this value is used if FDP12N=0)

    NNODEP

    Nodes for pressure loss calculation

    FDP12PL

    Method for calculation of part-load pressure loss

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Depending on mass flow

    =1: Depending on mass and volume flow

    =2: Constant at nominal value (calculated according to FDP12N)

    =3: Model based calculation

    =4: Adaptation function EDP12PL

    EDP12PL

    for FDP12PL=4 adaptation function part-load pressure drop in relation to nominal pressure drop.

    KS

    Equivalent sand roughness of absorber pipe inner surface (this value is used if FDP12N=1)

    ZETA

    Pressure loss coefficient for additional pressure losses not covered by the pressure loss model (used if FDP12N=1)

    Performance Parameters

    FOPT0

    Peak optical efficiency (related to the net aperture area LENGTH*AWIDTH*NRATIO at PHIINC=0)

    CLEANI

    Cleanliness of the mirrors as a ratio of actual reflectivity to nominal reflectivity assumed for FOPT0 (standard value is 1 indicating clean mirrors)

    FFOCUS

    Specification of focus state

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: By specification value FOCUS
    =1: By logical input 3

    FOCUS

    Focus state of the collector (0=not focused, 1=focused, linear in between)

    CORSHAD

    Factor to tune the result of the shading model (0= no shading,  1=based on model)

    FELOSS

    Method for calculation of optical end losses and end gains

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Optical end losses not considered (also end gains not considered)

    =1: End losses considered according to model and tuning factor CORELOS

    =2: End losses and end gains at inflow side of collector considered according to model and tuning factors CORELOS and COREGAI are used

    =3: End losses and end gains at outflow side of collector considered according to model and tuning factors CORELOS and COREGAI are used

    =4: End losses and end gains at both sides of collector considered according to model and tuning factors CORELOS and COREGAI are used

    CORELOS

    Factor to tune the optical end losses calculated from the end loss model (1=no correction of model)

    COREGAI

    Factor to tune the optical end gains calculated from the end gain model (1=no correction of model)

    FWIND

    Tuning factor for wind impact (1=no correction)

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Given by factor CORWIND  (Note: Specification values ​​SVWIND and AWIND remain grayed out, even if FSWIND switches over.)

    =1: Adaptation function EWIND

    CORWIND

    Factor to describe wind impact on optical performance

    Standard value = 1 (no influence)  
    or influence determined by EWIND)

    EWIND

    for FWIND=1 adaptation function wind impact:
    Result: 0..., default is 0.  (These specification values ​​are only used in kernel expressions that  the user has to write himself.)

    NNODE

    Number of nodes along the collector axis used for the calculation of the heat losses. The temperature difference along the collector is divided into NNODE intervals. Heat losses are calculated in the center of each section.

     Ambient and Irradiance Parameters

    FSPHI

    Definition of incident and transversal angle

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Given by parameters PHIINC and PHITRAN

    =1: Incident and transversal angle taken from calculation in superior SUN component with index ISUN

    =2: Incident and transversal angle calculated from collector orientation (given by parameters CAZIM and CSLOP) and sun position obtained from superior SUN component with index ISUN

    PHIINC

    Incident angle prescription (this value is used if FSPHI=0)

    PHITRAN

    Transversal angle prescription (this value is used if FSPHI=0)

    FSDNI

    Definition of direct normal irradiance

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Given by parameter DNI

    =1: Taken from superior SUN component with index ISUN

    DNI

    Direct normal irradiance (this value is used if FSDNI=0)

    FSTAMB

    Definition of ambient temperature

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Given by parameter TAMB

    =1: Taken from superior SUN component with index ISUN

    TAMB

    Ambient temperature (this value is used if FSAMB=0)

    FSWIND

    Definition of wind speed and wind direction (FWIND=1)

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Given by parameters VWIND and AWIND

    =1: Taken from superior SUN component with index ISUN

    (Note: FSWIND only specifies whether the wind speed and wind direction are specified individually for the solar collector
                are to be taken over by the "sun". It decides whether these values ​​are used Fwind.
                If FWIND = 0 (default value is CORWIND ", VWIND and AWIND must remain grayed out, even if FSWIND switches.)

    VWIND

    Wind speed (>0, this value is used if FSWIND=0)

    AWIND

    Wind direction (from south to north=0°, positive in east direction, values in the range of 0..360°, this value is used if FSWIND=0)

    ISUN

    Index of reference solar data component

    Note for VWIND und AWIND: Note for VWIND and AWIND: EBSILON does not calculate the wind influence and does not use these default values unless the user
                                                        defines it himself via EWIND.
                                                        Then he can use VWIND and AWIND, including the setting in FSWIND, i.e. H. VWIND and AWIND can also be taken from component 113. 

    Incident angle correction

    FIAM

    Method for calculation of incident angle correction

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Standard polynomial

    =1: Adaptation functions defined in EPHIINC and EPHITRAN

    =2: Table based values given by CIAMINC and CIAMTRAN

    IAMLA

    Coefficient for standard formulation (longitudinal) (A-Term)

    IAML0

    Coefficient for standard formulation (longitudinal) (constant Term)

    IAML1

    Coefficient for standard formulation (longitudinal) (linear Term)

    IAML2

    Coefficient for standard formulation (longitudinal) ( ^2-Term)

    IAML3

    Coefficient for standard formulation (longitudinal) ( ^3-Term)

    IAML4

    Coefficient for standard formulation (longitudinal) ( ^4-Term)

    IAML5

    Coefficient for standard formulation (longitudinal) ( ^5-Term)

    IAMLCOS

    Coefficient for standard formulation (longitudinal) (Cosine-Term)

    IAMT0

    Coefficient for standard formulation (transversal) (constant Term)          (FTYPE=1)

    IAMT1

    Coefficient for standard formulation (transversal) (linear Term)               (FTYPE=1)

    IAMT2

    Coefficient for standard formulation (transversal) ( ^2-Term)                   (FTYPE=1)

    IAMT3

    Coefficient for standard formulation (transversal) ( ^3-Term)                   (FTYPE=1)

    IAMT4

    Coefficient for standard formulation (transversal) ( ^4-Term)                   (FTYPE=1)

    IAMT5

    Coefficient for standard formulation (transversal) ( ^5-Term)                   (FTYPE=1)

    IAMTCOS

    Coefficient for standard formulation (transversal) (Cosine-Term)          (FTYPE=1)

    EPHIINC

    for FIAM=1 adaptation function for incident angle. Result: 0...90°

    EPHITRAN

    for FIAM=1 and FTYPE=1 adaptation function for transversal angle.
    Result: 0...90° or 270°...360°

     Heat loss parameters

    FQLOSS

    Method for heat loss calculation

    Like in Parent Profile (Sub Profile option only)

    Expression

    =0: Standard polynomial

    =1: Adaptation function defined in EQLOSS

    =2: Table based values given by QLOSSA and QLOSSB

    QLOSSA0

    Coefficient for standard formulation heat loss (no DNI dependency) (constant Term in dT)
    (for FQLOSS=0) 

    QLOSSA1

    Coefficient for standard formulation heat loss (no DNI dependency) (linear Term in dT)
    (for FQLOSS=0)

    QLOSSA2

    Coefficient for standard formulation heat loss (no DNI dependency) (^2 Term in dT)
    (for FQLOSS=0)

    QLOSSA3

    Coefficient for standard formulation heat loss (no DNI dependency) (^3 Term in dT)
    (for FQLOSS=0)

    QLOSSA4

    Coefficient for standard formulation heat loss (no DNI dependency) (^4 Term in dT)
    (for FQLOSS=0)

    QLOSSB0

    Coefficient for standard formulation heat loss (DNI dependency) (const. Term in dT)
    (for FQLOSS=0)

    QLOSSB1

    Coefficient for standard formulation heat loss (DNI dependency) (lin. Term in dT)
    (for FQLOSS=0)

    QLOSSB2

    Coefficient for standard formulation heat loss (DNI dependency) (^2 Term in dT)
    (for FQLOSS=0)

    QLOSSC1

    Coefficient for standard formulation heat loss (no DNI dependency) (lin. Term in T)
    (for FQLOSS=0)

    QLOSSC2

    Coefficient for standard formulation heat loss (no DNI dependency) (^2 Term in T)
    (for FQLOSS=0)

    QLOSSC3

    Coefficient for standard formulation heat loss (no DNI dependency) (^3 Term in T)
    (for FQLOSS=0)

    QLOSSC4

    Coefficient for standard formulation heat loss (no DNI dependency) (^4 Term in T)
    (for FQLOSS=0)

    QLOSSD1

    Coefficient for standard formulation heat loss (DNI dependency) (lin. Term in T)
    (for FQLOSS=0)

    QLOSSD2

    Coefficient for standard formulation heat loss (DNI dependency) (^2 Term in T)
    (for FQLOSS=0)

    EQLOSS

    for FQLOSS=1  adaptation function for receiver heat losses. Result: [W/m]

    Nominal parameters

    M1N                  

    Mass flow (nominal)

    VREFN            

    Specific volume at reference point (nominal)

    The parameters marked in blue are reference quantities for the off-design mode. The actual off-design values refer to these quantities in the equations used.

    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

     

     Result values

    RDNI

    Direct normal irradiance used for calculation

    RSHEIGHT

    Sun height angle used for calculation

    RSAZIM

    Sun azimuth angle used for calculation

    RPHIINC

    Incident angle used for calculation

    RPHITRAN

    Transversal angle used for calculation

    ETACOLL

    Collector efficiency QEFF / (RDNI*ANET)

    RFOCUS

    Value used for FOCUS

    QSOLAR

    Solar heat input

    QLOSS

    Heat losses collector

    QEFF

    Heat absorbed by fluid

    QLSOLAR

    Length specific solar heat input

    QLLOSS

    Length specific heat losses

    QLEFF

    Length specific effective heat

    QASOLAR

    Area specific solar heat input QSOLAR/ANET

    QALOSS

    Area specific heat losses QLOSS/ANET

    QAEFF

    Area specific effective heat QEFF/ANET

    ANET

    Net aperture area

    TAVER

    Average collector temperature 0.5*(T1+T2)

    KIA

    Incident angle modifier

    KIAINC

    Incident angle modifier (longitudinal part)

    KIATRAN

    Incident angle modifier (transversal part)

    ETASHAD

    Shading efficiency

    ETAENDL

    Endless efficiency

    ETASPILL

    Spillage efficiency

    DP12

    Pressure loss over collector

    RVWIND

    Wind speed used in calculation

    RAWIND

    Wind direction used in calculation

    RTAMB

    Ambient temperature used in calculation

     

    Physics used

    Equations

    General heat balance

    The heat input into the fluid flow is given by

     M1*(H2-H1) = QEFF .

    This equation is used for both parabolic trough and Linear Fresnel collectors. The effective heat input QEFF depends on the solar heat input QSOLAR and the thermal losses QLOSS.

     QEFF = QSOLAR - QLOSS  

    The solar input QSOLAR is determined by the equation

     QSOLAR = DNI * ANET * FOPT_0 * KIA * FOCUS * ETASHAD * ETAENDL * ETASPILL * ETA_CLEAN

    with the terms:

     DNI                            Direct normal irradiance in W/m**2

     ANET                         Net aperture area ANET=LENGTH*AWIDTH*NRATIO

     FOPT_0                     Peak optical efficiency (parameter FOPT0)

     KIA                            Incident angle correction (cosine losses already included)

     FOCUS                      Focus state of the collector

     ETASHAD                 Factor to include shading losses

     ETAENDL                  Factor to correct end loss effects determined from model

     ETASPILL                  Factor to include optical losses due to wind impact

     ETA_CLEAN              Factor to correct for actual mirror cleanliness ETA_CLEAN=CLEANI

    The heat losses of the collector to the surrounding are calculated based on the length-specific heat loss QLOSS by

     QLOSS = qloss * LENGTH

    The applied methods for the calculation of the terms are provided in the following section.

     

    Important hints for the user: 

    The performance of the collector depends on the product

     FOPT0 * LENGTH * AWIDTH * NRATIO .

    Since performance data found in literature do not have a unique structure please make sure that FOPT0 is always used together with the corresponding reference area. This might be the gross area as given by LENGTH*AWIDTH or the net area which is reduced by factor NRATIO. If the gross area is reported as the reference, NRATIO should be equal to 1 to come to correct results.

    For Linear Fresnel systems AWIDTH is considered as the width of the collector system. Due to the facetted structure, NRATIO is used to define the net aperture area. It is in the choice of the manufacturer if the net aperture area is defined as the area if all facettes are looking to the zenith or the one of the projected area at perpendicular irradiation. For correct results the aperture area definition should be consistent with the peak optical efficiency value and the incident angle correction values KIA.

    For linear Fresnel systems the peak optical efficiency might not be reached at perpendicular irradiance due to the specific optics of these systems. Two options for defining the parameters are open to the user:

     

    Solar input: QSOLAR

    The peak optical efficiency FOPT0 describes the optical efficiency of the collector under the assumptions

    Deviations from this ideal reference point are described by a number of factors that reduce the available heat. These are described in the following.

    Incident angle correction: KIA

    At non-perpendicular incident of the sun additional losses due to shading of collector structure elements, a longer optical path of the reflected sun rays and angle-dependent optical properties of mirrors and absorber tube occur. These optical effects are summarized in the incident angle correction KIA. Note that this factor already includes the cosine losses of the parabolic trough collector in order to allow the same methodology as for linear Fresnel systems.

    The incident angle correction KIA is calculated as

     KIA = KIAINC(RPHIINC) for parabolic trough systems,  

     KIA = KIAINC(RPHIINC) * KIATRAN(RPHITRAN) for linear Fresnel systems
    (with RPHITRAN=abs(PHITRAN)) .

    The user has three options to specify the relations between the angles RPHIINC, RPHITRAN and KIAINC and KIATRAN which are chosen by flag FSPHI.

    The standard function for parabolic troughs is defined as

     KIAINC = ( 1-IAMLA+IAMLA*cos(RPHIINC) ) * (IAMLCOS*cos(RPHIINC) + IAMLO + IAML1*RPHIINC + IAML2*RPHIINC**2 + IAML3*RPHIINC**3 + IAML4*RPHIINC**4 + IAML5*RPHIINC**5 )

    The structure of this function is chosen to be able to represent common formulations from literature. The terms in the first bracket are necessary if a polynomial-based relation for the incident angle correction has to be represented which does not already include the cosine of the incident angle.

    For the Linear Fresnel system the correlations

    KIAINC = IAMLO + IAML1*RPHIINC + IAML2*RPHIINC**2 + IAML3*RPHIINC**3 + IAML4*RPHIINC**4 + IAML5*RPHIINC**5  

    KIATRAN = IAMTO + IAMT1*RPHITRAN + IAMT2*RPHITRAN**2 + IAMT3*RPHITRAN**3 + IAMT4*RPHITRAN**4 + IAMT5*RPHITRAN**5  

    are used. In all cases the results of the functions are limited to a minimum value of 0. When using the adaptation function or table based method always check if the units (deg or rad) of RPHIINC, RPHITRAN fit to the values you specify.

     

    Shading losses: ETASHAD

    If the sun is near the horizon, parallel collector rows shade each other. This effect is considered by means of the term ETASHAD which is calculated based on geometric relation in dependence of the track angle (=transversal angle) of parabolic trough systems.

     ETASHAD=1 - min(1, CORSHAD * max( 0,1- ROWDIST * cos(RPHITRAN) / AWIDTH ) )

    The term min(...) describes the reduction of available energy as a fraction of the energy available if no shading occurs. If the sun is high above the horizon this term equals to 0 and the ETASHAD is 1. The user has the possibility to tune the model based shading effect by a tuning factor CORSHAD.

     

    Optical end loss effects: ETAENDL

    At incident angles <>0 some fraction of the reflected sun beams at the ends of the collector do not hit the absorber tube. This effect is called optical end loss and is a function of the incident angle RPHIINC. If the next collector is arranged in the same axis the lost sun beams of one collector can hit the absorber tube of the next collector. A fraction of the lost heat can be therefore be regained. This effect is called optical end gain. Since a collector can have a corresponding neighbour collector at one side or at both sides, the effective end gains depend on the position of the collector in the field. The user has the possibility to specify at which end of the collector end gains can be obtained. End gains are always less than the end losses. The flag FELOSS defines in which way the end effects are handled:

    The end loss effects ETA_ELOS are calculated based on the equation

     ETA_ELOS =1

    - CORELOS * min(1, kel * LFOCAL/LENGTH * tan(RPHIINC) )

    + COREGAI * max( 0, keg*min(1, kel * LFOCAL/LENGTH * tan(RPHIINC) ) - CDIST/LENGTH )

     

    where the term with CORELOS represents the end losses and the term with COREGAI represents the end gains. The parameters CORELOS and COREGAI are tuning factors with a default value of 1 that correct the calculated effects by a factor. The values kel and keg are used to include the user's choice in flag FELOSS:

    For the calculation based on FELOSS=2, 3 the sun position (azimuth angle SAZIM) is required to determine the sun position relative to the collector. The sun position has to be provided by the sun model with index ISUN. For the other options the sun position is not required to calculate the end losses.

     

    Wind effects on optical performance: ETASPILL

    Under wind loads the collector structure is deformed which reduces the optical efficiency. This effect is represented by the factor ETASPILL. There is no model or standard formulation for the spillage effect included, as this data on this effect are sparse. The user has two possibilities: 

     

    Cleanliness of mirrors: ETA_CLEAN

    The actual cleanliness of the mirrors relative to the ideally clean state can be specified by parameter CLEANI so that ETA_CLEAN=CLEANI.

     

    Heat losses: QLOSS

    Due to a temperature difference between the heat transfer fluid and the ambient air, heat losses occur in the collector. Heat losses are assumed to only depend on the temperature difference. The user has three options to specify the heat loss qloss per unit length of the collector:

    The predefined function is given as

    where  T-Tamb is the temperature difference between fluid and ambient air and hopt is defined as:

    hopt= KIA * FOCUS * ETASHAD * ETAENDL * ETASPILL * ETA_CLEAN .

    The formulation of the predefined function is chosen in a way that common formulations from literature can be represented. The terms with QLOSSAx model a simple heat loss term that only depends on the fluid temperature which is a simplifying assumption. In effect the heat loss is dependent on the temperature on the outside of the absorber tube. This temperature is thus dependent on the actual radial heat flow which goes in a first approximation linear with the effective irradiance RDNI*hopt. The terms with QLOSSBx are added to model this impact. The polynomial for calculating the heat losses has been expanded by such fractions that directly depend on the temperature (QLOSSC1 to QLOSSC4, QLOSSD1, QLOSSD2).

    The following table gives an overview on the units of the parameters.

    QLOSSA0

    W / m

    QLOSSA1

    W / (m K)

    QLOSSA2

    W / (m K**2)

    QLOSSA3

    W / (m K**3)

    QLOSSA4

    W / (m K**4)

    QLOSSB0

    m

    QLOSSB1

    m / K

    QLOSSB2

    m / K**2

    QLOSSC1

    W / (m °C)

    QLOSSC2

    W / (m °C**2)

    QLOSSC3

    W / (m °C**3)

    QLOSSC4

    W / (m °C**4)

    QLOSSD1

    m / °C

    QLOSSD2

    m / °C**2

      

    With the resulting formula two commonly used formulations can be represented. This is illustrated in the following:

     

    Example: Sandia efficiency correlation

    The correlation for collector efficiency suggested early by Sandia National Laboratories reads

     eta[%] = KIA * (73.3-0.007276*dT) - 0.496*dT / DNI - 0.0691 * dT**2 / DNI .

    In order to obtain the corresponding values for the EBSILON®Professional formulation the %-formulation and the sign have to be considered. Further the SANDIA formulation is based on the aperture area whereas the formulation here is length-specific. With an aperture width for the above given LS-2 collector, the parameters are obtained as:

    QLOSSA0, QLOSSA3, QLOSSA4 = 0

    QLOSSA1 = 0.496 W%/(m**2 K) * 5.0 m / 100% = 0.0248 W/ (m K)

    QLOSSA2 = 0.0691 W%/(m**2 K**2) * 5.0 m / 100% = 0.003455 W/ (m K**2)

    QLOSSB0, QLOSSB2 = 0

    QLOSSB1 = 0.007276 %/K * 5.0 m / 100% = 0.0003638

     

    Example: Eurotrough relation

    The correlation for the heat losses of a Eurotrough collector have been published as

     Qloss[W/m**2] = 0.00047 W/(m**2 K**) * dT**2 .

    This area-specific heat loss can simply be transformed into the length-specific heat loss formulation with the net aperture width of 5.45 m (values were obtained based on net area):

    QLOSSA0, QLOSSA1, QLOSSA3, QLOSSA4 = 0

    QLOSSA2 = 0.00047 W/(m**2 K**) * 5.45 m = 0.0025615 W/ (m K**2)

    QLOSSB0, QLOSSB1, QLOSSB2 = 0

     

    Heat loss / Pressure loss (NNODE / NNODEP)

    The heat loss the user can specify the number of sections to be taken as a basis in the calculation (specification value NNODE). This calculation by sections is also possible for the pressure drop, too (specification value NNODEP).

     

    Nominal pressure loss

    The user has the following options for the calculation of the nominal pressure loss:

     

    Part-load pressure loss

    The user has the following options for the calculation of the part-load pressure loss:

     

    Underlying pressure loss model

    To cover single and two-phase flow with one consistent model the two-phase pressure loss model by Friedel (VDI-Wärmeatlas) is used. This model is described in the following. The calculation of pressure drop coefficients is based on single phase Reynolds numbers RE_L and RE_G

    RE_L = MFLUX * DINNER / ETA_L

    RE_G = MFLUX * DINNER / ETA_G

    with the mass flow density MFLUX=M1/(pi/4*DINNER**2) and the dynamic viscosities of liquid (ETA_L) and gas (ETA_G). For two-phase flow, these Reynolds numbers are calculated as if each phase flows alone in the pipe. The single phase pressure loss coefficients are calculated as

    ZETA_L=( 0.86859 * ln( RE_L / (1.964 * ln(RE_L) - 3.8215) ) )**(-2) for RE_L > 1055

    ZETA_L=64 / RE_L for RE_L < 1055

    ZETA_G=( 0.86859 * ln( RE_G / (1.964 * ln(RE_G) - 3.8215) ) )**(-2) for RE_G > 1055

    ZETA_G=64 / RE_G for RE_G < 1055 .

    The term

    DP_L=ZETA_L * MFLUX**2 / (2*DINNER*RHO_L)

    describes the specific pressure loss (Pa/m), as if the whole mass were be present as liquid. A two-phase multiplier R is used to consider the impact of the gas phase

    DP_S = DP_L * R .

    For calculation of R the Weber number

    WE_L = MFLUX**2 * DINNER / RHO_L / SIGMA.

    and the Froude number

    FR_L = MFLUX**2 / ( 9.81 * DINNER * RHO_L**2 )

    of the liquid phase are required. With these numbers the terms

    A = (1-X)**2 + X**2 * (RHO_L * ZETA_G / RHO_G / ZETA_L)

    VV = (RHO_L/RHO_G)**0.8 * (ETA_G/ETA_L)**0.22 * (1-ETA_G/ETA_L)**0.89 * FR_L**(-0.047) * WE_L**(-0.0334)

    XX = 3.43 * X**0.685 * (1-X)**0.24

    R = A + XX * VV

    are calculated where X is the gaseous fraction of the flow (kg/s of gas per kg/s of total flow) and RHO and ETA are the density and the dynamic viscosity of gas and liquid.

    This model can be used in the whole range from single phase liquid flow over two-phase flow to single phase gas flow. This model does not consider pipe roughness.
    Since standard models for two-phase pressure drops in rough pipes are not available the following assumption is made. In the single phase regions the pressure loss in rough pipes is calculated by the correlation

    ZETA_LR = 0.25 * ( log10( KS /(3.7*DINNER) + 5.74 /RE_L**0.9 ) )**(-2)

    ZETA_GR = 0.25 * ( log10( KS /(3.7*DINNER) + 5.74 /RE_G**0.9 ) )**(-2)

    of Swamee and Jain (1976) which approximates the well known implicit equation by von Karman and Nikuaradse. The corresponding pressure losses are

    DP_LR = ZETA_LR * MFLUX**2 / (2*DINNER*RHO_L)

    DP_GR = ZETA_GR * MFLUX**2 / (2*DINNER*RHO_G) .

    Due to a lack of reliable models the pressure drop in rough pipes in the two-phase region is approximated as

    DP_R = DP_LR + (DP_GR-DP_LR) * X

    The final pressure drop is defined as the maximum of the smooth pipe pressure drop DP_S and the rough pipe pressure drop DP_R as

    DP = max ( DP_S, DP_R) .

    This approach has the advantage that the single phase pressure drops are represented very well. A common model for two-phase pressure drop is also available. With the approach defined here the pressure drop is continuous at the boundaries X=0 and X=1.

    Fluid properties (enthalpy, pressure, steam fraction, ...) for pressure loss calculation are taken from the center of the collector.

     

    Additional pressure drop defined by ZETA

    In addition to the pressure loss calculation for the pipe flow the user may enter a pressure loss coefficient ZETA to generate the additional pressure drop

    DPZETA = ZETA * RHO * ( MFLUX / RHO )**2

    In case of two-phase flow the mass averaged mixture density RHO is used. The density is evaluated in the center of the collector (H1+H2)/2 and (P1+P2)/2.

     


     

    Characteristic Lines 

    CIAMINC: Incident Angle Modifier (Longitudinal)

    Correction Factor = f(PHIINC)

     

    CIAMTRAN: Incident Angle Modifier (Transverse)

    Correction Factor = f(PHITRAN)

     

    CQLOSSA: Heat Loss (dT)

    Heat Loss = f(dT)

     

    CQLOSSB: Heat Loss (dT) / DNI

    Heat Loss = f(dT) / DNI

     


    Component Displays

    Display Option 1

    Example

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

     

    See Also