自然換気量の計算法 [多重開口の場合]

建物の開口部で生じる自然換気量は、下記の式[1]を用いて計算される。

Q = \alpha A\sqrt{\frac{2 \Delta P}{\rho}} \; \cdots \; [1]

ここで、Q [m³/s]は換気量、α [-]は流量係数、A [m²]は開口面積、ρ [kg/m³]は流体密度（空気の場合は、約1.2 kg/m³）、ΔP [Pa]は室内外圧力差である。

式[1]によると自然換気量は、室内外圧力差 ΔP（風圧差 ΔP_w か温度差 ΔT ）と、開口における流量係数 α（圧力損失係数 ξの平方根の逆数）に比例する形で計算される。以下では、流量係数 α、有効開口面積 αAについてそれぞれ概説し、風圧差 ΔP_wか温度差 ΔT によって発生する自然換気（風力換気／温度差換気）による自然換気量の計算法について説明する。

※ 多重開口ではなく、単一開口における自然換気量の計算法は計算原理が異なるため、下記のブログをご確認ください。
https://lee-lab.net/blog-contents-004

【1】流量係数（α）

空気が開口部を通過する際には、開口部の形状、材質、長さ、空気の流速などの要因により抵抗が生じ、その結果として空気の圧力が低下する。この現象は「圧力損失」と呼ばれる。この圧力損失の大きさを開口部を通過する気流速度（U_o）や密度（ρ）に依存しない形で比較するために、圧力損失（ΔP）を動圧（P_v）で割った無次元数が一般的に用いられる。これを「圧力損失係数」と呼び、式[2]に示すように 記号 ξ [-]で表される。

\xi = \frac{{\Delta}{P}}{P_{v}} = \frac{{\Delta}{P}}{ \frac{{1}}{{2}} \rho U_{o}^{2}} = \frac{{\Delta}{P}}{ \frac{{1}}{{2}} \rho \left ( {Q_{r}}/{A} \right )^{2}} \; \cdots \; [2]

ここで、P_v [Pa]は動圧、U_o [m/s]は開口部を通過する気流速度（代表速度）、Q_r [m³/s]は実際の流量または実測流量である。

一方、流量係数 αは、開口部を通過する実際の流量と、圧力損失がない理想状態で得られる理論流量と比（式[3]参照）として定義され、圧力損失係数 ξの平方根の逆数に相当する。

\alpha = \frac{{Q}_{r}}{A \sqrt {\frac{{2}\Delta{P}}{\rho}}} = \frac{{1}}{{\sqrt{\xi}}} \; \cdots \; [3]

図1に、様々な開口部における典型的な流量係数を示す。例えば、通常の窓は0.6～0.7程度であるが、ベルマウス形状のように入口形状がスムーズな場合、流れの縮流や剥離が抑えられ、流量係数は概ね1となる。また、ルーバーが設置される場合はその角度によって圧力損失が変わるため、角度によって異なる。開き窓の場合、開き角度によって変わる流量係数を図2に示す。また、網戸、カーテン、ブラインドなどを使用する場合にも圧力損失が生じるため、適正な流量係数を計測して適用する必要がある。

図1. 開口部における流量係数

図2. 窓の開き角度による流量係数
(Calculated by Akane Tsutsumi et al. [PDF])

【2】有効開口面積（αA）

有効開口面積（αA）は開口配置によって並列合成と直列合成に分けて計算でき、それぞれ式[7]と式[11]で計算される。

① αA の並列合成

有効開口面積（αA）の並列合成は図3に示したものであり、開口部1と2の前後に作用する全圧P_Ta とP_Tb が等しいとしてαA_合値を計算する。

図3. αAの並列合成

Q_{1} = \alpha_{1} A_{1} \sqrt{ \frac{{2}{\Delta}{P}}{\rho} } \; \cdots \; [4]

Q_{2} = \alpha_{2} A_{2} \sqrt{ \frac{{2}{\Delta}{P}}{\rho} } \; \cdots \; [5]

\def\arraystretch{1.2}\begin{aligned} Q & = Q_{1}+Q_{2} \\\ &= (\alpha_{1} A_{1}+\alpha_{2} A_{2}) \sqrt{ \frac{{2}{\Delta}{P}}{\rho} }
\; \cdots \; [6]
\end{aligned} 

\alpha A_{合} = \alpha_{1} A_{1} + \alpha_{2} A_{2}  \; \cdots \; [7]

② αA の直列合成

有効開口面積（αA）の直列合成は図4に示したものであり、開口部1と2を通過する風量 Q が等しいとして αA_合値を計算する。

図4. αA の直列合成

