A Solid-Based Approach for Modeling Simple Yield-Stress Fluids: Rheological Transitions, Overshoot and Relaxation

本論文は、 Carbopol 分散液などの単純な降伏応力流体の定常および過渡的なレオロジー挙動（特に始動せん断における応力過渡現象）を、Zener 型粘弾性固体要素と非線形粘性モデルなどを組み合わせた新しい構成方程式によって記述し、その過渡現象が異方性硬化や不均一な微細構造進化ではなく、垂直応力差による応力不変量の増大に起因する均一なメカニズムであることを明らかにしたものである。

要約

この論文は、**「カボポル（Carbopol）」**という特殊な液体の動きを、新しい数学のルール（モデル）を使って説明しようとした研究です。

カボポルは、歯磨き粉、ヘアジェル、リチウムイオン電池の電極スラリーなどに使われている「降伏応力流体」という物質です。一言で言うと、**「静かには固いスポンジのようだが、力を加え続けると液体のように流れる」**という不思議な性質を持っています。

これまでの研究では、この物質の「ゆっくり流れる様子」や「急激に動かす様子」を同時に正確に予測するのが難しかったのですが、この論文はそれを解決する新しい「レシピ」を見つけ出しました。

以下に、専門用語を避けて、身近な例え話で解説します。

---

1. 研究のゴール：なぜ新しいモデルが必要なのか？

これまでのモデルは、この物質を「液体」として扱うことが多かったのですが、実験を見ると、実は**「固体」**としての側面が強いことがわかりました。

• 実験の不思議な点：
  • スタートダッシュ（スタートアップ）： 急に動かそうとすると、最初は抵抗が大きくて「ピーク（山）」になり、その後落ち着いて流れる（応力過剰という現象）。
  • 止めた後（緩和）： 動かすのをやめても、完全に力が消えず、**「少しだけ力が残ったまま」**になる。
  • クリープ（じわじわ変形）： 軽い力をかけると固まったままだが、ある限界を超えると急に流れ出す。

これまでのモデルは、これらの複雑な動きを一度に説明できませんでした。そこで、著者たちは**「この物質は元々『固体（ゼリー）』で、そこに『水（溶媒）』が混ざっている」**と考え、新しいモデルを作りました。

2. 新しいモデルの仕組み：「ゼリーと水」の合体

著者たちが考えたのは、以下の 3 つのパーツを組み合わせた機械のようなイメージです。

1. ゼリー部分（Zener モデル）：
   • これは**「バネ」と「ゼリー（バネ＋ダッシュポット）」**の組み合わせです。
   • イメージ： 硬いスポンジの中に、少し柔らかいゼリーが入っている状態。
   • 役割： 物質の「骨格」を作ります。力を加えると変形しますが、元に戻ろうとする力（弾性）と、ゆっくり変形する力（粘性）を持っています。
2. 水部分（並列のダッシュポット）：
   • イメージ： ゼリーの隙間を埋めている**「水」**です。
   • 役割： 全体が動くとき、この水も一緒に動きます。これにより、摩擦のような抵抗（粘性）が生まれます。
3. 「記憶」する力（バックストレス）：
   • イメージ： 変形したゼリーが**「元に戻りたい」という記憶**を持つ力です。
   • 役割： 変形が大きいほど、元に戻ろうとする力が強くなります。これが「降伏（流れ始める瞬間）」を決める鍵になります。

3. このモデルが解明した「3 つの謎」

この新しい「ゼリー＋水」のモデルを使うと、実験で見つかった不思議な現象がすべて説明できました。

① なぜ「スタートダッシュ」でピークが出るのか？（応力過剰）

• 現象： 急に動かすと、最初は「ぐっと」力がかかり、その後「ふっと」力が抜けて安定する。
• 解説： 1 次元（単純な計算）では説明できませんでしたが、3 次元（実際の立体空間）で計算すると、ゼリーが歪む方向によって「内部の圧力」が変化します。
  • 例え： 硬いスポンジを横から押すと、最初は固く感じますが、ある角度で歪むと、内部の構造が崩れて急に柔らかくなる瞬間があります。この「歪みの方向による変化」が、ピーク（山）を生み出しました。
  • 重要： これまで「内部構造がバラバラになるから」とか「硬くなるから」と言われていましたが、このモデルは**「均一なゼリーが歪むだけでピークが出る」**ことを証明しました。

