SQU Journal for Science, 2018, 23(1), 1-7 DOI: http://dx.doi.org/10.24200/squjs.vol23iss1pp1-7 Sultan Qaboos University 1 A One-dimensional Model of Thermo-electroelasticity in Extended Thermodynamics Moustafa S. Abou-Dina and Ahmed F. Ghaleb* Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt. *Email: afghaleb@sci.cu.edu.eg ABSTRACT: The subject of thermo-electroelasticity involves many complications due to the multiple ways in which the mechanical, thermal and electric fields can interact, some of these involving non-linearities. In extended thermodynamics, an additional difficulty arises due to the requirement of finiteness of the speed of propagation of the thermal disturbances. This implies, as may be observed in the extensive literature on the subject, a re-visiting of the basic postulates of thermodynamics, ultimately leading to the desired generalization. There are only a few nonlinear models dealing with this subject. In order to consider general nonlinear models, it is necessary to study linear one s first, as they represent most of the basic features of the studied phenomena. This is particularly true when the problem is tackled numerically through iteration methods, in which case the starting field equations are linear. Here we study a one-dimensional system of equations of thermo-electroelasticity in extended thermodynamics and in the quasi-electrostatic regime. The nonlinear equations are given for reference only. The mixed character, parabolic-hyperbolic, of the associated linear system is established through the study of the characteristic curves. Two speeds of wave propagation are given in evidence, one for the usual coupled thermoelastic wave, and the other for a second sound. Parabolicity is due to the quasi-static distribution of the electric field. An example concerning the half- space is treated numerically by the Cranck-Nicolson method. The curves presented clearly illustrate the propagation of two types of waves, the usual coupled thermoelastic wave, and a faster wave generated by the second sound. It is hoped that the present study will clarify the basic features of the solution, as a prelude to tackling more sample, nonlinear equations. Keywords: Electroelasticity; Extended thermodynamics; Linear theory; Characteristics; Cranck-Nicolson method. الديناميكا الحرارية الموسعة إطار الكهروحرارية في نظرية المرونةل نموذج أحادي البعد *غالب ؤادف حمدأو أبو دينه ابرص صطفىم الميكانيكي، الحراري : بين المجاالتيتضمن موضوع المرونة الحرارية الكهربائية العديد من الصعوبات بسبب الطرق المتعددة للتفاعل :صالملخ في الديناميكا الحرارية الموسعة، تنشأ صعوبة إضافية نظرا لمحدودية سرعة انتشار االضطرابات . والكهربائي، وبعضها يتضمن ظواهر غير خطية زيارة للمفاهيم األساسية للديناميكا الحرارية، واسعة النطاق حول هذا الموضوع، إعادةال وهذا يعني، كما يمكن أيضا مالحظته في المرجعيات. الحرارية ومن أجل النظر في النماذج غير من النماذج غير الخطية التي تتناول هذا الموضوع، يوجد عدد قليل فقط . مما يؤدي في نهاية المطاف إلى التعميم المطلوب ويصدق ذلك بشكل خاص عندما تعالج . السمات األساسية للظواهر المدروسةالخطية العامة، فمن الضروري دراسة النماذج الخطية أوال ألنها تمثل معظم .المشكلة عدديا من خالل الطرق التكرارية، وفي هذه الحالة تكون معادالت خطوة البداية هى المعادالت الخطية وتعَطى .ستاتيكيارية الموسعة وفي النظام شبه الكهرووهنا ندرس نظاما أحادي البُعد لمعادالت نظرية المرونة الحرارية الكهربائية في الديناميكا الحر بواسطة من خالل النظام الخطي المرتبط "ناقصي –مكافئ " النظام المختلط الناتج من المعادالت التفاضلية يتم دراسةو .الخطية كمرجع فقط غيرالمعادالت ويرجع ذلك إلى .لما يُسمى بالصوت الثاني واآلخرى المقترنةواحدة للموجة المرنة الحرارية :انتشار الموجة المنحنيات المميزة، وتوضع سرعتان في وتوضح المنحنيات . نيكلسن -، ويتم التعامل مع مثال يتعلق بنصف فراغ رقميا بواسطة طريقة كرانك للمجال الكهربائي التوزيع شبه الكهروستاتيكي ويؤمل أن . الموجة المرنة الحرارية المعتادة المقترنة، وموجة سريعة تم إنشاؤها بواسطة الصوت الثاني :شار نوعين من الموجاتالمعروضة بوضوح انت .تمهيدا لمعالجة المعادالت غير الخطية األكثر وفرة ذه الدراسة السمات األساسية للحلتوضح ه .نيكلسن -كرانك طريقة -المنحنيات المميزة –النظرية الخطية –الديناميكا الحرارية الموسعة -المرونة الكهربائية :مفتاحيةالكلمات ال M.S. ABOU-DINA and A.F. GHALEB 2 1. Introduction he theory of thermo-electroelasticity in generalized thermodynamics may be of interest in the study of electrical disturbances which accompany the propagation of heat waves at low temperature (Rybalko et al. [1], Pashitskii and Ryabchenko [2], and Pashitskii et al. [3]). Dost [4] treats the case of thermoelastic dielectrics within the nonlinear theory. Ersoy [5, 6] considers electrically and thermally conducting magnetothermoelastic solids. His theory is based on the introduction of the electric current and the heat flux vectors as independent state variables. Montanaro [7] develops a model of nonlinear thermoelasticity in extended thermodynamics for electrically polarizable and finitely deformable, heat conducting elastic continua. Ghaleb [8] presents a fully nonlinear model for electrically polarizable, heat conducting elastic continuous media in the quasi-electrostatic approximation. Kuang [9] considers wave propagation in pyroelectrics and other materials with complex structure in extended thermodynamics within the linearized theory. Montanaro [10] treats the case of electrical continuous media within Green and Naghdi thermoelasticity theory. Chandrasekharaiah [11] develops a model for piezoelectrics with the heat flux as an independent state variable. Singh [12] investigates thermo-piezoelectric solids in extended thermodynamics using Green-Lindsay and Lord-Shulman theories. Zhou and Yang [13] investigate the propagation of plane waves in pyroelectric materials in the presence of viscous effects. Solutions of concrete problems are almost inexistant in the literature. Our attention was drawn lately to a recent publication concerning two-dimensional electro-magneto thermoelastic wave propagation in an electrically conducting cylinder, within the frame of extended thermodynamics with one relaxation time [14]. All of the above references illustrate a multitude of approaches, conceived to remove the paradox of infinite propagation of thermal disturbances and to incorporate the electric interactions. Most of them rely on rigorous thermodynamics, meaning that the basic laws of thermodynamics are respected, as well as the celebrated Clausius- Duhem inequality expressing the non-negativeness of the dissipation function. However, some of these approaches still need more solid experimental verification. A common ground to all of them is the need for an enrichment of the basic thermodynamical variables describing the electro-thermomechanical system. An interesting contribution discussing the difficulties facing the different models of extended thermodynamics may be found in [15]. For dielectric materials, which are the subject of the present investigation, we use the same approach as described in [7] and [8]. The basic set of thermodynamical parameters is enriched by adjoining the heat flow vector to the classical set involving strain and the electric field. That was shown to be consistent with the basic principles of thermodynamics in earlier work by Coleman et al. [16]. A new feature of the present work is the introduction of initial and boundary conditions for the heat flow vector, independently of those for temperature. In what follows, we study a one-dimensional system of linear equations of thermo-electroelasticity in extended thermodynamics and in the quasi-electrostatic regime. Thus, any magnetic contributions are disregarded from the outset. The original nonlinear equations are given for reference only. It is hoped that the obtained results will be helpful in clarifying the basic features of the solution, thus laying the background for the numerical treatment the nonlinear equations using an iterative method. Two speeds of wave propagation are put in evidence, one for the usual coupled thermoelastic wave, and the other for second sound. A numerical example is treated for the half-space. All the unknown functions are zeroed at a sufficiently large distance from the boundary of the half-space, as a result of which reflected waves are expected to arise. In order to avoid them, the computational grid includes so mewhat restricted time values. As noted above, future work is under progress to study nonlinear wave propagation in pyroelectric materials in two dimensions. 2. The nonlinear equations The following equations are a restriction to one spatial dimension of a more general model of electro- thermoelasticity based on the introduction of the heat flow vector as an additional state variable in the free energy density of the medium and are introduced in [8]. The equations are in dimensionless form, the velocity of elastic waves being taken as unity:   ,=1 3 2 1   EEUU xx (1) ),,()(=)32(1 311 2 txfUEUUUU xxxxxxxxtt   (2) ),,( 1 =) 2 1 ( 2 0 2 21 txrAQQQbUEU xxtxx     (3)   . 1 = 0 xt QQ   (4) Here, U denotes the mechanical displacement component,  is the absolute temperature as measured from a reference temperature 0 , Q is the heat flow vector and E is the electric field component. These four relations represent respectively: (i) the equation of electrostatics expressing the vanishing of the divergence of the electric T A ONE-DIMENSIONAL, MODEL OF THERMO-ELECTROELASTICITY 3 induction, after integration and setting the integration constant equal to zero; (ii) the equation of motion; (iii) the equation of heat conduction and (iv) the Vernotte-Cattaneo law which replaces the usual Fourier law for heat conduction. The coefficient A in eq.(3) is the multiplicative factor of the squared heat flow vector in the expression for the free energy of the system. The other coefficients appearing in the equations denote material constants having obvious meaning. In particular, the coefficients  and  express the dependence of the elastic moduli on strain. Many of these coefficients will drop out in the following sections, when considerations are restricted to the linear equations. It may be noticed that the form of the Vernotte-Cattaneo law used in (4) and proposed in [8] does not yield the classical Fourier law for heat conduction in the limit of small relaxation times, and hence it cannot be used under such a restrictive condition. Considering wave propagation in a half-space with initial uniform reference temperature 0 , the characteristic length 0 L , time 0 T and heat flow vector 0 Q used to deduce the dimensionless equations are taken as follows: ,==,=,= 0 0 00 0200 CV T CL Q CV K T CV K L      where  is the material density, C is the specific heat, K is the thermal conductivity and V is a characteristic velocity of propagation of the elastic waves. Tentatively, for common solids the different coefficients take on the following values in the SU system of units: 3 4 3 2 1 1 1 1 5 1 0 = 10 , = 10 , = 10 , = 237 , = 10 .K kgm C Jkg K K W m K V ms        In what follows we consider a problem for the half-space <<0 x . The formulation necessitates a boundary condition for the new variable of state Q , to be considered side by side with the boundary condition for temperature. 3. On the character of the system of linear equations Following [8] and using a well-known formalism based on the introduction of the heat flow vector as an additional state variable, one gets the linear equations of electro-thermoelasticity in extended thermodynamics. In the one-dimensional case, after dimension analysis to reduce the speed of the classical coupled thermoelastic waves to unity, the governing linear equations may be written as a system of first-order partial differential equations in six unknown functions },,,,,{ RPEQU  of the form: ,= 3        ER x U t E (5) ),,(= 1 txf x E xx P t P             (6) ),,(= 21 txr t E t R x Q t             (7) ,= Q xt Q        (8) ,= P x U t U      (9) ,= R x U t R       (10) The two unknown functions P and R are defined from the last two equations of the above system. The parameters  , are two positive small parameters introduced artificially for convenience. Subsequently, they will be made to decrease during the computations, starting from some initial values. The function f represents the external forces of non-electric origin acting on the medium, while r denotes the volume heat supply. The other parameters appearing in the equations are dimensionless quantities involving the physical parameters of the medium. Using standard analysis, it may be shown that all the characteristics of the considered linear system of equations are real, and do not depend on the two constants  , introduced earlier. Four of these characteristics yield the speeds of wave propagation, for the usual coupled thermoelastic wave and for second sound: . 1 =1,= 3,41,2   vv These are the same as for the purely thermoelastic case, showing that in this linear approximation, the electric field does not influence the velocities of wave propagation. The remaining two velocities are equal to zero. M.S. ABOU-DINA and A.F. GHALEB 4 Thus the equations under consideration form a system of partial differential equations of the first order of mixed, parabolic-hyperbolic type. It is believed that the parabolic element is due to the considered quasi-electrostatic approximation. 4. Numerical scheme In what follows, we use the finite difference Crank-Nicholson method [17] to solve the initial set of linear equations (5)-(10). For the computational work, consider finite intervals on the x and t axes. The domain in the ),( tx plane is discretized by a grid with step length hx = and time step .= kt Let NM , be natural numbers, and the coordinates of the mesh points are: ,...,2,1,0,=,=,...,2,1,0,=,= MmmktNnnhx mn where .<0,,/=,)/(= TtbxaMTtNabh  The numerical values of the variables QPUE ,,,,  and R at the grid point ),( mn tx is denoted, respectively, by m n m n m n m n m n QPUE ,,,,  and . m n R We use the following differences approximations: 1. For first order time derivative ).(=)( 1 kO k uu u m n m nm nt    (11) 2. For first order space derivative .)() 22 ( 2 1 =)( 2 1 1 1 111 hO h uu h uu u m n m n m n m nm nx         (12) 3. For source term ).( 2 1 =)( 1 m n m n m n fff   (13) These expressions will be used to approximate the partial derivatives of all the unknown functions QPUE ,,,,  and R in the proposed system. The above replacement changes this system into a linear algebraic system for the unknowns m n m n m n m n m n QPUE ,,,,  and m n R , and then the new system can be solved easily. Also, the local truncation error of this scheme is of the order ),( 2 khO  and it is well known that the Crank-Nicholson scheme is unconditionally stable when used to solve hyperbolic partial differential equations of the first order [17]. The system of equations is solved under the following initial-boundary conditions: 0,=,0)(=,0)(=,0)(=,0)(=,0)(=,0)( xRxQxxPxUxE  ,sin)cos(1=)(0,),cos0.5(1=)(0,,cos1=)(0, 10 tptptPttUttE  ),cos(1=)(0,)),(0,(=)(0,,cos1=)(0, 0 trtRtBitQtt   0,=),(=),(=),(=),(=),(=),( txRtxQtxtxPtxUtxE finalfinalfinalfinalfinalfinal  and for the following values of the material and geometrical constants .=0,=0.0012,=,= 43231 QfffERf   The unknowns P and R are not independent, but defined through the other unknowns. Thus, the constants 10 , pp and 0r cannot be arbitrary. Simple calculations show that   . 1 =, 2 1 =, 1 = 3 013 0 0          rp r p Taking 0.5,=0.1,=35,==,10=350,=0,= 1  MNTba 0.1.=0.98,=0.1,=0.0001,==0.5,=0.5,=0.1,= 321 Bi  All variables were zeroed at finalx , on the basis that waves have not reached this point at time finalt and therefore no reflected waves will appear during this period of time. The maximal value for time is 35finalt , while the values of the A ONE-DIMENSIONAL, MODEL OF THERMO-ELECTROELASTICITY 5 spatial coordinate x were allowed to run up to the value = 10 350final finalx t , in order to be able to view the fast wave travelling with velocity 3.16 1  . The solutions for the mechanical displacement, temperature, heat flow vector and electric field are represented as 3-D surfaces on Figures 1-4, in which one clearly notices two separate propagating waves generated by the boundary regimes. The computations were repeated many times for different values of  and  , and it was noted that all the figures remained unchanged as these two constants were made smaller by many orders of magnitude, thus justifying their introduction for convenience. Figure 1. Mechanical displacement component U. Figure 2. Temperature  . M.S. ABOU-DINA and A.F. GHALEB 6 Figure 3. Heat flow vector component Q. Figure 4. Electric field component E. 5. Conclusion A model of linear thermo-electroelasticity in extended thermodynamics and in the quasi-electrostatic approximation has been investigated using the method of characteristics, and then numerically for given initial and boundary conditions. The main ingredient in this model is a Cattaneo-type evolution equation, which effectively requires the heat flow vector to be considered as an additional state function, independent of temperature. This requires the introduction of initial and boundary conditions for the heat flow vector, independently of those for temperature, a fact that represents a fundamental difference from the thermoelasticity problems solved within the frame of different theories of extended thermodynamics. For these, the additions of temporal derivatives occur in the heat equation itself. Moreover, the problems are usually solved by Laplace transform, with the elimination of the heat flow vector in favor of the other unknowns of the problem, hence not requiring additional limiting conditions for this variable. The presently used boundary condition for the heat flow vector is a Robin thermal condition which may be controlled A ONE-DIMENSIONAL, MODEL OF THERMO-ELECTROELASTICITY 7 experimentally. Other types of boundary conditions for the heat flow vector may equally well be used for the computations. Dimension analysis and a study of the characteristic curves has revealed the mixed “parabolic-hyperbolic” character of the system of linear equations of electro-thermoelasticity. The parabolic element is due to the used quasi- electrostatic approximation. Moreover, within the linear theory, the electric field does not influence the speeds of propagation of the waves. This will certainly not be true anymore in the nonlinear case. A numerical application clearly shows two types of waves: The usual, coupled thermoelastic wave and the second sound, propagating in the medium. It is hoped that the presented results, together with the future extension to include nonlinear interactions, will help investigating interesting problems in thermo-electroelasticity. References 1. Rybalko, A.S., Rudavskii, E., Rubets, S., Tikhiy, V., Derkach, V. and Tarapov, S. Electric induction in He II. Journal of Low Temperature Physics, 2007, 148, 527-534. 2. Pashitskii, E.A., Ryabchenko, S.M. On the cause of electric activity of superfluid Helium upon excitation of a second sound wave and normal-component velocity oscillations in it. Journal of Low Temperature Physics, 2007, 33, 12-21. 3. Pashitskii, E.A., Tkachenko, O.M., Grygoryshyn, K.V. and Lev, B.I. On the nature of electrical activity in superfluid Helium at second sound excitations. Ukrainian Journal of Physics, 2009, 54, 89-93. 4. Dost, S. On generalized thermoelastic dielectrics. Journal of Thermal Stresses, 1981, 4(1), 51-57. 5. Ersoy, Y. A new nonlinear constitutive theory for conducting magnetothermoelastic solid. International Journal of Engineering Science. 1984, 22, 683-705. 6. Ersoy, Y. A new nonlinear constitutive theory of electric and heat conductions for magnetoelastothermoelectrical anisotropic solids. International Journal of Engineering Science, 1986, 24, 867-882. 7. Montanaro, A. On the constitutive relations for second sound in thermo-electroelasticity. Archives of Mechanics, 2011, 63(3), 225-254. 8. Ghaleb, A.F. Coupled thermoelectroelasticity in extended thermodynamics. Encyclopedia of Thermal Stresses (C), Ed. R.B. Hetnarski, Springer, 2014, 767-774. 9. Kuang, Z.B. Theory of Electroelasticity, Springer, 2014. 10. Montanaro, A.A Green-Naghdi approach for thermo-electroelasticity. Journal of Physics: Conference Series, 2015, 633, 012129. Doi:10.1088/1742-6596/633/1/012129. 11. Chandrasekharaiah, D.S. A generalized linear thermoelasticity theory for piezoelectric media. Acta Mechanica, 1988, 71, 39-49. 12. Singh, B. On the theory of generalized thermoelasticity for piezoelectric materials. Applied Mathematics and Computation, 2005, 171(1), 398-405. 13. Zhou, Z.D. and Yang, F.P. Plane waves in pyroelectrics with viscous effect. Acta Mechanica, 2014, 225, 509-521. Doi 10.1007/s00707-013-0962-7. 14. Sherief, H.H. and Allam, A.A. Electro-magneto interaction in a two-dimensional generalized thermoelastic solid cylinder. Acta Mechanica, 2017, 228, 2041-2062. Doi 10.1007/s00707-017-1814-7. 15. Li, S. and Cao, B.Y. On thermodynamics problems in the single-phase-lagging heat conduction model. Entropy, 2016, 18, 391-399. Doi:10.3390/e18110391. 16. Coleman, B.D., Fabrizio, M. and Owen, D.R. On the thermodynamics of second sound in dielectric crystals. Archive for Rational Mechanics Analysis, 1982, 80(2), 135-158. 17. Smith, G.D. Numerical Solution of Partial Differential Equations. (Second ed.), Oxford University Press, 1978. Received 4 April 2017 Accepted 13 August 2017