大気圏層構造シミュレーター 戻る
大気科学

大気圏層構造シミュレーター

国際標準大気(ISA)モデルで高度ごとの気圧・気温・密度・音速をリアルタイム計算。対流圏・成層圏・中間圏・熱圏の各層をグラフで比較し、大気科学の基礎を学べます。

プローブ上昇シミュレーション

km/s(描画用)
プリセット高度
対流圏(0〜11 km)
気圧
hPa
ライブ計算結果
高度
— km
現在の層
気温
— °C
密度 ρ
— kg/m³
音速 a
— m/s
地表比(気圧)
— %
動粘度 ν
— m²/s
平均自由行程
— nm
大気圏を上昇するプローブ
温度プロファイル(ISA)
気圧・密度プロファイル
理論・主要公式

対流圏(0〜11 km、気温減率 \(L=-6.5\) K/km):
\(T = 288.15 - 6.5\,h\) [K]、\(\;p = 101325\left(\frac{T}{288.15}\right)^{-g/(RL)}\) [Pa]

成層圏下部の等温層(11〜20 km):\(T = 216.65\) K、\(\;p = p_{11}\,e^{-g(h-h_{11})/(RT)}\)

密度と音速:\(\rho = p/(RT)\)、\(\;a = \sqrt{\gamma RT}\)。ここで \(R=287\) J/(kg·K)、\(g=9.80665\) m/s²、\(\gamma=1.4\)。

🎓 会話で学ぶ大気圏の構造

🙋
山に登ると寒くなりますよね。高度が高いほど太陽に近いのに、なぜ寒くなるんですか?太陽から近い方が暖かいはずじゃないんですか?
🎓
面白い疑問だ。太陽と地球の距離は1億5千万kmで、山の高さは最大でも8.8km。山の高さは太陽との距離に対してほぼ無視できる程度しか変わらない。大気が暖まるのは主に「地表から放射された赤外線を大気が吸収する」から。高度が上がるほど地表から遠くなり、大気密度も下がって吸収量が減る。だから対流圏は下が暖かい。
🙋
でも成層圏では逆に高くなるほど温度が上がるって言いましたよね。それはオゾン層のせいですか?具体的にはどんなメカニズムですか?
🎓
そう。オゾン(O₃)は高エネルギーの紫外線(UV-B: 280〜315nm, UV-C: 100〜280nm)を直接吸収して熱に変える。この吸収層が約15〜35kmに集中しているので、その高度では温度が上昇する。成層圏頂部(約50km)の気温は−3℃程度まで戻る。オゾン層が加熱源になることで成層圏は「暖かい蓋」になり、対流圏の水蒸気・雲が上に抜けにくくなる。
🙋
旅客機は高度10〜12kmを飛ぶとのことですが、そこって成層圏ですよね。なんでそんな高いところを飛ぶんですか?
🎓
主に燃費のため。大気密度が低いと空気抵抗が減る。高度10kmの密度は地表の約34%(気圧は約26%)。摩擦力は密度に比例するので抵抗が大幅に減る。一方、揚力も密度に比例して減るので速度を上げる必要があるが、それでも総合的に燃費が良い。また成層圏は乱気流が少なく飛行が安定する。ジェットエンジンは低温(−50°C付近)で熱効率が向上するメリットもある。
🙋
CAEの分野でこの大気モデルはどんな使われ方をするんですか?
🎓
航空・宇宙CFDで必須だよ。航空機の空力解析では高度ごとのISA気圧・気温・密度をCFDの境界条件として入力する。エンジン性能計算、翼のレイノルズ数計算、再突入カプセルの加熱解析(高層大気の密度)などに使われる。大気中の音速(\(a=\sqrt{\gamma RT}\))が変わるとマッハ数が変わるので、遷音速・超音速飛行のシミュレーションでは特に重要だ。

よくある質問

国際標準大気(ISA)とは何ですか?
ICAO(国際民間航空機関)が定める高度と気温・気圧・密度の標準モデルです。海面での基準値は気圧101325Pa、気温288.15K(15°C)、密度1.225kg/m³。対流圏(0〜11km)では6.5K/kmで気温が低下、成層圏下部(11〜20km)は等温216.65K、成層圏上部(20〜32km)は1K/kmで上昇、など複数の層で定義されます。
高度8848m(エベレスト)の気圧はどれくらいですか?
ISAモデルによると約314hPa(海面気圧の約31%)、気温は約−42°Cです。空気の酸素分圧も海面の約31%になるため、補助酸素なしでの登頂は非常に困難です。ヒラリーとテンジンが1953年に初登頂した際も酸素ボンベを使用しました。ISAは標準年の平均値なので実際の季節・緯度による変動があります。
音速は高度によってどう変わりますか?
音速は \(a = \sqrt{\gamma RT}\) で気温だけに依存し(γ=1.4、R=287 J/kgK)、気圧や密度には直接依存しません。対流圏では高度が上がるほど気温が下がり音速も低下します。海面340m/s → 成層圏下部(T=216.65K)で約295m/s。高速機のマッハ数はこの局所音速との比なので、同じ飛行速度でも高高度ほどマッハ数が大きくなります。
中間圏と熱圏の特徴を教えてください。
中間圏(50〜80km)は再び温度が低下し、80km付近が大気最低温度(約−90°C)になります。流星が燃える層でもあります。熱圏(80km以上)は太陽の極端紫外線(EUV)を吸収して1000°C以上になりますが、空気が極めて薄いため熱量は少ない。ISSが飛行するのも熱圏(高度400km程度)です。
ISAモデルを使った航空機CFD解析の方法は?
航空機CFDでは飛行高度のISA値(気圧p、気温T、密度ρ、動粘度μ)を計算域の遠方境界条件として設定します。動粘度はサザーランドの式 \(\mu = \mu_0(T/T_0)^{3/2}(T_0+S)/(T+S)\) で計算します。マッハ数M = V/aで決まる圧縮性流れの境界条件も高度によって変わります。OpenFOAMではfreestream境界条件で設定でき、ansysFluentでは「Operating Conditions」で高度による気圧補正ができます。