Q = \alpha_{1} A_{1} \sqrt{ \frac{{2}{\Delta}{P_{a} }}{\rho} } = \alpha_{2} A_{2} \sqrt{ \frac{{2}{\Delta}{P_{b} }}{\rho} } \; \cdots \; [8]

\def\arraystretch{1.5}\begin{aligned}
\Delta P & = \Delta P_{a} + \Delta P_{b} \\
& = \frac{1}{2} \rho \left(\frac{Q}{\alpha_{1} A_{1}} \right)^{2} +  \frac{1}{2} \rho \left(\frac{Q}{\alpha_{2} A_{2}} \right)^{2}
\; \cdots \; [9]
\end{aligned}

Q =\frac{1}{\sqrt {\left(\frac{1}{\alpha_{1} A_{1}} \right)^{2} + \left(\frac{1}{\alpha_{2} A_{2}} \right)^{2}}} \sqrt{ \frac{{2}{\Delta}{P}}{\rho} }  \; \cdots \; [10]

\alpha A_{合} = \frac{1}{\sqrt {\left(\frac{1}{\alpha_{1} A_{1}} \right)^{2} + \left(\frac{1}{\alpha_{2} A_{2}} \right)^{2}}} \; \cdots \; [11]

【3】風力換気

建物に風が吹き付けると、大気圧に風圧力が加算されて作用する。図5に示すように風上側に掛かって来る風圧 P_w [Pa]は、基準風速 （軒高風速が基準）U [m/s]と風圧係数 C_w [-]を用いて式[12]で計算される。

図5. 風力による圧力差

P_{\mathrm{w}} = C_{\mathrm{w}}\frac{1}{2}\rho_{\mathrm{o}}U^{2} \; \cdots \; [12]

ここで、P_w [Pa]は風圧、C_w [-]は風圧係数、ρ_o [kg/m³]は外気密度、U [m/s]は基準風速である。

開口部がある建物に風が吹き付けると、図6に示すように風上側に掛かって来る風圧 P_w1 [Pa]と風下側の風圧 P_w2 [Pa]の差 ΔP [Pa]は、基準風速 U [m/s]と風圧係数 C_w1 [-]と C_w2 [-]の差 (C_w1 – C_w2) [-]を用いて式[13]で計算できる。

図6. 風力による圧力差

\Delta P = \left ( C_{\mathrm{w1}}-C_{\mathrm{w2}} \right )\frac{1}{2}\rho_{\mathrm{o}}U^{2} \; \cdots \; [13]

ここで、ΔP [Pa]は風圧P_w1 [Pa]と風圧P_w2 [Pa]の差圧、C_w1 [-]は風上側の風圧係数、C_w2 [-]は風下側の風圧係数、ρ_o [kg/m³]は外気密度、U [m/s]は基準風速である。

式[13]を式[1]に代入して整理すると、風力による換気量 Qは式[14]で計算できる。

Q = \alpha A_{合} \cdot U \cdot \sqrt{C_{\mathrm{w1}}-C_{\mathrm{w2}} } \; \cdots \; [14]

一方、風圧係数は実測、模型を用いた風洞実験（図7参照）によって測定値を用いるか、それとも推定式、CFD解析などによる計算値を使う場合もある。図8、図9にAkinsら（1979）¹[1] Akins, R.E., J.A. Peterka, and J.E. Cermak. 1979. Averaged pressure coefficients for rectangular buildings. Wind Engineering: Proceedings of the Fifth International Conference, vol. 7, pp. 369-380.が高層ビルを対象として行った風洞実験の壁面、屋根面の結果値を示し、図10には低層ビルに対してSwamiら（1987）²[2] Swami, M.V., and S. Chandra. 1987. Procedures for calculating natural ventilation airflow rates in buildings. Final Report FSEC-CR-163-86. Florida Solar Energy Center, Cape Canaveral.が報告した風圧係数を示す。

図7. 模型を用いた風洞実験による風圧係数

図8. 壁面における風圧係数（Akins et al.）

図9. 屋根面における風圧係数（Akins et al.）

図10. 壁面における風圧係数（Swami et al.）

また、Walkerら（1994）³[3] Walker, I.S., and D.J. Wilson. 1994. Practical methods for improving estimates of natural ventilation rates. Proceedings of the 15th IEA Conference of the Air Infiltration and Ventilation Centre, Buxton, U.K., pp. 517-526.は、風の風向角による風圧係数の推定式を式[15]のように提案している。

