CI²MA - Publications | Preprints

Preprint 2014-12

Sebastiano Boscarino, Raimund Bürger, Pep Mulet, Giovanni Russo, Luis M. Villada:

Linearly implicit IMEX Runge-Kutta methods for a class of degenerate convection-diffusion problems


Multi-species kinematic flow models with strongly degenerate diffusive corrections give rise to systems of nonlinear convection-diffusion equations of arbitrary size. Applications of these systems include models of polydisperse sedimentation and multi-class traffic flow. Implicit-explicit (IMEX) Runge-Kutta (RK) methods are suitable for the solution of these convection-diffusion problems since the stability restrictions, coming from the explicitly treated convective part, are much less severe than those that would be deduced from an explicit treatment of the diffusive term. These schemes usually combine an explicit Runge-Kutta scheme for the time integration of the convective part with a diagonally implicit one for the diffusive part. In [R. Bürger, P. Mulet and L.M. Villada, SIAM J. Sci. Comput., 35 (2013), pp. B751–B777] a scheme of this type is proposed, where the nonlinear and non-smooth systems of algebraic equations arising in the implicit treatment of the degenerate diffusive part are solved by smoothing of the diffusion coefficients combined with a Newton-Raphson method with line search. This nonlinearly implicit method is robust but associated with considerable effort of implementation and possibly CPU time. To overcome these shortcomings while keeping the advantageous stability properties of IMEX-RK methods, a second variant of these methods is proposed, in which the diffusion terms are discretized in a way that more carefully distinguishes between stiff and nonstiff dependence, such that in each time step only a linear system needs to be solved still maintaining high order accuracy in time, which makes these methods much simpler to implement. In a series of examples of polydisperse sedimentation and multi-class traffic flow it is demonstrated that these new linearly implicit IMEX-RK schemes approximate the same solutions as the nonlinearly implicit versions, and in many cases these schemes are more efficient.

Download in PDF format PDF

This preprint gave rise to the following definitive publication(s):

Sebastiano BOSCARINO, Raimund BüRGER, Pep MULET, Giovanni RUSSO, Luis M. VILLADA: Linearly implicit IMEX Runge-Kutta methods for a class of degenerate convection-diffusion problems. SIAM Journal on Scientific Computing, vol. 37, no. 2, pp. B305-B331, (2015).