Department of Mathematics
In the last two decades, the field of geometric numerical integration and structure-preserving algorithms has focused on the design of numerical methods that preserve properties of Hamiltonian systems, evolution problems on manifolds, and problems with highly oscillatory solutions. In this paper, we show that a different geometric property, namely, the blow up of solutions in finite time, can also be taken into account in the numerical integrator, giving rise to geometric methods we call B-methods. We give a first systematic approach for deriving such methods for scalar and systems of semi- and quasi-linear parabolic and hyperbolic partial differential equations. We show both analytically and numerically that B-methods have substantially better approximation properties than standard numerical integrators as the solution approaches the blow-up time.
geometric integration, blow-up solutions, nonlinear partial differential equations, nonlinear systems of equations
Source Publication Title
SIAM Journal on Scientific Computing
Society for Industrial and Applied Mathematics
Link to Publisher's Edition
Beck, M., Gander, M., & Kwok, F. (2015). B-Methods for the Numerical Solution of Evolution Problems with Blow-Up Solutions Part I: Variation of the Constant. SIAM Journal on Scientific Computing, 37 (6), A2298-A3029. https://doi.org/10.1137/15M1011767