On linearly implicit IMEX Runge-Kutta methods for degenerate convection-diffusion problems modeling polydisperse sedimentation

Implicit-explicit (IMEX) Runge-Kutta (RK) methods are suitable for the solution of nonlinear, possibly strongly degenerate, 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. A particularly efficient variant of these schemes, so-called linearly implicit IMEX-RK schemes, arise from discretizing the diffusion terms 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. These schemes provide an efficient tool for the numerical exploration of sediment formation and composition under a strongly degenerate polydisperse sedimentation model.     

An embedded pair of exponentially fitted explicit Runge–Kutta methods

An embedded pair of exponentially fitted explicit Runge–Kutta (RK) methods for the numerical integration of IVPs with oscillatory solutions is derived. This pair is based on the exponentially fitted explicit RK method constructed in Vanden Berghe et al., and we confirm that the methods which constitute the pair have algebraic order 4 and 3. Some numerical experiments show the efficiency of our pair when it is compared with the variable step code proposed by Vanden Berghe et al. (J. Comput. Appl. Math. 125 (2000) 107).     

Linearly Implicit Approximations of Diffusive Relaxation Systems

Diffusive relaxation systems provide a general framework to approximate nonlinear diffusion problems, also in the degenerate case (Aregba-Driollet et al. in Math. Comput. 73(245):63–94, 2004; Boscarino et al. in Implicit-explicit Runge-Kutta schemes for hyperbolic systems and kinetic equations in the diffusion limit, 2011; Cavalli et al. in SIAM J. Sci. Comput. 34:A137–A160, 2012; SIAM J. Numer. Anal. 45(5):2098–2119, 2007; Naldi and Pareschi in SIAM J. Numer. Anal. 37:1246–1270, 2000; Naldi et al. in Surveys Math. Indust. 10(4):315–343, 2002). Their discretization is usually obtained by explicit schemes in time coupled with a suitable method in space, which inherits the standard stability parabolic constraint. In this paper we combine the effectiveness of the relaxation systems with the computational efficiency and robustness of the implicit approximations, avoiding the need to resolve nonlinear problems and avoiding stability constraints on time step. In particular we consider an implicit scheme for the whole relaxation system except for the nonlinear source term, which is treated though a suitable linearization technique. We give some theoretical stability results in a particular case of linearization and we provide insight on the general case. Several numerical simulations confirm the theoretical results and give evidence of the stability and convergence also in the case of nonlinear degenerate diffusion.     

On entropy stable schemes for degenerate parabolic multispecies kinematic flow models

Entropy stable schemes for the numerical solution of initial value problems of nonlinear, possibly strongly degenerate systems of convection–diffusion equations were recently proposed in Jerez and Parés's study. These schemes extend the theoretical framework of Tadmor's study to convection–diffusion systems. They arise from entropy conservative schemes by adding a small amount of viscosity to avoid spurious oscillations. The main condition for feasibility of entropy conservative or stable schemes for a given model is that the corresponding first‐order system of conservation laws possesses a convex entropy function and corresponding entropy flux, and that the diffusion matrix multiplied by the inverse of the Hessian of the entropy is positive semidefinite. As a new contribution, it is demonstrated in the present work, first, that these schemes can naturally be extended to initial‐boundary value problems with zero‐flux boundary conditions in one space dimension, including an explicit bound on the growth of the total entropy. Second, it is shown that these assumptions are satisfied by certain diffusively corrected multiclass kinematic flow models of arbitrary size that describe traffic flow or the settling of dispersions and emulsions, where the latter application gives rise to zero‐flux boundary conditions. Numerical examples illustrate the behavior and accuracy of entropy stable schemes for these applications.     

Second-order schemes for conservation laws with discontinuous flux modelling clarifier–thickener units

