In order to model the impact of incidents an alternative resolution method is proposed for the Lighthill-Whitham-Richards model. Is is based on an explicit handling of shock waves (generation, tracking and collision). Contrary to finite difference methods which impose a fixed discretization grid, the Wave Tracking method is event based. The different possible events are studied and they all lea...

In this paper we establish the stability condition of a general class of finite difference schemes applied to nonlinear complex reaction-diffusion equations. We consider the numerical solution of both implicit and semi-implicit discretizations. To illustrate the theoretical results we present some numerical examples computed with a semi-implicit scheme applied to a nonlinear equation.

We analyze the discretization of initial and boundary value problems with a stationary interface in one space dimension for the heat equation, the Schrödinger equation, and the wave equation by finite difference methods. Extending the concept of the elliptic projection, well known from the analysis of Galerkin finite element methods, to our finite difference case, we prove second-order error es...

Bounds are given for the error constant of stable finite-difference methods for first-order hyperbolic equations in one space dimension, which use r downwind and s upwind points in the discretization of the space derivatives, and which are of optimal order p = min(r + s,2r + 2,2s). It is known that this order can be obtained by interpolatory methods. Examples show, however, that their error con...

An adaptive finite difference method for first order nonlinear systems of ordinary differential equations subject to multipoint nonlinear boundary conditions is presented. The method is based on a discretization studied earlier by H. B. Keller. Variable order is provided through deferred corrections, while a built-in natural asymptotic estimator is used to automatically refine the mesh in order...

This paper examines formulation of the discretization schemes for diffusion process simulation that allow coarse grid spacings in the areas of exponentially varying concentrations and fluxes. The method of integral identities is used as a common framework for exponential fitting of both the finite difference and finite element schemes. An exponentially fitted finite difference scheme, with disc...

  In this article, vibration analysis of an Euler-Bernoulli beam resting on a Pasternak-type foundation is studied. The governing equation is solved by using a spectral finite element model (SFEM). The solution involves calculating wave and time responses of the beam. The Fast Fourier Transform function is used for temporal discretization of the governing partial differential equation into a se...

In this paper, we derive a discretized multi-group epidemic model with time delay by using a nonstandard finite difference (NSFD) scheme. A crucial observation regarding the advantage of the NSFD scheme is that the positivity and boundedness of solutions of the continuous model are preserved. Furthermore, we show that the discrete model has the same equilibria, and the conditions for their stab...

in this paper, a new numerical method for solving time-fractional diffusion equations is introduced. for this purpose, finite difference scheme for discretization in time and chebyshev collocation method is applied. also, to simplify application of the method, the matrix form of the suggested method is obtained. illustrative examples show that the proposed method is very efficient and accurate.

In this paper, we derive a highly accurate numerical method for the solution of one-dimensional wave equation with Neumann boundary conditions. This hyperbolic problem is solved by using semidiscrete approximations. The space direction is discretized by wavelet-Galerkin method and the time variable is discretized by using various classical finite difference schemes. The numerical results show t...