Abstract and Applied Analysis

Volume 2013 (2013), Article ID 614874, 14 pages

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

Research Article

A One Step Optimal Homotopy Analysis Method for Propagation of Harmonic Waves in Nonlinear Generalized Magnetothermoelasticity with Two Relaxation Times under Influence of Rotation

S. M. Abo-Dahab,1,2 Mohamed S. Mohamed,1,3 and T. A. Nofal1,4

¹Mathematics Department, Faculty of Science, Taif University, P.O. Box 888, Saudi Arabia

²Mathematics Department, Faculty of Science, South Valley University, Qena 83523, Egypt

³Mathematics Department, Faculty of Science, Al-Azhar University, Nasr City 11884, Cairo, Egypt

⁴Mathematics Department, Faculty of Science, Minia University, Minia, Egypt

Received 1 May 2013; Revised 2 June 2013; Accepted 4 June 2013

Academic Editor: Santanu Saha Ray

Copyright © 2013 S. M. Abo-Dahab et al. 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 aim of this paper is to apply OHAM to solve numerically the problem of harmonic wave propagation in a nonlinear thermoelasticity under influence of rotation, thermal relaxation times, and magnetic field. The problem is solved in one-dimensional elastic half-space model subjected initially to a prescribed harmonic displacement and the temperature of the medium. The HAM contains a certain auxiliary parameter which provides us with a simple way to adjust and control the convergence region and rate of convergence of the series solution. This optimal approach has a general meaning and can be used to get fast convergent series solutions of the different type of nonlinear fractional differential equation. The displacement and temperature are calculated for the models with the variations of the magnetic field, relaxation times, and rotation. The results obtained are displayed graphically to show the influences of the new parameters.

1. Introduction

In the past recent years, much attention has been devoted to simulate some real-life problems which can be described by nonlinear coupled differential equations using reliable and more efficient methods. Nonlinear partial differential equations are useful in describing various phenomena in disciplines. The nonlinear coupled systems of partial differential equations often appear in the study of circled fuel reactor, high-temperature hydrodynamics, and thermoelasticity problems, see [1–4]. From the analytical point of view, a lot of work has been done for such systems. With the rapid development of nanotechnology, there appears an ever increasing interest of scientists and researchers in this field of science. Nanomaterials, because of their exceptional mechanical, physical, and chemical properties, have been the main topic of research in many scientific publications. Wave generation in nonlinear thermoelasticity problems has gained a considerable interest for its utilitarian aspects in understanding the nature of interaction between the elastic and thermal fields as well as the system of PDEs for its applications. A lot of applications were paid on existence, uniqueness, and stability of the solution of the problem, see [5–7].

Recently, much attention has been devoted to numerical methods, which do not require discretization of space-time variables or linearization of the nonlinear equations, among the homotopy analysis methods. Since most of the nonlinear FDEs cannot be solved exactly, approximate and numerical methods must be used. Some of the recent analytical methods for solving nonlinear problems include the homotopy analysis method HAM [8–14]. The HAM, first proposed in 1992 by Liao [8], has been successfully applied to solve many problems in physics and science. This method is applied to solve linear and nonlinear systems.

The homotopy perturbation method HPM has the merits of simplicity and easy execution. The homotopy perturbation method was first proposed by He [15]. Unlike the traditional numerical methods, the HPM does not need discretization and linearization. Most perturbation methods assume that a small parameter exists, but most nonlinear problems have no small parameter at all. Many new methods have been proposed to eliminate the small parameter. Recently, the applications of homotopy theory among scientists appeared, and the homotopy theory became a powerful mathematical tool, when it is successfully coupled with perturbation theory.

Recently, Gepreel et al. [16] investigated the homotopy perturbation method and variational iteration method for harmonic waves propagation in nonlinear magnetothermoelasticity with rotation. Abd-Alla and Abo-Dahab [17] investigated the effect of rotation and initial stress on an infinite generalized magnetothermoelastic diffusion body with a spherical cavity. Abo-Dahab and Mohamed [18] studied the influence of magnetic field and hydrostatic initial stress on reflection phenomena of P (Primary) and SV (Shear Vertical) waves from a generalized thermoelastic solid half space. Abd-Alla and Mahmoud [19] investigated the magnetothermoelastic problem in rotating non-homogeneous orthotropic hollow cylinder under the hyperbolic heat conduction model. Abd-Alla et al. [20] studied the thermal stresses effect in a non-homogeneous orthotropic elastic multilayered cylinder. Abd-Alla et al. [21] studied the generalized magnetothermoelastic Rayleigh waves in a granular medium under the influence of a gravity field and initial stress. Abd-Alla and Abo-Dahab [22] investigated the time-harmonic sources in a generalized magnetothermoviscoelastic continuum with and without energy dissipation.

