Ibrahim A. Abbas1
1Mathematics Department, Faculty of Science and Arts – Khulais, King Abdulaziz University, Jeddah, Saudi Arabia
1Department of Mathematics, Faculty of Science, Sohag University, Sohag, Egypt
Journal of Measurements in Engineering, Vol. 2, Issue 4, 2014, p. 175-184.
Received 20 October 2014; received in revised form 30 November 2014; accepted 12 December 2014; published 31 December 2014
Copyright © 2014 JVE International Ltd.
The theory of generalized thermoelastic diffusion with relaxation times, in the context of Lord and Shulman theory, is used to investigate the thermoelastic diffusion problem in a transversely isotropic thermoelastic infinite medium with a cylindrical cavity. The surface of cavity is taken to be traction free and subjected to heating. The analytical solution in the Laplace transform domain is obtained by using the eigenvalue approach. The numerical results of displacement, temperature, mass concentration and stress as well as the chemical potential represented analytical and graphically. An appreciable effect of relaxation times is observed on various resulting quantities.
Keywords: eigenvalue approach, Laplace transform, Lord and Shulman theory, thermoelastic diffusion.
The generalized theory of thermoelasticity is one of the modified versions of the classical uncoupled and coupled theory of thermoelasticity and has been developed in order to remove the paradox of the physically impossible phenomena of infinite velocity of thermal signals in the classical coupled thermoelasticity. The first of such modeling is the extended thermoelasticity theory (LS) of Lord and Shulman , who introduced the concept of thermal relaxation time into the classical Fourier law of heat conduction. Subsequently, modifying the stress versus strain relationship as well as the entropy relationship with relaxation time, Green and Lindsay  proposed the temperature rate dependent thermoelasticity (GL) theory. The theory was extended for anisotropic body by Dhaliwal and Sherief . Recently, [4-13] have considered different problems by using the generalized thermoelasticity theories.
Diffusion can be defined as the random walk of an ensemble of particles, from regions of high concentration to regions of lower concentration. The studying of diffusion became increasingly important. This is due mainly to its many applications in geophysics and industrial applications. In integrated circuit fabrication, diffusion is used to introduce “dopants” in controlled amounts into the semiconductor substrate. In particular, diffusion is used to form the base and emitter in bipolar transistors, form integrated resistors, and form the source/drain regions in MOS transistors and dope poly-silicon gates in MOS transistors. The theory of thermoelastic diffusion is developed by [14-16]. In this theory, the coupled thermoelastic model is used.  developed the theory of generalized thermoelastic diffusion that predicts finite speeds of propagation for thermoelastic and diffusive waves. [18, 19] studied the theory of generalized thermoelastic diffusion in an infinite elastic body with a spherical cavity while [20, 21] investigated the diffusion in a generalized thermoelastic solid in an infinite body with a cylindrical cavity. The counterparts of our problem in the contexts of the thermoelasticity theories have been considered by using analytical and numerical methods [22-29].
In this paper, the thermoelastic interaction with mass diffusion in an infinite medium with a cylindrical cavity is investigated. Laplace transform techniques and eigenvalue approach are used. The physical explanations are given seriatim corresponding to the distributions of the considered physical variables obtained in this paper. An appreciable effect of relaxation times is observed on various resulting quantities.
Following Sherief et al.  and Abbas et al. , the governing equations for an anisotropic, homogeneous thermoelastic diffusion solid in the absence of body forces, heat source and mass diffusion source considered as:
– the equations of motion:
– the equation of heat conduction:
– the equation of mass diffusion:
– the constitutive equations:
where are the elastic constants, is the density of the medium, are the components of displacement vector, is the temperature, is the concentration of diffusive material in the elastic material, and are the tensor of thermal and diffusion moduli respectively and is the time. is the thermal relaxation time, is the thermal conductivity, is the specific heat at constant strain, are the components of strain tensor, is the reference temperature assumed to obey the inequality and is the measures of thermodiffusion effect. is the diffusion relaxation time, is the diffusion coefficient and is the measures of diffusive effects. is the Kronecker symbol, are the components of stress tensor and is the chemical potential.
Let us consider an elastic unbounded body with a cylindrical cavity occupying the region whose state can be expressed in terms of the space variable and the time . are taken as the cylindrical coordinates with -axis aligned along the cylinder axis. We note that due to symmetry, the only non-vanishing displacement component is the radial one . The strain tensor has the non-vanishing following components:
Substituting for , into equations (4) and (5) we get:
where and , in which , are the coefficients of linear thermal expansion along the radial and tangential directions while and , are the coefficients of linear diffusion expansion along the radial and tangential directions. Thus, the equations (1), (2) and (3) can be written as:
Now, we introduce the following non-dimensional variables:
where and .
Using the above dimensionless variables, the equations (8)-(13), take the form where we have dropped the primes for convenience:
where , , , , , , , .
The solution for the present problem is completed by the application of the initial and boundary conditions. The initial conditions can be expressed as:
It is assumed that the inner surface of cavity is taken to be traction free and subject to a temperature which is decaying with time. The chemical potential is also assumed to be a known function of time on the inner surface of cavity. Therefore, the boundary conditions may be expressed as:
where and are constants, is an exponent of the decayed heat flux and denotes the Heaviside unit step function.
Applying the Laplace transform for above equations define by the formula:
Hence, the equations (14)-(22) take the form:
Differentiating equation (25) and (26) with respect to r and using equation (24) we get:
Equations (24), (32) and (33) can be written in a vector-matrix differential equation as follows :
where , and .
Let us now proceed to solve equation (34) by the eigenvalue approach proposed by . The characteristic equation of the matrix takes the form:
The roots of the characteristic equation (34) which are also the eigenvalues of matrix may be written in the form:
The eigenvector , corresponding to eigenvalue can be calculated as:
From equations (36)-(39), we can easily calculate the eigenvector , corresponding to eigenvalue , . For further reference, we shall use the following notations:
The solution of equation (34) which is bounded as , is given by:
where , , , is the modified of Bessel’s function of order one, , and are constants to be determined from the boundary condition of the problem. Thus, the field variables can be written for and as:
To complete the solution we have to know the constants , and by using the boundary conditions (31) we can obtain:
In order to determine the conductive and thermal temperature, displacement and stress distributions in the time domain, the Riemann-sum approximation method is used to obtain the numerical results. In this method, any function in Laplace domain can be inverted to the time domain as:
where is the real part and is imaginary number unit. For faster convergence, numerous numerical experiments have shown that the value of satisfies the relation .
The material was chosen for purposes of numerical evaluations and the constants of the problem were taken as following :
The effect of the exponent of the decayed heat flux on the distribution of all quantities through the thickness of the medium is investigated. The numerical technique outlined above was used to obtain the temperature, displacement, concentration, radial stress, hoop stress and the chemical potential distributions respect to distance at for different values of the decayed heat flux . Figure 1 show the temperature decreases with the increase of the distance. It is noticed that the temperature always starts from different certain values according to the value of and terminates at the zero value. Of course, decreases with the increase of the exponent of the decayed heat flux . It is obvious from figure 2 that the displacement is negative at where its magnitude is maximum. The displacement increases from the negative value to a positive value. In the positive values, the displacement has a peak value that depends on the values of the decayed heat flux . Figure 3 show the distribution of concentration with distance . It is noticed that the concentration have a maximum magnitude at the boundary then decreases with the increasing of distance while it increase when increasing the decayed heat flux . Figure 4 display the variation of radial stress with distance for different values of . Its magnitude increases from zero to a maximum value after that decreases rapidly as increases and then the stress reduces to zero which agree with our theoretical boundary condition. Figure 5 display the variation of the hoop stress with the distance for different values of . Figure 6 show the variation of chemical potential with the distance . From figure 6, it is seen that the chemical potential starts with its maximum value, which satisfies the boundary conditions and decreases with the increasing of the distance .
Fig. 1. The variation of temperature with distance
Fig. 2. The variation of displacement with distance
Fig. 3. The variation of concentration with distance
Fig. 4. The variation of radial stress with distance
Fig. 5. The variation of hoop stress with distance
Fig. 6. The variation of chemical potential with distance