End-pinching and inertial-capillary reopening in viscoplastic ligaments at low Ohnesorge number

本論文は、低オネスゴ数における Herschel-Bulkley 流体の液柱再閉合メカニズムを解明し、せん断希薄化領域において粘性がゼロの極限でも端絞りによる破断ではなく、曲率誘起圧力勾配に支配される純粋な慣性 - 毛管力による再開口が生じることを示した。

要約

この論文は、**「粘り気のある液体の糸（リグメント）が、表面張力で縮むとき、どうやって切れるか（あるいは切れないか）」**という現象を、新しい視点から解明した研究です。

専門用語を避け、身近な例え話を使って解説します。

1. 物語の舞台：「液体の糸」と「魔法のハサミ」

想像してください。太い糸のような液体（例えば蜂蜜やケチャップ）が空中に浮かんでいます。この液体は、表面張力という「縮みたがる力」によって、自分自身を丸めようとしています。

• 普通の水（ニュートン流体）の場合：
  水のような液体は、縮むと先端が太くなり、その根元が細く絞られていきます。まるで**「首を絞められるように」細くなり、最後はポロッと一滴が飛び去ります。これを「エンド・ピンチング（端の絞め切り）」**と呼びます。

  • これまでの常識： 「水は粘度が低ければ低いほど、この『絞め切り』が起きやすくなる。粘度ゼロ（完全にサラサラ）なら、必ず切れるはずだ」と考えられていました。

• 今回の研究のテーマ：
  しかし、現実には「ケチャップ」や「歯磨き粉」のように、**「ある程度力を入れないと動かない（降伏応力）」かつ「動く速さによって粘度が変わる（せん断速度依存性）」という、複雑な性質を持つ液体（ヒュースキー・バルキー流体）があります。
  この論文は、「そんな複雑な液体の糸が縮むとき、どんなドラマが起きるのか？」**をシミュレーションで詳しく調べました。

2. 発見された 4 つの「結末」

研究チームは、液体の性質（「力を入れれば動くか」「動く速さで粘度がどう変わるか」）を変えながらシミュレーションを行い、4 つの異なる結末を見つけました。

① 定番の「絞め切り」（エンド・ピンチング）

• どんな時？ 粘度が低めで、降伏応力も小さい場合。
• 様子： 普通の水と同じように、先端が太くなり、根元が細くなってポロッと切れます。
• 例え： 飴玉を引っ張って細くし、最後にはちぎれるような感じ。

② 奇跡の「逃げ出し」（エンド・ピンチングからの脱出）

ここがこの論文の最大の発見です。液体が「絞め切り」を回避して、切れないで丸くなるパターンです。実は2 通りの逃げ方がありました。

• A. 「太ることで逃げる」場合（せん断増粘）

  • 仕組み： 液体が「細くなろうとすると、急に固くなる（粘度が上がる）」性質を持っている時です。
  • 例え： 糸の根元が細くなろうとすると、その部分の液体が**「固まってゴムのように硬くなる」**のです。すると、細くなるのを防いで、逆に渦を巻いて流れが戻り、糸が太り直します。
  • イメージ： 絞められそうになった瞬間、首に「硬い防具」ができて、絞め殺されずに生き延びる。

• B. 「圧力で押し戻す」場合（せん断減粘）

  • 仕組み： 液体が「細くなると、逆にサラサラになる（粘度が下がる）」性質を持っている時です。
  • 意外な事実： 以前は「粘度がゼロに近づくと、必ず切れる」と思われていました。しかし、この研究では**「粘度が極端に低い（ほぼ水に近い）状態でも、切れないで太り直す」**現象が見つかりました。
  • 例え： 糸の根元が細くなると、その部分の圧力が複雑に変化し、**「押し戻す力」**が働きます。まるで、細くなろうとする糸を、見えない手が「太く戻せ！」と押しているようです。
  • 重要： これは「粘性」のおかげではなく、「慣性」と「表面張力」のバランスだけで起きる現象で、これまで誰も気づいていなかった新しい法則です。

③ 「首なし」の丸まり（No-neck）

• どんな時？ 降伏応力（動くための最低限の力）が少し大きい場合。
• 様子： 糸の根元が細く絞られる（首ができる）前に、液体全体が均一に太くなり、ドロっとした玉になります。
• 例え： 細い糸が絞められる前に、**「太いロープ」**として丸まってしまいます。首が作られる前に、全体が「動かない部分」に覆われてしまうため、細く絞られずに済みます。