In the present paper, investigation is devoted for solving numerically the problem of harmonic wave propagation in a nonlinear thermoelasticity under influence of magnetic field, thermal relaxation times, and rotation. The problem is solved in one-dimensional elastic half-space model subjected initially to a prescribed harmonic displacement and the temperature of the medium.

The HAM contains a certain auxiliary parameter which provides us with a simple way to adjust and control the convergence region and rate of convergence of the series solution. The -curve of the third-order approximate solutions is displayed graphically to show the interval that the exact and approximate solutions take the same values. The displacement and temperature are calculated for the methods with the variations of the magnetic field and rotation. The results obtained are displayed graphically to show the influences of the new parameters.

2. A One-Step Optimal Homotopy Analysis Method for PDEs

To describe the basic ideas of the HAM, we consider the following general nonlinear differential equation:

where is a nonlinear operator for this problem, and denote independent variables, and is an unknown function.

By means of the HAM, one first constructs the following zero-order deformation equation:

where is the embedding parameter, is an auxiliary parameter, is an auxiliary function, is an auxiliary linear operator, and   is an initial guess. Obviously, when and , it holds

Liao [8, 9] expanded in Taylor series with respect to the embedding parameter , as follows:

where

Assume that the auxiliary linear operator, the initial guess, the auxiliary parameter , and the auxiliary function are selected such that the series (4) is convergent at ; then we have from (4)

Let us define the vector

Differentiating (2) times with respect to , then setting and dividing then by , we have the following th-order deformation equation:

where

Applying the integral operator on both sides of (8), we have

where the th-order deformation equation (8) can be easily solved, especially by means of symbolic computation software such as Mathematica, Maple, and MathLab. The convergence of the homotopy analysis method for solving these equations is discussed in [23].

Abbasbandy and Jalili [24] and Turkyilmazoglu [25–29] applied the homotopy analysis method to nonlinear ODEs and suggested the so-called optimization method to find out the optimal convergence control parameters by minimum of the square residual error integrated in the whole region having physical meaning. Their approach is based on the square residual error.

Let denote the square residual error of the governing equation (1) and express it as

where

the optimal value of is given by a nonlinear algebraic equation:

3. Application of HAM on the Nonlinear Magnetothermoelastic with Rotation Equations

In this section, we use the homotopy analysis method to calculate the approximate solutions of the following nonlinear magnetothermoplastic model with rotation equations

where , , , , , , , , and are arbitrary constants with the initial conditions

where is an arbitrary constant and the boundary conditions

To demonstrate the effectiveness of the method, we consider the system of nonlinear initial-value problem (14) with the initial conditions (15) and the boundary conditions (16) by choosing the linear operators

with the property , , where , are the integral constants and the nonlinear operators are defined as

Choosing for , the zeroth-order deformation equations are

where

Then, the th-order deformation equations become

where

For simplicity, we suppose ; the system (21) has the following general solutions:

In this case, where and are constants, the general solution of (23) is taking the following form:

The problems above can be readily solved by symbolic computation packages such as Mathematica. Thereupon, successive solving of these problems yields

Now we make calculations for the results obtained by the HAM using the Mathematica software package with the following arbitrary constants:

To investigate the influence of on the convergence of the solution series given by the HAM, we first plot the so-called -curves of and . According to the -curves, it is easy to discover the valid region of . We used terms in evaluating the approximate solution and . Note that the solution series contains the auxiliary parameter which provides us with a simple way to adjust and control the convergence of the solution series. In general, by means of the so-called -curve that is, a curve of a versus . As pointed by Liao [8] and Turkyilmazoglu [25], the valid region of is a horizontal line segment. Therefore, it is straightforward to choose an appropriate range for which ensures the convergence of the solution series (Tables 1 and 2). We sketch the -curve of and in Figure 1, which shows that the solution series is convergent when .

Table 1: The optimal values of at third-order approximate solutions of (14) when , for .

Optimal value of
Minimum value

1
0
−0.95623

0
1
−0.97543

1
0
−0.95623

0
1
−0.97543

Table 2: The optimal values of at 5th-order approximate solutions of (14) when , for .

Optimal value of
Minimum value

1
0
−0.990321

0
1
−0.995339

1
0
−0.990321

0
1
−0.995339

(a)

(b)

Figure 1: The -curve of the third-order approximate solutions of (14) when ,  ; for LS model when ,  ,  ,  ,  ,  , and .

4. Discussion

