A positivity preserving scheme for Poisson-Nernst-Planck Navier-Stokes equations and its error analysis

Ziyao Yu,Qing Cheng,Jie Shen,Changyou Wang

DOI: https://doi.org/10.48550/arXiv.2311.17349

2023-11-29

Abstract:We consider in this paper a numerical approximation of Poisson-Nernst-Planck-Navier- Stokes (PNP-NS) system. We construct a decoupled semi-discrete and fully discrete scheme that enjoys the properties of positivity preserving, mass conserving, and unconditionally energy stability. Then, we establish the well-posedness and regularity of the initial and (periodic) boundary value problem of the PNP-NS system under suitable assumptions on the initial data, and carry out a rigorous convergence analysis for the fully discretized scheme. We also present some numerical results to validate the positivity-preserving property and the accuracy of our scheme.

Numerical Analysis,Analysis of PDEs

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is to propose a new numerical approximation scheme for the Poisson - Nernst - Planck - Navier - Stokes (PNP - NS) system. Specifically, the authors aim to construct a discrete scheme to ensure the following three key properties:

1. **Positivity Preserving**: Ensure that the ion concentrations \(p\) and \(n\) are always positive during the calculation process.
2. **Mass Conserving**: Guarantee that the mass of ions in the system does not lose or increase with the change of time during the numerical simulation process.
3. **Unconditionally Energy Stable**: Ensure that the numerical scheme can maintain the energy dissipation characteristics of the system at any time step.

In addition, the authors also strive to establish the well - posedness and regularity results of the PNP - NS system and conduct a strict convergence analysis of the proposed fully - discrete scheme. To achieve these goals, the authors adopt the following techniques:

- Rewrite the PNP system in the form of Wasserstein gradient flow and use related techniques to maintain positivity and energy stability.
- Use the projection method to decouple the velocity and pressure.
- Introduce an additional \(O(\Delta t)\) term to explicitly handle the convection term in the PNP equation while maintaining stability.

Through these methods, the authors not only propose a new numerical scheme but also verify its theoretical properties and practical application effects, thus providing a powerful tool for research in fields such as electrochemical systems and biomembrane channels.