Instabilities of ring-rivulets: Impact of substrate wettability

本論文は、数値シミュレーションを用いて、環状バンドや径方向の接触角勾配といった基板の濡れ性パターンを調整することで、リング状リブレットの安定性、分裂ダイナミクス、および結果として生じる液滴の形態を精密に制御できることを示している。

要約

テーブルの上に、小さくて完璧な水の輪が置かれているところを想像してみてください。もしテーブルがただの普通の表面であれば、この水の輪は不安定です。それは、今にも弾けそうな風船のようなものですが、弾ける代わりに、自らを修復しようとします。輪の幅がその大きさに対してどの程度かによって、次の2つのどちらかの行動をとります。

1. 粒状に分裂する： もし輪が細く、幅が狭ければ、バラバラの水の滴へと分解され、真珠のネックレスのような形になります。
2. 水たまりへと崩壊する： もし輪が太く、幅が広ければ、内側へと縮まり、すべての水を中央の丸い水たまりへと引き寄せます。

この論文は、その水の輪を制御するための「レシピ本」のようなものです。研究者たちは、コンピューター・シミュレーションを用いて、水の下にあるテーブルの「質感（基板）」を変えると何が起こるかを調査しました。彼らは、テーブルの上に特定のパターンを描くことで、たとえ水が自然にそうしたかったとしても、望み通りの動きを強制できることを見出したのです。

以下に、簡単な比喩を用いてその方法を説明します。

1. 「ベルクロ・トラック」（環状バンド）

水の輪がトラックを走るランナーだと想像してください。普通のトラックでは、ランナーは疲れ果てて中央で止まってしまう（崩壊する）か、あるいはつまずいて散らばってしまう（分裂する）かもしれません。

研究者たちは、水のすぐ下に特別な「ベルクロ（マジックテープ）」の帯（非常に粘着性の高い表面の輪）を設置しました。

• 結果： この粘着性のある輪は、フェンスのような役割を果たします。これは、水が中央の水たまりへと崩壊するのを阻止します。輪がどれほど幅広であっても、ベルクロが水を押さえつけるため、水は必ず粒状に分解されます。
• 制御： この「ベルクロ」を、周囲のテーブルと比較して、より粘着性が高くしたり、逆に低くしたりすることで、形成される液滴の数を正確にコントロールできました。これは、ギターの弦の張力を調整して、特定の音の数を得るようなものです。

2. 「丘と谷」（放射状グラデーション）

次に、彼らは中心から外側に向かって、緩やかな丘や谷のように、テーブルの「粘着性」が徐々に変化するように変更しました。

• 「谷」（下向きの傾斜）： テーブルが中心に向かうほど粘着性が高くなると想像してください。これは滑り台のような役割を果たします。水の輪は中央に向かって強い引きを感じます。たとえ、通常なら粒状に分裂するはずの広い輪であっても、この「滑り台」によって、水は中央へと急いで流れ込み、単一の水たまりへと崩壊させられます。
• 「丘」（上向きの傾斜）： テーブルが中心に向かうほど粘着性が低くなる（あるいは外側に向かって粘着性が高くなる）と想像してください。これは、水が登らなければならない「丘」のようなものです。もし輪が内側へ崩壊しようとすると、抵抗の「壁」にぶつかります。これにより、崩壊が完全に阻止され、たとえ普通のテーブルの上では幅広で不安定な状態であっても、輪は安定して形を保ち続けます。

大きな展望

主な教訓は、これらの液体の輪の形状や振る舞いは、水そのものだけでなく、その下に置かれた表面に大きく影響を受けるということです。

• 均一な表面： 水は自身の自然な本能に従います（分裂するか、崩壊するか）。
• パターン化された表面： 研究者たちは、表面を「交通整理の警察官」として機能するように「プログラミング」することができます。彼らは水に対して、「10個の液滴に分かれなさい」とか、「一つの大きな輪のままでいなさい」、あるいは「中央へ急げ」といった指示を与えることができるのです。

化学的な「質感」を特定のパターンで変えるだけで、彼らは、水の輪がバラバラになるのか、崩壊するのか、あるいは安定したままなのか、そして正確にいくつの液滴を作るのかを、完全に制御することに成功したのです。

技術概要

技術要約：リング状リベットレットの不安定性：基材の濡れ性が与える影響

問題提起
本論文は、部分的に濡れ性の異なる基材上に堆積した、リング状のリベットレット（液状の筋）の脱濡れダイナミクスを扱っている。直線状のリベットレットの崩壊はレイリー・プラトー不安定性によって支配される既知の問題であるが、リング状のリベットレットは、非一様な曲率および内側と外側の接触線の異なる曲率に起因する追加の複雑さを呈している。これまでの理論および数値的研究により、リング状リベットレットの行方（単一の中央ドロップレットへと崩壊するか、あるいは複数のドロップレットへと分裂するか）は、主にアスペクト比（$\psi_0$：リベットレットの幅と半径の比）と基材接触角によって決定されることが確立されている。しかし、設計された基材の濡れ性パターンを通じて、これらの不安定性の経路や結果として生じる脱濡れ形態を能動的に制御できる可能性については、リング形状において未だ十分に探索されていない。本研究では、空間的に変化する接触角を用いることで、リング状リベットレットの安定性、分裂のタイムスケール、および最終的なドロップレット構成を操作できるかどうかを調査している。

