Adaptive time-stepping schemes for the solution of the Poisson-Nernst-Planck equations

DOI10.1016/J.APNUM.2021.01.018zbMATH Open1462.65116arXiv1703.10297OpenAlexW3126406745MaRDI QIDQ1995976FDOQ1995976

Authors: David Yan, M. C. Pugh, F. P. Dawson 

Publication date: 2 March 2021

Published in: Applied Numerical Mathematics (Search for Journal in Brave)

Abstract: The Poisson-Nernst-Planck equations with generalized Frumkin-Butler-Volmer boundary conditions (PNP-FBV) describe ion transport with Faradaic reactions, and have applications in a number of fields. In this article, we develop an adaptive time-stepping scheme for the solution of the PNP-FBV equations based on two time-stepping methods: a fully implicit (BDF2) method, and an implicit-explicit (SBDF2) method. We present simulations under both current and voltage boundary conditions and demonstrate the ability to simulate a large range of parameters, including any value of the singular perturbation parameter $epsilon$. When the underlying dynamics is one that would have the solutions converge to a steady-state solution, we observe that the adaptive time-stepper based on the SBDF2 method produces solutions that ``nearly converge to the steady state and that, simultaneously, the time-step sizes stabilize to a limiting size $d{t}_{i}nfty$. In the companion to this article cite{YPD_Part2}, we linearize the SBDF2 scheme about the steady-state solution and demonstrate that the linearized scheme is conditionally stable. This conditional stability is the cause of the adaptive time-stepper's behaviour. While the adaptive time-stepper based on the fully-implicit (BDF2) method is not subject to such time-step constraints, the required nonlinear solve yields run times that are significantly longer.

Full work available at URL: https://arxiv.org/abs/1703.10297

Recommendations

• A study of the numerical stability of an ImEx scheme with application to the Poisson-Nernst-Planck equations
• Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations
• An adaptive conservative finite volume method for Poisson-Nernst-Planck equations on a moving mesh
• A conservative discretization of the Poisson-Nernst-Planck equations on adaptive Cartesian grids
• Poisson-Nernst-Planck equations for simulating biomolecular diffusion-reaction processes. I: Finite element solutions

zbMATH Keywords

Poisson-Nernst-Planck equationsadaptive time-steppingsemi-implicit methodsImEx methods

Mathematics Subject Classification ID

PDEs in connection with optics and electromagnetic theory (35Q60) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Mesh generation, refinement, and adaptive methods for the numerical solution of initial value and initial-boundary value problems involving PDEs (65M50) Electrochemistry (78A57)

Cites Work

• Title not available (Why is that?)
• Spectral methods for incompressible viscous flow
• Numerical initial value problems in ordinary differential equations.
• Implicit-Explicit Methods for Time-Dependent Partial Differential Equations
• The Theoretical Accuracy of Runge–Kutta Time Discretizations for the Initial Boundary Value Problem: A Study of the Boundary Error
• Approximation theory and approximation practice
• Current-Voltage Relations for Electrochemical Thin Films
• Une méthode multipas implicite-explicite pour l'approximation des équations d'évolution paraboliques
• Variable step-size implicit-explicit linear multistep methods for time-dependent partial differential equations
• A fully implicit, fully adaptive time and space discretisation method for phase-field simulation of binary alloy solidification
• On the Removal of Boundary Errors Caused by Runge–Kutta Integration of Nonlinear Partial Differential Equations
• A BDF2 integration method with step size control for elasto-plasticity
• Electrochemical Thin Films at and above the Classical Limiting Current
• A highly accurate adaptive finite difference solver for the Black-Scholes equation
• \texttt{Chebfun}: A new kind of numerical computing
• Digital simulation in electrochemistry
• The analysis of the generalized-\(\alpha\) method for nonlinear dynamic problems
• A study of the numerical stability of an ImEx scheme with application to the Poisson-Nernst-Planck equations

Cited In (4)

• A study of the numerical stability of an ImEx scheme with application to the Poisson-Nernst-Planck equations
• A multi-scale simulation of retinal physiology
• Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations
• Convergence and superconvergence analysis for a mass conservative, energy stable and linearized BDF2 scheme of the Poisson-Nernst-Planck equations

Uses Software

• Matlab
• Chebfun
• GitHub

This page was built for publication: Adaptive time-stepping schemes for the solution of the Poisson-Nernst-Planck equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1995976)