Magneto-thermoelastic interactions in an initially stressed isotropic homogeneous elastic half-space with two temperatures are studied using mathematical methods under the purview of the L-S model of linear theory of generalized thermoelasticity. The formalism deals with the state space approach with the purpose of counteracting the difficulties of handling the displacement potential functions. Of specific concern here is the propagation of waves owing to ramp type increase in temperature and load. The medium is considered to be permeated by a uniform magnetic field. The expressions for different field parameters such as displacement, temperature, strain, and stress in the physical domain are obtained by applying a numerical inversion technique. Results of some earlier workers have been deduced from the present formulation. Numerical work is also performed for a suitable material with the aim of illustrating the results.

1. Introduction

The classical coupled thermoelasticity theory proposed by Biot [1] with the introduction of the strain-rate term in the Fourier heat conduction equation leads to a parabolic-type heat conduction equation, called the diffusion equation. This theory predicts finite propagation speed for elastic wave but an infinite speed for thermal disturbance. This is physically unrealistic. To overcome such an absurdity, generalized thermoelasticity theories have been propounded by Lord and Shulman [2] as well as Green and Lindsay [3] advocating the existence of finite thermal wave speed in solids. These theories have been developed by introducing one or two relaxation times in the thermoelastic process, either by modifying Fourier’s heat conduction equation or by correcting the energy equation and Neumann-Duhamel relation. According to these generalized theories, heat propagation can be visualized as a wave phenomenon rather than a diffusion one; in the literature, it is usually referred to as the second sound effect. These two theories are structurally different from one another, and one cannot be obtained as a particular case of the other. Various problems characterizing these theories have been investigated and have revealed some interesting phenomenon. Brief reviews of this topic have been reported by Chandrasekharaiah [4, 5].

The interplay of the Maxwell electromagnetic field with the motion of deformable solids is largely being undertaken by many investigators owing to the possibility of its application to geophysical problems and certain topics in optics and acoustics. Moreover, the earth is subject to its own magnetic field, and the material of the earth may be electrically conducting. Thus, the magnetoelastic nature of the earth’s material may affect the propagation of waves. Many authors have considered the propagation of electromagneto-thermoelastic waves in an electrically and thermally conducting solid. A comprehensive review of the earlier contributions to the subject can be found in the study by Paria [6]. Among the authors who considered the generalized magneto-thermoelastic equations are Nayfeh and Nemat-Nasser [7] who studied the propagation of plane waves in a solid under the influence of an electromagnetic field. They have obtained the governing equations in the general case and the solution for some particular cases. Choudhuri [8] extended these results to rotating media. Ezzat [9] has studied the problem of generation of generalized magneto-thermoelastic waves by thermal shock in a perfectly conducting half-space. Ezzat et al. [10] have established the model of two-dimensional equations of generalized magneto-thermoelasticity.

In dealing with classical or generalized thermoelastic problems in most situations, the displacement potential function approach is used. However, Bahar and Hetnarski [11] outlined several disadvantages of the potential function approach. These may be summarized in the fact that the boundary and initial conditions of the problem are not related directly to the potential function, as it has no physical meaning explicitly. Secondly, more stringent assumptions must be made on the behaviour of potential functions than on the actual physical quantities. Finally, 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. To get rid of these difficulties, Bahar and Hetnarski [11] introduced the state space formulation in thermoelastic problems. This state space approach has been further developed in Sherief [12] to include the effect of heat sources. He et al. [13] surveyed a one-dimensional thermal shock problem for a semi-infinite piezoelectric rod using state space approach. In a survey article, Ezzat [14] has discussed the developments in the theory of thermoelasticity and fluid mechanics. He has also reviewed the method of matrix exponential, which constitutes the basis of the state space approach in the same work. Youssef and El-Bary [15] put forward an analysis for a generalized thermoelastic infinite layer problem under three theories using state space approach. State space formulation to the vibration of gold nanobeam in femtoseconds scale was done by Elsibai and Youssef [16]. Recently, Deswal et al. [17] investigated the effects of magnetic field and initial stress on a fractional order generalized thermoelastic half-space.

