Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

Abstract

We have developed efficient numerical algorithms for solving 3D steady-state Poisson-Nernst-Planck (PNP) equations with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by a finite difference scheme and solved iteratively using the Gummel method with relaxation. The Nernst-Planck equations are transformed into Laplace equations through the Slotboom transformation. Then, the algebraic multigrid method is applied to efficiently solve the Poisson equation and the transformed Nernst-Planck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed, which reduces computational complexity from $O(N^2)$ to $O(N\log N)$, where $N$ is the number of grid points. Integrals involving the Dirac delta function are evaluated directly by coordinate transformation, which yields more accurate results compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for lithium-ion (Li-ion) batteries are shown to be in good agreement with the experimental data and the results from previous studies.

Keywords and phrases: Poisson-Nernst-Planck equations, classical density functional theory, algebraic multigrid method, fast Fourier transform, Li-ion battery. 26 pages, 8 figures

References in corpus (2)

• P3DFFT: a framework for parallel computations of Fourier transforms in three dimensions
• Ionic channels in biological membranes. Electrostatic analysis of a natural nanotube

Cited by in corpus (12)

• NWChem: Past, Present, and Future
• A Positive and Energy Stable Numerical Scheme for the Poisson-Nernst-Planck-Cahn-Hilliard Equations with Steric Interactions
• Modeling transport of charged species in pore networks: solution of the Nernst-Planck equations coupled with fluid flow and charge conservation equations
• A local approximation of fundamental measure theory incorporated into three dimensional Poisson-Nernst-Planck equations to account for hard sphere repulsion among ions
• Positivity-preserving third order DG schemes for Poisson--Nernst--Planck equations
• Structure-Preserving and Efficient Numerical Methods for Ion Transport
• Impedance resonance in narrow confinement
• High-order space-time finite element methods for the Poisson-Nernst-Planck equations: Positivity and unconditional energy stability
• Interactive Multiscale Modeling to Bridge Atomic Properties and Electrochemical Performance in Li-CO${}_2$ Battery Design
• Efficient, positive, and energy stable schemes for multi-D Poisson-Nernst-Planck systems
• A Modified Poisson--Nernst--Planck Model with Excluded Volume Effect: Theory and Numerical Implementation
• The role of correlation and solvation in ion interactions with B-DNA