`Mathematical Problems in EngineeringVolume 2011, Article ID 731876, 23 pageshttp://dx.doi.org/10.1155/2011/731876`
Research Article

## One-Dimensional Problem of a Conducting Viscous Fluid with One Relaxation Time

Department of Mathematics, Faculty of Girls, Ain Shams University, Cairo 11757, Egypt

Received 2 August 2010; Revised 28 December 2010; Accepted 23 February 2011

Academic Editor: Sergio Preidikman

Copyright © 2011 Angail A. Samaan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We introduce a magnetohydrodynamic model of boundary-layer equations for conducting viscous fluids. This model is applied to study the effects of free convection currents with thermal relaxation time on the flow of a viscous conducting fluid. The method of the matrix exponential formulation for these equations is introduced. The resulting formulation together with the Laplace transform technique is applied to a variety problems. The effects of a plane distribution of heat sources on the whole and semispace are studied. Numerical results are given and illustrated graphically for the problem.

#### 1. Introduction

The modification of the heat-conduction equation from diffusive to a wave type may be affected either by a microscopic consideration of the phenomenon of heat transport or in a phenomenological way by modifying the classical Fourier law of heat conduction.

Many authors have considered various aspects of this problem and obtained similarity solutions. Samaan [1] investigated steady oscillating magnetohydrodynamic flow in a circular pipe. Analytical and numerical methods for the momentum and energy equations of a viscous incompressible fluid along a vertical plate have been considered by Samaan [2]. Chamkha [3] studied the magnetohydrodynamic flow of a uniformly stretched vertical permeable surface in the presence of heat generation/absorption and chemical reaction. Ishak et al. [4] investigated theoretically the unsteady mixed convection boundary layer flow and heat transfer due to a stretching vertical surface in a quiescent viscous and incompressible fluid.

Many authors presented some mathematical results, and a good amount of references can be found in the papers by Liao and Pop [5] and Nazar et al. [6]. Further, the stagnation region encounters the highest pressure, the highest heat transfer, and the highest rate of mass deposition studied by Wang [7]. Singh et al. [8] investigated the problem of heat transfer in the flow of an incompressible fluid.

Samaan [9] investigated the heat and mass transfer over an accelerating surface with heat source in presence of suction and magnetic field. The flow of an unsteady, incompressible magnetohydrodynamics (MHDs) viscous fluid with suction is investigated by Muhammad et al. [10]. Heat and mass transfer over an accelerating surface with heat source in presence of magnetic field is derived by Samaan [11].

Recently, Samaan [12] studied the effects of variable viscosity and thermal diffusivity on the steady flow in the presence of the magnetic field. variable viscosity effects on hydrodynamic boundary layer flow along a continuously moving vertical plate were done by Mostafa [13]. Concerning the studies of state space formulation for MHDs and free convection flow with two relaxation times and the free convection effective perfectly conducting couple stress fluid, we may refer to Samaan [14], Ezzat et al. [15], and Hayat et al. [16]. Using differential transform method and Pade approximate for solving MHDs flow in a laminar liquid film from a horizontal stretching surfaces investigated by Rashidi et al. [17].

In the present work, we use a more general model of magnetohydrodynamic free convection flow, which also includes the relaxation time of heat conduction and the electric displacement current [18, 19]. An attempt to account for the time dependence of heat transfer, Cattaneo [20] and Vernotte [21] independently modified Fourier’s law to include the relaxation time of the system. Generalized thermoelasticity stands for a hyperbolic thermoelasticity in which a thermomechanical load applied to a body is transmitted in a wave-like manner throughout the body, only transient thermoelastic waves are included in the survey by Hertnarski and Ignaczak [22]. The inclusion of the relaxation time and the electric displacement current modifies the governing thermal and electromagnetic field equations, changing them from the parabolic to a hyperbolic type, and thereby eliminating the unrealistic result that thermal and electromagnetic disturbances are realized instantaneously within a fluid.

The solution is obtained using a state-space approach. The importance of state-space analysis is recognized in fields where the time behavior of physical process is of interest. The first writer to introduce the state space approach in magnetohydrodynamic free convection flow was Ezzat [23, 24]. His works dealt with free convection flow in the absence of the applied magnetic field or when there are no heat sources. The present work is an attempt to generalize these results to include the effects of heat sources. The results obtained are used to solve a problem for the whole space with a plane distribution of heat sources. The solutions obtained are utilized in combination with the method of images to obtain the solution for a problem with heat sources distributed situated inside a semispace whose surface bounded by an infinite vertical plate.

