Article history: Received 27 June 2013 Received in revised form 22 January 2014 Accepted 25 March 2014 Available online 2 April 2014

Fractional-order diffusion equations are viewed as generalizations of classical diffusion equations, treating super-diffusive flow processes. In this paper, in order to solve the fractional advection-diffusion equation, the fractional characteristic finite difference method is presented, which is based on the method of characteristics (MOC) and fractional finite difference (FD) procedures. The ...

Abslrczc-The 3-dimensional finite-difference time-domain method in non-orthogonal co-ordinates (non-standard FDTD) is used to calculate the frequencies of resonators. The numerical boundary conditions of the method are presented. The influences of boundary conditions and discrete meshes on the numerical accuracy are investigated. We present the nonstandard FDTD method using the boundary-orthogo...

Using classic differential quadrature formulae and uniform grids, this paper systematically constructs a variety of high-order finite difference schemes, and some of these schemes are consistent with the so-called boundary value methods. The derived difference schemes enjoy the same stability and accuracy properties with correspondent differential quadrature methods but have a simpler form of c...

High-order finite difference methods are efficient, easy to program, scale well in multiple dimensions and can be modified locally for various reasons (such as shock treatment for example). The main drawback has been the complicated and sometimes even mysterious stability treatment at boundaries and interfaces required for a stable scheme. The research on summation-byparts operators and weak bo...

and Applied Analysis 3 The grid function y(x, t) is a function defined at the grid points of g. we denote the nodal values of a grid function y(x, t) between time levels t 0 and t 0 as y (x, t) = y (x 1 , x 2 , t l,j i ) = y l,j n1 ,n2 , (11) for x ∈ ω i , i > 0, j = 0, . . . , m i . For x ∈ ω 0 we define y (x, t) = y (x 1 , x 2 , t l+1 0 ) = y l+1 n1 ,n2 . (12) δ x1 , δ x1 and δ x2 , δ x2 are ...

We formulate and analyze an unconditionally stable nonstandard finite difference method for a mathematical model of HIV transmission dynamics. The dynamics of this model are studied using the qualitative theory of dynamical systems. These qualitative features of the continuous model are preserved by the numerical method that we propose in this paper. This method also preserves the positivity of...

When one solves differential equations, modeling biological or physical phenomena, it is of great importance to take physical constraints into account. More precisely, numerical schemes have to be designed such that discrete solutions satisfy the same constraints as exact solutions. In this work, we introduce explicit finite difference schemes based on the nonstandard discretization method to a...