modelIsobaric
Extends from ThermofluidStream.Interfaces.SISOFlow (Base Model with basic flow eqautions for SISO).
Information
Isobaric process (p_in = p_out) suitable for modeling both heating/cooling (open system, steady-flow, systemSpec==Flow)
or periodic heat transfer in a closed cycle system (systemSpec==Cycle), see section Advanced.
The following quantities can be specified either as parameters or as time-varying input signals:
- Temperature difference:
dT = T_out - T_in - Outlet temperature:
T_out - Specific enthalpy difference:
dh = h_out - h_in - Outlet specific enthalpy:
h_out
During simulation, the physical positive direction of the heat flow rate is visualized on the icon layer by a red arrow.
Further assumptions for the isobaric process (open system, steady-flow):
-
Steady-state conditions:
dE_sys/dt = 0, dm_sys/dt = 0 -
No net external forces on the system: The control volume is not subject to acceleration as a rigid body, i.e.
\sum F_external = 0. -
Rigid boundary, no expansion work:
w_exp = 0 -
Negligible kinetic and potential energy changes: Differences between inlet and outlet are neglected, i.e.
g*z_2 + 1/2*c_2^2 ≈ g*z_1 + 1/2*c_1^2. -
Constant composition: No change in species mass fractions across the control volume, i.e.
X_in = X_out.
Advanced settings
Periodic closed cycle system
The isobaric process (p_in = p_out) can also be used to represent periodic heat transfer in a closed cycle system,
which is highlighted on icon layer by p-c (isobaric, closed) compared to p.
In this case Inlet amd Outlet refer to the initial and final state of the isobaric process.
Assumptions for the isobaric process (closed system, cyclic operation):
-
Cycle-averaged steady state: Over multiple consecutive cycles, the system satisfies
dE_sys/dt = 0anddm_sys/dt = 0. -
The mass flow rate
m_flowrepresents the cycle-averaged mass throughput of the system. -
No net external forces on the system: The control volume is not subject to acceleration as a rigid body, i.e.
\sum F_external = 0. -
Isobaric boundary work: The process involves a movable boundary with mean expansion work at constant pressure
p. -
Negligible kinetic and potential energy changes: Differences between inlet and outlet are neglected, i.e.
g*z_2 + 1/2*c_2^2 ≈ g*z_1 + 1/2*c_1^2. -
Constant composition: No change in species mass fractions across the control volume, i.e.
X_in = X_out.
Heat flow rate input/output
The heat flow rate can optionally be an EnergyFlowInput or an
EnergyFlowOutput; see
UsersGuide.EnergyFlow
for details.
The model supports four possible operating modes:
-
If the outlet state is set and the heat flow rate is not an input,
the model is balanced.
In this case, the heat flow rateQ_flowis computed from the given specific enthalpy differencedhand mass flow ratem_flow:Q_flow = m_flow*dh. -
If the outlet state is set and the heat flow rate is an input,
the model is locally overdetermined.
This condition is indicated by a plus symbol on the icon layer. The mass flow rate is then computed from the given heat flow rateQ_flowand specific enthalpy differencedh:m_flow = Q_flow/dh. In this case zero specific enthalpy differencedh = 0can cause the simulation to fail. To obtain an overall balanced system, a locally underdetermined component, such as Sink_free , must be connected. -
If the outlet state is not set and the heat flow rate is an input,
the model is balanced.
In this case, the specific enthalpy differencedhis calculated from the given mass flow ratem_flowand heat flow rateQ_flow:dh = Q_flow/m_flow. In this case zero mass flow ratem_flow = 0can cause the simulation to fail. -
If the outlet state is not set and the heat flow rate is not an input,
the model is not permissible.
In this case, both the heat flow rateQ_flowand the specific enthalpy differencedhare unknown, and the equationQ_flow = m_flow*dhcannot be balanced. This condition is indicated by a prominent red lightning symbol.
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| String | instanceName (from DropOfCommonsPlus) | getInstanceName() | Instance name |
| SI.AbsolutePressure | p_inf | 1e5 | Ambient pressure |
| Boolean | isCycle | systemSpec == ThermofluidStream.Idealized.Types.SystemModel.Cycle | = true, if systemSpec == Cycle; auxiliary variable |
| Boolean | iconIsHeating | iconType == ThermofluidStream.Idealized.Types.Icons.HeatTransfer.Heating | = true, if iconType == Heating; auxiliary variable |
| Medium.SpecificEnergy | eps_q | Modelica.Constants.eps | Regularization specific heat flow |
| SI.MassFlowRate | eps_m_flow | Modelica.Constants.eps | Regularization mass flow rate |
| SI.HeatFlowRate | eps_Q_flow | Modelica.Constants.eps | Regularization heat flow rate |
| Layout | |||
| Boolean | displayInstanceName (from DropOfCommonsPlus) | dropOfCommons.displayInstanceNames | = true, if instance name is displayed |
| Boolean | displayParameters (from DropOfCommonsPlus) | dropOfCommons.displayParameters | = true, if displaying parameters is enabled |
| Advanced | |||
| Boolean | considerInertance (from SISOFlow) | dropOfCommons.considerInertance | =true, if transient momentum (inertance) term is considered; disable only for advanced use |
| Utilities.Units.Inertance | L (from SISOFlow) | dropOfCommons.L | Inertance |
| StateSelect | m_flowStateSelect (from SISOFlow) | StateSelect.default | State selection for mass flow rate |
| Boolean | clip_p_out (from SISOFlow) | = false, if dr_corr=0 (correction of inertial pressure difference) | |
| Medium.AbsolutePressure | p_min (from SISOFlow) | dropOfCommons.p_min | Minimum steady-state output pressure |
| Initialization › Mass flow rate | |||
| InitializationMethods | initM_flow (from SISOFlow) | ThermofluidStream.Utilities.Types.InitializationMethods.none | Initialization method for mass flow rate |
| SI.MassFlowRate | m_flow_0 (from SISOFlow) | 0 | Initial value for mass flow rate |
| Utilities.Units.MassFlowAcceleration | m_acceleration_0 (from SISOFlow) | 0 | Initial value for derivative of mass flow rate |
| Specification | |||
| IconType | iconType | ThermofluidStream.Idealized.Types.Icons.HeatTransfer.Heating | Defines the initial icon prior to simulation |
| OutletSpecification | outletSpec | ThermofluidStream.Idealized.Types.OutletSpecification.Isobaric.TemperatureDifference | Quantity used to define the outlet state |
| ValueSpecification | outletValueSpec | ThermofluidStream.Types.ValueSpecification.Fixed | Specifies whether the quantity is fixed or prescribed |
| SI.TemperatureDifference | dT_fixed | 0 | Fixed temperature difference (dT = T_out - T_in) (OM-Bug) |
| Medium.Temperature | T_out_fixed | Medium.T_default | Fixed outlet temperature |
| SI.SpecificEnthalpy | dh_fixed | 0 | Fixed specific enthalpy difference (dh = h_out - h_in) |
| SI.SpecificEnthalpy | h_out_fixed | Medium.h_default | Fixed outlet specific enthalpy |
| Advanced specification | |||
| HeatFlowSignal | heatFlowSignal | ThermofluidStream.Idealized.Types.EnergyFlowSignalMode.Disabled | Heat flow signal configuration |
| SystemSpecification | systemSpec | ThermofluidStream.Idealized.Types.SystemModel.Flow | Select whether the system is steady-flow (open) or a closed cycle (periodic) |
| Boolean | specifyOutlet | true | = true to specify the outlet state is specified |
| Layout › Display parameters | |||
| Boolean | showOutletSpecification | true | = true to show the fixed outlet specification value (either dT_fixed, T_out_fixed, dh_fixed, h_out_fixed) |
| Boolean | showHeatFlowDirection | true | = true to show the actual heat flow direction |
| Warnings | |||
| AssertionLevel | assertionLevel | AssertionLevel.warning | Assertion level (pressure drop) |
Connectors
| Type | Name | Default | Description |
|---|---|---|---|
| Inlet | inlet (from SISOFlow) | ||
| Outlet | outlet (from SISOFlow) | ||
| Modelica.Blocks.Interfaces.RealInput | outletSpec_prescribed | Prescribed outlet specification [SI-units] | |
| EnergyFlow.Interfaces.EnergyFlowInput | Q_flow_in | Q_flow | Heat flow rate, dircted into the system [W] |
| EnergyFlow.Interfaces.EnergyFlowOutput | Q_flow_out | -Q_flow | Heat flow rate, directed out of the system [W] |
| EnergyFlow.Interfaces.EnergyFlowOutput | P_out | -P | Power (mean net expansion work for systemSpec == Cycle), directed out of the system [W] |
Components
| Type | Name | Default | Description |
|---|---|---|---|
| SI.MassFlowRate | m_flow (from SISOFlow) | inlet.m_flow | Mass flow rate |
| SI.Pressure | dr_corr (from SISOFlow) | Correction of inertial pressure difference | |
| SI.Pressure | dp (from SISOFlow) | Pressure difference | |
| Medium.AbsolutePressure | p_in (from SISOFlow) | Medium.pressure(inlet.state) | Inlet pressure |
| Medium.SpecificEnthalpy | h_in (from SISOFlow) | Medium.specificEnthalpy(inlet.state) | Inlet specific enthalpy |
| Medium.MassFraction[Medium.nXi] | Xi_in (from SISOFlow) | Medium.massFraction(inlet.state) | Inlet mass fractions |
| Medium.AbsolutePressure | p_out (from SISOFlow) | Outlet pressure | |
| Medium.SpecificEnthalpy | h_out (from SISOFlow) | Outlet specific enthalpy | |
| Medium.MassFraction[Medium.nXi] | Xi_out (from SISOFlow) | Outlet mass fractions | |
| Medium.Temperature | T_in | Medium.temperature(inlet.state) | Inlet state temperature |
| Medium.Temperature | T_out | Outlet state temperature | |
| SI.TemperatureDifference | dT | T_out - T_in | Temperature difference (T_out - T_in) |
| Medium.SpecificEnthalpy | dh | h_out - h_in | Difference in specific enthalpy (h_out - h_in) |
| SI.SpecificEnergy | q | Specific heat flow rate | |
| SI.HeatFlowRate | Q_flow | Heat flow rate | |
| Medium.SpecificEnergy | u_in | h_in - p/rho_in | Inlet specific internal energy |
| Medium.SpecificEnergy | u_out | h_out - p/rho_out | Outlet specific internal energy |
| SI.SpecificEnergy | du | u_out - u_in | Difference in specific internal energy(u_out - u_in) |
| Medium.Density | rho_in | Medium.density(inlet.state) | Inlet density |
| Medium.Density | rho_out | Outlet density | |
| SI.SpecificVolume | v_in | 1/rho_in | Inlet specific volume |
| SI.SpecificVolume | v_out | 1/rho_out | Outlet specific volume |
| SI.SpecificEnergy | w_exp | Specific expansion work | |
| SI.SpecificEnergy | w_amb | Specific ambient pressure work | |
| SI.SpecificEnergy | w_exp_net | Net specific expansion work | |
| SI.Power | P | Power (mean net expansion work) | |
| Real | singularityRegime | =+1.0 for dh:=Q_flow/m_flow -> infty, =-1 for m_flow:=Q_flow/dh -> infty, =0.0 else |
Revisions
-
2026-08, by Silvan Keim (silvan.keim@dlr.de):
Improved icon. -
2026, by Raphael Gebhart (raphael.gebhart@dlr.de):
Initial version.