Continuously operated clarifier–thickener (CT) units can be modeled by a non-linear, scalar conservation law with a flux that involves two parameters that depend discontinuously on the space variable. This paper presents two numerical schemes for the solution of this equation that have formal second-order accuracy in both the time and space variable. One of the schemes is based on standard total variation diminishing (TVD) methods, and is addressed as a simple TVD (STVD) scheme, while the other scheme, the so-called flux-TVD (FTVD) scheme, is based on the property that due to the presence of the discontinuous parameters, the flux of the solution (rather than the solution itself) has the TVD property. The FTVD property is enforced by a new nonlocal limiter algorithm. We prove that the FTVD scheme converges to a BV _t solution of the conservation law with discontinuous flux. Numerical examples for both resulting schemes are presented. They produce comparable numerical errors, while the FTVD scheme is supported by convergence analysis. The accuracy of both schemes is superior to that of the monotone first-order scheme based on the adaptation of the Engquist–Osher scheme to the discontinuous flux setting of the CT model (Bürger, Karlsen and Towers in SIAM J Appl Math 65:882–940, 2005). In the CT application there is interest in modelling sediment compressibility by an additional strongly degenerate diffusion term. Second-order schemes for this extended equation are obtained by combining either the STVD or the FTVD scheme with a Crank–Nicolson discretization of the degenerate diffusion term in a Strang-type operator splitting procedure. Numerical examples illustrate the resulting schemes.     

The fast scalar auxiliary variable approach with unconditional energy stability for nonlocal Cahn–Hilliard equation

Comparing with the classical local gradient flow and phase field models, the nonlocal models such as nonlocal Cahn–Hilliard equations equipped with nonlocal diffusion operator can describe more practical phenomena for modeling phase transitions. In this paper, we construct an accurate and efficient scalar auxiliary variable approach for the nonlocal Cahn–Hilliard equation with general nonlinear potential. The first contribution is that we have proved the unconditional energy stability for nonlocal Cahn–Hilliard model and its semi‐discrete schemes carefully and rigorously. Second, what we need to focus on is that the nonlocality of the nonlocal diffusion term will lead the stiffness matrix to be almost full matrix which generates huge computational work and memory requirement. For spatial discretizaion by finite difference method, we find that the discretizaition for nonlocal operator will lead to a block‐Toeplitz–Toeplitz‐block matrix by applying four transformation operators. Based on this special structure, we present a fast procedure to reduce the computational work and memory requirement. Finally, several numerical simulations are demonstrated to verify the accuracy and efficiency of our proposed schemes.     

An IMEX‐BDF2 compact scheme for pricing options under regime‐switching jump‐diffusion models

In this paper, an implicit‐explicit two‐step backward differentiation formula (IMEX‐BDF2) together with finite difference compact scheme is developed for the numerical pricing of European and American options whose asset price dynamics follow the regime‐switching jump‐diffusion process. It is shown that IMEX‐BDF2 method for solving this system of coupled partial integro‐differential equations is stable with the second‐order accuracy in time. On the basis of IMEX‐BDF2 time semi‐discrete method, we derive a fourth‐order compact (FOC) finite difference scheme for spatial discretization. Since the payoff function of the option at the strike price is not differentiable, the results show only second‐order accuracy in space. To remedy this, a local mesh refinement strategy is used near the strike price so that the accuracy achieves fourth order. Numerical results illustrate the effectiveness of the proposed method for European and American options under regime‐switching jump‐diffusion models.     

Nonstandard discretizations of the generalized Nagumo reaction‐diffusion equation

A competitive nonstandard semi‐explicit finite‐difference method is constructed and used to obtain numerical solutions of the diffusion‐free generalized Nagumo equation. Qualitative stability analysis and numerical simulations show that this scheme is more robust in comparison to some standard explicit methods such as forward Euler and the fourth‐order Runge‐Kutta method (RK4). The nonstandard scheme is extended to construct a semi‐explicit and an implicit scheme to solve the full Nagumo reaction‐diffusion equation. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 363–379, 2003.     

A large time‐step implicit moving mesh scheme for moving boundary problems

