Color-gradient lattice Boltzmann modeling of wetting boundary condition on curved solid boundaries

本論文は、ゴーストノード上のオーダーパラメータを更新することにより、カラー勾配格子ボルツマン法における曲面固体表面に対する濡れ境界条件を導入するものであり、この手法は、大きな密度および粘性コントラストを効果的に処理し、偽の流速（スプリアス電流）を最小限に抑えつつ、静的および動的な接触線の挙動の両方を正確に再現することをGPUハードウェア上で検証された。

要約

あなたは、水滴が曲面（葉っぱに落ちる雨粒や、曲面ガラスを滑り落ちる泡など）に当たったときにどのように振る舞うかをシミュレーションしようとしていると想像してください。これをコンピュータで行うために、科学者たちは**格子ボルツマン法（Lattice Boltzmann Method）**と呼ばれる手法を使用します。この手法は、コンピュータの画面を覆う巨大で目に見えない、小さなタイルのグリッドのようなものだと考えてください。各タイルには、流体の情報が少しずつ保持されており、コンピュータはこれらのタイルをステップごとに更新して、流体がどのように動くかを観察します。

難しい部分は、境界条件、特に流体が固体の壁に触れたときにどのように振る舞うかという点です。現実の世界では、水は壁に当たってピタッと止まるわけではありません。表面が濡れているか（きれいなガラスのように）、あるいは乾いているか（ワックスを塗った車の表面のように）によって、特定の角度（接触角と呼ばれます）を形成します。

問題点：機械の中に潜む「ゴースト」

コンピュータのシミュレーションにおいて、固体の壁は滑らかな線ではなく、正方形のグリッドタイルでできているため、ギザギザしています。計算を成立させるためには、コンピュータは、そこには流体が存在しないにもかかわらず、固体の壁の内部で流体が何をしているのかを知る必要があります。これらの壁の中にある想像上のスポットは、**「ゴーストノード（幽霊ノード）」**と呼ばれます。

これらのゴーストノードに何をさせるかを伝える従来の方法には、いくつかの欠点がありました。

• 流体が何らかの力なしに勝手に動いているように見える、偽の「ゴースト電流（スプリアス速度）」を作り出すことがありました。
• 曲面の扱いが苦手で、多くの場合、平らな壁のために設計されたかのように振る舞いました。
• 中立的な角度（水が広がったり、玉になったりもしない状態）を扱うために、特別な複雑な数学を必要とすることがありました。

解決策：ゴーストのための新しいルール

著者らは、これらのゴーストノードのための、よりシンプルで新しいルールを導入しました。

比喩： 流体には「気分（色で表現され、ガスは0、液体は1）」があると想像してください。現実の世界では、表面を横切るときに、この気分はガスから液体へと滑らかに変化します。

• 従来の方法： 壁の向きに立っている人の気分を、適当な推測で叫んで当てるようなものでした。
• 新しい方法： もし壁のすぐ外側（流体側）にいる人の「気分」を知っていれば、その滑らかな「気分」の曲線を、壁を通してゴーストノードまで数学的に延長できることに著者らは気づきました。彼らは単にこう問いかけます。「もし流体がここで45度の角度を作りたいのであれば、ゴーストノードの気分はどうあるべきか？」

この新しいルールは、継ぎ目のない架け橋のようなものです。それは、液滴の自然な形状を、固体の壁の直前、そしてわずかに壁の内側まで延長し、液滴が作る角度が科学者が要求した通りのものになるようにします。

検証内容

彼らの新しいルールが機能することを証明するために、彼らは非常に強力なコンピュータチップ（NVIDIA A100 GPU）を使用して、いくつかのシミュレーションを実行しました。

1. 静止した液滴： 平らなプレートと円柱の上に水滴を置き、液滴が指定した通りの角度で落ち着くかどうかを確認しました。
   • 結果： 彼らの新しいルールは、特に角度が非常に鋭い場合（玉になった液滴）や非常に平らな場合（広がった液滴）において、従来の最高の手法よりも正確でした。
2. 浮遊する粒子： 油と水の境界に浮かぶ円柱をシミュレートしました。
   • 結果： 彼らの手法は、水線の位置を以前よりも正確に算出しました。
3. 落下する液滴： 液滴が落下して円柱に当たり、跳ねたり広がったりする様子をシミュレートしました。
   • 結果： 液滴は現実的に振る舞い、新しいルールによって流体に奇妙な偽の動きが発生することもありませんでした。

主なポイント

• 精度： 新しい手法は曲面をはるかにうまく扱い、壁が平らであっても丸であっても、流体の角度を正しく保ちます。
• 安定性： 「偽のノイズ（スプリアス電流）」をほとんど生成しないため、シミュレーション内の流体はより自然に見えます。
• 簡潔さ： 接触角がちょうど90度（中立）の場合に、従来のメソッドでは厄介だった特別な複雑な数学を回避できます。
• スピード： 最新のコンピュータチップ（GPU）と特定のプログラミングスタイルを使用することで、シミュレーションを非常に高速に実行できるようにしました。彼らは、少し精度の低い数値形式（単精度）を使用することで、ほとんどのテストの結果を損なうことなく、コンピュータの実行速度を2倍にできることを見出しました。

要約すると、著者らは、コンピュータのシミュレーションが液体と固体の境界をどのように扱うかについての、より優れた「ルールブック」を構築しました。これにより、デジタル上の液滴が、曲面であっても、現実の液滴のように振る舞うようにしたのです。

