Adaptive Finite Element Methods with Inexact Solvers for the Nonlinear Poisson-Boltzmann Equation

Abstract

In this article we study adaptive finite element methods (AFEM) with inexact solvers for a class of semilinear elliptic interface problems. We are particularly interested in nonlinear problems with discontinuous diffusion coefficients, such as the nonlinear Poisson-Boltzmann equation and its regularizations. The algorithm we study consists of the standard SOLVE-ESTIMATE-MARK-REFINE procedure common to many adaptive finite element algorithms, but where the SOLVE step involves only a full solve on the coarsest level, and the remaining levels involve only single Newton updates to the previous approximate solution. We summarize a recently developed AFEM convergence theory for inexact solvers, and present a sequence of numerical experiments that give evidence that the theory does in fact predict the contraction properties of AFEM with inexact solvers. The various routines used are all designed to maintain a linear-time computational complexity.

Submitted to DD20 Proceedings

Cited by in corpus (6)

• Dielectric continuum methods for quantum chemistry
• Modeling Charge-Sign Asymmetric Solvation Free Energies With Nonlinear Boundary Conditions
• Coupling finite and boundary element methods to solve the Poisson--Boltzmann equation for electrostatics in molecular solvation
• Weak formulations of the nonlinear Poisson-Boltzmann equation in biomolecular electrostatics
• Enriched Gradient Recovery for Interface Solutions of the Poisson-Boltzmann Equation
• Reliable numerical solution of a class of nonlinear elliptic problems generated by the Poisson-Boltzmann equation