Abstract:This work considers charged systems described by the modified Poisson--Nernst--Planck (PNP) equations, which incorporate ionic steric effects and the Born solvation energy for dielectric inhomogeneity. Solving the steady-state modified PNP equations poses numerical challenges due to the emergence of sharp boundary layers caused by small Debye lengths, particularly when local ionic concentrations reach saturation. To address this, we first reformulate the steady-state problem as a constraint optimization, where the ionic concentrations on unstructured Delaunay nodes are treated as fractional particles moving along edges between nodes. The electric fields are then updated to minimize the objective free energy while satisfying the discrete Gauss's law. We develop a local relaxation method on unstructured meshes that inherently respects the discrete Gauss's law, ensuring curl-free electric fields. Numerical analysis demonstrates that the optimal mass of the moving fractional particles guarantees the positivity of both ionic and solvent concentrations. Additionally, the free energy of the charged system consistently decreases during successive updates of ionic concentrations and electric fields. We conduct numerical tests to validate the expected numerical accuracy, positivity, free-energy dissipation, and robustness of our method in simulating charged systems with sharp boundary layers.

Existence of C Singular Solutions to Euler-Nernst-Planck-Poisson System on R with Free-Moving Charges

Existence of $C^{1,\alpha}$ Singular Solutions to Euler-Nernst-Planck-Poisson System on $\mathbb{R}^3$ with Free-Moving Charges

Existence of C1, singular solutions to Euler-Nernst-Planck-Poisson system on R3 with free-moving charges

Existence of C1,α Singular Solutions to Euler–Nernst–Planck–Poisson System on R

An Approach to Numerically Solving the Poisson Equation

Global Existence of Large Solutions for the 3D incompressible Navier--Stokes--Poisson--Nernst--Planck Equations

Structure of singularities for the Euler-Poisson system of ion dynamics

Global existence and asymptotic behavior of the solutions to the three-dimensional bipolar Euler–Poisson systems

Singularity formation for the relativistic Euler‐Poisson equations with repulsive force and damping

GLOBAL EXISTENCE AND ASYMPTOTIC BEHAVIOR OF THE FULL EULER SYSTEM WITH DAMPING AND RADIATIVE EFFECTS IN R3

Formation of singularities in plasma ion dynamics

Global solutions to the Nernst-Planck-Euler system on bounded domain

Multiple Blowing-Up and Concentrating Solutions for Liouville-type Equations with Singular Sources under Mixed Boundary Conditions

Solutions for planar Kirchhoff-Schrödinger-Poisson systems with general nonlinearities

Stability of planar diffusion wave for three‐dimensional bipolar Euler–Poisson system with time‐dependent over‐damping

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

Existence of Euler-Korteweg equations with the non-monotone pressure

Existence of positive solution to Kirchhoff-Schrödinger-Poisson system with strong singular term

Nonexistence of multi-dimensional solitary waves for the Euler–Poisson system

Finite time blowup explicit solution for 3D incompressible Euler equations

Local Structure-Preserving Relaxation Method for Charged Systems on Unstructured Meshes