Document Type: Research Paper
Authors
^{1} University of Calcutta
^{2} Univresity of Calcutta
Abstract
Keywords

Mechanics of Advanced Composite Structures 6 (2019) 95 – 104


Semnan University 
Mechanics of Advanced Composite Structures journal homepage: http://MACS.journals.semnan.ac.ir 
ElastoThermodiffusive Response in a TwoDimensional Transversely Isotropic Medium
A. Sur ^{*}, M. Kanoria
Department of Applied Mathematics, University of Calcutta, 92 A. P. C. Road, Kolkata 700009, India
Paper INFO 

ABSTRACT 
Paper history: Received 20171226 Received in revised form 20180818 Accepted 20180923 
The present article investigates the elastothermodiffusive interactions in a transversely isotropic elastic medium in the context of thermoelasticity with one relaxation time parameter and two relation time parameters. The resulting nondimensional coupled equations are applied to a specific problem of a halfspace in which the surface is free of tractions and is subjected to timedependent thermal and chemical loadings. The analytical expressions for the displacement components, stresses, temperature, strain, mass diffusion, and chemical potential are obtained in the physical domain by employing the normal mode analysis as a tool. These expressions are calculated for a copperlike material and the results are depicted graphically. A comparative study of a diffusive medium and a thermoelastic medium show that diffusion has a significant effect on the thermophysical quantities. Furthermore, in the absence of the effect of thermodiffusion, the results agree with the existing literature.




Keywords: Generalized thermoelastic diffusion Finite wave speed Normal mode analysis Transversely isotropic material Green–Lindsay model 


© 2019 Published by Semnan University Press. All rights reserved. 
Thermoelastic diffusion, which is also known as elastothermodiffusion in elastic solids, deals with the coupling effects of the fields of temperature, mass diffusion, and strain, in addition to heat and mass exchange with the environment. It has extensive applications in geophysics and other industrial fields, including the extraction of oil from oil deposits. In recent years, the subject of thermoelastic diffusion has received serious attention. The theory of thermoelastic diffusion was first developed by Nowacki [1–4]. Gawinecki et al. [5] proved a theorem about the existence, uniqueness, and regularity of the solutions for a nonlinear parabolic thermoelastic diffusion problem. Gawinecki and Szymaniec [6] established a theorem about the global existence of the solution for the same problem. In the theory developed by Nowacki [1–4], the classical coupled thermoelastic model was used [7]. It should be mentioned that the theory of coupled dynamical thermoelasticity predicts an infinite speed for thermal signals, which is physically unrealistic. There has been an increased interest in the field of heat propagation to remove this unrealistic prediction, which has led to the development of wellestablished theories of generalized thermoelasticity. Generalized thermoelasticity theories involve hyperbolictype governing equations and predict the finite speed of thermal signals. For example, Lord and Shulman [8] proposed the generalized thermoelasticity theory, which is known as the LS model, involving one relaxation time. Green and Lindsay [9] developed the temperatureratedependent thermoelasticity model (GL model) involving two relaxation times.
Recently, Sherief et al. [10] developed a generalized thermoelastic diffusion theory with one relaxation time, which allowed waves to propagate at finite speeds. Sherief and Saleh [11] investigated the problem of a thermoelastic halfspace in the context of the generalized thermoelastic diffusion theory with one relaxation time. Singh [12] discussed the reflection wave phenomena from the free surface of an elastic solid with generalized thermodiffusion with one relaxation time and two relaxation times in a later study [13]. Aouadi [14] studied diffusion in an infinitely long solid cylinder and in an infinite elastic body with a spherical cavity [15]. The uniqueness and reciprocity theorems for the equations of a generalized thermoelastic diffusion problem in isotropic media were proven by Aouadi [16] on the basis of the Laplace transform method. Kumar and Gupta [17] studied the wave propagation at the boundary surface of an inviscid fluid under thermoelastic diffusion. Recently, Othman et al. [18] analyzed the effects of diffusion on a twodimensional problem of generalized thermoelasticity in the context of the Green–Naghdi theory. Deswal and Choudhary [19–21] also analyzed a twodimensional thermoelastic diffusion problem using the same theory. Kumar and Kansal [22] discussed the propagation of waves on the free surface of a transversely isotropic body under generalized thermoelastic diffusion. Kothari and Mukhopadhyay [23, 24] investigated thermoelastic diffusion inside a spherical shell under three different theories. Wang et al. [25] studied the thermoelastic dynamic solution of a multilayered spherically isotropic hollow sphere for spherically symmetric problems. Such a body is said to possess transverse isotropy about any radius vector drawn from the center to a given point of material. Recently, several researchers, including ElSayed [26], Karmakar and Kanoria [27], and Bhattacharya and Kanoria [28, 29], have used the elastothermodiffusive response to solve several problems. In addition, a few remarkable works on generalized thermoelastic diffusion have been published [30–33].
In the present analysis, we study the generalized thermoelastic diffusion in a transversely isotropic twodimensional thermoelastic medium subjected to a prescribed temperature and chemical loading in which the boundary is free of traction. The analysis compares the thermoelastic diffusion model with two relaxation times, also known as the Green–Lindsay model with diffusion (GLD), and the thermoelastic diffusion model with one relaxation time, also known as the Lord–Shulman model with diffusion (LSD) for a thermodiffusive medium. Introducing a normal mode analysis, the governing equations have been expressed and solved in terms of normal modes. The numerical estimates for the thermal stresses, temperature, mass concentration, and chemical potential have been computed for a copperlike material and depicted graphically; the most significant points arising from our analysis have also been highlighted. In the absence of thermodiffusion, the LSD and GLD have been compared with the Lord–Shulman (LS) heat transfer model and the Green–Lindsay (GL) heat transfer model.
We consider a transversely isotropic elastic medium in a twodimensional plane subjected to thermal and chemical loadings on the plane. The displacement components in the and directions are given as:
(1) 
The stress–strain temperature relations for the present problem are given as:
(2) 

