# The Effect of Magnetic Field and Initial Stress on Fractional Order Generalized Thermoelastic Half-Space

Authors: Sunita Deswal; Sandeep Singh Sheoran; Kapil Kumar Kalkal
Journal: Journal of Mathematics
 (2013)

---

## Abstract

The aim of this paper is to study magneto-thermoelastic interactions in an initially stressed isotropic homogeneous half-space in the context of fractional order theory of generalized thermoelasticity. State space formulation with the Laplace transform technique is used to obtain the general solution, and the resulting formulation is applied to the ramp type increase in thermal load and zero stress. Solutions of the problem in the physical domain are obtained by using a numerical method of the Laplace inverse transform based on the Fourier expansion technique, and the expressions for the displacement, temperature, and stress inside the half-space are obtained. Numerical computations are carried out for a particular material for illustrating the results. Results obtained for the field variables are displayed graphically. Some comparisons have been shown in figures to present the effect of fractional parameter, ramp parameter, magnetic field, and initial stress on the field variables. Some particular cases of special interest have been deduced from the present investigation.


---

## Body

## 1. Introduction

Biot [1] developed the coupled theory of thermoelasticity to overcome the paradox inherent in the uncoupled theory that elastic changes have no effect on temperature. In this theory, the equations of elasticity and heat conduction are coupled. But, it shares the defect of the uncoupled theory in which it predicts infinite speed of propagation for heat waves. Generalized thermoelastic theories have been developed with the objective of removing the defect of coupled theory. At present, mainly two different models of generalized thermoelasticity are being extensively used, one proposed by Lord and Shulman [2] and the other proposed by Green and Lindsay [3]. The Lord-Shulman theory only modified the Fourier heat conduction equation and suggested one relaxation time, whereas the Green-Lindsay theory modified both the energy equation and the equation of motion and suggested two relaxation times. Dhaliwal and Sherief [4] extended the Lord-Shulman theory to include anisotropic case. Hetnarski and Ignaczak [5] presented a survey article of various representative theories in the range of generalized thermoelasticity. Youssef [6] studied two-dimensional generalized thermoelasticity problem for half-space subjected to ramp type heating.The theory of magneto-thermoelasticity has received attention of many researchers due to its application in widely diverse fields such as geophysics for understanding the effect of earth's magnetic field on seismic waves, damping of acoustic waves, emission of electromagnetic radiations from nuclear devices, optics, and so forth. The theory of magneto-thermoelasticity was introduced by Knopoff [7] and Chadwick [8] and developed by Kaliski and Petykiewicz [9]. The theoretical outline of the development of magneto-thermoelasticity was discussed by Paria [10]. Paria studied the propagation of plane magneto-thermoelastic waves in an isotropic unbounded medium under the influence of a magnetic field acting transversely to the direction of propagation. Nayfeh and Nemat-Nasser [11] studied the propagation of plane waves in a solid under the influence of electro-magnetic field. Sherief and Ezzat [12] discussed a thermal shock problem in magneto-thermoelasticity with thermal relaxation. Sherief and Helmy [13] illustrated a two-dimensional half-space problem in the context of electromagneto-thermoelasticity theory subjected to a nonuniform thermal shock. Ezzat and Youssef [14] constructed a problem on generalized magneto-thermoelasticity in a perfectly conducting medium. Baksi et al. [15] examined magneto-thermoelastic problem with thermal relaxation and heat source in three-dimensional, infinite rotating elastic medium.The development of initial stress in the medium is due to many reasons such as the process of quenching, resulting from difference of temperatures, slow process of creep, differential external forces, and gravity variations. The earth is supposed to be under high initial stress. The researchers have shown much interest to study the effect of these stresses on the propagation of waves. Biot [16] solved the dynamic problem of elastic medium under initial stress. Chattopadhyay et al. [17] studied the reflection of elastic waves under initial stress at a free surface. Montanaro [18] studied the isotropic linear thermoelasticity with hydrostatic initial stress by using Biot's linearization of the constitutive law for stress. Othman and Song [19] investigated the reflection of plane waves from an elastic solid half-space under hydrostatic initial stress without energy dissipation. Singh [20] explored the effect of hydrostatic initial stress on waves in a thermoelastic half-space.The theory of fractional derivative and integral was established in the second half of nineteenth century. The first application of fractional derivative was given by Abel who applied fractional calculus in the solution of an integral equation that arises in the formulation of tautochrone problem. In the recent years, fractional calculus has been applied successfully in various areas to modify many existing models of physical processes such as heat conduction, diffusion, viscoelasticity, wave propagation, and electronics. Caputo and Mainardi [21, 22] and Caputo [23] have established the relation between fractional derivative and theory of linear viscoelasticity. The generalization of the concept of derivative and integral to a noninteger order has been subjected to several approaches, and some various alternative definitions of fractional derivatives appeared in [24–27]. One can refer to Podlubny [28] for a survey of applications of fractional calculus. Povstenko [29] has proposed a quasistatic uncoupled theory of thermoelasticity based on fractional heat conduction equation. Sherief et al. [30] introduced a new model of thermoelasticity using fractional calculus, proved a uniqueness theorem, and derived a reciprocity relation and a variational principle. In this model, heat conduction equation takes the form as
(1)qi+τ0∂α∂tαqi=-kijθ,j,
where qi are the components of the heat flux vector, θ is the temperature, τ0 is the thermal relaxation time parameter, kij is thermal conductivity tensor, and α is a fractional parameter such that 0<α≤1. The above heat conduction equation reduces to the Maxwell-Cattaneo law in the limiting case when α→1. It should be mentioned here that the Maxwell-Cattaneo law has been employed by Lord and Shulman [2] to develop first generalized theory of thermoelasticity. Youssef [31] constructed another model of thermoelasticity in the context of a new consideration of heat conduction with a fractional order and proved the uniqueness theorem. In this model, Youssef described different cases of conductivity: 0<α<1 corresponds to weak conductivity, α=1 corresponds to normal conductivity, and 1<α≤2 corresponds to superconductivity. Ezzat [32, 33] established a model of fractional heat conduction equation by using the new Taylor series expansion of time-fractional order developed by Jumarie [34]. El-Karamany and Ezzat [35] introduced two general models of fractional heat conduction law for a nonhomogeneous anisotropic elastic solid. Uniqueness and reciprocal theorems are proved, and the convolutional variational principle is established and used to prove a uniqueness theorem with no restriction on the elasticity or thermal conductivity tensors except symmetry conditions. The two-temperature dynamic coupled, Lord-Shulman and fractional coupled thermoelasticity theories, result as limit cases. For fractional thermoelasticity not involving two-temperatures, El-Karamany and Ezzat [36] established the uniqueness, reciprocal theorems and convolution variational principle. The dynamic coupled and the Green-Naghdi thermoelasticity theories result as limit cases. The reciprocity relation in case of quiescent initial state is found to be independent of the order of differintegration [35, 36]. Fractional order theory of a perfect conducting thermoelastic medium not involving two temperatures was investigated by El-Karamany and Ezzat [37]. Kothari and Mukhopadhyay [38] studied a half-space problem under fractional order theory of thermoelasticity and analyzed the effect of the fractional order parameter on the field variables.In the present paper, we study the effect of magnetic field and initial stress under fractional order theory of thermoelasticity proposed by Sherief et al. [30]. We employ a state space approach developed by Bahar and Hetnarski [39] on the formulation. The Laplace transform technique is used to obtain the general solution. The inverse Laplace transform is carried out using a numerical inversion method developed by Honig and Hirdes [40]. Finally, the effect of fractional parameter, ramp parameter, magnetic field, and initial stress on field variables is displayed graphically.