The Laplace transform techniques is applied to one-dimensional problem. The inversion of the Laplace transforms is carried out using a numerical technique [25].

#### 2. Formulation of the Problem

Let a constant magnetic field of strength act in the direction of the -axis. Due to the effect of this magnetic field, there arises in the medium an induced magnetic field and an induced electric field . All the considered functions will depend on and the time only.

The electromagnetic quantities satisfy Maxwell’s equations [26] where is the electric current density, and are the magnetic and electric permeabilities, respectively, and , are the magnetic and electric induction vectors, respectively, and dotte is the time variable.

These equations are supplemented by Ohm’s law where is velocity vector of the fluid and is the electric conductivity.

As mentioned above, the applied field and the induced magnetic field has the components The vector and will have nonvanishing components only in the -direction. That is, where from (2.5) The vector (2.1) and (2.2) reduce to the following scalar equations Eliminating between (2.8) and (2.9), we obtain The pondermotive force has one nonvanishing component in the -direction where is the electric conductivity.

We investigate the free convective heat transfer in an incompressible hydromagnetic flow past an infinite vertical plate. The -axis is taken along the plate in the direction of the flow and the -axis normal to it. Let be the component of the velocity in the direction. All the fluid properties are assumed constant, except that the influence of the density variation with temperature is considered only in the body force term. In the energy equations, terms representing viscous and Joule’s dissipation are neglected, as they are assumed to be very small in free convection flow [27]. Also, in the energy equation, the term representing the volumetric heat source is taken as a function of the space variables. In view of the assumptions, the equations that govern unsteady one-dimensional free convection flow in an incompressible conducting fluid through a porous medium bounded by an infinite nonmagnetic vertical plate in the presence of a constant magnetic field are (2.10)–(2.12), and the equations describing the flow in the boundary layer reduce to [2830] And the constitutive equation In these equations, is the permeability of the porous medium, is the Alfven velocity, is the acceleration due to gravity, is the coefficient of volume expansion, is the temperature distribution, is the temperature of the fluid away from the plate, is specific heat at constant pressure, is the thermal diffusivity, is the relaxation time, is the density of the fluid far from the surface, and is the density of the fluid.

Let us introduce the following nondimensional variables In view of these transformations, (2.13)–(2.14) become where is the Grashof number, is the speed of light given by , , is a measure of the magnetic viscosity, and is the Prandtl number. From now on, we will consider a heat source of the form where and are the Dirac delta function and the Heaviside unit step function, respectively, is the strength of the applied heat source, and is a constant.

We will assume that the initial state of the medium is quiescent. Taking the Laplace transform, defined by the relation of both sides of (2.16), we obtain Eliminating between (2.19)-(2.20), we get where

We will choose as state variables the temperature , the velocity component , and induced magnetic field and their gradients. Equations (2.22)–(2.23) can be written as where . The above equations can be written in matrix form as where In order to solve the system (2.26), we need first to find the form of the matrix .

The characteristic equation of the matrix has the form where The roots , , and of (2.28) satisfy the relations One of the roots, say , has a simple expression given by The other two roots satisfy the relation The Taylor series expansion of the matrix exponential has the form

Using the well-known Cayley-Hamilton theorem, the infinite series can be truncated to the following form: where is the unit matrix of order 6 and are some parameters depending on and . The characteristic roots , of the matrix must satisfy the equations The solution of this system of linear equations is given by where Substituting for the parameters from (2.36) into (2.34) and computing , , , and , we get the elements () of the matrix to be

It is worth mentioning here that (2.32) have been used repeatedly in order to write the above entries in the simplest possible form. We will stress here that the above expression for the matrix exponential is a formal one. In the actual physical problem, the space is divided into two regions accordingly as or . Inside the region , the positive exponential terms, not bounded at infinity, must be suppressed. Thus, for , we should replace each by and each by . In the region , the negative exponentials are suppressed instead.

We will now proceed to obtain the solution of the problem for the region . The solution for the other region is obtained by replacing each by .

