Journal of Applied Mathematics

Volume 2013 (2013), Article ID 586870, 7 pages

http://dx.doi.org/10.1155/2013/586870

## Numerical Simulations for the Space-Time Variable Order Nonlinear Fractional Wave Equation

^{1}Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt^{2}Department of Mathematics, Faculty of Science, Umm Al-Qura University, Makkah 21955, Saudi Arabia

Received 3 March 2013; Revised 29 April 2013; Accepted 3 May 2013

Academic Editor: Chein-Shan Liu

Copyright © 2013 Nasser Hassan Sweilam and Taghreed Abdul Rahman Assiri. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The explicit finite-difference method for solving variable order fractional space-time wave equation with a nonlinear source term is considered. The concept of variable order fractional derivative is considered in the sense of Caputo. The stability analysis and the truncation error of the method are discussed. To demonstrate the effectiveness of the method, some numerical test examples are presented.

#### 1. Introduction

It is well known that the fractional calculus definitions are extensions of the usual calculus definitions [1–8], where the orders need not to be positive integers. On the other hand, the variable order calculus is a natural extension of the constant order (integer or fractional) calculus. In this sense, the order may function in any variable such as time and space variables or a system of other variables [9, 10]. In general, one can say that this extension is introduced by Samko and Ross in [11], where Marchaud fractional derivative and Riemann-Liouville derivative are extended to the variable order cases the order in this case is a function in the space variable only. Many authors have introduced different definitions of variable order differential operators, each of these with a specific meaning to suit desired goals. These definitions such as Riemann-Liouville, Grünwald, Caputo, Riesz [3, 12–16], and some notes as Coimbra definition [17, 18].

Coimbra in [17] used Laplace transform of Caputo’s definition of the fractional derivative as the starting point to suggest a novel definition for the variable order differential operator. Because of its meaningful physical interpretation, Coimbra’s definition is better suited for modeling physical problems. The variable order differentials are an important tool to study some systems such as the control of nonlinear viscoelasticity oscillator (for more details see [17–19] and the references cited therein), where the order changes with respect to a parameter or more parameters.

In the following, we present the basic definition for the variable order fractional derivatives which we will use in this paper.

*Definition 1 (see [14]). *The Caputo space variable order derivative is defined as follows:
where .

The main aim of this work is to use the explicit finite difference method (EFDM) to study numerically the following nonlinear space-time variable order wave equation: subject to initial conditions and the following boundary conditions where is a constant, are smooth functions, and is a nonlinear scour term that satisfies the Lipschitz condition, that is, where the constant is called a Lipschitz constant for .

#### 2. Discretization for EFDM

In this section, EFDM is used to study the model problem (2), then the space-time solutions domain will be discretized. The discrete form for the pervious Caputo derivative can be written as follows: Then, Now, pick two positive integers , and define the step size of space and time by , , respectively, where and . Also we introduce the following notations: , , and . Then, By the same way, we have For simplicity, let us define Then, we can rewrite (2) in the following form: that is, The previous equation can be expressed in the following matrix form: and for where ,, and is a matrix with the following coefficients: for , and . Also, we note that then .

Lemma 2. *The coefficients and satisfy the following conditions:*(1)*, and ,*(2)*, and , for *

#### 3. The Stability Analysis and the Truncation Error

Let us consider and to be two different numerical solutions of (15) with initial values given by and , respectively.

Theorem 3. *The explicit method approximation defined by (15) to the variable order space-time wave equation (2) is unconditionally stable, that is,
*

*Proof. *Let us define .

From (15) we have
where
and .

Noting that , for any .

Let . From (20), we have , then
Now, we analyze the stability via mathematical induction [10]. From (14) we have , where is a constant.

Now, assume that , then from (22), we have
Then, the theorem holds.

Lemma 4. *Let
**
be a smooth function; then
*

*Proof. *In terms of standard centered difference formula, we have
By the integral mean value theorem, we have
where . Combining the pervious two formulae, we have
Now, by using Lemma 4, we can derive the truncation error of explicit finite difference scheme (14). It has a local truncation error of (from the left side) and (from the right side).

*Remark 5. *The pervious explicit method was shown to be stable. This method is consistent with a local truncation error which is . Therefore, according to the Lax Equivalence Theorem [2], it converges at this rate.

#### 4. Numerical Examples

*Example 1. *Consider the following variable-order linear fractional wave equation:
with , and
subjected to the following initial conditions:
where , , and .

The exact solution of this problem when is .

In Figure 1, a comparison between the numerical and the exact solutions when at is presented.

In Figures 2(a) and 2(b), we report the approximate solutions at and , respectively.

In Figures 3, 4, and 5, respectively, we report the approximate solutions in three dimensions, where the axis’s are , and , respectively.

*Example 2. *Consider the following variable-order nonlinear fractional wave equation:
with , and ,, where , and .