The development of initial stresses in the medium is due to many reasons, for example, resulting from variation of temperature, process of quenching, overburden layer, and slow process of creep, gravitation, weight, largeness, and so forth. The earth is supposed to be under high initial stresses. It is therefore of great interest to study the influence of these stresses on the propagation of stress waves. It was the achievement of Biot [18] to show that acoustic propagation under initial stresses would be fundamentally different from that under stress free state and could not be represented merely by introducing stress dependent elastic coefficients in the classical theory. He has obtained the velocities of longitudinal and transverse waves propagating along the coordinate axes only. Many authors have contributed their efforts to study the wave propagation in solids under initial stresses for various models. The study of reflection and refraction phenomenon of plane waves in an unbounded medium under initial stresses is due to Chattopadhyay et al. [19]. Montanaro [20] constructed the model of isotropic linear thermoelasticity with hydrostatic initial stress. Othman and Song [21] and Singh [22] studied the reflection of thermoelastic waves from a free surface under hydrostatic initial stress in the context of different theories of generalized thermoelasticity.

The theory of heat conduction in a deformable body, formulated by Chen and Gurtin [23] and Chen et al. [24, 25] depends on two different temperatures: the conductive temperature and the thermodynamical temperature . Chen et al. [24] have suggested that the difference between these two temperatures is proportional to heat supply. In absence of heat supply, these two temperatures are identical for time independent situation. However, for time dependent cases, particularly for problems related to wave propagation, two temperatures are in general different, regardless of heat supply. The two-temperature thermoelasticity theory has gained much attention of the researchers in the recent years. The existence, structural stability, convergence, and spatial behaviour of two-temperature thermoelasticity have been provided by Quintanilla [26]. Youssef [27] has developed a new model of generalized thermoelasticity that depends on two temperatures and , where the difference between the two temperatures is proportional to heat supply with a nonnegative constant (length2). Later on, Youssef and Al-Lehaibi [28] and Youssef and El-Bary [29] investigated various problems on the basis of two temperature thermoelasticity with relaxation times and showed that the obtained results are qualitatively different as compared to those in case of one-temperature thermoelasticity.

It is worthwhile mentioning here that, in most of the earlier studies, mechanical or thermal loading on the bounding surface is considered to be in the form of a shock. But the sudden change in the load is definitely an idealized situation as it is difficult to interpret physically a pulse represented mathematically by a step function; even a very rapid rise-time (of the order of 10−9 second) may be slow in continuum terms. This is particularly true in the case of second sound effects when the thermal relaxation times for typical metals are less than 10−9 second. It is thus felt that a finite time of rise of external load applied on the bounding surface should be considered while studying a practical problem of this nature. Considering this vital aspect of rise of time, Misra et al. [30] and Youssef with many authors investigated many applications in which the ramp type heating is used [15, 16].

The present investigation is intended to demonstrate the use of state space approach in analyzing the propagation of magneto-thermoelastic waves with hydrostatic initial stress in two-temperature theory of thermoelasticity. Laplace transform technique is employed to obtain the general solution. The inverse Laplace transform is carried out numerically using a method based on Fourier expansion technique in Honig and Hirdes [31]. Numerical values of the field quantities are calculated by considering illustrative examples, and their variations with respect to space coordinate are displayed graphically. Some particular cases have also been deduced from the present analysis.

2. Formulation of the Problem

Let us consider a perfectly conducting thermoelastic half-space with hydrostatic initial stress and two temperatures of a linear, homogeneous, and isotropic medium whose state can be expressed in terms of the space variable and the time variable . The initial conditions are assumed to be homogeneous. The adjacent free space is assumed to be permeated by a uniform magnetic field acting parallel to the boundary . This produces an induced magnetic field and induced electric field which satisfy the linearized equations of electromagnetism and are valid for slowly moving media:

The above equations are supplemented by the field equations and constitutive relations with initial stress and two temperatures according to Lord and Shulman [2], Montanaro [20], and Youssef [27] as follows. (1) Constitutive relations: (2) Heat conduction equation: (3) Equation of motion: where (4) Relation between conductive temperature and thermodynamical temperature:

The displacement components for one-dimensional medium have the form

The strain component takes the form

