A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries

本論文は、複雑な幾何学形状を有する境界面上での接触線ダイナミクスとヒステリシスを高精度にシミュレートするため、混合セルにおける厳密な質量保存と時間刻み制限の回避、および高度な接触角条件の適用を実現する、新しい3 次元幾何学的かつ保存的な VOF 法を提案するものである。

要約

水滴の「足跡」を完璧に追う：複雑な形をした世界での新しいシミュレーション技術

この論文は、**「水滴が複雑な形をした物体の上をどう動くか」**を、コンピュータの中で非常に正確に再現する新しい方法について書かれたものです。

想像してみてください。雨粒が丸いボールの上を転がったり、インクが細かいメッシュ（網目）の上で広がったりする様子。これらは日常的に見られる現象ですが、コンピュータで「完璧に」シミュレーションするのは、実はとても難しいのです。

この研究は、その難問を解決するための**「3 次元の新しいデジタル・カメラと、水滴の動きを制御する魔法のルール」**を開発しました。

---

1. 従来の問題：なぜ難しいのか？

これまでのシミュレーションでは、水滴の動きを計算する際に、2 つの大きな壁にぶつかっていました。

🧱 壁その 1：小さな「穴」にハマる

コンピュータは空間を小さな箱（セル）の集まりで表現します。水滴が丸いボールや複雑な形をした物体に当たると、その箱の端が「カクカク」した形になってしまいます。

• 昔の方法： 水滴が箱の端に少しだけ入っている場合、その「わずかな部分」を無視したり、計算が狂って水滴が勝手に増えたり減ったりしていました。まるで、**「小さな穴からこぼれた水を、掃除機で吸い取ろうとして逆に部屋中を汚してしまう」**ような状態でした。
• この研究の解決策： 「再分配（リディストリビューション）」という新しいルールを作りました。水滴が箱の端で少し余分になったり足りなくなったりしたら、それを隣の箱に「お裾分け」したり「借りたり」して、全体の水量が絶対に変わらないように調整します。これにより、計算が極端に遅くなるのを防ぎつつ、水滴の量（体積）を完璧に守れるようになりました。

📐 壁その 2：水滴の「足跡」の角度

水滴が固体に触れる部分（接触線）には、決まった角度（接触角）があります。しかし、物体の表面が曲がっていたり、斜めだったりすると、この角度を正しく設定するのが大変でした。

• 昔の方法： 水滴の形を「直線」で近似していました。これは、**「丸い山を直線で描こうとする」**ようなもので、特に角度が急な場合や、曲面の上では大きくズレてしまいました。
• この研究の解決策： **「放物線（パラボラ） fitting（当てはめ）」という新しい手法を使いました。水滴の形を「直線」ではなく、「滑らかな曲線（放物線）」で表現し、さらにその曲線が固体の表面とどう接するかを、事前に計算して補正する「下書き（プレフィッティング）」の工程を追加しました。これにより、「どんなに曲がった山の上でも、水滴が自然な角度で止まる」**ように再現できるようになりました。

---

2. 具体的な成果：どんなことができた？

この新しい技術を使って、いくつかの難しいシミュレーションを行いました。

• 回転する水滴： 水滴が回転する実験で、従来の方法だと水滴の形が歪んでしまいましたが、新しい方法では**「10 回転しても、元の形がほとんど崩れない」**ほど正確でした。
• 丸いボールの上： 水滴が丸いボールの上で広がる様子をシミュレーション。従来の方法だと水滴の輪っかが四角く歪んでいましたが、新しい方法では**「完璧な円」**を描きました。
• メッシュ（網目）の上： 水滴が細かい網の上を広がる様子。これは水滴が「登ったり、降りたり」する複雑な動きをしますが、この方法でも**「滑らかに、かつ安定して」**再現できました。
• 接触角のヒステリシス（記憶効果）： 水滴が「進もうとする時」と「戻ろうとする時」で、表面との接する角度が変わる現象（濡れにくさの違い）も、この方法なら正確に再現できます。

---

3. なぜこれが重要なのか？（日常への応用）

