Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
194 tokens/sec
GPT-4o
7 tokens/sec
Gemini 2.5 Pro Pro
45 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

Local Structure-Preserving Relaxation Method for Charged Systems on Unstructured Meshes (2310.05411v1)

Published 9 Oct 2023 in math.NA and cs.NA

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.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (37)
  1. A finite-element framework for a mimetic finite-difference discretization of Maxwell’s equations. SIAM J. Sci. Comput., 43:A2638–A2659, 2021.
  2. Simple and robust solver for the Poisson-Boltzmann equation. Phys. Rev. E, 80:016705, 2009.
  3. A stabilized finite element method for the numerical simulation of multi-ion transport in electrochemical systems. Comput. Methods Appl. Mech. Engrg., 223-224:199–210, 2012.
  4. Diffuse-charge dynamics in electrochemical systems. Phys. Rev. E, 70(2):021506, 2004.
  5. Steric effects in electrolytes: A modified Poisson-Boltzmann equation. Phys. Rev. Lett., 79:435–438, 1997.
  6. PUMIPic: A mesh-based approach to unstructured mesh Particle-In-Cell on GPUs. J. Parallel Distrib. Comput., 157:1–12, 2021.
  7. Positivity preserving finite difference methods for Poisson–Nernst–Planck equations with steric interactions: Application to slit-shaped nanopore conductance. J. Comput. Phys., 397:108864, 2019.
  8. B. Eisenberg. Ionic channels in biological membranes-electrostatic analysis of a natural nanotube. Contemp. Phys., 39(6):447–466, 1998.
  9. PNP equations with steric effects: a model of ion flow through channels. J. Phys. Chem. B, 116(37):11422–11441, 2012.
  10. A multi-scale approach to model K(+) permeation through the KcsA channel. Front. Mol. Biosci., 9:880660, 2022.
  11. J. Hu and X. Huang. A fully discrete positivity-preserving and energy-dissipative finite difference scheme for Poisson–Nernst–Planck equations. Numer. Math., 145:77–115, 2020.
  12. Steric effects in the dynamics of electrolytes at large applied voltages. II. Modified Poisson–Nernst–Planck equations. Phys. Rev. E, 75:021503, 2007.
  13. B. Li. Continuum electrostatics for ionic solutions with nonuniform ionic sizes. Nonlinearity, 22:811–833, 2009.
  14. B. Li. Minimization of electrostatic free energy and the Poisson-Boltzmann equation for molecular solvation with implicit solvent. SIAM J. Math. Anal., 40:2536–2566, 2009.
  15. A positivity-preserving, energy stable and convergent numerical scheme for the Poisson-Nernst-Planck system. Math. Comput., 90:2071–2106, 2021.
  16. H. Liu and W. Maimaitiyiming. Efficient, positive, and energy stable schemes for multi-D Poisson–Nernst–Planck systems. J. Sci. Comput., 87(3):1–36, 2021.
  17. J. Liu and B. Eisenberg. Molecular mean-field theory of ionic solutions: A Poisson-Nernst-Planck-Bikerman model. Entropy, 22:550, 2020.
  18. Modified Poisson–Nernst–Planck model with accurate Coulomb correlation in variable media. SIAM J. Appl. Math., 78:226–245, 2018.
  19. B. Lu and Y. Zhou. Poisson-Nernst-Planck equations for simulating biomolecular diffusion-reaction processes II: Size effects on ionic distributions and diffusion-reaction rates. Biophys. J., 100:2475–2485, 2011.
  20. Modified Poisson–Nernst–Planck model with Coulomb and hard-sphere correlations. SIAM J. Appl. Math., 81:1645–1667, 2021.
  21. A. C. Maggs and V. Rossetto. Local simulation algorithms for Coulomb interactions. Phys. Rev. Lett., 88:196402, 2002.
  22. Semiconductor equations. Springer Science & Business Media, 2012.
  23. Energetically stable discretizations for charge transport and electrokinetic models. J. Comput. Phys., 306:1–18, 2016.
  24. J. Nocedal and S. Wright. Numerical Optimization. Springer Verlag, New York, 1999.
  25. I. Pasichnyk and B. Dünweg. Coulomb interactions via local dynamics: A molecular-dynamics algorithm. J. Phys: Condens. Matter, 16:S3999, 2004.
  26. A positive and energy stable numerical scheme for the Poisson–Nernst–Planck–Cahn–Hilliard equations with steric interactions. J. Comput. Phys., 426:109908, 2021.
  27. Ion-flow field calculation of HVDC overhead lines using a high-order stabilisation technique based on Petrov–Galerkin method. IET Gener. Transm. Distrib., 12:1183–1189, 2018.
  28. A Maxwell-Ampère Nernst-Planck framework for modeling charge dynamics. SIAM J. Appl. Math., 83:374–393, 2023.
  29. Structure-preserving numerical method for Maxwell-Ampère Nernst-Planck model. J. Comput. Phys., 475:111845, 2023.
  30. A finite element framework for some mimetic finite difference discretizations. Comput. Math. Appl., 70(11):2661–2673, 2015.
  31. Transport phenomena in nanofluidics. Rev. Mod. Phys., 80:839–883, 2008.
  32. J. Shen and J. Xu. Unconditionally positivity preserving and energy dissipative schemes for Poisson–Nernst–Planck equations. Numer. Math., 148:671–697, 2021.
  33. M. Stynes. Steady-state convection-diffusion problems. Acta Numer., 14:445–508, 2005.
  34. P. N. Vabishchevich. Finite-difference approximation of mathematical physics problems on irregular grids. Comput. Methods Appl. Math., 5(3):294–330, 2005.
  35. S. Wang. Solving convection-dominated anisotropic diffusion equations by an exponentially fitted finite volume method. Comput. Math. Appl., 44(8):1249–1265, 2002.
  36. Structures of hard-sphere fluids from a modified fundamental-measure theory. J. Chem. Phys., 117(22):10156–10164, 2002.
  37. A class of finite element methods with averaging techniques for solving the three-dimensional drift-diffusion model in semiconductor device simulations. J. Comput. Phys., 458:111086, 2022.

Summary

We haven't generated a summary for this paper yet.

Ai Sparkles Streamline Icon: https://streamlinehq.com Summarize Now