② 止めても力が残る理由（非ゼロの完全緩和応力）

• 現象： 動かすのをやめても、力がゼロにならず、少し残る。
• 解説： 液体なら力がゼロになりますが、これは**「固体（ゼリー）」**だからです。
  • 例え： 伸びたゴムを放しても、完全に縮む前に止まることがあります。この物質は、液体ではなく「流れることができる固体」なので、力をやめても「元に戻ろうとする力」が少し残るのです。

③ 固いのか、液体なのか？（クリープ転移）

• 現象： 軽い力なら固まったまま、重い力なら流れ出す。
• 解説： このモデルには**「流れ始める限界の力（降伏応力）」**がはっきり定義されています。
  • 例え： 氷に軽い指を乗せれば固いですが、重い鍋を乗せれば溶けて流れます。このモデルは、その「境目」を正確に予測できます。

4. この研究のすごいところ

• 固体から見た視点： これまで「液体」として扱われていたものを、**「固体（ゼリー）」**の視点からモデル化した点が画期的です。
• シンプルで正確： 複雑なパラメータを大量に使う必要がなく、シンプルながら実験結果（ピーク、残留応力、流れの境目）をすべて再現できました。
• 実用性： このモデルを使えば、電池の製造工程や、3D プリンターでインクを押し出すような、複雑な流れのシミュレーションがより正確に行えるようになります。

まとめ

この論文は、**「カボポルという不思議な物質は、実は『記憶を持つゼリー』と『水』の組み合わせでできている」**という新しい視点を提供しました。

これまでの「液体」という考え方に固執せず、「固体」として捉え直すことで、なぜ急に動かすとピークが出るのか、なぜ止めても力が残るのか、といった長年の謎を、シンプルで美しい数学のルールで解き明かしました。これは、新しい材料の開発や製造プロセスの最適化に大きく貢献するはずです。

技術概要

以下は、提出された論文「A Solid-Based Approach for Modeling Simple Yield-Stress Fluids: Rheological Transitions, Overshoot and Relaxation（単純な降伏応力流体のモデル化のための固体ベースのアプローチ：レオロジー的遷移、オーバーシュート、緩和）」の技術的サマリーです。

1. 研究の背景と課題 (Problem)

降伏応力流体（Carbopol 分散液や電池用電極スラリーなど）は、自然現象から産業応用まで広く存在しますが、そのレオロジー挙動を記述する構成方程式の開発には依然として課題が残っています。

• 既存モデルの限界: 従来のモデル（Herschel-Bulkley 式など）は定常せん断流れを記述できますが、過渡的挙動（スタートアップせん断、クリープ、応力緩和）を統一的に予測することは困難です。
• 未解決の現象:
  • 応力オーバーシュート: スタートアップせん断流れで観測される応力のピーク現象を、単純な降伏応力流体（チキソトロピーが弱い Carbopol など）で説明する際、従来のモデルは「等方性硬化」や「空間的不均一な微構造進化」を仮定する必要があり、物理的メカニズムの解明が不十分でした。
  • 完全緩和応力: 応力緩和実験において、応力がゼロにならず有限の値で飽和する現象（固体特性）を、既存の流体ベースのモデルが正確に再現できていませんでした。
  • クリープ遷移: 応力大小によって、ひずみ飽和（固体挙動）から時間とともに増加する流動（流体挙動）へ遷移する現象の統一的な記述が求められていました。

2. 提案手法とモデル化 (Methodology)

本研究では、降伏応力流体を「流体」ではなく「粘弾性固体」の枠組みで記述する新しい構成方程式を提案しました。Carbopol 分散液の微構造（高濃度の水和マイクロゲル粒子と間隙溶媒）に基づき、以下の要素を組み合わせています。

• 基礎的な力学モデル:
  • ゼナー型粘弾性固体要素: ゲル相の応答を記述するために、バネ（弾性）とケルビン・フォイト要素（バネとダッシュポットの直列）を直列に接続。
  • 並列ダッシュポット: 間隙溶媒による粘性散逸を表現するため、ゼナー要素に並列に線形ダッシュポットを追加。
