粘性散逸
理論と物理
粘性散逸とは
先生、粘性散逸って「粘性で熱が出る」くらいのイメージですけど、正確にはどういう現象ですか？
流体の運動エネルギーが粘性力の作用で不可逆的に内部エネルギー（熱）に変換される現象だ。分子レベルで言えば、流体層間のせん断に伴う分子の摩擦で運動エネルギーが散逸する。
日常的な例を挙げると:
- タイヤのゴムが走行中に発熱する（粘弾性体の変形散逸）
- スペースシャトルの大気圏再突入時の高温（衝撃波内の粘性散逸）
- 高分子の射出成形でゲート部が局所的に加熱される
- ダムの放水路で水温がわずかに上昇する
ダムの落水でも温度が上がるんですか！
理論的にはそうだ。高さ $h = 100\,\text{m}$ から落下する水のポテンシャルエネルギーがすべて熱に変わったとすると:
これはJouleが19世紀に実測した値とほぼ一致する。
粘性散逸関数の導出
エネルギー方程式における粘性散逸項を正確に書くと:
非圧縮性ニュートン流体の場合、これは:
テンソル表記では:
ここで重要なのは、$\Phi \geq 0$ が常に成り立つことだ。これは熱力学第二法則の要請で、粘性散逸は常に運動エネルギーを熱に変換する一方通行の過程だ。
「不可逆過程」ということですね。
Brinkman数 — 粘性散逸の重要度指標
粘性散逸を無視できるかどうかを判断する無次元数がBrinkman数だ。
$\mu$ は粘度、$U$ は代表速度、$k$ は熱伝導率、$\Delta T$ は代表温度差。
| Br の値 | 解釈 | 例 |
|---|---|---|
| $\text{Br} \ll 1$ | 粘性散逸は無視可能 | 通常の水の管内流 |
| $\text{Br} \sim O(1)$ | 粘性散逸を考慮すべき | 高分子加工、潤滑油膜 |
| $\text{Br} \gg 1$ | 粘性散逸が支配的 | 超高速流、宇宙再突入体 |
Br数が大きい場合にだけエネルギー方程式に散逸項を入れればいいんですね。
そうだ。不必要に散逸項を入れると計算コストが増え、収束性も悪化するので、Br数による事前判断が重要だよ。
粘性散逸の発見史——ジュール加熱から流体摩擦熱へ（ストークスの貢献）
流体の粘性散逸（流れのエネルギーが熱に変換される現象）を流体力学的に定式化したのはジョージ・ガブリエル・ストークス（G.G. Stokes）だ。1845年の論文「On the Theories of the Internal Friction of Fluids in Motion」で粘性応力テンソルとエネルギー散逸の関係を数学的に記述し、N-S方程式の粘性項の物理的意味を明確にした。これ以前はジュール（1843年）が固体抵抗での電気-熱変換（ジュール加熱）を示しており、流体の摩擦熱はその流体力学版と位置づけられる。高速飛翔体（弾道ミサイル、宇宙船再突入）での空力加熱はこの粘性散逸が起源であり、ストークスの150年前の理論がアポロカプセルの熱シールド設計にまで直結している。
各項の物理的意味
- 時間項 $\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で粘性散逸を計算するとき、どの方程式をどう変更するんですか？
エネルギー方程式にソース項として粘性散逸関数 $\Phi$ を追加する。
有限体積法での離散化では、$\Phi$ はセル体積にわたる積分として評価される:
各セルの中心で速度勾配からせん断速度を計算し、$\Phi = \mu |\dot{\gamma}|^2$ を評価してエネルギー方程式のソース項に加える。
運動方程式は変更しなくていいんですか？
基本的にはそうだ。ただし、粘性散逸による温度上昇が粘度に影響する場合（温度依存粘度）は、運動方程式にもフィードバックが生じる。この場合、エネルギー方程式と運動方程式を連成して解く必要がある。
乱流場における粘性散逸
乱流場では、粘性散逸の扱いが層流とは異なる。RANS（Reynolds平均）の枠組みでは:
$\overline{\Phi}_{\text{mean}}$ は平均速度場による散逸、$\varepsilon$ は乱流エネルギーの散逸率（$k$-$\varepsilon$ モデルの $\varepsilon$ そのもの）だ。
実は $\varepsilon$ はほとんどの乱流問題でエネルギー方程式のソース項として自動的に考慮されている。乱流運動エネルギー方程式の散逸項がまさにこれだ:
$P_k$ は乱流エネルギーの生成項、$\varepsilon$ が散逸率。この $\varepsilon$ が最終的に熱になる。
$k$-$\varepsilon$ モデルの $\varepsilon$ が粘性散逸と直結しているんですね！
まさにその通り。乱流の全運動エネルギーのカスケードの最終段階が粘性散逸（$\varepsilon$）だ。Kolmogorovのスケーリングでは、散逸が起きる最小スケールは:
数値精度の注意点
粘性散逸の数値計算で注意すべき点をまとめよう。
| 項目 | 注意点 | 対策 | ||
|---|---|---|---|---|
| 速度勾配の精度 | $\Phi$ は速度勾配の2乗に比例するため、勾配精度がダイレクトに効く | 二次精度以上のスキーム使用 | ||
| メッシュ解像度 | 壁面近傍のせん断が最大となる領域 | 壁面メッシュを十分に細かく | ||
| 数値散逸との区別 | 風上差分の数値粘性も非物理的な散逸を生む | 高次精度スキームで数値散逸を低減 | ||
| 非ニュートン流体 | $\Phi = \eta(\dot{\gamma}) | \dot{\gamma} | ^2$ でモデル依存 | 粘度モデルと整合したΦの実装を確認 |
数値散逸と物理的な粘性散逸を混同しないように気をつける必要があるんですね。
特にLES（Large Eddy Simulation）ではSGS（Sub-Grid Scale）モデルが追加の散逸を導入する。物理的な散逸とモデル散逸の切り分けが精度に直結するんだ。
粘性散逸の数値的扱い——高速流体でのエネルギー保存精度確保
高Mach数流れや高粘性流体の解析では、速度勾配から熱が生じる「粘性散逸（Viscous Dissipation）項Φ=τ:∇u」をエネルギー方程式で正確に解くことが重要だ。この項は速度の2乗に比例するため、流速が高い領域で急激に大きくなる。数値的に問題になるのはメッシュが粗い領域で速度勾配が過小評価されると散逸も過少評価され、エネルギー保存が崩れること。実務的には「Brinkman数Br=μU²/(kΔT)」で散逸の重要性を事前評価し、Br>0.1なら散逸を必ず有効にする。また散逸項を含む計算では温度残差と速度残差が相互に影響するため、収束基準（残差1e-6以下）をより厳格に設定する必要がある。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
粘性散逸を考慮すべきケースの判断
実務で粘性散逸を有効にすべきかどうか、どう判断すればいいですか？
まずBrinkman数を概算する。それに加えて、以下のケースでは粘性散逸を考慮すべきだ。
| 応用分野 | 典型的なBr | 散逸の影響 |
|---|---|---|
| 高分子射出成形 | $1 \sim 100$ | ゲート部で50-100K以上の温度上昇 |
| 軸受・潤滑膜 | $0.1 \sim 10$ | 油膜温度が性能を支配 |
| 押出成形（ダイ内） | $1 \sim 50$ | スクリュー回転による発熱 |
| マイクロチャネル流れ | $0.01 \sim 1$ | チャネル径が小さいとBr増大 |
| 超音速/極超音速流 | $\gg 1$ | 衝撃波内の断熱温度上昇 |
| 一般的な水配管 | $\ll 0.01$ | 無視可能 |
高分子加工ではほぼ必須なんですね。
そうだ。射出成形ではゲート通過時にせん断速度が $10^4 \sim 10^5\,\text{s}^{-1}$ に達し、粘度が $10^2 \sim 10^3\,\text{Pa}\cdot\text{s}$ あるので、粘性散逸による温度上昇は数十度になる。この温度上昇を無視すると、充填パターンの予測が全く合わない。
Ansys Fluent
Fluentでの粘性散逸の有効化手順:
1. Models > Energy をONにする
2. Energy Dialog > Options で「Viscous Dissipation」にチェック
3. これだけでエネルギー方程式に $\Phi$ のソース項が追加される
TUIコマンドでは: /define/models/energy yes yes yes （最後のyesがviscous dissipation）
OpenFOAM
OpenFOAMでは使用するソルバーによって扱いが異なる。
- buoyantPimpleFoam: エネルギー方程式に粘性散逸が含まれている
- chtMultiRegionFoam: 流体領域のエネルギー方程式でサポート
- カスタムソルバーの場合は
fvm::Sp(mu*magSqr(symm(fvc::grad(U))), he)のような項を追加
STAR-CCM+
STAR-CCM+では:
1. Physics > Energy Model を有効化
2. Physics > Viscous Dissipation Model を有効化
3. 自動的にエネルギー方程式に組み込まれる
検証問題: Couette流れの粘性散逸
解析解がある問題で検証しよう。上壁が速度 $U$ で移動するCouette流れでの定常温度分布は:
上下壁とも温度 $T_0$ に固定した場合、最高温度は中央（$y = H/2$）で:
解析解と比較すれば、ソフトの粘性散逸の実装が正しいことを確認できますね。
まさにその通り。2D平行板でメッシュを十分に細かく（$H$ 方向に50セル以上）して計算し、解析解との一致を確認してから実問題に進もう。水（$\mu = 10^{-3}\,\text{Pa}\cdot\text{s}$, $k = 0.6\,\text{W/(m}\cdot\text{K)}$）で $U = 10\,\text{m/s}$, $H = 1\,\text{mm}$ だと $\Delta T \approx 0.02\,\text{K}$ で微小。高分子メルト（$\mu = 1000\,\text{Pa}\cdot\text{s}$, $k = 0.2$）なら $\Delta T \approx 62,500\,\text{K}$ と計算上は巨大になるが、実際にはせん断発熱で粘度が下がるフィードバックが働く。
粘性散逸でコーヒーが温まる——攪拌の意外な物理
コーヒーをスプーンでかき混ぜると、わずかに温度が上がります。これが粘性散逸の実体です。速度勾配が大きいほど散逸が大きく、熱に変わる。実際には回転数1000rpmの工業用ミキサーで粘性の高いペースト（η≈100 Pa·s）を攪拌すると、1分間で数℃上昇することもあります。食品・製薬業界のプロセス設計では「散逸熱による品質劣化」が深刻な問題で、攪拌翼形状の最適化は粘性散逸分布の制御そのものです。地味に見えてこんな場面で主役になる項です。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
粘性散逸対応ツールの比較
粘性散逸の計算に強いソフトウェアはどれですか？
粘性散逸自体はほぼすべてのCFDソルバーでサポートされているが、非ニュートン流体との連成や温度依存粘度の扱いで差が出る。
| ツール | 粘性散逸 | 非ニュートン連成 | 温度依存粘度 | 特記事項 |
|---|---|---|---|---|
| Ansys Fluent | チェック一つでON | 完全対応 | Arrhenius, WLF他 | UDFで完全カスタマイズ可能 |
| STAR-CCM+ | モデル選択でON | 完全対応 | 各種モデル内蔵 | Field Functionで柔軟な定義 |
| OpenFOAM | ソルバー依存 | 対応（要設定） | UDF相当で実装 | ソースコード修正で完全自由 |
| COMSOL | Weakフォームで追加 | 完全対応 | 任意式入力可 | マルチフィジックス連成が容易 |
| Ansys Polyflow | 標準搭載 | 粘弾性含む | 各種モデル | 高分子加工専用で最も充実 |
| Moldflow | 自動考慮 | 射出成形特化 | 金型DB連携 | 射出成形ではデフォルトON |
高分子加工分野での選定
粘性散逸が特に重要な高分子加工分野では、専用ソフトの利点が大きい。
Ansys Polyflow は粘弾性流体のdieスウェル（ダイ出口での膨張）や粘性散逸による温度分布を正確に計算できる。Free surface tracking、ブロー成形、熱成形にも対応。
Moldflow は射出成形に特化しており、粘性散逸はデフォルトでエネルギー方程式に含まれる。Cross-WLFモデルが標準で、材料データベースが充実している:
Moldflowの材料データベースは具体的にどのくらいの材料が入っているんですか？
約10,000種以上のプラスチック材料がプリセットされている。各材料のCross-WLFパラメータ、PVT特性、熱物性がすべて含まれているので、材料グレードを選ぶだけで粘性散逸を含む正確なシミュレーションが可能だ。
汎用CFDでの選定
潤滑問題や一般的な熱流体問題での粘性散逸なら、FluentやSTAR-CCM+で十分だ。判断基準をまとめよう。
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 射出成形 | Moldflow, Moldex3D | 材料DB、金型冷却連成 |
| 押出・ブロー成形 | Polyflow | 粘弾性、自由表面 |
| 潤滑・軸受 | Fluent, STAR-CCM+ | 薄膜流れの取り扱い |
| マイクロチャネル | COMSOL, Fluent | マルチフィジックス連成 |
| 極超音速流 | Fluent (Density-based) | 圧縮性流れ＋散逸 |
用途によって使い分けが明確なんですね。
粘性散逸モデルのツール実装差——FluentとCFX、デフォルト設定の落とし穴
ANSYSのFluent と CFXは同じ会社の製品でも粘性散逸の扱いに差がある。Fluentでは「Energy Equation」パネルの「Viscous Heating」オプションが非圧縮流体ではデフォルトOFFになっており、高速・高粘性流れで意図せず散逸を無視する設定ミスが起きやすい。CFXでは「Heat Transfer」設定でThermal Energy（散逸なし）とTotal Energy（散逸あり）を明示的に選ぶ仕組みになっており、選択ミスが少ない設計だ。OpenFOAMのbuoyantSimpleFoamは粘性散逸をデフォルトで含む実装だが、viscousフラグをfalseにするユーザーが多く、設定の意図確認が必要。どのツールでも設定後に単純なCouetteフロー（解析解既知）でエネルギーバランスを検証することが確実な方法だ。
選定で最も重要な3つの問い
- 「何を解くか」：粘性散逸に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
粘弾性流体における散逸
粘弾性流体の場合、粘性散逸はどう変わるんですか？
粘弾性流体では、エネルギー散逸が「粘性的部分」と「弾性的部分」に分かれる。弾性エネルギーは一時的に蓄えられ、後に散逸するか仕事として回収される。
Oldroyd-Bモデルでの散逸関数は:
Peters & Baaijens (1997) は粘弾性流体における散逸の正確な定式化を示し、弾性蓄積と不可逆散逸を分離する方法を提案した。射出成形のゲート部などで、弾性効果が温度場に与える影響が数十Kに達することがある。
弾性エネルギーの蓄積と散逸を区別するのが粘弾性流体のポイントなんですね。
極超音速流における粘性散逸
Ma > 5の極超音速流では、衝撃波を通過する際の温度上昇が数千Kに達する。
正常衝撃波での温度比は:
Ma = 10 の場合（空気）、$T_2/T_1 \approx 20$。入口 $T_1 = 250\,\text{K}$ なら衝撃波後に約5000 Kとなり、化学反応（分子の解離、電離）が始まる。
このレベルでは通常のNavier-Stokes方程式の枠組みでは不十分で:
- 高温気体効果: $c_p(T)$, $\mu(T)$, $k(T)$ の温度依存性
- 化学反応: $\text{N}_2 \rightleftharpoons 2\text{N}$, $\text{O}_2 \rightleftharpoons 2\text{O}$ の解離反応
- 放射伝熱: 高温ガスからの輻射
これらをすべて連成して解く必要がある。NASAのUS3DやDplr、Ansys Fluent(Density-based solver with chemical kinetics)が使われる。
大気圏再突入って、粘性散逸の極端なケースなんですね。
マイクロスケールでの粘性散逸
マイクロチャネル（$D_h < 100\,\mu\text{m}$）では、マクロスケールでは無視できた粘性散逸が重要になる場合がある。
Brinkman数は $\text{Br} \propto U^2 / \Delta T$ だが、マイクロチャネルでは流速が大きく、かつ外部からの加熱/冷却による $\Delta T$ が小さいため、Brが無視できなくなる。
実験的にも、マイクロチャネル内の摩擦係数がマクロスケールの理論値（$f = 64/\text{Re}$）からずれる原因の一つとして、粘性散逸による粘度変化が指摘されている。
研究動向
粘性散逸の研究の最前線はどのあたりですか？
古典的な概念なのに、まだまだ研究の余地があるんですね。
超音速エンジン内の粘性散逸——発熱量は想像を超える
スクラムジェットエンジン（極超音速燃焼ラムジェット）では、流入空気の速度がマッハ5〜10に達します。このとき粘性散逸による加熱は非常に大きく、燃焼室入口ですでに空気温度が1000K以上に跳ね上がる。粘性散逸項 Φ=μ(∂u_i/∂x_j+∂u_j/∂x_i)∂u_i/∂x_j を無視した解析では、入口温度を100〜200Kも低く見積もってしまい、燃料の着火タイミング設計が根本から狂う。先端の超音速流体解析では、粘性散逸は「細かい補正」ではなく、主要な物理です。
トラブルシューティング
粘性散逸計算の典型的トラブル
粘性散逸を有効にしたら計算がうまくいかなくなった、というケースを教えてください。
粘性散逸特有のトラブルパターンを整理しよう。
1. 温度が非物理的に上昇して発散
症状: 局所的に温度が数万Kになり発散する
原因: 粘性散逸 → 温度上昇 → 粘度低下（温度依存粘度の場合）→ せん断速度増大 → さらに散逸増大、という正のフィードバックループ
対策:
- 粘度の下限値を設定する（物理的に妥当な範囲で）
- 時間ステップを小さくする（非定常計算の場合）
- 緩和係数を下げる（エネルギー方程式: 0.7-0.8、粘度: 0.5-0.7）
- まず等温（粘度固定）で流れ場を収束させてからエネルギー方程式をONにする
2. 壁面近傍で温度スパイクが発生
症状: 壁面の一層目のセルで温度が異常に高い
原因: メッシュが粗く、壁面のせん断速度が過大評価されている。$\Phi \propto \mu (\partial u/\partial y)^2$ なので、壁面第一層で速度勾配を過大評価すると散逸が急増する。
対策:
- 壁面近傍のメッシュを十分に細かくする
- プリズム層の成長率を1.2以下に抑える
- $y^+$ を確認し、壁関数モデルと整合する値になっているか確認
メッシュ品質が散逸の精度を直接左右するんですね。
3. 散逸量が理論値より過大になる
症状: Couette流やPoiseuille流で、散逸による温度上昇が解析解の数倍になる
原因: 数値粘性（人工散逸）が物理的散逸に上乗せされている。特に一次風上差分（First Order Upwind）で顕著。
対策:
- 二次精度以上の離散化スキームを使用（Second Order Upwind, QUICK, Central Differencing）
- メッシュを十分に細分化してRichardson外挿で格子収束を確認
- OpenFOAMの場合、
fvSchemesでGauss linearやGauss limitedLinear 1を使用
4. 粘性散逸をONにしたら収束が大幅に悪化
症状: 散逸なしでは300反復で収束していたのに、ONにすると5000反復でも収束しない
原因: エネルギー方程式と運動方程式の連成が強い場合（高Br数）、セグリゲート法での交互求解が収束しにくい
対策:
- 連成ソルバー（Coupled Solver）を使用
- エネルギー方程式の緩和係数を段階的に上げる（最初は0.5、安定したら0.8-0.9）
- 擬似時間ステップ（Pseudo Transient）をONにする
高Br数ほど連成が強くなるから収束が難しいんですね。
チェックリスト
粘性散逸計算の品質保証チェック項目:
- Brinkman数を事前に概算し、散逸の影響度を把握したか
- 解析解のある問題（Couette流、Poiseuille流）で検証したか
- メッシュ収束を確認し、散逸量がメッシュ非依存になっているか
- 数値スキームの次数は二次以上か
- 温度依存粘度を使う場合、粘度の上下限は適切に設定されているか
- 断熱壁 vs. 等温壁の境界条件は問題に適合しているか
Br数の事前推定から始めて、検証問題で確認してから本問題に進む、という流れですね。
その通り。粘性散逸は物理的に正しい現象だが、数値的な人工散逸と混同しないよう、系統的な検証が不可欠だ。
CFDで温度が際限なく上昇し続ける——粘性散逸の二重計上とソルバの収束問題
圧縮性流れのCFDで「温度が反復ごとに際限なく上昇し続け収束しない」症状は、エネルギー方程式の粘性散逸項が二重計上されているサインの場合がある。Fluentでは「Viscous Heating」をONにしつつ「Compressibility Effect」も有効にすると、一部のソルバ設定で散逸の貢献が重複して加算される既知の実装上の問題がある。また低速非圧縮流体でViscous Heating ONにすると、理論的には無視すべき微小な散逸が蓄積してエネルギーバランスが崩れることがある。切り分け方は①エネルギー方程式をOFFにして速度場が収束するか確認する、②Brinkman数Brを計算し散逸の重要性を評価する、③マニュアルの「Energy Source Terms」実装仕様を確認する、の3ステップが有効だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——粘性散逸の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告