Advances in Materials Science and Engineering

Volume 2016 (2016), Article ID 5821604, 9 pages

http://dx.doi.org/10.1155/2016/5821604

## Application of Generalized Fractional Thermoelasticity Theory with Two Relaxation Times to an Electromagnetothermoelastic Thick Plate

Department of Mathematics, Alexandria University, Alexandria, Egypt

Received 19 February 2016; Revised 18 July 2016; Accepted 25 August 2016

Academic Editor: Santiago Garcia-Granda

Copyright © 2016 A. M. Abd El-Latief. 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

The fractional mathematical model of Maxwell’s equations in an electromagnetic field and the fractional generalized thermoelastic theory associated with two relaxation times are applied to a 1D problem for a thick plate. Laplace transform is used. The solution in Laplace transform domain has been obtained using a direct method and its inversion is calculated numerically using a method based on Fourier series expansion technique. Finally, the effects of the two fractional parameters (thermo and magneto) on variable fields distributions are made. Numerical results are represented graphically.

#### 1. Introduction

Fractional calculus (FC) is a very useful tool in describing the evolution of systems with memory, which typically are dissipative and complex systems such as glasses, biopolymers, biological cells, porous materials, amorphous semiconductors, and liquid crystals. Scaling laws and self-similar behavior are supposed to be fundamental features of complex systems. In recent decades the FC and in particular the fractional differential equations have attracted interest of researchers in several areas including mathematics, physics, chemistry, biology, engineering, and economics [1–4].

FC theory has been used successfully in thermoelasticity and thermoviscoelasticity, such that a quasi-static uncoupled theory of thermoelasticity based on the fractional heat-conduction equation was put forward by Povstenko [5]. The theory of thermal stresses based on the heat-conduction equation with the Caputo time-fractional derivative is used by Povstenko [6] to investigate thermal stresses in an infinite body with a circular cylindrical hole. Sherief et al. [7] introduced a fractional order theory of thermoelasticity. Raslan [8, 9] has solved 1D problems in the context of this theory and applied this theory to 2D problem of thick plate [10]. The fractional parameter effect of this theory on thermoelastic material with variable thermal variable thermal material properties has been studied in [11, 12]. Ezzat [13] established a model of fractional heat-conduction equation by using the Taylor series expansion of time-fractional order developed by Jumarie [14].

Hamza et al. established a new mathematical model of Maxwell’s equations in an electromagnetic field in [15] and derived a fractional model for thermoelasticity associated with two relaxation times in [16]. A model for unsteady thermoelectric magnetohydrodynamics (TEMHD) flow and heat transfer of two immiscible second-grade fluids with two fractional parameters was introduced in [17].

The previous work introduced by Hamza et al. [18] describes one-dimensional problems in the context of the theory [16] in which the electromagnetic field effects are ignored. The model in the previous work depends only upon one fractional parameter which does not have electric or magnetic affects. Another application of this theory introduced in [16] depends on two fractional parameters *α* and *β*. The first fractional parameter *α* appears only in the heat equation and is absent from both the equation of motion and the constitutive equations. In this work, we solve a 1D problem for a thick plate that was not solved in the context of both theories [15, 16].

The new model here depends on two fractional parameters *α* and *β*. The first fractional parameter *α* appears in the heat equation, the equation of motion, and the constitutive equations.

The effects of the fractional parameters corresponding to two models [15, 16] are discussed. The solution is obtained in Laplace transformed domain using a direct approach. All the studied field variables are represented graphically.

#### 2. The Mathematical Model

The governing fractional Maxwell’s equations in an electromagnetic field are given by [15] where *β* is the order of the fractional derivative in Caputo sense with respect to time such that , is the magnetic permeability, and is the electric permeability. and are the magnetic and electric field intensities, respectively. is the current density and *τ* is a positive constant.

Equation (2) is called “Fractional Faraday’s Law of Magnetic Induction” which was derived in [15]. The proof uses Faraday’s induction law and the fractional Taylor’s series expansion developed by Jumarie [14].

Ohm’s law for moving media states thatwhere is the electric conductivity of the medium (assumed to be infinite) and is the displacement vector.

Since is bounded and is infinite, it follows that

The governing equations for generalized fractional thermoelasticity associated with two relaxation times in the absence of external body forces and heat sources are given by [16] the following.

(i) The constitutive equations:where is the temperature of the medium, are the components of the stress tensor, the constants and are Lamé’s constants, and , where is the coefficient of linear thermal expansion and is the order of the time-fractional derivative. is a positive constant and are the components of the strain tensor.

(ii) The strain-displacement relations:

(iii) The equation of motion:where is the density assumed independent of time , is the displacement vector component, are the components of the stress tensor, and is the component of the Lorentz force given by

(iv) The fractional heat equation:where the specific heat at constant strain is and is the thermal conductivity. is a positive constant and is a reference temperature assumed to be such that . As usual, superimposed dots denote time derivatives. The convention of summing over repeated indices is used.

#### 3. Formulation of the Problem

We consider a magnetothermoelastic thick plate of perfect conductivity occupying the region in an initial magnetic field in direction at a uniform reference temperature . This produces an induced magnetic field in the direction of the -axis and an electric field in the -axis (perpendicular to and ).

The -axis is perpendicular to the surface of the plate. The upper surface () of the plate is taken to be traction-free and is subjected to a thermal shock that is a function of time. The lower surface of the plate is taken to be thermally isolated and laid on a rigid foundation.

It is assumed that all the state functions depend on and only.

Thus

In the context of generalized fractional thermoelasticity associated with two relaxation times, the constitutive equations, the equation of motion, and fractional heat equation can be expressed in our case as follows:

Now, (1)–(4b) yield From (13), (14), and (12), we obtainFrom (13), (15), and (8), we obtainWe shall use the following nondimensional variables

Now, using the above nondimensional variables, the system of equations of the problem will reduce towhere

Also, we assume that the medium is initially at rest and the undisturbed state is maintained at uniform reference temperature. Then we have

Now, taking the Laplace transform (denoted by an overbar) with parameter of both sides of the above system of equations, we getSolving (21) and (22), we obtainThe corresponding characteristic equation of (27) is where The solution of (27) compatible with (22) has the formSubstituting (30) and (31) into (26), we get Using (31), (23)–(25) become We assume that the boundary conditions have the formwhere is a known function of . Equations (30)–(32) give

#### 4. Inversion of the Laplace Transform

We shall now outline the method used to invert the Laplace transforms in the above equations. Let be the Laplace transform of a function . The inversion formula for Laplace transforms can be written as [19]where is an arbitrary real number greater than all the real parts of the singularities of . Expanding the function in a Fourier series in the interval , we obtain the approximate formula [19]: where

Two methods are used to reduce the total error. First, the “Korrektur” method is used to reduce the discretization error. Next, *ε*-algorithm is used to reduce the truncation error and therefore to accelerate convergence.

The Korrektur method uses the following formula to evaluate the function :

We shall now describe *ε*-algorithm that is used to accelerate the convergence of the series in (37). Let be an odd natural number and let be the sequence of partial sums of (37). We define *ε*-sequence by It can be shown that [19] the sequence converges to faster than the sequence of partial sums.

#### 5. Numerical Results and Discussion

In order to obtain the solutions for the field functions in the physical domain, we have applied the Laplace inversion formula mentioned in the above section. FORTRAN programming language was used on a personal computer. The accuracy maintained was 7 digits for the numerical program. For computational purposes, a copper-like material has been taken into consideration. The values of the material constants are taken as in Table 1.