(3) 

(4) 
where is the elastic coefficient, is the stress tensor, is the strain tensor, and are the tensors of thermal and diffusion moduli, respectively, and and are the thermal and diffusion relaxation times, respectively.
The equations of motion in the and directions are given by:
(5) 

(6) 
where is the density.
The heat conduction equation corresponding to the problem, introduced by some unified parameters, is defined as:

(7) 
where and are the thermal and diffusion relaxation times satisfying the relations and . For from Equation (7), we have the GLD model. However, if , we have the LSD model in the presence of thermodiffusion; is the relaxation time for the LS model.
The chemical potential is given by:
(8) 
where and are the measures of the thermodiffusion effect and the diffusive effect, respectively. The mass flux is given by:
(9) 
where is the diffusive constant. The diffusion equation is given by:
(10) 
The following nondimensional variables are as follows:
, , , , ,, , , , , , , ,
After removing the primes, the above equations can be written in a nondimensional form as:

(11) 

(12) 
(13) 
where

Therefore, the equations of motion in the and directions are given by:

(14) 

(15) 
The heat conduction equation is given by:

(16) 
where

The chemical potential is given by:

(17) 
The diffusion equation takes the form:
(18) 
where


First, we choose the following equation:
(19) 
where is the complex time constant and is the wave number in the direction. Therefore, employing the normal mode analysis, the above equations can be written as:
(20) 

(21) 

(22) 


(23) 
(24) 
When we eliminate from Equations (20)–(22) and (24), they can be simplified to the following:

(25) 
(26) 

(27) 
where













Furthermore, if we eliminate and from Equations (25)–(27), we obtain the following equation:
(28) 
where





The solution of Equation (28) is obtained as:
(29) 
where are the roots of the equation

(30) 
It can be shown that and satisfy the same equation as follows:
(31) 
Thus, the solutions are:

