大気境界層プロファイル計算 戻る
大気境界層 / CFD

大気境界層プロファイル計算

地表粗度を選択し、基準高度の風速を設定するだけで、対数則とべき乗則による風速鉛直プロファイルをリアルタイムに比較。乱流強度・風エネルギー密度も同時算出。

パラメータ設定
地表粗度クラス
基準高度風速 U_ref (m/s) 10.0
風力タービンハブ高度 (m) 80
ハブ高度での計算結果
U_log (m/s)
U_pow (m/s)
乱流強度 I
風力密度 (W/m²)

理論メモ

対数則:$U(z)=\dfrac{u_*}{\kappa}\ln\!\left(\dfrac{z}{z_0}\right)$
$u_* = \dfrac{U_{ref}\,\kappa}{\ln(z_{ref}/z_0)},\quad \kappa=0.41$

べき乗則:$U(z)=U_{ref}\!\left(\dfrac{z}{z_{ref}}\right)^{\!\alpha}$

乱流強度:$I(z)\approx\dfrac{1}{\ln(z/z_0)}$

風力密度:$E=\dfrac{1}{2}\rho U^3,\;\rho=1.225\,\text{kg/m}^3$

地表粗度イメージと風速矢印(矢印の長さが風速に対応)

大気境界層プロファイル計算とは

🧑‍🎓
「大気境界層」って何ですか?風力発電の設計でよく聞くけど、具体的にどうやって風速を予測するんですか?
🎓
ざっくり言うと、地表の摩擦で風速が遅くなる大気の一番下の層だよ。例えば、地面すれすれの風は弱いけど、高い風車のブレードの位置では強い風が吹くよね。このツールでは、上の「地表粗度クラス」を「市街地」に変えてみると、地表付近の風速がガクッと落ちるのがわかるよ。地面が凸凹だと摩擦が大きくなるからね。
🧑‍🎓
え、そうなんですか!でも、予測する式が「対数則」と「べき乗則」の2つあるみたいです。どっちを使えばいいんですか?
🎓
実務では用途で使い分けるんだ。対数則は物理的に厳密で、CFD(数値流体力学)のシミュレーションでよく使う。一方、べき乗則は計算がシンプルで、建築基準法など多くの工学基準で採用されているんだ。このシミュレーターで「基準高度風速」のスライダーを動かすと、2つの法則で描かれる曲線がどう変わるか、すぐに比較できるよ。
🧑‍🎓
「乱流強度」も出てきますね。これって何に影響するんですか?
🎓
風の「揺らぎ」の激しさだね。値が大きいと、風車のブレードに加わる力が変動して疲労寿命が短くなったり、建物が揺れたりする原因になる。例えば、地表粗度を「高層ビル街」に設定すると、低高度で乱流強度が大きくなるのがグラフで確認できるよ。風力発電のサイト選定では、平均風速だけでなくこの乱流強度も超重要なパラメータなんだ。

物理モデルと主要な数式

地表近くの中性大気境界層における風速の鉛直分布を記述する、物理ベースの「対数則」です。摩擦速度 $u_*$ と地表粗度長 $z_0$ が鍵となります。

$$U(z)=\frac{u_*}{\kappa}\ln\!\left(\frac{z}{z_0}\right)$$

$U(z)$: 高度 $z$ における平均風速 [m/s], $u_*$: 摩擦速度 [m/s], $\kappa\approx0.41$: カルマン定数, $z_0$: 地表粗度長 [m]。粗度長 $z_0$ が大きい(地表が粗い)ほど、同じ高度での風速 $U(z)$ は小さくなります。

工学的に広く用いられる経験則である「べき乗則」です。基準高度での風速から、任意の高度の風速を簡便に推定できます。

$$U(z)=U_{ref}\!\left(\frac{z}{z_{ref}}\right)^{\!\alpha}$$

$U_{ref}$: 基準高度 $z_{ref}$ での既知の風速 [m/s], $\alpha$: べき指数。べき指数 $\alpha$ は地表粗度や大気安定度に依存し、例えば平坦な草地では約0.14、市街地では約0.25〜0.4と大きくなります。このツールでは、対数則から換算された $\alpha$ が使われています。

実世界での応用

風力発電のサイトアセスメント:潜在的な設置場所の風況を評価するために使用されます。ハブ高度での平均風速と乱流強度を計算し、年間発電量の予測とタービン機種の選定、構造負荷の評価を行います。

建築物・構造物の耐風設計:建築基準法や各種設計指針では、べき乗則を用いて設計用基準風速を設定高度まで外挿します。特に高層ビルや長大橋など、風荷重が支配的な構造物の設計に不可欠です。

