Unconditional Stability for Multistep ImEx Schemes: Theory

This paper presents a new class of high order linear ImEx multistep schemes with large regions of unconditional stability. Unconditional stability is a desirable property of a time stepping scheme, as it allows the choice of time step solely based on accuracy considerations. Of particular interest are problems for which both the implicit and explicit parts of the ImEx splitting are stiff. Such splittings can arise, for example, in variable-coefficient problems, or the incompressible Navier-Stokes equations. To characterize the new ImEx schemes, an unconditional stability region is introduced, which plays a role analogous to that of the stability region in conventional multistep methods. Moreover, computable quantities (such as a numerical range) are provided that guarantee an unconditionally stable scheme for a proposed implicit-explicit matrix splitting. The new approach is illustrated with several examples. Coefficients of the new schemes up to fifth order are provided.Comment: 33 pages, 7 figure

A Modified Split Bregman Algorithm for Computing Microstructures Through Young Measures

The goal of this paper is to describe the oscillatory microstructure that can emerge from minimizing sequences for nonconvex energies. We consider integral functionals that are defined on real valued (scalar) functions $u(x)$ which are nonconvex in the gradient $\nabla u$ and possibly also in $u$. To characterize the microstructures for these nonconvex energies, we minimize the associated relaxed energy using two novel approaches: i) a semi-analytical method based on control systems theory, ii) and a numerical scheme that combines convex splitting together with a modified version of the split Bregman algorithm. These solutions are then used to gain information about minimizing sequences of the original problem and the spatial distribution of microstructure.Comment: 34 pages, 10 figure