In order to gain physical insight, the temperature and radial displacement have been discussed by assigning numerical values to the parameter encountered in the problem in which the numerical results are displayed with the graphical illustrations in 2D and 3D formats. The variations are shown in Figures 1–15, with the view of illustrating the theoretical results obtained in the preceding sections; a numerical result is calculated for the homotopy analysis method.

(a)

(b)

Figure 2: The -curve of the third-order approximate solutions of (14) when ,  ; for GL model when , ,  ,  ,  ,  , and .

(a)

(b)

Figure 3: Variations of the displacement and the temperature for various values of the -axis and time when , , , , , , and .

(a)

(b)

Figure 4: Variations of the displacement and the temperature for various values of the -axis and rotation when , , , , , , and .

(a)

(b)

Figure 5: Variations of the displacement and the temperature for various values of the -axis and magnetic field when , , , , , , and .

(a)

(b)

Figure 6: Variations of the displacement and temperature for various values of the -axis and magnetic field when , , , , and .

(a)

(b)

Figure 7: Variations of the displacement and the temperature for various values of the -axis and time when , , , , and .

(a)

(b)

Figure 8: Variations of the displacement and the temperature for various values of the -axis and time when , , , , and .

(a)

(b)

Figure 9: Variations of the displacement and the temperature for various values of the -axis and rotation when , , , , , and .

(a)

(b)

Figure 10: Variations of the displacement and the temperature for various values of the -axis and magnetic field when , , , , , and .

(a)

(b)

Figure 11: Variations of the displacement and the temperature for various values of the -axis and magnetic field when ,  , , , , and .

(a)

(b)

Figure 12: The displacement as a function of time and the temperature without and with rotation and magnetic field (LS) at ,  , , , , ,  , , , , , , (, —), and (, , , - - -).

(a)

(b)

Figure 13: The displacement as a function of time and the temperature without and with rotation and magnetic field (GL) at ,  , , , , ,  , , , , , , (, —), and (, , , - - -).

(a)

(b)

Figure 14: The displacement as a function of time and the temperature for two values of (LS) at , , ,  , , , , ,  , , , , , , , — , and - - - .

(a)

(b)

Figure 15: The displacement as a function of time and the temperature for two values of (GL) at , , ,  , , , , , , , , , , , , — , and - - - .

Figures 1 and 2 display the -curve of the third-order approximate solutions (14) when ,  ; it is concluded that the displacement and temperature increase with increasing the values of to their maxima and then decrease with the high values of ; also, it is shown that is convergent when and is convergent when .

Figures 2–7 show the variations of the radial displacement and temperature with respect to axial , respectively, for different values of the time , rotation , and sensitive parts of the magnetic fields and . In both figures, it is clear that the radial displacement and temperature have a zero value only in a bounded region of space. It is observed from Figure 3 that the displacement and the temperature start from their maximum values, decrease, and increase periodically with an increase of the coordinate ; also, it is obvious that their values take the minimum values and increases with the increasing values of the time . From Figure 4, one can see that and decrease with an increase of the rotation . It is shown that the components of the displacement and the temperature start from the minimum values near zero, increase, and then decrease periodically with the coordinate ; it is clear also that there is a slight increase with an increase of the sensitive parts of the magnetic field (see, Figures 5 and 6). It is shown that the increasing of the coordinate sensitive an increasing and decreasing on them periodically due to appearance of the pairs () in the initial condition and the approximate solutions; it is also clear that the components begin from their minimum values and increase absolutely with the variation of the time . The variations of the rotation and magnetic field tend to slightly affect the displacement and the temperature.

From Figures 7 and 8 (GL model), it is clear that the displacement component and temperature if the rotation and magnetic field are vanish, take larger values than the corresponding values with the rotation and magnetic field effects.

Figures 9, 10, and 11 show the variations of the displacement and temperature with respect to the time with LS and GL models; it is shown that the radial displacement and the temperature increase with an increase of that takes a slight change with the rotation if , if , and if .

Figures 12 and 13 show clearly the variations of the displacement and temperature in the presence and absence of the rotation and sensitive magnetic field; it is observed that and in presence of the parameters are smaller than the corresponding values in the absence of , , and , but there is a slight change for LS and GL models.

Finally, Figures 14 and 15 show the variations of the displacement and the temperature with respect to the time for different values of and with and without rotation and magnetic field effects, respectively. It is obvious that the radial displacement and the temperature increase with an increase of ; the displacement increases with an increase of parameter, but the temperature is not affected by . Is also seen that the radial displacement and the temperature take large values with the rotation and magnetic field effects. Also, it is concluded that takes large values for LS comparing with those in GL model, vice versa for the temperature.