Velocity‐based moving mesh methods update the mesh at each time level by using a velocity equation with a time‐stepping scheme. A particular velocity‐based moving mesh method, based on conservation, uses explicit time‐stepping schemes with small time steps to avoid mesh tangling. However, this can prove to be impractical when long‐term behavior of the solution is of interest. Here, we present a semi‐implicit time‐stepping scheme which manipulates the structure of the velocity equation such that it resembles a variable‐coefficient heat equation. This enables the use of maximum/minimum principle which ensures that mesh tangling is avoided. It is also shown that this semi‐implicit scheme can be extended to a fully implicit time‐stepping scheme. Thus, the time‐step restriction imposed by explicit schemes is overcome without sacrificing mesh structure. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 321–338, 2014     

Robust spectral method for numerical valuation of european options under Merton's jump‐diffusion model

We propose a novel numerical method based on rational spectral collocation and Clenshaw–Curtis quadrature methods together with the “” transformation for pricing European vanilla and butterfly spread options under Merton's jump‐diffusion model. Under certain assumptions, such model leads to a partial integro‐differential equation (PIDE). The differential and integral parts of the PIDE are approximated by the rational spectral collocation and the Clenshaw–Curtis quadrature methods, respectively. The application of spectral collocation method to the PIDE leads to a system of ordinary differential equations, which is solved using the implicit–explicit predictor–corrector (IMEX‐PC) schemes in which the diffusion term is integrated implicitly, whereas the convolution integral, reaction, advection terms are integrated explicitly. Numerical experiments illustrate that our approach is highly accurate and efficient for pricing financial options.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1169–1188, 2014     

On a Parabolic Sine-Gordon Model

We consider a parabolic sine-Gordon model with periodic boundary conditions. We prove a fundamental maximum principle which gives a priori uniform control of the solution. In the one-dimensional case we classify all bounded steady states and exhibit some explicit solutions. For the numerical discretization we employ first order IMEX, and second order BDF2 discretization without any additional stabilization term. We rigorously prove the energy stability of the numerical schemes under nearly sharp and quite mild time step constraints. We demonstrate the striking similarity of the parabolic sine-Gordon model with the standard Allen-Cahn equations with double well potentials.     

On convergence of finite volume schemes for one-dimensional two-phase flow in porous media

This paper deals with development and analysis of finite volume schemes for a one-dimensional nonlinear, degenerate, convection-diffusion equation having application in petroleum reservoir and groundwater aquifer simulation. The main difficulty is that the solution typically lacks regularity due to the degenerate nonlinear diffusion term. We analyze and compare three families of numerical schemes corresponding to explicit, semi-implicit, and implicit discretization of the diffusion term and a Godunov scheme for the advection term. L^∞ stability under appropriate CFL conditions and BV estimates are obtained. It is shown that the schemes satisfy a discrete maximum principle. Then we prove convergence of the approximate solution to the weak solution of the problem. Results of numerical experiments using the present approach are reported.     

Asset pricing for an affine jump‐diffusion model using an FD method of lines on nonuniform meshes

We present a novel numerical scheme for the valuation of options under a well‐known jump‐diffusion model. European option pricing for such a case satisfies a 1 + 2 partial integro‐differential equation (PIDE) including a double integral term, which is nonlocal. The proposed approach relies on nonuniform meshes with a focus on the discontinuous and degenerate areas of the model and applying quadratically convergent finite difference (FD) discretizations via the method of lines (MOL). A condition for observing the time stability of the fully discretized problem is given. Also, we report results of numerical experiments.     

Error estimates for a mixed finite element discretization of some degenerate parabolic equations

We consider a numerical scheme for a class of degenerate parabolic equations, including both slow and fast diffusion cases. A particular example in this sense is the Richards equation modeling the flow in porous media. The numerical scheme is based on the mixed finite element method (MFEM) in space, and is of one step implicit in time. The lowest order Raviart–Thomas elements are used. Here we extend the results in Radu et al. (SIAM J Numer Anal 42:1452–1478, 2004), Schneid et al. (Numer Math 98:353–370, 2004) to a more general framework, by allowing for both types of degeneracies. We derive error estimates in terms of the discretization parameters and show the convergence of the scheme. The features of the MFEM, especially of the lowest order Raviart–Thomas elements, are now fully exploited in the proof of the convergence. The paper is concluded by numerical examples.     

Polynomial spectral collocation method for space fractional advection–diffusion equation

