軸流圧縮機段
理論と物理
概要
先生、軸流圧縮機って何段も翼が並んでるやつですよね？ あれ1段1段でどのくらい圧力上がるんですか？
良い質問だね。亜音速の軸流圧縮機だと1段あたりの圧力比は概ね 1.15～1.4 程度。段数を重ねて全体で 10:1 以上の圧力比を達成するのがガスタービン用圧縮機の典型だよ。
そんなに少しずつなんですか。なんで一気に圧縮しないんですか？
圧縮過程では流れが減速するから、急な減速は境界層の剥離を招く。翼列1段あたりのディフュージョンファクタ（DF）を抑えることが安定運転の鍵なんだ。Liebleinの基準では DF < 0.6 が目安とされている。
支配方程式 ― Euler仕事方程式
段のエネルギー伝達を表す式って何ですか？
ターボ機械の根幹をなすEulerのターボ機械方程式だよ。
ここで $U$ は翼周速、$C_{\theta}$ は絶対流れの旋回成分。添字1が入口、2が出口だ。この式はエネルギー保存から厳密に導かれるので、圧縮機でもタービンでも成立する。
じゃあ旋回速度の変化が大きいほど仕事が増えるわけですね。
そのとおり。ただし旋回が増えすぎると翼の負荷が上がりすぎるから、速度三角形を描いて相対マッハ数や流れ角を確認するのが設計のイロハだ。反動度 $R$ を使って仕事の分配も管理する。
$R = 0.5$ の50%反動度段は動翼と静翼の速度三角形が対称になり、損失が最小化されやすい設計として広く使われる。
de Haller基準と拡散限界
翼列の減速はどこまで許容されるんですか？
de Haller数が実務でよく使われる簡便な指標だ。
$W$ は相対速度で、この比が0.72を下回ると境界層剥離のリスクが急増する。もう少し精密な評価にはLiebleinのディフュージョンファクタを用いる。
$\sigma$ はソリディティ（翼弦/翼ピッチ）。DF < 0.6 が設計限界の目安だ。
この2つの基準はCFDの前段階、1D設計の段階で使うんですか？
その通り。1Dの平均線解析（Mean-Line Analysis）で速度三角形を決めてからCFDに進むのが標準ワークフローだよ。
商用ツールと翼列設計
軸流圧縮機の翼列設計に使われるソフトって何がありますか？
代表的なものを挙げよう。
| ツール | 用途 | 特徴 |
|---|---|---|
| Ansys CFX / TurboGrid | 3D CFD＋専用翼列メッシャー | 構造格子で翼間流路を高品質にメッシュ化 |
| NUMECA FINE/Turbo | 3D CFD（AutoGrid5） | 多段ターボ機械に特化、Non-Matching Interface対応 |
| Concepts NREC (AXIAL) | 1D/2D予備設計 | Mean-Line＋Throughflow解析が一体 |
| AxSTREAM (SoftInWay) | 1D～3D統合設計 | 予備設計からCFDへのシームレス連携 |
TurboGridって聞いたことあります。あれは何が良いんですか？
ATM最適化（Automatic Topology and Meshing）で翼形状に合わせたH/J/C/L型トポロジを自動生成してくれる。翼前縁・後縁のO-gridも自動配置されるから、境界層の解像が格段に楽になるんだ。
速度三角形が描けなかった昔のエンジニアたち
軸流圧縮機の設計に不可欠な「速度三角形」——絶対速度・相対速度・翼周速の関係を表したあのベクトル図を、19世紀末のエンジニアたちは手計算と試行錯誤で解いていました。オイラー方程式そのものは1750年代に確立されていたのに、それを翼列設計に体系的に適用する手法（Mean-Line法）が整ったのは1940年代。ジェットエンジン開発の急速な需要が「速度三角形の実用化」を一気に加速させた歴史があります。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
回転系の定式化
圧縮機の動翼をCFDで解くとき、回転をどう扱うんですか？
最も一般的なのはMRF（Multiple Reference Frame）法、つまり回転座標系での定常解析だ。運動量方程式にコリオリ力と遠心力の体積力項が加わる。
$\mathbf{v}_r$ は回転座標系での相対速度、$\boldsymbol{\omega}$ は角速度ベクトルだ。
コリオリ力のところが追加分ですね。これって定常計算なら時間微分はゼロですか？
そう、定常MRFなら左辺第1項はゼロだ。動翼-静翼間のインターフェースにはMixing Plane（周方向平均）を置くことで、異なるピッチの翼列同士を定常的に接続できる。
乱流モデルの選択
圧縮機のCFDではどの乱流モデルが定番ですか？
実務上の標準は SST（Shear Stress Transport）k-omega モデルだ。翼面の逆圧力勾配による剥離をk-epsilon系より正確に捕らえられる。CFXでは「SST」を選ぶだけで壁近傍の自動切り替えが有効になる。
| 乱流モデル | 長所 | 短所 | 推奨用途 |
|---|---|---|---|
| SST k-omega | 逆圧力勾配に強い、壁近傍精度良好 | 遷移を捉えない | 設計点の定常解析 |
| Gamma-Theta遷移モデル | 層流-乱流遷移を予測可能 | 校正パラメータが必要 | 低Re翼型の性能評価 |
| SAS / DES | 大規模剥離の非定常構造を解像 | 計算コスト大 | 失速・サージ近傍の解析 |
遷移モデルってどういう場面で必要になりますか？
翼弦Reが $5 \times 10^5$ 以下の低速ファンや小型圧縮機では、翼面上の層流域がかなり長くなる。遷移位置の違いが損失に直結するから、Gamma-Theta（Langtry-Menter）モデルが有効だよ。
境界条件の設定
入口と出口にはどんな条件を設定するんですか？
典型的な圧縮機段の境界条件は以下だ。
- 入口: 全温 $T_0$、全圧 $p_0$、流れ角（旋回成分）、乱流強度（通常5%程度）
- 出口: 静圧 or 質量流量指定
- 翼面: No-slip、断熱壁が一般的
- ハブ/シュラウド: No-slip回転壁（動翼側）or 静止壁（静翼側）
- 周期面: 1ピッチの回転周期境界条件
出口を静圧指定にすると、流量が結果として出てくるわけですね。
そう。特性曲線を描くには、出口の背圧を段階的に上げてサージ寄りの運転点まで計算を繰り返す。質量流量が急減するあたりがストール限界だ。ただし定常計算では真のサージを捉えきれないから、限界付近では非定常解析が必要になる。
境界条件で圧縮機マップが変わる話
軸流圧縮機のCFDで「どの境界条件を使うか」は、実は結果に大きく影響します。出口に静圧を指定するか質量流量を指定するかで、サージ限界付近の挙動がまったく変わる。現場では「背圧指定のほうが収束しやすい」とよく言われますが、その代わり流量が結果任せになる。ターボ機械CAEの難しさは、こうした「何を固定して何を求めるか」の選択が技術者の腕の見せどころになる点です。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
翼列メッシュの基本方針
圧縮機のメッシュで一番気をつけるべきことは何ですか？
翼面境界層の解像だ。SSTモデルで壁関数を使わない「Low-Re解法」を選ぶなら、第一層セル高さは $y^+ \approx 1$ を狙う必要がある。
代表的な条件（翼弦 50mm、Re = $5 \times 10^5$）だと第一層高さは約 5～10 $\mu$m になることが多い。
そんなに薄いんですか。成長率はどのくらいにすればいいですか？
プリズム層の成長率は 1.1～1.2 が推奨だ。15～25層のプリズムレイヤーで翼面の境界層をしっかり覆う。TurboGridの場合、O-gridの厚さと分割数で直接制御できる。
メッシュ密度の目安
翼間全体のセル数としてはどのくらいが目安ですか？
1ピッチ1翼列あたりの目安を表にまとめよう。
| 解析目的 | セル数/ピッチ | 備考 |
|---|---|---|
| 予備評価（粗い） | 10～30万 | 全体傾向の把握、パラメトリック |
| 設計評価（標準） | 50～100万 | 効率・特性の定量予測 |
| 詳細評価（細かい） | 200～500万 | 二次流れ構造、チップ漏れの詳細 |
| LES/DES | 1000万～ | 非定常渦構造の解像 |
多段だとピッチごとに100万でも全体では大変な数になりますね。
そう。例えば10段で各段ロータ+ステータなら20翼列。1翼列80万でも1600万セルだ。だからMixing Planeで定常計算するのが実務で主流なんだよ。
チップクリアランスのモデル化
翼端のチップ隙間ってどうやってメッシュに入れるんですか？
TurboGridではチップ領域に専用のH-gridブロックが自動生成される。重要なのは隙間方向（径方向）に最低でも10～15セル入れること。チップ漏れ渦はセル数が不足すると数値拡散で消えてしまう。
実機でチップクリアランスってどのくらいですか？
大型ガスタービンの高圧圧縮機で翼高さの1～2%、小型機で2～3%程度。たったこれだけの隙間で段効率が1～3ポイント落ちるから、CFDでの正確な評価が非常に重要なんだ。
メッシュ収束性確認
メッシュが十分かどうかはどうやって判断しますか？
GCI（Grid Convergence Index）による3水準評価が基本だ。粗・中・細の3メッシュで段効率を比較し、Richardson外挿で格子非依存解を推定する。効率の変化が0.1ポイント以内に収まれば十分と判断することが多い。
y+ひとつで効率予測が1ポイント変わる
軸流圧縮機のCFD実務でよく起きるのが「メッシュは細かいのに効率が実験値と合わない」問題。原因の多くはy+の不一致です。壁関数モデルを使っているのにy+が1以下だと壁せん断力が過大評価され、効率が低めに出る。逆にLow-Reynolds型モデルでy+が10を超えると、境界層の解像が足りず摩擦損失が過小評価される。翼列CFDではy+の管理が1段の断熱効率を文字通り1ポイント以上左右することがあります。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
Mixing Planeの仕組み
多段圧縮機のCFDでMixing Planeが必須だと聞きましたが、具体的に何をしてるんですか？
動翼と静翼は翼枚数が異なるため、同時に定常計算するには界面で周方向の情報を平均化する必要がある。Mixing Planeでは上流側の出口流量を周方向に質量平均し、下流側の入口条件として渡す。
$\phi$ は全圧・全温・流れ角などの保存変数、$N_b$ は翼枚数だ。
平均化すると非定常の翼間干渉の情報は失われますよね？
その通り。ウェイク（後流）やポテンシャル干渉は捕えられない。それでも段効率や圧力比の予測精度は実験値と1～2ポイント以内で一致することが多いから、設計には十分実用的なんだ。
CFXでの多段設定
Ansys CFXで多段圧縮機を組む手順を教えてください。
大まかな流れはこうだ。
1. TurboGridで翼列ごとにメッシュ生成 - ロータ1、ステータ1、ロータ2… と個別に作成
2. CFX-Preでドメイン統合 - 各メッシュをインポートし、回転/静止を指定
3. Stage Interface設定 - 動翼-静翼間に「Stage（Mixing Plane）」を選択
4. Turbo Mode - CFX-Preの Turbo Mode で入口/出口/壁/周期面を自動検出
NUMECAだとどう違いますか？
FINE/TurboのAutoGrid5では、全段をまとめて1つのプロジェクト内でメッシュ生成できる。Mixing Planeの位置は「Row Interface」として自動配置される。段間にブリーディング（抽気）を入れる場合も専用GUIがあるから、多段の設定作業はFINE/Turboのほうが効率的な場面が多い。
特性曲線の取得
圧縮機マップを作るにはどうやって計算するんですか？
回転数一定で出口背圧を段階的に上げていく。各運転点で収束した結果から圧力比と断熱効率を抽出してプロットする。
| 運転点 | 出口背圧 | 質量流量 | 圧力比 | 断熱効率 |
|---|---|---|---|---|
| 1（チョーク近傍） | 低 | 最大 | 低 | 低～中 |
| 2（設計点近傍） | 中 | 設計値 | 設計値 | 最高 |
| 3（ストール近傍） | 高 | 低下 | ピーク付近 | 低下 |
収束しなくなる点がストール限界ですか？
定常計算ではそう解釈することが多いが、物理的にはまだストールしていない領域で数値的に収束しなくなることもある。正確なストールマージン評価には非定常計算（Sliding Meshや時間刻みMixing Plane）が必要だ。
NUMECAが航空エンジンに愛される理由
軸流圧縮機解析の世界でNUMECAのFINE/Turboが高い評価を受ける背景には、ベルギーのvon Karman Instituteとの長年の共同研究があります。翼列乱流モデルや遷音速流れの安定化スキームが、実際の圧縮機試験データを元に繰り返し検証されてきた。大手航空エンジンメーカーがベンチマークに使うレベルの精度がある一方で、セットアップの自動化機能も充実。「研究精度と実務効率の両立」というターボ機械CAEの永遠のテーマに正面から向き合ったツールといえます。
選定で最も重要な3つの問い
- 「何を解くか」：軸流圧縮機段に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
サージとストールの違い
サージとストールって何が違うんですか？ よく混同されてる気がして…
重要な区別だ。ストール（回転失速）は翼列の一部に局所的な剥離セルが発生し、回転方向に伝播する現象。サージは圧縮機全体の流れが周期的に逆流する系統不安定だ。
- 回転失速: 局所的・翼列内部の現象。圧力変動は翼通過周波数より低い
- サージ: 系統全体。流れが完全に逆転するディープサージでは振動・破損の危険
CFDでサージを予測するのは難しいですか？
非常に難しい。サージは圧縮機だけでなく下流のプレナム容積やダクト系の特性に依存するから、GreitzerのBパラメータモデルのような系統モデルとCFDを連成させるアプローチが研究されている。
$U$:翼周速、$a$:音速、$V_p$:プレナム容積、$A_c$:圧縮機流路面積、$L_c$:圧縮機長さ。$B > 0.7$ 程度でサージに移行しやすいとされる。
ケーシングトリートメント
サージマージンを改善する方法としてケーシングトリートメントがあると聞きました。
シュラウド側に周方向の溝や軸方向のスロットを設けて、チップ付近の低運動量流体を再エネルギー化する手法だ。CFDでは溝の形状・深さ・幅をパラメトリックに変えて最適化する。
| 手法 | 構造 | 効果 | 効率への影響 |
|---|---|---|---|
| 周方向溝 | シュラウドに環状溝 | ストールマージン+5～15% | 効率低下 0.5～1% |
| 軸方向スロット | 翼弦方向にスリット | ストールマージン+10～20% | 効率低下 1～2% |
| 自己循環ケーシング | ブリード＋リインジェクション | 大きなマージン改善 | 効率低下やや大 |
効率が少し落ちても安定性が大事ってことですね。
そう。航空エンジンでは鳥の吸い込みやエンジン加減速時にサージマージンが食われるから、トリートメントによる保険は重要なんだ。
最適化とAIの活用
翼形状の最適化にはどんな手法が使われていますか？
NUMECAにはそういう機能がありますか？
FINE/Designに最適化エンジンが組み込まれていて、AutoGrid5+FINE/Turbo+Design環境で翼形状パラメトリックスタディから多目的最適化まで一貫して行える。Ansys側ではoptiSLangとの連携が一般的だ。
ケーシングトリートメントの意外な歴史
ジェットエンジンのサージマージンを改善するケーシングトリートメント（シュラウドの溝・スロット加工）は、1970年代にNASAが偶然発見したとも言われる技術です。当初は「効率が少し落ちても、失速しにくくなるなら実用的」という経験的知見でした。そのメカニズムが翼端チップ漏れ流れの再エネルギー化である、とCFDで定量的に解明されたのは2000年代以降のこと。「現場の経験がCFDに証明してもらう」という流れが、ターボ機械の世界では今もよく起きています。
トラブルシューティング
典型的な収束不良パターン
圧縮機のCFDで計算が全然収束しないんですが、何を確認すればいいですか？
圧縮機CFDの収束不良は翼列特有の原因が多い。チェックリストを見てみよう。
1. 出口背圧が高すぎる
一番多い原因がこれだ。サージ限界を超えた背圧を設定すると物理的に定常解が存在しない。まずは背圧を下げて収束を確認し、段階的に上げていくこと。
2. Mixing Planeでの不整合
Mixing Planeが問題になることがあるんですか？
上流ロータの出口と下流ステータの入口でスパン方向のメッシュ分布が大きく異なると、Mixing Planeでの補間で数値的なスパイクが生じる。特にハブやシュラウド付近の曲率変化が急な箇所で起きやすい。
対策: 段間インターフェース前後のメッシュでスパン方向の節点分布を概ね揃える。
3. チップ隙間のメッシュ品質
チップクリアランスのメッシュが歪んでいると、高マッハ数のチップ漏れ流れで負の圧力が発生して発散する。
対策: チップ領域の最小直交品質が0.15以上、最大アスペクト比が1000以下になるようにメッシュを調整する。
非物理的な結果のチェック
収束はしたけど結果が怪しいときはどうすれば？
以下の項目を確認しよう。
| チェック項目 | 期待値 | 異常時の原因 |
|---|---|---|
| 質量流量の入口/出口差 | 0.01%以内 | 収束不十分、メッシュ漏れ |
| 翼面y+ | 0.5～3.0 | 境界層メッシュ不適切 |
| 出口旋回角分布 | スパンで滑らかに変化 | Mixing Plane不整合 |
| エントロピー生成分布 | 翼面・端壁に集中 | 数値拡散が支配的（メッシュ粗すぎ） |
特に効率が実験より異常に高いとか低いとか、ありがちですか？
よくある。メッシュが粗いと数値拡散で損失が過大評価される（効率低め）。逆に壁関数の適用が不適切で摩擦損失が過小評価されると効率が高く出る。y+の値を必ず確認して、選択したモデルの前提範囲に入っているかチェックすることが最重要だ。
並列計算のTips
多段で計算が重いとき、並列数はどう決めればいいですか？
CFXの場合、1パーティションあたり5万～10万セルが効率的だ。例えば800万セルなら80～160パーティション。ただしMixing Planeのインターフェースをまたぐ分割は通信コストが増えるから、ドメイン境界でのパーティショニングに注意が必要だ。
「収束した嘘」はベテランも引っかかる
軸流圧縮機CFDのあるある失敗談：残差も収束していて、コンター図もきれいで、でも実験値と全段効率で3ポイント以上ずれていた——。Mixing Planeのスパン方向補間で数値的なエネルギーが消えていたケースが現場では珍しくありません。「収束＝正確」ではない。ターボ機械CFDでは、全体のエネルギーバランス（入力仕事量とΔh_0の整合）を必ずチェックする習慣が、ベテランエンジニアを見分けるひとつの指標になっています。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——軸流圧縮機段の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告