手法
著者らは、Lattice Boltzmann法（LBM）を用いたSwalbeソルバーによる薄膜方程式（TFE）に基づく数値シミュレーションを用いている。支配方程式は、膜厚 $h(x,t)$ の進化を以下のように記述する：
$$\partial_t h(x,t) = \nabla \cdot (M_\delta(h) \nabla p)$$
ここで、圧力 $p$ は毛管力 ($-\gamma \nabla^2 h$) と、流体-固体相互作用を考慮し、空間的に変化する平衡接触角 $\theta(x)$ を許容する離散圧 $\Pi(h, \theta)$ を含む。移動度関数 $M_\delta(h)$ は、接触線の特異性を正則化するために有効スリップ長 $\delta$ を組み込んでいる。

初期条件は、原点を中心とするトーラス状のキャップ形状であり、これを平衡形状に緩和させた後、ガウスノイズで摂動を加えたものである。本研究では以下の項目を系統的に変化させている：

1. 一様な基材： 一定の接触角 $\theta_0$ およびアスペクト比 $\psi_0$。
2. 環状バンドパターン： 初期リベットレットの幅内に限定された低接触角 $\theta_a$ の領域と、それを囲む高接触角 $\theta_b$ の背景。
3. 径方向の勾配： 中心に向かって $\theta$ が減少する（内向き）または中心に向かって $\theta$ が増加する（外向き）線形接触角プロファイル。

主要な結果

• 一様な基材：

  • 分裂 vs 崩壊： 臨界アスペクト比 $\psi_0 \approx 0.2$ が、分裂と崩壊を判別する。$\psi_0 < 0.2$ の場合、リングは複数のドロップレット ($n_d > 1$) に分裂する。一方、$\psi_0 > 0.2$ の場合、リングは単一の中央ドロップレットへと崩壊する。
  • ドロップレット数： 分裂レジームにおいて、ドロップレット数 $n_d$ は一般に線形安定性解析（LSA）の予測 $n_d \approx \pi / (2\psi_0)$ に従う（特に小さな接触角の場合）。しかし、大きな接触角 ($\theta_0 > 10^\circ$) では、$n_d$ は $\psi_0$ に対してほぼ独立となり、このプロセスが基材の濡れ性の低下によって駆動されていることを示唆している。
  • タイムスケール： 分裂時間 ($\tau_b$) は、非常に小さな接触角では $\psi_0$ とともに増加するが、大きな角度では $\psi_0$ に依存しなくなる。崩壊時間 ($\tau_c$) は $\tau_c \sim \psi_0^{-2}$ とスケールし、これは現象論的に導出され、シミュレーションによって確認された関係である。

• 環状バンドパターン（低接触角バンド）：

  • 崩壊の抑制： 高接触角の背景の中に低接触角の環状バンドを配置してリベットレットを堆積させると、崩壊モードが事実上排除される。この場合、すべてのテストされたアスペクト比において、リングは分裂を起こす。
  • 形態の制御： バンドと背景の接触角のコントラストが制御パラメータとして機能する。ミスマッチを大きくすると、形成されるドロップレットの数は均質な予測から逸脱し、初期アスペクト比に依存しなくなる。
  • 安定性： このパターンは、一様な基材と比較してリベットレットの破断に対する安定性を高め、分裂時間を長くさせる。

• 径方向の勾配：

  • 内向き勾配（中心に向かって濡れ性が増加）： この構成はリングの後退を加速させる。急峻な勾配は急速な崩壊を誘発し、分裂を完全に阻止する場合がある。緩やかな勾配では、リングが縮小しながら分裂し、最終的に中央で合体するドロップレットへと至る。これにより、ドロップレット形成の過渡的な制御が可能となる。
  • 外向き勾配（中心に向かって濡れ性が減少）： この勾配は、幅広のリング ($\psi_0 > 0.2$) が持つ自然な崩壊傾向に対する安定化力として作用する。これは、リングの主半径を安定させ、穴の閉鎖を防ぐことで、リング形状のままより大きな液量を堆積することを可能にする。

意義および主張
本論文は、基材の濡れ性パターンを調整することで、リング状リベットレットの安定性と脱濡れ形態を精密に制御できると主張している。具体的には：

1. モード選択： （環状バンドによる）崩壊モードの除去、あるいは（外向き勾配による）分裂モードの抑制を選択的に行うことができ、それによってシステムが単一のドロップレットへと進化するか、あるいはドロップレット配列へと進化するかを決定できる。
2. 形態制御： ドロップレットの数と位置は、初期幾何学形状だけでなく、接触角のコントラストや勾配の方向によっても操作可能である。
3. タイムスケール工学： 濡れ性パターンは、後退を加速させる、あるいはリングを崩壊から安定させるなど、脱濡れの特性的なタイムスケールを変化させることができる。

著者らは、これらの知見が、円形のドロップレットパターンの前駆体であるリング状リベットレットの行方を制御するメカニズムを提供するものであると結論付けている。また、提示された数値的枠組みを超えて、時間依存の濡れ性パターンを探索することで、不安定性をさらに操作し、新たな脱濡れ経路を発見できる可能性があることを示唆している。