modelPressureIndependent

Model for an air damper whose mass flow is proportional to the input signal

Extends from Buildings.Fluid.Actuators.BaseClasses.PartialDamperExponential (Partial model for air dampers with exponential opening characteristics).

Information

Model for an air damper whose airflow is proportional to the input signal, assuming that at y = 1, m_flow = m_flow_nominal. This is unless the pressure difference dp is too low, in which case a kDam = m_flow_nominal/sqrt(dp_nominal) characteristic is used.

The model is similar to Buildings.Fluid.Actuators.Valves.TwoWayPressureIndependent, except for adaptations for damper parameters. Please see that documentation for more information.

Computation of the damper opening

The fractional opening of the damper is computed by

  • inverting the quadratic flow function to compute the flow coefficient from the flow rate and the pressure drop values (under the assumption of a turbulent flow regime);
  • inverting the exponential characteristics to compute the fractional opening from the loss coefficient value (directly derived from the flow coefficient).

The quadratic interpolation used outside the exponential domain in the function Buildings.Fluid.Actuators.BaseClasses.exponentialDamper yields a local extremum. Therefore, the formal inversion of the function is not possible. A cubic spline is used instead to fit the inverse of the damper characteristics. The central domain of the characteritics having a monotonous exponential profile, its inverse can be properly approximated with three equidistant support points. However, the quadratic functions used outside of the exponential domain can have various profiles depending on the damper coefficients. Therefore, five linearly distributed support points are used on each side domain to ensure a good fit of the inverse.

Note that below a threshold value of the input control signal (fixed at 0.02), the fractional opening is forced to zero and no more related to the actual flow coefficient of the damper. This avoids steep transients of the computed opening while transitioning from reverse flow. This is to be considered as a modeling workaround (avoiding the introduction of an additional state variable) to prevent control chattering during shut off operation where the pressure difference at the damper boundaries can vary between slightly positive and negative values due to outdoor pressure variations.

Parameters

TypeNameDefaultDescription
BooleanhomotopyInitialization (from PartialResistance)true= true, use homotopy method
Modelica.Units.SI.MassFlowRatem_flow_turbulent (from PartialResistance)Turbulent flow if |m_flow| >= m_flow_turbulent
Booleanuse_deltaM (from PartialDamperExponential)trueSet to true to use deltaM for turbulent transition, else ReC is used
RealdeltaM (from PartialDamperExponential)0.3Fraction of nominal mass flow rate where transition to turbulent occurs
Modelica.Units.SI.Velocityv_nominal (from PartialDamperExponential)(2/rho_default/k1*dpDamper_nominal)^0.5Nominal face velocity
Modelica.Units.SI.AreaA (from PartialDamperExponential)m_flow_nominal/rho_default/v_nominalFace area
BooleanroundDuct (from PartialDamperExponential)falseSet to true for round duct, false for square cross section
RealReC (from PartialDamperExponential)4000Reynolds number where transition to turbulence starts
RealkFixed (from PartialDamperExponential)if dpFixed_nominal > Modelica.Constants.eps then m_flow_nominal/sqrt(dpFixed_nominal) else Modelica.Constants.infFlow coefficient of fixed resistance that may be in series with damper, k=m_flow/sqrt(dp), with unit=(kg.m)^(1/2).
Assumptions
BooleanallowFlowReversal (from PartialTwoPort)true= false to simplify equations, assuming, but not enforcing, no flow reversal
Nominal condition
Modelica.Units.SI.MassFlowRatem_flow_nominal (from PartialTwoPortInterface)Nominal mass flow rate
Modelica.Units.SI.PressureDifferencedp_nominal (from PartialResistance)Pressure drop at nominal mass flow rate
Modelica.Units.SI.PressureDifferencedpDamper_nominal (from PartialDamperExponential)Pressure drop of fully open damper at nominal mass flow rate
Modelica.Units.SI.PressureDifferencedpFixed_nominal (from PartialDamperExponential)0Pressure drop of duct and resistances other than the damper in series, at nominal mass flow rate
Advanced
Modelica.Units.SI.MassFlowRatem_flow_small (from PartialTwoPortInterface)1E-4*abs(m_flow_nominal)Small mass flow rate for regularization of zero flow
Booleanfrom_dp (from PartialResistance)false= true, use m_flow = f(dp) else dp = f(m_flow)
Realn (from PartialResistance)2Flow exponent, n=1 for laminar, n=2 for turbulent
Booleanlinearized (from PartialResistance)false= true, use linear relation between m_flow and dp for any flow rate
Booleanuse_constant_density (from PartialDamperExponential)trueSet to true to use constant density for flow friction
Reall20.01Gain for mass flow increase if pressure is above nominal pressure
Realdeltax0.02Transition interval for flow rate
Advanced › Diagnostics
Booleanshow_T (from PartialTwoPortInterface)false= true, if actual temperature at port is computed
Dynamics › Actuator position
Booleanuse_strokeTime (from ActuatorSignal)trueSet to true to continuously open and close valve using strokeTime
Modelica.Units.SI.TimestrokeTime (from ActuatorSignal)120Time needed to fully open or close actuator
Modelica.Blocks.Types.Initinit (from ActuatorSignal)Modelica.Blocks.Types.Init.InitialOutputType of initialization (no init/steady state/initial state/initial output)
Realy_start (from ActuatorSignal)1Initial position of actuator
Damper coefficients
Reala (from PartialDamperExponential)-1.51Coefficient a for damper characteristics
Realb (from PartialDamperExponential)0.105*90Coefficient b for damper characteristics
RealyL (from PartialDamperExponential)15/90Lower value for damper curve
RealyU (from PartialDamperExponential)55/90Upper value for damper curve
Realk0 (from PartialDamperExponential)2*rho_default*(A/kDamMin)^2Loss coefficient for y=0 (pressure drop divided by dynamic pressure)
Realk1 (from PartialDamperExponential)0.45Loss coefficient for y=1 (pressure drop divided by dynamic pressure)
Reall (from PartialDamperExponential)0.0001Damper leakage, ratio of flow coefficients k(y=0)/k(y=1)