この技術は、単に「水滴のシミュレーションが上手になった」だけではありません。

• インクジェットプリンター： 紙や特殊な素材に、水滴（インク）を正確に落とす技術の向上。
• マイクロ流体デバイス： 小さなチップの中で薬液を運ぶ技術。
• 自己洗浄コーティング： 雨粒が汚れを落としていく仕組みの設計。
• 3D プリント： 溶けた材料が複雑な形の上にどう乗るか、の予測。

これらすべてに、**「複雑な形をした表面での、水滴の動きを正確に予測する力」**が必要です。

---

まとめ：この研究のすごいところは？

この論文は、**「複雑な形をした世界でも、水滴の動きを『質量保存（量が減らない）』と『正確な角度』の両立で、3 次元で完璧に再現する」**という、長年の難問を解決しました。

まるで、**「カクカクしたブロックの世界でも、滑らかな水滴の動きを、一滴たりともこぼさずに、自然な曲線で描き出す魔法のペン」**を手に入れたようなものです。これにより、将来のエンジニアリングや科学技術において、より効率的で美しい「濡れ」の制御が可能になるでしょう。

技術概要

論文の技術的サマリー：複雑な境界におけるヒステリシスを伴う移動接触線動力学のモデル化のための 3 次元鋭い界面かつ保存的な VOF 法

1. 背景と課題 (Problem)

移動接触線（液体/気体/固体の接点）の動力学は、マイクロ流体、インクジェット印刷、薄膜操作など、多くの工学的・科学的プロセスにおいて重要な役割を果たします。従来の研究は主に平坦な基板（規則的な境界）に限定されており、曲がった基板や任意の形状を持つ複雑な幾何学構造（埋め込み境界）上の接触線挙動を正確にシミュレーションすることは困難でした。

特に、複雑な固体境界を扱う際、以下の 2 つの主要な数値的課題が存在します：

1. 混合セル（Mixed Cells）における界面の移流と保存性: 固体境界と流体界面が交差するセル（混合セル）では、従来の幾何学的 VOF（Volume-of-Fluid）法を直接適用すると、局所的な質量保存が破綻し、非物理的な流体の蓄積や欠損が発生します。また、小さなカットセル（固体に切断されたセル）が存在すると、CFL 条件による時間刻みの厳格な制限が生じ、計算効率が著しく低下します。
2. 複雑な表面における接触角の正確な付与: 接触角条件を正しく課すことは、界面の曲率推定と表面張力力の評価に不可欠です。従来の高さ関数（Height Function, HF）法や線形補間法は、平坦な境界では機能しますが、複雑な形状や極端な接触角（鋭角・鈍角）を持つ境界では精度が低下し、接触線の非対称な挙動や誤った曲率評価を招きます。

2. 提案手法 (Methodology)

本研究では、Open-source ソルバー「Basilisk」を基盤とし、幾何学的 VOF 法と埋め込み境界法（Embedded Boundary Method, EBM）を結合した 3 次元数値枠組みを提案しています。

2.1. 混合セルにおける保存的な幾何学的 VOF 移流法

• 修正されたフラックス計算: 混合セルにおける体積分率の移流において、固体境界による体積の欠損を補正する「修正された移流幅（$u_n^*$）」を導入しました。これにより、界面の幾何形状を考慮したフラックス計算を行い、局所的かつ大域的な体積保存を厳密に保証します。
• 3 次元洪水充填アルゴリズム（Flood-fill Algorithm）: 不規則な凸多面体として表現される混合セル内の液体/気体界面を再構築するために、3 次元版の洪水充填アルゴリズムを開発しました。これにより、任意の形状の固体境界内でも PLIC（Piecewise Linear Interface Calculation）による界面再構築が可能になります。
• フラックス再分配戦略（Redistribution Strategy）: 小さなカットセルによる時間刻み制限を回避するため、過剰または不足した体積分率を流下方向の隣接セルへ再分配する手法を導入しました。これにより、CFL 条件の制限を完全に除去しつつ、局所・大域の体積保存を維持し、計算効率を大幅に向上させました。

