[論文レビュー] Particle-based simulations of steady-state mass transport at high Péclet numbers

Thomas Müller, Paolo Arosio+6 more|arXiv (Cornell University)|Oct 17, 2015

Groundwater flow and contamination studies被引用数 3

ひとこと要約

本稿では、初期濃度ではなくフラックス分布に従って粒子を初期化することにより、高Péclet数における定常状態の物質移動を粒子ベースでシミュレートする手法を提示する。これらのフラックス重み付き粒子をモンテカルロ経路に従って進める際、事前に定義された領域内を通過する粒子を検出することで、下流の濃度プロファイルを正確に予測する。実験的なマイクロフルイディクスデータ（蛍光コロイドを用いたもの）と非常に良好な一致を示し、粒子半径25 nmを用いた検証によって本手法の有効性が裏付けられている。

ABSTRACT

Conventional approaches for simulating steady-state distributions of particles under diffusive and advective transport at high Péclet numbers involve solving the diffusion and advection equations in at least two dimensions. Here, we present an alternative computational strategy by combining a particle-based rather than a field-based approach with the initialisation of particles in proportion to their flux. This method allows accurate prediction of the steady state and is applicable even at high Péclet numbers where traditional particle-based Monte-Carlo methods starting from randomly initialised particle distributions fail. We demonstrate that generating a flux of particles according to a predetermined density and velocity distribution at a single fixed time and initial location allows for accurate simulation of mass transport under flow. Specifically, upon initialisation in proportion to their flux, these particles are propagated individually and detected by summing up their Monte-Carlo trajectories in predefined detection regions. We demonstrate quantitative agreement of the predicted concentration profiles with the results of experiments performed with fluorescent particles in microfluidic channels under continuous flow. This approach is computationally advantageous and readily allows non-trivial initial distributions to be considered. In particular, this method is highly suitable for simulating advective and diffusive transport in microfluidic devices.

研究の動機と目的

• 高Péclet数において従来の粒子ベースモンテカルロ法が初期化が不十分なために計算効率が著しく低下し、失敗する問題を解決すること。
• Fokker-Planck方程式を解くことなく、フラックスに基づく初期条件から直接粒子軌道をシミュレートする戦略を開発すること。
• 特にポアウイユ流れを示すマイクロフルイディクスデバイスにおいて、複雑な幾何形状でも定常状態の濃度プロファイルを正確に予測できること。
• マイクロチャネル内での蛍光コロイドを用いた実験的検証を行い、シミュレートされたプロファイルと測定データを比較すること。
• 測定済みの初期確率分布を入力として用いることで、フィッティングパラメータを最小限に抑えた手法のロバストネスを示すこと。

提案手法

• 粒子はt=0およびx=0で、フラックス重み付き分布Φ(x=0,y,z) = C(x=0,y,z)·v(y,z)/v̄に従って初期化される。ここでCは初期濃度、vは局所流速、v̄は平均流速である。
• 各粒子の軌道は、流れ場に基づく確定的広義と、確率的熱運動を組み込んだLangevin方程式を用いてシミュレートされる。
• 定常状態の濃度プロファイルは、下流の事前に定義された検出領域に粒子が何回時間ステップを費やしたかを数えることで再構築される。
• 全粒子を一度の初期時刻ステップからシミュレートすることで、連続的な粒子供給の必要性を回避し、計算コストを削減する。
• マイクロチャネル内の放物型速度プロファイルを考慮し、流れが遅い壁付近での検出バイアスを排除する。
• 粒子半径を除き、第一原理の物理法則に基づくシミュレーションであり、初期分布の測定値と粒子サイズのみを自由パラメータとして用い、実験的初期分布で検証済みである。

実験結果

リサーチクエスチョン

• RQ1 従来のモンテカルロ法が失敗する高Péclet数条件下でも、粒子ベースのシミュレーション手法が定常状態の物質移動を正確に予測できるか？
• RQ2 初期濃度ではなくフラックス分布に従って粒子を初期化することで、対流拡散移動のシミュレーションにおける精度と効率が向上するか？
• RQ3 実験的に測定された初期確率分布を入力として用いる場合、フィッティングなしに下流の濃度プロファイルをどの程度正確に予測できるか？
• RQ4 マイクロフルイディクスチャネル内でのシミュレートされた粒子分布は、実験的蛍光顕微鏡データとどの程度一致するか？
• RQ5 シミュレートされた拡散的幅広がりと測定されたものとの間で、定量的な一致度はどの程度か？

主な発見

• 本手法は、マイクロフルイディクスチャネル内でのシミュレート済み濃度プロファイルと実験的測定値との間で、非常に良好な定量的一致を達成しており、25 nmの既知の粒子半径を用いた場合に最も良い一致が得られた。
• 20 mm下流における25 nm粒子シミュレーションの正規化平方誤差は380であり、実験的ノイズレベルのおよそ20倍以内の一致を示している。
• 10 nmまたは40 nmの粒子半径を用いたシミュレーションでは、正規化誤差が数千～数万に達し、正しい粒子半径が使用された場合の感度と精度が確認された。
• 本手法は、80 mmにわたる伝播距離における蛍光コロイドビームの拡散的幅広がりを正確に捉えており、複数の下流位置での実験画像と一致した。
• 計算コストは幾何形状や境界条件にほとんど依存しないため、本手法は計算的に効率的かつスケーラブルである。
• 検証結果から、シミュレーションはデータへのフィッティングではなく、測定済みの初期分布と粒子サイズのみを用いた第一原理的予測であることが示された。