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.     

Fourth‐order compact scheme with local mesh refinement for option pricing in jump‐diffusion model

The value of a contingent claim under a jump‐diffusion process satisfies a partial integro‐differential equation. A fourth‐order compact finite difference scheme is applied to discretize the spatial variable of this equation. It is discretized in time by an implicit‐explicit method. Meanwhile, a local mesh refinement strategy is used for handling the nonsmooth payoff condition. Moreover, the numerical quadrature method is exploited to evaluate the jump integral term. It guarantees a Toeplitz‐like structure of the integral operator such that a fast algorithm is feasible. Numerical results show that this approach gives fourth‐order accuracy in space. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

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     

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.     

A uniform error estimate in time for spectral Galerkin approximations of the magneto‐micropolar fluid equations

We consider Galerkin approximations for the equations modeling the motion of an incompressible magneto‐micropolar fluid in a bounded domain. We derive an optimal uniform in time error bound in the H¹ and L² ‐norms for the velocity. This is done without explicit assumption of exponential stability for a class of solutions corresponding to decaying external force fields. Our study is done for no‐slip boundary conditions, but the results obtained are easily extended to the case of periodic boundary conditions. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 689–706, 2012     

Least‐squares method for the Oseen equation

This article studies the least‐squares finite element method for the linearized, stationary Navier–Stokes equation based on the stress‐velocity‐pressure formulation in d dimensions (d = 2 or 3). The least‐squares functional is simply defined as the sum of the squares of the L² norm of the residuals. It is shown that the homogeneous least‐squares functional is elliptic and continuous in the norm. This immediately implies that the a priori error estimate of the conforming least‐squares finite element approximation is optimal in the energy norm. The L² norm error estimate for the velocity is also established through a refined duality argument. Moreover, when the right‐hand side f belongs only to , we derive an a priori error bound in a weaker norm, that is, the norm. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1289–1303, 2016     

A collocation method of lines for two‐sided space‐fractional advection‐diffusion equations with variable coefficients

We present the method of lines (MOL), which is based on the spectral collocation method, to solve space‐fractional advection‐diffusion equations (SFADEs) on a finite domain with variable coefficients. We focus on the cases in which the SFADEs consist of both left‐ and right‐sided fractional derivatives. To do so, we begin by introducing a new set of basis functions with some interesting features. The MOL, together with the spectral collocation method based on the new basis functions, are successfully applied to the SFADEs. Finally, four numerical examples, including benchmark problems and a problem with discontinuous advection and diffusion coefficients, are provided to illustrate the efficiency and exponentially accuracy of the proposed method.     

First passage time for multivariate jump‐diffusion processes in finance and other areas of applications

The first passage time (FPT) problem is an important problem with a wide range of applications in science, engineering, economics, and industry. Mathematically, such a problem can be reduced to estimating the probability of a stochastic process first to reach a boundary level. In most important applications in the financial industry, the FPT problem does not have an analytical solution and the development of efficient numerical methods becomes the only practical avenue for its solution. Most of our examples in this contribution are centered around the evaluation of default correlations in credit risk analysis, where we are concerned with the joint defaults of several correlated firms, the task that is reducible to a FPT problem. This task represents a great challenge for jump‐diffusion processes (JDP). In this contribution, we develop further our previous fast Monte Carlo method in the case of multivariate (and correlated) JDP. This generalization allows us, among other things, to evaluate the default events of several correlated assets based on a set of empirical data. The developed technique is an efficient tool for a number of financial, economic, and business applications, such as credit analysis, barrier option pricing, macroeconomic dynamics, and the evaluation of risk, as well as for a number of other areas of applications in science and engineering, where the FPT problem arises. Copyright © 2008 John Wiley & Sons, Ltd.     

Two grid finite element discretization method for semi‐linear hyperbolic integro‐differential equations

In this paper, we will investigate a two grid finite element discretization method for the semi‐linear hyperbolic integro‐differential equations by piecewise continuous finite element method. In order to deal with the semi‐linearity of the model, we use the two grid technique and derive that once the coarse and fine mesh sizes H, h satisfy the relation h = H² for the two‐step two grid discretization method, the two grid method achieves the same convergence accuracy as the ordinary finite element method. Both theoretical analysis and numerical experiments are given to verify the results.     

A uniformly optimal‐order error estimate for a bilinear finite element method for degenerate convection‐diffusion equations

We prove an optimal‐order error estimate in a degenerate‐diffusion weighted energy norm for bilinear Galerkin finite element methods for two‐dimensional time‐dependent convection‐diffusion equations with degenerate diffusion. In the estimate, the generic constants depend only on certain Sobolev norms of the true solution but not the lower bound of the diffusion. This estimate, combined with a known stability estimate of the true solution of the governing partial differential equations, yields an optimal‐order estimate of the Galerkin finite element method, in which the generic constants depend only on the Sobolev norms of the initial and right side data. Preliminary numerical experiments were conducted to verify these estimates numerically. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A Riesz–Feller space‐fractional backward diffusion problem with a time‐dependent coefficient: regularization and error estimates

In this paper, we consider a Riesz–Feller space‐fractional backward diffusion problem with a time‐dependent coefficient We show that this problem is ill‐posed; therefore, we propose a convolution regularization method to solve it. New error estimates for the regularized solution are given under a priori and a posteriori parameter choice rules, respectively. Copyright © 2016 John Wiley & Sons, Ltd.     

A nonconventional Eulerian‐Lagrangian single‐node collocation method with Hermite polynomials for unsteady‐state advection‐diffusion equations

