modelHotWaterStorage_constProp_L4

Temperature and heat flow rate based model of a stratified thermal storage with finite volume discretisation (1=top, n=bottom) and constant fluid properties, electric heating rods can be added

Extends from TransiEnt.Basics.Icons.ThermalStorageBasic (Icon for thermal storage model).

Information

1. Purpose of model

One dimensional fluid storage model with stratification. Intention of the model is to represent a hot water storage in a bigger system with more accurate outflow temperatures compared to a zero dimensional storage model. Several ports on the producer and grid side can be added as well as electric heating rods.

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

L4: Storage is divided in layered volumes. Each volume is ideally stirred. Between the fluid volumes, heat conduction, heat diffusion, heat losses and boyancy are considered.

Heat losses to the ambient are simplified as heat conduction through top, side wall and bottom.

3. Limits of validity

- The storage model includes just a vertical temperature distribution. No horizontal temperature distribution is modeled. Mixing effects due to the velocity of the fluid at inlets an outlets are not modeled.

- Losses to the ambient are modeled as linear dependent from temperature difference (no radiation or convection modeled)

- Thermodynamic properties of fluid are constant (no temperature dependency modeled)

- No pressure losses or levels modeled

- No change of gaseous state modeled

- Geometry is cylindric

- arbitrary fluid port numbers with individual heigths

4. Interfaces

Heat

heatLosses: ambient temperature and the collected heat flow to the ambient through top, side wall and bottom

port: heat port connected to each layer of the storage tank if fluid ports are disabled.

Fluid

waterPortIn_prod: fluid connection from producer, fluid flows to the storage

waterPortOut_prod: fluid connection to producer, fluid flows from the storage

waterPortIn_grid: fluid connection from heating grid, fluid flows to the storage

waterPortOut_grid: fluid connection to heating Grid, fluid flows from the storage

5. Nomenclature

(no remarks)

6. Governing Equations

Energy and mass or volume balance inside every volume segment. Heat losses due to one dimensional thermal conductance through top, bottom and side wall. Thermal conductance between volume segments. Modeled boyancy introducing heat flow from lower to higher segment if the lower segemnt has a higher temperature. Direct fluid connection between the volumes.

7. Remarks for Usage

The allowed minimum number of volume segements is two. The higher the number of segments the higher the number of equations.

If FluidPorts are disabled, the amount of heat that is transferred to each layer needs to be specified via the heat port 'port'.

8. Validation

The model is validated with hot water storage Vitocell 160E. The storage tank has a capacity of 1000 liters and an inner height of 1.88 metres (without insulation). The tank has multiple fluid inflow and outflow connections. The storage is used for climatisation and is installed at TUHH for research purposes. The model is validated against measurements and simulations from Harmsen.

9. References

(no remarks)

10. Version History

Model created by Tobias Ramm (tobias.ramm@tuhh.de) on Fri Mar 20 2015

Model revised by Lisa Andresen (andresen@tuhh.de), Jun 2015

Model revised and redesigned by Pascal Dubucq (dubucq@tuhh.de) on Wed August 24, 2016

- Switched order of control volumes to more intuitive counting (1=bottom, then counting up until N_cv=top)

- Renamed variables to Transient / ClaRa code conventions

- Changed diffusion modeling to second order finite difference approximation

- Removed TIL Media object that was used to calculate *some* of the thermodynamic properties (but not all of them). Now all properties are constant which reduces model complexity. Relevant properties (density, thermal conductivity, heat cacpacity) are mainly temperature dependent and storage operates in fairly narrow temperature range.

Added cost statistics on March 30, 2017

Model expanded to include several producer and grid ports on different heights and electric heating rods by Carsten Bode (c.bode@tuhh.de), Nov 2018

Modification by Anne Hagemeier (anne.hagemeier@umsicht.fraunhofer.de) in August 2021 (Enabled the model to be used without fluid ports)

Parameters

