Buildings.Fluid.Storage.BaseClasses
Package with base classes for Buildings.Fluid.Storage
Information
This package contains base classes that are used to construct the models in Buildings.Fluid.Storage.
Extends from Modelica.Icons.BasesPackage (Icon for packages containing base classes).
Package Content
| Name | Description |
|---|---|
 Buoyancy
|
Model to add buoyancy if there is a temperature inversion in the tank |
 IndirectTankHeatExchanger
|
Heat exchanger typically submerged in a fluid with a second fluid circulating through it |
 PartialStratified
|
Partial model of a stratified tank for thermal energy storage |
 ThirdOrderStratifier
|
Model to reduce the numerical dissipation in a tank |
 Examples
|
Examples for BaseClasses models |
Buildings.Fluid.Storage.BaseClasses.Buoyancy
Model to add buoyancy if there is a temperature inversion in the tank
Information
This model outputs a heat flow rate that can be added to fluid volumes in order to emulate buoyancy during a temperature inversion. For simplicity, this model does not compute a buoyancy induced mass flow rate, but rather a heat flow that has the same magnitude as the enthalpy flow associated with the buoyancy induced mass flow rate.
Extends from Modelica.Blocks.Icons.Block (Basic graphical layout of input/output block).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | Modelica.Media.Interfaces.Pa... | Medium model | |
| Volume | V | Volume [m3] | |
| Integer | nSeg | 2 | Number of volume segments |
| Time | tau | tau( ... | Time constant for mixing [s] |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package Medium | Medium model | |
| HeatPort_a | heatPort[nSeg] | Heat input into the volumes |
Modelica definition
model Buoyancy
"Model to add buoyancy if there is a temperature inversion in the tank"
extends Modelica.Blocks.Icons.Block;
replaceable package Medium = Modelica.Media.Interfaces.PartialMedium "Medium model";
parameter Modelica.SIunits.Volume V "Volume";
parameter Integer nSeg(min=2) = 2 "Number of volume segments";
parameter Modelica.SIunits.Time tau(min=0) "Time constant for mixing";
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a[nSeg] heatPort
"Heat input into the volumes";
Modelica.SIunits.HeatFlowRate[nSeg-1] Q_flow
"Heat flow rate from segment i+1 to i";
protected
parameter Medium.ThermodynamicState sta_default = Medium.setState_pTX(T=Medium.T_default,
p=Medium.p_default, X=Medium.X_default[1:Medium.nXi])
"Medium state at default properties";
parameter Modelica.SIunits.Density rho_default=Medium.density(sta_default)
"Density, used to compute fluid mass";
parameter Modelica.SIunits.SpecificHeatCapacity cp_default=Medium.specificHeatCapacityCp(sta_default)
"Specific heat capacity";
parameter Real k(unit="W/K") = V*rho_default*cp_default/tau/nSeg
"Proportionality constant, since we use dT instead of dH";
Modelica.SIunits.TemperatureDifference dT[nSeg-1]
"Temperature difference between adjoining volumes";
equation
for i in 1:nSeg-1 loop
dT[i] = heatPort[i+1].T-heatPort[i].T;
Q_flow[i] = k*noEvent(smooth(1, if dT[i]>0 then dT[i]^2 else 0));
end for;
heatPort[1].Q_flow = -Q_flow[1];
for i in 2:nSeg-1 loop
heatPort[i].Q_flow = -Q_flow[i]+Q_flow[i-1];
end for;
heatPort[nSeg].Q_flow = Q_flow[nSeg-1];
end Buoyancy;
Buildings.Fluid.Storage.BaseClasses.IndirectTankHeatExchanger
Heat exchanger typically submerged in a fluid with a second fluid circulating through it
Information
This model is a heat exchanger with a moving fluid on one side and a stagnant fluid on the other. It is intended for use when a heat exchanger is submerged in a stagnant fluid. For example, the heat exchanger in a storage tank which is part of a solar thermal system.
This component models the fluid in the heat exchanger, convection between the fluid and the heat exchanger, and convection from the heat exchanger to the surrounding fluid.
The model is based on Buildings.Fluid.HeatExchangers.BaseClasses.HACoilInside and Buildings.Fluid.HeatExchangers.BaseClasses.HANaturalCylinder.
The fluid ports are intended to be connected to a circulated heat transfer fluid while the heat port is intended to be connected to a stagnant fluid.
Extends from Buildings.Fluid.Interfaces.TwoPortFlowResistanceParameters (Parameters for flow resistance for models with two ports), Buildings.Fluid.Interfaces.LumpedVolumeDeclarations (Declarations for lumped volumes), Buildings.Fluid.Interfaces.PartialTwoPortInterface (Partial model transporting fluid between two ports without storing mass or energy).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package MediumHex | Modelica.Media.Interfaces.Pa... | Heat transfer fluid flowing through the heat exchanger | |
| replaceable package MediumTan | Modelica.Media.Interfaces.Pa... | Heat transfer fluid inside the tank | |
| replaceable package Medium | PartialMedium | Medium in the component | |
| Integer | nSeg | nSeg( ... | Number of segments in the heat exchanger |
| HeatCapacity | CHex | Capacitance of the heat exchanger [J/K] | |
| Volume | volHexFlu | Volume of heat transfer fluid in the heat exchanger [m3] | |
| Diameter | dExtHex | Exterior diameter of the heat exchanger pipe [m] | |
| Nominal condition | |||
| PressureDifference | dp_nominal | Pressure difference [Pa] | |
| MassFlowRate | m_flow_nominal | Nominal mass flow rate [kg/s] | |
| HeatFlowRate | Q_flow_nominal | Heat transfer at nominal conditions [W] | |
| Temperature | TTan_nominal | Temperature of fluid inside the tank at UA_nominal [K] | |
| Temperature | THex_nominal | Temperature of fluid inside the heat exchanger at UA_nominal [K] | |
| Real | r_nominal | 0.5 | Ratio between coil inside and outside convective heat transfer |
| Flow resistance | |||
| Boolean | computeFlowResistance | true | =true, compute flow resistance. Set to false to assume no friction |
| Boolean | from_dp | false | = true, use m_flow = f(dp) else dp = f(m_flow) |
| Boolean | linearizeFlowResistance | false | = true, use linear relation between m_flow and dp for any flow rate |
| Real | deltaM | 0.1 | Fraction of nominal flow rate where flow transitions to laminar |
| Dynamics | |||
| Equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Dynamics | massDynamics | energyDynamics | Type of mass balance: dynamic (3 initialization options) or steady state |
| Dynamics | energyDynamicsSolid | energyDynamics | Formulation of energy balance for heat exchanger solid mass |
| Real | mSenFac | 1 | Factor for scaling the sensible thermal mass of the volume |
| Initialization | |||
| AbsolutePressure | p_start | Medium.p_default | Start value of pressure [Pa] |
| Temperature | T_start | Medium.T_default | Start value of temperature [K] |
| MassFraction | X_start[Medium.nX] | Medium.X_default | Start value of mass fractions m_i/m [kg/kg] |
| ExtraProperty | C_start[Medium.nC] | fill(0, Medium.nC) | Start value of trace substances |
| ExtraProperty | C_nominal[Medium.nC] | fill(1E-2, Medium.nC) | Nominal value of trace substances. (Set to typical order of magnitude.) |
| Assumptions | |||
| Boolean | allowFlowReversal | true | = false to simplify equations, assuming, but not enforcing, no flow reversal |
| Advanced | |||
| MassFlowRate | m_flow_small | 1E-4*abs(m_flow_nominal) | Small mass flow rate for regularization of zero flow [kg/s] |
| Boolean | homotopyInitialization | true | = true, use homotopy method |
| Diagnostics | |||
| Boolean | show_T | false | = true, if actual temperature at port is computed |
| Modeling detail | |||
| Boolean | hA_flowDependent | true | Set to false to make the convective heat coefficient calculation of the fluid inside the coil independent of mass flow rate |
| Boolean | hA_temperatureDependent | true | Set to false to make the convective heat coefficient calculation of the fluid inside the coil independent of temperature |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package MediumHex | Heat transfer fluid flowing through the heat exchanger | |
| replaceable package MediumTan | Heat transfer fluid inside the tank | |
| replaceable package Medium | Medium in the component | |
| FluidPort_a | port_a | Fluid connector a (positive design flow direction is from port_a to port_b) |
| FluidPort_b | port_b | Fluid connector b (positive design flow direction is from port_a to port_b) |
| HeatPort_a | port[nSeg] | Heat port connected to water inside the tank |
Modelica definition
model IndirectTankHeatExchanger
"Heat exchanger typically submerged in a fluid with a second fluid circulating through it"
replaceable package MediumHex = Modelica.Media.Interfaces.PartialMedium
"Heat transfer fluid flowing through the heat exchanger";
replaceable package MediumTan = Modelica.Media.Interfaces.PartialMedium
"Heat transfer fluid inside the tank";
extends Buildings.Fluid.Interfaces.TwoPortFlowResistanceParameters;
extends Buildings.Fluid.Interfaces.LumpedVolumeDeclarations(
redeclare final package Medium = MediumHex);
extends Buildings.Fluid.Interfaces.PartialTwoPortInterface(
redeclare final package Medium = MediumHex,
final show_T=false);
parameter Integer nSeg(min=2) "Number of segments in the heat exchanger";
parameter Modelica.SIunits.HeatCapacity CHex
"Capacitance of the heat exchanger";
parameter Modelica.SIunits.Volume volHexFlu
"Volume of heat transfer fluid in the heat exchanger";
parameter Modelica.SIunits.HeatFlowRate Q_flow_nominal
"Heat transfer at nominal conditions";
final parameter Modelica.SIunits.ThermalConductance UA_nominal=
abs(Q_flow_nominal/(THex_nominal-TTan_nominal))
"Nominal UA value for the heat exchanger";
parameter Modelica.SIunits.Temperature TTan_nominal
"Temperature of fluid inside the tank at UA_nominal";
parameter Modelica.SIunits.Temperature THex_nominal
"Temperature of fluid inside the heat exchanger at UA_nominal";
parameter Real r_nominal(min=0, max=1)=0.5
"Ratio between coil inside and outside convective heat transfer";
parameter Modelica.SIunits.Diameter dExtHex
"Exterior diameter of the heat exchanger pipe";
parameter Modelica.Fluid.Types.Dynamics energyDynamics=Modelica.Fluid.Types.Dynamics.DynamicFreeInitial
"Formulation of energy balance for heat exchanger internal fluid mass";
parameter Modelica.Fluid.Types.Dynamics energyDynamicsSolid=energyDynamics
"Formulation of energy balance for heat exchanger solid mass";
parameter Modelica.Fluid.Types.Dynamics massDynamics=energyDynamics
"Formulation of mass balance for heat exchanger";
parameter Boolean homotopyInitialization = true "= true, use homotopy method";
parameter Boolean hA_flowDependent = true
"Set to false to make the convective heat coefficient calculation of the fluid inside the coil independent of mass flow rate";
parameter Boolean hA_temperatureDependent = true
"Set to false to make the convective heat coefficient calculation of the fluid inside the coil independent of temperature";
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a port[nSeg]
"Heat port connected to water inside the tank";
Buildings.Fluid.FixedResistances.PressureDrop res(
redeclare final package Medium = MediumHex,
final dp_nominal=dp_nominal,
final m_flow_nominal=m_flow_nominal,
final allowFlowReversal=allowFlowReversal,
final homotopyInitialization=homotopyInitialization,
final show_T=show_T,
final from_dp=from_dp,
final linearized=linearizeFlowResistance)
"Calculates the flow resistance and pressure drop through the heat exchanger";
Buildings.Fluid.MixingVolumes.MixingVolume vol[nSeg](
redeclare each package Medium = MediumHex,
each nPorts=2,
each m_flow_nominal=m_flow_nominal,
each V=volHexFlu/nSeg,
each energyDynamics=energyDynamics,
each massDynamics=massDynamics,
each p_start=p_start,
each T_start=T_start,
each X_start=X_start,
each C_start=C_start,
each C_nominal=C_nominal,
each final prescribedHeatFlowRate=false,
each final allowFlowReversal=allowFlowReversal) "Heat exchanger fluid";
Modelica.Thermal.HeatTransfer.Components.HeatCapacitor cap[nSeg](
each C=CHex/nSeg,
T(each start=T_start,
each fixed=(energyDynamicsSolid == Modelica.Fluid.Types.Dynamics.FixedInitial)),
der_T(
each fixed=(energyDynamicsSolid == Modelica.Fluid.Types.Dynamics.SteadyStateInitial))) if
not energyDynamicsSolid == Modelica.Fluid.Types.Dynamics.SteadyState
"Thermal mass of the heat exchanger";
protected
Buildings.Fluid.Sensors.MassFlowRate senMasFlo(
redeclare package Medium = MediumHex,
allowFlowReversal=allowFlowReversal)
"Mass flow rate of the heat transfer fluid";
Modelica.Thermal.HeatTransfer.Components.Convection htfToHex[nSeg]
"Convection coefficient between the heat transfer fluid and heat exchanger";
Modelica.Thermal.HeatTransfer.Components.Convection HexToTan[nSeg]
"Convection coefficient between the heat exchanger and the surrounding medium";
Modelica.Thermal.HeatTransfer.Sensors.TemperatureSensor temSenHex[nSeg]
"Temperature of the heat transfer fluid";
Modelica.Thermal.HeatTransfer.Sensors.TemperatureSensor temSenWat[nSeg]
"Temperature sensor of the fluid surrounding the heat exchanger";
Modelica.Blocks.Routing.Replicator rep(nout=nSeg)
"Replicates senMasFlo signal from 1 seg to nSeg";
Buildings.Fluid.HeatExchangers.BaseClasses.HACoilInside hAPipIns[nSeg](
each m_flow_nominal=m_flow_nominal,
each hA_nominal=UA_nominal/nSeg*(r_nominal + 1)/r_nominal,
each T_nominal=THex_nominal,
each final flowDependent=hA_flowDependent,
each final temperatureDependent=hA_temperatureDependent)
"Computation of convection coefficients inside the coil";
Buildings.Fluid.HeatExchangers.BaseClasses.HANaturalCylinder hANatCyl[nSeg](
redeclare each final package Medium = Medium,
each final ChaLen=dExtHex,
each final hA_nominal=UA_nominal/nSeg*(1 + r_nominal),
each final TFlu_nominal=TTan_nominal,
each final TSur_nominal=TTan_nominal - (r_nominal/(1 + r_nominal))*(
TTan_nominal - THex_nominal))
"Calculates an hA value for each side of the heat exchanger";
Modelica.Thermal.HeatTransfer.Sensors.TemperatureSensor temSenSur[nSeg]
"Temperature at the external surface of the heat exchanger";
equation
for i in 1:(nSeg - 1) loop
connect(vol[i].ports[2], vol[i + 1].ports[1]);
end for;
connect(rep.u,senMasFlo. m_flow);
connect(port,HexToTan. fluid);
connect(vol[1].ports[1],senMasFlo.port_b);
connect(cap.port,HexToTan.solid);
connect(vol.heatPort,htfToHex. fluid);
connect(htfToHex.solid,HexToTan. solid);
connect(temSenHex.T, hAPipIns.T);
connect(hAPipIns.hA, htfToHex.Gc);
connect(HexToTan.solid,temSenSur. port);
connect(temSenWat.port, port);
connect(temSenSur.T, hANatCyl.TSur);
connect(hANatCyl.TFlu, temSenWat.T);
connect(port_a, senMasFlo.port_a);
connect(vol[nSeg].ports[2], res.port_a);
connect(res.port_b, port_b);
connect(temSenHex.port, vol.heatPort);
connect(rep.y, hAPipIns.m_flow);
connect(hANatCyl.hA,HexToTan. Gc);
end IndirectTankHeatExchanger;
Buildings.Fluid.Storage.BaseClasses.PartialStratified
Partial model of a stratified tank for thermal energy storage
Information
This is a partial model of a stratified storage tank.
See the Buildings.Fluid.Storage.UsersGuide for more information.
Implementation
This model does not include the ports that connect to the fluid from the outside, because these ports cannot be used for the models that contain the Buildings.Fluid.Storage.BaseClasses.ThirdOrderStratifier.
Extends from Buildings.Fluid.Interfaces.PartialTwoPortInterface (Partial model transporting fluid between two ports without storing mass or energy).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | PartialMedium | Medium in the component | |
| Volume | VTan | Tank volume [m3] | |
| Length | hTan | Height of tank (without insulation) [m] | |
| Length | dIns | Thickness of insulation [m] | |
| ThermalConductivity | kIns | 0.04 | Specific heat conductivity of insulation [W/(m.K)] |
| Integer | nSeg | 2 | Number of volume segments |
| Time | tau | 1 | Time constant for mixing [s] |
| Nominal condition | |||
| MassFlowRate | m_flow_nominal | Nominal mass flow rate [kg/s] | |
| Assumptions | |||
| Boolean | allowFlowReversal | true | = false to simplify equations, assuming, but not enforcing, no flow reversal |
| Advanced | |||
| MassFlowRate | m_flow_small | 1E-4*abs(m_flow_nominal) | Small mass flow rate for regularization of zero flow [kg/s] |
| Diagnostics | |||
| Boolean | show_T | false | = true, if actual temperature at port is computed |
| Dynamics | |||
| Equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Formulation of energy balance |
| Dynamics | massDynamics | energyDynamics | Formulation of mass balance |
| Initialization | |||
| AbsolutePressure | p_start | Medium.p_default | Start value of pressure [Pa] |
| Temperature | T_start | Medium.T_default | Start value of temperature [K] |
| MassFraction | X_start[Medium.nX] | Medium.X_default | Start value of mass fractions m_i/m [kg/kg] |
| ExtraProperty | C_start[Medium.nC] | fill(0, Medium.nC) | Start value of trace substances |
Connectors
| Type | Name | Description |
|---|---|---|
| FluidPort_a | port_a | Fluid connector a (positive design flow direction is from port_a to port_b) |
| FluidPort_b | port_b | Fluid connector b (positive design flow direction is from port_a to port_b) |
| output RealOutput | Ql_flow | Heat loss of tank (positive if heat flows from tank to ambient) |
| HeatPort_a | heaPorVol[nSeg] | Heat port that connects to the control volumes of the tank |
| HeatPort_a | heaPorSid | Heat port tank side (outside insulation) |
| HeatPort_a | heaPorTop | Heat port tank top (outside insulation) |
| HeatPort_a | heaPorBot | Heat port tank bottom (outside insulation). Leave unconnected for adiabatic condition |
Modelica definition
model PartialStratified
"Partial model of a stratified tank for thermal energy storage"
extends Buildings.Fluid.Interfaces.PartialTwoPortInterface;
import Modelica.Fluid.Types;
import Modelica.Fluid.Types.Dynamics;
parameter Modelica.SIunits.Volume VTan "Tank volume";
parameter Modelica.SIunits.Length hTan "Height of tank (without insulation)";
parameter Modelica.SIunits.Length dIns "Thickness of insulation";
parameter Modelica.SIunits.ThermalConductivity kIns = 0.04
"Specific heat conductivity of insulation";
parameter Integer nSeg(min=2) = 2 "Number of volume segments";
////////////////////////////////////////////////////////////////////
// Assumptions
parameter Types.Dynamics energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial
"Formulation of energy balance";
parameter Types.Dynamics massDynamics=energyDynamics
"Formulation of mass balance";
// Initialization
parameter Medium.AbsolutePressure p_start = Medium.p_default
"Start value of pressure";
parameter Medium.Temperature T_start=Medium.T_default
"Start value of temperature";
parameter Medium.MassFraction X_start[Medium.nX] = Medium.X_default
"Start value of mass fractions m_i/m";
parameter Medium.ExtraProperty C_start[Medium.nC](
quantity=Medium.extraPropertiesNames)=fill(0, Medium.nC)
"Start value of trace substances";
// Dynamics
parameter Modelica.SIunits.Time tau=1 "Time constant for mixing";
////////////////////////////////////////////////////////////////////
// Connectors
Modelica.Blocks.Interfaces.RealOutput Ql_flow
"Heat loss of tank (positive if heat flows from tank to ambient)";
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a[nSeg] heaPorVol
"Heat port that connects to the control volumes of the tank";
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a heaPorSid
"Heat port tank side (outside insulation)";
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a heaPorTop
"Heat port tank top (outside insulation)";
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a heaPorBot
"Heat port tank bottom (outside insulation). Leave unconnected for adiabatic condition";
// Models
Buildings.Fluid.MixingVolumes.MixingVolume[nSeg] vol(
redeclare each package Medium = Medium,
each energyDynamics=energyDynamics,
each massDynamics=massDynamics,
each p_start=p_start,
each T_start=T_start,
each X_start=X_start,
each C_start=C_start,
each V=VTan/nSeg,
each m_flow_nominal=m_flow_nominal,
each final mSenFac=1,
each final m_flow_small=m_flow_small,
each final allowFlowReversal=allowFlowReversal) "Tank segment";
protected
parameter Medium.ThermodynamicState sta_default = Medium.setState_pTX(
T=Medium.T_default,
p=Medium.p_default,
X=Medium.X_default[1:Medium.nXi]) "Medium state at default properties";
parameter Modelica.SIunits.Length hSeg = hTan / nSeg
"Height of a tank segment";
parameter Modelica.SIunits.Area ATan = VTan/hTan
"Tank cross-sectional area (without insulation)";
parameter Modelica.SIunits.Length rTan = sqrt(ATan/Modelica.Constants.pi)
"Tank diameter (without insulation)";
parameter Modelica.SIunits.ThermalConductance conFluSeg = ATan*Medium.thermalConductivity(sta_default)/hSeg
"Thermal conductance between fluid volumes";
parameter Modelica.SIunits.ThermalConductance conTopSeg = ATan*kIns/dIns
"Thermal conductance from center of top (or bottom) volume through tank insulation at top (or bottom)";
BaseClasses.Buoyancy buo(
redeclare final package Medium = Medium,
final V=VTan,
final nSeg=nSeg,
final tau=tau) "Model to prevent unstable tank stratification";
Modelica.Thermal.HeatTransfer.Components.ThermalConductor[nSeg - 1] conFlu(
each G=conFluSeg) "Thermal conductance in fluid between the segments";
Modelica.Thermal.HeatTransfer.Components.ThermalConductor[nSeg] conWal(
each G=2*Modelica.Constants.pi*kIns*hSeg/Modelica.Math.log((rTan+dIns)/rTan))
"Thermal conductance through tank wall";
Modelica.Thermal.HeatTransfer.Components.ThermalConductor conTop(
G=conTopSeg) "Thermal conductance through tank top";
Modelica.Thermal.HeatTransfer.Components.ThermalConductor conBot(
G=conTopSeg) "Thermal conductance through tank bottom";
Modelica.Thermal.HeatTransfer.Sensors.HeatFlowSensor heaFloTop
"Heat flow at top of tank (outside insulation)";
Modelica.Thermal.HeatTransfer.Sensors.HeatFlowSensor heaFloBot
"Heat flow at bottom of tank (outside insulation)";
Modelica.Thermal.HeatTransfer.Sensors.HeatFlowSensor heaFloSid[nSeg]
"Heat flow at wall of tank (outside insulation)";
Modelica.Blocks.Routing.Multiplex3 mul(
n1=1,
n2=nSeg,
n3=1) "Multiplex to collect heat flow rates";
Modelica.Blocks.Math.Sum sum1(nin=nSeg + 2);
Modelica.Thermal.HeatTransfer.Components.ThermalCollector theCol(m=nSeg)
"Connector to assign multiple heat ports to one heat port";
equation
connect(buo.heatPort, vol.heatPort);
for i in 1:nSeg-1 loop
// heat conduction between fluid nodes
connect(vol[i].heatPort, conFlu[i].port_a);
connect(vol[i+1].heatPort, conFlu[i].port_b);
end for;
connect(vol[1].heatPort, conTop.port_a);
connect(vol.heatPort, conWal.port_a);
connect(conBot.port_a, vol[nSeg].heatPort);
connect(vol.heatPort, heaPorVol);
connect(conWal.port_b, heaFloSid.port_a);
connect(conTop.port_b, heaFloTop.port_a);
connect(conBot.port_b, heaFloBot.port_a);
connect(heaFloTop.port_b, heaPorTop);
connect(heaFloBot.port_b, heaPorBot);
connect(heaFloTop.Q_flow, mul.u1[1]);
connect(heaFloSid.Q_flow, mul.u2);
connect(heaFloBot.Q_flow, mul.u3[1]);
connect(mul.y, sum1.u);
connect(sum1.y, Ql_flow);
connect(heaFloSid.port_b, theCol.port_a);
connect(theCol.port_b, heaPorSid);
end PartialStratified;
Buildings.Fluid.Storage.BaseClasses.ThirdOrderStratifier
Model to reduce the numerical dissipation in a tank
Information
This model reduces the numerical dissipation that is introduced by the standard first-order upwind discretization scheme which is created when connecting fluid volumes in series.
The model is used in conjunction with Buildings.Fluid.Storage.Stratified. It computes a heat flux that needs to be added to each volume of Buildings.Fluid.Storage.Stratified in order to give the results that a third-order upwind discretization scheme (QUICK) would give.
The QUICK method can cause oscillations in the tank temperatures since the high order method introduces numerical dispersion. There are two ways to reduce the oscillations:
-
To use an under-relaxation coefficient
alphawhen adding the heat flux into the volume. -
To use the first-order upwind for
hOut[2]andhOut[nSeg]. Note: Using it requiresnSeg ≥ 4.
Both approaches are implemented in the model.
The model is used by Buildings.Fluid.Storage.StratifiedEnhanced.
Limitations
The model requires at least 4 fluid segments. Hence, set nSeg to 4 or higher.
Extends from Modelica.Blocks.Icons.Block (Basic graphical layout of input/output block).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | Modelica.Media.Interfaces.Pa... | Medium model | |
| MassFlowRate | m_flow_small | m_flow_small( ... | Small mass flow rate for regularization of zero flow [kg/s] |
| Integer | nSeg | nSeg( ... | Number of volume segments |
| Real | alpha | 0.5 | Under-relaxation coefficient (1: QUICK; 0: 1st order upwind) |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package Medium | Medium model | |
| HeatPort_a | heatPort[nSeg] | Heat input into the volumes |
| input RealInput | m_flow | Mass flow rate from port a to port b |
| input RealInput | H_flow[nSeg + 1] | Enthalpy flow between the volumes |
| FluidPort_a | fluidPort[nSeg + 2] | Fluid port, needed to get pressure, temperature and species concentration |
Modelica definition
model ThirdOrderStratifier
"Model to reduce the numerical dissipation in a tank"
extends Modelica.Blocks.Icons.Block;
replaceable package Medium = Modelica.Media.Interfaces.PartialMedium
"Medium model";
parameter Modelica.SIunits.MassFlowRate m_flow_small(min=0)
"Small mass flow rate for regularization of zero flow";
parameter Integer nSeg(min=4) "Number of volume segments";
parameter Real alpha(
min=0,
max=1) = 0.5 "Under-relaxation coefficient (1: QUICK; 0: 1st order upwind)";
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a[nSeg] heatPort
"Heat input into the volumes";
Modelica.Blocks.Interfaces.RealInput m_flow
"Mass flow rate from port a to port b";
Modelica.Blocks.Interfaces.RealInput[nSeg + 1] H_flow
"Enthalpy flow between the volumes";
Modelica.Fluid.Interfaces.FluidPort_a[nSeg + 2] fluidPort(redeclare each
package Medium = Medium)
"Fluid port, needed to get pressure, temperature and species concentration";
protected
Modelica.SIunits.SpecificEnthalpy[nSeg + 1] hOut
"Extended vector with new outlet enthalpies to reduce numerical dissipation (at the boundary between two volumes)";
Modelica.SIunits.SpecificEnthalpy[nSeg + 2] h
"Extended vector with port enthalpies, needed to simplify loop";
Modelica.SIunits.HeatFlowRate Q_flow[nSeg]
"Heat exchange computed using upwind third order discretization scheme";
// Modelica.SIunits.HeatFlowRate Q_flow_upWind
// "Heat exchange computed using upwind third order discretization scheme"; //Used to test the energy conservation
Real sig
"Sign used to implement the third order upwind scheme without triggering a state event";
Real comSig
"Sign used to implement the third order upwind scheme without triggering a state event";
equation
assert(nSeg >= 4,
"Number of segments of the enhanced stratified tank should be no less than 4 (nSeg>=4).");
// assign zero flow conditions at port
fluidPort[:].m_flow = zeros(nSeg + 2);
fluidPort[:].h_outflow = zeros(nSeg + 2);
fluidPort[:].Xi_outflow = zeros(nSeg + 2, Medium.nXi);
fluidPort[:].C_outflow = zeros(nSeg + 2, Medium.nC);
// assign extended enthalpy vectors
for i in 1:nSeg + 2 loop
h[i] = inStream(fluidPort[i].h_outflow);
end for;
// Value that transitions between 0 and 1 as the flow reverses.
sig = Modelica.Fluid.Utilities.regStep(
m_flow,
1,
0,
m_flow_small);
// at surface between port_a and vol1
comSig = 1 - sig;
// at surface between port_a and vol1
hOut[1] = sig*h[1] + comSig*h[2];
// at surface between vol[nSeg] and port_b
hOut[nSeg + 1] = sig*h[nSeg + 1] + comSig*h[nSeg + 2];
// Pros: These two equations can further reduce the temperature overshoot by using the upwind
// Cons: The minimum of nSeg hase to be 4 instead of 2.
hOut[2] = sig*h[2] + comSig*h[3];
// at surface between vol1 and vol2
hOut[nSeg] = sig*h[nSeg] + comSig*h[nSeg + 1];
// at surface between vol[nSeg-1] and vol[nSeg]
for i in 3:nSeg - 1 loop
// at surface between vol[i-1] and vol[i]
// QUICK method
hOut[i] = 0.5*(h[i] + h[i + 1]) - comSig*0.125*(h[i + 2] + h[i] - 2*h[i + 1])
- sig*0.125*(h[i - 1] + h[i + 1] - 2*h[i]);
// hOut[i] = 0.5*(h[i]+h[i+1]); // Central difference method
end for;
for i in 1:nSeg loop
// difference between QUICK and UPWIND; index of H_flow is same as hOut
Q_flow[i] = m_flow*(hOut[i + 1] - hOut[i]) - (H_flow[i + 1] - H_flow[i]);
end for;
// Q_flow_upWind = sum(Q_flow[i] for i in 1:nSeg); //Used to test the energy conservation
for i in 1:nSeg loop
// Add the difference back to the volume as heat flow. An under-relaxation is needed to reduce
// oscillations caused by high order method
heatPort[i].Q_flow = Q_flow[i]*alpha;
end for;
end ThirdOrderStratifier;