This model is an air damper with flow coefficient that is an exponential function
of the opening angle. The model is as in ASHRAE 825-RP.
A control signal of y=0
means the damper is closed, and y=1
means the damper
is open. This is opposite of the implementation of ASHRAE 825-RP, but used here
for consistency within this library.
For yL < y < yU
, the damper characteristics is
kd(y) = exp(a+b (1-y)).
Outside this range, the damper characteristic is defined by a quadratic polynomial that
matches the damper resistance at y=0
and y=yL
or y=yU
and
y=1
, respectively. In addition, the polynomials are such that
kd(y) is
differentiable in y and the derivative is continuous.
The damper characteristics kd(y) is then used to compute the flow coefficient k(y) as
k(y) = (2 ρ ⁄ kd(y))1/2 A,
where A is the face area, which is computed using the nominal
mass flow rate m_flow_nominal
, the nominal velocity
v_nominal
and the density of the medium. The flow coefficient k(y)
is used to compute the mass flow rate versus pressure
drop relation as
m = sign(Δp) k(y) √ Δp
with regularization near the origin.
ASHRAE 825-RP lists the following parameter values as typical:
opposed blades | single blades | |
---|---|---|
yL | 15/90 | 15/90 |
yU | 55/90 | 65/90 |
k0 | 1E6 | 1E6 |
k1 | 0.2 to 0.5 | 0.2 to 0.5 |
a | -1.51 | -1.51 |
b | 0.105*90 | 0.0842*90 |
P. Haves, L. K. Norford, M. DeSimone and L. Mei, A Standard Simulation Testbed for the Evaluation of Control Algorithms & Strategies, ASHRAE Final Report 825-RP, Atlanta, GA.
kDam_default
and kThetaSqRt_default
from initial algorithm
to the variable declaration, to avoid a division
by zero in OpenModelica.