Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 696597, 6 pages

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

## A Possible Generalization of Acoustic Wave Equation Using the Concept of Perturbed Derivative Order

^{1}Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa^{2}Department of Mathematics and Institute for Mathematical Research, University Putra Malaysia, 43400 Serdang, Malaysia

Received 18 February 2013; Accepted 18 March 2013

Academic Editor: Guo-Cheng Wu

Copyright © 2013 Abdon Atangana and Adem Kılıçman. 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 standard version of acoustic wave equation is modified using the concept of the generalized Riemann-Liouville fractional order derivative. Some properties of the generalized Riemann-Liouville fractional derivative approximation are presented. Some theorems are generalized. The modified equation is approximately solved by using the variational iteration method and the Green function technique. The numerical simulation of solution of the modified equation gives a better prediction than the standard one.

#### 1. Introduction

Acoustics was in the beginning the study of small pressure waves in air which can be detected by the human ear: sound. The possibility of acoustics has been extended to higher and lower frequencies: ultrasound and infrasound. Structural vibrations are now often included in acoustics. Also the perception of sound is an area of acoustical research. In our present paper we will limit ourselves to the original definition and to the propagation in fluids like air and water. In such a case acoustics is a part of fluid dynamics. A major problem of fluid dynamics is that the equations of motion are nonlinear. This implies that an exact general solution of these equations is not available. Acoustics is a first-order approximation in which nonlinear effects are neglected. The corresponding relative density fluctuations are considered very small [1]. The acoustic wave equation governs the propagation of acoustic waves through a material medium. The form of the equation is a second-order partial differential equation. The equation describes the evolution of acoustic pressure or particle velocity as a function of position and time . A simplified form of the equation that describes acoustic waves in only one spatial dimension is considered in this paper

A derivation of general linearized wave equations is discussed by Pierce and Goldstein [1, 2]. However, neglecting the nonlinear effects in this equation, may lead to inaccurate prediction of the propagation of acoustic wave through the medium. Therefore in order to include explicitly the effect of corresponding relative density fluctuations in the mathematical formulation, one needs to insert it in the partial differential equation that governs the propagation of the acoustic wave. Recently, the acoustic equation was extended to the concept of fractional order derivative in [3]. Therefore in this paper, our concern is the modification of the previous equation by perturbing the order of the first derivative by replacing the first order of the derivative with where is a positive small parameter. Also, when we consider diffusion process in porous medium, if the medium structure or external field changes with time, in this situation, the ordinary integer order and constant-order fractional diffusion equation model cannot be used to well characterize such phenomenon (see [3–9]).

#### 2. Definitions and Approximation

To describe the propagation of acoustic waves through a material medium with coordinate and time-dependent perturbed dimension, one must use Riemann-Liouville fractional order derivative that was introduced and used in a number of works (see [3, 6–8, 10]). These derivatives are defined as (see [3, 6–8])

Here, is the Euler gamma function; , where is the integer part of for , that is, and for . Following (2) we have that and . The integral operator defined previously for fractional exponents and depending on coordinates and time can be expressed in terms of ordinary derivative and integral [11] for . Here and are considered as the corresponding relative density fluctuations that vary slightly in time and space, respectively. For this matter, generalized Riemann-Liouville fractional derivatives satisfy the approximate relations:

The previous relations make it possible to describe the dynamic system, including the effect of the corresponding relative density fluctuations, by means of partial differential and integral equations.

#### 3. Some Properties of the Approximation

Let us examine some properties of the previous derivative operator.

(i) *Addition *

If , , and are differentiable in the opened interval **I** then,

* Proof. *We have

(ii) *Division *

If and are differentiable on the open interval **I **then

(iii) *Multiplication *

If , and are differentiable in the open interval **I** then

(iv) *Power *

If and are differentiable in the open interval **I** then

If and are two times differentiable in the open interval **I** then

##### 3.1. Clairaut’s Theorem for the Approximation