技術概要

技術要約：曲面固体境界における濡れ境界条件のカラー勾配格子ボルツマンモデリング

問題提起
濡れの現象を正確にモデリングすることは、多孔質媒体内の流れ、マイクロフルイディクス、コロイド膜の乾燥といった技術的に重要なプロセスをシミュレーションする上で極めて重要である。拡散界面を利用した格子ボルツマン法（LBM）フレームワークにおいて、曲面固体境界上に濡れ境界条件を実装することには、大きな課題が伴う。既存の手法は、人工的な質量源または吸収源の導入、界面付近での偽の（寄生）速度の発生、平面幾何学への制限、およびノンスリップ条件と濡れ条件を同時に課す際の不整合といった制限にしばしば直面する。さらに、多くの手法は中立接触角（$\theta = 90^\circ$）に対して特別な処理を必要としたり、誤差を導入する線形補外に依存したりしている。

手法
著者らは、二次元カラー勾配LBMフレームワーク内で曲面固体境界に特化して設計された濡れ境界条件を提案している。この手法の核心は、固体相の内部に位置する「ゴーストノード」における秩序パラメータ（カラー場、$\phi$）の更新にある。

• ゴーストノード更新則: 定数値を規定したり線形補外を行ったりする代わりに、本スキームは平衡カラープロファイルを固体相へと拡張する。ゴーストノード（$x_G$）におけるカラー場の値は、固体法線方向に距離 $\Delta x$ だけ離れた流体領域内の対応点（$x_F$）におけるカラー場に基づいて算出される。
• 数学的定式化: 秩序パラメータの分布が平衡状態に近いと仮定すると、固体法線方向における秩序パラメータの微分は、平衡プロファイルと規定された接触角 $\theta$ から導出される。この関係を積分することで、以下の非線形更新則が得られる：
  $$\phi_G = \frac{\phi_F}{\phi_F + (1 - \phi_F) \exp(\varepsilon)}$$
  ここで、$\varepsilon = -4\Delta x \cos \theta / W$ であり、$W$ は界面幅である。
• 主な特徴: この定式化により、ゴーストノードのカラー場が0から1の間に限定されることが保証され、人工的な質量移動が防止される。決定的なことに、このルールは、以前のスキームとは異なり、$\theta = 90^\circ$（$\varepsilon = 0$）の場合に特別な処理を必要としない。本手法は、カーネル融合と自動微分を利用したNVIDIA A100 GPU上で動作するJAXベースのフレームワークを用いて実装されている。

主要な貢献と結果
提案されたスキームは、いくつかのテストケースにおいて、解析解およびFakhariら[20]のベンチマークスキームに対して検証された。

1. 平板上の静止液滴: 本スキームは、$30^\circ$から$150^\circ$の範囲において、Fakhariらの手法と比較して、測定された接触角の誤差が一貫して低い値を示した。また、ベンチマークスキームと比較して、より小さな寄生流（$5.74 \times 10^{-7}$から$7 \times 10^{-7}$格子単位）を生成した。
2. 円柱上の液滴: 曲面円柱の上に載った液滴において、提案手法は液滴頂点の高さ（$h_{max}$）の予測において、特に極端な接触角（$30^\circ$および$150^\circ$）において精度が向上した（ベンチマークの誤差約$0.76\%$に対し、約$0.35\%$に減少）。
3. 液ガス界面に固定された粒子: 界面に捕捉された円柱の周期配列のシミュレーションにおいて、提案スキームは、極端な角度における界面の垂直変位の誤差を、ベンチマークの約$14\%$に対し、約$9\%$に低減させた。
4. 動的な液滴衝突: 本スキームは、高レイノルズ数および高ボンド数における円柱への液滴衝突を正常にシミュレートし、親水性表面での液滴の広がりや疎水性表面での分裂といった複雑な動的挙動を捉えた。結果は、膜厚に関する解析的なスケーリング則とよく一致した。
5. 計算性能: GPU実装により、単精度（FP32）計算は、精度を維持しつつ（静止および中程度の動的ケースにおいて）、倍精度（FP64）計算よりも少なくとも2倍高速であることが示された。しかし、高レイノルズ数／ウェーバー数での液滴衝突テストでは、FP32は不安定になったため、このような高エネルギーのシナリオにはFP64が必要となった。

意義と主張
本論文は、提案された濡れ境界条件が、曲面固体境界における非混和性二相流をシミュレートするための堅牢かつ正確な代替案であることを主張している。特定のゴーストノード更新則を通じて平衡プロファイルを固体相へと拡張することにより、本手法は以下の特性を持つ：

• 大きな密度および粘度コントラストを扱うモデルの能力を維持する。
• 比較的小さな寄生流を生成する（一般的な擬似ポテンシャルモデルよりも数桁小さい）。
• 他の手法に見られる中立接触角に関連する数値的な特異性を回避する。
• 曲面幾何学における静的および動的な接触線について、既存のスキーム（具体的にはFakhariら）よりも優れた精度を示す。

著者らは、現在の実装はシングルGPU実行および二次元に限定されているものの、基礎となる物理学はフェーズフィールドモデルや自由エネルギーモデルに直接適用可能であり、三次元へ拡張可能であると述べている。本研究は、特に接触線毛管力計算やヒステリシス・スキームの進展と組み合わせた際の、液ガス界面における固体粒子のシミュレーションに向けた基盤を提供するものである。