2.2. 3 次元高さ関数（HF）に基づく接触角付与法

• 放物面フィッティング（Paraboloid Fitting）: 従来の線形フィッティングに代わり、接触線近傍の界面形状をより正確に表現するため、放物面フィッティング手法を 3 次元に拡張しました。
  • 事前フィッティング（Pre-fitting）: 流体領域内の有効な高さ関数データ点を用いて初期の放物面を推定し、安定した補間点を選択します。
  • 制約付き最終フィッティング: 選択された補間点と、接触線位置での接触角条件（$\theta$）を制約として組み込み、非線形方程式系を解くことで、接触角を正確に満たす最終的な放物面を構築します。
• 垂直高さ関数（Vertical HF）: 固体表面が水平面に対してほぼ平行で接触角が極端な場合（$\theta < 45^\circ$ または $\theta > 135^\circ$）、従来の水平 HF が機能しない問題に対処するため、界面の向きに依存しない垂直 HF 法を開発しました。
• ヒステリシスモデル: 接触角ヒステリシス（進行角と後退角の範囲内でのピン留め現象）をモデル化し、体積保存の原則に基づいて接触線のピン留め・脱着を反復的に決定するアルゴリズムを実装しました。

3. 主要な成果 (Key Results)

提案手法の精度、頑健性、保存性を検証するために、一連のベンチマークテストを実施しました。

• 体積保存と移流精度:
  • 同心球モデルと円錐カット付き回転液殻のシミュレーションにおいて、提案した修正移流法と再分配戦略は、混合セル内での局所体積誤差を機械精度レベル（$O(10^{-15})$）にまで低減しました。
  • 従来の手法では混合セルで大きな形状誤差が生じましたが、提案手法は長時間のシミュレーションでも形状を維持し、CFL 制限を回避することで計算コストを削減しました。
• 接触角付与の精度:
  • 平坦および球面上での液滴の自発的拡散シミュレーションにおいて、提案した放物面フィッティング法は、理論解と極めて良く一致する平衡接触線半径を与えました。
  • 従来の線形フィッティング法で見られた接触線の非対称性（矩形や三角形の歪み）が解消され、任意の埋め込み深さや基板の傾きに対しても、接触角条件が正確に課されることが確認されました。
• 複雑な幾何学構造と動的挙動:
  • せん断流中の液滴: 接触角ヒステリシスを考慮したせん断流下での液滴の運動をシミュレーションし、ピン留め・脱着の非対称な挙動を正確に再現しました。
  • 正弦波状およびメッシュ状基板: 粗面や複雑なメッシュ構造上の液滴運動において、表面の山や谷での接触線のピン留め、液滴の捕捉、および曲率勾配による自発的移動（円錐面上の液滴の登攀）を安定してシミュレーションできました。
  • オリフィスへの液滴衝突: 丸い穴と鋭い縁を持つ穴への液滴衝突シミュレーションを行い、実験結果や既存の 2 次元結果と整合する、接触線のピン留めや部分的な閉じ込め現象を 3 次元で正確に捉えました。

4. 意義と貢献 (Significance)

本研究は、複雑な 3 次元幾何学構造上の移動接触線動力学をシミュレーションするための、完全に幾何学的かつ保存的な VOF ベースの枠組みを初めて確立した点に大きな意義があります。

1. 技術的ブレイクスルー: 混合セルにおける厳密な質量保存と、CFL 制限のない時間積分を両立させるアルゴリズムを開発しました。
2. 高精度な境界条件: 任意の形状の固体表面に対して、接触角条件を正確かつ頑健に付与する新しい 3 次元 HF 法を提案し、従来の線形近似の限界を克服しました。
3. 実用性: 液滴の衝突、拡散、自発的移動など、多様な物理現象を複雑な 3 次元形状で高精度にシミュレーション可能であり、基礎研究から実用的な工学応用（マイクロ流体デバイス設計、印刷技術など）まで幅広く貢献するツールとなります。

この研究は、Basilisk ソルバーに実装され、複雑な多相流問題に対する計算流体力学の能力を大幅に向上させるものです。