Assume that , , and are functions for which , , , and exist and are continuous over a domain then, and exist and are continuous over the domain D. If in addition then

*Proof. *If , , and are functions for which , ,, and exist and are continuous over a domain then
Now interchanging by we obtain .

If = then according to Clairaut’s theorem; thus replacing by in , we obtain that

##### 3.2. Chain-Rule for the Approximation

We have

##### 3.3. Rolle's Theorem for the Approximation

If a real-valued functions and are continuous on a closed interval [*a*, *b*], differentiable on the open interval , and , then there exist a in the open interval and a small parameter such that

*Proof. *Following Rolle's theorem, there exists a *c* in the open interval * *such that . For this we have that

If , , and are differentiable in an open interval **I**, then there exist and such that

*Proof. * Let ; then
But is very small such that and ; it follows that

It is important to observe that if , we recover the properties of normal derivatives.

#### 4. Modification of the Equation

In order to include explicitly the possible effect of the corresponding relative density fluctuations into the mathematical formulation, in this paper, we replace the classical version of the derivative of (1) by the modified Riemann-Liouville fractional derivative approximation (3) to obtain Making use of (3) and relation (8), the previous equation can be transformed to the following partial differential equation for : Omitting the terms of in the previous equation, we obtain the following:

Now since the small parameters representing the perturbation additions to unity are small, the right- and left-hand sides of (21) can be divided by . In this case, we obtain Here, We have,

It is easy to observe that (22) differs from (1) in three properties.

First the velocity of the sound in this case depends on time and coordinates due to the effect of the corresponding relative density fluctuations. Secondly the force appears due to the coordinate and time dependence of the corresponding relative density fluctuations, and such force was considered in [12]. Third there is a derivative-free term that depends on both time and space, and it is proportional to and characterizes, depending on the coefficient sign, the retardation or enhancement propagation of acoustic waves through a material medium. Therefore, even weak memory, which is taken into account by generalized Riemann-Liouville fractional derivatives and presents the characteristics of a corresponding relative density fluctuations, transforms constant-coefficient velocity to varying-coefficient velocity. Moreover, this memory is responsible for a force with which the corresponding relative density fluctuations act on a propagation of acoustic waves through a material medium. This force appears only if the propagating acoustic wave has memory depending on coordinates and time; that it “remembers” its trajectories and time. Those terms in (22) that involve the corresponding relative density fluctuations additions ( and ) to the time and space dimensions are very small. Also the exact analytical solution of this modified equation is not easy to be determined. Therefore, this equation can be solved approximately by changing the function to , which satisfies (1).

#### 5. Solutions of the Modified Equation

The modified equation can be reformulated as follows:

To solve (25) we need to give explicitly and . For example if one consider these functions to be where and are very significantly small such that . Since these parameters are very small then the solution of (25) can be sought in the form where is the solution of (1) and is proportional to and . Here and are frequencies characterizing variation in corresponding relative density fluctuations.

To solve (25) together with (26) we make use of two techniques including the Green function techniques and the variational iterative decomposition technique. Here we will start with the variational iteration techniques.

##### 5.1. Variational Iteration Method

Variational iteration method has been favourably applied to various kinds of nonlinear problems. The main property of the method is in its flexibility and ability to solve nonlinear equations accurately and conveniently. Very recently it was recognized that the variational iteration method [11, 13–20] can be an effective procedure for solution of various nonlinear problems without usual restrictive assumptions. In this paper we will make use of this iterative decomposition technique to solve the modified wave (26) together with (25). To solve (25) by means of variational iteration method, we put (25) in the form of

The correction functional for (5) can be approximately expressed for this matter as follows: where is a general Lagrange multiplier [21, 22], which can be recognized optimally by means of variation assumption [21, 22]; here , and are considered as constrained variations. Making the previous functional stationary, we obtain

Capitulates the next Lagrange multipliers, giving up to the following Lagrange multipliers for the case where and for . For these matters , we obtain the following iteration formula:

