Solution of Moore–Gibson–Thompson Equation of an Unbounded Medium with a Cylindrical Hole

In the current article, in the presence of thermal and diffusion processes, the equations governing elastic materials through thermodiffusion are obtained. The Moore–Gibson–Thompson (MGT) equation modifies and defines the equations for thermal conduction and mass diffusion that occur in solids. This modification is based on adding heat and diffusion relaxation times in the Green–Naghdi Type III (GN-III) models. In an unbounded medium with a cylindrical hole, the built model has been applied to examine the influence of the coupling between temperature and mass diffusion and responses. At constant concentration as well as intermittent and decaying varying heat, the surrounding cavity surface is traction-free and is filled slowly. Laplace transform and Laplace inversion techniques are applied to obtain the solutions of the studied field variables. In order to explore thermal diffusion analysis and find closed solutions, a suitable numerical approximation technique has been used. Comparisons are made between the results obtained with the results of the corresponding previous models. Additionally, to explain and realize the presented model, tables and figures for various physical fields are presented.


Introduction
The theory of thermoelasticity has received a lot of interest from researchers and scientists due to its many applications in different fields. The fields of architecture, structural features, plasma physics, geophysics, aeronautics, missiles, steam turbine generators, etc. are among the most important areas in which this theory is used. For the first time, an uncombined principle of thermoelasticity was foreseen, in that the elastic strain is independent of the heat transfer and conversely. However, this hypothesis became invalid in displaying the practical results of many concrete problems. Then, researchers provided alternative theories of coupled thermoelasticity, which became widely known as a general thermoelastic theory. This concept has been gaining fame in recent years because its aim is to solve the contradiction of the unlimited heat propagation rate. This description, however, has become appropriate to some extent to communicate practical performance, such as high-speed energy transportation and low-temperature and high-heat transfer engineering.
The classical coupled thermoelasticity theory (CTE) [1] suggested by Biot expects a theoretically unrealistic unlimited velocity of the spread of heat. In order to remove the contradiction of the extraordinary physical phenomenon of unlimited speed in the CTE theory, non-classical thermoelasticity models known as generalized thermoelasticity have been formulated. Lord-Shulman [2], Green-Lindsay [3], and Green-Naghdi [4][5][6] theories The thermo-diffusion theory and the combination of quasi-stationary diffusion problems for an elastic substrate have been investigated by several researchers [32,[39][40][41][42][43][44]. The effects of cross-effects caused by the mixture of temperature, mass diffusion, and strain were studied, resulting in additional mass concentration due to thermal excitation and producing additional temperature ranges [45][46][47][48][49][50].
It was found from previous studies that the theory of the third form (GN-III) by Green and Naghdi has a similar defect as the usual Fourier model and also forecasts the immediate spread of thermal waves. Giorgi et al. [45] explained that such phenomena relevant to constant states are outside the control of the theory of Type III, which should be seen in Fourier as a principle of thermal conduction rather than as a simple theory. This theory is true for both stationary and sluggish thermal phenomena. To solve this dilemma, it is, therefore, natural to amend this suggestion. A new theory of thermal diffusion is formulated in this present work, in which the Moore-Gibson-Thompson equation defines the heat conduction and diffusion formulae. This model was designed to explore the interaction between elasticity, heat, and the mechanisms of diffusion of elastic materials that allow the propagation of thermal waves at finite rates. Many studies have been performed on the inversion of the Laplace transform in the literature [51][52][53][54][55][56][57][58]. More details for them will be given in the next sections. Indeed, the governing system equations can be realized in the new model after the addition of two relaxation factors in the GN-III model. Discrete singular convolution is also applied for some heat and cavity problems [59,60]. A modified and correct solution technique is proposed using the integral boundary formulations of the heat equation by Chernov and Reinarz [61].
The topic of thermo-diffusion interactions in an unbounded isotropic homogeneous solid with a cylindrical hole has been researched in the sense of the generalized model of Moore-Gibson-Thompson thermo-diffusion (MGT-TD). The surface of the cylinder is free of traction and, respectively, undergoes time-dependent convection and chemical load. By means of the Laplace transform method, an accurate solution to the issue is obtained first. The Laplace transformations have been reversed numerically. To show the diffusion effects and different physical phenomena of these solids, the study findings have also been numerically measured and graphically portrayed. Numerical values have been provided in figures and tables to illustrate the comparisons between the physical fields in order to allow a distinction between the results we obtained and the corresponding results in other special models.

The Basic MGT Thermo-Diffusion Equations
Many mathematicians, physicists, and engineers, as well as industry, widely use Fourier and Vic formulas to explain the thermal conductivity and diffusion in elastic materials. In this section, we will derive a new paradigm that allows describing the phenomena of heat transfer as well as propagation within materials in a manner consistent with the physical and chemical aspects.
In a homogeneous isotropic, elastic solids, the governing equations, and constitutive relationships for generalized thermo-diffusion behavior are established by [33][34][35]. The Fourier's law: The Fick's law: The strain-displacement relations: The coupled conservation heat energy equation: The continuity equation [46]: The chemical potential P: where heat flux q, θ = T − T 0 denotes the temperature increment in which T is the absolute temperature, T 0 is the reference temperature, K indicates the thermal conductivity, e ij denotes the strain tensor and u i is the displacement vector. Additionally, in Equations (1)-(6), P denotes the chemical potential, η denotes the flow of the diffusing mass vector, D is the diffusion coefficient, where div u = e kk is the cubical dilatation, Q is the source of heat, ρ is the density of the medium, C e denotes the specific heat at constant strain, a is the measure of thermoelastic diffusion effect, C is the concentration of the diffusive material and β 1 = (3λ + 2µ)α t are the material constants (thermoelastic coupling), where α t is the coefficient of the thermal expansion.
Although the elastic thermal diffusion laws of Fourier and Fick have been well tested for most practical problems, they do not specify a temporary short-time temperature area (higher frequencies and smaller wavelengths). We obtain parabolic partial differential equations through the combination of Equation (1) with (4) and (2) with (5). They emit an infinite velocity of propagation as a result. This property, from a physical point of view, conflicts with physical phenomena. In the past 3 decades, the non-classical diffusion and thermal elasticity theories have been replaced by more general equations to solve the previous inconsistency issue. The generalized model of thermo-diffusion with a time of relaxation was developed by Sherief et al. [36], which allows thermal wave propagation to be reduced.
In 1948, to solve the infinitely rapid propagation in the heat equation, Cattaneo proposed an updated equation for the heat equation. The relationship (1) was replaced by where τ 0 is a thermal relaxation time parameter. Sherief et al. [36] postulated a similar mass flow equation in the same way and similar to Equation (7), given by: where τ 1 indicates the diffusion time relaxation parameter. This ensures that the C equation also forecasts the limited velocity of spread of the substance from medium to media. This is achieved with concentration C. The development of an alternative construction of heat diffusion was then suggested by Green and Naghdi [4][5][6] as an entirely new thermoelastic theory. The function ϑ is seen as a new constituent vector in the models of Green and Naghdi as a gradient of thermal displacement. In mechanical fields and in thermal fields, the scalar function ϑ is generally seen as the mechanical displacement equivalent. In addition, the function ϑ fulfills . ϑ = θ. The law of improved thermal conductivity of the GN-III theory is described as [5] In this instance, the parameter K * > 0 is a constant material property and is also mentioned as the heat conductivity rate. Abouelregal [39,40] suggested a similar relation for the vector of heat flow to Equation (9) as: Mathematics 2021, 9, 1536 5 of 21 In Equation (10), the chemical displacement rate υ is regarded as a new constitutive variable in which parameter D * is the diffusion rate factor. The chemical displacement function υ satisfies the relation . υ = P. The combination of Fourier's modified law (9) and energy Equation (4) has been found to lead to a series of elements in the point continuum, rendering the actual component infinitely dependent on the solutions. Equation (9) has the exact same flaw as the usual Fourier hypothesis that the waves of thermal conduction spread instantly. Therefore, this proposal was also generally updated, and Quintanilla [15,16] included a relaxation time to solve this problem, given the MGT equation. Quintanilla [15] made an adjustment to the proposed updated heat conduction equation (MGTE) after adding the relaxation factor in the Green-Naghdi form III model as follows: By differentiating Equation (11) and using the relationship . ϑ = θ, we have As in Equation (11), the mass flux Equation (10) is assumed to be similar: After differentiating the previous relationship with respect to time and using . υ = P, we obtain By combining Equations (4) and (12), a modified type of heat conduction equation can be obtained, which is developed based on the MGTE equation [15,16]. When we take a divergence operator to Equation (12) and use Equation (4), we obtain: We also obtain the modified equation of mass diffusion by taking the divergence of Equation (15) and the use of Equations (5) and (6) as Complementing the system of equations that govern the behavior of thermal, dynamic, and diffusion propagation inside bodies with normal properties and homogeneity, we have [36]: where and β 2 is a material constant (diffusion coupling) given by β 2 = (3λ + 2µ)α c , where α c is the coefficient of the linear diffusion equation, σ ij is the stress tensor, λ and µ are Lame's constants, and δ ij is the Kronecker delta function. Equation (15) occurs in viscous thermal soothing applications, with high intensity, medical and industrial applications such as lituosis, thermal therapy, or ultrasound washing. The model is also known as the Kelvin or Zener model. This equation is found in the viscoelasticity theory, the classic linear viscoelastic model, to clarify the conduct of certain viscoelastic materials, including the complex and viscous micro-structures of fluids [48].

Implementation of the Proposed Model
We will study an infinitely homogeneous thermo-diffusive solid with a cylindrical hole of radius R and the boundary initially in a non-turbulent state and at a uniform temperature T 0 and an initial concentration of C 0 . We will use the cylindrical coordinate system (r, Θ, z) because of the nature of the problem, considering that the z-axis lies along the axis of the cylindrical hole. The state of a solid can be represented in terms of state r and time t variables because of the symmetry of the problem so that the values of all physical field variables vanish at infinity. We also assumed that there are no sources of external body force and heat that influence the body, implying that in the beginning, the body is calm.
The displacement and strain tensor components are given by Based on Equation (20), we can obtain the cubical dilatation as e = e rr + e ϕϕ + e ΘΘ = ∂u ∂r From Equations (6) and (18), the radial and hoop stresses and chemical potential can be written as The modified equations for MGT heat conduction (15) and MGTE mass diffusion (16) can be written as follows: where ∇ 2 = ∂ 2 ∂r 2 + 1 r ∂ ∂r is the cylindrical operator of Laplace. Equation (17) can be expressed as Inserting Equations (22) and (23) in Equation (27), we obtain Multiplying the cylindrical div operator on both sides of Equation (28), then we have By using the following dimensionless quantities, the control equations can be more conveniently put into practice: Equations (22)- (26) and (29) in the non-dimensional variables can be expressed as follows: where We have lowered the primes for simplicity in Equations (31)-(36).

Initial and Boundary Conditions
The initial conditions can be assumed as The following boundary conditions can be considered where θ 0 and P 0 are constants, and H(t) is the Heaviside unit step function.

Solution in the Laplace Transform Domain
Applying the Laplace transform to Equations (31)-(36) under the initial conditions (38), described by we obtain where Eliminating e, C between Equations (43)-(45), we obtain where By the same manner, we can demonstrate that the functions e and C can fulfill the equations Introducing the parameters m i (i = 1, 2) into Equation (50), we have where m 2 1 , m 2 2 and m 2 3 are the solutions for the equation: The parameters m 1 , m 2 and m 3 can be determined as The solution of Equation (54), which is infinitely limited, is given by where B i , (i = 1, 2, 3) are the integral parameters and K 1/2 (·) denotes the second kind of Modified Bessel function of zero order. Likewise, we can write Inserting Equations (58) and (59) into Equations (43)- (45), the following relations can be obtained Then, Equations (58) and (59) may be expressed as After using the relationship between u and e, the solution of the radial displacement u will be Using the expressions (57), (61)-(63) and the well-known relation we can obtain the solutions of the thermal stresses and the chemical potential as The boundary conditions in Equations (39)-(41) at r = R in the transformed domain may be expressed as Using the boundary conditions (68) with Equations (57), (65), and (67), the following linear system of equations can be obtained: The unknown parameters B 1 , B 2 , and B 3 can be determined by resolving the system (69)-(71). By the setting of the constants B 1 , B 2 , and B 3 and substitution in the general solutions of the different functions, this means that we have finished solving the problem in the field of the Laplace Transform.

Generalized MGT Model of Thermoelasticity without DIFFUSION
In this case, we take C = P = 0. However, the functions θ and e may be obtained Laplace transform domain by solving the two differential equations: Eliminating e(r, s) or θ(r, s) in the above equation, one obtains the fourth-order differential equation: where The solutions to Equation (55) are given by where are integral constants and n 2 i are given by Based on these solutions, the quantities of other physical fields can be simplified in the same way as in the previous analysis: To complete the solution in this case, the boundary conditions in Equations (39) and (40) are sufficient to obtain the integral constants G i (i = 1, 2).

Special Cases of Thermoelasticity and Thermo-Diffusion Models
The model and equations derived and obtained in the second section of this paper are valid for many special cases that can be inferred from our constructed model. A system of equations containing thermal conduction equation (Equation (15)) and diffusion equation (Equation (16)) may include at least 10 generalized theories of thermal elasticity in the absence of diffusion as well as some generalized diffusion theories. The following can be summarized and documented in some special cases: Case I: When the diffusion effect is absent If the diffusion effect is missing (C = 0), it is possible to obtain the following models: • The traditional theory of thermoelasticity (CTE) [1] when K * = τ 0 = 0. • The Lord and Shulman generalized theory (LS) [2] by setting K * = 0.
• The Green and Naghdi model of Type II (GN-II) [6] when the first term on the righthand side of Equation (15) is disregarded and τ 0 = 0. • The Green and Naghdi theory of Type III (GN-III) [5] when τ 0 = 0.

•
The generalized thermoelastic diffusion model (GTD) [36] by setting K * = 0 and D * = 0. • The Green and Naghdi thermoelastic diffusion theory of Type II (GND-II) [49] when the first term on the right-hand side of Equation (15) is neglected with the addition of The Green and Naghdi thermoelastic diffusion theory of Type III (GND-III) [50] when thermal relaxation τ 0 = 0.

Inversion of the Laplace Transforms
Many problems have solutions that may be expressed in terms of a Laplace transform, which is therefore too difficult to invert using complex analysis techniques. The numerical evaluation of the Laplace inversion integral has been accomplished using a variety of approaches. Many papers have been written in the literature on the inversion of the Laplace transform. See, for example, [51][52][53][54], and a detailed bibliography can be found in [55,56].
Some methods work amazingly well with some image functions, while providing horrific results for others. Since no approach is best in all cases, we recommend using several alternative strategies to solve a particular reversal problem. We can be more confident and accept a numerically inverse Laplace transform if two or more methods give approximately the same result. In [57], Abate and Valkó have classified those algorithms into four categories according to the basic approach of the method as follows: (1) Fourier series expansion; (2) Laguerre function expansion; (3) combination of Gaver functionals; (4) deform the Bromwich contour.
This work presents the numerical inversion technique of Laplace transforms based on the Durbin and Fourier series expansion [58] and Honig and Hirdes [53]. The disadvantage of this type of inversion method, the remarkable dependence of estimation, and truncation error on free parameters is overcome by applying a procedure to reduce the estimation error, a method to speed up the Fourier series convergence, and a procedure that roughly computes the "best" option for all free parameters at the same time. For two reasons, this approach is effective in meeting the demands of standalone digital computing. Compared to other more complex numerical squaring methods, it is fast (economical) on the digital computers that are currently accessible. Second, because the inverse function is given as a Fourier cosine series with the coefficients being acceptable forward transform values, the technique is theoretically straightforward and requires minimal programming effort.
The Laplace transform and its inversion formula of the function f (t) are defined as follows: where ϑ > 0 is arbitrary but must be chosen so that it is greater than the real parts of all the singularities of F(r, s). For purposes of discussion, it is herein assumed that these integrals exist for Re(s) > ϑ > 0. Using the inversion formula which is equivalent to (2) with s = ϑ + iw, ϑ, w ∈ R, and a Fourier series expansion of e −ϑt f (r, t) in the interval [0, 2t 1 ], Durbin derived the approximation formula where the parameter N is a satisfactorily large integer and denotes the number of stops in the truncated Fourier series. The parameter N can be chosen so that it is e e ϑt Re e iNπt/t 1 g(c + iNπ/t 1 ) ≤ ε 1 where ε 1 is a persecuted small positive number that corresponds to the degree of accuracy to be achieved. A good choice of free parameters N and ϑt is important not only for the accuracy of the results but also for the application of the "Korrektur" method and convergence acceleration methods. These methods do not improve results if the parameters are badly chosen. For faster convergence, numerical experiments have shown that the value of the parameter ϑ that satisfies the above relation is given in terms of time as ϑ =≈ 4.7/t [13].

Numerical Example and Discussion
We will now present some empirical results numerically to explain the new model introduced in this work and on the basis of preceding section outcomes. In order to answer the issue more comprehensively and to see how the studied fields depend on diffusion and certain physical influences, we have also numerically calculated for a particular physical material. The physical values of copper were selected for the purpose of mathematical estimation. The properties of this substance are therefore given in the SI units as follows [39,40]: All of the parameters calculated t = 0.12, γ 1 = 1 and γ 2 = 0.5 will be used in calculations unless are determined otherwise. The numerical procedure defined in (84) has been applied to obtain the temperature θ, displacement u, radial stress σ rr , hoop stress σ ΘΘ , and concentration C as the chemical potential P distributions within the medium. We assume that the cylinder cavity, with its center, has a radius of R = 1 from the origin.
By comparing the basic Green and Naghdi theories of the second and third types with the modified thermal diffusion model, this section will provide some examples and applications for testing the validity of the given model. Plot activities are seen in the hollows cylinder's radial pathway. The tabulation of some of the results for future comparisons to other researchers is an important part of this study.
The effects of the field variables analyzed based on the thermo-diffusion theories CTED, CTED, GND-II, GND-III, and MGTED have been reported in Tables 1-6. In the radial direction of the cylinder, Figures 1-6 display further sample graphs of all models to assess the effect of various thermal diffusion models on the quantities of the physical fields.
In the presence of the diffusion effect, Table 1 and Figure 1 display the variance of the temperature vs. the radial distance r for the different theories. The temperature values θ reach the maximum values on the surface of the spherical cavity r = 1 as shown in the table and the figure, and then gradually decrease. The distribution of temperature satisfies the condition of the thermal boundary imposed upon the problem. The values of θ gradually decrease and ultimately decrease in the opposite direction of the heat-wave propagation to a value of zero as the radius r increases. the temperature vs. the radial distance for the different theories. The temperature values reach the maximum values on the surface of the spherical cavity = 1 as shown in the table and the figure, and then gradually decrease. The distribution of temperature satisfies the condition of the thermal boundary imposed upon the problem. The values of gradually decrease and ultimately decrease in the opposite direction of the heat-wave propagation to a value of zero as the radius increases.  A great difference in the temperature field has been found near the boundary surface of the body under all thermo-diffusion theories, considering the convergence of all models when increasing the radius within the body. It was also noted that in the CTED and GND-III models, higher temperature values were observed compared to the GN-II, GTD, and MGTED models. In the CTED and GND-III models, thermal waves do not fade easily, unlike in other thermal diffusion models in which it is clear that thermal waves spread at slow speeds inside the medium. This is because the heat waves in the CTED model propagate at an unlimited velocity. It is also noted that in the case of the MGTED model, the A great difference in the temperature field has been found near the boundary surface of the body under all thermo-diffusion theories, considering the convergence of all models when increasing the radius within the body. It was also noted that in the CTED and GND-III models, higher temperature values were observed compared to the GN-II, GTD, and MGTED models. In the CTED and GND-III models, thermal waves do not fade easily, unlike in other thermal diffusion models in which it is clear that thermal waves spread at slow speeds inside the medium. This is because the heat waves in the CTED model propagate at an unlimited velocity. It is also noted that in the case of the MGTED model, the numerical values are lower inside the body away from the spherical cavity of the surface than in the case of the GND-III model, and the waves also fade away faster. The explanation for the rapid decrease, as predicted, is the inclusion of relaxation times in the equations extracted.
The effect of all thermal diffusion models on the displacement u in the radial direction of the cylinder is shown in Table 2 and Figure 2. The displacement distribution begins at a negative value and then increases gradually as the distance increases until it reaches a peak value of approximately r = 1.2, after which it eventually decreases to zero with r increasing. From the table and figure, it can be seen that the relaxation times τ 0 and τ 1 that were used in the current thermo-diffusion model have a significant role to play in body deformation. at a negative value and then increases gradually as the distance increases until it reaches a peak value of approximately = 1.2, after which it eventually decreases to zero with increasing. From the table and figure, it can be seen that the relaxation times and that were used in the current thermo-diffusion model have a significant role to play in body deformation.  In the case of modified GTD, GND-II, and modified MGTED models, numerical values of disolactivity are also smaller compared to the findings of the CTED and GN-III models. In general, it seems that in various theories, there is great similarity and convergence in the action of displacement. Only in magnitude and peak points does the difference arise. It can be detected that when relaxation times are used, the displacement profile In the case of modified GTD, GND-II, and modified MGTED models, numerical values of disolactivity are also smaller compared to the findings of the CTED and GN-III models. In general, it seems that in various theories, there is great similarity and convergence in the action of displacement. Only in magnitude and peak points does the difference arise. It can be detected that when relaxation times are used, the displacement profile for the Green and Naghdi thermo-diffusion model (GND-III) is greater in the absence of relaxation times than for the MGTED thermo-diffusion model.
For various combinations of thermo-diffusion models, Tables 3 and 4 and Figures 3 and 4 are provided to investigate the effect of thermo-diffusion on the variance of thermal stresses σ rr and σ ΘΘ . The comparison between the various models and the various distributions of heat stress can be observed. It is also noted that, in comparison to the CTED model, the magnitudes of stress values in the case of thermal diffusion GN-III models are the largest possible.     The radial thermal stress starts at zero, filling the conditions of the proposed problem, then decreases rapidly until at the peak point it reaches its maximum value, and then decreases steadily until it disappears to zero (see Table 3 and Figure 3). With the hoop stress , the same previous behavior occurs except that it starts with a negative value different from zero at the surface of the spherical cavity. In behavior and principles, the variations in GTD-, GND-II-, and MGTED-generalized thermal diffusion models are similar together.
The existence of the relaxation parameters of generalized thermoelasticity (GND-III) and the diffusion that have been integrated into the current model (MGTED) have a major effect on the properties and actions of the vibrations of the different distributions. The radial thermal stress σ rr starts at zero, filling the conditions of the proposed problem, then decreases rapidly until at the peak point it reaches its maximum value, and then decreases steadily until it disappears to zero (see Table 3 and Figure 3). With the hoop stress σ ΘΘ , the same previous behavior occurs except that it starts with a negative value different from zero at the surface of the spherical cavity. In behavior and principles, the variations in GTD-, GND-II-, and MGTED-generalized thermal diffusion models are similar together.
The existence of the relaxation parameters of generalized thermoelasticity (GND-III) and the diffusion that have been integrated into the current model (MGTED) have a major effect on the properties and actions of the vibrations of the different distributions.
It can be attributed to the compressive essence of stress that atoms are suppressed together in the solids and voluminously induce compressive stress because of the dissolution of the substance. In both cases, the comparison of each figure shows that, as time pass away from the source, the magnitude of the various quantities considered decreases, thus showing the wave fronts.
For the five separate thermo-diffusion models, the numerical computations of the chemical potential P against the radial distance r are shown graphically in Figure 5 and in Table 5. Figure 5 and Table 5 show that at the beginning, the variance of the chemical potential P is great and that the vibrations of the transition decrease and fade as the r value increases. The chemical potential P profile has also been established, starting with the greatest surface value, meeting the boundary conditions imposed on the problem, and this confirms and validates the outcomes we have obtained. In Figure 5 and Table 5, we note that the variations in the case of CEDT and GND-III are similar to each other and have the highest values, while the GTD models GND-II and MGTED are smaller in size due to the effect of generalized thermal and diffusion relaxation parameters.
The comparison with space for the distribution of concentration C is shown in Figure 6 and Table 6 for different thermos-diffusion models. From the numerical values, it was found that there were some variations between the values of concentration C and the variation in the model. The concentration C curves and tables have the same behavior but vary only in quantity.    The comparison with space for the distribution of concentration is shown in Figure 6 and Table 6 for different thermos-diffusion models. From the numerical values, it was found that there were some variations between the values of concentration and the variation in the model. The concentration curves and tables have the same behavior but vary only in quantity.  The concentration begins with a positive value at the cavity of the cylinder and then rises steadily until = 1.2 reaches its highest local maximum value. After that, increasing the distance r within the medium steadily decreases until it dissipates to obtain its null value. In contrast with the CTED and GND-III models, the -field values can be inferred  The concentration begins with a positive value at the cavity of the cylinder and then rises steadily until r = 1.2 reaches its highest local maximum value. After that, increasing the distance r within the medium steadily decreases until it dissipates to obtain its null value. In contrast with the CTED and GND-III models, the C-field values can be inferred from the numeric values to be smaller in the case of the MGTED model. The GND-III model, however, gives the GND-II model a larger model.
The medium is strained in the vicinity of the cylinder and, with the passing of time, is greater. Cross-effects may result in the presence of tensile stress near the cavity surface due to temperature coupling, mass dissemination, and strain fields. Owing to these crosseffects, an additional accumulation of thermal stress arises from the thermal excitation. The primary difference between the GND-III model of thermoelastic diffusion and the MGTED-based thermal conductivity equations of thermoelastic diffusion is the additional delay times τ 0 and τ 1 .
These relaxation times provide an additional diffusion and propagation mechanism in the governing equation for the dissipation effect. The delay time τ 0 is known to monitor the propagation behavior of the thermal wave, slow down the speed of heat-wave propagation, and show the characteristics of the thermal wave. Thermal energy can be diffused, and diffusion wave decay characters can be produced in the MGTED heat transfer by the τ 1 effect.

Conclusions
A new generalized thermoelastic diffusion model has been derived in the present paper, connecting heat and mass flux in elastic materials. It is understood that, due to the coupling between temperature and mass diffusion as well as stress fields, thermoelastic diffusion takes place in a flexible solid material. Centered on the Moore-Gibson-Thompson equation, the constitutive equations, as well as the heat and mass diffusion equations, have been modified. The goal behind the implementation of this new model was to deal with the evident inconsistency between the infinite thermal and diffusion rates predicted by the classical theories of thermoelastic diffusion [1,[33][34][35] and the energy dissipation model of Green and Naghdi (GN-III) [5]. In the constructed modified model, the laws of Fourier and Fick have also been improved to incorporate the time needed to speed up the heat wave and the time derivative of the diffusive mass flow.
The primary distinction between the thermoelastic diffusion model GND-III and the thermoelastic diffusion thermal conductivity equations based on MGTED is the additional delay times τ 0 and τ 1 . These relaxation times provide an additional diffusion and propagation mechanism in the governing equation for the dissipation effect. The delay time τ 0 is known to monitor the propagation behavior of the thermal signal, slow down the speed of thermal wave propagation, and show the characteristics of the thermal wave.
Thermal energy can be diffused, and diffusion wave decay characters can be produced in the MGTED heat transfer by the τ 1 effect.
The problem of an unbounded solid with a cylindrical cavity with a permeable material shaped on the surface in contact with the cylindrical cavity has been studied to explain and examine the proposed model. In the area of Laplace transformation, the solution is achieved by a direct approach. Centered on Fourier expansion techniques, a computational approach is used to invert Laplace transforms to solve the problem in the physical field.
It was found from our observations that thermal and diffusion waves are propagated as a wave with a slight velocity rather than an endless velocity on the conventional thermoelastic diffusion model. The thermoelastic diffusion models can be deduced as special cases: Green-Naghdi Type II and III models, generalized thermo-diffusion with one relaxation. Furthermore, the GND-III and CTED models have a near-conduct but are distinct from MGTED, which shows that the model of thermoelastic diffusion was significant. The results in this article should be useful for science and engineering research as well as for those working on solid mechanics growth. A wide variety of thermodynamic concerns are related to the method used in this article. These theoretical conclusions will provide experimental experts and researchers who work on this subject with interesting knowledge.