[論文レビュー] An Immersed Boundary Method with Divergence-Free Velocity Interpolation

Yuanxun Bao, Aleksandar Donev+3 more|arXiv (Cornell University)|Jan 25, 2017

Lattice Boltzmann Simulation Studies参考文献 31被引用数 4

ひとこと要約

本稿では、流体-構造連成シミュレーションにおける体積保存性を著しく向上させるために、発散なし速度補間と力拡散を強制する、新しい浸漬め境界（IB）法を提案する。補間されたラグランジュ速度場がC¹かつ発散なしであることを保証し、随伴に基づく力拡散を用いることで、特に圧力勾配が大きい薄い境界においても、不自然な体積変化を大幅に低減する。計算コストの増加はわずかにとどまる。

ABSTRACT

The Immersed Boundary (IB) method is a mathematical framework for constructing robust numerical methods to study fluid-structure interaction in problems involving an elastic structure immersed in a viscous fluid. The IB formulation uses an Eulerian representation of the fluid and a Lagrangian representation of the structure. The Lagrangian and Eulerian frames are coupled by integral transforms with delta function kernels. The discretized IB equations use approximations to these transforms with regularized delta function kernels to interpolate the fluid velocity to the structure, and to spread structural forces to the fluid. It is well-known that the conventional IB method can suffer from poor volume conservation since the interpolated Lagrangian velocity field is not generally divergence-free, and so this can cause spurious volume changes. In practice, the lack of volume conservation is especially pronounced for cases where there are large pressure differences across thin structural boundaries. The aim of this paper is to greatly reduce the volume error of the IB method by introducing velocity-interpolation and force-spreading schemes with the property that the interpolated velocity field in which the structure moves is at least $\mathscr{C}^1$ and satisfies a continuous divergence-free condition, and the force-spreading operator is the adjoint of the velocity-interpolation operator. We confirm through numerical experiments in two and three dimensions that our new IB method equipped with divergence-free velocity interpolation and force spreading is able to achieve substantial improvement in volume conservation compared to other existing IB methods, at the expense of a modest increase in the computational cost.

研究の動機と目的

• 従来の浸漬め境界（IB）法において、特に圧力差が大きい薄い構造境界で顕著な体積保存性の悪さが生じるという、長年の問題に対処すること。
• 補間されたラグランジュ速度場が少なくともC¹であり、連続的な発散なし条件を満たすようにする速度補間スキームの開発。
• 速度補間演算子の随伴として構築された力拡散演算子の構築により、一貫性と保存性の性質を保証すること。
• 数値的安定性を損なわず、過度な計算コストを増加させることなく、体積保存性を著しく向上させること。

提案手法

• 速度補間にC¹の正則化されたデルタ関数カーネルを用いることで、流体速度をラグランジュ構造点に滑らかで発散なしに補間する。
• 速度補間演算子は、連続的意味で補間された速度場が発散なしであるように設計されており、体積の整合性を保つ。
• 力拡散演算子は、速度補間演算子の随伴として構築され、ラグランジュ力とオイラー流体力の間の一貫性を保証する。
• 変分式を用いて補間および拡散演算子を導出することで、数学的整合性と改善された保存性特性を確保する。
• 2次元および3次元の両方で実装され、離散的レベルでも発散なしの性質を維持する数値スキームが採用されている。
• 体積保存性と精度を評価するために、ベンチマークとしての流体-構造連成問題を用いて検証が行われた。

実験結果

リサーチクエスチョン

• RQ1 発散なし速度補間スキームは、特に薄い構造・高圧力勾配の状況において、浸漬め境界法における体積誤差を顕著に低減できるか？
• RQ2 補間カーネルにC¹正則性を課すことで、IBシミュレーションにおける流体体積保存性にどのような影響を与えるか？
• RQ3 随伴に基づく力拡散演算子の使用は、IB法の安定性および精度にどのような影響を与えるか？
• RQ4 提案手法は、標準的なデルタ関数カーネルを用いた従来のIB法と比較して、どの程度体積保存性を向上させるか？

主な発見

• 提案手法は、特に圧力差が大きい薄い構造を有するケースにおいて、従来のIB法と比較して顕著な体積保存性の向上を達成した。
• 補間された速度場は発散なしかつC¹であるため、標準的なIB形式でよく見られる不自然な体積変化を防止する。
• 随伴力拡散演算子の使用により、ラグランジュ力とオイラー流体力の間の一貫性が保証され、保存性特性が向上する。
• 2次元および3次元における数値実験により、本手法が高い精度を維持しながら、体積誤差を顕著に低減することが確認された。
• 計算コストの増加はわずかであり、実世界の流体-構造連成シミュレーションにおいて実用的である。
• 本手法は、さまざまなベンチマーク問題においても強固な性能を示し、有効性と一般性を確認した。