Eigen value approach with memory dependant derivative on homogeneous isotropic infinitely extended rotating plate of a finite thickness in absence of heat source

. Our present manuscript is an attempt to derive a model of generalized thermoelasticity with dual phase lag heat conduction by using the methodology of memory dependent derivative for a isotropic rotating plate subject to the prescribed boundary conditions with constant magnetic and electric intensities. Two integral transform such as Laplace transform for time variable and Fourier transform for space variable are employed to the governing equations to formulate vector-matrix differential equation which is then solved by eigenvalue approach methodology. The inversion of two integral transformations is carried out using suitable numerical techniques. Numerical computations for displacement, thermal strain and stress component


Introduction
At present, the theory of magneto thermo elasticity which deals with the interactions among the strain, temperature and magnetic field has attracted research personalities not only due to its extensive uses in different fields such as geophysics, for understanding the effects of the earth's magnetic field, seismic waves, emission of electromagnetic radiations from nuclear devices etc., but also for its applications in different fields such as engineering, high energy The classical coupled thermo elasticity theory of Biot [1] includes the strain-rate gradient in particle accelerators, nuclear reactor's design etc.
Fourier's law of heat conduction equation which is transformed to parabolic type diffusion equation predicting infinite speed of propagation. Whereas Cattaneo [4] and Vernotte [18] theory admits the existence of thermal waves which propagate at a finite speed. Tzou [17] made the generalization of thermo elasticity theory with dual phase lag and suggested two different phase lags (DPL) in the Fourier Law of heat conduction in which the first is for the heat flux vector and the second for the temperature gradient. Generalization of thermo elasticity theory with two relaxation time has been made by Green and Lindsay [14] who obtained the explicit versions of the constitutive equations. Ezzat [10] obtained independently these equations and the fundamental solutions for this theory. Ezzat and El-Karamany [12] established the uniqueness and reciprocity theorems for anisotropic media. Ezzat, Othman and El-Karamany [13] introduced a formulation of the boundary integral equation method for generalized thermo elasticity with two relaxation times. The heat conduction equation and all the equations of motion were modified by Green and Lindsay [14] without violating Fourier Law introducing two relaxation time parameters. Ibrahim A. Abbas et al. [20] [23] discussed fractional order thermoelasticity in infinite medium.
The theory of materials with voids or vacuous pores was first proposed by Jace W.Nunziato and Stephen C. Cowin [21] in 1979 and in a study [22] by Marin Marin, Mohamed I. A. Othman and Ibrahim A. Abbas in which the domain of influence theorem has been extended to cover the generalized thermoelasticity of anisotropic bodies with voids in the context of Lord-Shulman and Green-Lindsay theories where it was proved that for a finite time > 0 the displacement field, the change in volume fraction and the temperature generate no disturbance outside a bounded domain .
The G-N theory introduced by Green and Naghdi [15] is based on generalized thermo elasticity where thermal displacement gradient in the constitutive equations is considered.
So we are getting different types of problems which have been solved using different approaches: Solutions of some problems are not convergent while involving potential functions, whereas they are actually convergent in physical problems.
Because physical problems with given initial and boundary conditions are directly related to considered physical quantities (like stress, strain, heat, temperature etc.) and space variables (displacement, time etc.) but not to potential function. Bahar and Hetnarski [3] discussed the limitations of using potential functions.
Solutions are found by state-space approach where Cayley-Hamilton theorem is applied to the co-efficient matrix of the field variables.
Sometimes solutions are found by Eigen value approach, as in Das and Lahiri [6] in which using the basic equations a vector matrix differential equation is formed and the solutions for the resulting field equations are found by evaluating the Eigen values and Eigen vector of the corresponding coefficient matrix.
The use of memory dependent derivative (MDD) in the Fourier law of heat conduction means that the heat conduction equation and the constitutive equations are modified with kernel function and the time delay in which they can be chosen freely which provides more flexible approaches to describe a material's practical response.
According to Wang and Li [19], as a new alternative to the fractional order derivative, which has certain features that lead to difficulties when applied to real world problems, the concept of memory-dependent derivative is developed.
For a differentiable function ( ): is called the 'first order memory-dependent derivative' at t relative to the time delay > 0 , to denote memory dependence as, +̇= − ∇ , being the heat flux relative to the temperature gradient and the kernel function ( − ) being differentiable at and can be chosen freely as follows: where, and are constants.
The present study deals with the concept of "Memory dependent derivative" in heat transfer process in an elastic solid body in absence of heat source in the context of two-phase-lag model of generalized thermo elasticity. Two integral transform such as Laplace transform for time variable and Fourier transform for space variable are employed to the governing equations to formulate vector-matrix differential equation which is then solved by the eigen value approach methodology. The upper and lower surfaces of the plate are shear-stress free whereas a time dependent compression is given to the both surfacesand the entire elastic medium is also rotating with a uniform angular velocity. The inversion of Fourier transformation is carried out using suitable numerical techniques while numerical inversion of Laplace transformation is done by Bellman's method. Numerical computations for displacement, thermal strain and stress component, temperature distribution are evaluated and presented graphically under influences of different physical parameters.

