Abstract:In this work, we consider numerical methods for the Poisson–Nernst–Planck–Cahn–Hilliard (PNPCH) equations with steric interactions, which correspond to a non-constant mobility H−1 gradient flow of a free-energy functional that consists of electrostatic free energies, steric interaction energies of short range, entropic contribution of ions, and concentration gradient energies. We propose a novel energy stable numerical scheme that respects mass conservation and positivity at the discrete level. Existence and uniqueness of the solution to the proposed nonlinear scheme are established by showing that the solution is a unique minimizer of a convex functional over a closed, convex domain. The positivity of numerical solutions is further theoretically justified by the singularity of the entropy terms, which prevents the minimizer from approaching zero concentrations. A further numerical analysis proves discrete free-energy dissipation. Extensive numerical tests are performed to validate that the numerical scheme is first-order accurate in time and second-order accurate in space, and is capable of preserving the desired properties, such as mass conservation, positivity, and free energy dissipation, at the discrete level. Moreover, the PNPCH equations and the proposed scheme are applied to study charge dynamics and self-assembled nanopatterns in highly concentrated electrolytes that are widely used in electrochemical energy devices. Numerical results demonstrate that the PNPCH equations and our numerical scheme are able to capture nanostructures, such as lamellar patterns and labyrinthine patterns in electric double layers and the bulk, and multiple time relaxation with multiple time scales. The multiple time relaxation dynamics with metastability take long time to reach an equilibrium, highlighting the need for robust, energy stable numerical schemes that allow large time stepping. In addition, we numerically characterize the interplay between cross steric interactions of short range and the concentration gradient regularization, and their impact on the development of nanostructures in the equilibrium state.

Energy Dissipative and Positivity Preserving Schemes for Large-Convection Ion Transport with Steric and Solvation Effects.

Structure-Preserving and Efficient Numerical Methods for Ion Transport

Second-order, Positive, and Unconditional Energy Dissipative Scheme for Modified Poisson-Nernst-Planck Equations

Second-order, Positive, and Unconditional Energy Dissipative Scheme for Modified Poisson–Nernst–Planck Equations

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

Positivity Preserving Finite Difference Methods for Poisson-Nernst-Planck Equations with Steric Interactions: Application to Slit-Shaped Nanopore Conductance

A linear, second-order accurate, positivity-preserving and unconditionally energy stable scheme for the Navier–Stokes–Poisson–Nernst–Planck system

Efficient, positive, and energy stable schemes for multi-D Poisson-Nernst-Planck systems

A Positivity-Preserving, Energy Stable and Convergent Numerical Scheme for the Poisson-Nernst-Planck System

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

A Modified Poisson--Nernst--Planck Model with Excluded Volume Effect: Theory and Numerical Implementation

Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations

Convergence Analysis on a Structure-Preserving Numerical Scheme for the Poisson-Nernst-Planck-Cahn-Hilliard System

A Positivity-Preserving and Free Energy Dissipative Hybrid Scheme for the Poisson-Nernst-Planck Equations on Polygonal and Polyhedral Meshes

A positivity-preserving, linear, energy stable and convergent numerical scheme for the Poisson–Nernst–Planck (PNP) system

A stabilized finite volume element method for solving Poisson–Nernst–Planck equations

Structure-preserving numerical method for Maxwell-Ampère Nernst-Planck model

A Second Order Accurate, Positivity Preserving Numerical Method for the Poisson–Nernst–Planck System and Its Convergence Analysis

Unconditional positivity-preserving and energy stable schemes for a reduced Poisson-Nernst-Planck system

Ion Transport in Dipolar Medium I: A Local Dielectric Poisson–Nernst–Planck/Poisson–Boltzmann Model

Poisson-Nernst-Planck Systems for Ion Flow with a Local Hard-Sphere Potential for Ion Size Effects.