modelDiversionOpenLoop

Model illustrating the operation of diversion circuits with constant speed pump

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=0 are 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_bal is set to true. 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_bal switched from false to true
  • Valve characteristic: ThreeWayValve switched 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.

Diversion circuit bypass flow rate
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).

Diversion circuit direct flow rate
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.

Diversion circuit pump flow rate
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.

Diversion circuit heat flow rate fully open
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).

Diversion circuit heat flow rate 10% open
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

TypeNameDefaultDescription
Buildings.Fluid.HydronicConfigurations.Types.Controltyp (from PartialActivePrimary)Buildings.Fluid.HydronicConfigurations.Types.Control.HeatingLoad type
IntegernTer (from PartialActivePrimary)2Number of terminal units
RealkSizPum (from PartialActivePrimary)1.0Pump oversizing coefficient
Modelica.Units.SI.Pressurep_min (from PartialActivePrimary)200000Circuit minimum pressure
Modelica.Units.SI.TemperatureTLiqEnt_nominal (from PartialActivePrimary)if typ == Buildings.Fluid.HydronicConfigurations.Types.Control.Heating then 60 + 273.15 else 7 + 273.15Liquid entering temperature at design conditions
Modelica.Units.SI.TemperatureTLiqLvg_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.TemperatureTLiqEntChg_nominal (from PartialActivePrimary)60 + 273.15Liquid entering temperature in change-over mode
Modelica.Units.SI.TemperatureTLiqSup_nominal (from PartialActivePrimary)TLiqEnt_nominalLiquid primary supply temperature at design conditions
Modelica.Units.SI.TemperatureTLiqSupChg_nominal (from PartialActivePrimary)TLiqEntChg_nominalLiquid primary supply temperature in change-over mode
Modelica.Units.SI.PressureDifferencedpValve_nominaldpTer_nominalControl valve pressure drop at design conditions
Nominal condition
Modelica.Units.SI.MassFlowRatemTer_flow_nominal (from PartialActivePrimary)1Terminal unit mass flow rate at design conditions
Modelica.Units.SI.MassFlowRatem1_flow_nominal (from PartialActivePrimary)m2_flow_nominalMass flow rate in primary branch at design conditions
Modelica.Units.SI.MassFlowRatem2_flow_nominal (from PartialActivePrimary)nTer*mTer_flow_nominalMass flow rate in consumer circuit at design conditions
Modelica.Units.SI.PressureDifferencedpTer_nominal (from PartialActivePrimary)3E4Terminal unit pressure drop at design conditions
Modelica.Units.SI.PressureDifferencedpPip_nominal (from PartialActivePrimary)0.5E4Pipe section pressure drop at design conditions
Modelica.Units.SI.PressureDifferencedpPum_nominal (from PartialActivePrimary)Pump head at design conditions
Modelica.Units.SI.MassFlowRatemPum_flow_nominal (from PartialActivePrimary)m1_flow_nominalPrimary pump mass flow rate at design conditions
Modelica.Units.SI.PressureDifferencedpBal3_nominalif is_bal then dpTer_nominal else 0Bypass balancing valve pressure drop at design conditions
Dynamics › Conservation equations
Modelica.Fluid.Types.DynamicsenergyDynamics (from PartialActivePrimary)Modelica.Fluid.Types.Dynamics.FixedInitialType of energy balance: dynamic (3 initialization options) or steady state
Configuration
Booleanis_balfalseSet to true for a balanced bypass
Buildings.Fluid.HydronicConfigurations.Types.ValveCharacteristictypChaBuildings.Fluid.HydronicConfigurations.Types.ValveCharacteristic.EqualPercentageControl valve characteristic

Components

TypeNameDefaultDescription
Sources.Boundary_pTref (from PartialActivePrimary)Pressure and temperature boundary condition
Buildings.Fluid.HydronicConfigurations.Components.Pumppum (from PartialActivePrimary)Circulation pump
FixedResistances.PressureDropres1 (from PartialActivePrimary)Pipe pressure drop
Sensors.TemperatureTwoPortT1Ret (from PartialActivePrimary)Return temperature sensor
Sensors.TemperatureTwoPortT1Sup (from PartialActivePrimary)Supply temperature sensor
Buildings.Controls.OBC.CDL.Reals.SubtractdT1 (from PartialActivePrimary)Primary Delta-T
Delays.DelayFirstOrderdel1 (from PartialActivePrimary)Fluid transport delay
ActiveNetworks.DiversionconHydronic connection
Buildings.Fluid.HydronicConfigurations.ActiveNetworks.Examples.BaseClasses.LoadloaLoad
.Buildings.Controls.OBC.CDL.Reals.Sources.ConstantfraLoaLoad modulating signal
ActiveNetworks.Diversioncon1Hydronic connection
Buildings.Fluid.HydronicConfigurations.ActiveNetworks.Examples.BaseClasses.Loadloa1Load
FixedResistances.PressureDropres1bPipe pressure drop
Sensors.RelativePressuredpDifferential pressure
Sensors.RelativePressuredp1Differential pressure
.Buildings.Controls.OBC.CDL.Reals.Sources.TimeTableopeValve opening signal
Buildings.Controls.OBC.CDL.Integers.Sources.ConstantmodeOperating mode
Buildings.Controls.OBC.CDL.Integers.GreaterThresholdisEnaReturns true if enabled

Revisions

  • June 30, 2022, by Antoine Gautier:
    First implementation.