Mathematical model
As discussed in Chakraborty, Das and Lahiri [5], in the fields relating to generalized magneto thermoelasticity with memory-dependent derivative, the displacements and the thermal fields as well as the stress-strain-temperature relations for a linear isotropic medium take the following forms: The equation of motion: where, = × = 0 ( × ) represents the Lorentz's force. The generalized heat conduction equation with memory dependent derivative (MDD) in absence of heat source: The constitutive stress components are: The strain displacement components are: where, electro-magnetic field components are satisfying the Maxwell's equation which is given by: Conventionally, we can conclude the following four renowned theories such as: For classical thermoelasticity (CD) theory [1]: 1 = 2 = = 0 = 0 = 0.

Formulation of physical problem
Here, the problem is considered for a homogeneous and rotating isotropic infinitely extended thick plate of a finite thickness 2L along the -axis occupying region = {( , , ): − < < , −∞ < < ∞}, as shown in Fig. 1.
We get the displacement components as = ( ( , , ), ( , , ), 0) with cubical dilatation, electric and induced magnetic intensities are given by: The heat conduction equation with MDD is: and the constitutive stress components are: We now define the following non-dimensional variables: ).
Using the non-dimensional variables and for convenience omitting '*' symbol, we can get the non-dimensional form of the Eqs. (5-10) such as: where: where, and are constants.
For the solution of the vector-matrix differential equation (23), we now apply the method of Eigen value method as in Das, Chakraborty and Lahiri [9] also in Das and Lahiri [8]. The characteristic equation of the matrix ( , ) can be written as: The roots of the characteristic Eq. (24) are: = ± , where, = 1(1)3, which are also the eigenvalues of the matrix and the corresponding eigenvector corresponding to the eigenvalue can be calculated as: = [ ] = ± , = 1(1)3, where, the values of ' , = 1(1)5, are given in Appendix.
Considering the regularity condition at infinity, as in Das, Chakraborty and Lahiri [9] also in Das and Lahiri [8], the general solution of the vector-matrix differential Eq. (23) can be written as: and Eq. (16) also gives the expressions for stress components such as:

boundary conditions
The boundary conditions are described as follows:

Thermal boundary condition
The temperature gradient is zero to the both thermally insulated surfaces = ± of the plate that is:

Mechanical boundary condition
The upper and lower surface of the plate is shear-stress free whereas a time dependent JOURNAL OF ENGINEERING AND THERMAL SCIENCES. DECEMBER 2022, VOLUME 2, ISSUE 2 compression is given to the both surfaces i.e.: = − 0 ( , ) , = ± , = 0, = ± .
To get the values of ′ , we apply Eq. (11) and (19) to the Eqs. (27) and (28) also using Eqs. (25) and (26), we obtain a system of simultaneous equations satisfied by the arbitrary constants, which are given in Appendix.

Numerical results and discussions
The inversion procedure of Laplace and Fourier transform for the expressions of displacement components, temperature distribution and stress components in space-time domain are very much complicated. For this purpose, we prefer to use an efficient computer programme using suitable programming language for the inversion of two integral transforms. Inversion of Laplace transform for the time variable is carried out numerically using Bellman [2] method whereas the inversion of Fourier transform for the space variable is carried out by Gaussian formula for fixed values of and .
In the process of numerical calculation, we consider the homogeneous, isotropic, thermally and electrically conducting material medium as that of copper. The material constants (in SI unit) are given by = 7.76×10 10 ; = 3.    In Fig. 9, the variation of temperature (Θ) is found for the different . The magnitude of temperature decreases up to ≤ 0.2 and no significant changes occur in the region 0.2 ≤ ≤ 0.35.

Application
As a new alternative to the Integral Order Calculus and Fractional Order Calculus, having certain features that lead to difficulties when applied to real world fields, such as geophysics, for understanding the effects of the earth's magnetic field, seismic waves, emission of electromagnetic radiations from nuclear devices, engineering, high energy particle accelerators, nuclear reactor's designing etc., the concept of Memory-Dependent Derivative was developed. We get more accurate results than the others while using Laplace and Fourier Transformations with time parameters and kernel function with the time-delay parameter.
The results will guide us to slove a 2-D problem with dynamic response of memory-dependent derivatives in homogeneous isotropic magneto-thermoelastic medium with two temperatures which is advantageous to successful applications of memory dependence in heat conduction.

Conclusions
The main goal of this article is to introduce a unified new model of generalized magneto thermoelasticity theory with time delay (> 0) and the kernel function ( − ) using the definition for reflection of the memory effect. The form of the kernel function for time delaying memory effect can be chosen freely according to the various applications of the problems. Graphical representations are also explained the necessary importance for the stress distributions in the medium of the plate. Therefore, we can conclude that 1) Thermally and mechanically disturbances are very clear in the neighboring region of the upper and lower surfaces of the plate. No significant disturbance occurred in the middle plane or neighboring region of the middle plane.
2) Significant variations of temperature ( Θ ) are observed for different choices of kernel functions. To analyze the time delaying memory effect, it is more fruitful to us to deal with different value of kernel function and it is also more realistic and eventually more compatible to the physical aspects.
3) This problem may be reduced to the simple generalized magneto thermo-elastic model using more conventional Fourier's law of heat conduction equation defining the values of kernel function as constant and one can compared this result as existing literature.