The constitutive equation, heat conduction equation, and equation of motion assume the shape


Equation (10) can be rewritten as

The relation between the conductive temperature and thermodynamical temperature takes the form

In order to nondimensionalize (9), (12), and (13), let us define the following set of dimensionless variables: where and .

Introducing (14) in (9), (12), and (13), we obtain the following equations in nondimensional forms (dropping the primes): where , , , , and are the coupling parameters.

Equations (15) will change in the Laplace transform space to where , and denotes the Laplace transform variable (bar over a symbol indicates the quantity in Laplace transform domain).

Elimination of and from (16) and (17) yields where

Having chosen conductive temperature and stress component as the state variables, (18) may be recasted in the form where

The formal solution of the differential equation (20) can be expressed as


after cancelling the exponential having positive power to get bounded solution for large .

Now we have to find out the form of the matrix . The characteristic equation of the matrix is obtained as where and are the roots of (24) which satisfy the relations

Now, we can write the spectral decomposition of as follows, Charles [32] where is a matrix resolution of the identity; that is,

Then we have where

In particular, if , we have

Finally, we get

Now the Taylor series expansion of has the form

Using the well-known Cayley-Hamilton theorem, we can express and higher orders of the matrix in terms of and , where is the unit matrix of order 2. Thus, the infinite series in (32) can be truncated as where and are some parameters depending on and .

The characteristic equation of the matrix can be written as The roots of this equation, namely, and , take the forms

According to Cayley-Hamilton theorem, the characteristic roots and of the matrix must satisfy (33). Therefore, we get

The solution of this system of linear equations is given by

Plugging the above expressions along with and into (33), we obtain the components of the matrix as

Hence solution in (22) can be written as

Inserting the values of , , and into (39) and after straight-forward calculation, the expressions for conductive temperature and stress are evaluated as


Combination of (17) and (40) yields

3. Application

3.1. Problem I: Ramp Type Increase in Boundary Temperature of an Elastic Half-Space

We consider a homogeneous isotropic thermoelastic solid occupying the half-space . The boundary plane of the half space is subject to a ramp type increase in temperature. Mathematically, these boundary conditions can be written as where is a constant temperature and is a function of defined as where is a fixed moment of time.

Employing (14) and Laplace transform on (44) and (45), one can find

We can now use the values of and from (46) to finally get the solutions in the Laplace transform domain as where

Considering (10) along with (11) and using nondimensional variables and Laplace transform, the displacement component is evaluated as

Substituting from (48) into (52), we find

3.2. Problem II: Ramp Type Increase in Boundary Load of an Elastic Half-Space

In this case we assume that the plane boundary of the half-space is kept at the constant reference temperature and is subject to a ramp type increase in load. Mathematically, these boundary conditions can be written as where is a constant load and is a function as defined by (45). Laplace transform of the equations (54) containing boundary conditions yields (in non-dimensional form)

Adopting the same procedure, the expressions for all the considered physical quantities can be found as

4. Limiting Cases

4.1. Case I: With One Temperature

The expressions for different distributions in the above two sets of boundary conditions with one temperature, that is, the case when conductive temperature coincides with thermodynamical temperature, can be obtained from expressions (47)–(50), (53), and (56) by setting and consequently in the field equations.

4.2. Case II: Without Magnetic Effect

In order to discuss the problem of wave propagation in two-temperature thermoelasticity with initial stress, it is sufficient to set , and accordingly, we have . Taking into consideration the above-mentioned modifications, the corresponding expressions of the physical fields can be obtained from expressions (47)–(50), (53), and (56).

4.3. Case III: Without Initial Stress

The required expressions for the components of displacement, stress, strain, and temperature can be procured from the expressions (47)–(50), (53), and (56) in the absence of initial stress by putting and in the governing equation. In addition to this if we also neglect magnetic effect (i.e., ), then our results match with those obtained by Youssef and Al-Lehaibi [28].

5. Numerical Inversion of the Laplace Transform

