🙋
「レイリー流れ」って初めて聞きました。普通の管路の流れと何が違うんですか？
🎓
いい質問だね。レイリー流れというのは、一定断面の管にガスが流れていて、その途中で「熱だけ」を加えたり奪ったりする、という理想化された圧縮性流れのモデルなんだ。摩擦は無視する。これがジェットエンジンの燃焼器とか、ラムジェットの加熱ダクトの「教科書モデル」になっている。実機の燃焼器は強烈な燃焼熱が一気に流れに入るから、摩擦よりも熱の出入りが圧倒的に支配的で、レイリー流れで十分に直感がつかめるんだよ。
🙋
なるほど。でも、ガスに熱を加えたら「温度が上がる」だけじゃないんですか？シミュレーターで加熱量 q を上げると、出口マッハ数 M₂ が上がっていって驚いたんですが…
🎓
そこがレイリー流れ最大のキモなんだ。亜音速の流れに熱を加えると、ガスは温まって膨張する。管の断面積は一定だから、質量保存則 ρuA=一定 から、密度ρが下がった分だけ速度 u が上がる。マッハ数 M=u/√(γRT) は、温度上昇で音速も少し上がるけど、亜音速では「速度上昇」が「音速上昇」を上回って、M が 1 に向かってグングン上がっていく。「加熱で流れが加速する」というのは、初学者の直感と真っ向から衝突する有名な結果だよ。
🙋
じゃあ熱をどんどん加えていけば、ずっと加速し続けるんですか？
🎓
いや、そこに「熱閉塞（thermal choking）」という硬い天井があるんだ。ある熱量を加えた時点で、出口マッハ数がちょうど 1 に到達する。それ以上は同じ管に熱を「物理的に」入れられない。本ツールでは q_max として表示している。デフォルト条件（M₁=0.3、空気）だと q_max≈578 kJ/kg。これを超えるとどうなるかというと、上流の圧力が勝手に上がって流量が落ち、運転点が変わってしまう。実機の燃焼器では、設計マージンを取って q が q_max の 60〜70% を超えないように燃料投入を制限するのが普通だね。
🙋
出口マッハ数 M₂ って、解析的に解けるんですか？式を見ると T₀ について整理されているような…
🎓
いいところに気がついた。レイリー関数 T₀/T₀* = (γ+1)M²(2+(γ−1)M²)/(1+γM²)² は、M の 4 次式が分母分子に出てくるので、M について閉じた形では解けないんだ。だから本ツールでは、「加熱後の T₀₂ から決まる比 T₀₂/T₀*」を計算して、その値になる M を [M₁, 1] の区間で二分法で 60 回反復して解いている。亜音速領域ではレイリー関数は単調増加なので、二分法でしっかり収束する。古典的な気体力学の教科書では、これを表（Rayleigh table）として暗記しろと言われていたけど、数値計算で一瞬で解けるようになったわけだね。
🙋
勉強になります！最後に、これって実機のどこで使われるんでしょうか？
🎓
代表は航空機エンジンの燃焼器、ラムジェット、スクラムジェット、アフターバーナー、それから固体ロケットの後部ノズル前のダクトだね。どれも「一定断面に熱を一気に入れる」という設定が支配的。地上の工業炉でも、長いダクトに加熱コイルを入れて流量を制御するような場面では、レイリー流れの考え方で「これ以上加熱したら流量が決まる側に行く＝熱閉塞」という限界をすぐに見積もれる。圧縮性流体力学の中で最も実機に近い理想モデルの一つだよ。
レイリー流れとは何ですか？フェノ流れとどう違いますか？
レイリー流れは、一定断面の管路を流れる理想気体に「熱だけ」を出し入れし、摩擦は無視するという理想モデルです。加熱が支配的でジェットエンジンの燃焼器のような加熱を扱うのに最適です。一方フェノ流れは、断熱（熱の出入りなし）で「摩擦だけ」が働く流れを扱う双子のモデルで、長い管路の流れ損失を見るのに使います。両者とも 1次元・等断面の圧縮性流れですが、駆動するエネルギーが熱か摩擦かで使い分けます。
なぜ亜音速の流れに熱を加えると、流れが速くなるのですか？
熱を加えると気体は膨張して密度ρが下がります。一方、管断面積 A は一定で、質量保存則から ρ·u·A = 一定 なので、密度が下がった分だけ速度 u は上がります。マッハ数は M = u/√(γRT) で、加熱で T も上がりますが、亜音速領域では「密度の低下による速度上昇」が「音速の上昇」を上回るため、結果としてマッハ数は 1 に向かって上がっていきます。これがレイリー流れ最大の直観に反する結果です。
熱閉塞（thermal choking）とは何ですか？避ける方法は？
亜音速の流れに熱を加え続けると、出口マッハ数がぴったり 1 に到達した時点で、その管にはそれ以上熱を入れられなくなります。これが熱閉塞です。本ツールでは q_max（最大加熱量）として表示しています。回避策は (1) 上流の質量流量を下げる、(2) 管断面を広げる、(3) 上流圧力が上がって運転点が変わるのを受け入れる、の3つです。実機では燃料投入量を制限したり、燃焼器を広げて流速を落とす設計で対処します。
出口マッハ数 M2 はどのように求めていますか？
レイリー関数 T0/T0* = (γ+1)M²(2+(γ−1)M²)/(1+γM²)² は M=1 で 1 になる単調関数（亜音速領域）です。加熱後の全温 T02 を求め、その T02/T0* に等しくなる M を、区間 [M1, 1] の二分法で 60 回反復して解いています。これにより閉じた解析解のないこの式を数値的にミリ秒で求められます。加熱量が q_max を超える場合は target>1 となるため、M2=1 にクランプし「ほぼ熱閉塞」と表示します。
ジェットエンジン燃焼器： ターボジェット・ターボファンのコンバスタは、まさにレイリー流れがそのまま当てはまる場所です。圧縮機で加圧した空気に燃料を噴射・燃焼させ、流れに数百〜千 kJ/kg 程度の熱を加えます。設計時は q が q_max を超えないように燃焼器を太く（流速 M₁ を 0.2〜0.3 程度に抑えて）作り、熱閉塞による圧力損失と燃焼室出口の温度暴走を防ぎます。本ツールで M₁ を 0.5 以上に上げ、q を増やしてみると、q_max がどれほど厳しくなるか体感できます。
ラムジェット・スクラムジェット： 超音速で吸気し、燃焼器で熱を加えて推進力を得るエンジンでは、燃焼器内の流速管理が生死を分けます。ラムジェットは燃焼器を亜音速にして（インテイクで減速）レイリー流れに従わせ、スクラムジェットは超音速のまま燃焼させて熱閉塞を回避します。後者は加熱で「逆に減速して M が 1 に向かう」という超音速側の挙動を利用しています。
アフターバーナー（再加熱燃焼）： 戦闘機エンジンが瞬時に推力を 50% 増やすあの「青い炎」は、タービン後の排気にもう一度燃料を入れて熱を加える装置です。ここも完全にレイリー流れの世界。アフターバーナー使用時は熱閉塞ギリギリで運転されるため、燃料制御は q が q_max の 80〜90% を上下するシビアなフィードバック制御になっています。
気体力学・CAE 解析の事前検討： 本ツールのような 1 次元のレイリーモデルは、実機の燃焼器の詳細 CFD 解析（化学反応・乱流・3D 形状）の前の「当たりづけ」に使えます。CFD で出てきた出口マッハ数や圧力比が、レイリー解析と桁違いなら、境界条件や燃焼モデルの設定ミスを疑うサニティチェックになります。学部・修士の圧縮性流体力学の演習でも、レイリー線（h-s 線図上の）と合わせて必ず登場します。
最大の誤解は、「加熱したらガスは温度だけ上がる」 と思い込むことです。閉じた容器の中なら確かに温度だけ上がりますが、開いた管路の流れでは熱を入れると膨張→密度低下→速度上昇となり、亜音速ならマッハ数が 1 に向かって上がっていきます。逆に超音速のレイリー流れでは、加熱で減速して M が 1 に向かいます。「加熱はマッハ数を 1 に近づける」というのが、亜音速・超音速の両方をまとめる正しい言い方です。q を上げると M₂ が上がるのは亜音速のとき特有の現象です。
次に、「q_max を上回る加熱量を入れたらどうなる？」 という疑問。物理的にはそんな状態は安定して存在せず、上流の圧力が勝手に上がる（背圧として伝わる）か、流量が下がるか、衝撃波が立つかして、運転点が q ≤ q_max を満たすところに「自動的に」収束します。本ツールでは q > q_max のとき M₂ = 1 にクランプして「ほぼ熱閉塞」と表示しますが、実機ではポンピング（流量振動）や燃焼振動として現れる、極めて危険な状態です。設計時は q を q_max の 60〜70% 以下に抑えるのが鉄則です。
最後に、「レイリー流れは摩擦を無視している」 という点を忘れないことです。実際の長い燃焼器や加熱ダクトでは、摩擦と加熱が同時に効くため、出口マッハ数はレイリー解より少し高め（摩擦も M を 1 に向ける効果がある）に出ます。短い高熱負荷の燃焼器ならレイリーでほぼ十分ですが、長くて熱負荷が低めのダクトではフェノ流れ（断熱摩擦流）との重ね合わせや一般化された Fanno+Rayleigh 解析が必要です。本ツールはあくまで「加熱が支配的な短い区間」の理想モデルとして使ってください。
