Five-Structures Preserving Algorithm for charge dynamics model
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866911578547290112 |
|---|---|
| author | Sun, Haoran Wu, Wancheng Wang, Kun |
| author_facet | Sun, Haoran Wu, Wancheng Wang, Kun |
| contents | This paper develops a family of fast, structure-preserving numerical algorithms for the nonlinear Maxwell-Ampere Nernst-Planck equations. For the first-order scheme, the Slotboom transformation rewrites the Nernst-Planck equation to enable positivity preservation. The backward Euler method and centered finite differences discretize the transformed system. Two correction strategies are introduced: one enforces Gauss's law via a displacement correction, and the other preserves Faraday's law through potential reconstruction. The fully discrete scheme exactly satisfies mass conservation, concentration positivity, energy dissipation, Gauss's law, and Faraday's law, with established error estimates. The second-order scheme adopts BDF2 time discretization while retaining the same structure-preserving strategies, exactly conserving mass, Gauss's law, and Faraday's law. Numerical experiments validate both schemes using analytical solutions, confirming convergence orders and positivity preservation. Simulations of ion transport with fixed charges demonstrate exact preservation of Gauss's and Faraday's laws over long-time evolution, reproducing electrostatic attraction, ion accumulation, and electric field screening. The results fully support the theoretical analysis and the schemes' stability and superior performance. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2604_08228 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Five-Structures Preserving Algorithm for charge dynamics model Sun, Haoran Wu, Wancheng Wang, Kun Numerical Analysis This paper develops a family of fast, structure-preserving numerical algorithms for the nonlinear Maxwell-Ampere Nernst-Planck equations. For the first-order scheme, the Slotboom transformation rewrites the Nernst-Planck equation to enable positivity preservation. The backward Euler method and centered finite differences discretize the transformed system. Two correction strategies are introduced: one enforces Gauss's law via a displacement correction, and the other preserves Faraday's law through potential reconstruction. The fully discrete scheme exactly satisfies mass conservation, concentration positivity, energy dissipation, Gauss's law, and Faraday's law, with established error estimates. The second-order scheme adopts BDF2 time discretization while retaining the same structure-preserving strategies, exactly conserving mass, Gauss's law, and Faraday's law. Numerical experiments validate both schemes using analytical solutions, confirming convergence orders and positivity preservation. Simulations of ion transport with fixed charges demonstrate exact preservation of Gauss's and Faraday's laws over long-time evolution, reproducing electrostatic attraction, ion accumulation, and electric field screening. The results fully support the theoretical analysis and the schemes' stability and superior performance. |
| title | Five-Structures Preserving Algorithm for charge dynamics model |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2604.08228 |