EBSILON®Professional Online Documentation
In This Topic
    Component 119: Indirect Storage
    In This Topic

    Component 119: Indirect Storage


    Specifications

    Line Connections

    1

    Fluid inlet

    2

    Fluid outlet

    3

    Heat extraction from fluid to storage medium

     General       User Input Values       Physics Used       Characteristic Lines      Results       Displays       Example

     

    General

    This component provides fast access to the representation of non-steady state processes in Ebsilon. Although by principle, Ebsilon is limited to steady-state solutions because the component equations do not contain a time-dependency and Ebsilon thus solves a steady state equation system, this component allows the representation of non-steady state processes like e.g.

    by adding respective sources and sinks:

    The component "indirect storage" calculates the non-steady state heat exchange of the material with the fluid that flows through and around it, respectively. For this, the structure flown through is represented by a pipe model of similar mass, heat exchanger surface, and material characteristics. The influence of the kind of flow and the geometry of the actual structure on the heat transfer is taken into consideration by the heat transfer coefficient α at the heat transfer surface inside the pipe. The modeling as a pipe has the advantage that a large number of technical applications can be translated to the theory of the heat transfer in pipes. Another advantage of the pipe geometry consists in the cycle-symmetric temperature distribution and heat transfer, so that the calculation of the temperature field is reduced to a cross-section along the direction of flow.

    Based on an initial state of the temperature field, component 119 calculates the change of the temperature field of the pipe due to a change of the determining factors within a defined period of time. The determining factors are defined by the specification values of the component as well as by the state variables of the fluid at the component inlet. The component has to be used in combination with the time series calculation in which the calculation period is defined and by means of which the determining factors for defined points in time can be changed. During a calculation step in the time series calculation, the specification values and state variables of the model remain constant.


     

    User Input Values

    FMODE

    Flag: Calculation mode (design / off-design)

    =0:  Global
    = 1: Local off-design

    FINIT

     

    Flag: Initializing state

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

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

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

    FINST

     

    Flag: Determination of transient calculation modes

    = 0: Transient solution according to time series table
    = 1: Always stationary solution
    = 2: Transient solution as single calculation, time=TIMESING
    = 3: Transient solution as single calculation, use TIMEMAX from model options

    FALGINST

     

    Flag: Determination of transient calculation algorithms
    = 1: Crank-Nicolson-Algorithm
    = 4: Combined numerical and analytical solution

    FSTO

    Flag: Definition of storage geometry

    =0: By LSTO, ASTO and MSTO
    =1: By LSTO, DIAI and THSTO

    LSTO

    Flow length of storage

    DIAI

    Inner diameter of storage pipe

    THSTO

    Thickness of storage pipe

    ASTO

    Heat exchanging area of storage

    MSTO

    Mass of storage

    FVFLUID

    Determination of the fluid volume that defines the fluid mass at a given fluid density

    =0: Fluid volume defined by VFLUID
    =1: Fluid volume corresponds to the volume enclosed by the pipe

    VFLUID

    Flag: Switch for calculation of fluid volumes

    =0: Given from specification VFLUID
    =1: Given from geometry  -> pipe geometry with circular cross section

    FDATA

     

    Flag: Data source and interpolation

    = 1: Constant from specification values RHO, LAM, CP
    = 2: Linear from specification values RHO/DRHO, LAM/DLAM, CP/DCP
    = 3: From characteristic lines

    RHO

    Density of storage at TREF, or constant and independent of TREF in case of FDATA = 1

    DRHO

    Density change per degree defines the gradient of the linear function FDATA=2

    LAM

    Heat conductivity of storage at TREF, or constant and independent of TREF in case of FDATA = 1

    DLAM

    Heat conductivity change per degree defines the gradient of the linear function FDATA=2

    CP

    Spezific heat capacity of storage at TREF, or constant and independent of TREF in case of FDATA = 1

    DCP

    Change of specific heat capacity per degree defines the gradient of the linear function FDATA=2

    TREF

    Reference temperature for RHO, LAM, CP when FDATA=2

    THISO

    Thickness of insulation

    LAMISO

    Heat conductivity of insulation

    TAUADJ

    Correction factor for the time constant of the wall (reduced physical model only)

    LAMADJ

    Multiplication factor to 1/LAMBDA - the walls heat conductivity resistance (reduced physical model only)

    Setting LAMADJ=0 is equivalent to neglecting the walls heat conductivity resistance: either the wall thickness is infinitely small or lambda value is infinitely high

    Setting LAMADJ=1 is equivalent to computing the  walls heat conductivity with the original value of LAMBDA and wall thickness

    LAMADJ<1 leads to decreasing the walls heat conductivity resistance

    LAMADJ>1 leads to increasing the walls heat conductivity resistance

    FSPECM

    Flag: Handling of fluid mass

    = 1: Fluid mass neglectible
    = 2: Fluid mass considered, outlet equal inlet mass flow
    = 3: Fluid mass considered, outlet different from inlet mass flow
    = 4: Mass flows specified, pressure calculated

    FTTI

    Flag: Handling of temperature during time interval

    =0: Actual temperature at the end of time step
    =1: Average temperature for time step interval
    =2: Linear interpolation at each time step

    FTSTEPS

    Flag: Specification of (sub-) time steps

    =1: By specification value TISPEP
    =2: 0.2 of the stable theoretical time increment
    =3: 0.5 of the stable theoretical time increment
    =4: 1.0 of the stable theoretical time increment
    =5: 2.0 of the stable theoretical time increment
    =6: 5.0 of the stable theoretical time increment

    ISUBMAX

    Maximum number of time sub steps for initialization

    IERRMAX

    Maximum allowed error for initializing step

    TISTEP

    Time step

    NFLOW

    Number of points in x-direction (max. 100) in flow direction. When FSPECM=4 NX is fixed to =1

    NRAD

    Number of points in y-direction (max. 30) in radial direction

    TIMESING

    Integration time for single calculation when FINST=2

    FFREQ

     

    Flag: Frequency of transient calculation

    =1: At each iteration step
    =2: At each 2nd iteration step
    =3: At each 4th iteration step
    =4: At each 8th iteration step

    FSTART

     

    Flag: Specification of start temperature

    =1: From specification value TSTART
    =2: From stationary solution

    TSTART

    Start temperature when FSTART=1

    TMIN

    Lower limit for storage temperature

    TMAX

    Upper limit for storage temperature

    FSTAMB

     

    Flag: Definition of ambient temperature

    =0: Definition specification value (TAMB)
    =1: Defined from superior model

    TAMB

    Ambient temperature

    ISUN

    Index for solar parameters

    FALPHI

     

    Flag: Determination of alpha inside

    =0: From constant value ALPHI
    =1: From function EALPHI

    ALPHI

    Heat transfer coefficient alpha inside (wall to fluid) for design case FALPHI=0;
    In part load cases htc is given by: ALPHI/ALPHIN= (M1/M1N)**EX12. Caution! In case of M1=0 the heat transfer coefficient also becomes 0! To avoid errors it´s recommended to use the expression of EALPHI.

    EALPHI

    Function for alpha inside - In the expression, the values of the fluid temperature T_FLUID and the wall temperature T_WALL can be used as arguments (e.g. calculation of the radiant heat transfer)

    FALPHO

    Flag: Determination of alpha outside

    =0: From constant value ALPHO
    =1: From function EALPHO

    ALPHO

    Outer heat transfer coefficient (to ambient) when FALPHO=0

    EALPHO

    Function for alpha outside

    EX12

    Exponent for heat transfer coefficient
    RALPHI = ALPHI*(M1/M1N**EX12N)

    FVOL

    Flag: Part-load pressure drop

    =0: Only depending on mass flow
    =1: Mass flow and density dependent
    =2: Constant (equal nominal value)

    FDPNUM

    Pressure loss handling in the numerical solution

    =0: Using the average fluid pressure between inlet and outlet
    =1: Use of a linear pressure distribution between inlet and outlet, corresponding pressure values in the individual NFLOW fluid elements

    DP12N

    Pressure drop (nominal)

    TAVSTART

    Starting value for average medium temperature

    HAVSTART

    Starting value for average fluid enthalpy

    PAVSTART

    Starting value for average fluid pressure / starting pressure in case of pressure calculation

    FDIR

    Flag: Pipe direction

    =0: Normal
    =1: Reverted, direction of flow is going to be changed by reflecting the numerical grid vertically without changing the connectors.
          (Caution: Reverting is always carried out, when FDIR=1!)

    TIMETOT0

    Total time at start of calculation (Sum of previous time steps)

    M1N

    Mass flow (nominal)

    V1N

    Spezific volume at inlet (nominal)

    TM12N

    Mean temperature to calculate the alpha number (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

    TAVBEG

    Mean caloric temperature of the storage in the beginning of the time step

    TAVEND

    Mean caloric temperature of the storage at the end of the time step

    T2BEG

    Outlet temperature in the beginning of the time step

    T2END

    Outlet temperature at the end of the time step

    RTAMB

    Ambient temperature

    QSTO

    Amount of heat stored during time step, depending on the mode

    FSPECM =1: only the pipe walls are thermal storage
    FSPECM =2, 3 ,4. pipe masses and fluid masses are part of the thermal storage

    QAV

    Mean heat flow through the storage (QSTO/TIMEINT)

    QEND

    Heat flow through storage at the end of time step

    QAVI

    Mean heat flow from fluid to storage

    QENDI

    Heat flow from fluid to storage at the end of time step

    QAVO

    Mean heat flow from storage to ambient

    QENDO

    Heat flow from storage to ambient at the end of time step

    RALPHI

    Calculated heat transfer coefficient fluid-storage

    RALPHO

    Calculated heat transfer coefficient storage-ambient

    RASTO

    Area for heat transfer   

    RMSTO

    Mass of the storage (only pipe walls without fluid)

    RDIAI

    Inner diameter of pipe

    RTHSTO

    Wall thickness of storage pipe

    RDIAO

    Outer diameter of pipe

    RVFLUID

    Overall fluid volume in the storage

    MFLUID

    Overall fluid mass in storage

    RHOFLAV

    Mean fluid density

    PFLAV

    Mean fluid pressure

    HFLAV

    Mean fluid enthalpy

    TFLAV

    Mean fluid temperature

    BIOT

    Biot number

    FOUR

    Fourier number

    TIMEINT

    Integration time of the actual time step

    TIMETOT

    Total time at the end of calculation

    TIMESUB

    Time interval of the sub time steps

    INSFRAC

    Number of transient calculation steps

    ISUB

    Number of sub time steps

    TISUBREC

    Recommended time interval

    PREC

    Precision indicator


     

    Physics used / Equations

    1    Pipe geometry

    FSTO=0: geometry from LSTO, ASTO and MSTO

    In order to handle complex structures, the storage or pipe mass MSTO, the heat exchanging area ASTO of the pipe side facing the fluid  can be set by the user as well as the length of flow LSTO. The governing equations are:

    (1)
    (2)

    whereas Dout and Din mean outer and inner diameters, RHO represents the density of the pipe/storage material. The resulting geometry of the pipes is defined:

    (3)
    (4)

    The fluid containing volume is determined by the geometry of the pipe:

    (5)

    or the specific value of VSTO can be chosen, what shows advantages in case of complex structures, which deviate from the pipe model. This provides a flexible approach for many modeling applications.

    FSTO=1: via LSTO, DIAI and THSTO

    Alternative it is possible to fix the geometry of the storage by defining the following parameters:

    LSTO:      Pipe length
    DIAI:        Inner diameter of the pipe
    THSTO:    Thickness of the pipe wall

     

    2.1  Calculation of the heat transfer in the storage wall

    For a certain conductivity of temperature:

    (6)

    and given distances dX and dY to the neighbour points, the distribution of temperatures Θ in the walls is calculated by solving the differential equation for all discrete elements in the storage:

    (7)

     simplified to two dimensions the yield is:

    (8)

    Changes to the temperatures starting at an initial state always occur due to heat flow through the wall and the insulating surface (see also chapter 3).

    2.1.1  Crank-Nicolson-Algorithm

    For FALGINST=1 (2D-Grid Crank-Nicolson-Algorithm) the equation (8) is discretized and solved numerically. The number of grid elements in X direction (here the flow direction of fluid) can be manipulated by the user input value NFLOW. The number of grid elements in Y direction (here the direction normal to the tube wall) can be manipulated by the user input value  NRAD. The numerical solution results in a 2-dimensional Temperature field in the storage wall. This resulting field can be observed in matrices MXTSTO, RXTSTO (see below).

    2.1.2  Combined analytical and numerical method

    For FALGINST=4 (model using combined analytical and numerical methods) only 1 grid element is used in wall normal direction Y. For each element there is a mean temperature which is searched for. The heat conductivity in the flow direction is neglected.

    The storage temperature after the time step Dt can be computed as

    (8.1)
        (8.2)

    As no temperature gradients occur within a discretisation element, this approach can be integrated directly and we receive the form of the solution for the storage temperature as a function of the time step width shown in Equation (8.2). The two temperatures with the index “sto“ describe the condition of the storage tank before and after the time interval Dt, here Tinf as a driving force describes the steady-state final temperature of the storage wall element, and t is the storage time constant mc/kA

     

    2.2  Calculation of the media temperatures

    Heat flow between the walls and the fluids depends on position and time, what is described mathematically by linking the convective 9a to the conductive terms 9b:

    (9a)
    (9b)

    whereas  represents the temperature of the wall area facing the fluid.  is the temperature of the media and indicates that the direction of the heat flow is orthogonal to the flow direction of the fluid. Heat flow at outer wall will be calculated equivalently.

    Heat fluxes at the walls induce a change of the media temperatures, which also influences the whole 1-D grid of the fluid. Gradients orthogonal to the flow direction and the flow dynamics as well will be neglected in order to keep the model simple and save CPU-time. Calculation of time dependent changes of the temperatures will be done based on mass and energy balances for every volume element.

    (10a)
    (10b)

    The enthalpy of a fluid element is given by , introduced to eqn 10b the result is:

    (11)

    Left side of this equation shows the energy stored in the fluids, while the right side gives the resulting energy flows from the walls trough the boundaries. Mass of a fluid element is defined by:

    (12)

     All property calculations of the media are based on the EBSILON media library.

     

    2.3   Calculation modes according to FSPECM

    FSPECM=1: In case of continuous flow Min = Mout conditions and neglectible fluid mass in comparison to the mass of the wall, left side of eqn 11 is set to zero, what simplifies the following calculation to:

    (13)

    FSPECM=2: In case of continuous flow Min = Mout or a zero mass flow Mfluid = 0, there´s no change of the masses in the duct pipe. But if the second condition concerning the mass ratio of fluid and storage and coupled to this the energy content of the fluid has to be taken into account, it´s necessary to chance eqn 11 to the following:

    (14)

    Example:    Cooling down a stagnant thermo liquid in duct.

    FSPECM=3: Fluid mass inside the pipe is variable, what can be the result of changing density as a function of pressure and temperature. In assumption of a constant duct volume different mass flows at inlet and outlet can occur. To define the thermodynamical state of the fluid in the beginning of the time series the specification values TAVSTART, HAVSTART and PAVSTART can be used.

    Example:    Evaporation from a vessel after a pressure drop.

    FSPECM=4: Inlet and outlet mass flow need to be specified and don´t have to be equal. As a consequence the properties and also the fluid mass in the storage are changeable. The numerical grid is reduced to one point in flow direction (NX=1), different input will be ignored. Initial state of the fluid is given by TAVSTART, HAVSTART and PAVSTART, what also fixes the fluid mass in storage via the relation m = r(p,H) V.

     

    2.4 Calculation of the heat flows

    The calculation of the heat fluxes is done using eqn 11 combined with the settings given by switch FSPEM. In case of neglectible fluid mass the energy balance for both of the fluids appears as following:

    The mass of the storage, whose changes of temperature induces a heat flow (Q0(t)), only consists of the wall elements. The fluids are connected thermally via the following equation:

    (15)

    In the same way the calculation of heat losses to the ambient (blue arrows) is carried out and passed to the result values of QAVO and QENDO. The transient storage behaviour is realized by solving the volume integral of all temperature changes:

    (16)

    Balancing all heat flows (according to eqn 2.13) delivers the stored amount of thermal energy with further integration concerning the time step. The result is passed to the variable QSTO:

    (17)

    This value QSTO divided by the interval of the time step Dt shows the mean heat flow QAV between wall and FLUID. While this heat flux changes time dependent, there´s a variable QEND representing the value at the end of the time step. To get a closed energy balance the amount of heat received from or given to the fluid is stored in the variables QAVI (mean value) and QENDI, which is also passed to line 3. The balance shows the following relations:  

    Tip: Dependent on the magnitude of the temperature gradients and the choice of the time step, it is not possible to get a closed energy balance. If greater miss balances appear, it will be reasonable to reduce the time interval or do a refinement of the time step by adjusting FTSTEPS.

    If simulations are done with respect to the fluid mass (FSPECMXX = 2,3 and 4), balancing becomes different as shown below:

    Left side of eqn 2.8 is not zero anymore, hence fluid mass is considered in this equation. The walls and fluids as well are acting as transient energy storage.

     

    3. Initializing temperature profile in time series calculation

    In the beginning of a time series calculation the temperature profile has to be defined. This is done whenever the switch FINIT ist set to 1. According to switch FSTART there are two possibilities to get initial values for the storage temperature.

    FSTART=1: The whole temperature profile is set to the value TSTART.
    FSTART=2: A first calculation delivers the temperature distribution in the storage as a long
                    term asymptote using the specifications of ambient temperature, thickness of the
                    pipe, heat conductivity LAM, the heat transfer coefficients ALPHI and ALPHO. 

     

    4. Simulation using the time series dialogue and the EBSILON GUI

    The switch FINST generally controls the transient behaviour of the component.

    FINST=0: Intervals of the time dependent calculation steps are given by the time series dialogue, which takes control of the whole calculation. All transient terms including thermal storages in the equations are going to be solved.

    FINST=1: The calculation is always steady state. All fluid data will be passed directly to line 2 without any calculations. 

    FINST=2: An transient single calculation is done with the time step interval given by TIMESING. The simulation is started with the "Simulate" button.

    FINST=3: An transient single calculation is done with the time step interval given by TIMEMAX from the model options panel (Simulation > Transient > Time handling > calculate maximum) 
    The simulation is started with the "Simulate" button.

     

    Figure 2: Internal numerical grid of the storage device

    5. Global Initialization of Transient Components

    All transient components that possess the flag FINIT can be commonly controlled via one global flag.

    For this, the flag FINIT has been expanded by the position GLOBAL: 0.

    If it is set to this value, the control of the transient simulation will be handed over to the global variable “Transient mode“, which can be found under

    Extras \Model Options\Simulation\Transient\ Combo Box "Transient mode".

    This will then pass on the desired mode (first iteration or following iteration) to the components. This can be controlled from the time series dialog by means of the expression “@calcoptions.sim.transientmode“.

     


     

    Characteristic Lines and Matrices

    In order to set up temperature dependent material properties, three characteristic lines are introduced:

    All plots showing the temperature on their X-axes.

    All other "characteristic lines" form a circular buffer. The user doesn´t have to take care of them.

    Corresponding to this "characteristic lines", there are also result curves.

    Specification matrix MXTSTO and result matrix RXTSTO

    The matrix MXTSTO is linked to the output field RXTSTO in the same way as the characteristic curves and result curves mentioned above.

    The distribution of the values in the storage and the fluids is stored in both matrices (default matrix MXTSTO for time step t-1 and result matrix RXTSTO for time step t).

    For the structure of the matrices, see matrices of component 119.


     

    Displays

    Form 1


    Example

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

    See Also