modelDiversionOpenLoop
Extends from BaseClasses.PartialActivePrimary (Partial model of active primary network).
Information
This model represents a heating system where the configuration Buildings.Fluid.HydronicConfigurations.ActiveNetworks.Diversion is used to modulate the heat flow rate transmitted to a constant load. Two identical secondary circuits are connected to a primary circuit with a constant speed pump. The main assumptions are enumerated below.
- The model is configured in steady-state.
-
The design conditions at
time=0are defined without considering any load diversity. - Each circuit is balanced at design conditions.
-
The bypass of the diversion circuit is balanced at
design conditions if the parameter
is_balis set totrue. Otherwise no fixed flow resistance is considered in the bypass branch, only the variable flow resistance corresponding to the bypass port of the three-way valve.
When simulated with the default parameter values, this example shows the following points.
-
The overflow caused by the unbalanced bypass when the valve
is fully closed in the first circuit (see plot #2 at
time=100) creates a concomitant flow shortage in the second circuit with the valve fully open. However, the flow shortage (4%) is of a much lower amplitude than the overflow (30%). Indeed the equivalent flow resistance seen by the pump is lower than at design conditions, leading a shift of the operating point of the pump towards a flow rate value higher than design, which partly compensates for the overflow. - The impact on the heat flow rate transferred to the load (see plot #4) is of an even lower amplitude (2%). This is dependent on the emission characteristic of the terminal unit that would theoretically be steeper with a higher effectiveness at design conditions, although the order of magnitude would remain the same. For reference, the terminal unit in this example represents a recirculating air terminal such as a fan coil unit in heating mode with a water ΔT of 10 °C at design conditions, so an effectiveness ε = ΔTliq / (Tliq, inl - Tair, inl) = 0.25.
- The equal-percentage / linear characteristic of the control valve yields a relationship between the heat flow rate transferred to the load and the valve opening that is close to linear (see plot #4), with a Pearson correlation coefficient equal to 0.99.
Sensitivity analysis
Those observations are confirmed by a sensitivity study to the following parameters.
- Ratio of the terminal unit pressure drop to the pump head at design conditions (refer to the schematic in the documentation of Buildings.Fluid.HydronicConfigurations.ActiveNetworks.Diversion for the nomenclature): ψ = ΔpJ-A / Δppump varying from 0.1 to 0.4
- Ratio of the control valve authority: β = ΔpA-AB / ΔpJ-AB varying from 0.1 to 0.7
-
Balanced bypass branch:
is_balswitched fromfalsetotrue -
Valve characteristic:
ThreeWayValveswitched from equal percentage-linear (EL) to linear-linear (LL).
Direct and bypass mass flow rate
The overflow in the bypass branch when the valve is fully closed increases with ψ and decreases with β. It is close to 70% for ψ = 40% and β = 10%. However, the concomitant flow shortage in the other terminal unit with a valve fully open (see Figure 2) is limited to 12%. For a valve authority of β = 50% one may note that the flow shortage is below 5%, indicating that selecting the control valve with a suitable authority largely dampens the impact of an unbalanced bypass.
Figure 1. Bypass mass flow rate (ratio to design value) at fully closed conditions
as a function of
ψ for various valve authorities β (color scale),
and a bypass branch either balanced (right plot) or not (left plot).
Figure 2. Direct mass flow rate (ratio to design value) at fully open conditions
as a function of
ψ for various valve authorities β (color scale),
and a bypass branch either balanced (right plot) or not (left plot).
Primary mass flow rate
The total primary mass flow rate (or pump mass flow rate) is plotted on Figure 3. This helps assess the actual flow variation in "constant flow circuits", i.e., constant speed pump distribution systems with terminal units equipped with three-way control valves. The total pump flow can vary up to 50% when the bypass of the three-way valves is not balanced, whereas the flow variation is limited to about 20% when the bypass of the three-way valves is balanced. If the characteristic of the valves is equal percentage and linear, this variation is rather by higher values for an unbalanced bypass and by lower values for a balanced bypass. If the characteristic of the valves is linear and linear, this variation is always by higher values. Eventually, when the control valve authority is higher than 0.5 the flow variation is limited to about ±20% in all cases.
Figure 3. Pump mass flow rate (ratio to design value) as a function of
ψ for various valve authorities β (color scale),
a bypass branch either balanced (right plots) or not (left plots)
and either an equal-percentage / linear valve characteristic (top plots)
or a linear / linear valve characteristic (bottom plots).
Heat flow rate transferred to the load
The heat flow rate transferred to the load is presented
- on Figure 4 at full opening to illustrate the impact on the coil capacity of the flow shortage previously discussed,
- on Figure 5 at 10% opening to illustrate the linearity of the relationship between the transmitted heat flow rate and the valve opening at low load in the case of an equal percentage valve (where values of the heat flow rate close to 10% of the coil capacity indicate a good linearity).
The impact on the coil capacity of the flow shortage due to an unbalanced bypass is limited to about 5% and is less than 2% for an authority higher or equal to 0.5. A balanced bypass tends to disturb the linearity of the heat flow rate with the valve opening. But again, if the valve is selected with an authority higher or equal to 0.5 that disturbance is highly reduced.
Figure 4. Heat flow rate (ratio to design value) at fully open conditions
as a function of
ψ for various valve authorities β (color scale),
and a bypass branch either balanced (right plot) or not (left plot).
Figure 5. Heat flow rate (ratio to design value) at 10% open conditions
as a function of ψ for various valve authorities β (color scale),
and a bypass branch either balanced (right plot) or not (left plot).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Buildings.Fluid.HydronicConfigurations.Types.Control | typ (from PartialActivePrimary) | Buildings.Fluid.HydronicConfigurations.Types.Control.Heating | Load type |
| Integer | nTer (from PartialActivePrimary) | 2 | Number of terminal units |
| Real | kSizPum (from PartialActivePrimary) | 1.0 | Pump oversizing coefficient |
| Modelica.Units.SI.Pressure | p_min (from PartialActivePrimary) | 200000 | Circuit minimum pressure |
| Modelica.Units.SI.Temperature | TLiqEnt_nominal (from PartialActivePrimary) | if typ == Buildings.Fluid.HydronicConfigurations.Types.Control.Heating then 60 + 273.15 else 7 + 273.15 | Liquid entering temperature at design conditions |
| Modelica.Units.SI.Temperature | TLiqLvg_nominal (from PartialActivePrimary) | TLiqEnt_nominal + (if typ == Buildings.Fluid.HydronicConfigurations.Types.Control.Heating then -10 else +5) | Liquid leaving temperature at design conditions |
| Modelica.Units.SI.Temperature | TLiqEntChg_nominal (from PartialActivePrimary) | 60 + 273.15 | Liquid entering temperature in change-over mode |
| Modelica.Units.SI.Temperature | TLiqSup_nominal (from PartialActivePrimary) | TLiqEnt_nominal | Liquid primary supply temperature at design conditions |
| Modelica.Units.SI.Temperature | TLiqSupChg_nominal (from PartialActivePrimary) | TLiqEntChg_nominal | Liquid primary supply temperature in change-over mode |
| Modelica.Units.SI.PressureDifference | dpValve_nominal | dpTer_nominal | Control valve pressure drop at design conditions |
| Nominal condition | |||
| Modelica.Units.SI.MassFlowRate | mTer_flow_nominal (from PartialActivePrimary) | 1 | Terminal unit mass flow rate at design conditions |
| Modelica.Units.SI.MassFlowRate | m1_flow_nominal (from PartialActivePrimary) | m2_flow_nominal | Mass flow rate in primary branch at design conditions |
| Modelica.Units.SI.MassFlowRate | m2_flow_nominal (from PartialActivePrimary) | nTer*mTer_flow_nominal | Mass flow rate in consumer circuit at design conditions |
| Modelica.Units.SI.PressureDifference | dpTer_nominal (from PartialActivePrimary) | 3E4 | Terminal unit pressure drop at design conditions |
| Modelica.Units.SI.PressureDifference | dpPip_nominal (from PartialActivePrimary) | 0.5E4 | Pipe section pressure drop at design conditions |
| Modelica.Units.SI.PressureDifference | dpPum_nominal (from PartialActivePrimary) | Pump head at design conditions | |
| Modelica.Units.SI.MassFlowRate | mPum_flow_nominal (from PartialActivePrimary) | m1_flow_nominal | Primary pump mass flow rate at design conditions |
| Modelica.Units.SI.PressureDifference | dpBal3_nominal | if is_bal then dpTer_nominal else 0 | Bypass balancing valve pressure drop at design conditions |
| Dynamics › Conservation equations | |||
| Modelica.Fluid.Types.Dynamics | energyDynamics (from PartialActivePrimary) | Modelica.Fluid.Types.Dynamics.FixedInitial | Type of energy balance: dynamic (3 initialization options) or steady state |
| Configuration | |||
| Boolean | is_bal | false | Set to true for a balanced bypass |
| Buildings.Fluid.HydronicConfigurations.Types.ValveCharacteristic | typCha | Buildings.Fluid.HydronicConfigurations.Types.ValveCharacteristic.EqualPercentage | Control valve characteristic |
Components
Revisions
-
June 30, 2022, by Antoine Gautier:
First implementation.