In the previous section, we obtained the solutions of two problems in the Laplace transform domain . We shall now outline briefly the numerical inversion method used to find the solution in the physical domain. Let be the Laplace transform of a function . The inversion formula for Laplace transform can be written as Churchill [33] where is an arbitrary real number greater than all the real parts of the singularities of . Taking and using Fourier series in the interval , we obtain the approximate formula Honig and Hirdes [31] where

Two methods are used to reduce the total error. First, the “Korrecktur” method is used to reduce the discretization error. Next, the -algorithm is used to reduce the truncation error and hence to accelerate convergence. The details of these methods can be found in Honig and Hirdes [31]. The values of and are chosen according the criteria outlined in Honig and Hirdes [31].

6. Numerical Results and Discussions

With an aim to illustrate the contribution of initial stress, two temperatures, and magnetic field on the field quantities, a numerical analysis is carried out. The material chosen for the purpose of numerical computation is copper for which

The relation connecting Young’s modulus , the modulus of rigidity , and Poisson’s ratio for isotropic elastic solids does not hold well for real earthy materials, namely, sand, and soil. Weiskopf [34] has established a stress-strain relation for the same. He has investigated that due to the slipping of granules on each other, the resistance to shear is much less than in a solid and the resulting shearing deflection is much greater. For these materials . So the relation has been considered here. Also the other generalized Lame’s constant may be defined as where is defined as sandiness or initial stress parameter. corresponds to isotropic elastic medium with no initial stress.

The values of other parameters arising in the present analysis are taken as

With these numerical values of the parameters, the values of the field quantities are computed for the time . All the field quantities have been examined for the following four cases: (i) initially stressed magneto-thermoelastic solid with two temperatures (ISMTT, , , and ), (ii) initially stressed magneto-thermoelastic solid (ISMT, , , and ), (iii) magneto-thermoelastic solid with two temperatures (MTT, , , and ), and (iv) initially stressed thermoelastic solid with two temperatures (ISTT, , , , and ).

From application point of view, we have divided the graphs into three groups. The first group (Figures 14) represents the distribution of variables for the case of ramp type increase in boundary load. The variation of the field variables for the case of ramp type increase in boundary temperature is observed in second group (Figures 58). In this group, we have not considered the ISTT case (i.e., without magnetic effect) graphically keeping in view the resolution and clarity of plots. However, magnetic field has the same effect on the physical fields as in the case of ramp type increase in boundary load. The third group (Figures 912) displays the influence of ramp parameter on the considered physical variables for three values of which are taken as 0.1, 0.5, and 0.9.

6.1. Boundary Load

Figure 1 displays a comparison of displacement component in the context of four cases. The displacement field shows maximum negative value at the boundary of the half-space and then changes to zero value with increase in the spatial coordinate. Also the effects of two temperature, magnetic field, and initial stress parameter are very much pronounced. The two-temperature parameter has decreased the magnitude of the displacement field, while initial stress parameter and magnetic field have increased the magnitude. Figure 2 exhibits the spatial variation of the thermodynamical temperature . The thermodynamical temperature is enjoying the similar pattern in all the four cases except that the magnitudes in all the cases are different. In case of ISMT, it starts with zero value and further decreases continuously to a maximum negative value −1.8108 at . Afterwards, its value increases and attains the zero value. Significant difference is noticed in case of MTT. Also the two-temperature parameter is playing a vital role. It has both increasing and decreasing effects on the thermodynamical temperature. It can also be seen from the plot that magnetic field acts to increase the magnitude of thermodynamical temperature.

The variation of stress component with spatial coordinate has been observed in Figure 3. The stress component starts with positive values in all the cases, and it decreases continuously to attain its lowest value. Presence of initial stress, two temperature parameter, and magnetic field has appreciably increased the magnitude of the stress component. Another interesting observation is the nonvanishing behaviour of the stress component in ISMTT, ISMT, and ISTT cases, though the values are very small. This important phenomenon has happened owing to the presence of initial stress. Figure 4 is plotted to describe conductive temperature for all the four cases. The conductive temperature starts with a zero value for all cases and has decreasing behaviour for the interval . It attains its highest negative values −1.7235, −1.8108, −4.2750, and −1.6373 at for the ISMTT, ISMT, MTT, and ISTT cases, respectively. Initial stress parameter has strongly influenced the conductive temperature. It is also apparent from the figure that two temperature parameter and magnetic field have also produced significant difference in conductive temperature.

