[論文レビュー] Modeling electrochemical systems with weakly imposed Dirichlet boundary conditions

Sun-Gu Kim, Khanwale, Makrand A.+2 more|arXiv (Cornell University)|Oct 17, 2020

Nanopore and Nanochannel Transport Studies参考文献 43被引用数 4

ひとこと要約

本論文では、ポアソン-ナーランスト・プランク方程式にディリクレ境界条件を弱い形で強制するディリクレ・トゥ・ノイマン変換を用いた有限要素フレームワークを提示する。この手法により、結合されたナビエ-ストークス-ポアソン-ナーランスト・プランク方程式に従う電気化学的系における電流フラックスを正確に計算できる。本手法は、粗いメッシュでも薄いデバイ層を正確に解像可能であり、強制的実装と比較して優れた収束性と低コストな計算を実現している。これは、3次元および不安定な電気済み流れの複雑な系に対しても同様に有効である。

ABSTRACT

Finite element modeling of charged species transport has enabled analysis, design, and optimization of a diverse array of electrochemical and electrokinetic devices. These systems are represented by the Poisson-Nernst-Planck equations coupled with the Navier-Stokes equation, with a key quantity of interest being the current at the system boundaries. Accurately computing the current flux is challenging due to the small critical dimension of the boundary layers (small Debye layer) that require fine mesh resolution at the boundaries. We resolve this challenge by using the Dirichlet-to-Neumanntransformation to weakly impose the Dirichlet conditions for the Poisson-Nernst-Planck equations. The results obtained with weakly imposed Dirichlet boundary conditions showed excellent agreement with those obtained when conventional boundary conditions with highly resolved mesh we reemployed. Furthermore, the calculated current flux showed faster mesh convergence using weakly imposed conditions compared to the conventionally imposed Dirichlet boundary conditions. We illustrate the approach on canonical 3D problems that otherwise would have been computationally intractable to solve accurately. This approach substantially reduces the computational cost of model-ing electrochemical systems.

研究の動機と目的

• 高長さ比の境界層を有する電気化学的系において、薄いデバイ層を解像するための計算的課題に対処すること。
• マルチスケール電気済み流れ系における境界での正確な電流フラックスをシミュレートするための計算コストを低減すること。
• 境界付近で粗いメッシュを用いても精度を維持できる堅牢な数値フレームワークの開発

提案手法

• ポアソン-ナーランスト・プランク方程式に対して、ディリクレ・トゥ・ノイマン変換を用いてディリクレ境界条件を弱い形で強制する。
• ナビエ-ストークス方程式とPNP方程式を効率的に結合するため、ブロック反復解法を採用する。
• 境界付近に要素を集中させるために、進行比1.021の幾何的メッシュグレーディングを適用し、極めて細かい解像度を必要としない。
• 検証のため、マニュフェイクド・ソリューションおよび解析的ベンチマーク（例：電気オスマシスフロー）を用いる。
• 1次元、2次元、3次元問題に本手法を適用し、特に陽性選択性膜付近の電気済み不安定性を含む。
• 非次元化された形で、時間刻みΔt = 1×10⁻⁶で連成系を解く。

実験結果

リサーチクエスチョン

• RQ1 粗いメッシュでも、薄いデバイ層を有する電気化学的系において、弱い形でのディリクレ条件の実装が正確な電流フラックスを達成できるか？
• RQ2 境界付近でのメッシュの細分化に伴う収束速度について、弱い実装法と強い実装法とを比較すると、どちらが優れているか？
• RQ3 極めて細かいメッシュの細分化を最小限に抑えても、イオン濃度ポテンシャルの偏極や電気対流不安定性といった複雑な電気済み現象を正確にシミュレートできるか？
• RQ4 流動条件や境界層の厚さが変化しても、3次元幾何形状において本手法は精度を維持できるか？

主な発見

• 弱い実装法は、デバイ層にたった1〜2要素しか割り当てない場合（Λ = 1×10⁻³）でも、強い実装法の結果と優れた一致を示す。一方、強い実装法では顕著に細かいメッシュが必要となる。
• 境界付近でのメッシュ細分化に伴う収束速度は、強い実装法よりも本手法が速く、優れた数値的効率性を示している。
• 弱い実装法で計算された境界フラックスは、解析的解およびマニュフェイクド・ソリューションと比較して、誤差が著しく小さい（残差誤差〜10⁻⁸）。
• 粗いメッシュ（1280×180の四角形要素）を用いても、カチオン選択性膜付近の電気済み不安定性を正確に再現でき、先行研究のベンチマークで観測された指状の電荷密度構造を再現した。
• 脱塩用の3次元分岐マイクロチャネルのシミュレーションにより、本フレームワークのスケーラビリティおよび複雑な幾何形状におけるロバスト性を示した。
• 極めて細かい境界層メッシュを回避することで、計算コストを顕著に低減し、かつては計算不能であった3次元電気化学的シミュレーションを可能にした。