SIMIT 7 FLOWNET Reference Manual
SIMIT 7 FLOWNET Reference Manual
SIMIT 7
FLOWNET Library
Reference manual
s
Edition
January 2013
Siemens offers simulation software to plan, simulate and optimize plants and machines. The simulation- and optimization-
results are only non-binding suggestions for the user. The quality of the simulation and optimizing results depend on the
correctness and the completeness of the input data. Therefore, the input data and the results have to be validated by the user.
Trademarks
SIMIT® is a registered trademark of Siemens AG in Germany and in other countries.
Other names used in this document can be trademarks, the use of which by third-parties for their own purposes could violate
the rights of the owners.
Contents
1 PREFACE 1
1.2 Contents 1
1.3 Symbols 1
2 INTRODUCTION 3
3 FLOWNETS 4
Table of figures
Figure 3-1: Example flownet 4
Figure 3-2: Flownet graph 5
Figure 3-3: Flownet graph for example shown in Figure 3-1 5
Figure 3-4a,b: Minimal flownet graph and minimal flownet 6
Figure 3-5: Incomplete flownet branch 6
Figure 3-6: Error message for incomplete flownet branch 6
Figure 3-7: Error message for an isolated flownet component 7
Figure 3-8: Error message for components with different cycles 7
Figure 3-9: Branch with branch objects 8
Figure 4-1: Connector component types in the basic library 15
Figure 4-2: Topology flownet connector 15
Figure 4-3: Flownet connector on a branch 16
Figure 4-4: Flownet connector as joint 16
Figure 4-5: GENERAL library directory 17
Figure 4-6: Valve with drive 18
Figure 4-7: Application of the Characteristic component on valve 18
Figure 4-8: Parameters for the Valve component type 19
Figure 4-9: Additional parameters for the Valve component type 19
Figure 4-10: Symbolic representation of the valve position 19
Figure 4-11: Representation of the valve position and flow direction in symbol 20
Figure 4-12: Operating window for the valve component type 20
Figure 4-13: Parameters for the StopValve component type 21
Figure 4-14: Pump with drive 22
Figure 4-15: Pump characteristic 22
Figure 4-16: Extended characteristic 23
Figure 4-17: Parameters for the Pump component type 24
Figure 4-18: Additional parameters for the Pump component type 24
Figure 4-19: Representation of the pump operating status in symbol 24
Figure 4-20: Operating window for the Pump component type 24
Figure 4-21: Operating window for the component type Pnode 25
Figure 4-22: Extended operating window for the component type Pnode 25
Figure 4-23: Operating window for the Mnode component type 26
Figure 4-24: Extended operating window for the Mnode component type 27
Figure 4-25: Parameters for the NetParam component type 28
Figure 4-26: Additional parameters for the NetParam component type 29
Figure 4-27: Parameters for the BranchParam component type 29
Figure 4-28: Operating window for the Joint component type 30
Figure 4-29: Parameters for the JointParam component type 31
Figure 4-30: Additional parameters for the JointParam component type 31
Figure 4-31: Operating window for the JointParam component type 31
Figure 4-32: Library directory MEASURE 32
Figure 4-33: Indication of medium flow for the PipeMeasure component type 33
List of tables
Table 3-1: Flownet variables 8
Table 3-2: List of flownet parameters 14
Table 4-1: Signals of the Measure connection type 32
Table 5-1: FLN3 connection type signals 98
Table 5-2: Signals of connection type FLN5 101
Table 5-3: Signals of connection type FLN7 103
Table 5-4: Constants for flownet components 103
Table 5-5: State variables for water/steam 104
Table 5-6: State functions for water/steam 104
1 PREFACE
1.2 Contents
This manual describes the component types contained in the FLOWNET library and the
modelling approach on which the library is based. Section 2 is an introduction explaining the
basic features of the FLOWNET library.
Section 3 describes how piping networks can be simulated using flownets. The modelling
approach for flownets is described in detail.
The components contained in the FLOWNET library are described in section 4 and the
subsequent sections are useful for understanding their function in a flownet.
Section 5 explains how to create your own component types for flownets. The topological
aspects and the data exchange between the components and the flownet solver are
explained in detail. The information on flownets provided in section 3 is required in order to
understand this section.
1.3 Symbols
Particularly important information is highlighted in the text as follows:
NOTE
CAUTION
This means that the system will not respond as described if the specified
precautionary measures are not applied.
WARNING
STOP This means that the system may suffer irreparable damage or that data may be
lost if the relevant precautionary measures are not applied.
2 INTRODUCTION
The FLOWNET library is an extension of SIMIT, which provides component types for
creating simulations of piping networks. By connecting components in this library, a model of
a piping network, a flownet, is created, which simulates the thermodynamic processes in
piping networks. In conjunction with the FLOWNET library, a special solution method can be
applied in SIMIT that calculates the flow, pressure and specific enthalpy in the piping
network simulation.
Although the SIMIT flownets are based on a modelling approach that employs physical
balance equations, the aim is not to use the dynamic process simulations to facilitate the
design of system components or systems, but instead to provide a physically plausible
simulation of the thermodynamic variables in piping networks for virtual commissioning. This
simulation should be easy to create using components on a graphical interface and be stable
even in extreme situations. During the implementation of the component types in the
FLOWNET library, the focus was on simple parameterisation of the components and stable
behaviour in the flownet, rather than a detailed simulation of the physics.
The component types in the FLOWNET library can be used to create flownets for a variety of
media:
• water/steam,
• liquids or
• ideal gases.
The Component Type Editor (CTE) of SIMIT can be used to create your own flownet
components and thus extend your flownet library. The flownet solution method can also be
used by the components via FLOWNET specific connection types.
NOTE
When the simulation starts, it will check whether your SIMIT installation has a
licence for the FLOWNET library. If there are flownet components in the
simulation project, i.e. components that use the solution method for flownets,
then the simulation can only be run if you have a FLOWNET licence.
3 FLOWNETS
A flownet in SIMIT is a connection of flownet components used for simulating
thermodynamic processes in pipe networks. The simulation of flownets is based on a special
solution method, which is parameterised and configured via the flownet components. The
modelling approach described later restricts the flownet to homogenous media, but can be
used for liquids, ideal gases or water in physical condition that is either liquid or steam.
The FLOWNET library provides component types that can be used to configure flownets, i.e.
to model flownets. As usual in SIMIT, the diagram editor is used to model flownets. The
symbol used for the flownet component types, such as valves or pumps, is the same as that
is normally used in piping diagrams. This allows a flownet model, as shown in Figure 3-1, to
be easily constructed in the form of a piping diagram using the component type symbols.
The flownet topology for configuring the flownet solver is derived from the interconnection of
flownet components. While the simulation is running, the flownet solver and components
exchange data: calculated values or flownet parameters.
like the flownet. The graph represents the topology of the flownet on the level of pipelines
and joints. The graph in Figure 3-2 depicts seven nodes K i and nine branches Z i .
The example in Figure 3-1 results in a graph like the one in Figure 3-3 with three branches
and four nodes.
The boundary conditions for the flownet are set via the nodes K 1, K 3 and K 4 . For example,
the pressure on the connections of both tanks is defined at these external nodes. In contrast,
K 2 is an internal node, for which the relevant variables, such as pressure, are calculated
using the flownet solution method.
The solution method for flownets is based on the flow in the branches being dependant on
the pressure using the momentum balance and the balancing of the flow of matter and
enthalpy at the nodes. The state variables for such a system are the flow of mass in the
branches and the pressure and specific enthalpy at the nodes. Other variables, such as
density and temperature in the flownet can be derived according to the medium in question.
If the flownet components in a branch change neither the rate nor the enthalpy of flow in that
branch, the branch is removed from the flownet, i.e. both nodes are merged to form one
node.
The solution method for flownets is, as usual for SIMIT, a cyclical solution method with
equidistant cycles, which extends the standard solution method. Flownet component types
can also be used in addition to other component types, such as those in the basic library.
Specific connections are used for data exchange between flownet components and the
flownet solver. Through these connections, the components receive values calculated by the
flownet solver; these are then used to calculate variables that are sent back to the flownet
solver.
NOTE
It is clear that a flownet must consist of a minimum of one branch with two nodes. Figure
3-4a shows the minimal graph and corresponding minimal flownet. The nodes can be
external (as in Figure 3-4b) or internal.
If, as shown in Figure 3-5, a branch is not closed by two nodes, the error message "Isolated
branch component(s)“ shown in Figure 3-6, stating the components in the branch, appears
on start up of the simulation and the simulation will not run.
A similar error message (Figure 3-7) appears if the flownet only contains an isolated
component. The corresponding graph then consists of only one branch or node and does not
meet the minimum requirements for a flownet.
The flownet components and the flownet solver are also processed cyclically in SIMIT. The
flownet components must be assigned to a cycle. All components of a flownet must be
parameterised with the same cycle. Otherwise the simulation start up will be interrupted with
a message such as that shown in Figure 3-8.
This condition is always met made for incompressible media. For compressible media
no mass is stored in the branch if the branches have "no" volume. The density
decreases for compressible media in the direction of falling pressure, i.e. in the
direction of flow, the density change is negligible for slight throttling, so it is acceptable
to regard the density as constant.
Thus only the pressure on the connection points needs to be considered for branch objects.
The relationships between the pressures on either side of a branch object and the mass flow
in a branch, such as ρ ∆p ~ m 2 , are purely of an algebraic nature.
Momentum balance is applied to each branch k using pressure forces. Friction forces,
acceleration forces and gravity are ignored. With a uniform cross-section Ak on a branch of
length L k the following applies
dm k
Lk
dt
= Ak ∆p k − Ak ∑ ∆p
κ
k ,κ ,
dm k A ∆p k ∆p k ,κ kg
dt
= 10 5 k
Lk
bar −
∑κ
bar m s 2
.
The length and cross-section of a branch are generally unknown so a reasonable estimate
must be made. Taking a cross-section of 0.05 m2 for example and a length of 10 m, results in
a factor
Ak
Α k = 10 5
Lk
of 500 m. This factor is herein after reffered to as momentum factor.
Vi is the volume of the material balance envelope assigned to the node, ρ i is the density of
the medium within the material balance envelope and m κ are the inflows and outflows.
Inflows are positive ( m κ > 0 ), outflows are negative ( m κ < 0 ).
Using the equation of state ρ = ρ ( p , h ) and assuming an isenthalpic change of state for the
pressure in the node, the following applies:
−1
dp i dρ i
= Vi
dt dp i
∑ m κ .
h = const κ
The term
dρi
Vi
dpi
is a measure of the compressibility of the medium. Using the compression modulus K i
−1
dρ i
K i = M i Vi
dp i
the following applies for the mass balance:
dp i K i
dt
=
Mi
∑ m
κ
κ = ci ∑ m
κ
κ .
The default setting for the specific compression modulus ci = K i / M i in the flownet solution
method is the same for all nodes in the flownet. However, various factors result in increased
compressibility c i of gases and steam in contrast to liquids, and this is taken into
consideration for media with higher and lower densities (see section 3.6).
where M i is the mass and hi is the specific enthalpy of the medium within the material
balance envelope.
From
dhi 1
=
dt M i
∑ m (h
κ∈Z i
κ κ − hi ) ,
it follows that by using the difference in enthalpy only the inflows to the nodes i ( κ ∈ Z i ) need
to be summed.
Just like the mass balance factors c i , the default thermal factors mi = 1 / M i in the flownet
solution method are equal for all nodes, but have different values for media with higher and
lower densities.
Pa
p ⋅ 10 2
ρ= bar .
h
RS + 273,15 K
cp
In the flownet solution method, the specific gas constant has a value of 0.287 kJ/kgK
(specific gas constant for dry air).
dhi 1
=
dt M i κ∈Z i
∑
m κ (hκ − hi ) − ci (Ti − TEnv ) .
where T Env is the ambient temperature, T i is the temperature of the node and
ci = α i Ai
is the determining heat transfer factor, which is the product of the heat transfer coefficient α i
and the heat transfer surface A i .
and
The transition between both parameter values can also be set with a linear transfer function:
c L für ρ > ρ L
c = cG für ρ < ρ G
ρ G (ρ L − ρ)
c L + (cG − c L ) else
ρ(ρ L − ρ G )
and
m L für ρ > ρL
m = mG für ρ < ρG
mG − mL
mL + (v − v L ) else
vG − vL
The Topology connector can be used to create a topological flownet connection between two
or more branch components. Figure 4-3a shows two components connected by the CON_A
connector. The connection is functionally identical to the direct connection of both
components via a connecting line, as shown in Figure 4-3b. Figure 4-4a shows three
components connected by the NODE_A connector. This configuration is functionally identical
to the connection using a node shown in Figure 4-4b.
Symbol
Function
The Valve component type is used for simulating a control valve. Depending on the valve
position, the pressure drop across the valve is calculated by
2
∆p − m 2 sec
= 12960
bar kg h
cV2 ρ 3
m
Where
∆p = p B − p A is the pressure drop across the valve in bar,
m is the mass flow rate in kg/s,
ρ is the density of the medium in kg/m³ and
kV is the flow coefficient of the valve in m³/h.
The reference direction is defined for the mass flow m from connection A to connection B,
i.e. it is m > 0 for mass flow in the reference direction. The pressure drop is therefore a
negative value: ∆p < 0 .
At the Position connection, the H position of the valve drive is given as a percentage value,
to which the drive component types of the SIMIT basic library can be set (Figure 4-6), for
example.
The valve position value is limited to 0 ≤ H ≤ 100% .The valve position is mapped onto flow
coefficient cV using the valve characteristic. The following applies
cV 1 1 H
= + 1 − for a linear characteristic,
cV 100 SV SV 100%
2
cV 1 1 H
= + 1 − for a quadratic characteristic,
cV 100 SV SV 100%
H
cV −1
= SV 100% for an equal percentage characteristic.
cV 100
The position ratio
cV 100
SV =
cV 0
is applied here as a quotient consisting of the flow coefficient c V100 for a completely open
valve ( H = 100% ) and the flow coefficient c V0 for a completely closed valve ( H = 0 ).
Any valve characteristic can be created with the Characteristic component type from the
SIMIT standard library. The characteristic component is simply placed upstream of the
valve's Position input (Figure 4-7).
The linear characteristic is then set on the valve, so that the characteristic specified by the
characteristic component is only scaled by the factor (cV 100 − cV 0 ) / cV 0 and displaced by
cV 0 .
Parameter
The valve characteristic and flow coefficient can be adjusted via parameters.
• Characteristic
Set valve characteristics: linear, quadratic, equal-percentage; adjustable online
• Cvs
Flow coefficient c V100 with cV 100 ≥ cV 0 + 10 −6 m³ /h ; adjustable online
• Cv0
Flow coefficient c V0 with cV 0 ≥ 10 −6 m³/h
The parameters, their units and default values are shown in Figure 4-8.
Additional parameters
Additional parameters can be used to visualise the operating states of the components in
symbols (Figure 4-9).
When ShowFlow is set to True the valve position is shown as in Figure 4-10.
If ShowFlowDirection is also set to true, then the valve position and the direction of flow is
shown as in Figure 4-11.
Figure 4-11: Representation of the valve position and flow direction in symbol
Operating window
The pressure drop ∆p and the flow rate m are shown in the operating window (Figure 4-12).
Symbol
Function
The StopValve component type is used for simulating a non-return valve. Depending on the
direction of flow, the pressure drop across the non-return valve is given by
− m 2 sec
2
Parameter
Both flow coefficients for the valve can be set via parameters:
• Cvs
Flow coefficient c V100 with cV 100 ≥ cV 0 + 10 −6 m³ /h
• Cv0
Flow coefficient c V0 with cV 0 ≥ 10 −6 m³/h
The parameters, their units and default values are shown in Figure 4-13.
In order to implement the non-return function, the flow coefficient c V0 must be sufficiently
small.
Symbol
Function
The component type Pump calculates the pressure boost, which depends on the flow and
speed, according to
(
∆p = n 2 ∆p 0 + ∆p * − ∆p 0 ) m 2
für m > 0 .
(m )
* 2
Where
∆p = p B − p A is the pressure boost in bar,
m is the mass flow rate in kg/s,
The quadratic relation between pressure boost and flow defined above is illustrated in Figure
4-15 for the operation of the pump in normal range, i.e. for m > 0 ∆p > 0 .
(
∆p = n 2 ∆p 0 + ∆p * − ∆p 0 ) m 2
for m < 0
(k m )
* 2
when the flow is reversed ( m < 0 ). The throttling effect can be influenced by the flow
coefficient k. Figure 4-16 shows an extended characteristic for various values of k.
1,5
0,5
0
-2 -1,5 -1 -0,5 0 0,5 1 1,5 2
Pressure boost
k=1
-0,5
k=2
k=3
-1
K = 0,5
-1,5
-2
-2,5
-3
Flow
It is clear that the throttling effect is greater at smaller values of the flow coefficient k.
To stabilise a simulation that contains a component of this type, flow reversal ( ∆p < 0 ) is also
calculated using the above equations. For switched off pumps, i.e. with zero speed, this
results in pure throttling for both the flow in and against the pump direction:
∆p *
− ∆p (
m 2
)
for m < 0
∆p =
0
k m *
2
( )
∆p * − ∆p (
m 2
)
for m > 0
0
m * 2
( )
Parameter
The determining variables for the pump characteristic are set with parameters:
• ZeroFlowHead
Zero flow head ∆p 0 , ∆p 0 > ∆p * ; adjustable online
• NominalPressure
Nominal pressure boost ∆p * , ∆p * > 0 ; adjustable online
• NominalMassflow
Nominal mass flow m * , m * > 0 ; adjustable online
The parameters, their units and default values are shown in Figure 4-24.
Additional parameters
The additional parameter Throttling sets a positive flow coefficient k ( k > 0 ) (Figure 4-18).
If the additional parameter Showflow is set to true, then the operating status will be shown in
the pump symbol (Figure 4-19).
Operating window
The pressure boost ∆p and the flow rate m are shown in the operating window (Figure
4-20).
Symbol
Function
The Pnode component type defines values for the pressure p and specific enthalpy h at its
connection A. This type of component forms a boundary for the flownet. When the flownet is
represented as a graph, this corresponds to an (external) node, for which the pressure and
specific enthalpy are predefined.
Operating window
The values for pressure p and specific enthalpy h can be input digitally in the operating
window (Figure 4-21).
Figure 4-22: Extended operating window for the component type Pnode
Symbol
Function
The Mnode component type defined values for the mass flow rate m and specific enthalpy h
at its connection A. This type of component forms a boundary for the flownet. When the
flownet is represented as a graph, this corresponds to an in or outflow through an (internal)
node or branch. An internal node with defined inflow or outflow is added to the flownet for
each component of this type.
Operating window
The values for mass flow rate m and specific enthalpy h can be input digitally in the operating
window (Figure 4-23).
Figure 4-24: Extended operating window for the Mnode component type
Symbol
Function
The NetParam component type is used for parameterisation of the flownet. The component
can be added to any point on any branch of the flownet.
Parameter
The variables for the flownet can be specified as parameters for the components:
• Medium
The following can be selected as medium for the flownet "Water/Steam", "Liquid", or
"Ideal Gas"
• FactorMomentum
Momentum factor for the flow through a branch of the flownet
• sCompressionGas
Specific compression modulus for the "Water/Steam" medium with density
ρ < 500 kg/m³ or the "Ideal Gas" medium
• sCompressionLiquid
Specific compression modulus for the "Water/Steam" medium with density
ρ > 500 kg/m³ or the "Liquid" medium
• FactorThermalGas
Factor for the enthalpy balancing with "Water/Steam" medium with density
ρ < 500 kg/m³ or the "Ideal Gas" medium
• FactorThermalLiquid
Factor for the enthalpy balancing with "Water/Steam" medium with density
ρ > 500 kg/m³ or the "Liquid" medium
• DensityLiquid
Medium density; only applies to the "Liquid" medium
• sHeatCapGas
Specific heat capacity for gas‚ only applies to the "Ideal Gas" medium
• sHeatCapLiquid
Specific heat capacity for liquids; only applies to the "Liquid" medium
• GasConstant
Specific gas constant‚ only applies to the "Ideal Gas" medium
The parameters, their units and default values are shown in Figure 4-25.
Additional parameters
Initial values and specific characteristics for internal nodes of the flownet can be set via
additional parameters:
• PressureInit
Initialisation value for the pressure in the internal node of the flownet
• sEnthalpyInit
Initialisation value for the specific enthalpy in the internal node of the flownet
• SmoothTransition
When this additional parameter is set to true, the variables sCompression and
FactorThermal are set with a density dependant linear transfer function
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the medium in the internal nodes
and the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-26.
Symbol
Function
The BranchParam component type is used for parameterisation of the branches in a flownet.
The components are placed anywhere along the branch to be parameterised.
Parameter
The momentum balance factor can be defined in the branch that is to be parameterised.
• FactorMomentum
Momentum factor for the flow through a branch of the flownet
The default values for the parameter are shown in Figure 4-27.
Symbol
Function
The component type Joint can be used to join three branches connected at A, B and C in
one node. The Joint component adds an internal node to the flownet.
Operating window
The pressure p, specific enthalpy h, density r and the temperature T of the node are shown
in the operating window (Figure 4-28).
Symbol
Function
The component type JointParam can be used to join three branches connected at A, B and
C in one node. The JointParam component adds an internal node to the flownet.
Parameter
The following variables of the node associated with JointParam can be parameterised:
• sCompressionGas
Specific compression modulus for low density media (gas/steam, ρ < 500 kg/m³ )
• sCompressionLiquid
Specific compression modulus for high density media (liquids, ρ > 500 kg/m³ )
• FactorThermalGas
Factor for the enthalpy balancing for low density media (gas/steam, ρ < 500 kg/m³ )
• FactorThermalLiquid
Factor for the enthalpy balancing for high density media (liquids, ρ > 500 kg/m³ )
The parameters, their units and default values are shown in Figure 4-41.
Additional parameters
Initial values and specific characteristics for the node can be set via additional parameters:
• PressureInit
Initialisation value for the pressure in the node
• sEnthalpyInit
Initialisation value for the specific enthalpy in the node
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the medium in the internal nodes
and the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-30.
Operating window
The pressure p, specific enthalpy h, density r and the temperature T of the node are shown
in the operating window (Figure 4-31).
Connections from the Measure connection type are applied to the measuring components.
Descriptions of the individual signals for this type are given in Table 4-1.
Symbol
Function
The PipeMeasure component type creates a measuring point in the pipe. It is inserted with
its connections A and B at the required measuring point in the flownet. The measuring
process for the various variables is not simulated with suitable models, only the variables
calculated by the flownet solver are output.
The measurement variables are output at the Measure connection:
• absolute flow rate value m ,
Figure 4-33: Indication of medium flow for the PipeMeasure component type
Operating window
The variables output via the Measure connection are also displayed in the operating window
(Figure 4-34).
Values for the specific enthalpy and the density of the medium at the measuring point are
displayed in the extended operating window (Figure 4-35).
Figure 4-35: Extended operating window for the PipeMeasure component type
Symbol
Function
The MeasureAll component type provides the bundled values measured via the Measure
input as individual signals at its outputs: pressure, temperature, flow, level and weight.
This type of component can be connected to the pipe measuring point, for example and
thereby output measuring variables as individual signals (Figure 4-36).
Operating window
The output measuring variables are also displayed in the operating window for the
component type (Figure 4-37).
Symbol
Function
The FlowIndicator component type displays the defined flow at its Measure input.
Operating window
The flow rate value obtained via the Measure input is displayed in the operating window
(Figure 4-38).
Symbol
Function
The LevelIndicator component type displays the defined level at its Measure input.
Operating window
The level value obtained via the Measure input is displayed in the operating window (Figure
4-39).
Symbol
Function
The PressureIndicator component type displays the defined pressure at its Measure input.
Operating window
The pressure value obtained via the Measure input is displayed in the operating window
(Figure 4-40).
Symbol
Function
The TemperatureIndicator component type displays the defined temperature at its Measure
input.
Operating window
The temperature value obtained via the Measure input is displayed in the operating window
(Figure 4-41).
Symbol
Function
The WeightIndicator component type displays the defined weight at its Measure input.
Operating window
The weight value obtained via the Measure input is displayed in the operating window
(Figure 4-42).
Symbol
Function
The NetWS component type is used to parameterise a network for the water/steam medium.
The components can be added to any point on any branch of the flownet.
Parameter
The variables for the flownet can be specified as parameters for the components:
• FactorMomentum
Momentum factor for the flow through a branch of the flownet
• sCompressionSteam
Specific compression module for density ρ < 500 kg/m³ (steam)
• sCompressionWater
Specific compression module for density ρ > 500 kg/m³ (water)
• FactorThermalSteam
Factor for enthalpy balancing for density ρ < 500 kg/m³ (steam)
• FactorThermalWater
Factor for enthalpy balancing for density ρ > 500 kg/m³ (water)
The parameters, their units and default values are shown in Figure 4-44.
Additional parameters
Initial values and specific characteristics for internal nodes of the flownet can be set via
additional parameters:
• PressureInit
Initialisation value for the pressure in the internal node of the flownet
• sEnthalpyInit
Initialisation value for the specific enthalpy in the internal node of the flownet
• SmoothTransition
When this additional parameter is set to true, the variables sCompression and
FactorThermal are set with a density dependant linear transfer function
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the medium in the internal nodes
and the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-45.
Symbol
Function
The PnodeWS component type defines values for the pressure p and specific enthalpy h at
its connection A. This type of component forms a boundary for the flownet. When the flownet
is represented as a graph, this corresponds to an (external) node, for which the pressure
and specific enthalpy are predefined.
Operating window
The values for pressure p and specific enthalpy h can be input digitally in the operating
window (Figure 4-46).
Figure 4-47: Extended operating window for the component type PnodeWS
Function
The MnodeWS component type defines values for the mass flow rate m and specific
enthalpy h at its connection A. This type of component forms a boundary for the flownet.
When the flownet is represented as a graph, this corresponds to an in or outflow through an
(internal) node or branch. An internal node with defined inflow or outflow is added to the
flownet for each component of this type.
Operating window
The values for mass flow rate m and specific enthalpy h can be input digitally in the operating
window (Figure 4-48).
Figure 4-49: Extended operating window for the MnodeWS component type
Symbol
Function
The component type JointWS can be used to join three branches connected at A, B and C in
one node. The JointWS component adds an internal node to the flownet.
Operating window
The pressure p, specific enthalpy h, density r and the temperature T of the node are shown
in the operating window (Figure 4-50).
Symbol
Function
The component type JointParamWS can be used to join three branches connected at A, B
and C in one node. The JointParamWS component adds an internal node to the flownet.
Parameter
The node associated with JointParamWS can be parameterised:
• sCompressionSteam
Specific compression module for density ρ < 500 kg/m³ (steam)
• sCompressionWater
Specific compression module for density ρ > 500 kg/m³ (water)
• FactorThermalSteam
Factor for enthalpy balancing for density ρ < 500 kg/m³ (steam)
• FactorThermalWater
Factor for enthalpy balancing for density ρ > 500 kg/m³ (water)
The parameters, their units and default values are shown in Figure 4-51.
Additional parameters
Initial values and specific characteristics for the node can be set via additional parameters:
• PressureInit
Initial value for the pressure in the node
• sEnthalpyInit
Initial value for the specific enthalpy in the node
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the medium in the internal nodes
and the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-52.
Operating window
The pressure p, specific enthalpy h, density r and the temperature T of the node are shown
in the operating window (Figure 4-56).
Symbol
Function
The StorageTankWS component type provides the simulation with an open water tank, i.e. a
tank that is not sealed from the environment.
The medium flows in and out of the tank via STUBx. For each of the N connections, the
throttling effect is defined by
2
∆p − m 2 sec
= 12960 .
bar 2 kg h
cV ρ 3
m
The tank can have a minimum of one and a maximum of 16 connections. The connections
can be moved to any position on the outline of the component symbol by using the mouse
while simultaneously pressing ALT.
It is assumed that the water in the tank is immediately completely mixed i.e. it is an
homogenous medium with an overall consistent density and enthalpy.
The inflow and outflow of water is balanced across the N connections of the tank. The
balanced mass M of the N flow rates m i of the medium in the tank, is defined by
dM N
∑
= m i .
dt i =1
The water volume balancing via the volume flow results in the rate of change of density ρ
given by
dρ ρ N N
1
=
dt M
∑ m i − ρ ∑
i =1,i∈Z ρ i
m i ,
i =1,i∈Z
where only the inflows ( i ∈ Z ) need to be summed.
The specific enthalpy of the water is defined by the enthalpy balance equation
dh 1 N N
= ∑hi m i − h
dt M i =1,i∈Z
∑
m i
i =1,i∈Z
Again, only the inflows ( i ∈ Z ) need to be summed.
The dynamic behaviour of the water in the tank is described by these balances for mass M,
density ρ and specific enthalpy h.
The calculated variables for the level l of water in the tank, its mass M and temperature T as
well as pressure p 0 on the tank base are output at the Measure connection. The temperature
T is calculated from the density and specific enthalpy of the water using the equation of state
T = T (ρ , h ) .
The pressure p on the tank base is calculated as the sum of the weight pressure ρ g l of the
water and the atmospheric pressure p U according to
p = pU + ρ g l .
CAUTION
It is not checked whether the state variables (pressure, enthalpy) for the inflow
and tank contents always possess values for the state water.
In the simulation, corresponding specified inflows may result in a situation where
the state values for the tank contents describe a water-steam mixture or pure
steam, which has no physical meaning.
Parameter
All relevant geometric variables of the tank are set via parameters:
• Volume
Volume V of the tank; adjustable online
• Height
Height of the tank; adjustable online
• NbrOfStubs
Number N of connections
The parameters, their units and default values are shown in Figure 4-56 .
Additional parameters
The atmospheric pressure PressureOutside and initial value for the tank fill level (LevelInit),
the specific enthalpy (EnthalpyInit) and density (DensityInit) of the water are set via
additional parameters:
• PressureOutside
Atmospheric pressure p U ; adjustable online
• LevelInit
Initial value for the level
• EnthalpyInit
Initial value for the specific enthalpy in the contents
• DensityInit
Initial value for the density of the contents
• Cvs
A uniform flow coefficient c V for all connections
• Cv0
A uniform flow coefficient c V0 for strong throttling at the tank connections
• MinVolume
Minimal tank volume V min ; adjustable online
• MinVolumeHys
Hysteresis H min ; adjustable online
The additional parameters, their units and default values are shown in Figure 4-57.
Operating window
The level l, the pressure p, the temperature T and the mass M of the water in the tank is
displayed in the operating window (Figure 4-58).
Symbol
Function
The DrumWS component type provides the simulation with a drum, i.e. a closed container
that is water/steam tight.
The drum is assumed to be a cylindrical, horizontal or vertically standing container, where
media flow in and out via the STUBWx and/or STUBSx connections. For each of the
connections, the throttling effect is defined by
2
∆p − m 2 sec
= 12960
bar 2 kg h
cV ρ 3
m
At the N W connections STUBWx medium with the state variables of the saturated water can
be drawn out, likewise medium with the state variables of the saturated steam can be drawn
out at the N S connections STUBSx. The drum can have a minimum of one and a maximum
of eight connections of each type. The connections can be moved to any position on the
outline of the component symbol by using the mouse while simultaneously pressing ALT.
It can be assumed that the medium in the drum immediately separates into two homogenous
saturated phases: the saturated liquid phase and the saturated steam phase.
The inflow and outflow of water/steam are balanced via the N = N W + N S drum connections.
This results in a mass of water and steam in the drum, balanced via the flow rate, given by
dM N
∑
= m i
dt i =1
Given that
M = ρV
the rate of change in the mean density results in
N
dρ
V
dt
= ∑ m
i =1
i .
where h is the specific enthalpy (mean enthalpy) of water/steam in the drum and h i is the
specific enthalpy of the inflow and outflow at the i-th connection. The specific enthalpy of the
medium in the drum must be applied for outflows, i.e. h ′ for outflow of saturated water and
h ′′ for outflow of saturated steam. The balancing results from
N N N
dh
M = ∑ ∑ ∑
(hi − h )m i + (h ′ − h ) m i + (h ′′ − h ) m i
dt i =1,i∈Z i =1,i∈ A′ i =1,i∈ A′′
where the inflows Z , outflows of saturated water A′ and/or outflows of saturated steam A′′
are summed.
In order to model a heat exchange with an ideal insulated drum, the enthalpy balance
equation is extended as follows
N N N
dh
M = ∑ ∑ ∑
(hi − h )m i + (h ′ − h ) m i + (h ′′ − h ) m i + Aα(TT − TS ) .
dt i =1,i∈Z i =1,i∈ A′ i =1,i∈ A′′
Parameter
All relevant geometric variables of the drum are set via parameters:
• NbrOfStubsW
Number N W of connections STUBWx
• NbrOfStubsS
Number N S of connections STUBSx
• Position
Position of the drum: Vertically or Horizontally
• VolumeDrum
Volume V of the drum; adjustable online
• HeightOrLength
Height or length of the drum; adjustable online
• MassDrum
Mass M T of the drum; adjustable online
• SurfaceDrum
Inner surface A of the drum; adjustable online
• sHeatCapDrum
Specific heat capacity c T of the drum; adjustable online
• HeatTransCoe
Heat transfer coefficient α of the drum for water/steam ; adjustable online
The parameters, their units and default values are shown in Figure 4-61.
Additional parameters
Initial values and other variables can be set via additional parameters:
• VolumeInit
Initial value for the water volume
• TemperatureInit
Initial value for the temperature of water/steam (saturation temperature)
• Cvs
Uniform flow coefficient c V for all connections
• Cv0
Uniform flow coefficient c V0 for strong throttling at the drum connections
• MinVolume
Minimal drum volume V min ; adjustable online
• MinVolumeHys
Hysteresis H min ; adjustable online
The additional parameters, their units and default values are shown in Figure 4-62.
Operating window
The saturation pressure p S and the pressure p W at the STUBWx connections are displayed in
the operating window (Figure 4-63). The saturation pressure T S , the drum temperature T T ,
the mass M of water/steam and the water level l are also displayed.
The values for density r and specific enthalpy h of water/steam are shown in the extended
operating window (Figure 4-64). Mix indicates the mean variables, Water the variables for
saturated water and Steam the variables for saturated steam.
Figure 4-64: Extended operating window for the DrumWS component type
Symbol
Function
The ElectricalHeaterWS component type is used to simulate an electric heat exchanger. The
electric heating power P el in kW is set via the connection ElectricalPower.
Water/steam is directed via the connections TS_IN and TS_OUT as a heated medium. The
flow m is throttled according to the relation
2
∆p − m 2 sec
= 12900 .
bar kg h
cV2 ρ TS
m3
with the flow coefficient c V . The reference direction for the flow is chosen as from TS_IN to
TS_OUT.
For the heat exchange, a simple one tube model divided into segments is applied (Figure
4-65). The number N of segments can be set to a value between 4 and 16.
Figure 4-65: Diagram of the segments for the electrical heat exchanger
ElectricalHeaterWS
It is assumed that the supplied electrical energy is completely converted into heat. The heat
balance for a segment i of the tube with extremely simplified heat transfer is given by
dTT ,i
= aT Pel + bT (TTS ,i − TT ,i )
dt
where
1 A α
aT = , bT = TS TS .
M T cT M T cT
T T,i and T TS,i are the temperatures of the tube segments and the medium in the segments,
h T,i and h TS,i are the corresponding specific enthalpies.
Parameter
The segmentation and coefficients for the balance equations are defined by parameters:
• NbrOfSegments
Number N of segments: 4 ≤ N ≤ 16
• Cvs
Flow coefficient c V
• Volume
Volume V TS of the medium (tube inside volume)
• Surface
Surface A TS of the tube on the water/steam side
• HeatTransCoef
Heat transfer coefficient α TS from tube to water/steam
• sHeatCapTube
Specific heat capacity c T of the tube
• MassTube
Mass M T of the tube
The parameters, their units and default values are shown in Figure 4-66.
Operating window
The flow rate m , pressure drop ∆p and enthalpy difference ∆h of the heated medium as well
as the electrical heating power P el are displayed in the operating window (Figure 4-67). The
temperatures of the tube (Tube), and the heated medium (Tube side) for the first and last
segment (T 1 and T N ) are also displayed.
Symbol
Function
The HeatExchangerWS component type is used to simulate a heat exchanger for the media
water/steam on the tube side and shell side. The simulation is implemented for three types
• parallel-flow heat exchanger
• counter-flow heat exchanger and
• cross-flow heat exchanger.
Both media are routed on the tube side via the connections TS_IN and TS_OUT, and on the
shell side via the connections SS_IN and SS_OUT.
The flow m is throttled on both the tube and shell side according to the relation
2
∆p − m 2 sec
= 12900
bar 2 kg h
cV ρ 3
m
with the relevant flow coefficient c v . The chosen reference direction for the flow is from
connection _IN to _OUT on both the tube and shell side.
For the heat exchange, a simple one tube model divided into segments is applied (Figure
4-68). The number N of segments can be set to a value between 4 and 16.
The water/steam heat on the tube and shell side and the heat in the tube itself (Figure 4-69)
are balanced.
Both heat transfers from the shell side medium to the tube and from the tube to the tube side
medium are applied, in a very simplified form, as
Q SS −T = ASS q SS −T = ASS α SS (TSS − TT ) ,
dTT ,i
= aT (TSS ,i − TT ,i ) + bT (TTS ,i − TT ,i ) ,
dt
dhSS ,i
= a SS (hSS ,i −1 − hSS ,i ) + bSS (TT ,i − TSS ,i )
dt
apply with the following coefficients, which are valid for all segments:
N m TS A α N m SS A α A α A α
aTS = , bTS = TS TS , a SS = , bSS = SS SS , aT = SS SS , bT = TS TS .
ρ TS VTS ρ TS VTS ρ SS V SS ρ SS V SS M T cT M T cT
The values set for the media in the shell and tube side of the flownet apply to the density ρ TS ,
ρ SS and the specific heat capacity c TS , c SS .
For initialisation, the temperatures of the tube segments are set to the temperature of the
tube side medium, calculated from the pressure (input [Link]) and the specific
enthalpy (output [Link]).
Parameter
The segmentation and coefficients for the balance equations are defined by parameters:
• Type
Type of heat exchanger: ParallelFlow (parallel flow), CounterFlow (counter flow),
CrossFlow (cross flow)
• NbrOfSegments
Number N of segments: 4 ≤ N ≤ 16
• CvsSS
Shell side flow coefficient c V
• CvsTS
Tube side flow coefficient c V
• VolumeSS
Shell side volume V SS
• VolumeTS
Tube side volume V TS
• SurfaceSS
Shell side surface A SS of the tube (exterior surface of the tube)
• SurfaceTS
Tube side surface A TS of the tube (interior surface of the tube)
• HeatTransCoefSS
Heat transfer coefficient α SS on the tube exterior
• HeatTransCoefTS
Heat transfer coefficient α TS on the tube interior
• sHeatCapTube
Specific heat capacity c T of the tube
• MassTube
Mass M T of the tube
The parameters, their units and default values are shown in Figure 4-70.
Operating window
The flow rate m , pressure drop ∆p and enthalpy difference ∆h of the shell side and tube side
media are displayed in the operating window. The temperatures of the tube (Tube), and the
medium (Shell side and Tube side) for the first and last segment (T 1 and T N ) are also
displayed.
Symbol
Function
The NetLiquid component type is used for parameterisation of a flownet for liquids. The
components can be added to any point on any branch of the flownet.
Parameter
The variables for the flownet can be specified as parameters for the components:
• Density
Density of the liquid
• sHeatCapacity
Specific heat capacity of the liquid
• FactorMomentum
Momentum factor for the flow through a branch of the flownet
• sCompression
Specific compression modulus of the liquid
• FactorThermal
Factor for enthalpy balancing
The parameters, their units and default values are shown in Figure 4-73.
Additional parameters
Initial values and specific characteristics for internal nodes of the flownet can be set via
additional parameters:
• PressureInit
Initialisation value for the pressure in the internal nodes of the flownet
• TemperatureInit
Initialisation value for the temperature in the internal nodes of the flownet
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the medium in the internal nodes
and the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-74.
Symbol
Function
The PnodeLiquid component type outputs values for pressure and temperature at its
connection. It creates a boundary for the connected flownet. When the flownet is
represented as a graph, PnodeLiquid corresponds to an (external) node, for which the
pressure and temperature are predefined.
Operating window
The values for pressure and temperature can be input digitally in the operating window
(Figure 4-75).
Figure 4-76: Extended operating window for the PnodeLiquid component type
Symbol
Function
The MnodeLiquid component type outputs values for mass flow rate and temperature at its
connection. It creates a boundary for the connected flownet. When the flownet is
represented as a graph, MnodeLiquid corresponds to an in or outflow through an (internal)
node or branch. Thus an internal node with a defined in or outflow is added to the flownet.
Operating window
The values for mass flow rate and temperature can be input digitally in the operating window
(Figure 4-77).
Figure 4-78: Extended operating window for the MnodeLiquid component type
Symbol
Function
The JointLiquid component type can be used to join three connected branches in one node.
All connections are equal. The JointLiquid component adds an internal node to the flownet.
Operating window
The pressure p and temperature T of the node is displayed in the operating window (Figure
4-79).
Symbol
Function
The JointParamLiquid component type can be used to join three connected branches in one
node. All connections are equal. The JointParamLiquid component adds an internal node to
the flownet. This node can be parameterised.
Parameter
The node associated with JointParamLiquid can be parameterised:
• sHeatCapacity
Specific heat capacity of the liquid
• sCompression
Specific compression module of the liquid
• FactorThermal
Factor for enthalpy balancing
The parameters, their units and default values are shown in Figure 4-80.
Additional parameters
Initial values and specific characteristics for the node can be set via additional parameters:
• PressureInit
Initial value for the pressure in the node
• TemperatureInit
Initial value for the temperature in the node
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the liquid in the internal nodes and
the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-81.
Operating window
The pressure p and temperature T of the node is displayed in the operating window (Figure
4-82).
Function
The StorageTankLiquid component type provides the simulation with an open tank for liquid,
i.e. a tank that is not sealed from the environment.
The medium flows in and out of the tank via STUBx connections. For each of the N
connections, a throttling effect defined by
2
∆p − m 2 sec
= 12960
bar kg h
cV2 ρ 3
m
is assumed. The tank can have a minimum of one and a maximum of 16 connections. The
connections can be moved to any position on the outline of the component symbol by using
the mouse while simultaneously pressing ALT.
It is assumed that the liquid in the tank is immediately completely mixed i.e. it is an
homogenous medium with an overall consistent density and temperature.
The inflow and outflow is balanced across the N connections of the tank. The balanced mass
M of the N flow rates m i of the liquid in the tank, is defined by
dM N
∑
= m i .
dt i =1
The temperature of the liquid is balanced as the mixing temperature from the inflows ( i ∈ Z )
according to
dT 1 1 N N
=
dt M c L
∑ hi m i − T ∑ m . i
i =1,i∈Z i =1,i∈Z
The dynamic behaviour of the liquid in the tank is described by these balances for mass M,
and temperature T.
The calculated variables for the level l of liquid in the tank, its mass M and temperature T as
well as pressure p 0 on the tank base are output at the Measure connection. The pressure p
on the tank base is the sum of the weight pressure ρ g l of the liquid and the atmospheric
pressure p U :
p = pU + ρ g l .
Parameter
The following are set as geometric variables of the tank via parameters:
• Volume
Volume V of the tank; adjustable online
• Height
of the tank; adjustable online
• NbrOfStubs
Number N of connections
The parameters, their units and default values are shown in Figure 4-85.
Additional parameters
The initialisation values for the atmospheric pressure PressureOutside, the fill level
(LevelInit) of the tank and the temperature (TemperatureInit) of the liquid are set via
additional parameters:
• PressureOutside
Atmospheric pressure p U ; adjustable online
• LevelInit
Initial value for the level
• TemperatureInit
Initial value for the temperature of the contents
• Cvs
Uniform flow coefficient k V for all connections
• Cv0
Uniform flow coefficient k V0 for maximum throttling at the tank connections
• Density
Density ρ of the liquid in the tank
• sHeatCapacity
Specific heat capacity c G of the liquid in the tank
• MinVolume
Minimal tank volume V min ; adjustable online
• MinVolumeHys
Hysteresis H min ; adjustable online
The additional parameters, their units and default values are shown in Figure 4-86.
Operating window
The level l, the pressure p, the temperature T and the mass M of the liquid in the tank is
displayed in the operating window (Figure 4-87).
Function
The ElectricalHeaterLiquid component type is used to simulate an electric heat exchanger.
The electric heating power P el in kW is set via the connection ElectricalPower.
The liquid is directed via the connections TS_IN and TS_OUT as a heated medium. The flow
m is throttled according to the relation
2
∆p − m 2 sec
= 12900 .
bar 2 kg h
cV ρ 3
m
with flow coefficient c V . The reference direction for the flow is chosen as from TS_IN to
TS_OUT.
For the heat exchange, a simple one tube model divided into segments is applied (Figure
4-88). The number N of segments can be set to a value between 4 and 16.
It is assumed that the supplied electrical energy is completely converted into heat. The heat
balance for a segment i of the tube with extremely simplified heat transfer is given by
dTT ,i
= aT Pel + bT (TL ,i − TT ,i )
dt
where
1 Aα
aT = , bT = .
M T cT M T cT
T T,i and T L,i are the temperatures of the tube segments and the liquids in the segment i. The
flownet values apply for the density ρ and the specific heat capacity c G of the liquid.
Parameter
The segmentation and coefficients for the balance equations are defined by parameters:
• NbrOfSegments
Number N of segments: 4 ≤ N ≤ 16
• Cvs
Flow coefficient c V
• Volume
Volume V of the liquid (tube inside volume)
• Surface
Surface A of the tube interior
• HeatTransCoef
Heat transfer coefficient α from tube to liquid
• sHeatCapTube
Specific heat capacity c T of the tube
• MassTube
Mass M T of the tube
The parameters, their units and default values are shown in Figure 4-89.
Additional parameters
For initialisation of the heat exchanger, the temperatures of the liquid and the tube can be
set to the same value (Figure 4-90) via the additional parameter TemperatureInit.
Operating window
The flow rate m , pressure drop ∆p and temperature difference ∆T of the liquid as well as the
electrical heating power P el are displayed in the operating window (Figure 4-91). The
temperatures of the tube (Tube), and the liquid (Tube side) for the first and last segment (T 1
and T N ) are also displayed.
Symbol
Function
The HeatExchangerLiquid component type is used to simulate a heat exchanger for liquids
on the tube side and shell side. The simulation is implemented for three types
• parallel-flow heat exchanger
• counter-flow heat exchanger and
• cross-flow heat exchanger.
Both media are routed on the tube side via the connections TS_IN and TS_OUT, and on the
shell side via the connections SS_IN and SS_OUT.
The flow m is throttled on both the tube and shell side according to the relation
2
∆p − m 2 sec
= 12900 .
bar kg h
cV2 ρ 3
m
with the relevant flow coefficient c v . The chosen reference direction for the flow is from
connection _IN to _OUT on both the tube and shell side.
For the heat exchange, a simple one tube model divided into segments is applied (Figure
4-92). The number N of segments can be set to a value between 4 and 16.
The liquid heat on the tube and shell side and the heat in the tube itself (Figure 4-93) are
balanced.
Figure 4-93: Diagram of the heat balancing for the heat exchanger
Both heat transfers from the shell side liquid to the tube and from the tube to the tube side
liquid are applied, in a very simplified form, as
Q SS −T = ASS q SS −T = ASS α SS (TSS − TT ) ,
The values set for the liquid in the shell and tube side of the flownet apply to the density ρ TS ,
ρ SS and the specific heat capacity c TS , c SS .
For initialisation, the temperatures of the tube segments are set to the temperature of the
tube side liquid, calculated from the specific enthalpy ( [Link] input).
Parameter
The segmentation and coefficients for the balance equations are defined by parameters:
• Type
Type of heat exchanger: ParallelFlow (parallel flow), CounterFlow (counter flow),
CrossFlow (cross flow)
• NbrOfSegments
Number N of segments: 4 ≤ N ≤ 16
• CvsSS
Shell side flow coefficient c V
• CvsTS
Tube side flow coefficient c V
• VolumeSS
Shell side volume V SS
• VolumeTS
Tube side volume V TS
• SurfaceSS
Shell side surface A SS of the tube (exterior surface of the tube)
• SurfaceTS
Tube side surface A TS of the tube (interior surface of the tube)
• HeatTransCoefSS
Heat transfer coefficient α SS on the tube exterior
• HeatTransCoefTS
Heat transfer coefficient α TS on the tube interior
• sHeatCapTube
Specific heat capacity of the tube c T
• MassTube
Mass M T of the tube
The parameters, their units and default values are shown in Figure 4-94:
Operating window
The flow rate m , pressure drop ∆p and enthalpy difference ∆h of the shell side and tube side
liquid are displayed in the operating window (Figure 4-95). The temperatures of the tube
(Tube), and the liquid (Shell side and Tube side) for the first and last segment (T 1 and T N )
are also displayed.
Function
The NetGas component type is used for parameterisation of a flownet for gases. The
components can be added to any point on any branch of the flownet.
Parameter
The variables for the flownet can be specified as parameters for the components:
• GasConstant
Specific gas constant
• sHeatCapacity
Specific heat capacity of the gas
• FactorMomentum
Momentum factor for the flow through a branch of the flownet
• sCompression
Specific compression modulus of the gas
• FactorTherma
Factor for enthalpy balancing
The parameters, their units and default values are shown in Figure 4-97.
Additional parameters
Initial values and specific characteristics for internal nodes of the flownet can be set via
additional parameters:
• PressureInit
Initialisation value for the pressure in the internal nodes of the flownet
• TemperatureInit
Initialisation value for the temperature in the internal nodes of the flownet
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the gases in the internal nodes and
the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-98.
Symbol
Function
The PnodeGas component type outputs values for pressure and temperature at its
connection. It creates a boundary for the connected flownet. When the flownet is
represented as a graph, PnodeGas corresponds to an (external) node, for which the
pressure and temperature are predefined.
Operating window
The values for pressure and temperature can be input digitally in the operating window
(Figure 4-99).
Figure 4-100: Extended operating window for the component type PnodeGas
Symbol
Function
The MnodeGas component type outputs values for mass flow rate and temperature at its
connection. It creates a boundary for the connected flownet. When the flownet is
represented as a graph, Mnode corresponds to an in or outflow through an (internal) node or
branch. Thus an internal node with a defined in or outflow is added to the flownet.
Operating window
The values for mass flow rate and temperature can be input digitally in the operating window
(Figure 4-101).
Figure 4-102: Extended operating window for the component type MnodeGas
Symbol
Function
The JointGas component type can be used to join three connected branches in one node. All
connections are equal. The JointGas component adds an internal node to the flownet.
Operating window
The pressure p, temperature T and density r of the node are displayed in the operating
window (Figure 4-103).
Symbol
Function
The JointParamGas component type can be used to join three connected branches in one
node. All connections are equal. The JointParamGas component adds an internal node to
the flownet. This node can be parameterised.
Parameter
The node associated with JointParamGas can be parameterised:
• GasConstant
Specific gas constant
• sHeatCapacity
Specific heat capacity of the gas
• sCompression
Specific compression modulus of the gas
• FactorThermal
Factor for enthalpy balancing
The parameters, their units and default values are shown in Figure 4-104.
Additional parameters
Initial values and specific characteristics for the node can be set via additional parameters:
• PressureInit
Initialisation value for the pressure in the node
• TemperatureInit
Initialisation value for the temperature in the node
• TemperatureEnvironment
Ambient temperature
• FactorHeatExchangeEnv
Proportionality factor for heat exchange between the gas in the internal nodes and
the environment; at zero value no heat exchange occurs
The additional parameters, their units and default values are shown in Figure 4-105.
Operating window
The pressure p, temperature T and density r of the node are displayed in the operating
window (Figure 4-106).
Symbol
Function
The StorageTankGas component type provides the simulation with a gas tank.
The medium flows in and out of the tank via STUBx connections. For each of the N
connections, a throttling effect defined by
2
∆p − m 2 sec
= 12960
bar kg h
cV2 ρ 3
m
is assumed. The tank can have a minimum of one and a maximum of 16 connections. The
connections can be moved to any position on the outline of the component symbol by using
the mouse while simultaneously pressing ALT.
It is assumed that the gas in the tank is immediately completely mixed i.e. it is an
homogenous medium with an overall consistent density and temperature.
The inflow and outflow is balanced across the N connections of the tank. The balanced mass
M of the N flow rates m i of the gas in the tank, is defined by
dM N
∑
= m i .
dt i =1
The specific enthalpy h of the gas is balanced as mixed enthalpy from the inflows ( i ∈ Z ) as
per
dh 1 N N
= ∑hi m i − h
dt M i =1,i∈Z
∑ m i
i =1,i∈Z
The dynamic behaviour of the gas in the tank is described by these balances for mass M
and specific enthalpy h.
At the connection Measure, the variables for pressure p, temperature
h
T=
cG
and mass M of the gas are output. The pressure is calculated from the equation for ideal
gas:
R S (T + 273,15 K ) M
p= .
5 Pa
V ⋅ 10
bar
Parameter
The following are set as geometric variables of the tank via parameters:
• Volume
Volume V of the tank; adjustable online
• NbrOfStubs
Number N of connections
The parameters, their units and default values are shown in Figure 4-109 .
Additional parameters
Initial values and other variables can be predetermined via additional parameters:
• PressureInit
Initialisation value for the pressure; adjustable online
• TemperatureInit
Initialisation value for the temperature of the gases
• PressureMax
Maximum tank pressure p max ; adjustable online
• Cvs
Uniform flow coefficient c V for all connections
• Cv0
Uniform flow coefficient c V0 for maximum throttling at the tank connections
• GasConstant
Specific gas constant R S for the gas in the tank
• sHeatCapacity
Specific heat capacity c G of the gas in the tank
• MinMass
Minimum tank filling M min of the tank; adjustable online
• MinVolumeHys
Hysteresis H min ; adjustable online
The additional parameters, their units and default values are shown in Figure 4-110.
Operating window
The pressure p, temperature T and mass M of the gases in the tank are displayed in the
operating window (Figure 4-111).
Symbol
Function
The ElectricalHeaterGas component type is used to simulate an electric heat exchanger for
gas. The electric heating power P el in kW is set via the connection ElectricalPower.
The gas is directed via the connections TS_IN and TS_OUT as a heated medium. The flow
m is throttled according to the relation
2
∆p − m 2 sec
= 12900
bar 2 kg h
cV ρ 3
m
with flow coefficient c V . The reference direction for the flow is chosen as from TS_IN to
TS_OUT.
For the heat exchange, a simple one tube model divided into segments is applied (Figure
4-112). The number N of segments can be set to a value between 4 and 16.
It is assumed that the supplied electrical energy is completely converted into heat. The heat
balance for a segment i of the tube with extremely simplified heat transfer is given by
dTT ,i
= aT Pel + bT (TG ,i − TGT ,i )
dt
where
1 Aα
aT = , bT = .
M T cT M T cT
T T,i and T G,i are the temperatures of the tube segments and the gases in the segment i. The
flownet values apply for the density ρ and the specific heat capacity c G of the gases.
Parameter
The segmentation and coefficients for the balance equations are defined by parameters:
• NbrOfSegments
Number N of segments: 4 ≤ N ≤ 16
• Cvs
Flow coefficient c V
• Volume
Volume V of the gases (tube inside volume)
• Surface
Surface A of the tube interior
• HeatTransCoef
Heat transfer coefficient α from tube to gas
• sHeatCapTube
Specific heat capacity c T of the tube
• MassTube
Mass M T of the tube
The parameters, their units and default values are shown in Figure 4-113.
Additional parameters
For initialisation of the heat exchanger, the temperatures of the gas and the tube can be set
to the same value (Figure 4-114) via the additional parameter TemperatureInit.
Operating window
The flow rate m , pressure drop ∆p and temperature difference ∆T of the gas as well as the
electrical heating power P el are displayed in the operating window (Figure 4-115). The
temperatures of the tube (tube), and the gas (tube side) for the first and last segment (T 1
and T N ) are also displayed.
Symbol
Function
The HeatExchangerGas component type is used to simulate a heat exchanger for the gases
on the tube side and shell side. The simulation is implemented for three types
• parallel-flow heat exchanger
• counter-flow heat exchanger and
• cross-flow heat exchanger.
Both media are routed on the tube side via the connections TS_IN and TS_OUT, and on the
shell side via the connections SS_IN and SS_OUT.
The flow m is throttled on both the tube and shell side according to the relation
2
∆p − m 2 sec
= 12900
bar 2 kg h
cV ρ 3
m
with the relevant flow coefficient c v . The chosen reference direction for the flow is from
connection _IN to _OUT on both the tube and shell side.
For the heat exchange, a simple one tube model divided into segments is applied (Figure
4-116). The number N of segments can be set to a value between 4 and 16.
The gas heat on the tube and shell side and the heat in the tube itself (Figure 4-117) are
balanced.
Figure 4-117: Diagram of the heat balancing for the heat exchanger
Both heat transfers from the shell side gas to the tube and from the tube to the tube side gas
are applied, in a very simplified form, as
Q SS −T = ASS q SS −T = ASS α SS (TSS − TT ) ,
dTTS ,i
= aTS (TTS ,i −1 − TTS ,i ) + bTS (TT ,i − TTS ,i ) ,
dt
dTT ,i
= aT (TSS ,i − TT ,i ) + bT (TTS ,i − TT ,i ) ,
dt
dTSS ,i
= a SS (TSS ,i −1 − TSS ,i ) + bSS (TT ,i − TSS ,i )
dt
apply with the following coefficients, which are valid for all segments:
N m TS A α N m SS ASS α SS A α A α
aTS = , bTS = TS TS , a SS = , bSS = , aT = SS SS , bT = TS TS .
ρ TS VTS cTS ρ TS VTS ρ SS V SS c SS ρ SS V SS M T cT M T cT
The values set for the gas in the shell and tube side of the flownet apply to the density ρ TS ,
ρ SS and the specific heat capacity c TS , c SS .
For initialisation, the temperatures of the tube segments are set to the temperature of the
tube side gas, calculated from the specific enthalpy ( [Link] input).
Parameter
The segmentation and coefficients for the balance equations are defined by parameters:
• Type
Type of heat exchanger: ParallelFlow (parallel flow), CounterFlow (counter flow),
CrossFlow (cross flow)
• NbrOfSegments
Number N of segments: 4 ≤ N ≤ 16
• CvsSS
Shell side flow coefficient c V
• CvsTS
Tube side flow coefficient c V
• VolumeSS
Shell side volume V SS
• VolumeTS
Tube side volume V TS
• SurfaceSS
Shell side surface A SS of the tube (exterior surface of the tube)
• SurfaceTS
Tube side surface A TS of the tube (interior surface of the tube)
• HeatTransCoefSS
Heat transfer coefficient α SS on the tube exterior
• HeatTransCoefTS
Heat transfer coefficient α TS on the tube interior
• sHeatCapTube
Specific heat capacity of the tube c T
• MassTube
Mass M T of the tube
The parameters, their units and default values are shown in Figure 4-118.
Operating window
The flow rate m , pressure drop ∆p and enthalpy difference ∆h of the shell side and tube side
gas are displayed in the operating window (Figure 4-119). The temperatures of the tube
(Tube), and the gas (Shell side and Tube side) for the first and last segment (T 1 and T N ) are
also displayed.
If a connecting line is used to join topological connectors to each other in the diagram editor,
then the connectors are hidden and only visible as faint shadows (Figure 5-3).
They can also be joined by superimposing the connectors. Once linked, both connectors are
hidden (Figure 5-4).
Pressure boosts are applied in the reference direction, i.e. the pressure variation is defined
here as
∆p = p b − p a
For pressure drops in the reference direction (throttling), ∆p is therefore negative and for
pressure boosts in the reference direction ∆p is positive.
Figure 5-8: Data exchange between flownet components and flownet solver
Connection types FNL2 to FNL7 only provide connection to the flownet solver and thus must
be set as invisible on the component symbol. In the connection properties, the usage
"Property view only" (see Figure 5-9) or "CTE view only" must be set.
Figure 5-11 shows the input and output signals with the direction of data flow between the
components and the flownet solver.
The components receive four variables via the input signals for calculating the effect of the
branch object:
1. specific enthalpy h ([Link]),
2. mass flow m ([Link]),
3. pressure p a ([Link]),
4. densityρ ([Link]).
The pressure p a relates to the "from" connector (FROM) of the topological description, which
is the connector a in our example. Density ρ and enthalpy h are variables of the supplied
medium (Figure 5-12).
The output signals for the branch object are calculated in the component and sent to the
flownet solver:
1. change in specific enthalpy ∆h ([Link]),
2. pressure variation ∆p = p b – p a ([Link]),
3. derivative of pressure variation with respect to the mass flow d∆p/d M ([Link]).
The following applies: if the FN connector is defined with the direction OUT, then the
pressure and specific enthalpy are specified via this connector for the external node of the
flownet. If the connector is defined with the direction IN then the mass flow and the specific
enthalpy are specified for the external node.
The mass flow m is positive flow into the flownet ( m > 0 ) and for outflow it is negative
( m > 0 ). This does not depend on the direction choosen in the topology description.
Mass flow nodes operate in the flownet like internal nodes with mass inflow or outflow.
These nodes are therefore treated as internal nodes in the flownet, i.e. for an external node,
an internal node with additional inflow or outflow is created in the flownet. The pressure
returned to the components p e thus corresponds to the pressure in this internal node and
likewise for the specific enthalpy h e (Figure 5-16):
The variables calculated in the node i are provided via the connector inputs (see Figure
5-18):
1. pressure p i ([Link])
2. specific enthalpy h i ([Link])
3. density ρ i ([Link])
NOTE
The relations used for calculating the flownet variables, such as temperature, for
the relevant medium are defined via the MEDIUM parameter. This can only be
done during initialisation of the simulation. Changes made to the parameters
while the simulation is running have no effect.
Signal Description
MEDIUM Flownet medium:
MEDIUM := 0 for water/steam,
MEDIUM := 1 for ideal gas,
MEDIUM := 2 for liquid;
cannot be changed while the simulation is running
CG Specific compression modulus K / M for water/steam medium
with density ρ < 500 m³/kg and for ideal gas medium
Depending on the direction of the connector FN its signals are inputs or outputs. Depending
on the direction, flownet parameters can be set or used in components. If the connector is
defined by direction OUT, then the variables set in the components are provided to the
flownet as parameters. These parameters are available in the components for evaluation, if
the connector is defined with direction IN.
If the connector is defined with direction IN, then the momentum factor Α for the branch is
transferred to the branch object from the flownet solver.
The variables listed in Table 5-3 are transmitted to the flownet solver for the internal node.
Their meaning in the mathematical model of the flownet is explained in section 3.6.
Signal Description
CG Specific compression module K i / M i for water/steam medium with
density ρ < 500 m³/kg and for ideal gas medium
5.3.1 Constants
The available constants are listed in Table 5-4.
5.3.2 Functions
Functions for calculating state variables are available for flownet components with
water/steam medium. The available state variables and units are summarised in Table 5-5.
Variable Unit
p Pressure bar
pS Saturated pressure bar
T Temperature °C
TS Saturated temperature °C
h Specific enthalpy kJ/kg
h’ Specific enthalpy of saturated water kJ/kg
h’’ Specific enthalpy of saturated steam kJ/kg
ρ Density kg/m³
All state functions FNAME return a state variable ZVAL as an analog value. This is called up
by
ZVAL = _WaterSteam.FNAME(PARAM1 (, PARAM2)).
The available state functions are described in Table 5-6.
ρ rph p h ρ = ρ( p , h )
T trh ρ h T = T (ρ , h )
TS tvp pS - TS = T ( p S )
pS pvt TS - p S = p(TS )
ρ’ rsvp pS - ρ′ = ρ' ( p S )
The saturated pressure p S and saturated temperature T S are limited for the saturation
function as follows:
namespace Userdefined
{
public class FlownetFunctions
{
public static double trh(long m, double r, double h)
{
switch (m)
{
case -1:
// Calculations for media key number -1 …
return 0.0;
case -2:
// Calculations for media key number -2 …
return 0.0;
default:
// Should not occur!
return 0.0;
}
}
}
}
To assign a specific own defined media to a flownet simply set the appropriate negative
media key number for the signal MEDIUM of the connection type FLN5.
CAUTION
If you set a negative key number for the medium without providing the
necessary functions for the calculation of the density and temperature, the
flownet solver assignes a zero value to all densities and temeratures in the
flownet!
NOTE
As there are no variables available of the flownet solver during initialisation of the
components, it is advisable to pre-assign the inputs of the component types with suitable
values. Values that are consistent with the default values of the flownet solver are
recommended.