④ 完全な「凍結」（Motionless）

• どんな時？ 降伏応力が非常に大きい場合。
• 様子： 表面張力で縮もうとしても、液体が「固まりすぎて」全く動きません。
• 例え： 糸が**「コンクリート」**になってしまい、縮もうとしても動かない。表面張力という力が、固まった液体を動かすには弱すぎるのです。

3. この研究がなぜ重要なのか？

これまでの常識では、「液体がサラサラになればなるほど、糸は必ず切れる（エンド・ピンチングする）」と考えられていました。しかし、この研究は**「それは違う！粘度が極端に低い状態でも、液体は『太り直す』という別の道を選ぶことができる」**と証明しました。

• 実社会への応用：
  • インクジェット印刷： きれいな一滴を飛ばすには、糸が切れるタイミングを制御する必要があります。
  • 農業スプレー： 農薬が霧状になる仕組み。
  • 3D プリンティング： 柔らかい材料を積層する技術。

これらの技術において、「液体がどう振る舞うか」を正確に理解できれば、より精密な制御が可能になります。特に、ケチャップや化粧品、生体材料など「複雑な液体」を扱う際、この「太り直す」という新しい現象を考慮に入れることで、失敗を防げるかもしれません。

まとめ

この論文は、**「液体の糸が縮むとき、必ず切れるとは限らない」**という新しい世界を見せました。

• 固まる性質があれば、**「防具」**になって逃げられる。
• サラサラになりすぎる性質があれば、**「押し戻す力」**が働いて逃げられる。
• 固まりすぎれば、**「凍りついて」**動かない。

液体の振る舞いは、私たちが思っているよりもずっとドラマチックで、多様な結末があることを教えてくれました。

技術概要

この論文は、低オズノル数（Ohnesorge number, $Oh_K$）における粘塑性液柱（liquid ligament）の毛管再収縮（capillary retraction）における「端絞り（end-pinching）」現象と、その回避・再開口（reopening）のメカニズムを解明した研究です。ニュートン流体ではよく理解されている現象ですが、降伏応力とせん断速度依存性粘度を併せ持つヘルシュケル・バルクレー（Herschel-Bulkley）流体における挙動は未解明な部分が多く、本研究は数値シミュレーションを通じてそのダイナミクスを詳細に解析しています。

以下に、問題設定、手法、主要な貢献、結果、および意義について詳細にまとめます。

1. 問題設定と背景

• 背景: 液柱の毛管による再収縮は、インクジェット印刷、スプレーコーティング、農業用散布など、多くの産業プロセスにおいて重要な現象です。ニュートン流体では、低粘度領域において液柱の先端が細くなり、液滴が分離する「端絞り（end-pinching）」メカニズムが支配的であることが知られています。
• 課題: 既存の研究は、界面張力勾配によるマランゴニ効果や、単純な降伏応力流体（ビinghamモデル：降伏後の粘度が一定）に焦点を当てたものが中心でした。しかし、降伏応力とせん断速度依存性粘度（シアーシニング/シアーサンニング）の両方を持つヘルシュケル・バルクレー流体の挙動、特に低粘度（低 $Oh_K$）領域での挙動は十分に理解されていませんでした。
• 目的: ヘルシュケル・バルクレー流体の液柱が再収縮する際、降伏応力とせん断依存性粘度がどのように液滴の分離（ピンチオフ）を抑制し、あるいは新たな再開口メカニズムを生み出すかを解明すること。

2. 手法

• 数値シミュレーション: オープンソースソルバー「Basilisk」を用いた完全解像度の軸対称シミュレーションを行いました。
• 支配方程式: 非圧縮性 2 相ナビエ - ストークス方程式を「1-流体定式化（one-fluid formulation）」で解きました。
• 構成則: 液柱はヘルシュケル・バルクレー流体としてモデル化し、周囲はニュートン流体（空気）としました。
  • 有効粘度 $\mu_{eff}$ は、降伏応力 $\tau_y$ とせん断速度依存性（流挙動指数 $n$）を含む式で定義され、数値的な特異性を避けるため $\epsilon$-正則化（regularization）を適用しました。
  • 無次元パラメータとして、降伏応力と毛管力の比である「プラストキャピラリ数 $J$」と、慣性・毛管力対粘性力の比である「オズノル数 $Oh_K$」を定義しました。