\def\arraystretch{1.5}\begin{aligned}
{C}_{\textrm{p}}(\phi) = & 1/2 \{ [{{C}_{\textrm{p}}(1)}+{{C}_{\textrm{p}}(2)}](cos^{2}\phi)^{1/4} \\\ 
& +
[{{C}_{\textrm{p}}(1)}-{{C}_{\textrm{p}}(2)}](cos\phi)^{3/4} \\\ 
& +
[{{C}_{\textrm{p}}(3)}+{{C}_{\textrm{p}}(4)}](sin^{2}\phi)^{2} \\\ 
& +
[{{C}_{\textrm{p}}(3)}-{{C}_{\textrm{p}}(4)}]sin\phi \}
\end{aligned}
 \; \cdots \; [15]

ここで、C_p(1) は風向0°での風圧係数、C_p(2) は風向180°での風圧係数、C_p(3) は風向90°での風圧係数、C_p(4) は風向270°での風圧係数、φ は壁面に対して時計回り風向角である。

【4】温度差換気

図11に示すように室内上下に開口部があり、室内温度が室外温度より高い場合には、室内は上昇気流が生じる。これを温度差換気と言い、その換気量は以下の式[16]で室内外圧力差を求め、式[1]に代入することで計算できる。

\Delta P = \Delta \rho g h  \; \cdots \; [16]

ここで、ΔP [Pa]は開口部1における風圧P₁ [Pa]と開口部2における風圧P₂ [Pa]の差圧、Δρ [kg/m³]は室外の空気密度ρ_o [kg/m³]と室内の空気密度ρ_i [kg/m³]の密度差、g [m/s²]は重力加速度、h [m]は開口間の高さである。

図11. 温度差による圧力差

式[16]における室内外密度差Δρは、式[17]に示すシャルル法則（Charles’ law）によって室外温度T_i [K]と室内外温度差ΔT [K]に変換でき、式[18]のように表せられる。

\rho_{\mathrm{o}} T_{\mathrm{o}} =  \rho_{\mathrm{i}} T_{\mathrm{i}} \; \cdots \; [17]

\begin{aligned}
\Delta \rho & = \rho_{\mathrm{o}} - \rho_{\mathrm{i}} = \rho_{\mathrm{o}} - \frac{\rho_{\mathrm{o}}T_{\mathrm{o}}}{T_{\mathrm{i}}} \\\ 
& = \frac{\rho_{\mathrm{o}}(T_{\mathrm{i}} - T_{\mathrm{o}})}{T_{\mathrm{i}}}  = \frac{\rho_{\mathrm{o}}\Delta T}{T_{\mathrm{i}}} \; \cdots \; [18]
\end{aligned}

ここで、ΔT [K]は室内温度T_i [K]と室外温度T_o [K]の室内外温度差である。

式[16]に式[18]を代入して整理すると、式[19]のように表せられる。

\Delta P = \Delta \rho g h  = \frac{\rho_{\mathrm{o}}\Delta Tgh}{T_{\mathrm{i}}} \; \cdots \; [19]

式[19]を式[1]に代入して整理すると、温度差による換気量 Qは式[20]で計算できる。

Q = \alpha A_{合} \cdot \sqrt{\frac{2 \cdot \Delta Tgh}{T_{\mathrm{i}}} } \; \cdots \; [20]

※ 開口部が一箇所の建物では開口部での大気基準圧は0となり、それ以上では正圧（図12参照）となる。ここで大気基準圧とは、同一高度の大気圧を基準とした圧力表示を大気基準圧と言う。

図12. 開口部が一箇所の建物

また、上部にも開口部ができると、上部では正圧、下部では負圧となる分布が形成され、上部で室空気が流出し、下部から流入（図13参照）する。

図13. 開口部が二箇所の建物

一方、中性帯の高さは、開口部の有効開口面積（αA）が同じであれば中立位置に形成されるが、開口部の有効開口面積が異なるとαAが大きい方（図14参照）に近づく。

図14. 開口部の有効開口面積（αA）と中性帯高さの関係

★ 以上の自然換気量の計算法に基づき、例題を解くことで深く学習したい方は下記のブログをご確認ください。
https://lee-lab.net/blog-contents-005

[1] Akins, R.E., J.A. Peterka, and J.E. Cermak. 1979. Averaged pressure coefficients for rectangular buildings. Wind Engineering: Proceedings of the Fifth International Conference, vol. 7, pp. 369-380.【PDF LINK】
[2] Swami, M.V., and S. Chandra. 1987. Procedures for calculating natural ventilation airflow rates in buildings. Final Report FSEC-CR-163-86. Florida Solar Energy Center, Cape Canaveral.【PDF LINK】
[3] Walker, I.S., and D.J. Wilson. 1994. Practical methods for improving estimates of natural ventilation rates. Proceedings of the 15th IEA Conference of the Air Infiltration and Ventilation Centre, Buxton, U.K., pp. 517-526.【PDF LINK】

Written by Sihwan Lee
[Associate Professor, Tokyo University of Science]