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

Authors: S. M. Abo-Dahab; Mohamed S. Mohamed; T. A. Nofal
Journal: Abstract and Applied Analysis
 (2013)

---

## 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.


---

## Body

## 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. Theh-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:(1)N[u(x,t)]=0,
where N  is a nonlinear operator for this problem, x and t denote independent variables, and u(x,t) is an unknown function.By means of the HAM, one first constructs the following zero-order deformation equation:(2)(1-q)ℒ(ϕ(x,t;q)-u0(x,t))=qhH(t)N[ϕ(x,t;q)],
where q∈[0,1] is the embedding parameter, h≠0 is an auxiliary parameter, H(t)≠0 is an auxiliary function, ℒ is an auxiliary linear operator, and  u0(x,t) is an initial guess. Obviously, when q=0 and q=1, it holds
(3)ϕ(x,t;0)=u0(x,t),ϕ(x,t;1)=u(x,t).Liao [8, 9] expanded ϕ(x,t;q) in Taylor series with respect to the embedding parameter q, as follows:
(4)ϕ(x,t;q)=u0(x,t)+∑m=1∞‍um(x,t)qm,
where
(5)um(x,t)=1m!∂mϕ(x,t;q)∂qm|q=0.Assume that the auxiliary linear operator, the initial guess, the auxiliary parameterh, and the auxiliary function H(t) are selected such that the series (4) is convergent at q=1; then we have from (4)
(6)u(x,t)=u0(x,t)+∑m=1∞‍um(x,t).Let us define the vector(7)u→n(x,t)={u0(x,t),u1(x,t),u2(x,t),…,un(x,t)}.Differentiating (2) m  times with respect to q, then setting q=0 and dividing then by m!, we have the following mth-order deformation equation:
(8)ℒ(um(x,t)-ϰmum-1(x,t))=hH(t)ℛm(u→m-1),
where
(9)ℛm(u→m-1)=1(m-1)!∂m-1𝒩[ϕ(x,t;q)]∂qm-1|q=0,ϰm={0,m≤1,1,m>1.Applying the integral operator on both sides of (8), we have
(10)um(x,t)=ϰmum-1(x,t)+h∫0tH(t)ℛm(u→m-1)dt,where themth-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Δ(h) denote the square residual error of the governing equation (1) and express it as
(11)Δ(h)=∫Ω(N[u~n(t)])2dΩ,
where
(12)u~m(t)=u0(t)+∑k=1m‍uk(t);
the optimal value of h is given by a nonlinear algebraic equation:
(13)dΔ(h)dh=0.


## 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(14)(1+σ1)utt+Ωut-uxx(1-σ2+2γux+3δux2)-β1(1-iωτ2)θx-β2(θux)x=0,(θ(1-iωτ1)-aux(1-iωδτ1)-12bux2)t-[(1+αux)θx]x=0,
where σ1, σ2, Ω, γ, β1, β2, a, b, and α are arbitrary constants with the initial conditions
(15)u(x,0)=θ(x,0)=A(1-cos(x)),ut(x,0)=θt(x,0)=0,
where A is an arbitrary constant and the boundary conditions
(16)u(0,t)=θ(0,t)=0,ut(0,t)=θt(0,t)=0.
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
(17)ℒ1[ϕ1(x,t;q)]=∂2ϕ1(x,t;q)∂t2,ℒ2[ϕ2(x,t;q)]=∂ϕ2(x,t;q)∂t,
with the property ℒ1[c1+c2t]=0, ℒ2[c3], where ci, (i=1,2,3) are the integral constants and the nonlinear operators are defined as
(18)N1[ϕ1,ϕ2]=(1+σ1)∂2ϕ1∂t2+Ω∂ϕ1∂t-∂2ϕ1∂x2×(1-σ2+2γ∂ϕ1∂x+3δ(∂ϕ1∂x)2)-β1∂ϕ2∂x+iωτ2β1∂ϕ2∂x-β2∂∂x(ϕ2∂ϕ1∂x),N2[ϕ1,ϕ2]=∂∂t((1-iωτ1)ϕ2-a(1-iωτ1δ)∂ϕ1∂x-12b(∂ϕ1∂x)2-∂∂x(1+α∂ϕ1∂x)∂ϕ2∂x(1-iωτ1)ϕ2-a(1-iωτ1δ)∂ϕ1∂x).
Choosing Hi(t)=1 for i=1,2, the zeroth-order deformation equations are
(19)(1-q)ℒ1[ϕ1(x,t;q)-u0(x,t)]=qh1N1[ϕ1(x,t;q),ϕ2(x,t;q)],(1-q)ℒ2[ϕ2(x,t;q)-v0(x,t)]=qh2N2[ϕ1(x,t;q),ϕ2(x,t;q)],
where
(20)ϕ1(x,t;0)=u0(x,t),ϕ1(x,t;1)=u(x,t),ϕ2(x,t;0)=v0(x,t),ϕ2(x,t;1)=θ(x,t).
Then, the mth-order deformation equations become
(21)ℒ1[um(x,t)-ϰmum-1(x,t)]=h1ℛ1m(u→m-1,θ→m-1),ℒ2[θm(x,t)-ϰmθm-1(x,t)]=h2ℛ2m(u→m-1,θ→m-1),
where
(22)ℛ1m(u→m-1,θ→m-1)=∂2um-1∂t2+11+σ1(∑j=0m-1‍Ω(um-1)t-(um-1)xx+σ2(um-1)xx-2γ∑j=0m-1‍(uj)xx(um-1-j)x-3δ∑i=0m-1‍∑j=0i‍(uj)x(ui-j)x(um-1-i)xx-β2∑j=0m-1‍(θj)x(um-1-j)x-β1(θm-1)x-β2∑j=0m-1‍(θj)(um-1-j)xx+iωτ2β1θm-1),ℛ2m(u→m-1,θ→m-1)=∂θm-1∂t+11-iωτ1(∑j=0m-1‍-a(um-1)xthhhhhhhhhhhhhhhhhhh-12b∑j=0m-1‍(uj)xt(um-1-j)xthhhhhhhhhhhhhhhhhhh-(θm-1)xxhhhhhhhhhhhhhhhhhhh-α∑j=0m-1‍(θm-1-j)x(uj)xxhhhhhhhhhhhhhhhhhhh-α∑j=0m-1‍(θm-1-j)xx(uj)xhhhhhhhhhhhhhhhhhhh-(θm-1)xx+aiωτ1δ(um-1)xt∑j=0m-1‍).
For simplicity, we suppose h1=h2; the system (21) has the following general solutions:
(23)um(x,t)=ϰmum-1(x,t)+h∬0tℛ1m(u→m-1,θ→m-1)dtdt,θm(x,t)=ϰmθm-1(x,t)+h∫0t‍ℛ2m(u→m-1,θ→m-1)dt.
In this case, where u0 and θ0 are constants, the general solution of (23) is taking the following form:
(24)u(x,t)=u0(x,t)+∑m=1∞‍um(x,t),θ(x,t)=θ0(x,t)+∑m=1∞‍θm(x,t).
The problems above can be readily solved by symbolic computation packages such as Mathematica. Thereupon, successive solving of these problems yields
(25)u0(x,0)=A(1-cos(x)),θ0(x,0)=A(1-cos(x)),u1(x,t)=Aht24(1+σ1)((-2-3A2δ2Aβ2(-cos(x)+cos(2x))+cos(x)(-2-3A2δ+3A2δcos(2x)-4Aγsin(x)+2σ2A2)+2isin(x)β1(i+ωτ2)),θ1(x,t)=Ahtcos(x)(1+2Aαsin(x))-1+iωτ1,u2(x,t)=Aht24(1+σ1)×((-2-3A2δ+3A2δcos(2x)2Aβ2(-cos(x)+cos(2x))+cos(x)×(-2-3A2δ+3A2δcos(2x)h-4Aγsin(x)+2σ2-2-3A2δ+3A2δcos(2x)(x))+2isin(x)β1(i+ωτ2)(-2-3A2δ+3A2δcos(2x))+1384(1+σ1)2×(Ah2t4(-135A4δ2cos(5x)hhhhhhhhhhhh+240A3δγsin(4x)-16A2hhhhhhhhhhhh×(2cos(x)-5cos2(x)+3cos(3x))hhhhhhhhhhhh×β22-96Aγsin(2x)hhhhhhhhhhhh×(1+2A2δ-σ2)+3A2cos(3x)hhhhhhhhhhhh×(48δ+63A2δ2+32γ2-48δσ2)hhhhhhhhhhhh+8Aβ2((2+3A2δ)-2(2+3A2δ)cos(x)hhhhhhhhhhhhhhhhhh+10cos(2x)hhhhhhhhhhhhhhhhhh+A(-9sin(3x))3Aδ((4x)3cos(2x)hhhhhhhhhhhhhhhhhhhhhhhhhh+6cos(3x)hhhhhhhhhhhhhhhhhhhhhhhihhh-7cos(4x))hhhhhhhhhhhhhhhhhhhhhh-2γ(sin(x)hhhhhhhhhhhhhhhhhhhhhhhhhhh+6sin(2x)hhhhhhhhhhhhhhhhhhhhhhhhhhh-9sin(3x)))hhhhhhhhhhhhhhhhhh+2(2cos(x)hhhhhhhhhhhhhhhhhhhhhh-5cos(2x))σ2(2+3A2δ))hhhhhhhhhhhh-2cos(x)(8+A2hhhhhhhhhhhhhhhhhhhhh×(3δ(8+9A2δ)hhhhhhhhhhhhhhhhhhhhhhhh+16γ2)hhhhhhhhhhhhhhhhhhhhh+8σ2(-2-3A2δ+σ2))hhhhhhhhhhhh+4iβ1(β2+σ2)4sin(x)hhhhhhhhhhhhhhhhhh+A(-8γcos(2x)hhhhhhhhhhhhhhhhhhhhhi+3Aδ(sin(x)hhhhhhhhhhhhhhhhhhhhhhhhhhhh-3sin(3x))-8γcos(2x)-8γcos(2x)(2-8γcos(2x)-8γcos(2x)x))hhhhhhhhhhhhhhhhhh-4sin(x)(A(-1+2cos(x))hhhhhhhhhhhhhhhhhhhhhhhhhhh×β2+σ2))hhhhhhhhhhhh×(i+ωτ0)-135A4δ2cos(5x))Ah2t4)-124(1+σ1)2(i+ωτ1)×(Ah2t3(-2-3A2δ+3A2δcos(2x)hhhhhhhhhhhh-4Aγsin(x)+2σ2-2-3A2δ+3A2δcos(2x))hhhhhh×(i+ωτ1)+2Aβ2hhhhhh×(-2icos(2x)(1+Ω+σ1-i.Ωωτ1)iAαsin(x)(1+σ1)-3Aiαsin(3x)hhhhhhhhh×(1+σ1)+2Ωcos(x)(i+ωτ1)hhhhhhhhh-2icos(2x)(1+Ω+σ1-iΩωτ1))hhhhhh+4β1(i+ωτ2)hhhhhh×(-2Aαcos(2x)(1+σ1)hhhhhhhhh+sin(x)(1+Ω+σ1-iωΩτ1))Ah2t3(-2-3A2δ+3A2δcos(2x)),θ2(x,t)=Ah(1+h)tcos(x)(1+2Aαsin(x))-1+iωτ1+A2h2t3α6(1+σ1)(i+ωτ1)×(-12i(2sin(2x)(1+3A2δ-σ2)2Aγcos(x)-6Aγcos(3x)gggggHHggg-9A2δsin(4x)hhhhhhhhhh+2A(sin(x)+sin(2x)hhhhhhhhhhhhhhh-3sin(3x))β2hhhhhhhhhh+2sin(2x)(1+3A2δ-σ2))hhhh+cos(2x)β1(i+ωτ2)hhhh+18(i+ωτ1)hhhh×(12iAh2t2hhhhhhh×(a1+σ1hhhhhhhhhh×((i+ωτ2)-4sin(x)hhhhhhhhihhhh+A(8νcos(2x)hhhhhhhhhhhihhhhh-3Aδ(sin(x)hhhhhhhhhhhhhhhhhhhhhhh-3sin(3x)))hhhhhhhhihhhh+4sin(x)(×β2+σ2A(-1+4cos(x))hhhhhhhhhhhhhhhhhihhhh×β2+σ2)hhhhhhhhhihiiii+4cos(x)β1(1-iωτ2))hhhhhhhhhh+iaωδ11+σ1((i+ωτ2)4sin(3x)hhhhhhhhhhhhhhhhhh+A(-8νcos(x)hhhhhhhhhhhhhihihhhhhhh+3Aδ(sin(x)hhhhhhhhhhhhhhhhhhhhhhhhhhhh-3sin(3x)))hhhhhhhhhhhhhhhhhh-4sin(x)(A(-1hhhhhhhhhhhhhhhhhhhhhhhhhihhhh+4cos(x))hhhhhhhhhhhhhhhhhhhhhhhhhhhh×β2+σ2)hhhhhhhhhhhihihhhhh+4icos(x)β1(i+ωτ2))hhhhhhhhhii×τ1-1i+ωτ1hhhhhhhhhii×4i((1+A2hhhhhhihhhhhhhhh+α2cos(x)+Aαhhhhhhihhhhhhhhh×(-3A×αcos(3x)hhhhhhihhhhhhhhhhh+5sin(2x))(1+A2)(1+A2)1i+ωτ1).Now we make calculations for the results obtained by the HAM using the Mathematica software package with the following arbitrary constants:(26)a=0.5,A=0.001,b=0.5,α=1,β1=0.5,β2=0.05,γ=1,δ=0.8,Ω=1,σ1=0.2,σ2=0.1.
To investigate the influence of h on the convergence of the solution series given by the HAM, we first plot the so-called h-curves of u(1,1) and θ(1,1). According to the h-curves, it is easy to discover the valid region of h. We used 3 terms in evaluating the approximate solution u(x,t)≅∑i=02ui(x,t) and θ(x,t)≅∑i=02θi(x,t). Note that the solution series contains the auxiliary parameter h 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 h-curve that is, a curve of a versus h. As pointed by Liao [8] and Turkyilmazoglu [25], the valid region of h is a horizontal line segment. Therefore, it is straightforward to choose an appropriate range for h which ensures the convergence of the solution series (Tables 1 and 2). We sketch the h-curve of u(1,1) and θ(1,1) in Figure 1, which shows that the solution series is convergent when -1.45<h<-0.5.Table 1
The optimal values ofh at third-order approximate solutions of (14) when x=1.0,t=1.0, for δ=1,σ1=0.2,σ2=0.1,Ω=0.1,andω=0.02.
τ
1
τ
2
Optimal value ofh
Minimum value
u
(
1,1
)
1
0
−0.95623
1.37312
×
10
-
3
0
1
−0.97543
7.9675
×
10
-
7
θ
(
1,1
)
1
0
−0.95623
1.3242
×
10
-
2
0
1
−0.97543
2.1432
×
10
-
5Table 2
The optimal values ofh at 5th-order approximate solutions of (14) when x=1.0,t=1.0, for δ=1,σ1=0.2,σ2=0.1,Ω=0.1,andω=0.02.
τ
1
τ
2
Optimal value ofh
Minimum value
u
(
1,1
)
1
0
−0.990321
1.37312
×
10
-
7
0
1
−0.995339
7.9675
×
10
-
9
θ
(
1,1
)
1
0
−0.990321
1.3242
×
10
-
4
0
1
−0.995339
2.1432
×
10
-
6Theh-curve of the third-order approximate solutions of (14) when x=1.0,  t=1.0; for LS model when τ2=0,  τ1=1,  δ=1,  σ1=0.2,  σ2=0.1,  Ω=0.1, and ω=0.02.
(a)
(b)


## 4. Discussion

In order to gain physical insight, the temperatureT and radial displacement u 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.Theh-curve of the third-order approximate solutions of (14) when x=1.0,  t=1.0; for GL model when τ2=0.2, τ1=0.1,  δ=0,  σ1=0.2,  σ2=0.1,  Ω=0.1, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and time t when τ2=0, τ1=0.1, δ=1, σ1=0.2, σ2=0.1, Ω=0.1, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and rotation Ω when t=0.1, τ2=0, τ1=0.1, δ=1, σ1=0.2, σ2=0.1, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and magnetic field σ1when t=0.1, τ2=0, τ1=0.1, δ=1, σ2=0.1, Ω=0.1, and ω=0.02.
(a)
(b)Variations of the displacementu and temperature θ for various values of the x-axis and magnetic field σ2 when t=0.1, τ2=0, τ1=0.1, δ=1, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and time t when Ω=σ1=σ2=0, τ2=0.2, τ1=0.1, δ=1, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and time t when Ω=σ1=σ2=0, τ2=0.2, τ1=0.1, δ=0, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and rotation Ω when σ1=σ2=0, t=0.1, τ2=0.2, τ1=0.1, δ=0, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and magnetic field σ1 when Ω=σ2=0, t=0.1, τ2=0.2, τ1=0.1, δ=0, and ω=0.02.
(a)
(b)Variations of the displacementu and the temperature θ for various values of the x-axis and magnetic field σ2 when Ω=σ1=0,  t=0.1, τ2=0.2, τ1=0.1, δ=0, and ω=0.02.
(a)
(b)The displacement as a function of time and the temperature without and with rotation and magnetic field (LS) atx=100,  a=0.5,  β1=0.5, A=0.001, b=0.5, α=1,  β2=0.05, γ=1, δ=0.8, τ2=0, τ1=0.1, δ=1, ϖ=0.02 (Ω=σ1=σ2=0, —), and (Ω=0.1, σ1=0.2, σ2=0.1, - -  -).
(a)
(b)The displacement as a function of time and the temperature without and with rotation and magnetic field (GL) atx=100,  a=0.5,  β1=0.5, A=0.001, b=0.5, α=1,  β2=0.05, γ=1, δ=0.8, τ2=0.2, τ1=0.1, δ=0, ϖ=0.02 (Ω=σ1=σ2=0, —), and (Ω=0.1, σ1=0.2, σ2=0.1, - -  -).
(a)
(b)The displacement as a function of time and the temperature for two values ofβ1 (LS) at  x=100, a=0.5,  A=0.001,  b=0.5, α=1, β2=0.05, γ=1, δ=0.8,  Ω=0.1, σ1=0.2, σ2=0.1, τ2=0, τ1=0.1, δ=1, ϖ=0.02, — β1=0.5, and - -  -  β1=0.005.
(a)
(b)The displacement as a function of time and the temperature for two values ofβ1 (GL) at x=100, a=0.5,  A=0.001,  b=0.5, α=1, β2=0.05, γ=1, δ=0.8,  Ω=0.1, σ1=0.2, σ2=0.1, τ2=0.2, τ1=0.1, δ=0, ϖ=0.02, — β1=0.5, and - - - β1=0.005.
(a)
(b)Figures1 and 2 display the h-curve of the third-order approximate solutions (14) when x=1.0,  t=1.0; it is concluded that the displacement u(1,1) and temperature θ(1,1) increase with increasing the values of h to their maxima and then decrease with the high values of h; also, it is shown that u(1,1) is convergent when -0.9<h<-0.6 and θ(1,1) is convergent when -0.6<h<-0.4.Figures2–7 show the variations of the radial displacement and temperature with respect to axial x, respectively, for different values of the time t, rotation Ω, and sensitive parts of the magnetic fields σ1 and σ2. 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 u and the temperature θ start from their maximum values, decrease, and increase periodically with an increase of the coordinate x; also, it is obvious that their values take the minimum values and increases with the increasing values of the time t. From Figure 4, one can see that u and θ decrease with an increase of the rotation Ω. It is shown that the components of the displacement u and the temperature θ start from the minimum values near zero, increase, and then decrease periodically with the coordinate x; 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 x sensitive an increasing and decreasing on them periodically due to appearance of the pairs (cos,sin) 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 t. The variations of the rotation and magnetic field tend to slightly affect  the displacement and the temperature.From Figures7 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.Figures9, 10, and 11 show the variations of the displacement and temperature with respect to the time t with LS and GL models; it is shown that the radial displacement and the temperature increase with an increase of t that takes a slight change with the rotation Ω if σ1=σ2=0, σ1 if Ω=σ2=0, and σ2 if Ω=σ1=0.Figures12 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 u and θ in presence of the parameters are smaller than the corresponding values in the absence of σ1, σ2, and Ω, but there is a slight change for LS and GL models.Finally, Figures14 and 15 show the variations of the displacement and the temperature with respect to the timet for different values of β=0.5,0.005 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 t; 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 parameterh  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 Figures1–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.



---

*Source: 614874-2013-09-29.xml*