It is obvious that if the rotation and the sensitive part of the magnetic field are neglected, the approximate solutions obtained by HAM agree with the results obtained by Sweilam and Khader [1], taking into consideration VIM. Finally, it is obvious that the displacement takes large values if there are no rotation, thermal relaxation times, and sensitive part of the magnetic field parameters compared with the corresponding value with the rotation and magnetic fields parameters.

The results indicate that the effect of the rotation and the magnetic field on the radial displacement and the temperature is very pronounced.

5. Conclusion

Due to the complicated nature of the governing equations of the magnetothermoelastic, the finished works in this field are unfortunately limited. The method used in this study provides a quite successful in dealing with such problems. This method gives numerical solutions in the elastic medium without any restrictions on the actual physical quantities that appear in the governing equations of the considered problem. Important phenomena are observed in these computations.

(i)The homotopy analysis method has been successfully applied to obtain the numerical solutions of the nonlinear equation with initial conditions. The reliability of this method and reduction in computations give this method a wider applicability. HAM contains a certain auxiliary parameter which provides us with a simple way to adjust and control the convergence region and rate of convergence of the series solution. It was also demonstrated that the Adomian decomposition method, homotopy perturbation method, and variational iteration method are specialcases of. HAM is clearly a very efficient and powerful technique for finding the numerical solutions of the proposed equation. It therefore provides more realistic series solutions that generally converge very rapidly in real physical problems. HAM provides us with a convenient way of controlling the convergence of approximation series, which is a fundamental qualitative difference in analysis between HAM and other methods. The illustrative examples suggest that HAM is a powerful method for nonlinear problems in science and engineering. Mathematica has been used for computations in this paper. (ii)It was found that for large values of time the large and the generalization give numerical results. The case is quite different when we consider small values of rotation and magnetic field. The coupled theory predicts infinite speeds of wave propagation. The solutions obtained in the context of generalized thermoelasticity theory, however, exhibit the behavior of finite speeds of wave propagation.(iii)By comparing Figures 1–15 for thermoelastic medium with presence and absence of the rotation and magnetic field, it was found that they have the same behavior in both media. The effect of rotation and sensitive parts of the magnetic field is strongly effective in the displacement and temperature of the propagation of the harmonic waves propagation in nonlinear thermoelasticity.(iv)The results presented in this paper will be very helpful for researchers concerned with material science, designers of new materials, and low-temperature physicists, as well as for those working on the development of a theory of hyperbolic propagation of hyperbolic thermoelastic. The study of the phenomenon of rotation and magnetic field is also used to improve the conditions of oil extractions.

References

N. H. Sweilam and M. M. Khader, “Variational iteration method for one dimensional nonlinear thermoelasticity,” Chaos, Solitons & Fractals, vol. 32, no. 1, pp. 145–149, 2007.

N. H. Sweilam, “Harmonic wave generation in non linear thermoelasticity by variational iteration method and Adomian's method,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 64–72, 2007.

H. N. Sweilam, M. M. Khader, and F. R. Al-Bar, “On the numerical simulation of population dynamics with density-dependent migrations and the Allee effects,” Journal of Physics, vol. 96, no. 1, Article ID 012008, 10 pages, 2008.

N. H. Sweilam, M. M. Khader, and R. F. Al-Bar, “Nonlinear focusing Manakov systems by variational iteration method and adomian decomposition method,” Journal of Physics, vol. 96, no. 1, Article ID 012164, 7 pages, 2008.

A. N. Abd-Alla, A. F. Ghaleb, and G. A. Maugin, “Harmonic wave generation in nonlinear thermoelasticity,” International Journal of Engineering Science, vol. 32, no. 7, pp. 1103–1116, 1994.

C. A. De Moura, “A linear uncoupling numerical scheme for the nonlinear coupled thermo-dynamics equations,” in Numerical Methods, V. Pereyra and A. Reinoze, Eds., vol. 1005 of Lecture Notes in Mathematics, pp. 204–211, Springer, Berlin, Germany, 1983.

M. Slemrod, “Global existence, uniqueness, and asymptotic stability of classical smooth solutions in one-dimensional nonlinear thermoelasticity,” Archive for Rational Mechanics and Analysis, vol. 76, no. 2, pp. 97–133, 1981.

S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problem [Ph.D. thesis], Shanghai Jiao Tong University, 1992.

S. J. Liao, “An approximate solution technique not depending on small parameters: a special example,” International Journal of Non-Linear Mechanics, vol. 30, no. 3, pp. 371–380, 1995.