It is worth noting that if the zeroth component is defined, then the remaining components can be completely determined such that each term is determined by using the previous terms, and the series solutions are thus entirely determined. Finally, the solution is approximated by the truncated series

Here we choose the first term to be zero meaning and the second term can be determined as

Our next concern is to define ; that is, we first need to provide the solution of (1) which is found in the literatures [23]. The following solutions are obtained by separation of variables in different coordinate systems. They are phasor solutions, meaning that they have an implicit time-dependence factor of where is the angular frequency. The explicit time dependence is given by (33)

Introducing the previous expression in , we obtain the second expression. In this matter two components of the decomposition series were obtained of which was evaluated to have the following expansion:

The next figures show the graphical representation of the approximated solution of the modified acoustic wave equation and the exact solution of the standard version of acoustic wave equation as function of space and time (contour plot and density plot of both solution). A contour plot gives essentially a topographic map of a solution and the density plot shows the values of the function at a regular erray of points and lighter region of the contour plot is higher.

From the next Figures 1, 2, 3, and 4, respectively, one can see that there are more details with the solution of the modified acoustic wave equation than in the standard solution, meaning that the details left out by neglecting the small effect of the correspondent relative density fluctuations are very important when one needs to observe the propagation of the acoustic wave through the material medium.The approximate solution of (25) has been depicted in Figures 1 and 3 and the exact solution in Figures 2 and 4.

##### 5.2. Green Function Method

To solve (25) together with (26) by means of Green function technique, one needs first to construct a suitable green function.

If G is the green function to be constructed, then G must satisfy the following equation:

We are lucky enough, because the Green function to be constructed here is the green function of the wave equation and is given later in the case of the closed forms for the Green function for the infinite one-dimensional domain [23]: where

Following the Green function technique, the general solution of the modified acoustic wave equation is given later as where remains the same as defined earlier in Section 5.1.

#### 6. Conclusion

In this paper, an acoustic wave equation was extended to the concept of the modified Riemann-Liouville fractional order derivative. We presented in detail some properties of the generalized Riemann-Liouville fractional order derivative approximation. We presented the analysis of the generalized equation. We highlighted the three differences between the generalized equation and the standard one. First the velocity of the sound in this case depends on time and coordinates due to the effect of the corresponding relative density fluctuations. Second the force appears due to the coordinate and time dependence of the corresponding relative density fluctuations. Third there is a derivative-free term that depends on both time and space which is proportional to and characterizes, depending on the coefficient sign, the retardation or enhancement propagation of acoustic waves through a material medium. The modified equation is approximately solved by using the variational iteration method and the Green function technique. The solution of the modified equation gives better prediction than the standard one.

#### References

