modelPackedBedControlVolume_L4

An array of packed bed cells with a single energy equation for particles and fluid and heat connector at in- and outlet
Diagram of PackedBedControlVolume_L4

Extends from ClaRa.Basics.Icons.Volume_L4, ClaRa.Basics.Icons.ComplexityLevel (Displays the complexity level inside model icon ).

Information

1. Purpose of model

The model describes an axial-flow packed bed used for thermal energy storage. Gaseous heat transfer fluid and solid media storage material are generally assumed. It has been validated with air as heat transfer fluid and natural rock as storage material. Main focus is the

  • transient behavior of the temperature field and
  • fluid pressure loss.

2. Level of detail, physical effects considered, and physical insight

  • The packed bed is discretized in one dimension, which is the main fluid flow direction
  • A single energy equation is used for each finite volume, meaning particle and fluid temperature are not distinct, but a mean packed-bed temperature is used. This generally holds for ideal heat transfer between the storage material and heat transfer fluid. Neverthess, the effect of a limited heat transfer between both can be still taken into account by adoption of the effective packed bed thermal conductivity correlation for example using the approach of Vortmeyer (1974).
  • An additional heat port is added at the sides. It can be used to account for other means of heat transport such as conduction, radiation or natural convection besides the advective energy transport into and out of the packed bed.
  • An additional body force according to the Darcy-Forchheimer equation is added to the dynamic momentum balance to account for the packed bed flow resistance.
  • A storage material medium model is required, which has an additional state variable for the specific internal energy in order to account for a temperature variant specific heat capacity.


3. Limits of validity

  • The one-dimensional spatial representation leads to the plug-flow assumption, meaning no lateral temperature and velocity variations are taken into account
  • A horizontal air flow direction is assumed, thus no gravitational force is taken into account in the dynamic momentum balance
  • Natural convection is not taken into account


4. Interfaces

  1. Hot Air Inlet/Outlet
  2. Cold Air Inlet/Outlet


5. Nomenclature

(no remarks)

6. Governing Equations

(no remarks)

7. Remarks for Usage

  • Packed bed correlations are replaceable
  • The momentum balance is dynamic to allow very small mass flows
  • The mean sphericity describes the ratio of the surface of a set of monodisperse spheres with the same number and overall volume as the particle set to the particles set surface. It thus is not solely depended on the particles shape, but also on the particle size distribution.

8. Validation

The model is validated with two experimental setup of Siemens Gamesa Renewable Energy in Hamburg-Altenwerder (6 MWh_th) and -Bergedorf (130 MWh_th), Germany.

9. References

Abschlussbericht zum Teilprojekt der TUHH im Verbundforschungsprojekt Future Energy Solution (FES) (BMWI 03ET6072C) (2021)

Electric Thermal Energy Storage based on Packed Beds for Renewable Energy Integration, Dissertation, Hamburg University of Technology, Michael von der Heyde (2021)

10. Version History

First Version in 04.2020 for the research project Future Energy Solution (FES) by Michael von der Heyde (heyde@tuhh.de)

Parameters

TypeNameDefaultDescription
SI.DensityMassSpecific[geo.N_cv]rho_nomTILMedia.GasFunctions.density_pTxi(medium, p_nom, T_nom, xi_nom)Nominal density
Fundamental Definitions
TILMedia.GasTypes.BaseGasmediumsimCenter.airModelMedium of heat transfer fluid
BooleanfrictionAtInletfalseTrue if pressure loss between first cell and inlet shall be considered
BooleanfrictionAtOutletfalseTrue if pressure loss between last cell and outlet shall be considered
BooleanshowDatafalse|Summary and Visualisation||True, if a data port containing p,T,h,s,m_flow shall be shown, else false
Nominal Values
SI.Pressure[geo.N_cv]p_nom1e5*ones(geo.N_cv)Nominal pressure
SI.Temperature[geo.N_cv]T_nom293.15*ones(geo.N_cv)Nominal temperature
SI.MassFraction[medium.nc - 1]xi_nom{0.01, 0, 0.1, 0, 0.74, 0.13, 0, 0.02, 0}Nominal gas composition
SI.MassFlowRatem_flow_nom100Nominal mass flow
SI.PressureDelta_p_nom1e4Nominal pressure loss
Initialisation
IntegerinitOption0Type of initialisation
SI.Temperature[:]T_start293.15*ones(geo.N_cv)Initial temperature
SI.Pressure[:]p_start1e5*ones(geo.N_cv)Initial pressure
SI.MassFlowRate[geo.N_cv + 1]m_flow_startones(geo.N_cv + 1)*10Initial mass flow rate
SI.MassFraction[medium.nc - 1]xi_start{0.01, 0, 0.1, 0, 0.74, 0.13, 0, 0.02, 0}Initial gas composition
Initialisation › Model Settings
BooleanuseHomotopysimCenter.useHomotopytrue, if homotopy method is used during initialisation
Geometry
Reald_v_m0.02Mean Volume Equivalent Diameter
Realporosity0.4Porosity
Realsphericity0.8Sphericity

Connectors

TypeNameDefaultDescription
ClaRa.Basics.Interfaces.GasPortIninletInlet port
ClaRa.Basics.Interfaces.GasPortOutoutletOutlet port
ClaRa.Basics.Interfaces.HeatPort_a[geo.N_cv]heatExternal
ClaRa.Basics.Interfaces.HeatPort_a[geo.N_cv]heatInlet
ClaRa.Basics.Interfaces.HeatPort_a[geo.N_cv]heatOutlet

Components

TypeNameDefaultDescription
Stringcomplexity (from ComplexityLevel)"??"
ClaRa.SimCentersimCenter
TransiEnt.Basics.Media.SolidWithTemperatureVariantHeatCapacity[geo.N_cv]rock
PressureLosspressureLossPressure loss model
HeatTransferPB2WallheatTransferExternalExternal heat transfer model
HeatTransferPB2AirheatTransferOutletheat transfer model for hot side
HeatTransferPB2AirheatTransferInletheat transfer model for cold side
ThermalConductivitythermalConductivityThermal Conductivity model
TILMedia.Gas_pTfluidInletGas object at inlet port
TILMedia.Gas_pTfluidOutletGas object at outlet port
Summarysummary
TransiEnt.Storage.Heat.PackedBedStorage_L4.Basics.Records.IComPackedBedControlVolume_L4iCom
Geometrygeo
SI.TemperatureT_inletInlet temperature of component
SI.TemperatureT_outletOutlet temperature of component

Contents

NameDescription
Summary
medium_rock
PressureLoss
HeatTransferPB2Wall
HeatTransferPB2Air
ThermalConductivity
Geometry