Self-similar breakup of a liquid ligament with a solid particle

本研究は数値シミュレーションと解析的モデリングを通じて、固体粒子が伸長する液体リグメントにおいて普遍的かつ自己相似的なピンチオフダイナミクスを誘起し、その後の破断が粒子サイズに依存せず、リグメントの伸長とレイリー・プレトー不安定性の相互作用によって支配されることを示している。

要約

想像してみてください。空中に長い細い蜂蜜や濃厚なシロップの糸が垂れ下がっている様子を。その両端を引っ張ると、糸は次第に細くなり、最終的にはちぎれて別々の滴に分かれます。これは、葉から落ちる雨滴から、インクジェットプリンターが微細なドットを噴射するまで、自然や技術においてよく見られる光景です。

通常、このちぎれは液体が本質的に不安定であるために起こります。エネルギーを節約するために球体（滴）になろうとするからです。しかし、もしその粘り気のある糸の中に、砂粒やほこり粒子のような小さな固体の粒が閉じ込められていたらどうなるでしょうか？

この論文は、まさにそのシナリオを調査しています。研究者たちは、コンピュータシミュレーションと数学を用いて、単一の固体粒子が伸びる液体の糸のちぎれ方をどのように変化させるかを確認しました。

以下に、彼らの発見を簡単な概念に分解して物語ります。

設定：結び目のある伸びる紐

液体の糸を蜂蜜でできた伸縮性のある長い紐だと考えます。研究者たちはこの紐の両端を一定の速度で引っ張りました。紐の中には、単一の固体の球（粒子）を配置しました。

最初は紐が太く、球は紐の幅に比べて小さく見えます。太い庭用ホースの中にビー玉が入っているようなものです。ビー玉は何もせず、紐は伸びるにつれて予測可能なパターンに従って次第に細くなっていきます。

転換点：「ホース」が「ビー玉」のサイズまで縮むとき

紐が伸び続け、細くなっていきます。やがて、紐は内部のビー玉の表面にほぼ触れるほどに細くなります。

これが決定的な瞬間です。論文では、粒子のサイズと紐のサイズの比率が 1 に近づくときを指します。突然、ビー玉は紐の中の「結び目」や「膨らみ」として機能します。紐があまりに細いため、この膨らみは局所的な擾乱を生み出します。

驚き：「普遍的」なちぎれ

ここが発見の中で最も興味深い部分です。研究者たちは、異なるサイズのビー玉（小さいものから大きいものまで）でこの実験を行いました。

• ちぎれる前： 大きいビー玉は、小さいビー玉よりも早く紐をちぎらせました。これは理にかなっています。より大きな障害物はより早く問題を引き起こすからです。
• ちぎれる瞬間： 紐がビー玉に触れるほどに細くなると、不思議なことが起こりました。最終的なちぎれが起こる速度は、ビー玉が小さかろうと大きかろうと、完全に同じになったのです。

研究者たちはこれを**「自己相似」**的な振る舞いと呼んでいます。紐が障害物に触れるほどに細くなると、障害物の具体的なサイズはもはや重要ではなくなるかのように、液体は粒子の大きさを「忘れ」、ちぎれに至る普遍的で予測可能な経路をたどります。

比喩：交通渋滞

道路（液体の糸）を想像してください。車同士が互いに離れて走行し、交通が広がっている（伸びている）状態です。

• 初期段階： 道路の真ん中に小さな穴（小さな粒子）や大きな岩（大きな粒子）があっても、道路が広いのであまり問題になりません。
• 後期段階： 道路が一本の車線に狭まると、穴も岩もどちらも巨大な障害物となります。
• ちぎれ： 交通が障害物に衝突するほどに圧迫された瞬間、交通が麻痺して停止する（「分断」する）様子は、穴と岩の両方に対して全く同じように起こります。障害物のサイズは、最終的な渋滞のタイミングを変えなくなります。重要なのは、そこに何かがあるという事実だけなのです。

数学と物理学

研究者たちは単にこの現象を観察しただけではありません。いつちぎれるかを正確に予測する数学的な式を立てました。

• 彼らは、ちぎれる時間は伸び（紐を引っ張る力）と粘性（液体がどのくらい「太い」または粘り気があるか）という二つの力のせめぎ合いに依存することを発見しました。
• 蜂蜜のような粘り気の強い液体では、「粘り気」が支配的になります。
• 彼らの式は、コンピュータシミュレーションと完全に一致する形で、ちぎれる時間を正確に予測することに成功しました。

結論

この論文は、粒子が（粒子の近くで紐をより早く細くさせることで）分断がいつ始まるかを変化させる一方で、紐が粒子に触れるほどに細くなると、最終的なちぎれの行為は普遍的な法則に従うと結論付けています。

この特定の「粘り気のある」領域において、液体の糸は、ちぎれるための決まったプログラムを持つ機械のように振る舞います。粒子がそのプログラムをトリガーするほどに近づくと、粒子のサイズは無関係となり、糸は毎回予測可能で自己相似的な方法でちぎれます。

技術概要

技術的サマリー：固体粒子を有する液体リグメントの自己相似的破断

