Mohamed Abdelsabour Fahmy
1Department of Mathematics,Jamoum University College,Umm Al-Qura University,Alshohdaa,25371,Jamoum,Saudi Arabia
2Department of Basic Sciences,Faculty of Computers and Informatics,Suez Canal University,New Campus,Ismailia,41522,Egypt
Abstract: The main purpose of the current article is to develop a novel boundary element model for solving fractional-order nonlinear generalized porothermoelastic wave propagation problems in the context of temperaturedependent functionally graded anisotropic (FGA) structures.The system of governing equations of the considered problem is extremely very difficult or impossible to solve analytically due to nonlinearity,fractional order diffusion and strongly anisotropic mechanical and physical properties of considered porous structures.Therefore, an efficient boundary element method (BEM)has been proposed to overcome this difficulty,where,the nonlinear terms were treated using the Kirchhoff transformation and the domain integrals were treated using the Cartesian transformation method (CTM).The generalized modified shift-splitting(GMSS)iteration method was used to solve the linear systems resulting from BEM,also,GMSS reduces the iterations number and CPU execution time of computations.The numerical findings show the effects of fractional order parameter,anisotropy and functionally graded material on the nonlinear porothermoelastic stress waves.The numerical outcomes are in very good agreement with those from existing literature and demonstrate the validity and reliability of the proposed methodology.
Keywords: Boundary element method; fractional-order; nonlinear generalized porothermoelasticity; wave propagation; functionally graded anisotropic structures; Cartesian transformation method
The fractional order calculus (FOC) is the branch of mathematical analysis dealing with non-integer order calculus and its applications.The essential viewpoints are sketched out for fractional calculus theory in [1] and for fractional calculus applications in [2–6].FOC is nowadays extremely popular due to its applications in different fields such as diffusion equation, quantum mechanics, nanotechnology, solid mechanics, continuum mechanics, biochemistry, wave propagation theory, polymers, robotics and control theory, finance and control theory, electrochemistry,electrical engineering, fluid dynamics, signal and image processing, biophysics, electric circuits,viscoelasticity, electronics, field theory, group theory, etc.
Several researchers have contributed to the background of fractional calculus [7–9].Recently,Yu et al.[10] introduced new definitions of fractional derivative in the context of thermoelasticity.Research on generalized thermo-elasticity theories [11] has attracted much attention from many scientists, among which are research in magneto-thermoelasticity [12], viscothermoelasticity [13,14] and micropolar-thermoelasticity [15,16].
Because of computational complexity in solving complex fractional thermoelasticity problems not having any general analytical solution, computational techniques should be used to solve such problems.Among these computational techniques are the boundary element method (BEM) that has been used for magneto thermoviscoelasticity [17,18], computerized engineering models [19,20],and design sensitivity and optimization [21,22] and nonlinear problems [23–26].The BEM presents an attractive alternative numerical method to the domain methods for the investigation of thermoelastic wave propagation problems, like finite element method (FEM) [27–29] and finite volume method (FVM) [30–32].The main feature of BEM over the domain type methods is that it requires boundary-only discretization of the domain under consideration.This feature has significant importance for solving complex thermoelastic problems with fewer elements, and requires very little computational cost, much less preparation of input data, and therefore easier to use.
In the present paper, we introduce a new boundary element model for solving fractionalorder nonlinear generalized porothermoelastic wave propagation problems.The nonlinear terms are treated using the Kirchhoff transformation.The domain integrals were treated using the Cartesian transformation method.In the proposed BEM technique, the temperature and displacement distributions were calculated using a partitioned semi-implicit predictor–corrector coupling algorithm.Then, we can obtain the propagation of porothermoelastic stress waves in temperaturedependent FGA structures.Numerical results demonstrate the validity, accuracy and efficiency of our proposed model and technique.
The geometry of the considered problem is depicted in Fig.1.The governing equations for fractional-order nonlinear generalized porothermoelastic wave propagation problems in the context of FGA structures can be written as [33]
whereσijis the mechanical stress tensor,ρis the bulk density,ρFis the fluid density,Fiis the bulk body forces,φis the porosity,uiis the solid displacement andviis the fluid–solid displacement.
whereζis the variation of the fluid volume per unit reference volume,qis the instantaneous flux andis the source term.
The fractional nonlinear heat conduction equation can be expressed in non-dimensionless form as
in which
in which the heat source functionh(X,T,t)can be written as
whereTis the temperature,λis the thermal conductivity,Cijklis the constant elastic moduli,Ais the Biot’s effective stress coefficient,pis the fluid pressure,βijis the stress-temperature coefficients,is the permeability,T0is the reference temperature, Å is a unified parameter that introduces all generalized thermoelasticity theories into a unified system of equations,QandRare solid–fluid coupling parameters,τ0,τ1, andτ2are relaxation times,ρ0=ηφρFandηis the shape factor.
Figure 1:Geometry of the considered problem
According to finite difference scheme of Caputo at times(f+1)ΔτandfΔτ, we obtain [34]
where
On the basis of Eq.(7), the fractional heat conduction Eq.(3) can be expressed as
whereJ=1,2,...,Fandf=0,1,2,...,F.
By using the transformation of KirchhoffEq.(3) can be written as [35]
The decomposition of the right-hand side of (10) into linear and nonlinear sections, yields
The nonlinear section can be written as
Based on [24], we can write (11) into the following form
where
Now, by using the fundamental solution of (9), we can write the boundary integral equation corresponding to (13) as [36]
By substituting ofΘ(P,tn+1)=2ΘP,tn+(1/2)−Θ(P,tn)in (15), we get
Now, the domain integrals in Eq.(16) can be computed using CTM.Thus, the unknown boundary values can be calculated from the following system
whereΘΓandQΓareM′dimension vectors, andHandGareM′×M′dimension matrices.
Thus, the unknown internal values can be calculated from the following system
If we have assumed that the time step size is constant, then,H,G,andcan be computed at all time steps.Also,F,FNl,, andbe computed at all time steps using CTM.
Now, we are considering the following regular domain integral [37,38]
Based on Khosravifard et al.[39], we can write the domain integral (19) as follows
where
By applying the composite Gaussian quadrature method to (19), we obtain
which can be written as
By implementing the radial point interpolation method (RPIM) [40], we can write
where M equals the summation of boundary nodesM′and internal pointsM′′.
Based on [40], the functionp(x1,x2)may be described as
To build the RPIM shape functions, we applied the following Gaussian radial basis function
whereαiandbjare unknown coefficients which can be computed from the following system
and the followingconditions
By using Eqs.(27) and (28), we can expressαiandbjas
Thus, based on [40], and using (29), we can write Eq.(25) in the following form
Thus, we have
which can be written as
where p contains boundary and internalpvalues.
We now consider the following domain integrals that appear in the integral Eq.(16)whereEi(x)=−0.57721566+−ln(x)
According to [25], the weakly singular in (33) can be regularized to obtain
where
and
Also, the domain integral in (34) can be regularized to obtain
where
Jack Zipes offers a few theories behind the name s selection in The Great Fairy Tale Tradition: The derivation of this name is not certain. It was known through legend that there was a fool in Alexandria who was referred to as Carabas by the inhabitants of the city, and they mocked him by treating him as if he were a king. The Turkish word Carabag designates a beautiful place in the mountains where the sultans and princes would spend the summer months. Perrault might have come upon this term in the Dictionaire oriental, edited by Barthelemy d Herbelot during this time. (Zipes 2001, 398). Return to place in story.
and
Hence, from (18) we get
where a is an unknown matrix, while X and b are known matrices.
Based on the weighted residual technique, we can write Eqs.(1) and (2) as follows
where
On using integration by parts for the first term of Eqs.(42) and (43), we get
Based on Fahmy [24], elastic stress can be expressed as
where
Now, we consider the following definitions
Substituting above definitions into (47), we get
which after integration can be written as
where
Now, we can write (50) as
which can be expressed as follows
where the vectors Q, P, i, and j are displacements, tractions, pore pressure, and pore pressure gradients, respectively.
Substituting the boundary conditions into (54), we obtain the following system of equations
in which A represents unknown matrix, while X and B represent known matrices.
According to Breuer et al.[41], a robust and efficient partitioned semi-implicit predictor–corrector coupling algorithm was implemented with GMSS [42] for solving the resulting linear Eqs.(41) and (54) arising from the boundary element discretization, where poro-thermo-elastic coupling is considered instead of fluid-structure-interaction coupling.
The proposed BEM technique which is based on the coupling algorithm [41], should be applied to a wide variety of fractional-order nonlinear porothermoelastic wave propagation problems.
In the present paper, we considered the temperature-dependent properties of anisotropic porous copper material, where the specific heat and density are tabulated in Tab.1 [43].
Table 1:Temperature-dependent specific heat and density of porous copper material
The thermal conductivity is given by
The domain boundary of the current problem has been discretized into 42 boundary elements and 68 internal points as depicted in Fig.2.
Figure 2:Boundary element model of the considered problem
Figs.3–5 illustrate the propagation of nonlinear thermal stress wavesσ11,σ12, andσ22for different values(a=0.4,0.7 and 1.0)of the fractional order parameter (FOP).It can be seen from these figures that the FOP has a great influence on the nonlinear thermal stress waves of FGA porous structures.
Figure 3:Propagation of the nonlinear thermal stress σ11 waves with time t for different values of the fractional-order parameter
Figure 4:Propagation of the nonlinear thermal stress σ12 waves with time t for different values of the fractional-order parameter
Figure 5:Propagation of the nonlinear thermal stress σ22 waves with time t for different values of the fractional-order parameter
According to the relationship of elastic constants for anisotropic, isotropic, and orthotropic materials [44].We therefore considered these three materials in the current study.
Figs.6–8 show the propagation of nonlinear thermal stress wavesσ11,σ12, andσ22for anisotropic, isotropic and orthotropic functionally graded porous structures.It can be shown from these figures that the effects of anisotropy are very pronounced.
Figure 6:Propagation of the nonlinear thermal stress σ11 waves with time t for isotropic,orthotropic and anisotropic porous materials
Figure 7:Propagation of the nonlinear thermal stress σ12 waves with time t for isotropic,orthotropic and anisotropic porous materials
Figure 8:Propagation of the nonlinear thermal stress σ22 waves with time t for isotropic,orthotropic and anisotropic porous materials
Figs.9–11 display the propagation of nonlinear thermal stress wavesσ11,σ12, andσ22for homogeneous(m= 0)and functionally graded(m=0.4 and 0.7)porous structures.It can be shown from these figures that the effect of functionally graded material is very pronounced.
Figure 9:Propagation of the nonlinear thermal stress σ11 waves with time t for homogeneous and functionally graded porous materials
Figure 10:Propagation of the nonlinear thermal stress σ12 waves with time t for homogeneous and functionally graded porous materials
Figure 11:Propagation of the nonlinear thermal stress σ22 waves with time t for homogeneous and functionally graded porous materials
The effectiveness of our proposed approach has been established through the use of the GMSS which doesn’t need the entire matrix to be stored in the memory and converges quickly without the need for complicated calculations.During our treatment of the considered problem,we implemented GMSS, Uzawa-HSS, and regularized iteration methods [45].Tab.2 displays the number of iterations (IT), processor time (CPU), relative residual (RES), and error (ERR) of the considered methods computed for different fractional order values.It can be noted from Tab.2 that the GMSS needs the lowest IT and CPU times, which means that GMSS method has better performance than Uzawa-HSS and regularized methods.
For comparison purposes with other methods, we only considered the one-dimensional special case.Therefore, the time distribution results of the nonlinear thermal stressσ11are plotted in Fig.12 for the proposed BEM and compared with the FDM results obtained by Awrejcewicz et al.[46] and FEM results obtained by Shakeriaski et al.[47], it can be shown from Fig.12 that the BEM outcomes are in very good agreement with the FDM and FEM outcomes.Thus, the validity, accuracy, and usefulness of the proposed BEM have been demonstrated.
Table 2:Numerical results for the tested iteration methods
Figure 12:Propagation of the nonlinear thermal stress σ11 waves with time t for a special case and different methods
The main objective of the current paper is to develop a new boundary element model for solving fractional-order nonlinear generalized porothermoelastic wave propagation problems in FGA structures, which are difficult or impossible to solve analytically.Therefore, an efficient numerical procedure based on BEM has been proposed to overcome this challenge.The Kirchhoff transformation is first used to treat the nonlinear terms.Then, the Cartesian transformation method (CTM) has been applied to transform the domain integration into boundary integration,As a result, the computational complexity of integration and CPU computing time are significantly reduced.The memory requirements and Processing time are also reduced by applying the GMSS method which does not need that the entire matrix is stored in the memory, and it is rapidly converging without the need for complicated calculations.The numerical outcomes are presented graphically to show the effects of fractional parameter, anisotropy, and functionally graded material on the nonlinear thermal stress waves.The numerical outcomes also show very good agreement with the earlier work in the literature as a special case.These outcomes also confirm the validity, accuracy, and effectiveness of the proposed methodology.
Funding Statement:The author received no specific funding for this study.
Conficts of Interest:The author declares that he has no conflicts of interest to report regarding the present study.
Computers Materials&Continua2021年7期