TypeNameDefaultDescription
IntegerN_cv5Number of finite control volumes
SI.TemperatureT_startfill(273.15 + 60, N_cv)Start temperatures of control volumes
SI.TemperatureT_max_ref90 + 273.15Maximum reference temperature (SOC=1)
SI.TemperatureT_min_ref60 + 273.15Minimum reference temperature (SOC=0)
SI.Timetau_buoyancy1Time constant of buoyancy model
SI.Heightdxh/N_cvHeight of finite volumes
SI.MassmV*rho
SI.ThermalConductanceGcat(1, {A_top*U_top + A_wall/N_cv*U_wall}, fill(A_wall/N_cv*U_wall, N_cv - 2), {A_bottom*U_bottom + A_wall/N_cv*U_wall})
SI.Diameterdsqrt(4*V/h/Modelica.Constants.pi)
Integer[N_cv]nPortsBase.getPortCount(N_cv, i_prodIn, i_prodOut, i_gridIn, i_gridOut)
RealG_diffA_top*k/dx
Integer[4,max({n_prodIn, n_prodOut, n_gridIn, n_gridOut})]nPortSubIdxBase.getPortSubIndex(N_cv, i_prodIn, i_prodOut, i_gridIn, i_gridOut)
BooleanuseElHeaterif n_elHeater > 0 then true else falsetrue if electric heater should be included in the model
Fluid Definition
TILMedia.VLEFluidTypes.BaseVLEFluidmediumsimCenter.fluid1Medium to be used
Geometry
BooleanuseFluidPortstrue
SI.VolumeV1e3Volume of tank
Real[N_cv]p_Volumefill(1/N_cv, N_cv)Proportion of the total volume for the parts of the tank
SI.Heighth1Height of tank
Integern_prodIn1Number of inlet ports on producer side
Integern_prodOutn_prodInNumber of outlet ports on producer side
Integern_gridIn1Number of inlet ports on grid side
Integern_gridOutn_gridInNumber of outlet ports on grid side
SI.AreaA_topV/htop tank area
SI.AreaA_bottomA_topbottom tank area
SI.AreaA_wallModelica.Constants.pi*d*houtside tank area
Heating Units
SI.Height[n_prodIn]h_prodInfill(h, n_prodIn)Height of inlet of primary heat source
SI.Height[n_prodOut]h_prodOutfill(0, n_prodOut)Height of outlet of primary heat source
SI.Height[n_gridIn]h_gridInfill(0, n_gridIn)Height of inlet of grid side
SI.Height[n_gridOut]h_gridOutfill(h, n_gridOut)Height of outlet of grid side
Integer[n_prodIn]i_prodInTransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.getIdxFromHeight(N_cv, h, h_prodIn)Index of control volume containing heater(s) (1=top, n=bottom)
Integer[n_prodOut]i_prodOutTransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.getIdxFromHeight(N_cv, h, h_prodOut)Index of control volume containing heater(s) (1=top, n=bottom)
Integer[n_gridIn]i_gridInTransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.getIdxFromHeight(N_cv, h, h_gridIn)Index of control volume containing grid side inflow(s) (1=top, n=bottom)
Integer[n_gridOut]i_gridOutTransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.getIdxFromHeight(N_cv, h, h_gridOut)Index of control volume containing grid side outflow(s) (1=top, n=bottom)
Integer[n_elHeater]i_elHeaterTransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.getIdxFromHeight(N_cv, h, h_elHeater)Index of control volume containing electric heater(s) (1=top, n=bottom)
Thermodynamics
SI.CoefficientOfHeatTransferU_wall0.5Coefficient of heat transfer from wall to ambient
SI.CoefficientOfHeatTransferU_topU_wallCoefficient of heat transfer from top to ambient
SI.CoefficientOfHeatTransferU_bottomU_wallCoefficient of heat transfer from bottom to ambient
SI.ThermalConductivityk0.6Thermal conductivity of fluid in storage
SI.Densityrho1e3Density of fluid in storage
SI.SpecificHeatCapacityc_v4.185e3Heat capacity of fluid in storage
Electric Heater
Integern_elHeater0true if power port should be shown
BooleanusePowerPortfalsetrue if power port should be shown
SI.Height[n_elHeater]h_elHeaterfill(h/2, n_elHeater)Height of electric heater(s)
SI.HeatFlowRate[n_elHeater]Q_flow_n_elHeaterfill(100e3, n_elHeater)Nominal thermal power
SI.Efficiencyeta_elHeater0.95