CFD(数値流体力学)シミュレーションの境界条件設定:都市風況や汚染物質拡散の解析では、計算領域の入口境界に現実的な風速プロファイルを与える必要があります。対数則に基づくプロファイルがよく用いられます。

環境アセスメント・大気拡散評価:工場からの排ガスや粉塵の拡散範囲を予測するモデルにおいて、風速の鉛直分布は拡散係数に直接影響する重要な入力パラメータとなります。

よくある誤解と注意点

このツールを使い始めるとき、いくつか気をつけてほしいポイントがあるよ。まず、「地表粗度」の選択はシミュレーションの命だ。例えば、同じ「市街地」でも、低層住宅地と超高層ビルが密集した地域では風の通り方は全く違う。ツールのクラス分けはあくまで代表値だから、実際の現場の写真や土地利用データを見て、最も近いクラスを慎重に選ぶクセをつけよう。間違えると、想定より大幅に風速が低く(または高く)出て、設計が台無しになることもある。

次に、「べき乗則の指数αは万能ではない」という点。このツールでは対数則から換算したαを使っているけど、実はこのαは高度範囲によって少しずつ変化するんだ。例えば、地上10mから100mまでをひとつのαで近似すると、特に地表粗度が大きい場合、誤差が大きくなる。建築基準法では「高さ方向の区分ごとに異なるαを規定」しているのはそのためだ。ツールの結果を設計に直接使う前には、適用する規格や指針の定義を必ず確認してね。

最後に、「このプロファイルは“中性”という特殊な状態」だと覚えておいて。現実の大気は、日中の太陽加熱で不安定になったり、夜間に安定層ができたりする。例えば、晴天の昼間は熱対流で低高度の風速分布が対数則からずれることがある。このツールはあくまで基礎となる「基準状態」を計算するもの。実務では、季節や時刻による大気安定度の影響を別途考慮する必要があるんだ。

関連する工学分野

この風速プロファイルの計算は、思っているよりずっと幅広い分野の根底にある基礎技術なんだ。例えば、「都市気候学」や「環境工学」では、街中の熱環境や汚染物質の拡散を予測するために必須だ。排ガスがどの高さでどう広がるかをCFDでシミュレーションするとき、計算領域の入口にこの対数則プロファイルを設定する。そうしないと、現実とはかけ離れた空気の流れを計算してしまうことになる。

また、「航空工学」、特に離着陸の分野とも深く関わる。飛行機が滑走路に近づくとき、この大気境界層の中を通過する。パイロットは「ウィンドシア」と呼ばれる風速・風向の急激な変化に注意するけど、その原因の一つが地表粗度の変化による境界層プロファイルの乱れなんだ。空港周辺の地形や建物の影響を評価するのにも同じ理論が使われているよ。

さらには、「サウンドエンジニアリング」(音響工学)の世界でも応用されている。風車や建築物から発生する騒音の伝播は、風速と乱流の鉛直分布に大きく左右される。音が風下にどう運ばれ、減衰するかを正確に予測するためには、ここで計算している風速と乱流強度のプロファイルが重要な入力データになるんだ。

発展的な学習のために

このツールの背後にある理論をもっと深く理解したいなら、まずは「レイノルズ平均」の概念から入るのがおすすめだ。我々が扱っている平均風速 $U(z)$ は、瞬間的な風速の激しい揺らぎ(乱流)を時間平均したもの。この「平均化」の考え方が、対数則を導く出発点になる。教科書だと「乱流境界層」の章を読んでみよう。

次にステップアップとして、「大気安定度」の影響を学ぶと視界が一気に広がる。ここで扱った対数則は「中性」という条件付き。現実には、$$ \frac{U(z)}{u_*} = \frac{1}{\kappa} \left[ \ln\left(\frac{z}{z_0}\right) - \Psi_M\left(\frac{z}{L}\right) \right] $$ のような補正関数 $\Psi_M$ を導入する($L$はモニン・オブコフ長)。この補正を理解すれば、日射や地表面温度が風のプロファイルをどう変えるかがわかるようになる。

実践的な次の一歩としては、「CFDソフトウェアで実際に境界条件として設定してみる」ことを強く勧める。無料のオープンソースソルバーでも、流入境界の設定項目に「Log Law Profile」や「Power Law Profile」の入力欄があるはずだ。ツールで計算したパラメータ($z_0$ や $u_*$)を入力し、単純な流路内の流れをシミュレーションしてみよう。数式が実際の「流れ場」として可視化される体験は、理解を飛躍的に深めてくれるよ。