A finite difference method and analysis for 2D nonlinear Poisson-Boltzmann equations

Zhilin Li, C. V. Pao, Zhonghua Qiao

Research output: Journal article publicationJournal articleAcademic researchpeer-review

11 Citations (Scopus)


A fast finite difference method based on the monotone iterative method and the fast Poisson solver on irregular domains for a 2D nonlinear Poisson-Boltzmann equation is proposed and analyzed in this paper. Each iteration of the monotone method involves the solution of a linear equation in an exterior domain with an arbitrary interior boundary. A fast immersed interface method for generalized Helmholtz equations on exterior irregular domains is used to solve the linear equation. The monotone iterative method leads to a sequence which converges monotonically from either above or below to a unique solution of the problem. This monotone convergence guarantees the existence and uniqueness of a solution as well as the convergence of the finite difference solution to the continuous solution. A comparison of the numerical results against the exact solution in an example indicates that our method is second order accurate. We also compare our results with available data in the literature to validate the numerical method. Our method is efficient in terms of accuracy, speed, and flexibility in dealing with the geometry of the domain.
Original languageEnglish
Pages (from-to)61-81
Number of pages21
JournalJournal of Scientific Computing
Issue number1
Publication statusPublished - 1 Jan 2007
Externally publishedYes


  • Exterior problem
  • Immersed interface method
  • Monotone iterative method
  • Nonlinear Poisson-Boltzmann equation

ASJC Scopus subject areas

  • Theoretical Computer Science
  • Software
  • Engineering(all)
  • Computational Theory and Mathematics

Cite this