• 変形分解:
  • 1 次元モデル: 加法的ひずみ分解（全ひずみ = 弾性ひずみ + 塑性ひずみ）を使用。
  • 3 次元モデル: 有限変形と物質フレーム不変性の原則を満たすため、Kröner–Lee 分解（変形勾配の乗法的分解：$F = F^e \cdot F^p$）を採用。
• 構成則の要素:
  • 流動則: 塑性ひずみ速度を記述する非線形粘性モデル（Eyring 粘性モデルまたは Carreau-Yasuda 粘性モデル）を採用。降伏応力を超えた際のせん断希釈性を表現。
  • 背応力（Back Stress）の進化: 塑性変形に伴う硬化/軟化を記述するために、Armstrong-Frederick 理論に基づく背応力の進化方程式を導入。
  • 応力の構成: 全応力は、ゼナー要素による「ゲル応力」と、並列ダッシュポットによる「溶媒応力」の和として定義。

3. 主要な成果と結果 (Key Results)

1 次元モデルの結果

• 定常状態: 低せん断速度域で Herschel-Bulkley 型の挙動を再現し、高せん断速度域では第二ニュートン領域への遷移を予測。
• 応力緩和: 応力がゼロにならず、有限の「完全緩和応力」に収束することを示した。
• 限界: 1 次元モデルでは、スタートアップせん断における応力オーバーシュートを再現できなかった。これは、粘性が応力差の単調関数であるため、弾性負荷項と緩和項の競合が単調に収束してしまうため。

3 次元モデルの結果（本論文の核心）

3 次元への拡張により、以下の重要な現象を定量的・定性的に再現することに成功しました。

1. 応力オーバーシュートの発生メカニズムの解明:
   • 3 次元モデルでは、せん断流れ下でも法線応力差が発生し、これが応力不変量（$\xi$）に寄与します。
   • この法線応力差が粘性モデル（せん断希釈性）を加速させ、結果として塑性変形速度（緩和項）が一時的に弾性負荷項を上回ります。
   • 重要点: このオーバーシュートは、等方性硬化や微構造の空間的不均一性を仮定せず、均質な粘弾性固体のテンソル構造から自然に導かれることを示しました。
2. 応力緩和挙動の改善:
   • 3 次元モデルでは、完全緩和応力が前せん断速度に依存して減少する現象を再現し、実験事実（Vinutha et al., 2024）と一致しました。
3. クリープ遷移の予測:
   • 印加応力が Herschel-Bulkley 式で得られる降伏応力（$\tau_y$）を超えると、ひずみ率がゼロに減衰せず有限値に収束し、固体から流体への遷移が明確に予測されました。
4. 定常せん断応力:
   • 低せん断速度域では Herschel-Bulkley 則に従い、高せん断速度域ではモデル固有の曲率変化を示すことを確認しました。

4. 学術的・実用的意義 (Significance)

• 理論的革新: 降伏応力流体を「流体」ではなく「固体（降伏して流動する固体）」としてモデル化するアプローチの有効性を示しました。特に、応力オーバーシュートが「均質なメカニズム（法線応力差による粘性変化）」から生じることを初めて示唆し、従来の複雑な仮定（微構造不均一性など）を不要にしました。
• 汎用性の向上: Kröner-Lee 分解に基づく 3 次元構成方程式は、有限要素法（FEM）などの数値シミュレーションへの適用が可能であり、電池電極スラリーの塗布やダイレクトインク書き出し（DIW）など、複雑な変形モードを伴う産業プロセスの設計・解析に貢献します。
• 実験との整合性: Carbopol 分散液の実験データ（スタートアップ、クリープ、応力緩和）の多様な特徴を、単一の構成方程式で統一的に説明できることを実証しました。

結論

本研究は、降伏応力流体の複雑な過渡的挙動（特に応力オーバーシュートと完全緩和応力）を、微構造に基づいた「粘弾性固体＋並列溶媒散逸」のモデルから導き出すことに成功しました。このアプローチは、従来の流体ベースのモデルの限界を克服し、より物理的に妥当で実用的な構成方程式の枠組みを提供するものです。