Next Article in Journal
Discrete Optimization: The Case of Generalized BCC Lattice
Previous Article in Journal
MooFuzz: Many-Objective Optimization Seed Schedule for Fuzzer
Previous Article in Special Issue
Dissipativity of Fractional Navier–Stokes Equations with Variable Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Fractional Time Derivative on Two-Dimension Porous Materials Due to Pulse Heat Flux

1
Mathematics Department, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department Mathematics, Faculty of Science, Sohag University, Sohag 82524, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(3), 207; https://doi.org/10.3390/math9030207
Submission received: 21 December 2020 / Revised: 6 January 2021 / Accepted: 14 January 2021 / Published: 20 January 2021

Abstract

:
In the present article, the generalized thermoelastic wave model with and without energy dissipation under fractional time derivative is used to study the physical field in porous two-dimensional media. By applying the Fourier-Laplace transforms and eigenvalues scheme, the physical quantities are presented analytically. The surface is shocked by heating (pulsed heat flow problem) and application of free traction on its outer surface (mechanical conditions) by the process of temperature transport (diffusion) to observe the full analytical solutions of the main physical fields. The magnesium (Mg) material is used to make the simulations and obtain numerical outcomes. The basic physical field quantities are graphed and discussed. Comparisons are made in the results obtained under the strong (SC), the weak (WC) and the normal (NC) conductivities.

1. Introduction

Porous media appear in many forms of environmental, natural, and synthetic implementations and in several technologies. To overcome the first insufficiency in the decoupled thermoelasticity theorem, in 1956, Biot [1] presented the coupled thermoelasticity theorem to control the first insufficiency in the decoupled thermoelastic model, which prognosticates two phenomena not suitable for physical observation. Firstly, the thermal conductivity equation is parabolic, presenting an infinite propagation speed for thermal waves. Secondly, the thermal conductivity equation of this model does not contain an elastic term. Rosencwaig et al. [2] investigated the local thermoelastic deformation of the model caused by excitations.
Biot developed poroelasticity models [1,3,4] for a high–low-frequency range and built upon the coupled thermoelasticity hypothesis to overcome the inconsistency in the uncoupled hypothesis [1]. The heat conduction and elasticity equations in this theory are coupled. However, it includes a drawback of the uncoupled hypothesis in which the heat wave propagates with an infinite velocity that is impractical in nature. Then, to solve the problem of the coupled hypothesis, generalized thermoelasticity models were expanded. It is recognized that there are several generalizations of the thermoelasticity hypothesis, such as that presented by Lord-Shulman [5]. Green-Naghdi [6,7,8] formulated three types of models (GNI, GNII, and GNIII). The constitutive equations of the GN models are linearized, where the first type is similar to the classical coupled thermoelastic model, the second type demonstrates the propagations of thermal signals with finite velocity without energy dissipations, and the third type proposes the finite speed of propagations with energy dissipations. During the second half of the 19th century, it can be said that the complete model of fractional integrals and derivatives was determined. In the context of generalized thermoelastic models, Youssef [9,10] established generalized fractional-order thermoelastic models under strong, normal, and weak conductivities. Sherief et al. [11] presented a new theory by using the thermal conduction law. Ezzat and Elkaramany [12] established another theory for a generalized fractional-order thermoelasticity model based on Taylor’s expansions of time-fractional order. Ezzat et al. [13] studied modeling in a generalized thermoelastic model using a memory-dependent derivative. Marin [14] discussed basic models in elastostatic micropolar media. Saeed et al. [15] used the finite-element scheme to investigate thermoelastic interactions in poroelastic media using the GL model.
Ouyang et al. [16] proposed three-equation models considering local thermal nonequilibrium states. Many authors [17,18,19,20,21,22,23] have discussed the solutions to some problems using various generalized thermoelastic models. Hussein [24] presented a mathematical model for a spherical thermoelastic porous material. Hobiny and Abbas [25] discussed the effect of a dual-phase lag model in two-dimensional porous media. Biswas [26] investigated surface waves in nonlocal, orthotropic thermoelastic porous media. Carini and Zampoli [27] studied porous matrices with three delay times. Abbas and Marin [28] presented the analytical solution of a two-dimensional generalized thermodiffusion problem subjected to a laser pulse. Shekhar [29] studied deformations caused by thermal shock in a porous thermoelastic medium with properties dependent on the reference temperature. Calin et al. [30] investigated improvements in the rigidity of circular-plate composites. Abd-Elaziz et al. [31] discussed the Thomson impacts on thermoelastic porous materials based on Green and Naghdi’s model. Abbas and Kumar [32] investigated deformations induced by a heating source in the micropolar plane under generalized thermoelasticity using the finite-element scheme.
Ellahi et al. [33,34] studied the solutions of different problems under several boundary conditions in porous media. Singh [35] has investigated the wave propagations in a medium with voids under thermoelasticity models. Palani and Abbas [36] discussed the free convections magnetohydrodynamics flow with thermal radiations. Villatoro et al. [37] have applied the perturbation approach upon Laplace transform to get the solution of the heating equation in porous media which consists of gas and solid phases. The discontinuous front was observed in gas temperature because of incompatibility between the initial and boundary conditions, leading to the constant speed for thermal propagation. However, there was a smooth front at the solid temperature using an internal layer of asymptotic approximation. Abbas [38] discussed the natural frequency of a poro-elastic hollow cylinders. Alzahrani et al. [39] used the eigenvalues approaches to investigate the effects of thermal relaxation times in two-dimension porous media under strong, normal and weak thermal conductivities.
The objective of this work is to study the effects of strong, weak, and normal conductivities in a two-dimensional porous medium by using the eigenvalue scheme. By using the eigenvalue technique with Fourier–Laplace transforms on numerical and analytical methods, basic formulations are presented. The nondimensional temperature, displacements, stresses, and volume fraction are obtained and represented graphically. In the calculations, the impacts of strong, normal, and weak thermal conductivities on the considered variables are investigated and compared.