## 2. Governing Equations

The governing equations in the context of fractional order theory of generalized thermoelasticity with initial stress and magnetic field for isotropic and homogeneous elastic medium are considered as(i)
the equation of motion(2)ρu¨i=σji,j+Fi,whereF→=μ0J→×H→,(ii)
heat conduction equation(3)kθ,ii=ρcE(1+τ0∂α∂tα)∂θ∂t+βT0(1+τ0∂α∂tα)∂ui,i∂t,(iii)
constitutive relations(4)σij=-p(δij+ωij)+2μ-eij+λ-eδij-βθδij,eij=12(ui,j+uj,i),ωij=12(uj,i-ui,j).We take the linearized Maxwell equations governing the electromagnetic field for a perfectly conducting medium as(5)curlh→=J→+ε0∂E→∂t,curlE→=-μ0∂h→∂t,E→=-μ0(∂h→∂t×H→),divh→=0,
where ui are the components of displacement vector u→, θ=T-T0, T is the absolute temperature, T0 is the reference temperature assumed to obey the inequality |θ/T0|≪1, τ0 is the thermal relaxation time, σij are the components of the stress tensor, eij are the components of strain tensor, δij is the Kronecker delta function, e is the cubical dilation, ρ is the density of the medium, cE is the specific heat, k is the thermal conductivity, β=(3λ-+2μ-)αt, αt is the coefficient of linear thermal expansion, λ- and μ- are the lame constants, p is the initial stress, Fi are the components of Lorentz's body force vector F→, μ0 is the magnetic permeability, ε0 is the electric permittivity, H→ is the applied magnetic field, h→ is the induced magnetic field, E→ is the induced electric field, and J→ is the current density vector.


