Jehanzeb H. Chaudhry

Jehanzeb H. Chaudhry, J. B. Collins, John N. Shadid

Implicit-Explicit (IMEX) schemes are widely used for time integration methods for approximating solutions to a large class of problems. In this work, we develop accurate a posteriori error estimates of a quantity of interest for approximations obtained from multi-stage IMEX schemes. This is done by first defining a finite element method that is nodally equivalent to an IMEX scheme, then using typical methods for adjoint-based error estimation. The use of a nodally equivalent finite element method allows a decomposition of the error into multiple components, each describing the effect of a different portion of the method on the total error in a quantity of interest.

Jehanzeb H. Chaudhry

The Poisson-Boltzmann equation (PBE) models the electrostatic interactions of charged bodies such as molecules and proteins in an electrolyte solvent. The PBE is a challenging equation to solve numerically due to the presence of singularities, discontinuous coefficients and boundary conditions. Hence, there is often large error in the numerical solution of the PBE that needs to be quantified. In this work, we use adjoint based a posteriori analysis to accurately quantify the error in an important quantity of interest, the solvation free energy, for the finite element solution of the PBE. We identify various sources of error and propose novel refinement strategies based on a posteriori error estimates.