The formal solution of system (2.26) can be written in the form Evaluating the integral in (2.39) using the integral properties of the Dirac delta function, we obtain where Equation (2.40) expresses the solution of the problem in the Laplace transform domain for in terms of the vector representing the applied heat source and the vector representing the conditions at the plate . To evaluate the components of this vector, we note first, due to the symmetry of the problem, that the velocity component and induced magnetic field component vanish at the plane source of heat, thus Gauss’s divergence theorem will now be used to obtain the thermal condition at the plane source. We consider a short cylinder of unit base, whose axis is perpendicular to the plane source of heat and whose bases lie on opposite sides of it. Taking the limit as the height of the cylinder tends to zero and noting that there is no heat flux through the lateral surface, we get We will use the generalized Fourier’s law of heat conduction in the nondimensional form [31], namely, Taking the Laplace transform of both sides of this equation and using (2.43), we obtain Equation (2.42) and (2.45) give three components of the . To obtain the remaining three components, we substitute on both sides of (2.40) getting a system of linear equations whose solution gives Inserting the values from (2.42) and (2.46) into the right-hand side of (2.40) and performing the necessary matrix operations, we obtain where Also, substituting from (2.49) into (2.21) the induced electric field is given by Equations (2.47)–(2.51) determine completely the state of the fluid for . We mention in passing that these equations give also the solution to a semispace problem with a plane source of heat on its boundary that constitutes a rigid base. As mentioned before, the solution for the whole space when is obtained from (2.47)–(2.51), by taking the symmetries under considerations.

We will show that the solution obtained above can be used as a set of building blocks from which the solutions to many interesting problems can be constructed. For future reference, we will write down the solution to the problem in the case when the source of heat is located in the plane , instead of the plane . In this case, we have where the upper (plus) sign denotes the solution in the region , while the lower (minus) sign denotes the solution in the region .

#### 3. Applications

We will now consider the problems of a semispace with a plane source of heat located inside the medium at the position and subject to the following boundary conditions. (i)The shearing stress and the induced magnetic field are vanishing at the wall (), (ii)The temperature is kept at a constant value , which means that the temperature increment satisfies This problem can be solved in a manner analogous to the outlined above though the calculations become quite messy. We will instead use the reflection method together with the solution obtained above for the whole space. This approach was proposed by Nowacki in the context of coupled thermoelasticity [31].

The boundary conditions of the problem can be satisfied by locating two heat sources in an infinite space, one positive at and the other negative of the same intensity at . The temperature increment is obtained as a superposition of the temperature for both plane distribution. Thus, , where is the temperature due to the positive heat source, given by (2.52) and is the temperature due to the negative heat source and is obtained from (2.52) by replacing with , and noting that for all points of the semispace, we have . Thus, is given by Combining (2.52) and (3.2), we obtain Clearly, this distribution satisfies the boundary condition (3.2). We turn now to the problem of finding the distributions velocity, the induced magnetic field, and the induced electric field. Unfortunately, the above procedure of superposition cannot be applied to these fields as in the temperature fields. We define the scalar stream function by the relation By integration (2.53) and using (3.5), we obtain the stream function due to the positive heat source at the position as where the upper sign is valid for the region and the lower sign is valid for the region . Similarly, the stream function for the negative heat source at is given by Since is a scalar field, we can use superposition to obtain the stream function for the semispace problem as Using (3.8) and (3.5), we obtain the velocity distribution Differentiating (3.9) and using the resulting expressions together with (2.28), we obtain Substituting (3.10) and (2.21), we can get Clearly, in agreement with (3.1).

#### 4. Inversion of the Laplace Transform

In order to invert the Laplace transforms in the above equations, we will use a numerical technique based on Fourier expansions of functions.

Let be the Laplace transform of a given function . The inversion formula of Laplace transforms states that where is an arbitrary positive constant greater than all the real parts of the singularities of . Taking , we get This integral can be approximated by Taking , we obtain For numerical purposes, this is approximated by the function where is a sufficiently large integer chosen such that where is a preselected small positive number that corresponds to the degree of accuracy to be achieved, Formula (3.11) is the numerical inversion formula valid for [22]. In particular, we choose , getting

#### 5. Numerical Results

The constants of the problem were taken as , , , , and . All constants are given in SI units. The computations were carried out for the two different values of time, namely, and 1. The functions , , , and are evaluated. The results are shown in Figures 1, 2, 3, and 4. In these figures, solid lines represent the solution corresponding to using the generalized Fourier equation of heat conduction (), while dashed lines represent the solution corresponding to using the classical Fourier heat equation ().

