粘性散逸
粘性散逸の理論基礎
粘性散逸とは
先生、粘性散逸って「粘性で熱が出る」くらいのイメージですけど、正確にはどういう現象ですか？
流体の運動エネルギーが粘性力の作用で不可逆的に内部エネルギー（熱）に変換される現象だ。分子レベルで言えば、流体層間のせん断に伴う分子の摩擦で運動エネルギーが散逸する。
日常的な例を挙げると:
- タイヤのゴムが走行中に発熱する（粘弾性体の変形散逸）
- スペースシャトルの大気圏再突入時の高温（衝撃波内の粘性散逸）
- 高分子の射出成形でゲート部が局所的に加熱される
- ダムの放水路で水温がわずかに上昇する
ダムの落水でも温度が上がるんですか！
理論的にはそうだ。高さ $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年前の理論がアポロカプセルの熱シールド設計にまで直結している。
粘性散逸の数値計算手法
エネルギー方程式への粘性散逸の組み込み
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以下）をより厳格に設定する必要がある。
粘性散逸の実務適用
粘性散逸を考慮すべきケースの判断
実務で粘性散逸を有効にすべきかどうか、どう判断すればいいですか？
まず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ソルバーでサポートされているが、非ニュートン流体との連成や温度依存粘度の扱いで差が出る。
| ツール | 粘性散逸 | 非ニュートン連成 | 温度依存粘度 | 特記事項 |
|---|---|---|---|---|
| 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フロー（解析解既知）でエネルギーバランスを検証することが確実な方法だ。
粘性散逸の先端研究
粘弾性流体における散逸
粘弾性流体の場合、粘性散逸はどう変わるんですか？
粘弾性流体では、エネルギー散逸が「粘性的部分」と「弾性的部分」に分かれる。弾性エネルギーは一時的に蓄えられ、後に散逸するか仕事として回収される。
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ステップが有効だ。
関連トピック
なった
詳しく
報告