6.2. Boundary Temperature

Figure 5 is drawn to show the variation of the displacement component with distance . The displacement component starts with a positive maximum value and then decreases monotonically with attaining its lowest value at the location and then gradually diminishes to zero. Here the location at which the extreme values are attained is independent of the initial stress and two temperature parameter. In fact the effect of the stress and two temperature parameter is very tiny, and this effect vanishes with the non-dimensional location. Figure 6 represents the change in the thermodynamical temperature with location . The values of the thermodynamical temperature are found to be small in magnitude in MTT case. It can also be seen from the figure that the trend of thermodynamical temperature is similar in all the three cases.

Figure 7 gives the stress variation in the direction of wave propagation. At , the stress reduces to zero for the MTT case which agrees with the boundary condition. In this case also, the stress is having nonvanishing behaviour with very small magnitudes for the ISMTT and ISMT cases. The influence of initial stress parameter is very much pronounced on the stress component. In Figure 8, we have depicted the profile of conductive temperature with location . It is behaving like an increasing function in the entire range for all the three cases. Also the contribution of additional parameters like initial stress, two-temperature, and magnetic field reduces to zero with increase in the spatial coordinate.

6.3. Influence of Ramp Parameter (Boundary Load)

Figures 912 display the values of the considered physical variables in the direction of wave propagation for three values of ramp parameter which are taken to be 0.1, 0.5, and 0.9. We observe from these figures that all the variables behave in the same manner due to the change in the values of ramp parameter with some difference in magnitudes. Increase in the value of ramp parameter is making the magnitudes of the field variables small. Hence the ramp parameter has a salient effect on all the variables.

7. Conclusions

A mathematical study is made to discuss the problem of magneto-thermoelastic interactions in two-temperature theory with hydrostatic initial stress. The theoretical and numerical results reveal that all the four parameters, namely, initial stress, two temperature, magnetic field, and ramp parameter have significant effect on the considered physical variables. It is concluded the following.(1)The phenomenon of finite speed of propagation is manifested in all the figures except the stress component which is also having very small value.(2)With the increase in distance, the results are quite close to each other, which is in agreement with the generalized theories of thermoelasticity.(3)It is apparent from the figures that initial stress parameter strongly affects the physical quantities.(4)The nature of the source (boundary load or temperature) highly influences the physical quantities. The body is deformed to a large extent when a boundary load is applied. The displacement component exhibits increasing behaviour in case of boundary load, while in case of boundary temperature, it initially decreases to its maximum negative value and then gradually approaches to zero. Similarly, the other physical quantities behave significantly different when different type of sources are applied.(5)All the physical variables are observed to be sensitive (more or less) towards the two temperature parameter.(6)Presence of magnetic field has affected all the considered physical variables. It acts to increase the magnitude of all the field variables. (7)The ramp parameter has significantly affected all the studied fields. Increment in the value of ramp parameter causes decrement in the values of all the fields which is very obvious from the figures.(8)All the functions except stress vanish identically outside some finite interval. This is not the case in coupled thermoelasticity where the considered functions have non-vanishing values for all due to the infinite speed of propagation of heat waves.

This research demonstrates how state space approach can be extended to study the problem involving magneto-thermoelastic interactions in two temperature theory with initial stress which has been rarely employed in analyzing such problems. The inclusion of magnetic and initial stress parameters provides a more realistic model of the generalized thermoelasticity theory.


:Counterparts of Lame’s parameters
:Ramp parameter
:Initial pressure
:Initial stress parameter
:Young’s modulus
:Poisson’s ratio
:Two temperature parameter
:Coefficient of linear thermal expansion
:Thermodynamical temperature
:Temperature of the medium in its natural state assumed to be
:Conductive temperature
:Components of the stress tensor
:Components of the displacement vector
:Density of the medium
:Components of the strain tensor
:Cubical dilatation
:Specific heat at constant strain
:Thermal conductivity
:Thermal relaxation time
:Magnetic permeability
:Electric permittivity
:Lorentz force
:Kronecker delta function.