Figure 1: Variation of with , (−) for , (− ·) for , and at , , , , and representing onset stationary convection.
Figure 2: Variation of with , (−) for , (− ·) for , and at , , , , and representing onset stationary convection.
Figure 3: Variation of with , (−) for , (− ·) for , and at , , , , and representing onset stationary convection.
Figure 4: Variation of with , (−) for , (− ·) for , and at , , , , and representing onset stationary convection.

The important phenomenon observed in all computations is that the solution of any of the considered functions vanishes identically outside a bounded region of space surrounding the heat source at a distance from it equal to , and say that is a particular value of depending only on the choice of and is the location of the wave front. This demonstrates clearly the difference between the solution corresponding to using classical Fourier heat equation () and to using the non-Fourier case (). In the first and older theory, the waves propagate with infinite speeds, so the value of any of the functions is not identically zero (though it may be very small) for any large value of . In the non-Fourier theory, the response to the thermal and mechanical effects does not reach infinity instantaneously but remains in a bounded region of space given by for the semispace problem and by Min for the whole space problem.

We notice that results for all functions considered in the semispace problem when the relaxation time is appeared in heat equation are distinctly different from those when the relaxation time disappeared.

We also notice that for small values of time, the solution is localized in a finite region near the plane of heat sources. This region grows with increasing time until it fills the whole boundary-layer region.

Also, we observe from Figure 5. that the effect of heating by free convention currents when . In this case, it is noticed that as Alfven velocity increases, the velocity is found to decrease. This is mainly due to the fact that the effect of the magnetic field corresponds to a term signifying a positive force that tends to decelerate the fluid particles. Also, it is noticed that the velocity increases as increases, while it decreases when the Prandtl number increases.

Figure 5: Velocity distribution with , for different values of , , and , with , , , and representing onset stationary convection.

We notice that results for the temperature distribution when the relaxation time is appears in the heat equation are distinctly different from those when the relaxation time is not mentioned in the heat equation. This is due to the fact that thermal waves in the Fourier theory of heat equation travel with an infinite speed of propagation as opposed to finite speed in the non-Fourier case. We also notice that for small values of time, the solution is localized in a finite region near the plane surface. This region grows with increasing time until it fills the whole boundary-layer region. As time increases, results for both theories of heat equation almost coincide which is expected, since magnetohydrodynamic free convection flow effects are short lived. At all values of time, the velocity distributions for both theories coincide.

The effect of the relaxation time on the variation temperature and velocity distribution for this problem is shown in Figures 1 and 2, respectively.

The values of Grashof number have been chosen as they are interesting from physical point of view.

#### 6. Concluding Remarks

Many metallic materials are manufactured after they have been redefined sufficiently in the molten state. Therefore, it is a central problem in metallurgical chemistry to study the free convection effects on conducting fluid metal. For instance, liquid sodium Na (100°C) exhibit very small electrical resistivity.

The effects of Grashof number, Alfven velocity, Prandtl number, and relaxation time on the oneset temperature and velocity distribution are discussed. A discussion is provided for the effects of heating on a viscous conducting fluid is given.

The importance of state space analysis is recognized in fields where the time behavior of any physical process is of interest.

The state-space approach is more general than the classical Laplace and Fourier transform techniques. Consequently, state space is applicable to all systems that can be analyzed by integral transforms in time and is applicable to many systems for which transform theory breaks down [32].

Owing to the complicated nature of the governing equations for the unsteady magnetohydrodynamic with gradient pressure flow, few attempts have been made to solve problems in this field. These attempts utilized approximate methods valid for only a specific range of some parameters.

In this work, the method of direct integration by means of the matrix exponential, which is a standard approach in modern control theory and developed in detail in many texts such as Ogata [33], and Ezzat et al. [34, 35], is introduced in the field of magnetohydrodynamics and applied to two specific problems in which the temperature and velocity are coupled. This method gives exact solutions in the Laplace transform domain without any assumed restrictions on either the applied magnetic field or the velocity, temperature distributions, and viscoelastic parameter.

The method used in the present work is applicable to a wide range of problems. It can be applied to problems which are described by the linearized Navier-Stokes equations. The same approach was used quite successfully in dealing with problems in thermoelasticity theory [14, 15, 19, 34, 35].

