KMS Of Academy of mathematics and systems sciences, CAS
An inverse averaging finite element method for solving three-dimensional Poisson-Nernst-Planck equations in nanopore system simulations | |
Zhang, Qianru1,2; Wang, Qin1,2; Zhang, Linbo1,2; Lu, Benzhuo1,2 | |
2021-11-21 | |
Source Publication | JOURNAL OF CHEMICAL PHYSICS
![]() |
ISSN | 0021-9606 |
Volume | 155Issue:19Pages:13 |
Abstract | The Poisson-Nernst-Planck (PNP) model plays an important role in simulating nanopore systems. In nanopore simulations, the large-size nanopore system and convection-domination Nernst-Planck (NP) equations will bring convergence difficulties and numerical instability problems. Therefore, we propose an improved finite element method (FEM) with an inverse averaging technique to solve the three-dimensional PNP model, named inverse averaging FEM (IAFEM). At first, the Slotboom variables are introduced aiming at transforming non-symmetric NP equations into self-adjoint second-order elliptic equations with exponentially behaved coefficients. Then, these exponential coefficients are approximated with their harmonic averages, which are calculated with an inverse averaging technique on every edge of each tetrahedral element in the grid. Our scheme shows good convergence when simulating single or porous nanopore systems. In addition, it is still stable when the NP equations are convection domination. Our method can also guarantee the conservation of computed currents well, which is the advantage that many stabilization schemes do not possess. Our numerical experiments on benchmark problems verify the accuracy and robustness of our scheme. The numerical results also show that the method performs better than the standard FEM when dealing with convection-domination problems. A successful simulation combined with realistic chemical experiments is also presented to illustrate that the IAFEM is still effective for three-dimensional interconnected nanopore systems. |
DOI | 10.1063/5.0066194 |
Indexed By | SCI |
Language | 英语 |
Funding Project | National Key R&D Program of China[2019YFA0709600] ; National Key R&D Program of China[2019YFA0709601] ; National Natural Science Foundation of China[11771435] ; National Natural Science Foundation of China[22073110] |
WOS Research Area | Chemistry ; Physics |
WOS Subject | Chemistry, Physical ; Physics, Atomic, Molecular & Chemical |
WOS ID | WOS:000719748700001 |
Publisher | AIP Publishing |
Citation statistics | |
Document Type | 期刊论文 |
Identifier | http://ir.amss.ac.cn/handle/2S8OKBNM/59575 |
Collection | 中国科学院数学与系统科学研究院 |
Corresponding Author | Lu, Benzhuo |
Affiliation | 1.Chinese Acad Sci, Acad Math & Syst Sci, NCMIS, CEMS,LSEC, Beijing, Peoples R China 2.Univ Chinese Acad Sci, Sch Math Sci, Beijing, Peoples R China |
Recommended Citation GB/T 7714 | Zhang, Qianru,Wang, Qin,Zhang, Linbo,et al. An inverse averaging finite element method for solving three-dimensional Poisson-Nernst-Planck equations in nanopore system simulations[J]. JOURNAL OF CHEMICAL PHYSICS,2021,155(19):13. |
APA | Zhang, Qianru,Wang, Qin,Zhang, Linbo,&Lu, Benzhuo.(2021).An inverse averaging finite element method for solving three-dimensional Poisson-Nernst-Planck equations in nanopore system simulations.JOURNAL OF CHEMICAL PHYSICS,155(19),13. |
MLA | Zhang, Qianru,et al."An inverse averaging finite element method for solving three-dimensional Poisson-Nernst-Planck equations in nanopore system simulations".JOURNAL OF CHEMICAL PHYSICS 155.19(2021):13. |
Files in This Item: | There are no files associated with this item. |
Items in the repository are protected by copyright, with all rights reserved, unless otherwise indicated.
Edit Comment