## 3. Problem Formulation

We consider a perfectly conducting isotropic homogeneous and fractional order generalized thermoelastic half-space with hydrostatic initial stress subjected to a constant magnetic fieldH→(0,H0,0) which produces an induced magnetic field h→(0,h2,0) and induced electric field E→(0,0,E3). We assume one-dimensional motion for which all the field quantities are functions of x and t.The displacement components take the form(6)ux=u(x,t),uy=uz=0.The strain component becomes(7)e=exx=∂u∂x.The components of magnetic field vectors are(8)Hx=0,Hy=H0,Hz=0.The electric intensity vectorE→ is parallel to current density vector J→. Hence, components of E→ and J→ are given as
(9)Ex=Ey=0,Ez=E3,Jx=Jy=0,Jz=J3.Now, the Maxwell equation (5) provides the following results:
(10)E3=-μ0H0∂u∂t,h2=-H0e,J3=ε0μ0H0∂2u∂2t.Using (8) and (10) into the relation F→=μ0J→×H→, we obtain
(11)Fx=-ε0μ02H02u¨,Fy=Fz=0.The governing equations for one-dimensional case become(12)ρu¨=∂σ∂x+Fx,(13)k∂2θ∂x2=ρcE(1+τ0∂α∂tα)∂θ∂t+βT0(1+τ0∂α∂tα)∂e∂t,(14)σ=σxx=-p+(λ-+2μ-)e-βθ.Now, we will use the following nondimensional variables:(15)x′=c0η0x,u′=c0η0u,t′=c02η0t,τ0′=c02η0τ0,θ′=θT0,σ′=σλ-+2μ-,p′=pβT0,
where
(16)η0=ρcEk,c02=λ-+2μ-ρ.Expressing (12)–(14) in terms of the nondimensional variables given by (15) and dropping the prime sign for convenience, we have the following forms:
(17)∂2σ∂x2=α1e¨,∂2θ∂x2=(1+τ0*∂α∂tα)(θ˙+β1e˙),σ=e-ε1θ-δ1,
where
(18)α1=1+ε0μ02H02ρ,β1=βρcE,ε1=βT0λ-+2μ-,δ1=ε1p,τ0*=τ0(c02η0)1-α.These equations will be supplemented with appropriate boundary conditions relevant to the particular application under consideration.Taking the Laplace transform of (17) by using homogeneous initial conditions defined and denoted as
(19)f-(s)=∫0∞‍e-stf(t)dt,s>0,
we obtain
(20)∂2σ-∂x2=α1s2e-,(21)∂2θ-∂x2=(s+τ0*sα+1)(θ-+β1e-),(22)σ-=e--ε1θ--δ1s.Eliminating the value ofe-from (20) and (21) by using (22), we obtain
(23)∂2θ-∂x2=L1θ-+L2σ-+L3,∂2σ-∂x2=M1θ-+M2σ-+M3,
where
(24)L1=(1+β1ε1)(s+τ0*sα+1),L2=β1(s+τ0*sα+1),L3=β1δ1(1+τ0*sα),M1=α1ε1s2,M2=α1s2,M3=α1δ1s.


## 4. State-Space Formulation

Now, choosing the temperature of heat conductionθ and the stress component σ in x-direction as state variables, one can write (23) in matrix form as
(25)d2dx2V-(x,s)=A(s)V-(x,s)+C(x,s),
where
(26)A(s)=[L1L2M1M2],V-(x,s)=[θ-(x,s)σ-(x,s)],C(x,s)=[L3M3].The formal solution of (25) can be written in the form
(27)V-(x,s)=exp[-A(s)x][V-(0,s)+C(0,s)A-1(s)]-C(x,s)A-1(s),
where
(28)V-(0,s)=[θ-(0,s)σ-(0,s)]=[θ-0σ-0].In the above solution, we have cancelled the part of exponential having positive power to get bounded solution for largex.Now, we will use the Cayley-Hamilton theorem to find the form of the matrixexp[-A(s)x], and, for this, we proceed as follows.The characteristic equation of the matrixA(s) is obtained as
(29)λ2-(L1+M2)λ+(L1M2-M1L2)=0,
where the roots of (29), namely, λ1, and λ2 satisfy the following relations:
(30)λ1+λ2=L1+M2,λ1λ2=L1M2-L2M1.Now, we write the spectral decomposition of matrixA(s) as
(31)A(s)=λ1E+λ2F,
where E and F are called the projectors of A(s) and satisfy the following conditions(32a)E+F=I,(32b)EF=FE=O,(32c)E2=E,F2=F.Then, we have(33)A(s)=λ1E+λ2F,
where
(34)E=-1λ1-λ2[λ2-L1-L2-M1λ2-M2],F=1λ1-λ2[λ1-L1-L2-M1λ1-M2],I=[1001],O=[0000].Finally, we get(35)B(s)=A(s)=1λ1+λ2[λ1λ2+L1L2M1λ1λ2+M2].The characteristic equation of matrixB(s) can be written as
(36)k2-k(λ1+λ2)+λ1λ2=0,
where the roots of (36), namely, k1 and k2, can be written as
(37)k1=λ1,k2=λ2.Now, the Taylor series expansion ofexp[-A(s)x] yields
(38)exp[-A(s)x]=exp[-B(s)x]=∑n=0∞‍[-B(s)x]nn!.Using the Cayley-Hamilton theorem, we can expressB2 and higher orders of matrix B in terms of I and B, where I is the identity matrix of second order, Bahar and Hetnarski [39].Therefore, the infinite series in (38) can be expressed as
(39)L(x,s)=exp[-B(s)x]=b0I+b1B(s),
where b0 and b1 are coefficients depending on x and s. By the Cayley-Hamilton theorem, the characteristic roots k1 and k2 of the matrix B must satisfy (39), so we have
(40)exp[-k1x]=b0+b1k1,exp[-k2x]=b0+b1k2.By solving the above system of equations and using (37), we get
(41)b0=λ1e-λ2x-λ2e-λ1xλ1-λ2,b1=e-λ1x-e-λ2xλ1-λ2.Plugging the values ofb0 and b1 in (39), we have
(42)exp[-B(s)x]=L(x,s)=[lij],i,j=1,2,
where the entries lij(x,s) are given as
(43)l11=(λ1-L1)e-λ2x-(λ2-L1)e-λ1xλ1-λ2,l12=L2(e-λ1x-e-λ2xλ1-λ2),l21=M1(e-λ1x-e-λ2xλ1-λ2),l22=(λ1-M2)e-λ2x-(λ2-M2)e-λ1xλ1-λ2.The solution of (25) can be written in the following form:
(44)V-(x,s)=[lij][V-(0,s)+A-1(s)C(0,s)]-A-1(s)C(x,s).Substituting the values ofV-(x,s), A-1(s), and C(x,s) into (44) and performing the necessary matrix operations, we obtain the values of θ-(x,s) and σ-(x,s) as
(45)θ-(x,s)=1λ1-λ2×[e-λ1x{(λ1-L1)(θ-0+η1)-L2(σ-0+η2)}e-λ2x-{(λ2-L1)(θ-0+η1)-L2(σ-0+η2)}e-λ1x]-η1,(46)σ-(x,s)=1λ1-λ2×[e-λ1x{(λ1-M2)(σ-0+η2)-M1(θ-0+η1)}e-λ2x-{(λ2-M2)(σ-0+η2)-M1(θ-0+η1)}×e-λ1x]-η2,
where
(47)η1=M2L3-L2M3L1M2-M1L2,η2=L1M3-M1L3L1M2-M1L2.Now, considering (12) along with (15), and the Laplace transform, the displacement component is evaluated as
(48)u-(x,s)=1α1s2∂σ-∂x.Substituting (46) into (48), we get
(49)u-(x,s)=1α1s2(λ1-λ2)×[λ1{(λ2-M2)(σ-0+η2)-M1(θ-0+η1)}e-λ1x-λ2{(λ1-M2)(σ-0+η2)-M1(θ-0+η1)λ2}e-λ2xλ1].


