[論文レビュー] Making sessile drops easier

Amir H. Fatollahi, Maryam Hajirahimi|arXiv (Cornell University)|Apr 22, 2013

Surface Modification and Superhydrophobicity参考文献 10被引用数 3

ひとこと要約

本稿では、ヤング＝ラプラス方程式から直接関係式を導出し、従来の多パラメータ手法の複雑さを低減することで、軸対称なドロップの平衡形状を数値的に解くための一パラメータのシューティング法を提示する。この手法は、実験データおよびADSA-D手法と比較して、Bond数15まで高い精度を達成しており、実用的な用途を想定したMathematica実装が提供されている。

ABSTRACT

Using an identity, directly derived from the Young-Laplace equation, the problem of the equilibrium shape of an axisymmetric sessile drop is reduced to a one-parameter shooting method problem. Based on the method the numerical solutions for drops with Bond number up to 15 are plotted. The agreement between the method and the ADSA-D method as well as the experimental data is tested. A Mathematica code based on the method is presented.

研究の動機と目的

• パラメータ空間の削減により、軸対称なドロップの平衡形状の数値解法を簡素化すること。
• ヤング＝ラプラス方程式から導出された新しい恒等式を導出し、一パラメータのシューティングアプローチを可能にすること。
• 広範なBond数（最大15）の範囲で、既存の数値的手法（ADSA-D）および実験データと比較して、本手法の妥当性を検証すること。
• 再現可能でアクセス可能なドロッププロファイルの計算を可能にする、実用的でオープンソースのMathematicaコードを提供すること。

提案手法

• ヤング＝ラプラス方程式から、曲率と重力圧力の関係を示す重要な恒等式を導出し、境界値問題を単一のシューティングパラメータに還元可能にする。
• 頂点の曲率を唯一の調整可能なパラメータとして用いることで、軸対称なドロップ問題を一パラメータのシューティング法に変換する。
• MathematicaのNDSolveを用いて、頂点および接触線における境界条件を満たすように、得られた常微分方程式を数値的に解く。
• 物理的整合性を保証するため、頂点における傾きの絶対値が0.1未満である条件を収束基準として適用する。
• ドロッププロファイルを角度θのパラメトリック表現r(θ)として記述することで、等周点および頂点の位置を正確に追跡可能にする。
• プロファイルにおける最大半径を求めるために、根の探索手順を用いて等周半径および高さを特定する。

実験結果

リサーチクエスチョン

• RQ1 複数のパラメータではなく、1つのシューティングパラメータのみを用いて、ドロップの平衡形状を正確に計算できるか？
• RQ2 ヤング＝ラプラス方程式から導出された恒等式は、ドロップ形状問題の数値的取り扱いをどのように簡素化するか？
• RQ3 さまざまなBond数における実験データおよびADSA-D手法と比較して、一パラメータシューティング法の精度はどの程度か？
• RQ4 本手法は、鋭角および鈎角の接触角を有するドロップ、特に高Bond数のケースに対しても、どの程度耐性を示すか？

主な発見

• 一パラメータシューティング法は、実験データと非常に良好に一致しており、複数の液体-表面系において、接触半径および高さの相対誤差が5%未満である。
• 水をカーボン鋼上に滴下した場合（θ = 72°, V = 6.75×10⁻³ cm³）、実験値と比較してρ₀ = 0.1741 cm（0.4%誤差）およびh = 0.1199 cm（4.4%誤差）を計算した。
• 水銀をガラス上に滴下した場合（θ = 136°, V = 1.27×10⁻³ cm³）、h = 0.1119 cmを予測し、1.5%の誤差を示しており、高接触角ドロップに対しても安定性が確認された。
• 本手法は、重力と表面張力が釣り合う小・大規模ドロップをカバーするBond数15まで、ドロッププロファイルを正確に計算可能である。
• 提供されたMathematicaコードにより、体積、表面張力、密度、接触角といったパラメータを調整可能な状態で、直接的にドロッププロファイルを計算可能である。
• 多パラメータ最適化の必要性を排除することで、計算の複雑さを低減し、収束性を向上させ、従来のシューティング法よりも効率的である。