大気圏層構造シミュレーターとは

大気圏層構造シミュレーターの物理モデルは、国際標準大気(ISA)に基づき、高度 \( h \) における気温 \( T \)、気圧 \( P \)、密度 \( \rho \)、音速 \( a \) を層ごとに計算します。対流圏(高度0~11 km)では気温が高度とともに直線的に減少し、その勾配は \( \Gamma = -6.5 \times 10^{-3} \, \text{K/m} \) です。気温は \( T(h) = T_0 + \Gamma h \) で与えられ、気圧は静水圧平衡と理想気体の状態方程式から \( P(h) = P_0 \left( \frac{T(h)}{T_0} \right)^{-\frac{g}{R \Gamma}} \) と導出されます。ここで \( g \) は重力加速度、\( R \) は乾燥空気の気体定数です。成層圏では気温が一定または増加し、中間圏では再び減少、熱圏では急上昇します。密度は \( \rho = \frac{P}{R T} \) から、音速は \( a = \sqrt{\gamma R T} \)(\( \gamma \) は比熱比)から求められます。これらの式により、各層の物理量をリアルタイムで計算し、グラフ上で比較可能としています。

実世界での応用

産業での実際の使用例
航空機メーカー(例:ボーイング社)では、旅客機「787ドリームライナー」のエンジン設計時に、巡航高度約10kmの成層圏における気圧・密度データを本シミュレーターで取得。燃焼効率や推力特性の最適化に活用しています。また、気象観測機器メーカーがラジオゾンデの高度補正アルゴリズム検証に使用し、製品精度向上に貢献しています。

研究・教育での活用
大学の大気科学科では、対流圏から熱圏までの気温変化をリアルタイム可視化し、学生が各層の物理特性を直感的に理解する教材として採用。研究分野では、成層圏オゾン層の季節変動解析や、中間圏の大気重力波伝搬シミュレーションの初期条件設定に利用され、学術論文のデータ検証にも役立っています。

CAE解析との連携や実務での位置付け
本ツールはCFD(数値流体力学)解析の前処理段階で、高度別の標準大気条件を迅速に提供。例えば、ロケットノズル設計では、熱圏の低密度・高音速環境を入力値として燃焼ガス流れを解析。実務では、実験データ取得が困難な高高度域の境界条件設定を効率化し、設計サイクル短縮とコスト削減に寄与しています。

よくある誤解と注意点

「高度が上がるほど気温は下がり続ける」と思いがちですが、実際は成層圏で気温が上昇に転じます。これはオゾン層が紫外線を吸収するためで、対流圏から成層圏への遷移域である対流圏界面では気温が約−56.5℃で一定になる点に注意が必要です。

「国際標準大気は現実の大気を正確に再現する」と思いがちですが、実際は緯度・季節・天候の影響を平均化した理想モデルです。特に熱圏では太陽活動による温度変動が大きく、シミュレーション結果と実際の観測値が数百℃乖離することもあるため、あくまで標準的な傾向を把握するための基準値として利用する必要があります。

「音速は高度とともに単調に減少する」と思いがちですが、実際は気温に依存するため、成層圏で気温が上昇すると音速も増加します。特に成層圏界面付近では音速が対流圏より速くなる現象が発生し、航空機のマッハ数計算に影響を与える点に注意が必要です。

使い方

  1. 高度は0〜80 kmの範囲で、スライダーと数値入力のどちらも0.5 km刻みです。
  2. 層境界は getLayerInfo と同じく、対流圏0〜11 km、成層圏11〜47 km、成層圏頂47〜51 km、中間圏51〜80 kmとして表示します。
  3. 表示単位は気圧hPa、気温°C、密度kg/m³、音速m/sです。

計算例

高度5 kmでは T=255.65 K(-17.5°C)、p=54,020 Pa(540.2 hPa)です。高度47 km(成層圏頂付近)では T=270.65 K(-2.5°C)、p≈110.9 Pa(約1.11 hPa)となり、大気は非常に希薄です。

実務での注意点