3次元MPS法のための新しい自由表面同定法

動く水面を追跡することが重要な理由

堤防にぶつかる波や宇宙船燃料タンク内のスロッシングまで、液体の挙動とそれが及ぼす力を予測するために、エンジニアはますますコンピュータシミュレーションに依存しています。これらのシミュレーションの多くは、水を剛体格子ではなく仮想粒子の群れで表現します。しかし微妙な問題があります。コンピュータは常にどの粒子が液体の外側の表面にあり、周囲の空気の圧力がかかる場所なのかを判断し続けなければなりません。その判定がわずかにでも誤ると、壁や構造物に作用する予測圧力がノイズを含んだり誤解を招いたりします。本研究は、三次元シミュレーションでこれらの表面粒子を見つける新しい方法を紹介し、仮想の水をより実物に近づけます。

水と空気の境界を引く別の方法

移動粒子半陰的法（MPS）では、水は無数の小さな粒子としてモデル化され、粒子同士が移動・相互作用します。従来の手法では、粒子が自由表面にあるかどうかを単に近傍粒子の数で判定することが多いです。近傍が少ないと表面付近にいるとみなされます。しかしこの経験則は、乱流や鋭い角付近で粒子が不均一に配置される場合に失敗することがあります。著者らは代わりに、近傍粒子が空間にどのように配置されているかに着目し、相対位置発散（relative position divergence, RPD）と呼ばれる量を使います。これはある粒子の周囲に近傍がどれだけ「広がっている」かを測るもので、近傍が一方に偏っていると発散の値が変化し、それが表面である可能性を示します。

表面判定に方向性を加える

改良された発散量は理論的期待によく一致しますが、著者らは三次元ではこの値に単純な閾値を設けるだけでは、特に鋭い縁や空洞内部で内部粒子と表面粒子を完全に分離できないことを示します。これを解決するために、周囲の粒子から局所的な表面方向、すなわち法線ベクトルを推定する追加ステップを導入します。この方向を用いて、各候補粒子から外向きに向かう円錐状の領域を定義します。その円錐内に一定距離以内の他の粒子が存在しなければ、その粒子は自由表面上にあると判定します。この組み合わせアプローチ（RPD+NV）は、近傍の配置の情報と空間の空きの位置の両方を利用し、液体と空気の真の界面をより確実に描きます。

複雑な形状と静水条件での検証

研究者らはまず静止した純粋に幾何学的な形状で手法を試験します：ドームを載せた立方体、ボウル状の凹みをもつ立方体、内部に空洞を持つS字管などです。これらは膨らむ面やくぼむ面の両方でアルゴリズムに挑戦します。検出された表面粒子を既知の真の表面と比較することで、新しい手法が凸部や凹部を正しく捉え、隠れた空洞を正確に解像することを示しています。次に単純な静水容器に適用し、粒子間隔とタイムステップを合理的に選べば、底部での予測圧力が理論値と非常によく一致することを検証します。こうした静穏な条件下では、新しい表面検出器は従来法より圧力の振動が少なく、計算コストはほぼ変わりません。

波、衝撃、スロッシングの再現

次にチームは非常に動的な状況で手法を試します。三次元ダムブレークのシミュレーションでは、水柱が崩壊して壁に衝突し、急変する自由表面、飛沫、空気ポケットが生じます。近傍数や発散のみを基にした従来基準と比べて、RPD+NV法はより連続的で細部のある表面を追跡し、より滑らかで現実的な圧力場を生み出します。重要なのは、壁面での最大衝撃圧力が実験室測定にはるかに近づく点です。減衰したスロッシングのケースでは、タンク内の波が徐々に収まる過程で新手法が壁に作用する力の最も滑らかな変化をもたらし、数値的ノイズが低減されていることを示します。最後に、左右に揺れる環状タンク内の燃料スロッシングを模した試験では、シミュレーションで得られる表面形状と圧力ピークがハイスピードカメラ画像や圧力センサデータと良く一致し、表面での不自然な圧力スパイクも生じませんでした。

実世界のシミュレーションへの意義

計算流体力学の外の読者に向けた実務的メッセージは、粒子ベースのモデルにおいてコンピュータが水と空気の境界をどのように決めるかが改善されたということです。近傍が粒子をどのように取り囲んでいるかを洗練して測る指標と、方向性を持つ「空の円錐」テストを組み合わせることで、新しいRPD+NV法は表面粒子の誤判定を大幅に減らします。その結果、静水から激しい衝撃、複雑な三次元スロッシングまで広範な流れで、よりクリーンで安定した圧力予測が得られます。計算コストの増加が小さいため、沿岸構造物、船舶運動、貯蔵タンクなど、動く水面の正確な捕捉が重要な工学シミュレーションにとって実用的な改良手段を提供します。

引用: Geng, C., Wang, Wh., Heng, My. et al. A new free surface identification method for 3D MPS method. Sci Rep 16, 13829 (2026). https://doi.org/10.1038/s41598-026-44218-9

キーワード: 粒子ベースの流体シミュレーション, 自由表面検出, 数値的波衝撃, タンク内のスロッシング, 移動粒子法