HIGH-ORDER ASYMPTOTIC-PRESERVING METHODS FOR FULLY NONLINEAR RELAXATION PROBLEMS

We study solutions to nonlinear hyperbolic systems with fully nonlinear relaxation terms in the limit of, both, infinitely stiff relaxation and arbitrary late time. In this limit, the dynamics is governed by effective systems of parabolic-type which may contain degenerate and/or fully nonlinear diffusion terms. For this class of problems, we develop an implicit-explicit method based on Runge-Kutta discretization in time, and we apply this method to the investigation of several examples of interest in fluid dynamics. Importantly, we impose here a realistic stability condition on the time step and we demonstrate that solutions in the hyperbolic-to-parabolic regime can be computed numerically with high robustness and accuracy, even in the presence of fully nonlinear relaxation.