modelPackedBedStorage_L4

Model for thermal energy storage in packed beds with gaseous heat transfer fluid

Extends from TransiEnt.Basics.Icons.PackedBedThermalStorage (Icon for packed-bed thermal storage model).

Information

1. Purpose of model

The model describes a sensible high-temperature packed-bed thermal energy storage unit. It is composed of the packed-bed, fluid in and outlet spaces, and thermal insulation. 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 in the packed bed,
  • fluid pressure loss and
  • transmission heat losses to the environment.

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 body force according to the Darcy-Forchheimer equation is added to the dynamic momentum balance to account for the packed bed flow resistance.
  • Enclosed fluid volumes on the hot and cold side of the storage unit are included using a single control volume. They can be used to model additional pressure and heat loss at the fluid in- and outflow.
  • An additional heat port is added between the packed bed and air in and outlet space. 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.
  • The thermal insulation model is replaceable and can therefore consider static or dynamic heat transfer with various heat paths and insulation layers.
  • 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
  • No gravitational force is taken into account in the dynamic momentum balance, as a horizontal air flow direction has been used at all tested plants so far.
  • Natural convection inside the packed bed is not taken into account
  • The convective heat transfer between packed bed and insulation depends on the lateral packed bed temperature profile. This profile generally requires a dynamic model for the heat transfer coefficient. Nevertheless, several correlations are implemented, but the user should be aware of this.
  • The thermal capacity and flow resistance of the grating used to hold the storage material at in and outlet is not taken into account, but added to the packed bed volume


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

  • The heat transfer fluid momentum balance is dynamic to allow very small mass flows without numerical errors
  • 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 has been validated with two experimental setups of Siemens Gamesa Renewable Energy in Hamburg-Altenwerder (6 MWh_th) and -Bergedorf (130 MWh_th), Germany.

9. References

[1] M. von der Heyde, Abschlussbericht zum Teilprojekt der TUHH im Verbundforschungsprojekt Future Energy Solution (FES), BMWI 03ET6072C, 2021

[2] M. von der Heyde, Electric Thermal Energy Storage based on Packed Beds for Renewable Energy Integration, Dissertation, Hamburg University of Technology, 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.TemperatureT_start_bed_intif linearInit then linspace(T_start_hot, T_start_cold, N_cv) else T_start_bed
SI.HeatCapacityMassSpecificcp_nom(TILMedia.GasFunctions.specificIsobaricHeatCapacity_pTxi(medium_air, simCenter.p_amb_start, T_c, medium_air.xi_default) + TILMedia.GasFunctions.specificIsobaricHeatCapacity_pTxi(medium_air, simCenter.p_amb_start, T_d, medium_air.xi_default))/2
SI.EnthalpyMassSpecifich_nomTILMedia.GasFunctions.specificEnthalpy_pTxi(medium_air, p_nom, T_nom, medium_air.xi_default)Nominal enthalpy
Nominal Values
SI.PowerE_flow_nomm_flow_nom*cp_nom*(T_c - T_d)
SI.Pressurep_nomsimCenter.p_amb_startNominal pressure (for nominal pressure loss calculation)
SI.TemperatureT_nomsimCenter.T_amb_startNominal Temperature (for nominal pressure loss calculation)
SI.MassFlowRatem_flow_nom10Nominal mass flow rate
SI.PressureDelta_p_nom_bed2e3Nominal pressure loss (used for regularization around zero)
SI.TemperatureT_ref273.15Reference Temperature (for nominal storage capacity calculation)
SI.TemperatureT_c273.15 + 600Nominal Charge Temperature (for nominal storage capacity calculation)
SI.TemperatureT_d273.15 + 20Nominal Discharge Temperature (for calculation of nominal rate of charge and discharge)
SI.TemperatureT_stop_c273.15 + 100Outlet Temperature to stop charge operation (for usable capacity ratio calculation)
SI.TemperatureT_stop_d273.15 + 520Outlet Temperature to stop discharge operation (for usable capacity ratio calculation)
Initialisation
BooleanlinearInittrueTrue for linear bed temperature initialisation
SI.TemperatureT_start_bedfill(293.15, N_cv)Packed Bed Initial Temperature
SI.TemperatureT_start_hotsimCenter.T_amb_startHot Side Initial Temperature
SI.TemperatureT_start_coldsimCenter.T_amb_startCold Side Initial Temperature
SI.Pressurep_start_hotsimCenter.p_amb_startInitial hot side pressure
SI.Pressurep_start_coldsimCenter.p_amb_startInitial cold side pressure
SI.MassFlowRatem_flow_start0Initial mass flow rate
Geometry
Reald_v_m0.03Mean volume-equivalent diameter of storage material particles
Realporosity0.4Packed-bed porosity
Realsphericity0.8Packed-bed mean sphericity (see definition in documentation section)
SI.Lengthpipe_thickness0.005Wall thickness at enclosed air spaces
IntegerN_cv10Number of finite control volumes in flow direction
Fundamental Definitions
TILMedia.GasTypes.BaseGasmedium_airsimCenter.airModelMedium of heat transfer fluid

Connectors

TypeNameDefaultDescription
ClaRa.Basics.Interfaces.GasPortInHotAirPort
ClaRa.Basics.Interfaces.GasPortOutColdAirPort

Components

TypeNameDefaultDescription
TransiEnt.SimCentersimCenter
Summarysummary
TransiEnt.Storage.Heat.PackedBedStorage_L4.Basics.PackedBedControlVolume_L4packedBed
Insulation_bedinsulation_bed
Insulation_hotinsulation_hot
Insulation_coldinsulation_cold
ClaRa.Basics.ControlVolumes.GasVolumes.VolumeGas_L2_advancedair_hot
ClaRa.Basics.ControlVolumes.GasVolumes.VolumeGas_L2_advancedair_cold
ClaRa.Components.Adapters.Scalar2VectorHeatPortscalar2VectorHeatPort
ClaRa.Components.Adapters.Scalar2VectorHeatPortscalar2VectorHeatPort1
ClaRa.Basics.ControlVolumes.SolidVolumes.ThinPlateWall_L4thinPlateWall_L4_1
ClaRa.Basics.ControlVolumes.SolidVolumes.ThinPlateWall_L4thinPlateWall_L4_2

Contents

NameDescription
Summary
PackedBedGeometry
ColdAirGeometry
HotAirGeometry
medium_rock
medium_pipe
Insulation_bed
Insulation_hot
Insulation_cold
PressureLossPB
PressureLossHotAir
PressureLossColdAir
ThermalConductivityPB
HeatTransferPB2Wall
HeatTransferPB2Air
HeatTransferAir2Wall
HeatTransferAir2PB

Revisions