#### Nomenclature

 𝑡: Time 𝑥, 𝑦: Coordinates system 𝑣: =(𝑢,0,0) velocity vector of the fluid 𝑇: Temperature distribution 𝑇𝑜: Temperature of the plate 𝑇∞: Temperature of the fluid away from the plate 𝑝𝑟: Prandtl number 𝐺𝑟: Grashof number 𝑐𝑝: Specific heat at constant pressure 𝑔: Acceleration due to gravity 𝜌: Density 𝜆: Thermal diffusivity ℎ: Induced magnetic field 𝒬𝑜: Constant 𝒬: The strength of the applied heat source 𝐸: Electric field 𝑐: Spead of light 𝐹: Magnetic viscosity 𝐽: Electric current density 𝜇𝑜: Magnetic permeability 𝛽: Coefficient of volume expansion 𝜀𝑜: Electric permeability 𝐻𝑜: Magnetic intensity vector 𝜐0: Relaxation time 𝛿(𝑦): Dirac delta function 𝐻(𝑡): Heaviside unit step function 𝐾: Permeability of the porous medium 𝛼: Alfven velocity 𝜎𝑜: Electric conductivity.

#### Acknowledgment

The author is very grateful to the referees for their helpful suggestions which improved the paper.

#### References

1. A. A. Samaan, “Oscillating in a circular pipe with adverse pressure gradient,” International Journal of Dynamics of Fluids, vol. 5, pp. 119–127, 2009.
2. A. A. Samaan, “Analytical and numerical methods for the momentum, species concentration and energy equations of a viscous incompressible fluid along a vertical plate,” Nuovo Cimento della Società Italiana di Fisica B, vol. 125, no. 4, pp. 439–455, 2010.
3. A. J. Chamkha, “MHD flow of a uniformly streched vertical permeable surface in the presence of heat generation/absorption and a chemical reaction,” International Communications in Heat and Mass Transfer, vol. 30, no. 3, pp. 413–422, 2003.
4. A. Ishak, R. Nazar, and I. Pop, “Unsteady mixed convection boundary layer flow due to a stretching vertical surface,” The Arabian Journal for Science and Engineering B, vol. 31, no. 2, pp. 165–182, 2006.
5. S. J. Liao and I. Pop, “Explicit analytic solution for similarity boundary layer equations,” International Journal of Heat and Mass Transfer, vol. 47, no. 1, pp. 75–85, 2004.
6. R. Nazar, N. Amin, and I. Pop, “Unsteady boundary layer flow due to a stretching surface in a rotating fluid,” Mechanics Research Communications, vol. 31, no. 1, pp. 121–128, 2004.
7. C. Y. Wang, “Stagnation flow towards a shrinking sheet,” International Journal of Non-Linear Mechanics, vol. 43, no. 5, pp. 377–382, 2008.
8. K. R. Singh, V. K. Agrawal, and A. Singh, “Heat transfer in the flow of a non-Newtonian second-order fluid between two enclosed counter torsionally oscillating DISC,” International Journal of Dynamics of Fluids, vol. 4, pp. 145–157, 2008.
9. A. A. Samaan, “Heat and mass transfer over an accelerating surface with heat source in presence of Suction,” International Journal of Theoretical and Applied Mechanics, vol. 3, pp. 195–206, 2008.
10. R. Muhammad, K. Fakhar, and A. Farooq, “Unsteady magneto-fluid-dynamics fluid and heat flow,” International Journal of Dynamics of Fluids, vol. 4, pp. 1–12, 2008.
11. A. A. Samaan, “Heat and mass transfer over an accelerating surface with heat source in presence of magnetic field,” International Journal of Theoretical and Applied Mechanics, vol. 4, pp. 279–291, 2009.
12. A. A. Samaan, “The effects of variable viscosity and thermal diffusivity on the steady flow in the presence of the magnetic field,” Nuovo Cimento della Societa Italiana di Fisica B, vol. 124, no. 9, pp. 931–946, 2009.
13. A. A. Mostafa, “A note on variable viscosity and chemical reaction effect on mixed convection heat and mass transfer along a seme-infinite vertical plate,” Mathematical Problems in Engineering, vol. 10, pp. 1–7, 2007.
14. A. A. Samaan, “State space formulation for magnetohydrodynamic free convection flow with two relaxation times,” Applied Mathematics and Computation, vol. 152, no. 2, pp. 299–321, 2004.
15. M. A. Ezzat, M. Zakaria, A. A. Samaan, and A. Abd Elbary, “Free convection effective perfectly conducting couple stress fluid,” Journal of Technical Physics, vol. 47, pp. 5–30, 2006.
16. T. Hayat, H. Mambili-Mamboundou, and F. M. Mahomed, “Unsteady solutions in a third-grade fluid filling the porous space,” Mathematical Problems in Engineering, vol. 2008, Article ID 139560, 13 pages, 2008.
17. M. M. Rashidi and M. Keimanesh, “Using differential transform method and padé approximant for solving mhd flow in a laminar liquid film from a horizontal stretching surface,” Mathematical Problems in Engineering, vol. 2010, Article ID 491319, 2010.
18. M. A. Ezzat and M. Z. Abd-Elaal, “Free convection effects on a viscoelastic boundary layer flow with one relaxation time through a porous medium,” Journal of the Franklin Institute B, vol. 334, no. 4, pp. 685–706, 1997.
19. M. A. Ezzat, A. S. El-Karamany, and A. A. Samaan, “The dependence of the modulus of elasticity on reference temperature in generalized thermoelasticity with thermal relaxation,” Applied Mathematics and Computation, vol. 147, no. 1, pp. 169–189, 2004.
20. C. Cattaneo, “On the conduction of heat,” Comptes Rendus de l'Académie des Sciences, vol. 4, pp. 247–431, 1958.
21. P. Vernotte, “Les paradoxes de la théorie continue de l'équation de la chaleur,” Comptes Rendus de l'Académie des Sciences, vol. 246, pp. 3154–3155, 1958.
22. R. B. Hetnarski and J. Ignaczak, “Generalized thermoelasticity,” Journal of Thermal Stresses, vol. 22, no. 4-5, pp. 451–476, 1999.
23. M. A. Ezzat, “State space approach to unsteady two-dimensional free convection flow through a porous medium,” Canadian Journal of Physics, vol. 72, no. 5-6, pp. 311–317, 1994.
24. M. A. Ezzat, “Fundamental solution in thermoelasticity with two relaxation times for cylindrical regions,” International Journal of Engineering Science, vol. 33, no. 14, pp. 2011–2020, 1995.
25. G. Honig and U. Hirdes, “A method for the numerical inversion of Laplace transforms,” Journal of Computational and Applied Mathematics, vol. 10, no. 1, pp. 113–132, 1984.
26. A. H. Nayfeh and S. Nemat-Nasser, “Electromagneto- thermoelastic plane waves in solids,” Journal of Applied Mechanics E, vol. 39, pp. 108–113, 1972.
27. J. P. Holman, Heat Transfer, Mc Graw-Hill/Kogahusha, Tokyo, Japan, 1976.
28. M. A. Ezzat, “Free convection effects on perfectly conducting fluid,” International Journal of Engineering Science, vol. 39, no. 7, pp. 799–819, 2001.
29. D. W. Beard and K. Walters, “Elastico-viscous boundary-layer flows. I. Two-dimensional flow near a stagnation point,” Proceedings of the Cambridge Philosophical Society, vol. 60, pp. 667–671, 1964.
30. J. P. Goedbloed and S. Poedts, Principles of Magnetohydrodynamics, Cambridge University Press, Cambridge, UK, 2004.
31. W. Nowacki, “Some dynamic problems of thermoelasticity,” Archiwum Mechaniki Stosowanej, vol. 11, pp. 259–283, 1959.
32. D. Wibery, Theory and Problems of State Space and Linear System, Schaum’s Outline Series in Engineering, McGraw-Hill, New York, NY, USA, 1971.
33. K. Ogata, State Space Analysis Control System, chapter 6, Prentice-Hall, Englewood Cliffs, NJ, USA, 1967.
34. M. A. Ezzat, A. S. El-Karamany, and A. A. Samaan, “State-space formulation to generalized thermoviscoelasticity with thermal relaxation,” Journal of Thermal Stresses, vol. 24, no. 9, pp. 823–846, 2001.
35. M. A. Ezzat, M. I. Othman, and A. A. Smaan, “State space approach to two-dimensional electromagneto-thermoelastic problem with two relaxation times,” International Journal of Engineering Science, vol. 39, no. 12, pp. 1383–1404, 2001.