This article discusses the spectral collocation method for numerically solving nonlocal problems: one‐dimensional space fractional advection–diffusion equation; and two‐dimensional linear/nonlinear space fractional advection–diffusion equation. The differentiation matrixes of the left and right Riemann–Liouville and Caputo fractional derivatives are derived for any collocation points within any given bounded interval. Several numerical examples with different boundary conditions are computed to verify the efficiency of the numerical schemes and confirm the exponential convergence; the physical simulations for Lévy–Feller advection–diffusion equation and space fractional Fokker–Planck equation with initial δ‐peak and reflecting boundary conditions are performed; and the eigenvalue distributions of the iterative matrix for a variety of systems are displayed to illustrate the stabilities of the numerical schemes in more general cases. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 514–535, 2014     

Arbitrarily high-order accurate and energy-stable schemes for solving the conservative Allen–Cahn equation

In this paper, three high-order accurate and unconditionally energy-stable methods are proposed for solving the conservative Allen–Cahn equation with a space–time dependent Lagrange multiplier. One is developed based on an energy linearization Runge–Kutta (EL–RK) method which combines an energy linearization technique with a specific class of RK schemes, the other two are based on the Hamiltonian boundary value method (HBVM) including a Gauss collocation method, which is the particular instance of HBVM, and a general class of cases. The system is first discretized in time by these methods in which the property of unconditional energy stability is proved. Then the Fourier pseudo-spectral method is employed in space along with the proofs of mass conservation. To show the stability and validity of the obtained schemes, a number of 2D and 3D numerical simulations are presented for accurately calculating geometric features of the system. In addition, our numerical results are compared with other known structure-preserving methods in terms of numerical accuracy and conservation properties.     

IMEX Runge–Kutta schemes for reaction–diffusion equations

A fundamental research is carried out into convergence and stability properties of IMEX (implicit–explicit) Runge–Kutta schemes applied to reaction–diffusion equations. It is shown that a fully discrete scheme converges if it satisfies certain conditions using a technique of the B-convergence analysis, developed by Burrage, Hundsdorfer and Verwer in 1986. Stability of the schemes is also examined on the basis of a scalar test equation, proposed by Frank, Hundsdorfer and Verwer in 1997.     

An asymptotically stable Particle-in-Cell (PIC) scheme for collisionless plasma simulations near quasineutrality

We propose a new Particle-in-Cell scheme for the Vlasov–Poisson equation. This scheme remains stable when the Debye length and plasma period tend to zero without any restriction on the size of the time and length step. It relies on a semi-implicit integration of the particle trajectories. The numerical integration cost is that of the standard explicit method thanks to the use of a reformulation of the Poisson equation. To cite this article: P. Degond et al., C. R. Acad. Sci. Paris, Ser. I 343 (2006).     

Numerical Simulation of Two-Phase Flow through Heterogeneous Porous Media

A mixed finite element method is combined to finite volume schemes on structured and unstructured grids for the approximation of the solution of incompressible flow in heterogeneous porous media. A series of numerical examples demonstrates the effectiveness of the methodology for a coupled system which includes an elliptic equation and a nonlinear degenerate diffusion–convection equation arising in modeling of flow and transport in porous media.     

Convergence Analysis of a Numerical Scheme for the Porous Medium Equation by an Energetic Variational Approach

The porous medium equation (PME)is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al., J. Comput. Phys., 385 (2019), pp. 13–32], where the trajectory equation can be obtained and two numerical schemes have been developed based on different dissipative energy laws. It is also proved that the nonlinear scheme, based on $f$log $f$ as the total energy form of the dissipative law, is uniquely solvable on an admissible convex set and preserves the corresponding discrete dissipation law. Moreover, under certain smoothness assumption, we have also obtained the second order convergence in space and the first order convergence in time for the scheme. In this paper, we provide a rigorous proof of the error estimate by a careful higher order asymptotic expansion and two step error estimates. The latter technique contains a rough estimate to control the highly nonlinear term in a discrete $W$^1,∞norm and a refined estimate is applied to derive the optimal error order.