(32) 

(33) 
Substituting Equations (29), (32). and (33) into Equations (25) and (26), we obtain:


The mass concentration is given by:

(34) 
where
(35) 
Furthermore, the chemical potential is given by:

(36) 
where

(37) 
Therefore, by substituting Equations (29) and (32)–(34) into Equations (11)–(13), the stress components are given by:
(38) 

(39) 


(40) 
The problem is to solve subjected to the following boundary conditions. The bounding plane is subjected to a thermal loading as follows:

(41) 
The mechanical boundary conditions on the bounding plane are given by:

(42) 

(43) 
The bounding plane is subjected to a chemical loading as follows:

(44) 
Employing the normal mode analysis on the boundary conditions, we derive the following equations from Equations (41)–(44):

(45) 

(46) 

(47) 

(48) 
The aim of this section is to present the analytical numerical results obtained in the preceding sections. For the numerical computations, we have considered a copperlike material. Since is a complex time constant, we have and .
The values of the material constants are given as follows [27, 34, 35]:














































Figs. 1–7 have been plotted to study the effect of thermodiffusion on the thermophysical quantities when and for the thermoelastic diffusion model with two relaxation times (GLD) and the thermoelastic diffusion model with one relaxation time (LSD). In these figures, the continuous lines correspond to a thermodiffusive medium (WD) and the dotted lines correspond to a without thermodiffusive medium (WOD) for the GL model of heat transfer and the LS model of heat transfer.
Fig. 1 depicts the variation of the displacement component with respect to for at time . The figure shows that the displacement attains its maximum magnitude at ; as increases, the magnitude of decreases and reaches zero for both the LS and GL models. In addition, the magnitude of is larger for the GL model than the LS model. Furthermore, as seen in the figure, the magnitude of corresponding to the GLD and LSD models is greater than that of the GL and LS models.
Fig. 1. Variation of with respect to for and
Fig. 2. Variation of v with respect to y for t=0.4 and x=0.2.
Fig. 3. Variation of θ with respect to y for t=0.4 and x=0.2
Fig. 2 shows the variation of the displacement component with respect to y when and for both models. From the figure, we can observe that the vertical displacement attains its maximum magnitude on the plane ; asincreases, the magnitude of decreases. Furthermore, due the presence of thermodiffusion, the magnitude of is greater for the GLD and LSD models than for the GL and LS models. For both the diffusive medium and the elastic medium, the magnitude of for the GL model is greater than that of the LS model.
Fig. 3 depicts the variation of the temperature with respect to y for both the LS and GL models for WD and WOD when and . The figure shows that the magnitude of on the plane satisfies the thermal boundary condition of our problem as given in Equation (41). Furthermore, in the LS model, the magnitude of increases for the interval to attain its maximum value at in the diffusive medium, and then decreases sharply as increases. However, in the GL model, the smoothness in the profile of is revealed. The presence of thermodiffusion has a tendency to increase the magnitude of the profile of the temperature field.
Fig. 4 shows the variation of the stress component with respect to distance when and for the thermoelastic diffusion model with one relaxation time (LSD) and the model corresponding to two relaxation times (GLD) against the LS and GL heat transfer models. It is observed that has a value of zero on the plane for both models, which satisfies the mechanical boundary condition of the problem given in Equation (43). The figure also shows that the magnitude of is greater for the diffusive medium than for the elastic medium. The decay of the magnitude of is also faster for the elastic medium (WOD) than the diffusive medium (WD).
Fig. 5 depicts the variation of the stress component against distance for and . As shown in the figure, attains its maximum magnitude on the plane where the pressure is given and the magnitude of decreases as increases. The decay of is faster for the elastic medium (WOD) than for the diffusive medium (WOD).
Fig. 6 shows the variation of the chemical potential against the distance for and for both the LS and GL models in the presence and absence of thermodiffusion. As shown in the figure, for the LS model, the chemical potential increases for the interval and then decreases sharply as increases. However, the decrease of is slower in the GL model than in the LS model. Furthermore, because of the presence of thermodiffusion (GLD and LSD), the magnitude of becomes less in the elastic medium (WOD) than in the diffusive medium (WOD) for both the GL and LS heat transfer models.
Fig. 4. Variation of with respect to for and
Fig. 5. Variation of with respect to for and
Fig. 6. Variation of with respect to for and
Fig. 7 shows the variation of the mass concentration for and . The presence of thermodiffusion has a tendency to decrease the magnitude of the profile of for both the LS and GL models.
Figs. 8 and 9 depict the variation of the shearing stress and temperature for the LS model in a diffusive medium (WD) for different values of and when . It can be observed that as time increases, the magnitudes of the profiles of the shearing stress and the temperature distribution also increase, which supports the physical fact.
In the present analysis, the classical Fick’s diffusion law is replaced by a generalized expression that involves two relaxation times. It allows a delayed response between the relative mass flux vector and the potential gradient. A twodimensional transversely isotropic thermodiffusive medium has been considered in the context of the LSD model and the GLD model of generalized thermoelasticity. All the figures exhibit the different peculiarities that occur during the propagation of waves. The conclusions may be summarized as follows.
Fig. 7. Variation of 𝐶 with respect to 𝑦 for 𝑡=0.4 and 𝑥=0.2
Fig. 8. Profile of with respect to and when
Fig. 9. Profile of with respect to and when
We are grateful to Prof. S. C. Bose of the Department of Applied Mathematics, University of Calcutta for his valuable suggestions and guidance during the preparation of this paper.
References
[1] Nowacki W. Dynamical problems of thermoelastic diffusion in solids I. Bull. Acad. Pol. Sci. Ser. Sci. Tech. 1974; 22: 5564.
[2] Nowacki W. Dynamical problems of thermoelastic diffusion in solids II. Bull. Acad. Pol. Sci. Ser. Sci. Tech. 1974; 22: 129135.
[3] Nowacki W. Dynamical problems of thermoelastic diffusion in solids II. Bull. Acad. Pol. Sci. Ser. Sci. Tech. 1974: 22: 257266.
[4] Nowacki W. Dynamical problems of thermoelastic diffusion in elastic solids. Proceedings of vibration problems 1974; 15: 105128.
[5] Gawinecki J, Kacprzyk P, BarYoseph P. Initial boundary value problem for some coupled nonlinear parabolic system of partial differential equations appearing in thermoelastic diffusion in solid body. J. Math. Anal. Appl 2000; 19: 121130.
[6] Gawinecki J, Szymaniec A. Global solution of the Cauchy problem in nonlinear thermoelastic diffusion in solid body. Proc. Appl. Math. Mech 2002; 1: 446447.
[7] Biot MA. Thermoelasticity and irreversible thermodynamics. J. Appl. Phys 1956; 27: 240253.
[8] Lord HW, Shulman Y. A generalized dynamical theory of thermoelasticity. Journal of mechanics and physics of solids 1967; 15: 299309.
[9] Green AE, Lindsay KA. Thermoelasticity. Journal of Elasticity 1972; 2: 17.
[10] Sherief HH, Hamza FA, Saleh HA. The theory of generalized thermoelastic diffusion. International journal of engineering science 2004; 42: 591608.
[11] Sherief HH. Saleh HA. A halfspace problem in the theory of generalized thermoelastic diffusion. International journal of solids and structures 2005; 42: 44844493.
[12] Singh B. Reflection of P and SV waves from free surface of an elastic solid with generalized thermodiffusion. Journal of earth system science 2005; 114: 159168.
[13] Singh B. Reflection of SV waves from the free surface of an elastic solid in generalized thermoelastic diffusion. Journal of sound and vibration 2006; 291: 764778.
[14] Aouadi M. A generalized thermoelastic diffusion problem for an infinitely long solid cylinder. International journal of mathematics and mathematical sciences 2006; 115.
[15] Aouadi M. A problem for an infinite elastic body with a spherical cavity in the theory of generalized thermoelastic diffusion. International journal of solids and structures 2007; 44. 57115722.
[16] Aouadi M. Uniqueness and reciprocity theorems in the theory of generalized thermoelastic diffusion. Journal of thermal stresses 2007; 31: 665678.
[17] Kumar R, Gupta V. Wave propagation at the boundary surface of inviscid fluid halfspace and thermoelastic diffusion solid halfspace with Dualphaselag models. J. Solid. Mech 2015; 7(3): 312326.
[18] Othman MIA, Atwa SY, Farouk RM. The effect of diffusion on two dimensional problem of generalized thermoelasticity with Green Naghdi theory. International communications in heat and mass transfer 2009; 36: 857864.
[19] Deswal S, Choudhary S. Twodimensional interactions due to moving load in generalized thermoelastic solid with diffusion. Appl. Math. Mech. Engl. Ed 2008; 29: 207221.
[20] Deswal S, Choudhary S. Impulsive effect on an elastic solid with generalized thermodiffusion. J. Eng. Math 2009; 63: 7994.
[21] Deswal S, Choudhary S. Mechanical loads on a generalized thermoelastic medium with diffusion. Meccanica 2010; 45: 401413.
[22] Kumar R, Kansal T. Propagation of Rayleigh waves on free surface of transversely isotropic generalized thermoelastic diffusion. Applied mathematics and MechanicsEngland Edition 2008; 29(11): 14511462.
[23] Kothari S, Mukhopadhyay S. A study of influence of diffusion inside a spherical shell under thermoelastic diffusion with relaxation times. Mathematics and mechanics of solids 2012; 18(7):722737.
[24] Kothari S, Mukhopadhyay S. On the representations of solutions in the theory of generalized thermoelastic diffusion. Math. Mech. Solids 2011; 17(2): 120130.
[25] Wang HM, Ding HJ, Chen YM. Thermoelastic dynamic solution of a multilayered spherically isotropic hollow sphere for spherically symmetric problems. Acta Mechanica 2004; 173: 131145.
[26] ElSayed AM. A twodimensional generalized thermoelastic diffusion problem for a halfspace. Mathematics and mechanics of solids 2014; 52: 3743.
[27] Karmakar R, Kanoria M. Elastothermodiffusive response in a spherically isotropic hollow sphere. Journal of thermal stresses 2015; 38: 427446.
[28] Bhattacharaya D, Kanoria M. The influence of two temperature generalized thermoelastic diffusion inside a spherical shell. International journal of engineering and technical research 2014; 2: 151159.
[29] Bhattacharaya D, Kanoria M. The influence of twotemperature fractional order generalized thermoelastic diffusion inside a spherical shell. International journal of application or innovation in engineering and management 2014; 3(8): 96108.
[30] Ezzat MA, ElKaramany AS, ElBary AA. Modeling of memorydependent derivative in generalized thermoelasticity. European Physics Journal Plus 2016; 131:372.
[31] ElKaramany AS, Ezzat MA. Thermoelastic diffusion with memory dependent derivative. Journal of thermal stresses. 2016; 39(79): 10351050.
[32] Othman MIA, Eraki EEM. Effect of gravity on generalized thermoelastic diffusion due to laser pulse using dualphaselag model. Multidiscipline modeling in materials and structures 2017;14(3): 457481.
[33] Sur A, Pal P, Kanoria M. Modeling of memory dependent derivative in a fiberreinforced plate under gravitational effect. Journal of thermal stresses 2018; 41(8):973992.
[34] Sur A, Kanoria M. Propagation of thermal waves in a functionally graded thick plate. Mathematics and mechanics of solids 2015; 22(4), 718736.
[35] Sur A, Kanoria M. Fibrereinforced magnetothermoelastic rotating medium with fractional heat conduction. Procedia Engineering. 2015; 127:605612.