問題提起
薄化する液体リグメントの破断は、スプレー形成やインクジェット印刷からマイクロ流体滴生成に至るまで、多様な応用における基本的な過程である。外部の伸長がない場合、この破断はレイリー・プレート不安定によって支配され、表面の擾乱が表面張力によって増幅する。そのダイナミクスは通常、粘性力と慣性力および表面張力との比を表すオネスオーゲ数（$Oh$）によって特徴づけられる。清浄なリグメントの挙動はよく理解されているが、多くの実用的なシステムには固体粒子や不純物を含む懸濁液が含まれる。先行研究では多数の粒子を有するジェットや高密度懸濁液を対象とした粒子含有ジェットが検討され、変化した破断パターンが示されている。しかし、単一の固体粒子が局所的な擾乱として作用し、伸長するリグメントの破断ダイナミクスおよびピンチオフ時間をどのように変化させるかという基本的なメカニズムは、依然として十分に理解されていない。具体的には、リグメントの伸長、粘性、および局所的な粒子誘起擾乱との相互作用が、最終的な破断領域をどのように決定するかは不明である。

手法
著者らは、この問題を調査するために数値シミュレーションと解析的モデリングの組み合わせを採用している。

• 数値設定: 本研究は、二相非混和流体に対するカラーグラデーションアプローチに基づく格子ボルツマン法（LBM）を利用している。システムは、初期半径$R_i$および長さ$L$（アスペクト比$L/R_i = 28.8$）の円柱状液体リグメントを空気中にモデル化している。半径$r_p$の固体球形粒子が、すべりなし境界条件の下でリグメントの中心に配置される。
• 伸長メカニズム: 軸方向の伸長は、リグメントの端に線形速度プロファイル（$u_x = \frac{2}{L}(x - L/2)u_{max}$）を適用することで課され、これにより質量が領域から流出し、リグメント半径が時間とともに減少する。伸長率は、無次元流出率$\tilde{\dot{\epsilon}}$によって特徴づけられる。
• パラメータ: シミュレーションは、$Oh \sim 1.0$の粘性領域に焦点を当てている。粒子サイズは初期リグメント半径に対する相対値（$r_p/R_i$）で変化させられ、流出率は、伸長時間スケールを表面張力 - 粘性時間スケール（$\tau_v = \mu R_i / \sigma$）に対して制御するために調整される。
• 解析的アプローチ: 著者らは、リグメント半径を質量保存によって支配される非擾乱平均半径と、粒子によって誘起される局所的な擾乱振幅の和としてモデル化することで、ピンチオフ時間に対する解析式を導出した。彼らは、$Oh \sim 1$の領域で粘性効果が支配的であるレイリー・プレート不安定のダイナミクスに従って擾乱が増幅すると仮定している。

主要な結果

1. 2段階ダイナミクス: 破断過程は2つの明確な段階を示す。初期には、リグメント半径$R(t)$が粒子半径（$r_p/R(t) < 1$）よりも著しく大きい場合、リグメントは質量保存に従って単調に薄化し、粒子の影響は無視できる。
2. 自己相似性の開始: 伸長が進みリグメント表面が粒子に接近すると（$\beta(t) = r_p/R(t) \sim 1$）、局所的な擾乱が誘起される。この接近が達成されると、その後の擾乱の増幅および最終的なピンチオフダイナミクスは粒子サイズに依存しない自己相似的なものとなる。
3. 普遍的なピンチオフ時間: 破断までの総時間（$t_p$）は粒子サイズに依存する（より大きな粒子はより早期の破断を誘起する）が、擾乱増幅の開始から破断までの期間として定義される自己相似的ピンチオフ時間（$t_{ps} = t_p - t_{exp}$）は、与えられた伸長率において異なる粒子サイズに対して単一の値に収束する。
4. 解析的検証: 著者らは、伸長率と粘性表面張力不安定増幅率（$\lambda = \sigma / \mu R(t)$）との競合に基づいて、ピンチオフ時間（$t_{p,a}$）に対する解析式を導出した。この式、$t_{p,a} = \frac{1}{\lambda_0 e^{\dot{\epsilon} t_{p,a}} + \dot{\epsilon} \ln(R_i/a_0)}$は、様々な粒子サイズおよび流出率にわたる数値シミュレーションデータと優れた定量的な一致を示している。

意義と主張
本論文は、粘性領域において固体粒子を含むリグメントの破断を制御する局所的な擾乱に関する普遍的なメカニズムを明らかにすると主張している。主な貢献は、リグメント表面が固体粒子に接近すると、粒子の具体的なサイズが最終的な破断ダイナミクスに関係不再となり、システムが普遍的な自己相似領域に入ることを実証した点にある。

著者らは、この普遍性が、粒子誘起擾乱が持続し増幅することを可能にする粘性時間スケール（$t_\mu = \mu R / \sigma$）に起因すると強調している。これは、そのような擾乱が結果を著しく変化させる前に破断が起きる慣性領域（$Oh \ll 1$）とは対照的である。提案された解析的枠組みは、$Oh \sim 1$に対してこの普遍的なスケーリング領域を成功裡に捉え、伸長率および流体物性に依存するが、擾乱が誘起された後は特定の粒子サイズに依存しないピンチオフ時間の予測ツールを提供する。本研究は、対称的で静止した粒子配置に限定されており、これらの効果を分離するための最小構成として機能しており、より複雑な輸送粒子相互作用は将来の研究に委ねられている。