2. Basic Equations

For an isotropic 2D elastic porous material, the basic formulations based on [9,35] models without body force and the heating resources are expressed as:
( λ + μ ) u j , i j + μ u i , j j + b φ , i γ t Θ , i = ρ 2 u i t 2 ,
α φ , j j b u j , j ζ 1 φ ω o φ t + m Θ = ρ ψ 2 φ t 2 ,
( ( K * + K t ) I ϵ 1 Θ , j ) , j = t ( ρ c e Θ t + m T o φ + γ t T o u j , j t ) , 0 < ϵ 2 .
where the integral operator of fractional derivative is expressed as [9]
I ϵ 1 g ( r , t ) = 1 Γ ( ) 0 t ( t τ ) ϵ 1 g ( r , τ ) d τ ,   { 0 < ϵ < 1 , w e a k   c o n d u c t i v i t y ϵ = 1 , n o r m a l   c o n d u c t i v i t y 1 < ϵ 2 , s t r o n g   c o n d u c t i v i t y ,
where Γ ( ) is the Gamma function. The stress-displacement equations are defined by [35]
σ i j = μ ( u i , j + u j , i ) + ( λ u k , k + b φ γ t Θ ) δ i j .
where m is the coefficient of thermo-void, ψ is the equilibrated inertia, b is the measure of diffusion effects, ω o ,   α , ζ 1 are the parameters of voids material, c e is the specific heat, K * is the model characteristic material constant, Θ = T T o , T o is the reference temperature, ρ is the material density, σ i j are the stress components, K is the thermal conductivity coefficient, μ , λ are the Lame’s parameters, t is the time, u i are the components of displacement, i , j , k = 1 , 2 , 3 , γ t = ( 3 λ + 2 μ ) α t , α t is the linear coefficient of thermal expansions.

3. Formulations of the Problem

We consider a two-dimension porous material fills the region 0 x ,   y . By using the Cartesian co-ordinates ( x , y , z ) , the governing equations with the components of the displacement ( u , v , 0 ) can be given by
( λ + 2 μ ) 2 u x 2 + ( λ + μ ) 2 v x y + μ 2 u y 2 + b φ x γ t Θ x = ρ 2 u t 2
( λ + 2 μ ) 2 v y 2 + ( λ + μ ) 2 u x y + μ 2 v x 2 + b φ y γ t Θ y = ρ 2 v t 2 ,
α ( 2 φ x 2 + 2 φ y 2 ) ω o φ t ζ 1 φ b ( u x + v y ) + m Θ = ρ ψ 2 φ t 2 ,
( K * + K t ) I ϵ 1 ( 2 Θ x 2 + 2 Θ y 2 ) = t ( ρ c e Θ t + m T o φ + γ t T o t ( u x + v y ) ) ,
σ x x = λ v y + ( λ + 2 μ ) u x + b φ γ t Θ ,   σ x y = μ ( u y + v x ) .
The problem initial conditions are defined by
φ = φ t = 0 , Θ = Θ t = 0 ,   u = u t = 0 ,   v = v t = 0 , t = 0 .
While, the problem boundary conditions are given by
σ x x = 0 , σ x y = 0 , ( K * + K t ) Θ ( x , y , t ) x = q o t 2 e t t p 16 t p 2 H ( a | y | ) ,   φ x = 0
where q o is a constant, t p is the time of the flux pulse heating characteristics, and H is the function unit step. For appropriateness, the nondimensional variables can be taken as
Θ = Θ T o ,   φ = ψ η 2 c 2 φ ,   ( σ x x , σ x y ) = ( σ x x ,   σ x y ) ( λ + 2 μ ) , ( t , t p ) = η c 2 ( t , t p ) ,   ( u , v , x , y ) = η c ( u , v , x , y ) ,
where η = ρ c e k and c = λ + 2 μ ρ . In these nondimensional terms of the variables in Equation (13), the basic formulations can be written by (the dashes have been neglected for convenience)
2 u t 2 = 2 u x 2 + ( 1 r 1 ) 2 v x y + r 1 2 u y 2 + r 2 φ x r 3 Θ x ,
2 v t 2 = 2 v y 2 + ( 1 r 1 ) 2 u x y + r 1 2 v x 2 + r 2 φ y r 3 Θ y ,
r 4 2 φ t 2 = 2 φ x 2 + 2 φ y 2 r 5 φ t r 6 φ r 7 ( u x + v y ) + r 8 Θ ,
( r 9 + t ) ( 2 Θ x 2 + 2 Θ y 2 ) = ϵ t ϵ ( Θ t + r 10 φ + r 11 t ( u x + v y ) ) ,
σ x x = ( 1 2 r 1 ) v y + u x + r 2 φ r 3 Θ ,   σ x y = r 1 ( v x + u y ) .
σ x x = σ x y = 0.0 ,   ( r 9 + t ) Θ ( x , y , t ) x = q o t 2 e t t p 16 t p 2 H ( a | y | ) ,   φ x = 0 ,
where r 1 = μ ρ c 2   , r 2 = b ρ ψ η 2 c 4 , r 3 = γ t T o ρ c 2 , r 4 = ρ c 2 ψ α , r 5 = ω o α η , r 6 = ζ 1 α η 2 c 2 , r 7 = b ψ α , r 8 = m ψ T o α , r 9 = K * ρ c e c 2 , r 10 = m ψ η 4 c 4 k , r 11 = γ t ρ c e .
Now, the Laplace transforms for any function f ( x , y , t ) , are given by
f ¯ ( x , y , s ) = 0 f ( x , y , t ) e s t d t ,
however, the Fourier transform for any function f ¯ ( x , y , s ) can be expressed by
f ¯ * ( x , q , s ) = f ¯ ( x , y , s ) e i q y d y ,
Thus, the governing formulations are expressed to obtain the ordinary differential equations with the boundary conditions as follow
d 2 u ¯ * d x 2 = ( s 2 + r 1 q 2 ) u ¯ * i q ( 1 r 1 ) d v ¯ * d x r 2 d φ ¯ * d x + r 3 d   Θ ¯ * d x ,
d 2 v ¯ * d x 2 = ( s 2 + q 2 ) r 1 v ¯ * r 2 i q r 1 φ ¯ * + r 3 i q r 1   Θ ¯ * i q ( 1 r 1 ) r 1 d u ¯ * d x ,
d 2 φ ¯ * d x 2 = r 7 i q v ¯ * + ( s 2 r 4 + q 2 + r 6 + r 5 s ) φ ¯ * r 8   Θ ¯ * + r 7 d u ¯ * d x ,
d 2   Θ ¯ * d x 2 = r 11 i q s + 1 r 9 + s v ¯ * + r 10 s r 9 + s φ ¯ * + ( q 2 + s + 1 r 9 + s )   Θ ¯ * + s + 1   r 11 r 9 + s d u ¯ * d x ,
σ ¯ x x * = i q ( 1 2 r 1 ) v ¯ * + d u ¯ * d x + r 2 φ ¯ * r 3   Θ ¯ * .
σ ¯ x y * = r 1 ( d v ¯ * d x + i q u ¯ * ) ,
σ ¯ x x * = σ ¯ x y * = 0 , d φ ¯ * d x = 0 , d   Θ ¯ * d x = q o t p 8 ( r 9 + s ) ( s t p + 1 ) 3 2 π s i n ( q a ) q ,
Now, the vector-matrix differential equations (Equations (22)–(25)) are written as
d M d x = A M ,
where A and M are defined as shown in Appendix A. By using the eigenvalues method as in [17,40,41,42,43], the exact solutions for the Equation (29) are obtained. So, the matrix characteristic equation for A are taken the form as follow:
ξ 8 f 1 ξ 6 + f 2 ξ 4 + f 3 ξ 2 + f 4 = 0 ,
where f 1 , f 2 , f 3 and f 4 are determined as in Appendix B. To obtain the solution of Equation (29), the eigenvalue of matrix A and its eigenvectors must be computed, where ξ 1 ,   ξ 2 ,   ξ 3 , ξ 4 , ξ 1 ,   ξ 2 ,   ξ 3 and   ξ 4 are the eigenvalues which have the corresponding eigenvectors as in Appendix C. Thus, the analytical solutions of Equation (29) are written as:
M ( x , q , s ) = i = 1 4 B i Y i e ξ i x ,
The positive exponentials can be discarded, which, due to the conditions of the regulations of the solution at infinity and B 1 ,   B 2 ,   B 3 , and B 4 , are constants that are computed using the conditions of the problem boundary. Now, for each function f ¯ * ( x ,   q ,   s ) , the inverse Fourier transform can be expressed by
f ¯ ( x , y , s ) = 1 2 π f ¯ * ( x , q , s ) e i q y d q ,
Finally, to get the general solution for the displacement, the stresses components, the changes in volume fraction field of void distributions φ and the variations of temperature with respect to the distances x , y at any time t , Stehfest [44] numerical inversion schemes were chosen. In these schemes, the Laplace transform inverse for f ¯ ( x , y , s ) are defined by
f ( x , y , t ) = l n ( 2 ) t n = 1 N V n f ¯ ( x , y , n l n ( 2 ) t ) ,
where
V n = ( 1 ) ( N 2 + 1 ) p = n + 1 2 m i n ( n , N 2 ) ( 2 p ) ! p ( 1 + N 2 ) p ! ( n p ) ! ( N 2 p ) ! ( 2 n 1 ) ! ,
where N is the term number.

4. Numerical Results

For numerical examples, the magnesium mediums can be taken to object of numerical calculations. The parameters values of (Mg) are taken from [45]
α = 3.688 × 10 5 ( N ) ,   ζ 1 = 1.475 × 10 10 ( N ) ( m 2 ) , ω o = 0.0787 × 10 3 ( N ) ( m 2 ) ( s 1 ) , μ = 3.278 × 10 10 ( N ) ( m 2 ) ,   λ = 2.17 × 10 10 ( N ) ( m 2 ) , ρ = 1740 ( kg ) ( m 3 ) ,   t = 0.3 , T o = 298 ( k ) ,   c e = 1040   ( J ) ( kg 1 ) ( k 1 ) ,   a = 0.25 ,   ψ = 1.753 × 10 15 ( m 2 ) , α t = 1.98 × 10 6 ( k 1 ) , K = 1.7   ( W ) ( k 1 ) ( m 1 ) ,   β = 2.68 × 10 6 ( N ) ( m 2 ) ( k 1 ) , b = 1.13840 × 10 10 ( N ) ( m 2 ) , m = 2 × 10 6 ( N ) ( m 2 ) ( k 1 ) .
The above data have been used to study the strong (SC), the normal (NC) and the weak (WC) conductivities in 2D porous materials by the eigenvalue method. The voids changes in volume fraction field distribution φ , The variations of temperature Θ , the stresses σ x x ,   σ x y and the of displacement components u ,   v are studied. The material is considered to be homogeneous two-dimensions media. Figure 1 displays the change in the volume fraction field of void distributions φ along x . It is clear that it reduces with the rising x till attaining zeros. Figure 2 displays the temperature variations via x . It is noticed that it starts from heights value according to the application of boundary condition and reduces with the rising x to come to zeros. Figure 3 depicts the variations of vertical displacement with respect to x which have maximum value on x = 0 and reduces with the rising x . Figure 4 illustrates the horizontal displacement variations u via x . It is indicated that it attains maximums negative value and gradually rises till it attains peak value at specific locations in close nearness to x = 0 and after that reduces to come to zeros. As seen in Figure 4, the displacement changes continuously from negative to positive and after that goes down to zero, which is caused by the combined effect of the traction free bounding surface, thermal expansion and finite heat speed. Figure 5 and Figure 6 show the components of stress variations σ x y and σ x x along x . It is noticed that the magnitudes of components, permanently begin from zeros which obeyed the boundary condition.
Figure 7 and Figure 8 show the changes in volume fraction field of voids variations φ and the variations of temperature Θ via the distance y when x = 0.4 . It is indicated that the variations of changes in volume fraction field of void and the variations of temperature have maximum values at the length of heating surfaces ( ( | y | 0.4 ) ) after that begin to reduces totally near the edges ( ( | y | 0.4 ) ) where they reduce and reach to zeros values. Figure 9 displays the vertical displacement variation v via y . It is noticed that it starts the rising at the start and end of the thermal surfaces ( | y | 0.4 ) , and have small values at the center of the thermal surfaces ( | y | 0.4 ) , after that it begins the rising and reach ultimate value totally near the edges ( y = ± 0.4 ) , then it reduces to come to zeros. Figure 10 shows the variation of horizontal displacement u along x . It is observed that the horizontal displacement has a maximum value at the length of the heating surface ( | y | 0.4 ) , and then it starts to decrease totally near the edges ( y = ± 0.4 ) , and then reduces to zeros values. The stresses σ x x and σ x y with respect to y are presented in Figure 11 and Figure 12.
Finally, Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 explain the variations of all the studied fields along the distance y and the distance x at t = 0.3 . These figures display the predict curves during the strong, the weak and the normal conductivities. Unsurprisingly, you can find that the stages of the strong, normal and weak conductivities have great influences on the values of variables. According to the fractional-order generalized thermoelastic model, we have to construct new classifications for all mediums according to their fractional parameter where these parameters become new indicators of their power to conduct thermal energy.

5. Conclusions

In this work, we studied the impacts of strong, weak, and normal thermal conductivities in porous materials under the generalized fractional-order thermoelastic model. The resulting nondimensional equations were solved by the Fourier–Laplace transformation method and, following that, applying the eigenvalue approach. The significant impacts of the strong, normal, and weak thermal conductivities were discussed for all physical quantities. Accordingly, generalized thermoelastic fractional-order models are considered as an advancement in the study of porous elastic materials.

Author Contributions

T.S. and I.A.A. have equally contributed to the elaboration of this manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external founding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The matrix A = [ a i j ] ,   i , j = 1 8 , and a i j = 0 excepting
a 82 = r 11 i q s + 1 r 9 + s ,   a 83 = r 10 s r 9 + s ,   a 84 = q 2 + s + 1 r 9 + s ,   a 85 = s + 1   r 11 r 9 + s , a 72 = s 7 i q ,   a 73 = s 2 r 4 + q 2 + r 6 + r 5 s ,   a 74 = r 8 ,   a 75 = r 7 , a 62 = ( s 2 + q 2 ) r 1 ,   a 63 = r 2 i q r 1 ,   a 64 = r 3 i q r 1 ,   a 65 = i q ( 1 r 1 ) r 1 ,   a 51 = s 2 + r 1 q 2 ,   a 56 = i q ( 1 r 1 ) ,   a 57 = r 2 , a 58 = r 3 ,   a 15 = a 26 = a 37 = a 48 = 1 , and   M = [ u ¯ * v ¯ * φ ¯ * Θ ¯ * d u ¯ * d x d v ¯ * d x d φ ¯ * d x d Θ ¯ * d x ] T .

Appendix B

f 1 = a 58 a 85 + a 56 a 65 + a 57 a 75 + a 51 + a 62 + a 84 + a 73 , f 2 = a 56 a 65 a 84 + a 51 a 73 a 74 a 83 + a 58 a 62 a 85 a 57 a 74 a 85 + a 51 a 84 + a 57 a 75 a 84 a 56 a 64 a 85 + a 58 a 73 a 85 + a 62 a 73 + a 56 a 65 a 73 a 58 a 75 a 83 + a 57 a 62 a 75 a 56 a 63 a 75 + a 51 a 62 a 63 a 72 a 57 a 65 a 72 a 64 a 82 a 58 a 65 a 82 + a 62 a 84 + a 73 a 84 , f 3 = a 56 a 64 a 73 a 85 + a 57 a 65 a 72 a 84 + a 56 a 63 a 75 a 84 + a 57 a 62 a 74 a 85 a 56 a 65 a 73 a 84 a 57 a 62 a 75 a 84 a 58 a 63 a 75 a 82 a 58 a 65 a 72 a 83 + a 51 a 64 a 82 + a 64 a 73 a 82 + a 58 a 65 a 73 a 82 + a 51 a 74 a 83 + a 62 a 74 a 83 a 63 a 74 a 82 a 57 a 65 a 74 a 82 + a 63 a 72 a 84 a 56 a 64 a 75 a 83 a 51 a 62 a 84 + a 57 a 64 a 75 a 82 a 64 a 72 a 83 a 51 a 62 a 73 + a 58 a 63 a 72 a 85 a 57 a 64 a 72 a 85 a 58 a 62 a 73 a 85 a 51 a 73 a 84 a 62 a 73 a 84 + a 56 a 65 a 74 a 83 + a 58 a 62 a 75 a 83 a 56 a 63 a 74 a 85 + a 51 a 63 a 72 , f 4 = a 51 a 64 a 72 a 83 + a 51 a 62 a 73 a 84 a 51 a 64 a 73 a 82 + a 51 a 63 a 74 a 82 a 51 a 62 a 74 a 83 a 51 a 63 a 72 a 84 .

Appendix C

Y 4 = a 58 ξ ( ξ 6 + a 51 ( a 63 a 72 ξ 4 + ( ξ 2 a 73 ) a 62 + a 73 ξ 2 ) ξ 2 ( ( a 57 a 65 + a 63 ) a 72 + a 73 ξ 2 + a 57 a 75 ξ 2 a 62 ( a 57 a 75 ξ 2 + a 73 ) + a 56 ( a 75 a 63 + ( ξ 2 a 73 ) a 65 ) ) ) , Y 3 = ξ a 58 ( ( ( ξ 2 + a 62 ) ( ξ 2 + a 51 ) ξ 2 a 65 a 56 ) a 74 + ( a 72 ( ξ 2 a 51 ) + ξ 2 a 75 a 56 ) a 64 + a 58 ξ 2 ( ( ξ 2 a 62 ) a 75 + a 72 a 65 ) ) , Y 2 = ξ ( a 63 a 58 a 64 a 57 ) ( ( a 51 + ξ 2 ) a 74 + ξ 2 a 75 a 58 ) ξ ( a 64 ( ξ 2 a 51 ) + ξ 2 a 65 a 58 ) ( a 57 a 74 + a 58 ( ξ 2 a 73 ) ) , Y 1 = ( ξ a 58 ( a 62 ξ 2 ) ξ a 64 a 56 ) ( ξ a 58 ( a 73 ξ 2 ) ξ a 74 a 57 ) + ξ 2 ( a 57 a 64 + a 58 a 63 ) ( a 56 a 74 + a 58 a 72 ) , Y 5 = Y 1 ξ ,   Y 6 = Y 2 ξ , Y 7 = Y 3 ξ , Y 8 = Y 4 ξ ,

References

  1. Biot, M.A. Thermoelasticity and irreversible thermodynamics. J. Appl. Phys. 1956, 27, 240–253. [Google Scholar] [CrossRef]
  2. Rosencwaig, A.; Opsal, J.; Willenborg, D.L. Thin-film thickness measurements with thermal waves. Appl. Phys. Lett. 1983, 43, 166–168. [Google Scholar] [CrossRef]
  3. Biot, M.A. General solutions of the equations of elasticity and consolidation for a porous material. J. Appl. Mech. 1956, 23, 91–96. [Google Scholar]
  4. Biot, M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range. J. Acoust. Soc. Am. 1956, 28, 179–191. [Google Scholar] [CrossRef]
  5. Lord, H.W.; Shulman, Y. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  6. Green, A.E.; Naghdi, P.M. A re-examination of the basic postulates of thermomechanics. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1991, 432, 171–194. [Google Scholar]
  7. Green, A.; Naghdi, P. On undamped heat waves in an elastic solid. J. Therm. Stresses 1992, 15, 253–264. [Google Scholar] [CrossRef]
  8. Green, A.; Naghdi, P. Thermoelasticity without energy dissipation. J. Elast. 1993, 31, 189–208. [Google Scholar] [CrossRef]
  9. Youssef, H.M. Theory of fractional order generalized thermoelasticity. J. Heat Transf. 2010, 132, 061301. [Google Scholar] [CrossRef]
  10. Youssef, H.M.; Al-Lehaibi, E.A. Variational principle of fractional order generalized thermoelasticity. Appl. Math. Lett. 2010, 23, 1183–1187. [Google Scholar] [CrossRef] [Green Version]
  11. Sherief, H.H.; El-Sayed, A.M.A.; Abd El-Latief, A.M. Fractional order theory of thermoelasticity. Int. J. Solids Struct. 2010, 47, 269–275. [Google Scholar] [CrossRef]
  12. Ezzat, M.A.; El Karamany, A.S. Theory of fractional order in electro-thermoelasticity. Eur. J. Mech. A Solids 2011, 30, 491–500. [Google Scholar] [CrossRef]
  13. Ezzat, M.; El-Karamany, A.; El-Bary, A. Modeling of memory-dependent derivative in generalized thermoelasticity. Eur. Phys. J. Plus 2016, 131, 372. [Google Scholar] [CrossRef]
  14. Marin, M. Some basic theorems in elastostatics of micropolar materials with voids. J. Comput. Appl. Math. 1996, 70, 115–126. [Google Scholar] [CrossRef] [Green Version]
  15. Saeed, T.; Abbas, I.; Marin, M. A GL Model on Thermo-Elastic Interaction in a Poroelastic Material Using Finite Element Method. Symmetry 2020, 12, 488. [Google Scholar] [CrossRef] [Green Version]
  16. Ouyang, X.L.; Xu, R.N.; Jiang, P.X. Three-equation local thermal non-equilibrium model for transient heat transfer in porous media: The internal thermal conduction effect in the solid phase. Int. J. Heat Mass Transf. 2017, 115, 1113–1124. [Google Scholar] [CrossRef]
  17. Abbas, I.A. The effects of relaxation times and a moving heat source on a two-temperature generalized thermoelastic thin slim strip. Can. J. Phys. 2015, 93, 585–590. [Google Scholar] [CrossRef]
  18. Zenkour, A.M.; Abbas, I.A. Magneto-thermoelastic response of an infinite functionally graded cylinder using the finite element method. J. Vib. Control 2014, 20, 1907–1919. [Google Scholar] [CrossRef]
  19. Abbas, I.A. Nonlinear transient thermal stress analysis of thick-walled FGM cylinder with temperature-dependent material properties. Meccanica 2014, 49, 1697–1708. [Google Scholar] [CrossRef]
  20. El-Naggar, A.M.; Kishka, Z.; Abd-Alla, A.M.; Abbas, I.A.; Abo-Dahab, S.M.; Elsagheer, M. On the initial stress, magnetic field, voids and rotation effects on plane waves in generalized thermoelasticity. J. Comput. Theor. Nanosci. 2013, 10, 1408–1417. [Google Scholar] [CrossRef]
  21. Sur, A.; Mondal, S. A generalized thermoelastic problem due to nonlocal effect in presence of mode I crack. J. Therm. Stresses 2020, 43, 1277–1299. [Google Scholar] [CrossRef]
  22. Othman, M.I.; Mondal, S. Memory-dependent derivative effect on wave propagation of micropolar thermoelastic medium under pulsed laser heating with three theories. Int. J. Numer. Methods Heat Fluid Flow 2019, 30, 1025–1046. [Google Scholar] [CrossRef]
  23. Abbas, I.A.; El-Amin, M.; Salama, A. Effect of thermal dispersion on free convection in a fluid saturated porous medium. Int. J. Heat Fluid Flow 2009, 30, 229–236. [Google Scholar] [CrossRef]
  24. Hussein, E.M. Mathematical Model for Thermoelastic Porous Spherical Region Problems. Comput. Therm. Sci. Int. J. 2020, 12, 233–248. [Google Scholar] [CrossRef]
  25. Hobiny, A.; Abbas, I. Generalized thermoelastic interaction in a two-dimensional porous medium under dual phase lag model. Int. J. Numer. Methods Heat Fluid Flow 2020, 30, 4865–4881. [Google Scholar] [CrossRef]
  26. Biswas, S. Surface waves in porous nonlocal thermoelastic orthotropic medium. Acta Mech. 2020, 231, 2741–2760. [Google Scholar] [CrossRef]
  27. Carini, M.; Zampoli, V. On porous matrices with three delay times: A study in linear thermoelasticity. Mathematics 2020, 8, 371. [Google Scholar] [CrossRef] [Green Version]
  28. Abbas, I.A.; Marin, M. Analytical Solutions of a Two-Dimensional Generalized Thermoelastic Diffusions Problem Due to Laser Pulse. Iran. J. Sci. Technol. Trans. Mech. Eng. 2018, 42, 57–71. [Google Scholar] [CrossRef]
  29. Shekhar, S. Study of deformation due to thermal shock in porous thermoelastic material with reference temperature dependent properties. Mech. Based Des. Struct. Mach. 2020. [Google Scholar] [CrossRef]
  30. Itu, C.; Öchsner, A.; Vlase, S.; Marin, M.I. Improved rigidity of composite circular plates through radial ribs. Proc. Inst. Mech. Eng. Part L J. Mater. Des. Appl. 2019, 233, 1585–1593. [Google Scholar] [CrossRef]
  31. Abd-Elaziz, E.M.; Marin, M.; Othman, M.I. On the effect of Thomson and initial stress in a thermo-porous elastic solid under GN electromagnetic theory. Symmetry 2019, 11, 413. [Google Scholar] [CrossRef] [Green Version]
  32. Abbas, I.A.; Kumar, R. Deformation due to thermal source in micropolar generalized thermoelastic half-space by finite element method. J. Comput. Theor. Nanosci. 2014, 11, 185–190. [Google Scholar] [CrossRef]
  33. Sheikholeslami, M.; Ellahi, R.; Shafee, A.; Li, Z. Numerical investigation for second law analysis of ferrofluid inside a porous semi annulus: An application of entropy generation and exergy loss. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 1079–1102. [Google Scholar] [CrossRef]
  34. Ellahi, R.; Sait, S.M.; Shehzad, N.; Ayaz, Z. A hybrid investigation on numerical and analytical solutions of electro-magnetohydrodynamics flow of nanofluid through porous media with entropy generation. Int. J. Numer. Methods Heat Fluid Flow 2019, 30, 834–854. [Google Scholar] [CrossRef]
  35. Singh, B. Wave propagation in a generalized thermoelastic material with voids. Appl. Math. Comput. 2007, 189, 698–709. [Google Scholar] [CrossRef]
  36. Palani, G.; Abbas, I. Free convection MHD flow with thermal radiation from an impulsively started vertical plate. Nonlinear Anal. Model. Control 2009, 14, 73–84. [Google Scholar] [CrossRef] [Green Version]
  37. Villatoro, F.R.; Pérez, J.; Santander, J.L.G.; Borovsky, M.A.; Ratis, Y.L.; Izzheurov, E.A.; de Córdoba, P.F. Perturbation analysis of the heat transfer in porous media with small thermal conductivity. J. Math. Anal. Appl. 2011, 374, 57–70. [Google Scholar] [CrossRef] [Green Version]
  38. Abbas, I. Natural frequencies of a poroelastic hollow cylinder. Acta Mech. 2006, 186, 229–237. [Google Scholar] [CrossRef]
  39. Alzahrani, F.; Hobiny, A.; Abbas, I.; Marin, M. An Eigenvalues Approach for a Two-Dimensional Porous Medium Based Upon Weak, Normal and Strong Thermal Conductivities. Symmetry 2020, 12, 848. [Google Scholar] [CrossRef]
  40. Hobiny, A.; Alzahrani, F.; Abbas, A.; Marin, M. The Effect of fractional time derivative of bioheat model in skin tissue induced to laser irradiation. Symmetry 2020, 12, 602. [Google Scholar] [CrossRef]
  41. Das, N.C.; Lahiri, A.; Giri, R.R. Eigenvalue approach to generalized thermoelasticity. Indian J. Pure Appl. Math. 1997, 28, 1573–1594. [Google Scholar]
  42. Abbas, I.A.; Alzahrani, F.S.; Elaiw, A. A DPL model of photothermal interaction in a semiconductor material. Waves Random Complex Media 2019, 29, 328–343. [Google Scholar] [CrossRef]
  43. Saeed, T.; Abbas, I. Thermomechanical response in a two-dimension porous medium subjected to thermal loading. Int. J. Numer. Methods Heat Fluid Flow 2019, 30, 4103–4117. [Google Scholar] [CrossRef]
  44. Stehfest, H. Algorithm 368: Numerical inversion of Laplace transforms [D5]. Commun. ACM 1970, 13, 47–49. [Google Scholar] [CrossRef]
  45. Othman, M.I.; Marin, M. Effect of thermal loading due to laser pulse on thermoelastic porous medium under GN theory. Results Phys. 2017, 7, 3863–3872. [Google Scholar] [CrossRef]
Figure 1. The change in volume fraction field of void variations φ via x and y = 0.4 for strong normal and weak conductivities.
Figure 1. The change in volume fraction field of void variations φ via x and y = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g001
Figure 2. The variation of temperature Θ via x and y = 0.4 for strong normal and weak conductivities.
Figure 2. The variation of temperature Θ via x and y = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g002
Figure 3. The variations of vertical displacement v via x with y = 0.4 for strong normal and weak conductivities.
Figure 3. The variations of vertical displacement v via x with y = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g003
Figure 4. The variations of horizontal displacement u via x with y = 0.4 for strong normal and weak conductivities.
Figure 4. The variations of horizontal displacement u via x with y = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g004
Figure 5. The variation of stress σ x y via x with y = 0.4 for strong normal and weak conductivities.
Figure 5. The variation of stress σ x y via x with y = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g005
Figure 6. The variation of stress σ x x via x and y = 0.4 for strong normal and weak conductivities.
Figure 6. The variation of stress σ x x via x and y = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g006
Figure 7. The change in volume fraction field of void variations φ via y and x = 0.4 for strong normal and weak conductivities.
Figure 7. The change in volume fraction field of void variations φ via y and x = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g007
Figure 8. The variations of temperature Θ via y and x = 0.4 for strong normal and weak conductivities.
Figure 8. The variations of temperature Θ via y and x = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g008
Figure 9. The variations of vertical displacement v via y with x = 0.4 for strong normal and weak conductivities.
Figure 9. The variations of vertical displacement v via y with x = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g009
Figure 10. The variations of horizontal displacement u via y with x = 0.4 for strong normal and weak conductivities.
Figure 10. The variations of horizontal displacement u via y with x = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g010
Figure 11. The variation of stresses component σ x y versus y when x = 0.4   for strong normal and weak conductivities.
Figure 11. The variation of stresses component σ x y versus y when x = 0.4   for strong normal and weak conductivities.
Mathematics 09 00207 g011
Figure 12. The variation of stresses components σ x x versus y when x = 0.4 for strong normal and weak conductivities.
Figure 12. The variation of stresses components σ x x versus y when x = 0.4 for strong normal and weak conductivities.
Mathematics 09 00207 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saeed, T.; A. Abbas, I. The Effect of Fractional Time Derivative on Two-Dimension Porous Materials Due to Pulse Heat Flux. Mathematics 2021, 9, 207. https://doi.org/10.3390/math9030207

AMA Style

Saeed T, A. Abbas I. The Effect of Fractional Time Derivative on Two-Dimension Porous Materials Due to Pulse Heat Flux. Mathematics. 2021; 9(3):207. https://doi.org/10.3390/math9030207

Chicago/Turabian Style

Saeed, Tareq, and Ibrahim A. Abbas. 2021. "The Effect of Fractional Time Derivative on Two-Dimension Porous Materials Due to Pulse Heat Flux" Mathematics 9, no. 3: 207. https://doi.org/10.3390/math9030207

APA Style

Saeed, T., & A. Abbas, I. (2021). The Effect of Fractional Time Derivative on Two-Dimension Porous Materials Due to Pulse Heat Flux. Mathematics, 9(3), 207. https://doi.org/10.3390/math9030207

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop