# Influence of Rotation and Magnetic Fields on a Functionally Graded Thermoelastic Solid Subjected to a Mechanical Load

Authors: Ankush Gunghas; Rajesh Kumar; Sunita Deswal; Kapil Kumar Kalkal
Journal: Journal of Mathematics
 (2019)

---

## Abstract

The current manuscript is presented to study two-dimensional deformations in a nonhomogeneous, isotropic, rotating, magneto-thermoelastic medium in the context of Green-Naghdi model III. It is assumed that the functionally graded material has nonhomogeneous mechanical and thermal properties in thex-direction. The exact expressions for the displacement components, temperature field, and stresses are obtained in the physical domain by using normal mode technique. These are also computed numerically for a copper-like material and presented graphically to observe the variations of the considered physical variables. Comparisons of the physical quantities are shown in figures to depict the effects of angular velocity, nonhomogeneity parameter, and magnetic field.


---

## Body

## 1. Introduction

In the classical dynamical coupled theory of thermoelasticity formulated by Biot [1], thermal signals are propagating with infinite speed, which is not physically acceptable. To remove this drawback of the classical coupled dynamical theory of thermoelasticity, several theories of generalized thermoelasticity were developed. The first two generalized thermoelastic theories are Lord-Shulman (L-S) [2] theory and Green-Lindsay (G-L) [3] theory. In L-S theory, one thermal relaxation time parameter is introduced in the classical Fourier’s law of heat conduction, whereas in the G-L theory, two thermal relaxation times are introduced in the constitutive relations for force stress tensor and entropy equation. Green and Naghdi [4–6] speculated three new thermoelasticity theories that permit treatment of a much wider class of heat flow problems, labeled as G-N models I, II, and III. The nature of the constitutive equations in these three models is such that when the respective theories are linearized, the model I is same as the classical heat conduction theory, model II predicts the finite speed of heat propagation involving no energy dissipation, and model III indicates the propagation of thermal signals with finite speed.The theory of magneto-thermoelasticity is concerned with the effect of magnetic field on elastic and thermoelastic deformations of a solid body and has received the attention of many researchers due to its extensive use in various fields like optics, geophysics, and acoustics. The problem of distribution of thermal stresses and temperature in a perfectly conducting half-space, in contact with a vacuum, permeated by an initial magnetic field was studied by Ezzat [7]. A model of two-dimensional equations of generalized magneto-thermoelasticity with one relaxation time in a perfectly conducting medium was established by Ezzat and Othman [8]. The impacts of fractional order parameter, hydrostatic initial stress, and gravity field on the plane waves in a fiber-reinforced isotropic thermoelastic medium were investigated by Othmanet al. [9]. Deswal and Kalkal [10] employed Laplace and Fourier transforms technique to study the phenomenon of wave propagation in a fractional order micropolar magneto-thermoelastic half-space. In the frame of fractional order theory of generalized thermoelasticity, Deswalet al. [11] studied magneto-thermoelastic interactions in an initially stressed, isotropic, homogeneous half-space. The effects of initial stress and magnetic field on thermoelastic interactions in an isotropic, thermally and electrically conducting half-space, whose surface is subjected to mechanical and thermal loads, were explored by Othman and Eraki [12]. Xiong and Guo [13] investigated the electro-magneto-thermoelastic diffusive plane waves in a half-space with variable material properties under fractional order thermoelastic theory.Since the large bodies like the earth, the moon, and other planets have an angular velocity, it appears more realistic to study the thermoelastic problems in a rotating medium. The propagation of elastic waves in a rotating, homogeneous, and isotropic medium was investigated by Schoenberg and Censor [14]. Some results in thermoelastic rotating medium are due to Roy Choudhuri and Debnath [15] and Roy Choudhuri and Mukhopadhyay [16]. The effect of rotation in generalized thermoelastic solid under the influence of gravity with an overlying infinite thermoelastic fluid was analyzed by Ailawalia and Narah [17]. Abouelregal and Zenkour [18] used fractional order theory of thermoelasticity to scrutinize the effect of angular velocity on fiber-reinforced generalized thermoelastic medium whose surface is subjected to a Mode-I crack problem. Kumaret al. [19] discussed the propagation of plane waves at the free surface of thermally conducting micropolar elastic half-space with two temperatures. Othmanet al. [20] considered the dual-phase lag model to study the influence of the rotation on a two-dimensional problem of micropolar thermoelastic isotropic medium with two temperatures. Saidet al. [21] used normal mode technique to study thermodynamical interactions in a micropolar magneto-elastic medium with rotation and two-temperature. Abouelregal and Abo-Dahab [22] investigated a two-dimensional problem in the context of dual-phase-lag model with fiber-reinforcement and rotation using normal mode analysis. The effect of angular velocity on Rayleigh wave propagation in a fiber-reinforced, anisotropic magneto-thermo-viscoelastic media was discussed by Hussien and Bayones [23].Over the last few decades, some structural materials such as functionally graded materials have been rapidly developed and used in many engineering applications. In functionally graded materials (FGMs), material properties vary gradually with a location within the body. FGMs are usually designed to be used under high-temperature environments. So, FGMs can easily eliminate or control thermal stresses, when sudden heating or cooling happens. These types of material are broadly used in important structures such as body materials in the aerospace field and nuclear reactors. A thermoinelastic response of functionally graded composites was studied by Aboudiet al. [24]. Abd-Allaet al. [25] analyzed radial vibrations in a functionally graded orthotropic elastic half-space subjected to rotation and gravity field. The electro-magneto-thermoelastic response of an infinite functionally graded cylinder was studied by Abbas and Zenkour [26], by using finite element method. The problem of generalized thermoelasticity in a thick-walled functionally graded cylinder with one relaxation time was considered by Abbas [27]. In this problem, the effects of temperature-dependent properties, volume fraction parameter, and thermal relaxation time on thermophysical quantities are estimated.The aim of the present contribution is to consider two-dimensional disturbances in an infinite, isotropic, nonhomogeneous, rotating, magneto-thermoelastic medium in the context of G-N model III. All the mechanical and thermal properties of the FGM under consideration are supposed to vary as an exponential power of the space-coordinate. The numerical results for the physical quantities have been obtained for a copper-like material and presented graphically to estimate and highlight the effects of different parameters considered in this problem.


## 2. Basic Equations

Following Green-Naghdi [5] and Roy Choudhuri and Debnath [15], the field equations and stress-strain-temperature relations in a rotating thermoelastic medium in the presence of body forces Fi areConstitutive Law (1) σ i j = 2 μ e i j + δ i j λ e r r - β θ ,where(2)eij=12ui,j+uj,i.Stress Equation of Motion (3) σ j i , j + F i = ρ ∂ 2 u i ∂ t 2 + Ω → × Ω → × u → i + 2 Ω → × u → ˙ i .Equation of Heat Conduction (4) K ∗ ∇ 2 θ + K ∇ 2 θ ˙ = ρ C e θ ¨ + β T 0 e ¨ ,where λ, μ are Lame’s elastic constants, β=(3λ+2μ)αt, αt is the coefficient of linear thermal expansion, σij are the components of stress, eij are the components of strain, ρ is the reference mass density, Ce is the specific heat at constant strain, K is the thermal conductivity, K∗ is the material constant characteristic for this theory, u→ is the displacement vector, θ=T-T0, T is the absolute temperature, T0 is the reference temperature of the medium in its natural state assumed to be θ/T0≪1, e=err is the cubical dilatation, Ω→ is the rotation vector, and δij is the Kronecker delta.For a nonhomogeneous medium, the parametersλ, μ, β, K, K∗, and ρ are no longer constant but become space-dependent. Hence we replace λ, μ, β, K, K∗, and ρ by λ0f(x→), μ0f(x→), β0f(x→), K0f(x→), K0∗f(x→), and ρ0f(x→), respectively, where λ0, μ0, β0, K0, K0∗, and ρ0 are supposed to be constants and f(x→) is a given nondimensional function of the space variable x→=(x,y,z). Using these values of parameters, (1), (3), and (4) take the following form:(5)σij=fx→2μ0eij+δijλ0err-β0θ,(6)σji,j+Fi=ρ0fx→∂2ui∂t2+Ω→×Ω→×u→i+2Ω→×u→˙i,(7)K0∗fx→θ,i,i+K0fx→θ˙,i,i=fx→ρ0Ceθ¨+β0T0e¨,where(8)β0=3λ0+2μ0αt.Here, the superposed dot denotes derivative with respect to time and the comma denotes derivative with respect to space variable.


## 3. Mathematical Model

Consider a nonhomogeneous, isotropic, magneto-thermoelastic half-space under the purview of G-N model III. Rectangular Cartesian coordinates are introduced having the surface of the half-space as the planex=0, with x-axis pointing vertically downwards into the medium. The medium is rotating with an angular velocity Ω→. Thus the displacement equation of motion in the rotating plane has two extra terms: Ω→×(Ω→×u→), which is the centripetal acceleration due to time varying motion only and 2Ω→×u→˙, which is the Coriolis acceleration. The present formulation is restricted to xy-plane and thus all the field variables are independent of the space variable z. So the displacement vector u→ and angular velocity Ω→ will have the components:(9)u→=u1,u2,0,Ω→=0,0,Ω.It is also assumed that material properties are graded only in x-direction. So we take f(x→) as f(x). By virtue of (9), the stresses arising from (5) can be expressed as(10)σxx=fxλ0+2μ0∂u1∂x+λ0∂u2∂y-β0θ,(11)σyy=fxλ0+2μ0∂u2∂y+λ0∂u1∂x-β0θ,(12)σxy=fxμ0∂u1∂y+∂u2∂x.Due to the application of an initial magnetic field H→(0,0,H0), an induced magnetic field h→, an induced electric field E→, and a current density J→ are developed in the considered medium. The simplified linear equations of electrodynamics of a slowly moving medium for a nonhomogeneous, isotropic and thermally conducting elastic solid are given by, neglecting Thomson’s effect [28]:(13)∇×h→=J→,∇×E→=-μm∂h→∂t,∇·h→=0,∇·E→=0,(14)E→=-μm∂u→∂t×H→,h→=∇×u→×H→,where μm is the magnetic permeability of the medium.From the above expressions, one can obtain(15)E→=-μmH0∂u2∂t,μmH0∂u1∂t,0,h→=0,0,-H0∂u1∂x+∂u2∂y,J→=∂h3∂y,-∂h3∂x,0,where(16)h3=-H0∂u1∂x+∂u2∂y.By virtue of above expressions and replacing μm by μm0f(x), the components of the Lorentz force F→=μm(J→×H→) are given by(17)F1=fxμm0H02∂e∂x,F2=fxμm0H02∂e∂y,F3=0.Utilizing the components of Lorentz force into stress equation of motion along with the consideration of two-dimensional problem, the field equations (6) and (7) yield(18)fxλ0+2μ0∂2u1∂x2+λ0+μ0∂2u2∂x∂y+μ0∂2u1∂y2-β0∂θ∂x+∂fx∂xλ0+2μ0∂u1∂x+λ0∂u2∂y-β0θ+fxμm0H02∂e∂x=fxρ0∂2u1∂t2-Ω2u1-2Ω∂u2∂t,(19)fxλ0+2μ0∂2u2∂y2+λ0+μ0∂2u1∂x∂y+μ0∂2u2∂x2-β0∂θ∂y+μ0∂fx∂x∂u1∂y+∂u2∂x+fxμm0H02∂e∂y=fxρ0∂2u2∂t2-Ω2u2+2Ω∂u1∂t,(20)K0∗fx∇2θ+∂fx∂x∂θ∂x+K0fx∇2θ˙+∂fx∂x∂θ˙∂x=fxρ0Ceθ¨+β0T0e¨,where(21)∇2=∂2∂x2+∂2∂y2.


## 4. Exponential Variation of Nonhomogeneity

By assumingf(x)=e-nx, where n is a dimensionless parameter, one can conclude that the mechanical and thermal properties of the material vary exponentially along the x-direction. The governing equations can be recast in the dimensionless form by introducing the following dimensionless parameters:(22)x′,y′=w∗c1x,y,u1′,u2′=ρ0w∗c1β0T0u1,u2,t′=w∗t,θ′=θT0,σij′=σijβ0T0,Ω′=Ωw∗,where(23)w∗=ρ0Cec12K0,c12=λ0+2μ0ρ0.Now, in terms of the dimensionless parameters given in (22), (10)-(12) and (18)-(20) transform to(24)σxx=e-nx∂u1∂x+A11∂u2∂y-θ,(25)σyy=e-nx∂u2∂y+A11∂u1∂x-θ,(26)σxy=e-nxA12∂u1∂y+∂u2∂x,(27)1+A22∂2u1∂x2+A21+A22∂2u2∂x∂y+A12∂2u1∂y2-∂θ∂x-n∂u1∂x+A11∂u2∂y-θ=∂2u1∂t2-Ω2u1-2Ω∂u2∂t,(28)1+A22∂2u2∂y2+A21+A22∂2u1∂x∂y+A12∂2u2∂x2-∂θ∂y-nA12∂u1∂y+∂u2∂x=∂2u2∂t2-Ω2u2+2Ω∂u1∂t,(29)∇2θ-n∂θ∂x+A31∇2θ˙-n∂θ˙∂x=A32∂θ˙∂t+A33∂e˙∂t,where(30)A11=λ0ρ0c12,A12=μ0ρ0c12,A21=A11+A12,A22=μm0H02ρ0c12,A31=K0w∗K0∗,A32=ρ0Cec12K0∗,A33=β02T0ρ0K0∗.


## 5. Solution Methodology

In this section, the normal mode method is employed, which has the advantage of finding the exact solutions without any assumed constraints on the field variables. In this approach, the solution of the physical variables is decomposed in terms of normal modes and one gets exact solution without any assumed restrictions on the actual physical quantities that appear in the governing equations of the problem considered. Normal mode analysis is, in fact, to look for the solution in Fourier transform domain. It is assumed that all the functions are sufficiently smooth on the real line such that the normal mode analysis of these functions exists. So, the solution for the considered physical variables can be decomposed in terms of normal modes in the following form:(31)u1,u2,θ,σijx,y,t=u1∗,u2∗,θ∗,σij∗xexp⁡ωt+ιmy,where u1∗, u2∗, θ∗, and σij∗ are the amplitudes of the functions, ω is the angular frequency, ι is the imaginary unit, and m is the wave number in y-direction.Introducing expression (31) in (27)-(29), we get(32)B11D2-nD+B12u1∗+B13D+B14u2∗-D-nθ∗=0,(33)B13D-B15u1∗+A12D2-B16D+B17u2∗+B18θ∗=0,(34)B19Du1∗+B21u2∗+B22D2-B23D-B24θ∗=0,where(35)B11=1+A22,B12=Ω2-A12m2-ω2,B13=A21+A22ιm,B14=2Ωω-nA11ιm,B15=2Ωω+nA12ιm,B16=nA12,B17=Ω2-ω2-m21+A22,B18=-ιm,B19=-A33ω2,B21=-B18B19,B22=1+A31ω,B23=nB22,B24=m2+A31m2ω+A32ω2.The condition for the existence of a nonzero solution of the system of (32)-(34) provides us(36)D6+F1D5+F2D4+F3D3+F4D2+F5D+F6u1∗,u2∗,θ∗=0,where(37)F1=E11E18+E12E17-A12E24-E14E23E11E17-A12E23,F2=E11E19+E12E18+E13F17-E25A12-E14E24-E15E23E11E17-A12E23,F3=E11E21+E12E19+E13E18-A12E26-E14E25-E15E24-E16E23E11E17-A12E23,F4=E11E22+E12E21+E13E19-E14E26-E15E25-E16E24E11E17-A12E23,F5=E12E22+E13E21-E15E26-E16E25E11E17-A12E23,F6=E13E22-E16E26E11E17-A12E23,E11=B11B18+B13,E12=-nB18+B14+nB13,E13=B12B18+nB14,E14=-B15+nA12,E15=B13B18+B16+nB15,E16=B14B18-nB16,E17=A12B22,E18=-B15B22+B23A12,E19=B16B22+B15B23-A12B24,E21=B15B24-B16B23,E22=-B16B24+B18B21,E23=B22B13,E24=-B13B23+B17B22,E25=B17B23-B13B24-B19B18,E26=B17B24.The general solution of (36) which is bounded as x→∞ is given by(38)u1∗,u2∗,θ∗x=∑i=131,H1i,H2iMim,ωe-λix,forRe⁡λi>0,where Mi(m,ω)(i=1,2,3) are parameters, depending upon m and ω, and(39)H1i=E11λi2-E12λi+E13A12λi3-E14λi2+E15λi-E16,H2i=B21H1i-B19λiB24-B23λi-B22λi2.In view of solution (38), stress components (24)-(26) take the form(40)σxx∗,σxy∗,σyy∗x=∑i=13H3i,H4i,H5iMim,ωe-λix-nx,forRe⁡λi>0,where(41)H3i=-λi+ιmA11H1i-H2i,H4i=A12ιm-λiH1i,H5i=-A11λi+ιmH1i-H2i.


## 6. Application: Mechanical Load on the Surface of the Half-Space

A nonhomogeneous, rotating, magneto-thermoelastic medium, occupying the half-spacex≥0, has been considered. The surface x=0 of the half-space is acted upon by a mechanical load as shown in Figure 1. So, the boundary conditions are given by the following.Figure 1
Mechanical load over a functionally graded rotating magneto-thermoelastic medium.(i) Mechanical Boundary Conditions (i)
Normal stress component obeys(42)σxx0,y,t=-py,t, wherep(y,t) is a given function of y and t. (ii)
Tangential stress component vanishes at the surfacex=0,i.e.,(43)σxy0,y,t=0.(ii) Thermal Boundary ConditionSince, the plane boundary surfacex=0 is taken to be isothermal, so the thermal boundary condition is the vanishing of temperature θ,i.e.,(44)θ0,y,t=0.Application of nondimensional parameters and normal mode technique defined in (22) and (31) respectively, transforms the above boundary conditions to the form:(45)σxx∗x=-p∗,σxy∗x=0,θ∗x=0atx=0.Taking into account the nondimensional expressions for temperature and stresses from (38) and (40), the above boundary conditions reduce to a nonhomogeneous system of three equations, which can be written in matrix form as(46)H31H32H33H41H42H43H21H22H23M1M2M3=-p∗00.Solution of system (46) provides us the values of Mi(i=1,2,3) as follows:(47)M1=Δ1Δ,M2=Δ2Δ,M3=Δ3Δ,where(48)Δ=H31L1-H32L2+H33L3,Δ1=-p∗L1,Δ2=p∗L2,Δ3=-p∗L3,L1=H42H23-H22H43,L2=H41H23-H21H43,L3=H41H22-H21H42.Substitution of (47) into expressions (38) and (40) provides us the following expressions of field variables(49)u1∗,u2∗,θ∗x=1Δ∑i=131,H1i,H2iΔie-λix,forRe⁡λi>0,(50)σxx∗,σxy∗,σyy∗x=1Δ∑i=13H3i,H4i,H5iΔie-λix-nx,forRe⁡λi>0.


## 7. Notable Cases

### 7.1. Neglecting Rotational Effect

In the absence of rotation (i.e., Ω=0), we shall be left with the relevant problem in a nonhomogeneous, isotropic, magneto-thermoelastic medium in the context of GN theory III. In this limiting case, we get the corresponding expressions of the physical quantities from (49) and (50).


### 7.2. Neglecting Nonhomogeneity Effect

By settingn=0 in (24)-(29), one can get required expressions for different distributions from (49) and (50). In this limiting case, our results coincide with those of Abo-Dahabet al. [29] with appropriate changes in loading and boundary conditions.


## 7.1. Neglecting Rotational Effect

In the absence of rotation (i.e., Ω=0), we shall be left with the relevant problem in a nonhomogeneous, isotropic, magneto-thermoelastic medium in the context of GN theory III. In this limiting case, we get the corresponding expressions of the physical quantities from (49) and (50).


## 7.2. Neglecting Nonhomogeneity Effect

By settingn=0 in (24)-(29), one can get required expressions for different distributions from (49) and (50). In this limiting case, our results coincide with those of Abo-Dahabet al. [29] with appropriate changes in loading and boundary conditions.


## 8. Numerical Results and Discussion

With an aim to illustrate the obtained theoretical results in the preceding section, we now present some numerical results. The following relevant physical constants are taken from Abo-Dahabet al. [29] for a copper-like material:(51)λ0=7.76×1010Nm-2,μ0=3.86×1010Nm-2,K0=0.6×10-2Wm-1K-1,Ce=383.1Jkg-1K-1,ρ0=8954kgm-3,T0=293K,αt=1.78×10-5K-1.Sinceω is complex quantity, we can write ω=ω0+ιω1 so that eωt=eω0t[cos(ω1t)+ιsin(ω1t)]. So for small values of time we can assume ω as real (i.e.,ω=ω0). The other parameters for numerical computation are taken as ω=2, m=2, t=0.1, p∗=8, and y=0.2.Figures2–5 analyze the effect of rotation on the distribution of field variables by considering three different values of angular velocity as Ω=0.0 (solid line), Ω=0.5 (dashed line), and Ω=0.9 (dot-dashed line) with n=1.0 and H0=105. Figure 2 explains the spatial variation of the normal displacement component for different values of Ω. The figure shows that the distribution of normal displacement follows a similar trend for all the values of Ω and dissimilarity lies on the ground of magnitudes. Figure 3 is plotted to depict the variation of normal stress with location x for three different values of Ω. The figure shows that increase in the value of Ω results in an increase in the numerical values of normal stress. Therefore, angular velocity is having an increasing effect on the profile of normal stress. Variations in tangential stress distribution with spatial coordinate x have been displayed in Figure 4. These variations are having a common starting point of zero magnitude, which is in quite good agreement with the boundary conditions. The figure shows that the tangential stress increases in the beginning and starts decreasing near the point x=0.23 and thereafter converges to zero as x increases. Moreover, with increasing Ω, there is an increase in the magnitude of tangential stress distribution. Figure 5 is drawn to observe the effect of angular velocity on the pattern of temperature distribution. As expected, the temperature distribution is having a coincident starting point of zero magnitude for all the values of Ω, which agrees completely with the boundary conditions. It is also manifested from the figure that increasing values of Ω are having a decreasing effect on the magnitude of temperature variations.Figure 2
Effect of angular velocity on normal displacement distribution.Figure 3
Effect of angular velocity on normal stress distribution.Figure 4
Effect of angular velocity on tangential stress distribution.Figure 5
Effect of angular velocity on temperature distribution.Figures6–9 illustrate the effect of the nonhomogeneity parameter on the distribution of field variables by setting three different values of nonhomogeneity parameter as n=0.0 (solid line), n=0.5 (dashed line), and n=1.0 (dot-dashed line) with Ω=0.5 and H0=105. In Figure 6, we have shown the spatial variation of normal displacement for different values of n. From this figure, it is noted that all the curves have distinct starting points. It can also be noticed from the plot that the displacement distribution is strongly affected by the presence of nonhomogeneity parameter. Figure 7 illustrates the variation of normal stress with distance x for different values of nonhomogeneity parameter n. For three different values of nonhomogeneity parameter, σxx starts with value 8.9 (in magnitude). Normal stress distribution exhibits significant sensitivity towards the nonhomogeneity parameter and it is compressive for homogeneous medium. Figure 8 is plotted to show the variations of tangential stress σxy with distance x. The plot indicates that the tangential stress field is having a coincident starting point of zero magnitude for all the three cases, which signifies that the boundary conditions are satisfied. It can be seen from the plot that, for homogeneous medium, the behaviour of σxy is totally opposite to that in the nonhomogeneous medium. The difference in magnitudes becomes indistinct along with the passage of time. Figure 9 shows that the temperature starts with a value zero which is completely in agreement with the boundary conditions. It increases in the beginning and starts decreasing in the neighbourhood of x=0.3 and converges to zero as x increases. It is also manifested from the figure that the presence of nonhomogeneity is having a decreasing effect on the magnitude of temperature distribution.Figure 6
Effect of nonhomogeneity parameter on normal displacement distribution.Figure 7
Effect of nonhomogeneity parameter on normal stress distribution.Figure 8
Effect of nonhomogeneity parameter on tangential stress distribution.Figure 9
Effect of nonhomogeneity parameter on temperature distribution.The influence of magnetic field on various field variables is examined in Figures10–13. Figure 10 shows the transient effect of applied mechanical load on normal displacement distribution, in the medium with a magnetic field (solid line) and without magnetic field (dashed line). Solution curves for both the cases follow a similar pattern of variations with the difference in magnitudes. It can be noted that magnetic field has a noticeable impact on displacement distribution. Figure 11 describes the variation of normal stress with location x. The plot shows that stress σxx is having a coincident initial point for both the cases and diminution of the magnitude takes place as the distance from the boundary increases. Figure 12 represents the spatial variation of tangential stress. It can be seen from the plot that the values of tangential stress for both the cases increase in the beginning and start decreasing in the neighbourhood of x=0.2 and approach to zero with increasing x. Figure 13 displays the variation of temperature field with distance x. In this figure, all the curves have coincident beginning point with value zero that leads to satisfy the boundary conditions. The amplitude of vibrations of solution curves for temperature field gets suppressed in the absence of magnetic field.Figure 10
Effect of magnetic field on normal displacement.Figure 11
Effect of magnetic field on normal stress distribution.Figure 12
Effect of magnetic field on tangential stress distribution.Figure 13
Effect of magnetic field on temperature distribution.The 3D plots showing normal displacement distribution, normal stress distribution, tangential stress distribution, and temperature distribution are shown in Figures14–17 for a wide range of x(0≤x≤4) and for a wide range of dimensionless time t(0≤t≤0.4). Figure 14 describes the variation of normal displacement with distance x and with time t. From this figure, it can be seen that the dimensionless time t plays an important role on the distribution of normal displacement. Figure 15 represents the variation in the values of normal stress for a wide range of x. Numerical values of normal stress decrease as the distance x increases, while for time t, an increase in the values takes place. Figure 16 has been plotted to show the profile of tangential stress distribution. The tangential stress starts with a zero value which is completely in agreement with the boundary conditions prescribed. Figure 17 represents the distribution of temperature with distance x and with time t. The temperature distribution is behaving like an increasing function in the range 0.0≤x≤0.2 and for the rest of the domain it decreases and reaches a steady state about the point x=2.7.Figure 14
Profile of normal displacement atn=1.0, H0=105, and Ω=0.5.Figure 15
Profile of normal stress distribution atn=1.0, H0=105, and Ω=0.5.Figure 16
Profile of tangential stress distribution atn=1.0, H0=105, and Ω=0.5.Figure 17
Profile of temperature distribution atn=1.0, H0=105, and Ω=0.5.


## 9. Concluding Remarks

The investigation under consideration provides a mathematical model to obtain the behaviour of normal displacement, stresses, and temperature in a nonhomogeneous, isotropic, rotating, magneto-thermoelastic medium within the framework of G-N model III, by using normal mode technique. Theoretical and numerical results reveal that the parameters, namely, rotation, nonhomogeneity parameter, and magnetic field, have significant effects on the considered physical variables. Analysis of graphs permits the following concluding remarks:(i)
As expected, the values of all the physical quantities converge to zero as the distancex increases and from the distribution of all physical quantities, it can be found that wave type heat propagates through the medium.(ii)
It is apparent from figures that the rotational speed has an increasing effect on the profiles of normal stress and tangential stress but it has a decreasing effect on the profile of temperature.(iii)
The presence of nonhomogeneity is having a decreasing effect on the magnitude of temperature distribution while it has a mix effect on the remaining field variables. Therefore, while designing FGMs, the effect of nonhomogeneity should be taken into consideration.(iv)
The presence of magnetic field plays a significant role in the distribution of all physical variables.(v)
The method adopted here is applicable to a wide range of problems in thermodynamics and thermoelasticity. It can be employed to boundary-layer problems which are described by linearized Navier-Stokes equations in electro-hydrodynamics.The above study is of geophysical interest and finds applications in mechanical engineering, industrial sectors, and seismology. The results presented in this paper will prove useful for scientists in material science and designers of new materials as well as for those working on the development of magneto-thermoelasticity. Electro-magneto composite materials have applications in sensors, actuators, ultrasonic imaging devices, and many other emerging components. Magneto-thermoelasticity has drawn the attention of many engineers, because of its wide use in diverse areas, especially, geophysics for determining the effect of earth’s magnetic field on seismic waves, development of a highly sensitive superconducting magneto-meter, electrical power engineering, etc. FGMs are broadly used in the biomedical application as medical implants, because they are designed to mimic the human organs, which are FGMs in nature. These types of materials are also used in pressure vessels and pipes.



---

*Source: 1016981-2019-06-10.xml*