Connectors

TypeNameDefaultDescription
TransiEnt.Basics.Interfaces.Thermal.FluidPortIn[n_prodIn]waterPortIn_prod
TransiEnt.Basics.Interfaces.Thermal.FluidPortOut[n_prodOut]waterPortOut_prod
TransiEnt.Basics.Interfaces.Thermal.FluidPortIn[n_gridIn]waterPortIn_grid
TransiEnt.Basics.Interfaces.Thermal.FluidPortOut[n_gridOut]waterPortOut_grid
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_aheatPortAmbientHeat port connecting control volumes with ambient temperature to model energy losses (connect ambient temperature)
TransiEnt.Basics.Interfaces.Thermal.HeatFlowRateIn[n_elHeater]Q_flow_setSetpoint for thermal heat, should be negative
ElectricPowerPorteppChoice of power port
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a[N_cv]portHeat port connecting control volumes with ambient temperature to model energy losses (connect ambient temperature)
TransiEnt.Basics.Interfaces.Thermal.HeatFlowRateOutQ_flow_unload
TransiEnt.Basics.Interfaces.Thermal.HeatFlowRateOutQ_flow_load
TransiEnt.Basics.Interfaces.Thermal.HeatFlowRateOutQ_flow_load_total
TransiEnt.Basics.Interfaces.Thermal.HeatFlowRateOutQ_flow_unload_total

Components

TypeNameDefaultDescription
TransiEnt.SimCentersimCenter
TransiEnt.ModelStatisticsmodelStatistics
TransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.IncompressibleFluidVolumecontrolVolume
Base.IncompressibleFluidVolume_noFluidPortscontrolVolume1
Modelica.Thermal.HeatTransfer.Components.ThermalConductorcondAmbientThermal conductance through side wall of storage
TransiEnt.Components.Statistics.Collectors.LocalCollectors.StorageCostcollectStorageCosts
TransiEnt.Producer.Heat.Power2Heat.ElectricBoiler.ElectricBoiler[n_elHeater]electricHeater
TransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.BuoyancybuoyancyModels buoyancy in tank columes
TransiEnt.Storage.Heat.HotWaterStorage_constProp_L4.Base.Diffusiondiffusion
Modelica.Thermal.HeatTransfer.Components.ThermalCollectorthermalCollectorCollects the thermal lossen from top, sidewall and bottom
Modelica.Thermal.HeatTransfer.Sensors.TemperatureSensortemperatureSensor
Modelica.Blocks.Sources.RealExpressionrealExpression3
Modelica.Blocks.Sources.RealExpressionrealExpression2
Modelica.Blocks.Sources.RealExpressionrealExpression4
SI.Temperature[N_cv]TtemperatureSensor.T
SI.TemperatureT_meansum(temperatureSensor.T)/N_cv
RealSOC(T_mean - T_min_ref)/(T_max_ref - T_min_ref)
SI.EnergyE(sum(temperatureSensor.T) - T_min_ref)*c_v*m/N_cv
SI.HeatFlowRateQ_flow_loss2amb-heatPortAmbient.Q_flow
Modelica.Blocks.Sources.RealExpressionrealExpression1

Contents

NameDescription
ProducerCosts_elHeater
CostStatisticsModel
ElectricPowerPort
PowerBoundary