Connectors

TypeNameDefaultDescription
Modelica.Fluid.Interfaces.FluidPort_aport_a (from PartialTwoPort)Fluid connector a (positive design flow direction is from port_a to port_b)
Modelica.Fluid.Interfaces.FluidPort_bport_b (from PartialTwoPort)Fluid connector b (positive design flow direction is from port_a to port_b)
Modelica.Blocks.Interfaces.RealInputy (from ActuatorSignal)Actuator position (0: closed, 1: open)
Modelica.Blocks.Interfaces.RealOutputy_actual (from ActuatorSignal)Actual actuator position

Components

TypeNameDefaultDescription
Modelica.Units.SI.MassFlowRatem_flow (from PartialTwoPortInterface)port_a.m_flowMass flow rate from port_a to port_b (m_flow > 0 is design flow direction)
Modelica.Units.SI.PressureDifferencedp (from PartialTwoPortInterface)port_a.p - port_b.pPressure difference between port_a and port_b
Medium.ThermodynamicStatesta_a (from PartialTwoPortInterface)if allowFlowReversal then Medium.setState_phX(port_a.p, noEvent(actualStream(port_a.h_outflow)), noEvent(actualStream(port_a.Xi_outflow))) else Medium.setState_phX(port_a.p, noEvent(inStream(port_a.h_outflow)), noEvent(inStream(port_a.Xi_outflow)))Medium properties in port_a
Medium.ThermodynamicStatesta_b (from PartialTwoPortInterface)if allowFlowReversal then Medium.setState_phX(port_b.p, noEvent(actualStream(port_b.h_outflow)), noEvent(actualStream(port_b.Xi_outflow))) else Medium.setState_phX(port_b.p, noEvent(port_b.h_outflow), noEvent(port_b.Xi_outflow))Medium properties in port_b
Medium.Densityrho (from PartialDamperExponential)Medium density
RealkDam (from PartialDamperExponential)Flow coefficient of damper, k=m_flow/sqrt(dp), with unit=(kg.m)^(1/2)
Realk (from PartialDamperExponential)Flow coefficient of damper plus fixed resistance, k=m_flow/sqrt(dp), with unit=(kg.m)^(1/2)
Realphil + y_internal*(1 - l)Ratio actual to nominal mass flow rate of damper, phi=kDam(y)/kDam(y=1)

Contents

NameDescription
basicFlowFunction_dp_m_flowprotectedInverse of flow function that computes that computes the square inverse of flow coefficient
exponentialDamper_invprotectedInverse function of the exponential damper characteristics

Revisions

  • July 30, 2026, by Michael Wetter:
    Corrected argument for computing spline derivatives. Note that older models without this correction will fail to simulate due to an assertion in the spline derivative calculations.
    This is for IBPSA #2130.
  • August 11, 2021, by Michael Wetter:
    Reformulated initial equation section to avoid warning in OPTIMICA about variable array index.
    This is for IBPSA #1513.
  • June 10, 2021, by Michael Wetter:
    Changed implementation of the filter and changed the parameter order to a constant as most users need not change this value.
    This is for IBPSA #1498.
  • April 6, 2020, by Antoine Gautier:
    Added the computation of the damper opening.
  • December 23, 2019 by Antoine Gautier:
    Refactored as the model can now extend directly Buildings.Fluid.Actuators.BaseClasses.PartialDamperExponential.
    This is for IBPSA #1188.
  • March 21, 2017 by David Blum:
    First implementation.