A new mimetic finite difference scheme for solving the acoustic wave equation is presented. It combines a novel second order tensor mimetic discretizations in space and a leapfrog approximation in time to produce an explicit multidimensional scheme. Convergence analysis of the new scheme on a staggered grid shows that it can take larger time steps than standard finite difference schemes based o...

This paper presents a 3D model for pricing defaultable bonds with embedded call options. The pricing model incorporates three essential ingredients in the pricing of defaultable bonds: stochastic interest rate, stochastic default risk, and call provision. Both the stochastic interest rate and the stochastic default risk are modeled as a square-root diffusion process. The default risk process is...

We establish the global existence of L∞ solutions for a model of polytropic gas flow with diffusive entropy. The result is obtained by showing the convergence of a class of finite difference schemes, which includes the Lax– Friedrichs and Godunov schemes. Such convergence is achieved by proving the estimates required for the application of the compensated compactness theory.

A class of finite difference schemes for solving a fractional anti-diffusive equation, recently proposed by Andrew C. Fowler to describe the dynamics of dunes, is considered. Their linear stability is analyzed using the standard Von Neumann analysis: stability criteria are found and checked numerically. Moreover, we investigate the consistency and convergence of these schemes.

Several numerical schemes utilizing central difference approximations have been developed to solve the Goursat problem. However, in a recent years compact discretization methods which leads to high-order finite difference schemes have been used since it is capable of achieving better accuracy as well as preserving certain features of the equation e.g. linearity. The basic idea of the new scheme...

This article is devoted to the numerical study of various finite difference approximations to the stochastic Burgers equation. Of particular interest in the one-dimensional case is the situation where the driving noise is white both in space and in time. We demonstrate that in this case, different finite difference schemes converge to different limiting processes as the mesh size tends to zero....

11

In this paper we propose a new method based on adjoint systems for parametric sensitivity analysis of DAE’s. This method is employed in a series of experiments, and the results are compared with the estimation of sensitivity by a finite differences method, widely used because of the simplicity of its implementation. In addition, we propose a method extension based on the use of estimation funct...

This article presents a Monte Carlo methodology for solving the stochastic optimal timber harvest problem modeled as a recurrent American call option. A detailed description of the proposed methodology is given, and the Monte Carlo technique is contrasted with finite difference methods typically used to find solutions of the optimal harvest problem with stochastic prices. The use of the methodo...

A system of singularly perturbed ordinary differential equations of first order with given initial conditions is considered. The leading term of each equation is multiplied by a small positive parameter. These parameters are assumed to be distinct and they determine the different scales in the solution to this problem. A Shishkin piecewise–uniform mesh is constructed, which is used, in conjunct...