K. M. Hemida and M. S. Mohamed, “Numerical simulation of the generalized Huxley equation by homotopy analysis method,” Journal of Applied Functional Analysis, vol. 5, no. 4, pp. 344–350, 2010.

K. Hemida and M. S. Mohamed, “Application of homotopy analysis method to fractional order generalized Huxley equation,” Journal of Applied Functional Analysis, vol. 7, no. 4, pp. 367–372, 2012.

M. S. Mohamed, “Analytic approximations for fractional-order Newton-Leipnik system,” Journal of Advanced Research in Scientific Computing, vol. 3, no. 2, pp. 56–69, 2011.

M. S. Mohamed and H. A. Ghany, “Analytic approximations for fractional-order hyperchaotic system,” Journal of Advanced Research in Dynamical and Control Systems, vol. 3, no. 3, pp. 1–12, 2011.

H. A. Ghany and M. S. Mohammed, “White noise functional solutions for Wick-type stochastic fractional KdV-Burgers-Kuramoto equations,” Chinese Journal of Physics, vol. 50, no. 4, pp. 619–627, 2012.

J.-H. He, “Asymptotology by homotopy perturbation method,” Applied Mathematics and Computation, vol. 156, no. 3, pp. 591–596, 2004.

K. A. Gepreel, S. M. Abo-Dahab, and T. A. Nofal, “Homotopy perturbation method and variational iteration method for harmonic waves propagation in nonlinear magneto-thermoelasticity with rotation,” Mathematical Problems in Engineering, vol. 2012, Article ID 827901, 30 pages, 2012.

A. M. Abd-Alla and S. M. Abo-Dahab, “Effect of rotation and initial stress on an infinite generalized gagneto-thermoelastic diffusion body with a spherical cavity,” Journal of Thermal Stresses, vol. 35, pp. 892–912, 2012.

S. M. Abo-Dahab and R. A. Mohamed, “Influence of magnetic field and hydrostatic initial stress on wave reflection from a generalized thermoelastic solid half-space,” JVC/Journal of Vibration and Control, vol. 16, no. 5, pp. 685–699, 2010.

A. M. Abd-Alla and S. R. Mahmoud, “Magneto-thermoelastic problem in rotating non-homogeneous orthotropic hollow cylinder under the hyperbolic heat conduction model,” Meccanica, vol. 45, no. 4, pp. 451–462, 2010.

A. M. Abd-Alla, A. N. Abd-Alla, and N. A. Zeidan, “Thermal stresses in a nonhomogeneous orthotropic elastic multilayered cylinder,” Journal of Thermal Stresses, vol. 23, no. 5, pp. 413–428, 2000.

A. M. Abd-Alla, S. M. Abo-Dahab, H. A. Hammad, and S. R. Mahmoud, “On generalized magneto-thermoelastic Rayleigh waves in a granular medium under the influence of a gravity field and initial stress,” Journal of Vibration and Control, vol. 17, no. 1, pp. 115–128, 2011.

A. M. Abd-Alla and S. M. Abo-Dahab, “Time-harmonic sources in a generalized magneto-thermo-viscoelastic continuum with and without energy dissipation,” Applied Mathematical Modelling, vol. 33, no. 5, pp. 2388–2402, 2009.

M. A. F. Araghia and A. Fallahzadeh, “On the convergence of the Homotopy analysis method for solving the schrodinger equation,” Journal of Applied Sciences Research, vol. 2, pp. 6076–6083, 2012.

S. Abbasbandy and M. Jalili, “Determination of optimal convergence control parameter value in homotopy analysis method,” Numerical Algorithms, 2013.

M. Turkyilmazoglu, “Purely analytic solutions of magnetohydrodynamic swirling boundary layer flow over a porous rotating disk,” Computers and Fluids, vol. 39, no. 5, pp. 793–799, 2010.

M. Turkyilmazoglu, “Numerical and analytical solutions for the flow and heat transfer near the equator of an MHD boundary layer over a porous rotating sphere,” International Journal of Thermal Sciences, vol. 50, no. 5, pp. 831–842, 2011.

M. Turkyilmazoglu, “Analytic approximate solutions of rotating disk boundary layer flow subject to a uniform suction or injection,” International Journal of Mechanical Sciences, vol. 52, no. 12, pp. 1735–1744, 2010.

M. Turkyilmazoglu, “Analytic approximate solutions of rotating disk boundary layer flow subject to a uniform vertical magnetic field,” Acta Mechanica, vol. 218, no. 3-4, pp. 237–245, 2011.

M. Turkyilmazoglu, “A note on the homotopy analysis method,” Applied Mathematics Letters, vol. 23, no. 10, pp. 1226–1230, 2010.