This problem has the following exact solution, when

In Figure 6(a), we report the numerical solution when , are variables at and the exact solutions when .

In Figures 6(b) and 6(c), we report the approximate solution at and , respectively.

Figures 7, 8, and 9 show the approximate solution in three dimensions, where the axes are , and , respectively.

#### 5. Conclusions

In this paper, numerical studies using a simple explicit FDM for solving the variable order space-time wave equation are presented. The stability analysis and the truncation error of the proposed method are proved. Some test examples are given, and the results obtained by the method are compared with the exact solutions in integer order cases. Several figures are presented to simulate the solutions behaviors when the variable orders change with respect to space and time. The comparison certifies that FDM gives good results. Summarizing these results, we can say that the finite difference method in its general form gives reasonable calculations, easy to use, and can be applied for the variable order differential equations in general form. All results were obtained by using MATLAB version 7.6.0 (R2008a).

#### References

- I. Podlubny,
*Fractional Differential Equations*, Academic Press, 1999. View at MathSciNet - R. D. Richtmyer and K. W. Morton,
*Difference Methods for Initial-Vlue Problems*, Krieger, Malabar, Fla, USA, 2nd edition, 1994. View at MathSciNet - R. Lin, F. Liu, V. Anh, and I. Turner, “Stability and convergence of a new explicit finite-difference approximation for the variable-order nonlinear fractional diffusion equation,”
*Applied Mathematics and Computation*, vol. 212, no. 2, pp. 435–445, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. F. Lorenzo and T. T. Hartley, “A Solution to the fundamental linear fractional order differential equation,”
*NASA/TP*1998-208693, 1998. View at Google Scholar - C. F. Lorenzo and T. T. Hartley, “The vector linear fractional initialization problem,”
*NASA/TP*1999-208919, 1999. View at Google Scholar - N. H. Sweilam, M. M. Khader, and R. F. Al-Bar, “Numerical studies for a multi-order fractional differential equation,”
*Physics Letters A*, vol. 371, no. 1-2, pp. 26–33, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. H. Sweilam and M. M. Khader, “A Chebyshev pseudo-spectral method for solving fractional-order integro-differential equations,”
*The ANZIAM Journal*, vol. 51, no. 4, pp. 464–475, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. H. Sweilam, M. M. Khader, and A. M. Nagy, “Numerical solution of two-sided space-fractional wave equation using finite difference method,”
*Journal of Computational and Applied Mathematics*, vol. 235, no. 8, pp. 2832–2841, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. H. Sweilam, M. M. Khader, and H. M. Almarwm, “Numerical studies for the variable-order nonlinear fractional wave equation,”
*Fractional Calculus and Applied Analysis*, vol. 15, no. 4, pp. 669–683, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - S. Shen and F. Liu, “Error analysis of an explicit finite difference approximation for the space fractional diffusion equation with insulated ends,”
*The ANZIAM Journal*, vol. 46, pp. C871–C887, 2005. View at Google Scholar · View at MathSciNet - S. G. Samko and B. Ross, “Integration and differentiation to a variable fractional order,”
*Integral Transforms and Special Functions*, vol. 1, no. 4, pp. 277–300, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Ingman, J. Suzdalnitsky, and M. Zeifman, “Constitutive dynamic-order model for nonlinear contact phenomena,”
*Journal of Applied Mechanics, Transactions ASME*, vol. 67, no. 2, pp. 383–390, 2000. View at Publisher · View at Google Scholar · View at Scopus - D. Ingman and J. Suzdalnitsky, “Control of damping oscillations by fractional differential operator with time-dependent order,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 193, no. 52, pp. 5585–5595, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Liu, P. Zhuang, V. Anh, and I. Turner, “A fractional-order implicit difference approximation for the space-time fractional diffusion equation,”
*The ANZIAM Journal*, vol. 47, pp. C48–C68, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - C. F. Lorenzo and T. T. Hartley, “Initialized fractional calculus,”
*International Journal of Applied Mathematics*, vol. 3, no. 3, pp. 249–265, 2000. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. F. Lorenzo and T. T. Hartley, “Variable order and distributed order fractional operators,”
*Nonlinear Dynamics*, vol. 29, no. 1–4, pp. 57–98, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. F. M. Coimbra, “Mechanics with variable-order differential operators,”
*Annalen der Physik*, vol. 12, no. 11-12, pp. 692–703, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. F. M. Coimbra and L. E. S. Ramirez, “A variable order constitutive relation for viscoelasticity,”
*Annalen der Physik*, vol. 16, no. 7-8, pp. 543–552, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. F. M. Coimbra, C. M. Soon, and M. H. Kobayashi, “The variable viscoelasticity oscillator,”
*Annalen der Physik*, vol. 14, no. 6, pp. 378–389, 2005. View at Publisher · View at Google Scholar · View at Scopus