粘性散逸 — 先端技術と研究動向

粘弾性流体では、エネルギー散逸が「粘性的部分」と「弾性的部分」に分かれる。弾性エネルギーは一時的に蓄えられ、後に散逸するか仕事として回収される。

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}$）からずれる原因の一つとして、粘性散逸による粘度変化が指摘されている。

• 非ニュートン流体の粘性散逸最適化: 射出成形でのゲート設計にadjoint法を適用し、散逸による温度上昇を最小化
• 乱流散逸の直接計測: PIV/PLIFの組み合わせで乱流散逸率を実験的に高精度測定し、LES/DNSの検証データを構築
• エントロピー生成最小化: 粘性散逸をエントロピー生成率として定量化し、熱交換器の最適設計に適用（Bejan's entropy generation minimization）
• 量子流体力学: 超流動ヘリウムでの散逸メカニズム（量子渦のリコネクションによるエネルギーカスケード）