Structure-preserving numerical method for Ampere-Nernst-Planck model

04/25/2022
by   Zhonghua Qiao, et al.
0

Charge dynamics play essential role in many practical applications such as semiconductors, electrochemical devices and transmembrane ion channels. An Ampére-Nernst-Planck (ANP) model that describes charge dynamics via concentrations and the electric displacement is able to take effects beyond mean-field approximations into account. To obtain physically faithful numerical solutions, we develop a structure-preserving numerical method for the ANP model whose solution has several physical properties of importance. By the Slotboom transform with entropy-mean approximations, a positivity preserving scheme with Scharfetter-Gummel fluxes is derived for the generalized Nernst-Planck equations. To deal with the curl-free constraint, the dielectric displacement from the Ampére equation is further updated with a local relaxation algorithm of linear computational complexity. We prove that the proposed numerical method unconditionally preserves the mass conservation and the solution positivity at the discrete level, and satisfies the discrete energy dissipation law with a mild time-step restriction. Numerical experiments verify that our numerical method has expected accuracy and structure-preserving properties. Applications to ion transport with large convection, arising from boundary-layer electric field and Born solvation interactions, further demonstrate that the ANP formulation with the proposed numerical scheme has attractive performance and can effectively describe charge dynamics with large convection of high numerical cell Péclet numbers.

READ FULL TEXT

page 14

page 15

page 16

research
04/21/2021

A convergent numerical scheme for a model of liquid crystal dynamics subjected to an electric field

We present a convergent and constraint-preserving numerical discretizati...
research
09/29/2020

Numerical approximations of chromatographic models

A numerical scheme based on modified method of characteristics with adju...
research
05/03/2022

A dynamic mass transport method for Poisson-Nernst-Planck equations

A dynamic mass-transport method is proposed for approximately solving th...
research
12/11/2020

Trefftz Functions for Nonlocal Electrostatics

Electrostatic interactions in solvents play a major role in biophysical ...
research
02/22/2020

A Positive and Energy Stable Numerical Scheme for the Poisson-Nernst-Planck-Cahn-Hilliard Equations with Steric Interactions

We consider numerical methods for the Poisson-Nernst-Planck-Cahn-Hilliar...
research
11/18/2019

A structure preserving numerical scheme for Fokker-Planck equations of neuron networks: numerical analysis and exploration

In this work, we are concerned with the Fokker-Planck equations associat...
research
10/08/2019

Charge-conserving, variational particle-in-cell method for the drift-kinetic Vlasov-Maxwell system

This paper proposes a charge-conserving, variational, spatio-temporal di...

Please sign up or login with your details

Forgot password? Click here to reset