## 5. Application

Problem: the ramp type boundary temperature of an elastic half-space.We consider a homogeneous isotropic thermoelastic solid occupying the half-spacex≥0. The boundary of half-space x=0 is affected by ramp type heating. In mathematical notations, the boundary conditions can be denoted as
(50)σ(0,t)=σ0=-p,θ(0,t)=θ0=θ*h(t),
where θ* is constant temperature and h(t) is defined as
(51)h(t)={0,t≤0tt0,0<t≤t01,t>t0,
where t0 is ramping parameter.We employ nondimensional variables given in (15) on (50) and taking the Laplace transform, we get
(52)σ-(0,s)=σ-0=-δ1s,θ-(0,s)=θ-0=θ*(1-est0t0s2).Substituting the values ofσ-0 and θ-0 from (52) into (45)-(46) and (49), we find
(53)θ-(x,s)=1λ1-λ2×[{(λ1-L1)(θ*η*+η1)-L2(η2-δ1s)}×e-λ2x-{(η2-δ1s)(λ2-L1)(θ*η*+η1)-L2(η2-δ1s)}e-λ1x{(λ1-L1)(θ*η*+η1)-L2(η2-δ1s)}]-η1,σ-(x,s)=1λ1-λ2×[{(λ1-M2)(η2-δ1s)-M1(θ*η*+η1)}×e-λ2x-{(λ2-M2)(η2-δ1s)-M1(θ*η*+η1)(η2-δ1s)}e-λ1x(η2-δ1s)]-η2,u-(x,s)=1α1s2(λ1-λ2)×[λ1{(λ2-M2)(η2-δ1s)-M1(θ*η*+η1)(η2-δ1s)}e-λ1x-λ2{(λ1-M2)(η2-δ1s)-M1(θ*η*+η1)(η2-δ1s)}e-λ2x(η2-δ1s)],
where
(54)η*=1-e-st0t0s2.


## 6. Limiting Cases

(1)
The mathematical expressions for the field variables studied in the context of generalized magneto-thermoelasticity theory with initial stress under applied boundary condition can be obtained by applyingα=1 in (13).
(2)
Neglecting initial stress effect, by substitutingp=0 in (14) and η=1 in (58), we obtain the expressions of field variables under fractional order generalized magneto-thermoelasticity.
(3)
The expressions for studied fields in the context of the fractional order generalized thermoelasticity theory with initial stress can be deduced by settingH0=0 in (12).


## 7. Numerical Inversion of the Laplace Transforms

We will now outline the numerical inversion method used to find the solution in the physical domain. The inversion formula of the Laplace transform is defined as(55)f(t)=12πι∫c-ι∞c+ι∞‍estf-(s)ds,
where f-(s) is the Laplace transform of function f(t).In order to invert the Laplace transforms in the equations given in Section5, we apply a numerical inversion method based on the Fourier series expansion explained by Honig and Hirdes [40]. In this method, the inverse transform f(t) of the Laplace transform f-(s) is approximated by the relation as
(56)f(t)=ectt1[12f-(c)+Re∑k=1N‍eιkπt/t1f-(c+ιkπt1)],0≤t<2t1,
where N is a sufficiently large integer representing the number of terms in the truncated Fourier series chosen such that
(57)ectRe[eιNπt/t1f-(c+ιNπt1)]≤ε1,
where ε1 is a prescribed small positive value that corresponds to the degree of accuracy to be achieved and c is a positive constant and must be greater than the real parts of all the singularities of f-(s). The optimal choice of c was obtained according to the criteria described by Honig and Hirdes [40].


## 8. Numerical Results and Discussions

To illustrate and compare the theoretical results obtained in the Section5, we now present some numerical results which depict the variations of displacement, temperature, and stress component. The material chosen for the purpose of numerical evaluations is copper, for which we take the following values of the different physical constantsE
=
36.9
×
1
0
10 kgm−1 s−2, ε0=(10-9/36π) Fm−1, k=386 Wm−1 K−1, T0=293 K, ρ=8954 kgm−3, τ0=0.02 s, v=0.33, αt=1.78×10-5 K−1, cE=383.1 Jkg−1 K−1, H0=(107/4π) Am−1, μ0=4π10-7 Hm−1.The general Lame constantsλ- and μ- are given as
(58)λ-=Evη(1+v)(1-2v),μ-=E2η(1+v),
where η is the initial stress parameter, E is Young's modulus; and v is Poisson ratio. For isotropic elastic medium with no initial stress, we take η=1.The computations are carried out fort=0.1, t0=0.5, α=0.5, and θ*=1. The numerical technique, outlined in previous section, was used to invert the Laplace transform in (53), providing the displacement, temperature, and stress distributions in the physical domain. The results are represented graphically for different positions of x.Figures1, 2, and 3 exhibit the space variations of the field quantities in the context of fractional order theory of thermoelasticity with magnetic field and initial stress for different values of fractional parameter α.Figure 1
Displacement distributionu for different values of α at t=0.1.Figure 2
Temperature distributionθ for different values of α at t=0.1.Figure 3
Stress distributionσ for different values of α at t=0.1.Figure1 displays the variations of displacement component u for different values of α, and it is noticed that the magnitude of displacement component u decreases with the increase in the value of fractional parameter α. In both the cases (i.e., α=0.5 and α=1.0), the displacement component attains maximum value at the boundary of half-space, and then continuously decreases to zero. Hence, displacement component u has similar trend for both the values of α.Figure2 depicts the variations of temperature θ with distance x for different values of α, and it is noticed that in both the cases (i.e., α=0.5 and α=1.0), maximum value of θ is 2 which is on the boundary of half-space. We observe from the figure that the difference is negligible in the beginning, and with the increase in x, the difference is much pronounced up to x≤1.4; both the series approach to zero. The trends of both the series are alike only up to x=0.8.Figure3 shows the variations of stress component σ with distance x for different values of α. It is evident from the figure that both the series have similar trend, that is, first increases to a maximum value and then decreases to a minimum value. The value of σ increases with the increase of fractional parameter when 0≤x≤0.8, while the trend of change is adverse when 0.8≤x≤6. The difference is significant in the range 0≤x≤2.8.Figures4, 5, and 6 depict the variations of displacement component, temperature, and stress component for three different cases defined as (a) ISMT (initially stressed magneto-thermoelasticity), η=2.5, p=1; (b) MT (magneto-thermoelasticity), η=1, p=0; and (c) IST (initially stressed thermoelasticity), η=2.5, p=1, H0=0. In all the above cases (a), (b), and (c), the value of fractional order parameter is taken to be α=0.5.Figure 4
Displacement distributionu at α=0.5, t=0.1.Figure 5
Temperature distributionθ at α=0.5, t=0.1.Figure 6
Stress distributionσ at α=0.5, t=0.1.Figure4 shows the variations of displacement component under the three cases discussed in (a), (b), and (c). It is observed that the displacement component in case of ISMT and IST has maximum positive values but in case of MT, the displacement component has maximum negative value on the boundary of half-space. The trend for displacement component in case of ISMT and IST is same, that is, both the series are continuously decreasing to zero, but in case of MT, the series rapidly approaches to zero. Figure 4 also exhibits that in the absence of magnetic field (IST) and initial stress (MT) decreases the value of displacement. It is also observed that as compared to ISMT, the difference is significant in case of IST but much pronounced in case of MT.Figure5 displays the variations of temperature under the ISMT, MT, and IST theories. It is observed that the trend of the series in case of ISMT, MT, and IST is similar, and the difference is significant. It is found that the absence of initial stress (MT) increases the value of temperature component, but the absence of magnetic field (IST) decreases the value of temperature as compared to the general case (ISMT).Figure6 exhibits the variations of stress under ISMT, MT, and IST theories. We found that the behavior of stress component in all the three cases is alike. It also observed that as compared to ISMT theory the stress component has large values in IST theory, but has small values in MT theory. The difference is significant in both IST and MT theories compared with ISMT theory, but the difference is much pronounced in MT theory as compared to IST theory.Figures7, 8, and 9 exhibit the variations of displacement component, temperature, and stress component in the context of the fractional order theory of thermoelasticity with magnetic field and initial stress for different values of ramp parameter. The values of ramp parameter are taken as 0.1, 0.3, and 0.5, and the value of fractional order parameter α is taken to be 0.5 for all the three cases.Figure 7
Displacement distributionu for different values of ramp parameter at α=0.5, t=0.1.Figure 8
Temperature distributionθ for different values of ramp parameter at α=0.5, t=0.1.Figure 9
Stress distributionσ for different values of ramp parameter at α=0.5, t=0.1.Figure7 demonstrates that as we increase the value of ramp parameter, then the value of displacement component also increases. It is also observed that trend for all the cases is similar and effect of change in ramp parameter is significant in first two cases and much pronounced in last two cases.Figure8 depicts that as we increase the value of ramp parameter, then the value of the temperature rapidly decreases. It is also observed that trend for the series in all the three cases is similar and the effect of change in ramp parameter is much pronounced in the initial range of distance 0≤x≤1.0.Figure9 displays that as we increase the value of ramp parameter, then the value of the stress component decreases. It is also observed that trend for the series in all the three cases is similar (i.e., all the series have same initial value, and first increases to a maximum value then decreases to a minimum value) and the effect of change in ramp parameter is much pronounced. It is also apparent from the figure that as we increase the value of ramp parameter from 0.1 to 0.3, the values of stress component decrease rapidly as compared to increase in the value of ramp parameter from 0.3 to 0.5.


## 9. Summary

We consider a perfectly conducting elastic homogeneous half-space in the context of fractional order generalized thermoelasticity theory with magnetic field and initial stress. The method of the matrix exponential, which constitutes the basis of the state space approach of modern theory, is applied to the nondimensional equations. The importance of state space approach is recognized in the fields, where the time behaviour of physical process is of interest. The state space approach is more general than the classical Laplace and Fourier transform technique. Consequently, state space is applicable to all systems that can be analyzed by integral transform in time and also is applicable to many systems for which transform theory breaks down [41]. The potential function approach is often used to solve problems of thermoelasticity theory. This, however, has several disadvantages as outlined in [39]. These may be summarized in the fact that the boundary conditions for physical problems are related directly to the physical quantities under consideration not to the potential functions. Secondly, more stringent assumptions must be made on the behaviour of potential functions than on the actual physical quantities. Last of all, it was found that many integral representations of physical quantities are convergent in the classical sense, while their potential function representations only converge in the mean. All these reasons have led many authors to avoid the use of potential functions. Among the alternatives is the state space formulation. This approach enables one to use the methodology of modern control theory in solving problems of thermoelasticity.The main conclusions due to the influence of magnetic field, initial stress, fractional parameter, and ramp parameter can be summarized as follows:(1)
the phenomenon of finite speed of propagation is preserved for all the field variables except for stress component due to the presence of initial stress, and all results are in agreement with the generalized theory of thermoelasticity;(2)
the fractional parameter has a significant effect on all the studied fields. The displacement component decreases with the increase in the value of the fractional parameter;(3)
the thermodynamic temperature first decreases then increases, but the stress component first increases then decreases to a minimum value with the increase in the value of fractional parameter;(4)
the effect of magnetic field is silent as compared to initial stress on studied fields. Also, all the field variables for the cases ISMT, MT, and IST behave alike except the displacement component for MT;(5)
in the context of fractional order generalized thermoelasticity theory with magnetic field and initial stress, the increase in the value of ramp parameter increases the magnitude of displacement component but decreases the magnitude of temperature distribution and stress component;(6)
displacement componentu, temperature distribution θ, and stress component σ show almost similar pattern for different values of fractional parameter α and ramp parameter t0.



---

*Source: 489863-2013-04-09.xml*