[論文レビュー] Simple and efficient representations for the fundamental solutions of Stokes flow in a half-space

Zydrunas Gimbutas, Leslie Greengard+1 more|arXiv (Cornell University)|Apr 1, 2015

Lattice Boltzmann Simulation Studies参考文献 30被引用数 39

ひとこと要約

本稿では、スカラーハーモニック関数に基づく単一のパプコビッチ＝ニューバー補正項と、反射されたストークスレットを組み合わせることで、半空間内におけるストークス流れの基本解を簡略化・効率化した表現を提示する。この手法は、1つのイメージ力と1つのハーモニック補正項のみを用いて、正確なスリップなし境界条件を達成し、古典的なブレイク型の複雑な公式とは対照的に、より直感的で計算に適した代替手法を提供する。

ABSTRACT

We derive new formulas for the fundamental solutions of slow, viscous flow, governed by the Stokes equations, in a half-space. They are simpler than the classical representations obtained by Blake and collaborators, and can be efficiently implemented using existing fast solvers libraries. We show, for example, that the velocity field induced by a Stokeslet can be annihilated on the boundary (to establish a zero slip condition) using a single reflected Stokeslet combined with a single Papkovich-Neuber potential that involves only a scalar harmonic function. The new representation has a physically intuitive interpretation.

研究の動機と目的

• 古典的なブレイク型の公式とは対照的に、半空間内ストークス流れの基本解をより簡潔かつ効率的に表現すること。
• 複数のディポールおよびクアドリポール項を含む従来のイメージ法による計算複雑性を軽減すること。
• 最小限の構成要素で、スリップなし境界条件を物理的に直感的かつ明確に実装するフレームワークを提供すること。
• 自由空間の基本解およびハーモニックポテンシャルに対する既存の高速ソルバーを活用して、効率的な実装を可能とすること。
• ストレステット、ロットレット、ストークスダブルットに対しても、統一的かつコンactな表現を用いてこの手法を拡張すること。

提案手法

• パプコビッチ＝ニューバー表現を用い、ハーモニック関数 φ(x) をストークス解に写像する。具体的には、u(x) = x₃∇φ(x) − [0,0,φ(x)]ᵀ および p(x) = 2∂φ/∂x₃ を用いる。
• 境界 x₃ = 0 上の接線方向速度成分をキャンセルするため、画像点 yI = (y₁,y₂,−y₃) に反射されたストークスレットを用いる。
• 境界上の法線方向速度成分をキャンセルするため、yI における源およびディポールから導かれるハーモニック関数 φ(x) を用いた、単一のパプコビッチ＝ニューバー補正項を適用する。
• ストレステットおよびダブルットの場合、補正ポテンシャル φ は、yI におけるディポールおよびクアドリポールのハーモニック源から構成され、力および方向ベクトルから導かれる係数を含む。
• 自由空間のグリーン関数（ストークスレット、ディポール、クアドリポール）およびハーモニックポテンシャルの評価に既存の高速アルゴリズムを活用し、FMMベースのソルバーへの統合を可能にする。
• 2次元ストークス流れおよび線形弾性に対しても同様の定式化を導出し、本手法の広範な適用可能性を示す。

実験結果

リサーチクエスチョン

• RQ1 ブレイクの解が複雑な複数のディポール構造を必要とするのに対し、半空間内ストークスレットに対しては、より単純なイメージ表現が得られるか？
• RQ2 複数の高次特異項を用いるのではなく、単一のハーモニックポテンシャル補正項のみでスリップなし境界条件を満たすことは可能か？
• RQ3 パプコビッチ＝ニューバー表現は、半空間内におけるすべての基本解（ストークスレット、ストレステット、ロットレット、ダブルット）に対して、物理的に直感的かつ計算に適した定式化を可能にするか？
• RQ4 既存の自由空間基本解用高速ソルバーを活用して、新規定式化を効率的に実装できるか？
• RQ5 新規定式化は、ブレイクのオリジナルの公式と比較して、複雑さおよび計算コストの点で優位性を示すか？

主な発見

• 提案手法は、ブレイクの複雑なイメージ構造（複数のディポールおよびクアドリポール項を含む）を、スカラーハーモニック関数に基づく単一のパプコビッチ＝ニューバー補正項に置き換える。
• ストークスレットの場合、補正ポテンシャル φ(x) は、画像点 yI における源およびディポールの線形結合から構成され、x₃ = 0 で u₃ = −φ(x) となる。これにより法線方向速度がキャンセルされる。
• 本手法は、1つの反射ストークスレットと1つのパプコビッチ＝ニューバー補正項のみを用いて、境界 x₃ = 0 における正確なスリップなし条件（u = 0）を満たす。
• ストレステットの補正には、yI におけるディポールおよびクアドリポールからなるハーモニックポテンシャル φT(x) を用い、係数は νI·gI および y₃ に依存する。
• ロットレットの補正には、係数が −2ν₃^I および 2g₃^I の2つのディポールを用い、ストークスダブルットの補正は、対称および反対称成分を1つの φD(x) に統合し、3つの項を含む。
• 本定式化は、2次元ストークス流れおよび線形弾性へ自然に拡張可能であり、2次元ハーモニックポテンシャルを用いた同様のパプコビッチ＝ニューバー補正が適用可能である。