Droplets sitting on thin elastic sheets: A study with the boundary element method

技術要約：境界要素法を用いた薄い弾性シート上の液滴に関する研究

問題提起
本論文は、液滴と薄い弾性シートとの間の弾性毛管相互作用を調査している。剛体基板とは異なり、柔軟な固体は毛管力の影​​響下で変形し、液滴がシートに沈み込んでレンズ状の構造を形成するという複雑な平衡形状を生じさせる。これまでの研究では、波状のパターン（wrinkling patterns）、デュロタキシス（軟らかい領域への移動）、および応力下での異方的な形状などが探索されてきたが、様々な境界条件や張力プロトコルの下で液滴の平衡形状をシミュレートするための統一された数値フレームワークが必要とされていた。具体的には、著者らは、シートの厚さ、等方的な張力、および一軸方向の引き伸ばしが、液滴の形態、応力分布、および光学特性にどのように影響するかを定量化することを目的としている。

手法
著者らは、液滴と弾性シートの結合ダイナミクスをシミュレートするために、以前開発した境界要素法（BEM）を拡張した。この手法の核心は以下の通りである：

1. 動的方程式： システムは、液滴表面の速度場（$\dot{R}_0$）と弾性シートの速度場（$\dot{R}_{sh}$）が結合した過減衰ダイナミクスを用いてモデル化される。支配方程式は、線形システム $G\dot{R} = -\nabla_R F$ である。ここで、$G$ は摩擦行列、$F$ は全自由エネルギーを表す。
2. 自由エネルギー定式化：
   • 液滴 ($F_0$)： 液体のガス表面張力と、ヤングの接触角によって制御される基板・液体界面に対する仕事を含む。ラインテンション項も含まれるが、検討されている弾性シートにおいては無視できるため、ゼロとして設定されている。
   • シート ($F_{sh}$)： 面内のせん断と面積膨張（超弾性）のためのスカラック（Skalak）モデル、および面外曲げのためのヘルフィッツ（Helfrich）モデルを用いてモデル化される。弾性係数（せん断 $\kappa_S$、面積 $\kappa_A$、および曲げ $\kappa_b$）は、バルクの材料特性（ヤング率 $E$ およびポアソン比 $\nu$）とシートの厚さ $d$ から導出される。
   • 制約 ($F_{cons}$)： 液滴の一定の体積を維持し、界面における液滴メッシュとシートメッシュの幾何学的連結性を確保するために、ラグランジュ乗数を用いて強制される。
3. 摩擦行列： 行列 $G$ は、粘性せん断応力（オーセンドテンソルから導出）、接触線摩擦（Moffattの結果）、およびナビエの滑り条件によって制御される液体・基板摩擦を組み込んでいる。
4. 数値実装： システムは gmsh によって生成された三角形メッシュを用いて離散化される。著者らは、液滴とシートの異なるメッシュ間で力と速度をマッピングするために、重み付き補間スキーム（Surrogates.jl を使用）を採用している。シミュレーションは、接触線半径 $R_0$ と、Cox-Voinov則から導出される特性時間スケールを用いて無次元化されている。

主要な結果
本研究では、拡張されたBEMを以下の3つの特定のプロトコルに適用している：

1. 固定された弾性シート：

   • 液滴の平衡形状はレンズ状であり、シートの厚さが増すにつれて接触線の平面に対してますます非対称になる。
   • 厚さ依存性： 厚さ ($d$) が増加するにつれて、液滴の下部 ($h_l$) は減少し、上部 ($h_u$) は増加するが、全高は比較的一定に保たれる。
   • 接触角： 接触線における微視的な接触角はヤングの角と一致しており、モデルの妥当性を検証している。しかし、マクロな接触角（レンズ全体の形状によって決定される）はヤングの角とは大きく異なり、シートの厚さに依存する。
   • デカップリング（分離）： 液滴の下半分の形状は主にシートの弾性によって決定され、平衡接触角 ($\theta_{eq}$) にはほとんど依存しない。一方で、$\theta_{eq}$ は主に液滴の上半分に影響を与える。

2. 等方的に引き伸ばされたシート：

   • 等方的な張力を加えると、液滴が持ち上がり、シートの変形が抑制され、高い伸長時（$\epsilon \approx 0.1$）には再び平面に戻る。
   • 応力プロファイル： 著者らは、方位角方向（$T_{\phi\phi}$）および半径方向（$T_{rr}$）の弾性応力の径方向プロファイルを定量化している。応力は伸長とともに増加し、接触線は応力曲線の変曲点の近くに位置する。
   • 光学的なチューニング： レンズ状の液滴は液体光学レンズとして機能する。レンズメーカーの公式を介して計算される焦点距離は、加えられた張力を変化させることで、約15%可逆的に調整可能である。この範囲は、現代のスマートフォンのレンズのチューニング範囲に匹らである。

3. 一軸方向に引き伸ばされたシート：

   • シートを一方向に引き伸ばすと、液滴は長楕円形になる。伸長の程度は、側方の境界条件（固定、制御された内側への移動、または自由）に依存する。
   • 形態： 側方の端が自由な場合、最も顕著な伸長が生じ、大きな垂直変位が発生し、シートに上方への折り畳みやディンプル（窪み）が生じる。接触線は非平面となり、側方方向には高く、引き伸ばし方向には低くなる。
   • 異方性： 接触線のフットプリントのアスペクト比と、マクロな接触角の異方性が定量化されている。この伸長は、液滴が引き伸ばし方向と側方方向で異なる実効的なシート剛性を経験することに起因する。

意義と主張
本論文は、拡張されたBEMが、流体ダイナミクスと弾性変形の相互作用を捉えることに成功した、弾性毛管システムをシミュレートするための堅牢な数値ツールであることを主張している。主な貢献は以下の通りである：

• 方法論的貢献： 液滴表面のダイナミクスと、面内および曲げエネルギーを含む弾性シートの変形を結合させた、定式化されたBEM。
• 物理的洞察： シートの厚さが液滴の非対称性を制御する重要なパラメータであること、および、マクロな変形にもかかわらず微視的な接触角がヤングの法則と一致していることを実証した。
• 応用： 弾性シートに張力を加えることで、液体の焦点距離を動的かつ可逆的に調整できることを明示的に示した。これは、マイクロ・オプトフルイディック・デバイスにおける光学特性を制御するための潜在的なメカニズムを提供する。
• 今後の展望： 著者らは、この手法が波状パターンの研究、張力勾配におけるデュロタキシス、および複数の液滴間の弾性毛管相互作用の研究への道を開くものであると示唆している。特に、一軸方向の張力は、異方的な相互作用ポテンシャルや複雑なダイナミクスをもたらす可能性がある。