• 検証: ニュートン流体の端絞り現象、ビingham流体、ヘルシュケル・バルクレー流体の定常流（ポアズイユ流）の解析解との比較により、ソルバーの精度と正則化パラメータの影響を十分に検証しました。

3. 主要な貢献と発見

本研究は、ヘルシュケル・バルクレー流体の再収縮において、以下の4 つの異なるダイナミクス領域を特定し、その境界を $(n, J)$ 平面にマッピングしました。

A. 端絞り回避の 2 つのメカニズム

ニュートン流体では低粘度で端絞りが起こると考えられてきましたが、本研究では以下の 2 つの異なるメカニズムにより端絞りを回避（再開口）することが示されました。

1. シアーサンニング（Shear-thickening, $n > 1$）領域:

   • メカニズム: 頸部（ネック）が細くなる過程で局所的なせん断速度が増加し、見かけの粘度が急激に上昇します。これにより、曲率の高い頸部から強い渦（vorticity）が剥離し、界面に沿った逆流（backflow）が発生します。
   • 結果: この逆流が毛管特異性（capillary singularity）に抵抗し、頸部の再開口を引き起こします。このメカニズムは、ビingham流体には見られない、粘度の空間的変化に起因するものです。

2. シアーシニング（Shear-thinning, $n < 1$）領域およびニュートン極限 ($Oh_K \to 0$):

   • メカニズム: 予想に反し、強いシアーシニング流体や、極めて低粘度のニュートン流体（$Oh_K \to 0$）においても端絞りは回避されます。これは、渦の生成による粘性抵抗ではなく、**曲率誘起圧力勾配（curvature-induced pressure gradients）**に起因します。
   • 詳細: 頸部の細化に伴い、慣性力と毛管力のバランスが変化し、軸方向の圧力勾配が非単調になります。これにより、頸部付近で流体が押し戻され、再開口が発生します。
   • 画期的な発見: 従来の研究では $Oh_K \approx 10^{-4}$ まで端絞りが起こると考えられていましたが、本研究は $Oh_K \to 0$ の極限において、粘性が支配的になる前に「慣性 - 毛管再開口（inertial-capillary reopening）」が起こることを初めて示しました。これは、無粘性極限でも端絞りが唯一の帰結ではないことを意味します。

B. その他の領域

• ノネック（No-neck）領域: 降伏応力が中程度（$J \approx 0.05$）の場合、頸部が完全に降伏せず、局所化が抑制されます。その結果、液柱は頸部を形成せずに滑らかに単一の液滴へと収縮します。
• 静止（Motionless）領域: 降伏応力が非常に大きい（$J \gtrsim 0.5$）場合、毛管力が降伏応力を上回ることができず、液柱は完全に静止します。

4. 結果の定量的分析

• 相図の作成: $n$（流挙動指数）と $J$（プラストキャピラリ数）の関数として、端絞り、回避、ノネック、静止の 4 つの領域を明確に区分けする相図を作成しました。
• スケーリング則: 端絞り回避の閾値や、再開口時の頸部半径の時間発展を解析し、降伏応力と幾何学的形状（曲率）のバランスから理論的に説明しました。
• 数値的頑健性: 格子解像度、正則化パラメータ、時間刻みを変化させても結果が変化しないことを確認し、発見された再開口現象が数値アーティファクトではないことを証明しました。

5. 意義と結論

• 理論的意義: 従来のビinghamモデル（定数粘度）では捉えられなかった、ヘルシュケル・バルクレー流体特有の「粘度の空間的変化」が流体力学的挙動に与える質的な影響を初めて明らかにしました。
• ニュートン流体への示唆: 低オズノル数極限におけるニュートン流体の挙動に関する従来の見解（端絞りが漸近的な結果である）を修正し、慣性 - 毛管バランスに基づく再開口メカニズムが存在することを示しました。
• 応用: インクジェット印刷やスプレー技術において、粘性流体の液滴形成を制御する新たなパラメータ（降伏応力とせん断依存性）の重要性を指摘しています。特に、端絞りを抑制して液滴の分離を防ぐための戦略として、流体のレオロジー特性を設計する道を開きました。

総じて、この論文は、降伏応力流体の複雑なレオロジーが、自由表面流における特異な現象（端絞り）をどのように変容させるかを包括的に解明し、流体ダイナミクスとレオロジーの交差点における重要な知見を提供するものです。