非ニュートン流体
理論と物理
非ニュートン流体とは
先生、非ニュートン流体って「変な粘度の流体」くらいのイメージなんですけど、もう少しちゃんと教えてもらえますか？
ニュートン流体ではせん断応力 $\tau$ とせん断速度 $\dot{\gamma}$ が線形関係 $\tau = \mu \dot{\gamma}$ にある。この線形関係が成り立たない流体が非ニュートン流体だ。身の回りにも多い。血液、塗料、ケチャップ、シャンプー、高分子溶液、セメントスラリーなど。
非ニュートン流体は大きく以下に分類される。
| 種類 | 特徴 | 例 |
|---|---|---|
| Shear-thinning（擬塑性） | せん断速度↑ → 粘度↓ | 塗料、血液、高分子溶液 |
| Shear-thickening（ダイラタント） | せん断速度↑ → 粘度↑ | コーンスターチ水溶液 |
| Bingham塑性 | 降伏応力以下で変形しない | 歯磨き粉、チョコレート |
| 粘弾性流体 | 粘性と弾性の両方を持つ | 高分子メルト、DNA溶液 |
| チキソトロピー | 時間とともに粘度低下 | ペンキ、ヨーグルト |
ケチャップが出にくいのも非ニュートンなんですね！ 振ると粘度が下がるから出やすくなる。
構成方程式（粘度モデル）
非ニュートン流体の核心は、せん断速度に依存する粘度関数 $\eta(\dot{\gamma})$ の定式化にある。主要なモデルを紹介しよう。
1. Power-law（べき乗則）モデル
最もシンプルなモデルだ。
$K$ は整合性係数（consistency index）、$n$ はべき乗指数。$n < 1$ でshear-thinning、$n > 1$ でshear-thickening。$n = 1$ でニュートン流体に帰着する。
2. Carreauモデル
Power-lawの欠点（低・高せん断速度での非現実的な振る舞い）を改善したモデル。
$\eta_0$ はゼロせん断粘度、$\eta_\infty$ は無限せん断粘度、$\lambda$ は緩和時間。
3. Herschel-Bulkleyモデル
降伏応力を持つ流体（Bingham塑性の一般化）。
$\tau_y$ が降伏応力。$\tau < \tau_y$ では流動しない（剛体的挙動）。$n = 1$ でBinghamモデルに帰着する。
降伏応力のある流体って、数値的に扱いにくそうですね。
その通り。降伏面の位置が未知なので、正則化（regularization）手法がよく使われる。Papanastasiouモデルでは:
パラメータ $m$ を大きくすると理想的なBingham挙動に近づくが、数値的には剛性が増す。$m = 100\text{--}1000\,\text{s}$ が典型的な範囲だ。
非ニュートン流体論の礎——ビンガムとOstwald-de Waele（1906〜1929年）
非ニュートン流体の体系的研究は20世紀初頭に始まった。Eugene Bingham(1916)は泥漿・ペーストの「降伏応力（Yield Stress）」を測定し、降伏後は線形粘性に従う「ビンガム流体（Bingham Plastic）」モデルを提唱した。一方Ostwald(1925)とde Waele(1923)は独立に「べき乗則（Power Law）流体（τ=K·γ̇ⁿ）」を定式化し、剪断減粘性（n<1：血液、高分子溶液）と剪断増粘性（n>1：コーンスターチ水溶液）を統一的に記述した。これら100年前のモデルは今もCFDの材料モデルライブラリの基本として実装されており、より洗練されたCross、Casson、Herschel-Bulkleyモデルへの出発点となっている。
各項の物理的意味
- 時間項 $\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で解くとき、ニュートン流体と何が違うんですか？
根本的な違いは、粘度が速度場に依存することだ。Navier-Stokes方程式の粘性項が非線形になる。
$\mathbf{D} = \frac{1}{2}(\nabla\mathbf{u} + \nabla\mathbf{u}^T)$ はひずみ速度テンソルだ。
解法としてはPicard反復（逐次代入法）が基本になる。
1. 前のステップの速度場から $\dot{\gamma}$ を計算
2. $\eta(\dot{\gamma})$ を更新
3. 更新された粘度でN-S方程式を解いて新しい速度場を得る
4. 収束するまで1-3を繰り返す
ニュートン流体に比べて反復が1つ増えるんですね。
そうだ。この外側の反復ループが収束しにくいのが非ニュートン流体計算の難しさだ。特にshear-thinning度が強い（$n$ が小さい）場合や降伏応力がある場合は要注意。
一般化レイノルズ数
非ニュートン流体ではレイノルズ数の定義も変わる。Power-law流体の管内流ではMetzner-Reed一般化レイノルズ数が使われる。
層流-乱流遷移は $\text{Re}_{\text{MR}} \approx 2100$ で起きる（ニュートン流体とほぼ同じ臨界値）。ただし遷移後の乱流挙動はニュートン流体とは大きく異なる。
一般化Re数の式は結構複雑ですね。
有限体積法での離散化
CFDで非ニュートン流体を解く際の離散化の注意点を説明しよう。
粘度のフェース値補間が重要になる。セル中心で計算した粘度をフェースに補間する際:
- 調和平均: 粘度が急変する界面で推奨。$\eta_f = \frac{2\eta_L \eta_R}{\eta_L + \eta_R}$
- 算術平均: 粘度変化が緩やかな場合。$\eta_f = \frac{\eta_L + \eta_R}{2}$
shear-thinning流体では、壁面近傍と流路中心で粘度が数桁異なることがある。例えば高分子メルトの射出成形では $\eta$ が $10^1 \sim 10^4\,\text{Pa}\cdot\text{s}$ の範囲で変化する。この急変に対応するため、壁面近傍のメッシュを十分に細かくする必要がある。
| Power-law $n$ | 粘度比 $\eta_{\text{center}}/\eta_{\text{wall}}$ | メッシュ要求 |
|---|---|---|
| 0.8 | 約3 | 緩やかな細分化で十分 |
| 0.5 | 約30 | 壁面近傍に10層以上推奨 |
| 0.2 | 約1000 | 非常に細かいメッシュが必要 |
$n$ が小さいほどメッシュ要求が厳しくなるんですね。収束も難しくなりそう。
その通り。緩和係数を0.3-0.5に下げ、粘度の更新にも緩和を入れるのが実務的なテクニックだ。
非ニュートン流体CFDの数値安定性——降伏応力流体（ビンガム流体）の特異点処理
ビンガム流体（降伏応力τ_y>0）のCFDでは「降伏応力以下では固体、以上では液体」という物性の不連続性が数値不安定の原因になる。剪断速度γ̇→0の極限でビンガムモデルの粘度が無限大に発散するため、ニュートン法が収束しない。実用的な解決策はPapanastasiou(1987)の正則化近似：τ = (τ_y/γ̇)(1-exp(-mγ̇)) + ηγ̇ で、パラメータm（典型値10³〜10⁵）で収束性と精度のバランスを取る。mが大きすぎると数値剛性で発散し、小さすぎると降伏応力の効果を過少評価する。適切なm値はτ_y/η（粘性時間の逆数）の100倍程度を目安に、メッシュ精細化とセットで確認することが推奨されている。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
粘度モデルのパラメータ決定
非ニュートン流体のシミュレーションを始めるとき、粘度モデルのパラメータはどうやって決めるんですか？
レオメータ（粘度計）で測定した「せん断速度 vs. 粘度」のデータに、モデル式をフィッティングする。具体的な手順はこうだ。
1. 回転式レオメータでせん断速度 $\dot{\gamma} = 0.01 \sim 10^4\,\text{s}^{-1}$ の範囲で粘度を測定
2. 両対数プロットで $\eta$ vs. $\dot{\gamma}$ をプロット
3. 直線部分からPower-lawパラメータ（$K$, $n$）を決定
4. プラトー領域があればCarreauモデル（$\eta_0$, $\eta_\infty$, $\lambda$, $n$）にフィット
5. 降伏応力があればHerschel-Bulkleyモデルに移行
代表的な流体のパラメータ例:
| 流体 | モデル | パラメータ例 |
|---|---|---|
| 血液 | Carreau | $\eta_0 = 0.056\,\text{Pa}\cdot\text{s}$, $\eta_\infty = 0.00345\,\text{Pa}\cdot\text{s}$, $\lambda = 3.31\,\text{s}$, $n = 0.357$ |
| ポリエチレンメルト | Power-law | $K = 10^4\,\text{Pa}\cdot\text{s}^n$, $n = 0.4$ |
| セメントスラリー | Herschel-Bulkley | $\tau_y = 10\,\text{Pa}$, $K = 0.5\,\text{Pa}\cdot\text{s}^n$, $n = 0.8$ |
| ケチャップ | Herschel-Bulkley | $\tau_y = 15\,\text{Pa}$, $K = 8\,\text{Pa}\cdot\text{s}^n$, $n = 0.25$ |
実際の物性値を知っておくと、結果の妥当性判断に役立ちそうですね。
解析フロー
非ニュートン流体のCFD解析フローを整理しよう。
Step 1: レオロジーデータの準備
- レオメータ測定データの取得
- 適切な粘度モデルの選択とフィッティング
- 温度依存性がある場合はWLFモデルやArrheniusモデルも設定
Step 2: メッシュ設計
- 壁面近傍を十分に解像（Power-law流体の速度分布はニュートン流体よりフラット）
- せん断速度の大きい領域を特定して局所細分化
Step 3: ソルバー設定
- まずニュートン流体（$\eta_0$ の値）で定常解を得る
- その解を初期条件にして非ニュートンモデルをON
- 緩和係数を下げる（粘度の緩和も0.5-0.8に設定）
Step 4: 収束確認
- 残差だけでなく、壁面せん断応力や圧力損失のモニターが安定しているか確認
ニュートン流体で初期解を作ってから非ニュートンに切り替えるのがコツなんですね。
検証用の解析解
結果の妥当性確認に使える解析解を紹介しよう。
Power-law流体の管内層流速度分布:
$n = 1$ でHagen-Poiseuilleの放物線分布に帰着する。$n < 1$ では中心部がフラットになり、$n > 1$ では尖った分布になる。
$n = 0.5$ だと速度分布がかなり平たくなるわけですね。
そうだ。この解析解とCFDの結果を比較して、粘度モデルの実装が正しいことを確認してから、複雑な形状の解析に進むのがベストプラクティスだ。
チョコレート製造の非ニュートン流体CFD——テンパリング工程の粘度制御
食品工業でのCFD応用例として、チョコレートのテンパリング工程が挙げられる。融解チョコレートはカカオバターの結晶構造により粘度が温度と剪断速度に強く依存する典型的な非ニュートン流体（Herschel-Bulkleyモデルが近似に使われる）だ。テンパリング機の撹拌パイプ内でChocolate masaが所定の温度・剪断履歴を受けることで、最も安定なV型結晶が生成され光沢のある口溶けの良い製品になる。CFD（非ニュートンモデル）を使って撹拌器内の温度-剪断の均一性を最適化することで、不均一な結晶生成による「ブルーミング」（白斑）を低減した事例が食品工学分野の論文で報告されている。CFDが食品品質管理にも貢献していることを示す面白い応用例だ。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
主要CFDツールの非ニュートン流体対応
非ニュートン流体に対応しているCFDソフトを教えてください。
主要ツールごとの対応状況をまとめよう。
Ansys Fluent
Fluentは非ニュートン流体のサポートが充実している。
- 組み込みモデル: Power-law, Carreau, Cross, Herschel-Bulkley, Bingham
- 設定場所: Materials > Viscosity > Non-Newtonian で選択
- UDF対応:
DEFINE_PROPERTYマクロで任意の粘度モデルを実装可能 - 温度依存性: 各モデルにArrhenius型の温度依存を追加可能
- 粘弾性: Oldroyd-B, Giesekusモデル（Fluent内蔵ではないが、UDFで実装例あり）
Simcenter STAR-CCM+
STAR-CCM+も同等の機能を持つ。
- 組み込みモデル: Power-law, Carreau-Yasuda, Cross, Herschel-Bulkley, Bingham
- 特徴: Carreau-Yasudaモデルが標準搭載（$a$ パラメータによる遷移領域の調整が可能）
- Field Function: 任意のせん断速度依存粘度をJava APIまたはField Functionで定義
- 粘弾性: Viscoelasticモデルモジュールあり（Oldroyd-B, FENE-P）
OpenFOAM
OpenFOAMは transportModels ライブラリで非ニュートンモデルを提供する。
- 利用方法:
constant/transportPropertiesでtransportModelを指定 - 対応モデル:
powerLaw,CrossPowerLaw,BirdCarreau,HerschelBulkley - カスタムモデル: C++で
viscosityModelクラスを継承して実装 - ソルバー:
nonNewtonianIcoFoam（層流）、pimpleFoam+ 非ニュートンtransportModel（乱流）
COMSOL Multiphysics
COMSOLはGUIでのモデル設定が直感的だ。
- CFDモジュール: 非ニュートンモデルを物性値として直接設定
- 特徴: せん断速度依存の粘度を数式入力できる自由度が高い
- 粘弾性: Polymer Flowモジュール（Oldroyd-B, Giesekus, PTT）
- 連成解析: 熱・化学反応との連成が容易
専門ソフトウェア
汎用CFD以外に、非ニュートン流体専用のソフトってあるんですか？
射出成形や押出成形の分野では専用ソフトが主流だ。
| ソフトウェア | 開発元 | 用途 |
|---|---|---|
| Moldflow | Autodesk | 射出成形シミュレーション |
| Moldex3D | CoreTech System | 射出成形（3D解析） |
| Polyflow | Ansys | 高分子加工全般（押出、ブロー成形） |
| FLOW-3D | Flow Science | 自由表面を伴う非ニュートン流動 |
Ansys Polyflowは粘弾性流体の定番ツールだ。ViscoelasticモデルやDifferential/Integral型の構成方程式を幅広くサポートしている。ダイ設計や押出成形のシミュレーションでは、Fluentよりも適している場合が多い。
用途によって使い分けが必要なんですね。射出成形ならMoldflow、押出しならPolyflow、一般的なCFDならFluentやSTAR-CCM+という感じでしょうか。
その理解で正しい。選定のポイントは、必要な構成方程式のサポート状況とワークフローの効率だ。
非ニュートン流体CFDのツール対比——Fluent vs Polyflow vs OpenFOAMの高分子向け機能
高分子溶融体や粘弾性流体（ポリマー押出、射出成形）のCFDには専用ツールの選択が重要だ。ANSYS Polyflowは粘弾性流体（Oldroyd-B、Giesekus、FENE-Pモデル）と自由表面（ダイスウェル）解析に特化しており、高分子加工業界では事実上の標準ツールだ。一方ANSYSのFluentも非ニュートン粘性モデル（Power-Law、Carreau等）は実装しているが粘弾性モデルはPolyflowより限定的。OpenFOAMはviscoelasticFluidFoamソルバがあり、研究用途では粘弾性モデルの自由な実装・拡張が可能。実務では「粘弾性・自由表面が重要→Polyflow」「剪断粘度依存性のみ→Fluent/OpenFOAM」という棲み分けが効率的な選択だ。
選定で最も重要な3つの問い
- 「何を解くか」：非ニュートン流体に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
粘弾性流体の数値解析
Shear-thinning以上に難しい粘弾性流体の計算について教えてください。
粘弾性流体では応力テンソルが独自の発展方程式（構成方程式）に従う。代表的なモデルを見てみよう。
Oldroyd-Bモデル:
ここで $\overset{\triangledown}{\boldsymbol{\tau}}$ は上対流Maxwell微分、$\lambda_1$ は緩和時間、$\eta_p$ はポリマー粘度。
高Weissenberg数問題（HWNP）: Weissenberg数 $\text{Wi} = \lambda_1 \dot{\gamma}$ が大きくなると、応力の指数的増大によりソルバーが発散する。これは粘弾性CFDの長年の難問だ。
対策としては:
- DEVSS法: 粘性項を追加し、余分な項を構成方程式側で引き去る安定化手法
- Log-conformation法: 応力テンソルの対数を変数とすることで指数的増大を抑制。Fattal & Kupferman (2004) の提案
- sBETA法: ポリマー応力のブレンド安定化
Log-conformation法は対数変換で発散を防ぐんですね。数学的に美しい。
チキソトロピーのモデリング
チキソトロピー流体は粘度が「せん断の履歴」に依存する。構造パラメータ $\xi$（0〜1）を導入して表現する。
$a$ が構造回復速度、$b$ が構造破壊速度だ。コンクリートやインクジェットインクの解析で使われる。
時間の概念が入ると非定常計算が必須になりますね。
非ニュートン流体の乱流モデリング
非ニュートン流体の乱流は、ニュートン流体とは質的に異なる挙動を示す。
特に重要なのがドラッグリダクション（Drag Reduction）現象だ。微量のポリマー（ppm レベル）を添加するだけで、管路の摩擦損失が最大80%も減少する。Tom'sの効果として知られている。
この現象のDNS（直接数値シミュレーション）では、FENE-Pモデルを使った粘弾性乱流の計算が主流だ。ポリマーが渦構造を変化させ、ストリーク間隔が広がることで乱流摩擦が減少するメカニズムが解明されている。
ppmレベルの添加で80%も摩擦が減るんですか！ パイプラインの省エネに大きいですね。
アラスカ横断パイプラインでは実際にポリマー添加によるドラッグリダクションが使われている。CFDでの予測が設計に活用される好例だ。
最新の研究動向
非ニュートン流体CFDの最前線を紹介しよう。
- 機械学習による構成方程式の発見: レオメータデータからニューラルネットワークで構成方程式を自動生成（GENERIC framework + NN）
- 粒子法（SPH）による自由表面非ニュートン流れ: 溶岩流、食品加工のシミュレーション
- マイクロ流体デバイスにおける粘弾性不安定性: Wi数が臨界値を超えると弾性乱流が発生する新しい混合メカニズムの解析
- 生体流体力学: 赤血球の変形を考慮した血流シミュレーション
非ニュートン流体のCFDは工学だけでなく生体や地球科学にも広がっているんですね。
血流の非ニュートン性——赤血球凝集がもたらす剪断減粘性とCFD動脈解析
血液は剪断速度が低い（<10 s⁻¹）と赤血球が数珠つなぎ（Rouleaux形成）になり見かけ粘度が急上昇する典型的な剪断減粘性流体だ。Carreau-Yasudaモデルや修正Crossモデルで記述され、動脈瘤や冠動脈分岐部のCFD解析では非ニュートン効果が壁面剪断応力（WSS）予測に15〜20%の差をもたらすことが報告されている。低WSS領域(<0.4 Pa)が動脈硬化プラーク形成と相関するため、血流CFDの精度は医学的診断価値に直結する。最先端の「4D Flow MRI」（MRIで3次元時系列速度場を計測）データとCFDの融合による個人特化型血流解析が臨床応用の研究フロンティアだ。
トラブルシューティング
非ニュートン流体計算の典型的トラブル
非ニュートン流体の計算で、よくハマるポイントを教えてください。
非ニュートン特有のトラブルをパターン別に整理しよう。
1. 粘度が発散して計算が止まる
症状: 粘度が $10^{30}$ のような非現実的な値になり発散
原因: Power-lawモデルで $\dot{\gamma} \to 0$ のとき $\eta \to \infty$（$n < 1$ の場合）。流れの停滞領域で発生する。
対策:
- 粘度上下限を設定: Fluentでは
Minimum ViscosityとMaximum Viscosityを指定。一般に $\eta_{\text{min}} = \eta_\infty$、$\eta_{\text{max}} = \eta_0$ に設定 - Carreauモデルへ変更: ゼロせん断粘度 $\eta_0$ で自動的に上限が設定される
- 背景粘度の追加: Power-law粘度に微小な定数粘度を加える
Power-lawモデルの欠点がここで出てくるんですね。
2. 降伏応力モデルの収束困難
症状: Herschel-Bulkleyモデルで残差が振動し収束しない
原因: 降伏面（$\tau = \tau_y$）の位置が反復ごとに変動し、「固体⇔液体」の遷移が安定しない
対策:
- Papanastasiou正則化パラメータ $m$ を段階的に増やす: $m = 10 \to 100 \to 1000$ と徐々に上げる
- Bi-viscosityモデル: 降伏前に非常に高い粘度（$10^3 \sim 10^5\,\text{Pa}\cdot\text{s}$）を仮定
- 緩和係数を下げる: 粘度の under-relaxation を0.3-0.5に設定
3. せん断速度の計算精度不足
症状: 壁面近傍で粘度が不正確、圧力損失が理論値と合わない
原因: メッシュが粗く、壁面のせん断速度 $\dot{\gamma}$ が正確に計算できていない
対策:
- 壁面第一層を十分に細かくする（管内流なら $R/\Delta r > 20$ が目安）
- 二次精度以上の勾配スキームを使用
- 壁面に直交するプリズム層を配置
せん断速度の精度が粘度の精度に直結するんですね。
4. 温度依存粘度での熱暴走
症状: 粘性散逸により温度が上昇、粘度が下がりさらにせん断速度が増大するフィードバックループで発散
原因: 高粘度流体（高分子メルトなど）の高速せん断で粘性散逸が支配的になる
対策:
- エネルギー方程式のViscoeus Dissipation項を有効にする
- 時間刻みを十分に小さくする
- 粘度の温度依存係数を確認（Arrhenius型: $\eta = \eta_0 \exp[E_a/(RT)]$）
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——非ニュートン流体の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告