We developed a nonconventional Eulerian‐Lagrangian single‐node collocation method (ELSCM) with piecewise‐cubic Hermite polynomials as basis functions for the numerical simulation to unsteady‐state advection‐diffusion transport partial differential equations. This method greatly reduces the number of unknowns in the conventional collocation method, and generates accurate numerical solutions even if very large time steps are taken. The method is relatively easy to formulate. Numerical experiments are presented to show the strong potential of this method. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 271–283, 2003.     

The backward euler anisotropic a posteriori error analysis for parabolic integro‐differential equations

We derive two optimal a posteriori error estimators for an implicit fully discrete approximation to the solutions of linear integro‐differential equations of the parabolic type. A continuous, piecewise linear finite element space is used for the space discretization and the time discretization is based on an implicit backward Euler method. The a posteriori error indicator corresponding to space discretization is derived using the anisotropic interpolation estimates in conjunction with a Zienkiewicz‐Zhu error estimator to approach the error gradient. The error due to time discretization is derived using continuous, piecewise linear polynomial in time. We use the linear approximation of the Volterra integral term to estimate the quadrature error in the second estimator. Numerical experiments are performed on the isotropic mesh to validate the derived results.© 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1309–1330, 2016     

A two‐grid characteristic finite volume element method for semilinear advection‐dominated diffusion equations

A two‐grid finite volume element method, combined with the modified method of characteristics, is presented and analyzed for semilinear time‐dependent advection‐dominated diffusion equations in two space dimensions. The solution of a nonlinear system on the fine‐grid space (with grid size h) is reduced to the solution of two small (one linear and one nonlinear) systems on the coarse‐grid space (with grid size H) and a linear system on the fine‐grid space. An optimal error estimate in H¹ ‐norm is obtained for the two‐grid method. It shows that the two‐grid method achieves asymptotically optimal approximation, as long as the mesh sizes satisfy h = O(H²). Numerical example is presented to validate the usefulness and efficiency of the method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Moving collocation method for a reaction‐diffusion equation with a traveling heat source

We consider a reaction‐diffusion equation with a traveling heat source on an unbounded domain. The numerical simulation of the problem is difficult because of the moving singularity, the blow‐up phenomenon, and the delta function in the equation. Because we are only interested in the solution behavior near the heat source, we choose a bounded moving domain which contains the heat source and has the same speed as the source. Local absorbing boundary conditions are constructed on the boundaries of the moving domain. Then, we transform the moving domain to a fixed one. At last, a special moving collocation method is adopted. The new method is much simpler than the existing moving finite difference methods. Moreover, numerical experiments illustrate the accuracy and efficiency of our moving collocation method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Combination of standard Galerkin and subspace methods for the time‐dependent Navier‐Stokes equations with nonsmooth initial data

We consider a combination of the standard Galerkin method and the subspace decomposition methods for the numerical solution of the two‐dimensional time‐dependent incompressible Navier‐Stokes equations with nonsmooth initial data. Because of the poor smoothness of the solution near t = 0, we use the standard Galerkin method for time interval [0, 1] and the subspace decomposition method time interval [1, ∞). The subspace decomposition method is based on the solution into the sum of a low frequency component integrated using a small time step Δt and a high frequency integrated using a larger time step pΔt with p > 1. From the H¹‐stability and L²‐error analysis, we show that the subspace decomposition method can yield a significant gain in computing time. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2009     

A priori error estimates of an extrapolated space‐time discontinuous galerkin method for nonlinear convection‐diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection‐diffusion equation. We employ a combination of the discontinuous Galerkin finite element (DGFE) method for the space as well as time discretization. The linear diffusive and penalty terms are treated implicitly whereas the nonlinear convective term is treated by a special higher order explicit extrapolation from the previous time step, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the L^∞(L²) –norm and the L²(H¹) –seminorm with respect to the mesh size h and time step τ. Finally, we present an efficient solution strategy and numerical examples verifying the theoretical results. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1456–1482, 2010     

A posteriori error estimates for the generalized Schwarz method of a new class of advection‐diffusion equation with mixed boundary condition

In this paper, a posteriori error estimates for the generalized Schwarz method with mixed boundary condition on the interfaces for advection‐diffusion equation with second‐order boundary value problems are proved using theta time scheme combined with Galerkin spatial method. Furthermore, a asymptotic behavior in Sobolev norm is deduced using Benssoussan‐Lions' algorithm.     

Stability and error analysis for a spectral Galerkin method for the Navier‐Stokes equations with H2 or H1 initial data

In this article we consider a spectral Galerkin method with a semi‐implicit Euler scheme for the two‐dimensional Navier‐Stokes equations with H² or H¹ initial data. The H²‐stability analysis of this spectral Galerkin method shows that for the smooth initial data the semi‐implicit Euler scheme admits a large time step. The L²‐error analysis of the spectral Galerkin method shows that for the smoother initial data the numerical solution u exhibits faster convergence on the time interval [0, 1] and retains the same convergence rate on the time interval [1, ∞). © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

An application of exp‐function method to approximate general and explicit solutions for nonlinear Schrödinger equations

In this work, we implement a relatively new analytical technique, the exp‐function method, for solving nonlinear equations and absolutely a special form of generalized nonlinear Schrödinger equations, which may contain high‐nonlinear terms. This method can be used as an alternative to obtain analytical and approximate solutions of different types of fractional differential equations, which is applied in engineering mathematics. Some numerical examples are presented to illustrate the efficiency and the reliability of exp‐function method. It is predicted that exp‐function method can be found widely applicable in engineering. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1016–1025, 2011