BuildSysPro.IBPSA.Airflow.Multizone.DoorOpen

# .BuildSysPro.IBPSA.Airflow.Multizone.DoorOpen

Door model for bi-directional air flow between rooms

## Information

Model for bi-directional air flow through a large opening such as a door.

In this model, the air flow is composed of two components, a one-directional bulk air flow due to static pressure difference in the adjoining two thermal zones, and a two-directional airflow due to temperature-induced differences in density of the air in the two thermal zones. Although turbulent air flow is a nonlinear phenomenon, the model is based on the simplifying assumption that these two air flow rates can be superposed. (Superposition is only exact for laminar flow.) This assumption is made because it leads to a simple model and because there is significant uncertainty and assumptions anyway in such simplified a model for bidirectional flow through a door.

#### Main equations

The air flow rate due to static pressure difference is

V̇ab,p = CD w h (2/ρ0)0.5 Δpm,

where V̇ is the volumetric air flow rate, CD is the discharge coefficient, w and h are the width and height of the opening, ρ0 is the mass density at the medium default pressure, temperature and humidity, m is the flow exponent and Δp = pa - pb is the static pressure difference between the thermal zones. For this model explanation, we will assume pa > pb. For turbulent flow, m=1/2 and for laminar flow m=1.

The air flow rate due to temperature difference in the thermal zones is V̇ab,t for flow from thermal zone a to b, and V̇ba,t for air flow rate from thermal zone b to a. The model has two air flow paths to allow bi-directional air flow. The mass flow rates at these two air flow paths are

ṁa1 = ρ0 (+V̇ab,p/2 + V̇ab,t),

and, similarly,

V̇ba = ρ0 (-V̇ab,p/2 + V̇ba,t),

where we simplified the calculation by using the density ρ0. To calculate V̇ba,t, we again use the density ρ0 and because of this simplification, we can write

ṁab,t = -ṁba,t = ρ0 V̇ab,t = -ρ0 V̇ba,t,

from which follows that the neutral height, e.g., the height where the air flow rate due to flow induced by temperature difference is zero, is at h/2. Hence,

V̇ab,t = CD ∫0h/2 w v(z) dz,

where v(z) is the velocity at height z. From the Bernoulli equation, we obtain

v(z) = (2 g z Δρ ⁄ ρ0)1/2.

The density difference can be written as

Δρ = ρa-ρb ≈ ρ0 (Tb - Ta) ⁄ T0,

where we used ρa = p0 /(R Ta) and Ta Tb ≈ T02. Substituting this expression into the integral and integrating from 0 to z yields

V̇ab,t = 1⁄3 CD w h (g h ⁄ (R T0 ρ0))1/2 Δp1/2.

The above equation is equivalent to (6) in Brown and Solvason (1962).

#### Main assumptions

The main assumptions are as follows:

- The air flow rates due to static pressure difference and due to temperature-difference can be superposed.

- For buoyancy-driven air flow, a constant density can be used to convert air volume flow rate to air mass flow rate.

From these assumptions follows that the neutral height for buoyancy-driven air flow is at half of the height of the opening.

#### Notes

For a more detailed model, use IBPSA.Airflow.Multizone.DoorDiscretizedOpen.

#### References

- Brown, W.G. and K. R. Solvason. Natural Convection through rectangular openings in partitions - 1. Int. Journal of Heat and Mass Transfer. Vol. 5, p. 859-868. 1962. doi:10.1016/0017-9310(62)90184-9.
  Also available at https://nrc-publications.canada.ca/eng/view/ft/?id=081c0ace-7c31-449c-9b3b-e6c14864b196.

## Revisions

- January 22, 2020, by Michael Wetter:
  Revised buoyancy-driven flow based on Brown and Solvason (1962).
- January 19, 2020, by Klaas De Jonge:
  Revised influence of stack effect.
- October 6, 2020, by Michael Wetter:
  First implementation for #1353.

Generated at 2026-05-09T18:18:56Z by OpenModelicaOpenModelica 1.26.7 using GenerateDoc.mos