- K. Aki and P. Richards,
*Quantitative Seismology*, Theory and Methods: Freeman, New York, NY, USA, 1980. - A. D. Pierce,
*Acoustics: An Introduction to Its Physical Principles and Applications*, McGraw-Hill Book Company, New York, NY, USA, 1981. - W. Chen and S. Holm, “Fractional Laplacian time-space models for linear and nonlinear lossy media exhibiting arbitrary frequency dependency,”
*Journal of the Acoustical Society of America*, vol. 115, no. 4, pp. 1424–1430, 2004. View at Publisher · View at Google Scholar - L. E. S. Ramirez and C. F. M. Coimbra, “On the selection and meaning of variable order operators for dynamic modeling,”
*International Journal of Differential Equations*, vol. 2010, Article ID 846107, 16 pages, 2010. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. E. S. Ramirez and C. F. M. Coimbra, “On the variable order dynamics of the nonlinear wake caused by a sedimenting particle,”
*Physica D*, vol. 240, no. 13, pp. 1111–1118, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Ross and S. Samko, “Fractional integration operator of variable order in the Hölder spaces ${H}^{\lambda (x)}$,”
*International Journal of Mathematics and Mathematical Sciences*, vol. 18, no. 4, pp. 777–788, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Umarov and S. Steinberg, “Variable order differential equations with piecewise constant order-function and diffusion with changing modes,”
*Zeitschrift für Analysis und ihre Anwendungen*, vol. 28, no. 4, pp. 431–450, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. L. Magin, O. Abdullah, D. Baleanu, and X. J. Zhou, “Anomalous diffusion expressed through fractional order differential operators in the Bloch-Torrey equation,”
*Journal of Magnetic Resonance*, vol. 190, no. 2, pp. 255–270, 2008. View at Publisher · View at Google Scholar - A. Atangana and J. F. Botha, “Generalized groundwater flow equation using the concept of variable order derivative,”
*Boundary Value Problems*, vol. 2013, p. 53, 2013. View at Publisher · View at Google Scholar - H. G. Sun, W. Chen, and Y. Q. Chen, “Variable order fractional differential operators in anomalous diffusion modeling,”
*Physica A*, vol. 388, no. 21, pp. 4586–4592, 2009. View at Publisher · View at Google Scholar - G. C. Wu, “New trends in the variational iteration method,”
*Communications in Fractional Calculus*, vol. 2, no. 2, pp. 59–75, 2011. View at Google Scholar - L. Ya. Kobelev, L. Ya. Kobelev, and L. Yu. Klimontovicht " ", “Anomalous diffusion with time and coordinate-dependent memory,”
*Dokiady Physics*, vol. 48, no. 6, pp. 264–268, 2003. View at Publisher · View at Google Scholar - G. C. Wu and D. Baleanu, “Variational iteration method for the Burgers' flow with fractional derivatives-New Lagrange multipliers,”
*Applied Mathematical Modelling*, vol. 37, no. 9, pp. 6183–6190, 2013. View at Publisher · View at Google Scholar - G. C. Wu and D. Baleanu, “Variational iteration method for fractional calculus—a universal approach by Laplace transform,”
*Advances in Difference Equations*, vol. 2013, p. 18, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - S. A. El-Wakil, E. M. Abulwafa, M. A. Zahran, and A. A. Mahmoud, “Time-fractional KdV equation: formulation and solution using variational methods,”
*Nonlinear Dynamics*, vol. 65, no. 1-2, pp. 55–63, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Atangana, A. Ahmed, and N. Bildik, “A generalized version of a low velocity impact between a rigid sphere and a transversely isotropic strain-hardening plate supported by a rigid substrate using the concept of non-integer derivatives,”
*Abstract Applied Analysis*, vol. 2013, Article ID 671321, 9 pages, 2013. View at Publisher · View at Google Scholar - A. Atangana and A. Secer, “Time-fractional coupled- the korteweg-de vries equations,”
*Abstract Applied Analysis*, vol. 2013, Article ID 947986, 8 pages, 2013. View at Publisher · View at Google Scholar - J. J. S. Duan, R. Rach, D. Bulean, and A. M. Wazwaz, “A review of the Adomian decomposition method and its applications to fractional differential equations,”
*Communications in Fractional Calculus*, vol. 3, no. 2, pp. 73–99, 2012. View at Google Scholar - N. Bildik and A. Konuralp, “The use of variational iteration method, differential transform method and Adomian decomposition method for solving different types of nonlinear partial differential equations,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 7, no. 1, pp. 65–70, 2006. View at Google Scholar - Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 7, no. 1, pp. 27–34, 2006. View at Google Scholar - M. Inokuti, H. Sekine, and T. Mura, “General use of the Lagrange multiplier in non-linear mathematical physics,” in
*Variational Method in the Mechanics of Solids*, S. Nemat-Nasser, Ed., pp. 156–162, Pergamon Press, Oxford, UK, 1978. View at Google Scholar - A. Atangana, “New class of boundary value problems,”
*Information Science Letters*, vol. 1, no. 2, pp. 67–76, 2012. View at Google Scholar - P. M. Morse and H. Feshbash,
*Methods of Theoretical Physics*, McGraw-Hill, New York, NY, USA, 1953.