
Citation: Oleg S. Sukharev. Economic crisis as a consequence COVID-19 virus attack: risk and damage assessment[J]. Quantitative Finance and Economics, 2020, 4(2): 274-293. doi: 10.3934/QFE.2020013
[1] | Jue Er Amanda Lee, Linda May Parsons, Leonie M. Quinn . MYC function and regulation in flies: how Drosophila has enlightened MYC cancer biology. AIMS Genetics, 2014, 1(1): 81-98. doi: 10.3934/genet.2014.1.81 |
[2] | Tú Nguyen-Dumont, Jenna Stewart, Ingrid Winship, Melissa C. Southey . Rare genetic variants: making the connection with breast cancer susceptibility. AIMS Genetics, 2015, 2(4): 281-292. doi: 10.3934/genet.2015.4.281 |
[3] | Ke-Sheng Wang, Xuefeng Liu, Daniel Owusu, Yue Pan, Changchun Xie . Polymorphisms in the ANKS1B gene are associated with cancer, obesity and type 2 diabetes. AIMS Genetics, 2015, 2(3): 192-203. doi: 10.3934/genet.2015.3.192 |
[4] | Michael J. Murray . Drosophila models of metastasis. AIMS Genetics, 2015, 2(1): 25-53. doi: 10.3934/genet.2015.1.25 |
[5] | Mohammad Hashemi, Fatemeh Bizhani, Hiva Danesh, Behzad Narouie, Mehdi Sotoudeh, Mohammad Hadi Radfar, Mehdi Honarkar Ramezani, Gholamreza Bahari, Mohsen Taheri, Saeid Ghavami . MiR-608 rs4919510 C>G polymorphism increased the risk of bladder cancer in an Iranian population. AIMS Genetics, 2016, 3(4): 212-218. doi: 10.3934/genet.2016.4.212 |
[6] | Carina Ladeira, Lenka Smajdova . The use of genotoxicity biomarkers in molecular epidemiology: applications in environmental, occupational and dietary studies. AIMS Genetics, 2017, 4(3): 166-191. doi: 10.3934/genet.2017.3.166 |
[7] | Noel Pabalan, Neetu Singh, Eloisa Singian, Caio Parente Barbosa, Bianca Bianco, Hamdi Jarjanazi . Associations of CYP1A1 gene polymorphisms and risk of breast cancer in Indian women: a meta-analysis. AIMS Genetics, 2015, 2(4): 250-262. doi: 10.3934/genet.2015.4.250 |
[8] | Helena E. Richardson . Drosophila models of cancer. AIMS Genetics, 2015, 2(1): 97-103. doi: 10.3934/genet.2015.1.97 |
[9] | Tin Tin Su . Non-autonomous consequences of cell death and other perks of being metazoan. AIMS Genetics, 2015, 2(1): 54-69. doi: 10.3934/genet.2015.1.54 |
[10] | Valeriia Mikhailova, Valeriia Gulaia, Vladlena Tiasto, Stanislav Rybtsov, Margarita Yatsunskaya, Alexander Kagansky . Towards an advanced cell-based in vitro glioma model system. AIMS Genetics, 2018, 5(2): 91-112. doi: 10.3934/genet.2018.2.91 |
Thermoelasticity, which examines the interaction between thermal and mechanical properties in materials, has long been a cornerstone in the fields of materials science and engineering. Traditional thermoelastic theories, largely founded on Fourier's law of heat conduction, have laid the groundwork for the analysis of thermal stress and energy dissipation [1,2]. Nevertheless, these classical models often fall short in accurately depicting the behavior of contemporary materials, particularly those with microstructural anomalies such as voids, cavities, or granular structures. These characteristics, frequently present in advanced composites, porous materials, and microelectromechanical systems (MEMS), create localized thermal gradients, contribute to wave scattering, and induce size-dependent phenomena, thereby challenging the validity of classical theoretical models [3].
To overcome these shortcomings, researchers have formulated generalized thermoelastic models, including the two-phase-lag (TPL) model [4,5,6], which accounts for time delays in the correlation between heat flux and temperature gradients. Moreover, Lord and Shulman [7] put forward a theory based on a revised heat conduction law that includes the rate of heat flux. At the same time, Green and Lindsay [8] established a theory that factors in the rate of temperature change as one of the constitutive variables. These generalized theories suggest that heat propagation should be understood as a wave-like phenomenon instead of solely a diffusion process, a concept frequently discussed in academic literature as the "second sound effect".
Although these cutting-edge models improve upon traditional methods, they frequently face challenges when addressing materials with nonlocal interactions, where both thermal and mechanical responses are shaped by spatial and temporal histories [9]. Furthermore, materials that display memory-dependent behavior require more sophisticated modeling approaches to accurately represent their intricate dynamics. These issues underscore the need for continued advancements in thermoelasticity to provide reliable predictions relevant to real-world applications. Nonlocal elasticity theories [10,11,12,13,14] have emerged as a notable breakthrough in materials science, particularly for analyzing the behavior of materials at the nanoscale. These theories introduce an internal length scale, which is essential for capturing size-dependent effects that classical elasticity models often fail to account for.
In recent years, the concept of spatiotemporal nonlocality has gained significant attention in physics and engineering, particularly in the study of nanoscale materials. Traditional local theories, which assume that a material's behavior at a given point is determined solely by its immediate surroundings, often fall short in accurately describing phenomena where size effects play a critical role [15]. Spatiotemporal nonlocality addresses these limitations by accounting for interactions that occur across both space and time, providing a more comprehensive framework to analyze complex material behaviors [16,17].
This approach is especially relevant for modeling phenomena such as wave propagation, thermal conduction, and mechanical deformation in nanostructured materials. By incorporating spatiotemporal nonlocality, researchers can better understand how microstructural features and boundary conditions significantly influence the overall behavior of materials at smaller scales [18]. This methodology captures the intricate interplay between distant regions and historical effects, which are often overlooked in classical models but are crucial for nanoscale systems [19,20].
The integration of spatiotemporal nonlocality into modeling frameworks allows for more accurate predictions of material responses under various loading conditions and environmental influences. This advancement not only enhances fundamental understanding but also opens new avenues for designing and optimizing advanced materials and structures [21]. As a result, it facilitates breakthroughs in diverse fields, including nanotechnology, aerospace engineering, and biomedical applications, among others, enabling the development of innovative technologies and solutions tailored to modern challenges.
Fractional derivatives generalize the concept of classical differentiation to non-integer orders, offering a robust mathematical framework for modeling processes that involve memory and hereditary effects. Among the various definitions of fractional derivatives, the Caputo derivative is one of the most commonly used, especially valued for its effectiveness in physical and engineering applications [22]. The emergence of non-singular fractional derivatives has attracted considerable interest in recent times, especially with the advent of Caputo–Fabrizio (CF) [23,24] and Atangana–Baleanu (AB) fractional derivatives [25,26]. These groundbreaking methods provide versatile and powerful tools for accurately modeling intricate physical phenomena. The CF derivative, introduced by Caputo and Fabrizio, operates at a fractional order and is designed to address the limitations of traditional fractional derivatives that often involve singular kernels [23,24]. The AB fractional derivative [25,26] is based on the Mittag-Leffler kernel, which allows for non-singular behavior.
The introduction of tempered fractional derivatives (TF derivatives) marks an advancement in fractional calculus (FC), expanding traditional concepts to include exponential factors. This enhancement enables a more refined modeling of diverse physical phenomena, especially in systems where standard fractional derivatives may fall short in accurately capturing dynamic processes [27]. TF derivatives and integrals are formulated by incorporating an exponential factor into fractional derivatives and integrals. This adjustment brings about a tempering effect that modifies the behavior of these derivatives, making them applicable to a broader spectrum of scenarios [28]. A key application of TF derivatives lies in tempered fractional diffusion systems. In these settings, the conventional second-order derivative concerning spatial coordinates is substituted with a TF derivative, enabling the modeling of diffusion processes that showcase nonlocal and memory-dependent behaviors [29]. The parameters associated with irregular walk models featuring exponentially tempered power-law barrier distributions can be effectively characterized using TF derivatives, thereby improving the models' ability to capture intricate diffusion patterns. TF time derivatives, derived from the tempered power law for future instances, have demonstrated their usefulness across various scientific domains, including geophysics [30]. They support the modeling of phenomena such as wave propagation, heat transfer, and material behavior in atypical conditions.
In this paper, a novel model is introduced to address a critical gap in the existing literature on nonlocal thermal conduction in porous thermoelastic materials, offering both innovative improvements and substantial contributions to the field. Unlike previous approaches, the proposed framework integrates the dual-phase lag (TPL) theory with both spatial and temporal nonlocality, employing Caputo-tempered fractional derivatives to capture behavior that is simultaneously memory-dependent and reflective of microscopic scale effects. These advanced elements form a unified methodology for modeling thermomechanical phenomena in porous thermoelastic materials, overcoming challenges in accurately predicting energy dissipation, wave propagation, and thermomechanical coupling.
By effectively addressing the complexities of memory effects, the model provides predictions of material behavior that closely mirror real-world observations, thereby bridging theoretical concepts with practical applications. Furthermore, it strikes a critical balance between mathematical simplicity and physical accuracy, marking a significant advancement in the study and application of porous thermoelastic materials. This achievement not only deepens understanding but also enhances the design and optimization of engineering systems, facilitating a wide range of applications and paving the way for future research in this pivotal area.
The remainder of this paper is organized as follows: Section 2 presents the mathematical formulation of the fractional Caputo-tempered two-phase-lag (FCT-TPL) model. In Section 3, direct application of the proposed model is provided, along with the corresponding problem formulation. Section 4 describes the mathematical techniques utilized to solve the governing equations, while Section 5 introduces a numerical approach for computing the inverse transformations of the physical fields. Subsequently, Section 6 delivers numerical results, accompanied by insights and discussions regarding the impact of the fractional order, the smoothing parameter, and the specific geometric configurations. Finally, Section 7 concludes the paper by summarizing the key findings, exploring the practical implications of the results, and identifying potential directions for future research in this field.
Understanding the behavior of porous materials containing voids under thermal and mechanical loading is crucial and necessitates the use of fractional nonlocal thermoelastic models. These models provide a more accurate depiction of the complex responses exhibited by such materials, which cannot be captured by classical elasticity theories. By incorporating nonlocal effects and fractional derivatives, these models offer a comprehensive analysis of the interactions between thermal, mechanical, and void-related phenomena [31].
The governing equations for the linear homogeneous two-phase thermoelastic model with voids, which considers the displacement field ui, the void volume fraction ϕ, and the temperature variation field θ in the absence of external forces, are expressed as follows [32]:
g=ξ1θ−ξ2ϕ−𝓋e | (1) |
where ξ1 is a constitutive parameter related to thermal response, ξ2 is the coupling parameter between voids and thermal effects, g denotes the void-related body force, e=εkk refers to cubic dilation, and 𝓋 is the mechanical response parameter due to voids [33]:
hi=pϕ,i, | (2) |
where p is a constitutive parameter related to the volume fraction and hi are the components of the equilibrated stress tensor [34]:
σij=𝓋ϕδij−βθδij+λεkkδij+2μεij, | (3) |
where β=(3λ+2μ)αt connects thermal and elastic properties, λ and μ are the Lamé parameters, αt is the linear thermal expansion coefficient, and δij is the Kronecker delta.
εij=12(ui,j+uj,i), | (4) |
where εij are the strain tensor components.
Nonlocal elasticity, which considers the influence of surrounding material points on a material point's behavior, is particularly relevant for nanostructured materials where size effects are significant. The fundamental equations of nonlocal elasticity include spatial integrals that account for contributions from surrounding points, providing a more comprehensive understanding of material behavior at the nanoscale.
In this section, we establish the constitutive relations and field equations for nonlocal thermoelastic materials with voids based on fractional derivative heat transfer. The nonlocal force stress tensor Tkl(P), equilibrated stress components Hi(P), and equilibrated body force G(P) are defined as [35,36]:
Tkl(P)=∫σkl(P′)ℏ(|P−P′|)dΩ(P), | (5) |
Hi(P)=∫hi(P′)ℏ(|P−P′|)dΩ(P′), | (6) |
G(P)=∫g(r′)ℏ(|P−P′|)dΩ(P′), | (7) |
where ℏ(|P−P′|) is the kernel function that incorporates nonlocal influences, and |P−P′| is the Euclidean distance between points P and P′.
Inspired by Eringen's nonlocal elastic model and the Boltzmann superposition integral model, stress and strain are treated as convolution functions of time and space. This means the stress and strain at a reference point depend on historical data and the state of the material at all points within the domain. The mathematical representation of this concept is as follows [20,21]:
Tkl(P,t)=∫t−∞∫Ωℏ(|P−P′|,t−t′)σkl(P′,t′)dΩ(P′)dt′. | (8) |
The kernel function ℏ(|P−P′|,t−t′) determines the nature and extent of nonlocal interactions. It typically exhibits properties ensuring that the influence of distant points diminishes with increasing distance and time.
In the same vein, the nonlocal equilibrated stress components Hi(P) and the equilibrated body force G(P) can be derived as follows [37]:
Hi(P,t)=∫ℏ(|P−P′|,t−t′)hi(P′,t′)dΩ(P′)dt′, | (9) |
G(P,t)=∫ℏ(|r−r′|,t−t′)g(P′,t′)dΩ(P′)dt′. | (10) |
To further develop the nonlocal elasticity theory, we suggest employing the Klein-Gordon (KG) operator as a means to capture both spatial and temporal nonlocalities within the framework of nonlocal elasticity. The KG differential operator, denoted as Δ, is defined in the following manner [21,38]:
Δ=1−η2∇2+ζ2∂2∂t2. | (11) |
Here, η represents the internal length scale parameter that is linked to spatial nonlocality, while ζ serves as the characteristic time scale that captures temporal nonlocality. Additionally, ∇2 denotes the Laplacian operator, which addresses spatial variations within the model.
When the KG operator Δ is applied to the nonlocal stress tensor τkl, the equilibrated stress components Hi, and the equilibrated body force G, we arrive at the governing equations for isotropic materials within the framework of KG-type nonlocal elasticity [39,40]:
(1−η2∇2+ζ2∂2∂t2)Tkl=σkl=(𝓋ϕ−βθ+λεkk)δij+2μεij, | (12) |
(1−η2∇2+ζ2∂2∂t2)G=g=−ξ2ϕ+ξ1θ−𝓋e, | (13) |
(1−η2∇2+ζ2∂2∂t2)Hi=hi=α1ϕ,i. | (14) |
According to Newton's second law, the fundamental equation of motion can be expressed as:
τij,j+Fi=ρ∂2ui∂t2. | (15) |
By inserting Eq (12) into Eq (15), the equation of motion transforms into:
μui,jj+(λ+μ)uj,ij−βθ,i+𝓋ϕ,i=ρ(1−η2∇2+ζ2∂2∂t2)∂2ui∂t2. | (16) |
In the absence of any external equilibrated body force, the volume fraction field equation is represented as [41]:
G=ρχ1¨ϕ−Hi,i, | (17) |
where χ1 refers to the equilibrated inertia. By incorporating G from Eq (13) and Hi from Eq (14), the volume fraction field equation is modified to [42]:
pϕ,ii−ξ2ϕ+ξ1θ−𝓋e=χ1ρ(1−η2∇2+ζ2∂2∂t2)¨φ. | (18) |
The equation that represents the coupling between entropy (S), temperature (θ), and voids (ϕ) is expressed as [43]:
ρT0S=ξ1T0ϕ+ρCEθ+βe, | (19) |
where CE represents specific heat at constant strain. The energy equation defines the regulation of heat balance by including factors like internal heat generation (S0) and the components of heat flux (qi) as follows [44]:
ρT0˙S=−qi,i+S0. | (20) |
The TPL thermoelastic model builds upon Fourier's law by adding phase lags for heat flux (τq) and the temperature gradient (τθ), which can be expressed as [4,5,6]:
(1+τq∂∂t)qi=−K(1+τθ∂∂t)θ,i, | (21) |
where K denotes the thermal conductivity.
To increase accuracy, fractional derivatives are incorporated into the TPL model, allowing for the representation of memory effects and interactions that are nonlocal in nature. The adjusted DPL equation that includes fractional derivatives is expressed as [45]:
(1+ταqDαt)qi=−K(1+ταθDαt)θ,i. | (22) |
In this context, Dαt denotes a fractional derivative of order α that falls within the interval (0,1), extending the notion of differentiation to include effects related to memory and heredity.
Various formulations of fractional derivatives are employed based on the specific requirements of each application. In this study, we will focus solely on two key examples of fractional operators. The Caputo fractional derivative of order α (0<α<1) for Y(t) is defined as [46]:
C0D(α)0Y(t)=1Γ(1−α)t∫01(t−𝓈)αdd𝓈Y(𝓈)d𝓈. | (23) |
The Caputo-tempered fractional derivative of order α (0<α<1) with tempering parameter σ>0 is defined as [47,48]:
CT0Dα,σtY(t)=e−σtΓ(1−α)t∫0dd𝓈(eσ𝓈Y(𝓈))(t−𝓈)−αd𝓈. | (24) |
The Caputo-tempered fractional derivative incorporates a tempering parameter (σ) to adjust long-range memory effects. When σ is set to 0, it simplifies to the Caputo derivative. If σ is greater than 0, it adds an exponential tempering effect, which is perfect for systems that display declining memory characteristics. The CT derivative proves especially beneficial for materials that demonstrate both short-term and long-term memory behaviors.
Combining the fractional TPL model with the integrated thermal-mechanical framework allows us to derive the fractional heat conduction equation as follows:
(1+ταqDαt)[ρCE∂θ∂t+βT0∂εkk∂t+ξ1T0∂ϕ∂t−∂Q∂t]=K(1+ταθDαt)θ,ii. | (25) |
This study focuses on a nonlocal, homogeneous, isotropic thermoelastic medium described as infinite, containing a spherical cavity of radius R, where R<r and r tends to infinity, filled with empty pores (see Figure 1). At the start, the surface of the cavity is held at a constant temperature T0, with neither stress nor deformation present. The cavity's center is designated as the origin in the spherical polar coordinate system (r,Θ,Φ). In this coordinate setup, the body exhibits symmetrical characteristics, and the displacement components are simply defined as ur=u(r,t), while both uΘ and uΦ are zero. This expression arises due to the spherical symmetry of the problem. This configuration is based on the assumption that there are no body forces, external heat sources, or outside equilibrium forces acting on the system.
The components of the strain tensor can be expressed as follows:
εrr=∂u∂r,εΘΘ=εΦΦ=ur,εrΦ=εrΘ=εΘΦ=0. | (26) |
Based on this arrangement, the cubic dilatation e is characterized as:
e=εkk=∂u∂r+2ur=1r2∂(r2u)∂r. | (27) |
The nonlocal constitutive relationships (12) in the spherical polar coordinate system can be expressed as:
(1−η2∇2+ζ2∂2∂t2)Trr=σrr=2μ∂u∂r+λe+𝓋ϕ−βθ,(1−η2∇2+ζ2∂2∂t2)TΘΘ=σΦΦ=2μur+λe+𝓋ϕ−βθ. | (28) |
Here, in spherical coordinates, the operator ∇2 is defined as:
∇2=∂2∂r2+2r∂∂r. | (29) |
Consequently, the motion equation in spherical coordinates can be defined as:
Trr,r+2r(Trr−Tψψ)=ρ∂2u∂t2. | (30) |
Moreover, the volume fraction field equation (18) linked to the voids can be rephrased as:
p∇2ϕ−𝓋e−ξ2ϕ+ξ1θ=ρχ1(1−η2∇2+ζ2∂2∂t2)∂2ϕ∂t2. | (31) |
By substituting Eq (28) into Eq (30), we derive:
(λ+2μ)(∂2u∂r2+1r∂u∂r−2ur2)+(𝓋∂ϕ∂r−β∂θ∂r)=ρ(1−η2∇2+ζ2∂2∂t2)∂2u∂t2. | (32) |
Taking the divergence of Eq (32) results in:
(λ+2μ)∇2e+(𝓋∇2ϕ−β∇2θ)=ρ(1−η2∇2+ζ2∂2∂t2)∂2e∂t2. | (33) |
The revised fractional TPL heat transfer equation can be articulated as:
K(1+ταθDαt)∇2θ=(1+ταqDαt)[ρCE∂θ∂t+ξ1T0∂ϕ∂t+βT0∂e∂t]. | (34) |
Assuming the medium starts from rest, the following initial conditions are set:
ϕ(r,0)=0,∂ϕ(r,0)∂t=0,θ(x,0)=0,∂θ(r,0)∂t=0,τij(r,0)=0,∂τij(r,0)∂t=0,u(r,0)=0,∂u(r,0)∂t=0. | (35) |
It is assumed that the surface of the medium is free of traction, with no voids crossing the boundary. Thus, the following mechanical boundary conditions are established:
Trr(x,t)=0,∂ϕ(r,t)∂r=0 at r=R. | (36) |
The surface of the cavity at r=R is expected to experience a pulsed heat flux. Therefore, the thermal boundary condition can be formulated as [49]:
−K∂θ(r,t)∂r=q0t216t2p exp (−ttp) at r=R, | (37) |
where q0 indicates the peak intensity of the applied heat flux, and tp serves as a time parameter designed to regulate the rise and fall rates of the pulse.
Utilizing the modified Fourier's law (23), the thermal condition can be represented as:
−K(1+ταθDαt)∂θ(r,t)∂r=q0(1+ταqDαt)[t216t2pe(−ttp)] at r=R. | (38) |
To create a non-dimensional formulation for the thermoelastic porous half-space, we introduce non-dimensional variables that simplify the equations and highlight the crucial parameters influencing the system's behavior, particularly how voids affect thermal and mechanical interactions. The following non-dimensional variables will be defined based on the characteristic scales of the problem:
(r′,u′,ℓ′)=υ1η1(r,u,ℓ),(t′,τ′q,τ′θ,τ′)=υ21η1(t,τq,τθ,τ),φ′=χ1υ21η21φ,τ′ij=τijυ21ρ,θ′=βθυ21ρ,η1=ρCEK,p′=pυ21χ1ρ,υ21=λ+2μρ. | (39) |
The governing Eqs (29), (31), (33), and (34) can thus be expressed in non-dimensional form as follows:
(1−η2∇2+ζ2∂2∂t2)Trr=Υ1∂u∂r+Υ2e+R1ϕ−θ,(1−η2∇2+ζ2∂2∂t2)TΘΘ=Υ1ur+Υ2e+R1ϕ−θ, | (40) |
∇2e+R1∇2ϕ−∇2θ=(1−η2∇2+ζ2∂2∂t2)∂2e∂t2, | (41) |
p∇2ϕ−R5e−R4ϕ+R6θ=(1−η2∇2+ζ2∂2∂t2)∂2ϕ∂t2, | (42) |
(1+ταθDαt)∇2θ=(1+ταqDαt)[R2∂e∂t+R3∂ϕ∂t+∂θ∂t], | (43) |
where
R1=𝓋υ41η21χ1ρ,R2=T0β2υ21η1,R3=ξ1T0βυ41η31χ1ρK,R4=ξ2υ41η21χ1ρ,R5=𝓋υ21ρ,R6=mβ,Υ1=2μλ+2μ,Υ2=λλ+2μ. | (44) |
The solution method involves converting the system of coupled partial differential equations (PDEs) into ordinary differential equations (ODEs) by utilizing the Laplace transform. The Laplace transform for a function ℊ(x,t), which varies with time t, is defined as:
−ℊ(x,s)=L[ℊ(x,t)]=∫∞0e−stℊ(x,t)dt. | (45) |
Applying the Laplace transform to the governing Eqs (40) through (43), while incorporating the initial conditions specified in Eq (35), leads to the derivation of the following transformed equations:
(1−λ0∇2)−Trr=Q1d−udr+Q0−ur+Q3−ϕ−Q2−θ,(1−λ0∇2)−TΘΘ=Q0d−udr+Q1−ur+Q3−ϕ−Q2−θ, | (46) |
(∇2−Q4)−e=Q5∇2−θ−Q6∇2ϕ, | (47) |
(d2dx2−Q0)−θ=Q7−e+ϖ0R3−ϕ, | (48) |
(∇2−Q9)−ϕ=Q8−e−Q10−θ, | (49) |
where
Q0=Υ21+ζ2s2,Q1=Υ11+ζ2s2,Q2=11+ζ2s2,Q3=R11+ζ2s2,Q4=s2(1+ζ2s2)1+η2s2,Q5=11+η2s2,Q6=R11+η2s2,Q7=ϖ0δ2,ϖ0=s(1+τqF(s))(1+τθF(s)),Q8=R5p+η2s2,Q9=s2(1+ζ2s2)+R4p+η2s2,Q10=R6p+η2s2,λ0=η21+ζ2s2. | (50) |
The transformed fractional operator F(s) is defined based on the type of fractional derivative used:
M(s)={sαfor C fractional operotor,(s+σ)αfor CT fractional operotor. | (51) |
By combining Eqs (47) to (49), we obtain the following sixth-order DE:
(∇6−L1∇4+L2∇2−L3)−G=0, | (52) |
where −G={−e,−θ,−ϕ} and coefficients L1, L2, and L3 are defined as:
L1=1Q7[p5+Q7p1+Q8p2],L2=1Q7[p1p5+Q7Q4Q0+p2p4+Q8p3],L3=1Q7[Q4Q0p5+p3p4]. | (53) |
Here,
p1=Q4+ϖ0+Q5Q7,p2=ϖ0R3−Q6Q7,p3=Q4ϖ0δ3,p4=Q8ϖ0+Q7Q10,p5=Q9Q7−Q8ϖ0R3. | (54) |
The DE (52) can be expressed as the product of three second-order factors:
(∇2−m21)(∇2−m22)(∇2−m23){−e,−θ,−ϕ}=0, | (55) |
where m21, m22, and m23 are roots with positive real parts of the characteristic equation:
m6−L1m4+L2m2−L3=0. | (56) |
The formulation of the solution to Eq (55), under the condition that the variables remain finite as r approaches infinity, can be expressed as follows:
{−θ,−ϕ,−e}=1√r∑3i=1{1,Ψi,Γi}CiK1/2(mir). | (57) |
Here, Ci represent constants that are established based on the boundary conditions, while Ψi and Γi are outlined as follows:
Ψi=m4i−p1m2i+Q4M0p2m2i−p3,Γi=m2i(Q5−Q4Ψi)m2i−Q4,i=1,2,3. | (58) |
By transforming Eq (27) using the Laplace transform and leveraging Eq (59), we can perform an integration with respect to r to derive the displacement ˉu as follows:
ˉu=−1√r∑3i=1ΓimiCiK3/2(mir). | (59) |
In the process of deriving Eq (59), the following relation for the modified Bessel functions of the second kind has been utilized:
ddrK1/2(mir)=12rK1/2(mir)−miK3/2(mir). | (60) |
As Y becomes large, the modified Bessel functions of the second kind behave asymptotically as follows:
K1/2(Y) e−Y√π2Y,K32(Y)=e−Y√π2Y(1+1Y). | (61) |
When the asymptotic relations (50) are substituted into Eqs (46) and (49), we derive:
{−θ,−ϕ,−e}=1√r∑3i=1{1,Ψi,Γi}Cie−mir√π2mir, | (62) |
−u=−1√r∑2i=1ΓimiCie−mir√π2mir(1+12mir). | (63) |
Ultimately, by inserting Eqs (62) and (63) into Eq (46), we are able to determine the distributions of thermal stress.
Once the Laplace transform has been applied to the boundary conditions specified in Eqs (36) and (38), they are reformulated accordingly:
−τrr(r,s)=0,d−ϕ(r,s)dr=0 at r=R, | (64) |
d−θ(r,t)dr=−δ0q0tp8(stp+1)3=−−G(s) at r=R. | (65) |
To ascertain the values of constants C1, C2, and C3, we substitute the functions −τrr, −ϕ, and −θ into the boundary conditions and solve the resulting system of equations. With this approach, we have successfully obtained the solutions across the entire domain within the Laplace transform space.
The Gaver-Stehfest method serves as an efficient numerical technique widely utilized in engineering and applied sciences for inverting the Laplace transform. This approach facilitates the transition from the Laplace domain into the physical domain, addressing both spatial and temporal variables. By combining numerical techniques with the essential features of the Laplace transform, this method provides a reliable way to compute inverse transforms, particularly for functions that are difficult to invert analytically.
The Gaver-Stehfest method can be employed to estimate the values of ℊ(r,t) at different time intervals through the following formula [50]:
ℊ(r,t)≈ln(2)t∑Nk=1ωk−g(r,ntln(2)),k≥1,t>0. | (66) |
In this scenario, ωk signifies the coefficients linked to the individual terms within the series expansion, and N denotes the total count of terms that are incorporated into the expansion; a number that is required to be even. The determination of the coefficients ωk is based exclusively on the quantity of expansion terms k, and they are defined in the manner described as [51]:
ωk=(−1)k+N2∑min{k,N/2}j=k+12jN2(2j)!(N2−j)!(k−j)!(2j−k)!(j−1)!j!,k≥1,1≤N≤k. | (67) |
The Gaver-Stehfest method is particularly advantageous for its use in numerical analysis due to several key benefits. Its primary advantage is its efficiency. The method uses a straightforward summation formula, bypassing the more intricate processes of contour integration or symbolic inversion, thereby substantially decreasing the time required for computations. Furthermore, it is highly adaptable, being able to handle a wide array of Laplace-transformed functions that are commonly found in areas such as heat conduction, diffusion, and fluid dynamics. When it comes to smooth functions in the Laplace domain that are free from singularities or discontinuities, the Gaver-Stehfest method is capable of providing highly precise outcomes [52]. The simplicity of its formula also makes it accessible and easy to incorporate into numerical algorithms.
Nevertheless, the method is not without its drawbacks. A primary issue is numerical instability. As the parameter N is raised, the factorial elements within the coefficients ωk can increase sharply, leading to a decrease in precision. Additionally, if the Laplace-domain function −ℊ(r,s) possesses singularities or discontinuities, the method's accuracy may be compromised. It is crucial to carefully select the value of N. While a higher N can improve precision for smooth functions, extremely large values of N might exacerbate numerical instability, potentially leading to difficulties in achieving dependable results.
This section focuses on the numerical simulation of the thermoelastic behavior of an infinite porous magnesium elastic body with a spherical cavity. The main goal is to validate the nonlocal thermoelastic model, which integrates spatial and temporal effects, fractional derivatives, and phase delays. Through this analysis, we aim to enhance our understanding of the interactions among different field variables and examine how the key model parameters affect the dynamic and transient responses of the material within the medium. For the simulations, we have utilized the unique material properties of porous elastic magnesium, a material commonly employed in mechanical and thermal engineering applications. The details of the material parameters are presented in Table 1:
Property | Symbol | Value | Units |
First Lamé parameter | λ | 2.17×1010 | Nm−2 |
Shear modulus | μ | 3.278×1010 | Nm−2 |
Density | ρ | 1740 | kg/m3 |
Void-strain coupling parameter | 𝓋 | 1.13849×1010 | Nm−2 |
Void coupling parameter | ξ2 | 1.475×1010 | Nm−2 |
Thermal conductivity | K | 1.7×102 | Wm−2K−1 |
Thermal coupling coefficient | β | 2.68×106 | Nm−2K−1 |
Thermal expansion coefficient | αt | 2.5×10−5 | K−1 |
Specific heat capacity | Ce | 1.04×1010 | J kgK−1 |
Reference temperature | T0 | 298 | K |
Void-thermal coupling parameter | ξ1 | 2×106 | Nm−2K−1 |
Void density parameter | χ1 | 1.753×10−15 | m2 |
Void interaction parameter | p | 3.688×10−5 | N |
Wave speed parameter | υ21 | 5.92×106 | m2s−2 |
Thermal diffusivity | η1 | 1.065×1011 | s−1 |
In this study, the parameters governing the system's behavior are crucial for understanding the material's response to thermal and mechanical stimuli. The fractional order (α) is a pivotal factor that determines the extent of memory effects in the material. Lower α values enhance long-term memory, while higher values diminish these effects. The tempering parameter (σ) introduces exponential decay to the material's response, with higher values limiting long-range interactions and wave scattering, and lower values amplifying nonlocal interactions and wave dispersion.
Nonlocal parameters, such as spatial nonlocality (η) and temporal nonlocality (ζ), are essential for capturing size-dependent mechanical behavior and transient thermal effects. The spatial nonlocal parameter (η) dictates the localization of the material's response—smaller η values yield more localized effects, whereas larger values produce smoother, broader responses. The temporal nonlocal parameter (ζ) governs the material's reaction speed, influencing wave propagation rates.
Lastly, phase lags (τq and τθ) are critical for describing time delays in the material's thermal and heat flux responses. These lags contribute to the relaxation behavior, reducing vibration amplitudes; larger phase lags lead to longer relaxation times. Together, these parameters are integral to the model's precision and its capacity to predict the material's behavior under varying conditions.
In this section, we analyze the numerical results to compare the behaviors of physical fields under two types of fractional derivatives: the fractional-Caputo derivative (C) and the tempered-Caputo fractional derivative (CT). By varying the fractional order (α) and the exponential tempering parameter (σ), we aim to clarify the impact of memory effects on the thermoelastic responses of porous materials. The results, presented in Figures 2 to 5, focus on critical physical quantities such as temperature (θ), radial displacement (u), nonlocal thermal stresses (Trr and TΘΘ), and the volume fraction field (ϕ), as influenced by the chosen fractional operators.
Figure 2 specifically illustrates the variation of the temperature field (θ) based on the type of fractional derivative, the fractional order (α), and the tempering parameter (σ). One notable observation is the effect of the fractional order: when employing the Caputo fractional derivative (C-TPL), reducing α from 1.0 to 0.5 causes a significant decrease in temperature values across all radial distances (r). This trend indicates that lower fractional orders enhance memory effects, resulting in more gradual thermal diffusion and lower peak temperatures. Additionally, the tempering parameter (σ) plays a critical role. With the tempered-Caputo fractional derivative (CT-TPL), increasing σ correlates with a further reduction in temperature values. For instance, at α=0.5, as σ increases from 0 to 20, the temperature drops considerably, suggesting that the effect of exponential tempering accelerates thermal memory decay, leading to faster attenuation of temperature profiles.
Moreover, when comparing the C-TPL and CT-TPL models, it becomes evident that the tempered-Caputo fractional derivative consistently results in lower temperature values than the traditional Caputo fractional derivative for the same fractional order (α). This observation highlights the additional damping effect introduced by the tempering parameter (σ), especially pronounced at larger radial distances (r). Lastly, the dependence of temperature (θ) on radial distance (r) is notable, showing a consistent decrease as r increases for all scenarios. However, the decay rate is affected by both the fractional order (α) and tempering parameter (σ). Lower values of α and higher values of σ lead to steeper declines in temperature, reflecting stronger memory effects and quicker thermal relaxation. These observations underscore the critical role of fractional calculus and tempering parameters in modeling thermoelastic responses. The ability to tune memory effects and thermal relaxation through α and σ provides a powerful tool for designing advanced materials with optimized thermoelastic properties. This study highlights how fractional operators, especially the tempered-Caputo derivative, enhance the accuracy and flexibility of modeling thermoelastic behaviors in porous materials.
Figure 3 illustrates the radial displacement (u) and how it is influenced by the type of fractional derivative, as well as the fractional order (α) and the exponential tempering parameter (σ). The examination of these parameters reveals significant trends in how displacement behaves within the porous medium. As the fractional order (α) is adjusted, variations in displacement patterns become evident, suggesting that the nature of memory effects plays a pivotal role in how the material responds. When the tempering parameter (σ) is held constant at 0.0, a reduction in the fractional order (α) from 1.0 to 0.5 leads to a marked decrease in the magnitude of radial displacement (u) across all radial distances (r). For instance, at r = 1, u changes from -0.590416 (for α = 1.0) to -0.511017 (for α = 0.5). This trend suggests that lower fractional orders (α < 1) enhance memory effects, resulting in slower changes in displacement and a greater retention of previous mechanical states. In contrast, as α approaches 1, the displacement stabilizes more rapidly, reflecting a behavior similar to classical elasticity with diminished memory influences.
With a fixed fractional order (α), increasing the tempering parameter (σ) leads to a reduction in the magnitude of radial displacement (u). For example, at r = 1 and α = 0.5, u decreases from -0.511017 (for σ = 0.0) to -0.383079 (for σ = 20). The exponential tempering factor e−σx speeds up the decay of memory effects, resulting in faster stabilization of displacement. This effect enhances the material's response by making it more localized over time and space, which is beneficial for applications that require quick recovery or dynamic stability.
The interplay between the fractional order (α) and tempering parameter (σ) allows for precise tuning of the material's mechanical response. For example, a combination of low α and high σ creates a highly damped response with rapid stabilization, ideal for applications like vibration damping or shock absorption. Conversely, high α combined with low σ results in a more traditional elastic reaction with slower stabilization, beneficial for scenarios needing prolonged deformation.
Initially, at smaller radial distances (r = 1), the displacement values are negative, indicating inward movement due to thermal and mechanical loading near the cavity. The radial displacement (u) shows a clear correlation with radial distance (r). At smaller radial distances (r = 1), displacement values manifest the greatest magnitude due to proximity to the spherical cavity, where thermal and mechanical interactions are most intense. As r increases, the displacement consistently decreases, suggesting that the influence of the cavity lessens with distance. The decay rate is determined by both the fractional order (α) and tempering parameter (σ), with lower α values and higher σ values resulting in steeper declines in displacement, highlighting stronger memory effects and quicker stabilization.
Manipulating the fractional order (α) and the tempering parameter (σ) provides a powerful approach to customizing the mechanical characteristics of porous materials. By fine-tuning these parameters, engineers can develop materials specifically suited for various applications. For instance, in dynamic systems, higher σ values allow for quick stabilization, which is essential for effective vibration damping or dynamic thermal management. Meanwhile, lower α values can improve energy absorption capabilities, making the materials particularly effective for resisting impacts. Additionally, materials designed with lower α can demonstrate reduced heat dissipation, making them well-suited for applications requiring thermal insulation or heat retention.
Figure 4 presents the radial nonlocal thermal stresses (Trr) and examines how these stresses depend on the type of fractional derivative, alongside the fractional order (α) and the exponential tempering parameter (σ). This analysis reveals critical insights into the behavior of thermal stresses in porous materials under varying conditions. Examining the stress distribution specifically near the cavity surface (at r = 1), it is notable that Trr remains zero across all scenarios. This observation aligns with the established boundary condition for a traction-free cavity surface, where no thermal stress is exerted. As one moves just outside of the cavity (r > 1), thermal stress begins to emerge, peaking at r = 1.4 before gradually tapering off with increasing radial distance. This pattern emphasizes the critical role of proximity to the cavity in determining the distribution of thermal stresses within the material. The radial thermal stress (Trr) shows a distinct correlation with radial distance (r). At smaller distances, particularly at r = 1.4, the stress values are at their highest, indicating the significant impact of thermal and mechanical disturbances close to the spherical cavity. As the distance from the cavity increases, stress levels decrease steadily, illustrating that the influence of the cavity weakens with greater separation.
When the tempering parameter (σ) is held constant at 0.0, lowering the fractional order (α) from 1.0 to 0.5 results in a reduction of the radial thermal stress (Trr) across all radial distances (r). For instance, at r = 1.4, Trr decreases from -0.220284 (when α = 1.0) to -0.171524 (when α = 0.5). This trend indicates that lower fractional orders (α < 1) tend to enhance memory effects, which contribute to slower stress relaxation and a longer retention of thermal stresses. As α approaches 1, the decay of stress occurs more swiftly, reflecting behavior characteristics of classical mechanics, where memory effects become less significant.
With a fixed fractional order (α), an increase in the tempering parameter (σ) leads to a notable decrease in the magnitude of radial thermal stress (Trr). For example, at r = 1.4 and α = 0.5, Trr drops from -0.171524 (for σ = 0.0) to -0.072994 (for σ = 20). The exponential tempering factor e−σx plays a critical role in accelerating the decay of memory effects, resulting in faster stress relaxation and a more localized response in the material. This characteristic makes the material particularly advantageous for applications that necessitate quick thermal recovery or minimized residual stresses.
The interaction between α and σ enables precise adjustment of the material's thermoelastic behavior. For instance, utilizing a combination of low α and high σ leads to a significantly damped response characterized by swift stress relaxation. This configuration is particularly beneficial for applications like thermal shock resistance or vibration damping, where quick recovery from mechanical or thermal loads is essential. On the other hand, employing a combination of high α and low σ yields a more traditional elastic response, marked by slower stress relaxation. This setup is ideal for scenarios that demand sustained stress retention, making it suitable for applications where stability and long-term performance are critical.
Figure 5 illustrates the variation of the volume fraction field (ϕ) in relation to the radius within a porous infinite elastic material featuring a spherical cavity. This representation highlights the dependence of the volume fraction on several factors, including the type of fractional derivative, the fractional order (α), and the exponential smoothing factor (σ). The volume fraction field (ϕ) clearly varies with radial distance (r). At smaller distances, specifically at r = 1, the volume fraction ϕ reaches its peak values, highlighting the significant impact of thermal and mechanical disturbances near the spherical cavity. As one moves further away and r increases, there is a steady decline in the volume fraction ϕ, signifying that the effect of the cavity weakens with increasing distance.
When the tempering parameter (σ) is held constant at 0.0, decreasing the fractional order (α) from 1.0 to 0.5 results in a reduction of the volume fraction field (ϕ) at all radial distances (r). For instance, at r = 1, ϕ drops from 0.177517 (with α = 1.0) to 0.13496 (with α = 0.5). This trend indicates that lower fractional orders (α < 1) strengthen memory effects, which leads to a slower evolution of the volume fraction field ϕ and a longer retention of previous conditions. As α nears 1, the volume fraction field stabilizes more rapidly, showcasing behavior akin to classical mechanics, where memory effects play a diminished role.
With a fixed fractional order (α), an increase in the tempering parameter (σ) markedly decreases the volume fraction field (ϕ). For example, at r = 1 and α = 0.5, ϕ reduces from 0.13496 (for σ = 0.0) to 0.096955 (for σ = 20). The exponential tempering factor e−σx accelerates the fading of memory effects, resulting in a quicker stabilization of the volume fraction field. This characteristic lends the material a more localized temporal and spatial response, making it particularly suited for applications that require swift recovery or dynamic stability.
The internal length scale (η) serves as a vital metric that indicates how the surrounding material points influence a specific location within a substance. This parameter is particularly significant in the context of wave dispersion and size-dependent phenomena, which are commonly observed in nanostructures. Here, the microstructural characteristics greatly shape the material's overall response to external stimuli. The characteristic time scale (ζ), on the other hand, captures the essence of time-delayed interactions, accounting for the lag between the application of loads and the resulting responses of the material. This aspect becomes crucial when modeling dynamic situations, such as the heat fluxes generated by laser pulses, where transient effects play a dominant role.
In the classical framework where both η and ζ are set to zero, the system defaults to the local two-phase lag (FTPL) thermoelastic model that incorporates fractional derivatives. In this scenario, the responses are strictly determined by the local thermal and mechanical properties, effectively ignoring the influences of spatial nonlocality and time scale effects. Stress and strain fields rely solely on nearby temperature gradients and deformations. However, this approach has its limitations, as it fails to account for wave dispersion and size-dependent effects, overlooking the critical transient delays in heat conduction and stress propagation, especially pertinent in nanoporous materials and scenarios involving laser-induced heat fluxes. In contrast, the nonlocal framework, characterized by non-zero values for both η and ζ, introduces a more nuanced perspective. It leverages the KG operator to incorporate both spatial and temporal nonlocality, enabling responses to be shaped by not just neighboring points but also by historical states of the material. This approach effectively captures size-dependent phenomena, such as stress concentrations near voids and wave scattering, while also modeling time-delayed thermal and mechanical effects. As a result, this framework provides improved accuracy in dynamic scenarios, making it a more robust choice for understanding complex material behaviors.
Figure 6 illustrates the impact of the intrinsic length scale (η) and characteristic time scale (ζ) on the non-dimensional temperature distribution (θ) as a function of radial distance (r). The figure demonstrates the interaction between spatial nonlocality (η) and temporal nonlocality (ζ), showing how these parameters influence the temperature distribution within a porous medium.
In the standard scenario, where both η and ζ are set to zero, the system behaves in line with classical local models that exclude any spatial or temporal nonlocal effects. Here, the highest temperature values emerge because heat diffusion operates without the complexities introduced by size or time-dependent memory effects. Conversely, in nonlocal scenarios where η is greater than zero and ζ is non-negative, there is a noticeable decline in temperature readings. Increasing the internal length scale (η), which introduces spatial nonlocality, along with the characteristic time scale (ζ), which accounts for temporal nonlocality, leads to diminished temperature values. This reduction arises from the interplay of spatial averaging and the delayed propagation of heat, resulting in slower heat diffusion and lower peak temperatures.
Radially, the temperature consistently decreases as the radial distance increases. Close to the cavity, at a radial distance of r = 1, the temperature is significantly higher due to concentrated thermal interactions. However, as one moves outward to larger distances, such as r = 5, the temperature experiences substantial decay, reflecting smaller values. Specifically examining the impact of η and ζ, in the absence of nonlocal effects (η = 0.000, ζ = 0.000), the temperature follows a predictable decreasing trend with respect to radial distance, illustrating classical heat conduction behavior. As η increases, the temperature consistently decreases across all radial distances, especially pronounced near the cavity. For instance, at r = 1, the temperature drops from 0.796883 (η = 0.000) to 0.739311 when η is increased to 0.004 and ζ to 0.003. In a similar vein, higher ζ values further reduce temperature, notably beyond r = 1.4, as evidenced by a decrease from 0.0163379 (η = 0.000, ζ = 0.000) to 0.0128245 (η = 0.004, ζ = 0.005) at r = 3.
The intrinsic length scale (η) plays a pivotal role in introducing spatial nonlocality, affecting each location's temperature by factoring in the influence of surrounding material points. This leads to increased spatial averaging, effectively smoothing temperature gradients and lowering peak values. Most prominently, this effect is observed near the cavity, where localized thermal conditions are spread over a broader area. As η escalates from 0.002 to 0.006, the temperature decreases across all radial distances, with the most significant reductions appearing close to the cavity.
In terms of the characteristic time scale (ζ), this parameter brings temporal nonlocality into the model, taking into account the delay in heat propagation over time. Higher ζ values result in slower thermal diffusion, which manifests as lower temperature readings because it takes longer for heat to reach an equilibrium state. This impact is particularly observable at greater radial distances, such as those beyond r = 2, where the delays play a significant role. As ζ progresses from 0.000 to 0.005, a decrease in temperature becomes evident, particularly in areas farther from the cavity.
The combined influence of η and ζ accentuates the nonlocal characteristics of the model. Elevated η values work to decrease the temperature by spatially averaging the thermal profile, while increased ζ values contribute to this decline by introducing delays in heat diffusion. For instance, at r = 1, the temperature descends from 0.796883 (η = 0.000, ζ = 0.000) to 0.625518 (η = 0.004, ζ = 0.005), illustrating the compounded effects of spatial and temporal nonlocality.
The findings hold several practical implications. For thermal insulation, materials characterized by higher values of η and ζ display slower heat diffusion and reduced peak temperatures, making them particularly suitable for insulation purposes. In dynamic thermal management systems, the delays in heat propagation associated with higher ζ values could be advantageous for applications requiring gradual heat dissipation, such as in electronics or nanostructured components. Lastly, incorporating nonlocal parameters such as η and ζ is crucial in accurately modeling the behavior of nanoporous materials, where size and time-dependent effects are critical considerations in their thermal responses.
Figure 7 examines the impact of the intrinsic length scale (η) and characteristic time scale (ζ) on the dimensionless displacement (u) within an infinite nanoelastic medium that contains a spherical cavity. The curves in Figure 7 illustrate the behavior of u as a function of the radial distance (r) for various combinations of η and ζ. The radial displacement (u) exhibits complex patterns that vary with the radial distance (r), the intrinsic length scale (η), and the characteristic time scale (ζ). At smaller radial distances (around r = 1), the displacement is negative, indicating inward movement due to thermal and mechanical influences near the cavity. As r increases, u becomes positive, signifying outward movement caused by the same loading effects. The magnitude of u diminishes with increasing r, showing that the cavity's influence lessens at greater distances.
In the standard scenario (η = 0, ζ = 0), the system behaves like a classical local model, displaying the largest displacement magnitudes. Near r = 1, significant negative values are observed due to the absence of spatial or temporal nonlocal effects, meaning displacement (u) depends only on local mechanical properties. However, when η and ζ are non-zero, the displacement magnitude decreases both near the cavity and at larger distances. This reduction is due to the smoothing effect of spatial nonlocality (η) and the damping effect introduced by temporal nonlocality (ζ).
The intrinsic length scale (η) plays a critical role in shaping displacement (u) behavior. As η increases, it promotes spatial averaging, reducing the localization of mechanical responses and resulting in a smoother displacement field (u), especially near the cavity. For example, at r = 1, as η increases from 0.000 to 0.006 (with ζ = 0.001), the inward displacement (u) decreases from −0.383079 to −0.338451. At larger distances (r > 3), the influence of η diminishes as displacement values become inherently small. For instance, at r = 4.6, the displacement u decreases slightly from 0.00159686 (η = 0.000) to 0.00141083 (η = 0.006, ζ = 0.001).
Similarly, the characteristic time scale (ζ) introduces time-delay effects that further dampen mechanical responses. As ζ increases, the stabilization of the displacement field (u) slows down, particularly at larger radial distances. For example, near r = 1, increasing ζ reduces inward displacement (u) from −0.355403 (η = 0.004, ζ = 0.003) to −0.320056 (η = 0.004, ζ = 0.005). At greater distances, such as r = 5, the damping effect of ζ becomes more pronounced, with u decreasing from 0.00130603 (η = 0.004, ζ = 0.003) to 0.00117 (η = 0.004, ζ = 0.005).
The interaction between η and ζ amplifies the overall nonlocal effects in the model. Higher η values redistribute mechanical responses spatially, while higher ζ values introduce delays, enhancing displacement dampening. For instance, at r = 1, the displacement (u) decreases from −0.383079 (classical case, η = 0.000, ζ = 0.000) to −0.320056 (η = 0.004, ζ = 0.005), illustrating the combined effects of these parameters.
These findings have significant practical implications. Materials with higher η and ζ values exhibit reduced displacement magnitudes and faster stabilization, making them ideal for applications requiring effective vibration damping or shock absorption. Additionally, the influence of η and ζ is crucial for modeling nanostructured materials, where size- and time-dependent effects dominate. The ability to control temporal nonlocality (ζ) also allows for precise management of time-delayed mechanical responses, enabling the design of materials with enhanced dynamic stability for various advanced applications.
Figure 8 illustrates the nonlocal, non-dimensional thermal stress distribution (Trr) for various values of the intrinsic length scale (η) and characteristic time scale (ζ). The figure highlights the influence of spatial and temporal nonlocality on the size- and time-dependent stress behavior of materials with nanoscale characteristics. From the analysis presented in Figure 8, it is clear that at a radial distance of r = 1, the thermal radial stress (Trr) remains at zero, regardless of the configurations of η and ζ. This observation aligns perfectly with the condition of a free-traction boundary, where no stress is exerted at that surface. Beyond this point, Trr takes on negative values, which signify compressive stress arising from thermal expansion effects. Moreover, as the radial distance increases, the magnitude of Trr decreases, illustrating the diminishing impact of the cavity with distance.
In the benchmark scenario where both η and ζ equal zero, the results in Figure 8 display adherence to the classical local model, revealing that the thermal stress Trr achieves its peak value in the absence of nonlocal effects. When delving into nonlocal scenarios where η is positive and ζ is non-negative, Figure 8 indicates a noticeable reduction in the magnitude of Trr, both close to the cavity at r = 1.4 and further out at r = 5. This decrease is caused by the smoothing effects related to spatial nonlocality (η) and the damping effects associated with temporal nonlocality (ζ).
As η increases, the stress gradients illustrated in Figure 8 become more uniform, leading to a reduction in stress concentrations near the cavity's surface. For example, at r = 1.4, the thermal stress Trr diminishes from -0.141646 (when η = 0.000) to -0.111186 (at η = 0.006 and ζ = 0.001). Even at a greater radial distance of r = 5, the stress Trr decreases from -0.000533963 (η = 0.000) to -0.000419138 (η = 0.006, ζ = 0.001), though the effect is less pronounced. The characteristic time scale represented by ζ introduces temporal nonlocality, which accounts for delays in stress propagation. In Figure 8, it is depicted that larger values of ζ lead to lower stress magnitudes, indicating a delay in the system's response to thermal and mechanical loads. For instance, at r = 1.4, the thermal stress Trr falls from -0.125144 (with η = 0.004 and ζ = 0.003) to -0.10307 (with η = 0.004 and ζ = 0.005). This damping effect becomes even more evident at r = 5, where the stress Trr reduces from -0.000471757 (η = 0.004, ζ = 0.003) to -0.000388546 (η = 0.004, ζ = 0.005).
Additionally, Figure 8 highlights the interactive effects of η and ζ on the thermal stress Trr. Higher values of η contribute to the spatial redistribution of thermal effects, while elevated ζ values further diminish stress through delayed mechanical responses. This interaction is particularly significant near the cavity, such as at r = 1.4 but begins to wane at extended distances. For example, the thermal stress Trr decreases from -0.141646 (in the absence of nonlocal effects) to -0.10307 (when incorporating both η = 0.004 and ζ = 0.005), illustrating how spatial and temporal nonlocality cooperatively smooth and weaken the stress distributions.
These insights carry substantial practical implications. Materials with increased values of η and ζ tend to show diminished stress concentrations and smoother stress distributions, significantly enhancing resistance to thermal shocks. In the realm of nanostructured materials, recognizing the contributions of spatial and temporal nonlocality is essential for accurately modeling stress behavior, especially in contexts where size and time-dependent phenomena are critical. Moreover, the influence of temporal nonlocality (ζ) facilitates controlled stress relaxation, thus boosting the dynamic stability of materials subjected to sudden loading conditions (an important attribute for numerous engineering applications).
Figure 9 illustrates the impact of the intrinsic length scale (η) and characteristic time scale (ζ) on the distribution of the non-dimensional volume fraction field (ϕ) within an infinite elastic medium containing a spherical cavity. This analysis investigates how spatial nonlocality (η) and temporal nonlocality (ζ) affect the evolution of ϕ, which represents the void or porosity field. The volume fraction field (ϕ) is a critical parameter in porous materials, as it indicates the distribution of voids under thermoelastic loading conditions.
Close to the cavity (at r = 1), the volume fraction field (ϕ) reaches its highest value due to the significant thermal and mechanical interactions occurring near the spherical cavity. However, as the parameters η and ζ increase, ϕ begins to decline, demonstrating the dampening effects of both spatial and temporal nonlocality. Moving farther from the cavity (when r > 1), ϕ consistently decreases as the radial distance increases, indicating that the impact of the cavity lessens with distance. The rate at which ϕ decays is influenced by both η and ζ, as larger values of these parameters promote smoother and more subdued distributions.
In the standard scenario (η = 0.000, ζ = 0.000), the results illustrate a classical local model where ϕ is maximized across all radial distances since no nonlocal effects are considered. Here, the behavior is predominantly determined by local thermal and mechanical properties, showcasing sharp gradients near the cavity. When η and ζ are greater than zero, the introduction of spatial (η) and temporal (ζ) nonlocality leads to a noticeable reduction in ϕ values throughout all radial distances, especially in close proximity to the cavity. Increasing η and ζ smooths the spatial layout of ϕ, reducing peak values near r = 1 and signifying a more gradual evolution of voids due to the influences of nonlocal interactions.
Understanding the physical interpretation of η shows that it introduces spatial nonlocality, where the state at a particular point is affected not only by local conditions but also by the neighboring areas. An increase in η results in smoother distributions, as the void fraction field becomes more averaged out. Near r = 1, it is noted that an increase in η leads to a decrease in ϕ. Specifically, ϕ diminishes from 0.13496 when η = 0.000 to 0.105938 with η = 0.006 and ζ = 0.001. As the distance from the cavity increases (r > 3), the influence of η diminishes; for example, at r = 3, ϕ decreases from 0.00302298 with η = 0.000 to 0.00237291 with η = 0.006 and ζ = 0.001.
The influence of ζ introduces temporal nonlocality, accounting for the delays that occur in the evolution of voids. Higher values of ζ effectively dampen the evolution of ϕ, leading to a slower response to thermal and mechanical stimuli. As observed near the cavity (at r = 1), increasing ζ decreases ϕ even further; for instance, it shifts from 0.119238 (η = 0.004, ζ = 0.003) to 0.0982058 (η = 0.004, ζ = 0.005). At greater distances (specifically r = 5), the damping effect of ζ is particularly significant, with ϕ declining from 0.000166467 (η = 0.004, ζ = 0.003) to 0.000137104 (η = 0.004, ζ = 0.005).
The interplay between η and ζ yields a powerful impact. Spatial nonlocality (η) reduces the highest values of ϕ and smooths out its distribution, while temporal nonlocality (ζ) adds further damping by inducing time delays. This combination is most pronounced near the cavity, where the evolution of the void fraction field occurs more gradually due to the synergistic effects of spatial and temporal influences. For example, at r = 1, ϕ decreases from 0.13496 (with η = 0.000 and ζ = 0.000) to 0.0982058 (with η = 0.004 and ζ = 0.005).
The insights gained from these findings offer substantial practical implications. With larger values of η and ζ, the resulting void fraction distributions become smoother and more damped. This characteristic is crucial for controlling localized void growth and enhancing the structural integrity of porous materials. Additionally, the dampening implications attributed to η and ζ help minimize peak values of ϕ, reducing stress concentrations and bolstering the overall stability of materials under thermal and mechanical loading. In the design of nanostructures, recognizing the role of spatial and temporal nonlocality is essential for accurately modeling void evolution, especially in scenarios where size and time-dependent effects play a dominant role.
This study examined the Caputo-tempered two-phase lag (FCT-TPL) heat conduction model and its influence on thermoelastic vibrations in materials containing spherical cavities and voids. By integrating fractional calculus and tempered derivatives, along with spatial and temporal nonlocality, the model effectively bridges the gap between classical and generalized thermoelastic theories, offering a robust framework for accurately analyzing size-dependent phenomena and wave scattering.
The results demonstrate that the distinction between the fractional Caputo (C) and tempered Caputo (CT) derivatives highlights the versatility of fractional calculus in modeling thermoelastic behaviors. Additionally, it is revealed that the fractional-order parameter and the stiffness parameter in differential actuators serve as effective tools for adjusting memory effects, thermal relaxation, and mechanical stability in porous materials. Notably, it is observed that employing low-fractional-order values in combination with high-stiffness parameter values leads to rapid attenuation of responses, making this configuration particularly advantageous for applications requiring dynamic stability and efficient vibration damping. Conversely, utilizing high-fractional-order values alongside low stiffness parameter values results in prolonged response durations, which is ideal for applications involving thermal insulation and long-term material performance. This capability to fine-tune material behaviors through these parameters underscores the immense potential of fractional calculus in real-world engineering and material design applications.
From the study and analysis, it is evident that the integration of the parameters η and ζ is crucial for accurately modeling size- and time-dependent phenomena in porous nanomaterials, where localized effects are highly significant. Higher values of spatial nonlocality are shown to enhance spatial averaging, resulting in smoother gradients (e.g., temperature, displacement, and stress) while simultaneously reducing peak magnitudes around localized features such as cavities. Moreover, increasing values of temporal nonlocality introduces time-dependent damping, which slows heat diffusion, stabilizes displacement, and influences void evolution.
The findings conclude that the interaction between spatial nonlocality and temporal nonlocality substantially improves the precision of material modeling under thermoelastic loading. These parameters effectively capture size-dependent behaviors and time-delayed responses, establishing them as indispensable tools in the design of materials with enhanced thermal, mechanical, and structural properties.
While the Caputo-tempered two-phase lag (FCT-TPL) heat conduction model effectively addresses size- and time-dependent phenomena, it faces notable challenges. These include high computational complexity from fractional calculus and nonlocal parameters, sensitivity to parameter calibration, and limited applicability to heterogeneous or anisotropic materials due to its assumptions of homogeneity and isotropy. Additionally, insufficient experimental validation across diverse materials and conditions limits its practical use.
Future advancements should focus on reducing computational demands through efficient numerical methods or machine learning, extending the model to handle irregular geometries and anisotropy, and integrating it with additional physical phenomena. Comprehensive experimental validation, refined parameter identification techniques, and exploration of applications in emerging technologies are essential. Incorporating nonlinear effects could further enhance its accuracy and relevance, particularly under extreme conditions, solidifying its value in engineering applications.
Kareem Alanazi and Ahmed E. Abouelregal contributed to the conceptualization, methodology, and investigation of the study; Kareem Alanazi focused on data analysis, visualization, and preparing the initial manuscript draft; Ahmed E. Abouelregal provided oversight, critical revisions, and theoretical insights to enhance the interpretation of results. Both authors reviewed and approved the final manuscript. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors declare no conflict of interest.
[1] |
Alonso-Carrera J, Raurich X (2018) Labor mobility, structural change and economic growth. J Macroecon 56: 292-310. doi: 10.1016/j.jmacro.2018.03.002
![]() |
[2] |
Beutels P, Jia N, Zhou QY, et al. (2009) The economic impact of SARS in Beijing, China. Trop Med Int Health 14: 85-91. doi: 10.1111/j.1365-3156.2008.02210.x
![]() |
[3] |
Brancaccio E, Garbellini N, Giammetti R (2018) Structural labour market reforms, GDP growth and the functional distribution of income. Struct Change Econ Dyn 44: 34-45. doi: 10.1016/j.strueco.2017.09.001
![]() |
[4] |
Chambers J, Barker K, Rouse A (2012) Reflections on the UK's approach to the 2009 swine flu pandemic: conflicts between national government and the local management of the public health response. Health Place 18: 737-745. doi: 10.1016/j.healthplace.2011.06.005
![]() |
[5] |
Cutrini E (2019) Economic integration, structural change, and uneven development in the European Union. Struct Change Econ Dyn 50: 102-113. doi: 10.1016/j.strueco.2019.06.007
![]() |
[6] |
Freire C (2019) Economic diversification: A model of structural economic dynamics and endogenous technological change. Struct Change Econ Dyn 49: 13-28. doi: 10.1016/j.strueco.2019.03.005
![]() |
[7] |
Guerard JB, Markowitz H, Xu G (2015) Earnings forecasting in a global stock selection model and efficient portfolio construction and management. Int J Forecasting 31: 550-560. doi: 10.1016/j.ijforecast.2014.10.003
![]() |
[8] | Gates B (2020) Responding to COVID-19-a once in-a-century pandemic? New Engl J Med. |
[9] |
Iamsiraroj S (2016) The foreign direct investment-economic growth nexus. Int Rev Econ Financ 42: 116-133. doi: 10.1016/j.iref.2015.10.044
![]() |
[10] |
Landesmann MA, Stöllinge R (2019) Structural change, trade and global production networks: An 'appropriate industrial policy' for peripheral and catching-up economies. Struct Change Econ Dyn 48: 7-23. doi: 10.1016/j.strueco.2018.04.001
![]() |
[11] |
Nicoll A, Ammon A, Amato A, et al. (2010) Experience and lessons from surveillance and studies of the 2009 pandemic in Europe. Public Health 124: 14-23. doi: 10.1016/j.puhe.2009.12.001
![]() |
[12] |
Nkeki CI (2018) Optimal investment risks and debt management with backup security in a financial crisis. J Comput Appl Math 338: 129-152. doi: 10.1016/j.cam.2018.01.032
![]() |
[13] |
Pasquini-Descomps H, Brender N, Maradan D (2017) Value for money in H1N1 influenza: a systematic review of the cost-effectiveness of pandemic interventions. Value Health 20: 819-827. doi: 10.1016/j.jval.2016.05.005
![]() |
[14] |
Samaniego RM, Sun JY (2016) Productivity growth and structural transformation. Rev Econ Dyn 21: 266-285. doi: 10.1016/j.red.2015.06.003
![]() |
[15] |
Singhal T (2020) A Review of Coronavirus Disease-2019 (COVID-19). Indian J Pediatrics, 281-286. doi: 10.1007/s12098-020-03263-6
![]() |
[16] |
Sukharev OS (2020) Structural analysis of income and risk dynamics in models of economic growth. Quant Financ Econ 4: 1-18 doi: 10.3934/QFE.2020001
![]() |
[17] |
Sukharev OS (2019) The restructuring of the investment portfolio: the risk and effect of the emergence of new combinations. Quant Financ Econ 3: 390-411. doi: 10.3934/QFE.2019.2.390
![]() |
[18] | Zhao S, Chen H (2020) Modeling the epidemic dynamics and control of COVID-19 outbreak in China. Quant Biology. |
[19] |
Vu KM (2017) Structural change and economic growth: Empirical evidence and policy insights from Asian economies. Struct Change Econ Dyn 41: 64-77. doi: 10.1016/j.strueco.2017.04.002
![]() |
[20] | Xie B, He D, Mercer T, et al. (2020) Global health crises are also information crises: A call to action. J Assoc Inf Sci Technol. |
[21] | Yu KDS, Aviso KB (2020) Modelling the Economic Impact and Ripple Effects of Disease Outbreaks. Process Integration Optim Sustainability. |
1. | Zhongjian Chen, Xiancong Huang, Yun Gao, Su Zeng, Weimin Mao, Plasma-metabolite-based machine learning is a promising diagnostic approach for esophageal squamous cell carcinoma investigation, 2021, 11, 20951779, 505, 10.1016/j.jpha.2020.11.009 | |
2. | Maxim Khotimchenko, Vladlena Tiasto, Aleksandra Kalitnik, Maria Begun, Rodion Khotimchenko, Ekaterina Leonteva, Igor Bryukhovetskiy, Yuri Khotimchenko, Antitumor potential of carrageenans from marine red algae, 2020, 246, 01448617, 116568, 10.1016/j.carbpol.2020.116568 | |
3. | Shiyao Zheng, Nan Lin, Qing Wu, Hongxin He, Chunkang Yang, Prognostic model construction and validation of esophageal cancer cellular senescence-related genes and correlation with immune infiltration, 2023, 10, 2296-875X, 10.3389/fsurg.2023.1090700 | |
4. | Claudia Ceci, Maria Grazia Atzori, Pedro Miguel Lacal, Grazia Graziani, Role of VEGFs/VEGFR-1 Signaling and Its Inhibition in Modulating Tumor Invasion: Experimental Evidence in Different Metastatic Cancer Models, 2020, 21, 1422-0067, 1388, 10.3390/ijms21041388 | |
5. | Masaki Shimizu, Yu-ichiro Koma, Hiroki Sakamoto, Shuichi Tsukamoto, Yu Kitamura, Satoshi Urakami, Kohei Tanigawa, Takayuki Kodama, Nobuhide Higashino, Mari Nishio, Manabu Shigeoka, Yoshihiro Kakeji, Hiroshi Yokozaki, Metallothionein 2A Expression in Cancer-Associated Fibroblasts and Cancer Cells Promotes Esophageal Squamous Cell Carcinoma Progression, 2021, 13, 2072-6694, 4552, 10.3390/cancers13184552 | |
6. | Hamed Gheibollahi, Amirreza Dehghanian, Negar Taheri, Saeid Tavanafar, Seyede Sona Mousavi, Hamidreza Abbasi, Mohammadreza Sasani, Kissing tumors; the concurrent diagnosis of leiomyosarcoma and squamous cell carcinoma of the esophagus, 2021, 21, 1471-2482, 10.1186/s12893-020-01031-z | |
7. | Eduardas Cicinskas, Maria A. Begun, Vladlena A. Tiasto, Andrei S. Belousov, Valeria V. Vikhareva, Valeria A. Mikhailova, Alexandra A. Kalitnik, In vitro antitumor and immunotropic activity of carrageenans from red algae Chondrus armatus and their low‐molecular weight degradation products , 2020, 108, 1549-3296, 254, 10.1002/jbm.a.36812 | |
8. | Jamielyn DC Cruz, David Paculdo, Divya Ganesan, Meredith Baker, Rebecca J Critchley-Thorne, Nicholas J Shaheen, Sachin Wani, John W Peabody, Clinical variation in surveillance and management of Barrett’s esophagus: A cross-sectional study of gastroenterologists and gastrointestinal surgeons, 2022, 101, 1536-5964, e32187, 10.1097/MD.0000000000032187 | |
9. | Xiaoxiao Zeng, Ling Ye, Mingying Luo, Danli Zeng, Yang Chen, Prognostic value of pretreatment systemic immune-inflammation index in Chinese esophageal squamous cell carcinoma patients receiving radical radiotherapy: A meta-analysis, 2023, 102, 0025-7974, e34117, 10.1097/MD.0000000000034117 | |
10. | Lihong Wang, Xiaosa Wang, Yangyang Wu, Jingjing Wang, Wenping Zhou, Jianyao Wang, Haoran Guo, Na Zhang, Lufang Zhang, Xuanyu Hu, Yuanyuan Zhao, Jinxin Miao, Zifang Zhang, Louisa S. Chard Dunmall, Danhua Zhang, Nicholas R. Lemoine, Zhenguo Cheng, Yaohe Wang, A novel microenvironment regulated system CAR-T (MRS.CAR-T) for immunotherapeutic treatment of esophageal squamous carcinoma, 2023, 03043835, 216303, 10.1016/j.canlet.2023.216303 | |
11. | Başak Can, Fuzail Majoo, Arda Öztürkcan, CONSUMER ACCEPTANCE, ATTITUDE AND KNOWLEDGE STUDIES ON ALTERNATIVE PROTEIN SOURCES: INSIGHT REVIEW, 2024, 1300-3070, 682, 10.15237/gida.GD24030 | |
12. | Guangkun Pei, Zhuoran Liang, Bianli Gu, Linlin Shi, Ze‐Xian Liu, Shegan Gao, Molecular subtyping and precision therapy for esophageal cancer, 2025, 2832-6245, 10.1002/INMD.20240105 |
Property | Symbol | Value | Units |
First Lamé parameter | λ | 2.17×1010 | Nm−2 |
Shear modulus | μ | 3.278×1010 | Nm−2 |
Density | ρ | 1740 | kg/m3 |
Void-strain coupling parameter | 𝓋 | 1.13849×1010 | Nm−2 |
Void coupling parameter | ξ2 | 1.475×1010 | Nm−2 |
Thermal conductivity | K | 1.7×102 | Wm−2K−1 |
Thermal coupling coefficient | β | 2.68×106 | Nm−2K−1 |
Thermal expansion coefficient | αt | 2.5×10−5 | K−1 |
Specific heat capacity | Ce | 1.04×1010 | J kgK−1 |
Reference temperature | T0 | 298 | K |
Void-thermal coupling parameter | ξ1 | 2×106 | Nm−2K−1 |
Void density parameter | χ1 | 1.753×10−15 | m2 |
Void interaction parameter | p | 3.688×10−5 | N |
Wave speed parameter | υ21 | 5.92×106 | m2s−2 |
Thermal diffusivity | η1 | 1.065×1011 | s−1 |
Property | Symbol | Value | Units |
First Lamé parameter | λ | 2.17×1010 | Nm−2 |
Shear modulus | μ | 3.278×1010 | Nm−2 |
Density | ρ | 1740 | kg/m3 |
Void-strain coupling parameter | 𝓋 | 1.13849×1010 | Nm−2 |
Void coupling parameter | ξ2 | 1.475×1010 | Nm−2 |
Thermal conductivity | K | 1.7×102 | Wm−2K−1 |
Thermal coupling coefficient | β | 2.68×106 | Nm−2K−1 |
Thermal expansion coefficient | αt | 2.5×10−5 | K−1 |
Specific heat capacity | Ce | 1.04×1010 | J kgK−1 |
Reference temperature | T0 | 298 | K |
Void-thermal coupling parameter | ξ1 | 2×106 | Nm−2K−1 |
Void density parameter | χ1 | 1.753×10−15 | m2 |
Void interaction parameter | p | 3.688×10−5 | N |
Wave speed parameter | υ21 | 5.92×106 | m2s−2 |
Thermal diffusivity | η1 | 1.065×1011 | s−1 |