理論メモ
べき乗則:$\tau = K\dot{\gamma}^n$ ビンガム:$\tau = \tau_0 + \mu\dot{\gamma}$ ($\tau > \tau_0$)ハーシェル-バルクレー:$\tau = \tau_0 + K\dot{\gamma}^n$
アレニウス温度補正:$\mu(T) = \mu_0 \exp\!\left(\frac{E_a}{RT}\right)$
ニュートン流体からビンガム体・べき乗則流体まで、4種のレオロジーモデルの粘度曲線とクエット流れ速度プロファイルをリアルタイム可視化。温度補正・レイノルズ数計算も対応。
レオロジーモデルの基本は、せん断応力 $\tau$ とせん断速度 $\dot{\gamma}$ の関係式(流動曲線)で表されます。ニュートン流体は最も単純な比例関係です。
$$\tau = \mu \dot{\gamma}$$ここで、$\mu$ は粘度 [Pa·s] で一定値です。$\tau$ はせん断応力 [Pa]、$\dot{\gamma}$ はせん断速度 [1/s] です。
非ニュートン流体を代表する「べき乗則(Power Law)モデル」は、粘度がせん断速度に依存することを表します。
$$\tau = K \dot{\gamma}^n$$$K$ は整合性指数 [Pa·sn]、$n$ はべき乗指数(流動指数)[-] です。$n < 1$ のときずり流動化(せん断速度が上がると見かけ粘度が下がる)、$n > 1$ のときずり濃化を示します。$n=1$のときニュートン流体に帰着します。
食品・化粧品工業:マヨネーズ、ケチャップ、歯磨き粉などの製品設計では、容器からの出やすさや塗布時の感触が重要です。ビンガムモデルやハーシェル・バルクレーモデルを用いて、降伏応力(流れ始めるのに必要な最小応力)を評価・制御します。
高分子・プラスチック加工:プラスチックを金型に射出成形する際、溶融樹脂は強いせん断を受けるため、その粘度変化(ずり流動化)を正確に把握する必要があります。べき乗則モデルはこのような高分子溶融体の流動解析の基礎データとして広く用いられています。
塗料・インク開発:塗料は、刷毛やローラーで塗る時は低粘度で広がりやすく(ずり流動化)、塗布後は高い粘度で垂れないことが求められます。このような「チキソトロピー」性の評価にレオロジー測定とモデル化が不可欠です。
CFD(数値流体力学)シミュレーション:血液、泥水、セメントスラリーなどの非ニュートン流体の流れを予測するため、本ツールで確認できるようなレオロジーモデルとそのパラメータが、OpenFOAMや商用CFDソフトへの入力データとして使用されます。レイノルズ数の計算は流れの状態を予測する第一歩です。
このツールを使い始めるとき、いくつかハマりやすいポイントがあるから気をつけてね。まず一つ目は、「べき乗則の指数nだけで流体の性質が全部わかる」と思ってしまうこと。確かにn<1ならずり流動化だけど、それはあくまで「せん断速度が増えると粘度が下がる傾向がある」というだけ。実際の材料、例えば高分子溶融体は、せん断速度が極端に低い領域や高い領域ではべき乗則から外れることが多いんだ。ツールでn=0.3とn=0.8を比べてみて? 傾向は似てるけど、粘度の低下の仕方は全然違うだろ? 実務では、測定した全データポイントにフィットする単一のモデルは稀で、領域ごとにモデルを切り替えることもあるんだ。
二つ目はレイノルズ数の解釈。ツールではパイプ流れを想定して計算してるけど、この値は「目安」でしかない。例えば、レイノルズ数が2300を超えたからといって、必ず乱流になるわけじゃない。流路の形状や入口条件で大きく変わるからね。あくまでCFDのメッシュ分割や解析手法を選ぶための初期判断材料として使おう。
最後に、「ビンガム流体=ケチャップ」と短絡的に結びつけないこと。確かにケチャップは降伏応力を持つけど、時間依存性(チキソトロピー)も強い。ツールのビンガンモデルは「一度流れ始めたらニュートン流体のように振る舞う」超理想化されたモデルだ。実物はもっと複雑で、かき混ぜてから時間が経つとまた固まってくる。ツールで学ぶのは「降伏応力」という概念の入り口。実際の製品設計では、もっと複雑なモデルが必要になることを頭の片隅に入れておいて。
このシミュレーターで扱っている粘度やレオロジーの概念は、思っている以上に幅広い分野の根幹を支えているんだ。例えばバイオメカニクスでは、血液を非ニュートン流体(具体的にはカッソンモデルなど)としてモデル化して、動脈瘤の血流シミュレーションや人工心臓の設計に応用されている。血管の細い部分では血液の見かけ粘度が上がる(ずり流動化)現象が、血栓リスクの評価に重要なんだ。
もう一つは地盤工学や土木。地震時の液状化現象や、土石流の挙動予測では、泥水や砂-水混合物をビンガム流体のように扱うことがある。降伏応力があって、それを超える力が加わると一気に流動化するからだ。防災対策のための数値シミュレーションの基礎データとして、レオロジーパラメータの推定が行われているよ。
3Dプリンティング(アディティブマニュファクチャリング)も熱い応用分野だ。材料を積層する方式(例えば材料押出し)では、ノズルから押し出される時の材料の流動性と、積層後の形状保持性が天秤になる。これはまさに「ずり流動化(押出し時は低粘度)」と「降伏応力(積層後は高粘度で垂れない)」の制御問題なんだ。ツールでいろんなモデルをいじる感覚が、この最適化プロセスに直結する。
ツールに慣れて「もっと知りたい」と思ったら、次のステップに進んでみよう。まずは数学的な背景を少し深堀りするのがおすすめ。ツールで出てきた「せん断速度」は、速度の空間微分 $\dot{\gamma} = \frac{du}{dy}$ で定義される。これは流体の中の速度の「勾配」だ。この勾配が応力を生み、その比例「係数」が粘度$\mu$。非ニュートン流体では、この係数が勾配の関数 $\mu(\dot{\gamma})$ になる、と理解するとすっきりするよ。
次に挑戦したいのは「一般化ニュートン流体(GNF)」という概念の理解だ。ツールのべき乗則やビンガムモデルは全てGNFの一種。要は「粘度がせん断速度のスカラー関数で表せる流体」という大きな枠組みがあるんだ。この概念を押さえると、CFDソフトの材料設定画面で「Non-Newtonian」のサブメニューに並ぶたくさんのモデルが、全て同じ原理で動いていることがわかってくる。
最終的には、時間依存性(チキソトロピー・レオペクシー)や粘弾性というさらに広い世界に踏み込もう。歯磨き粉は力を加えるとサラサラになるが、力を止めると元の硬さに戻る(チキソトロピー)。これはツールの静的なモデルでは表現できない「経過時間」の因子が入る。粘弾性は、ゴムのように「固体の弾性」と「液体の粘性」を併せ持った挙動で、緩和時間などのパラメータが必要になる。これらの現象は、製品の「使い心地」や「加工性」を決める極めて重要な要素で、レオロジー解析の醍醐味とも言える深い領域だ。