In this paper, we investigate the semi-analytical solution of Kaup-Kupershmidt equations with the help of a modified method known as the new iteration transformation technique. This method combines the Yang transform and the new iteration technique. The nonlinear terms can be calculated straightforwardly by a new iteration method. The numerical simulation results have been presented to demonstrate the reliability and validity of the proposed approach. The result confirms that the suggested technique is the best tool for dealing with any nonlinear problems arising in technology and science. In addition, in terms of figures for varying fractional order, the physical behavior of new iteration transformation technique solutions has been shown and the numerical simulation is also exhibited. The solutions of the new iteration transformation technique reveal that the projected technique is reliable, competitive and powerful for studying complex nonlinear fractional type models.
Citation: M. Mossa Al-Sawalha, Rasool Shah, Kamsing Nonlaopon, Imran Khan, Osama Y. Ababneh. Fractional evaluation of Kaup-Kupershmidt equation with the exponential-decay kernel[J]. AIMS Mathematics, 2023, 8(2): 3730-3746. doi: 10.3934/math.2023186
[1] | Muhammad Suhail Aslam, Mohammad Showkat Rahim Chowdhury, Liliana Guran, Isra Manzoor, Thabet Abdeljawad, Dania Santina, Nabil Mlaiki . Complex-valued double controlled metric like spaces with applications to fixed point theorems and Fredholm type integral equations. AIMS Mathematics, 2023, 8(2): 4944-4963. doi: 10.3934/math.2023247 |
[2] | Afrah Ahmad Noman Abdou . Chatterjea type theorems for complex valued extended b-metric spaces with applications. AIMS Mathematics, 2023, 8(8): 19142-19160. doi: 10.3934/math.2023977 |
[3] | Asifa Tassaddiq, Jamshaid Ahmad, Abdullah Eqal Al-Mazrooei, Durdana Lateef, Farha Lakhani . On common fixed point results in bicomplex valued metric spaces with application. AIMS Mathematics, 2023, 8(3): 5522-5539. doi: 10.3934/math.2023278 |
[4] | Jamshaid Ahmad, Abdullah Eqal Al-Mazrooei, Hassen Aydi, Manuel De La Sen . Rational contractions on complex-valued extended b-metric spaces and an application. AIMS Mathematics, 2023, 8(2): 3338-3352. doi: 10.3934/math.2023172 |
[5] | Afrah Ahmad Noman Abdou . A fixed point approach to predator-prey dynamics via nonlinear mixed Volterra–Fredholm integral equations in complex-valued suprametric spaces. AIMS Mathematics, 2025, 10(3): 6002-6024. doi: 10.3934/math.2025274 |
[6] | Saif Ur Rehman, Iqra Shamas, Shamoona Jabeen, Hassen Aydi, Manuel De La Sen . A novel approach of multi-valued contraction results on cone metric spaces with an application. AIMS Mathematics, 2023, 8(5): 12540-12558. doi: 10.3934/math.2023630 |
[7] | Ahmed Alamer, Faizan Ahmad Khan . Boyd-Wong type functional contractions under locally transitive binary relation with applications to boundary value problems. AIMS Mathematics, 2024, 9(3): 6266-6280. doi: 10.3934/math.2024305 |
[8] | Khaled Berrah, Abdelkrim Aliouche, Taki eddine Oussaeif . Applications and theorem on common fixed point in complex valued b-metric space. AIMS Mathematics, 2019, 4(3): 1019-1033. doi: 10.3934/math.2019.3.1019 |
[9] | Sumaiya Tasneem Zubair, Kalpana Gopalan, Thabet Abdeljawad, Nabil Mlaiki . Novel fixed point technique to coupled system of nonlinear implicit fractional differential equations in complex valued fuzzy rectangular b-metric spaces. AIMS Mathematics, 2022, 7(6): 10867-10891. doi: 10.3934/math.2022608 |
[10] | Fatima M. Azmi, Nabil Mlaiki, Salma Haque, Wasfi Shatanawi . Complex-valued controlled rectangular metric type spaces and application to linear systems. AIMS Mathematics, 2023, 8(7): 16584-16598. doi: 10.3934/math.2023848 |
In this paper, we investigate the semi-analytical solution of Kaup-Kupershmidt equations with the help of a modified method known as the new iteration transformation technique. This method combines the Yang transform and the new iteration technique. The nonlinear terms can be calculated straightforwardly by a new iteration method. The numerical simulation results have been presented to demonstrate the reliability and validity of the proposed approach. The result confirms that the suggested technique is the best tool for dealing with any nonlinear problems arising in technology and science. In addition, in terms of figures for varying fractional order, the physical behavior of new iteration transformation technique solutions has been shown and the numerical simulation is also exhibited. The solutions of the new iteration transformation technique reveal that the projected technique is reliable, competitive and powerful for studying complex nonlinear fractional type models.
Contemporary theoretical mechanics and engineering technology are mainly concerned with either the isotropic or the anisotropic behavior of elastic materials and structures [1,2]. Specifically, the classical theory of linear elasticity is considered to be an indivisible branch of the more general field of non-linear elasticity [3], wherein the mathematical description of the physical quantities, associated with such media [4,5], incorporates the fundamental characteristics of the corresponding physical problems. Therein, the relationship between the strain and the stress of any solid body, subjected to external forces, is considered to be linear. On the other hand, isotropic elasticity, being set in conjunction with problems of elasticity in the linear regime with [6] or without [7,8] body forces present, is an extensively developed area of continuum mechanics, embodying solid analysis that requires both analytical and numerical attention [9]. Even though the most attractive area for developing new methodologies appears to be the anisotropic linear elasticity [10,11], it is evident that there still exist open problems in the isotropic spectrum. In fact, the continuing interest in elaborating with such kind of aspects is still effective and the necessity in producing purely analytical solutions towards this direction, ready to accept proper numerical handling, stands in the frontline of the current research. Indeed, such a complete and comprehensive survey in linear isotropic elasticity brings insight to new elements of applied mathematics, which lead to the next step of studying the wave propagation [12,13] and, in general, the theory of scattering [14] in elastic materials.
Linear isotropic behavior of elastic media has an inherent mathematical interest due to the fact that even though the related theory is much simplified, many applications can accept the isotropic character without loss of robustness. In view of this aspect, this work is involved with the production of closed-type analytical formulae for the determination of the fundamental field in solid mechanics, i.e. the displacement. Doing so, we initially utilize the Hook's law, which provides the second order stress tensor in terms of the strain dyadic and the stiffness tetratic, while, in the sequel, we invoke the result into the Newton's law so as to obtain the generalized constitutive equation in elasticity. Thereafter, our purpose is twofold, that is we give special attention to elastostatics, neglecting the temporal derivatives and we exclude any external body forces, since any cause of disturbance can be entered into the boundary conditions of the particular physical problem. Moreover, we assume the specific expression for the stiffness tensor, which incorporates the appropriate components that inherit the isotropic character of the material. That way, we arrive at the known Navier partial differential equation for the displacement field, which accepts an analytical compact solution in the form of a partial differential operator, acting on harmonic functions with vector and scalar character, namely the Papkovich differential representation. This can be derived directly from the Naghdi-Hsu general solution [15] and since the spherical harmonics form a complete system, it also provides complete representations for the elastostatic fields. Hence, Papkovich representation of the solution of Navier equation is complete, that is any solution can be expressed in this form. The utility of such solutions is significant, considering the fact that they provide handy analytical expressions, for instance let us refer to a series of articles [16,17,18], which deal with elastic wave scattering at low frequencies around ellipsoidal solid bodies or cavities, using the Papkovich representation. Otherwise, the representation theory could be applied to inverse elastic scattering [19] or to more complicated mixed-type boundary value problems [20] in linear elasticity. Finally, let us not ignore a major advantage of Papkovich representation, which can offers us a certain degree of freedom, since the potentials needed to describe the field itself are redundant, thus the extra terms can be used conveniently, e.g. to compensate the force term in the case wherein it is included into the Navier equation.
It is the purpose of this research to provide progress in this interesting mathematical aspect, hence, using the spherical coordinate system [21], we obtain connection formulae, which relate the spherical harmonics that lead, through the Papkovich representation, to the displacement in linear isotropic elasticity. In other words, we calculate the displacement field, generated by the harmonic eigenfunctions [22], through the Papkovich representation, and then we provide the necessary analytical background so as to solve the inverse problem of identifying those harmonic eigenfunctions, which generate the same displacement fields. Henceforth, in the aim of solving sufficiently general interior and exterior problems, we use the internal and external solid spherical harmonic eigenfunctions in their complex form [22]. This way, we demonstrate the importance of attaining such kind of ready-to-use mathematical tools, by showing some interesting example problems inside and outside a sphere. The background of this mathematical procedure comes from the low-frequency scattering in elasticity [23] and shows how one can obtain a solution basis of a finite dimensional subspace. This reduces the calculation of the solution to the calculation of a finite number of scalar coefficients. However, our work leads to building blocks for constructing the solution of any relative problem.
Notwithstanding the existence of adequate and simultaneously convenient computational codes in solving problems in elasticity, we should not overlook the fact that pure analytical techniques are the backbone of numerical analysis. Hence, the important advantage of the performed mathematical analysis is based upon its ability to understand the physical background and to verify the credibility of numerical methods or other more sophisticated analytical models. On the other hand, the idea of building any solution from ready-to-use eigensolutions goes back to the classic references of Rayleigh [24], Kelvin [25], Maxwell [26], Sommerfeld [27] and Neuber [28]. In the present work, we tried to extend this idea to the theory of elasticity, wherein the ample literature survey of the fundamental classical references [29,30,31,32,33,34,35,36,37,38,39,40] verifies the necessity of such kind of analysis. Finally, by virtue of the representation theory, it is obvious that spherical geometry approximates sufficiently well most basic problems in linear isotropic elasticity. Nevertheless, the extension to spheroidal, ellipsoidal or even more complicated geometries [21,41] provides a challenging area for future investigation.
Let us introduce an arbitrarily defined smooth, either bounded or unbounded, three-dimensional elastic domain Ω(R3), which could be designated as interior or interior as the case may be. Then, each field within Ω(R3) is written in terms of its position vector r=x1ˆx1+x2ˆx2+x3ˆx3, expressed via the Cartesian basis ˆxj, j=1,2,3 in Cartesian coordinates (x1,x2,x3). On the other hand, the current investigation excludes the dependence on time, since we operate according to the steady state status of the particular situation, while the necessary fundamental information, which are adequate for this work can be found collected in [11].
The physical interpretation of linear isotropic elasticity is involved with the displacement field u, which comprises the measure of deformation of an elastic material. By means of the gradient ∇ and the Laplacian Δ operators, the displacement satisfies the well-known Navier equation in the presence of body forces, i.e.
μΔu(r)+(λ+μ)∇[∇⋅u(r)]+f(r)=0forr∈Ω(R3), | (1) |
wherein f is an external applied force that renders expression (1) non-homogeneous, while λ,μ∈R are the elastic parameters of the isotropic theory, being known as the Lamé constants. Papkovich proposed a differential representation of the solution for the homogeneous Navier equation [15] with no body forces, by considering f=0 within (1), which expresses the displacement field in differential form, in terms of two harmonic functions, one vector A and one scalar B, that is
u(r)=A(r)−λ+μ2(λ+2μ)∇[r⋅A(r)+B(r)],whereΔA(r)=0,ΔB(r)=0forr∈Ω(R3), | (2) |
which is known to be complete. It is not hard to prove that (2) satisfies (1), bearing in mind the vector identity Δ(r⋅A)=Δr⋅A+r⋅ΔA+2(∇⊗r)T:(∇⊗A)=2˜IT:(∇⊗A)=2˜I:(∇⊗A)=2∇⋅A and the interchange Δ∇=∇Δ, since it readily holds that Δr=ΔA=0 and ΔB=0, noting that ˜I is the unit dyadic, "T" denotes transposition, "⊗" refers to the tensor product and ":" stands for the double inner product. The general differential solution (2) provides a powerful analytical tool for solving the homogeneous and time-independent linearized equation of classical dynamic elasticity.
At this point we have to mention that even though representation (2) is complete, it is not unique, since we have four potentials (three components of the vector A and one component, which is the scalar B) to determine the three-dimensional vector displacement field. Actually, this is more than an advantage rather than a disadvantage, considering the fact that we are equipped with certain degrees of freedom to deploy conveniently the proper manner, which depends on the boundary value problem at hand. On the other end, this flexibility of the representation theory could be used to cancel the difficulty when the source term f is not absent (see equation (1) for instance). More precisely, in the same sense, we could use a similar differential solution like (2), in order to obtain a general solution for the non-homogeneous Navier equation (1). In fact, we can keep the same form for the displacement field, potential A could be harmonic, but then potential B should satisfy the mixed-type equation Δ(∇B)≡∇(ΔB)=2f/(λ+μ), in order for (1) to be satisfied. However, obtaining B is quite a difficult task.
This section includes the main results of our work, which are focused on the direct connection between the homogeneous Navier equation and the harmonic kernels of the implicated potentials A and B, using the spherical geometry. Before we proceed to the analysis in spherical coordinates, we take advantage of the already discussed flexibility of the differential general solution (2), thus and without loss of generality, we may suppose that the vector harmonic potential A is adequate enough to provide us with a complete solution for the displacement field u via (2). Consequently, assuming that B≡0, the representation (2) becomes
u(r)=A(r)−λ+μ2(λ+2μ)∇[r⋅A(r)],whereΔA(r)=0forr∈Ω(R3), | (3) |
which constitutes a complete solution for the Navier equation (1) if f=0. Taking the curl on both sides of (3) and defining a new arbitrary function Φ, we have
∇×u(r)=∇×A(r)⇒u(r)−A(r)=∇Φ(r)⇒A(r)=u(r)−∇Φ(r)forr∈Ω(R3). | (4) |
Next, we apply the div on (3), which leads us to
∇⋅u(r)=∇⋅A(r)−λ+μ2(λ+2μ)Δ[r⋅A(r)]=∇⋅A(r)−λ+μλ+2μ∇⋅A(r)=μλ+2μ∇⋅A(r)⇒∇⋅A(r)=λ+2μμ∇⋅u(r)forr∈Ω(R3). | (5) |
Next, we take the divergence of relation (4) and use simultaneously the outcome (5), to arrive at
∇⋅A(r)=∇⋅u(r)−ΔΦ(r)⇒ΔΦ(r)=−λ+μμ∇⋅u(r)forr∈Ω(R3), | (6) |
which is a Poisson equation with respect to the unknown function Φ, whose solution is the result of implying the fundamental solution of the Laplace's operator [20], obtaining
Φ(r)=λ+μ4πμ∭Ω(R3)∇r′⋅u(r′)|r−r′|dV(r′)forr∈Ω(R3). | (7) |
So, if we insert Φ from (7) into (4), we recover potential A as
A(r)=u(r)−λ+μ4πμ∇∭Ω(R3)∇r′⋅u(r′)|r−r′|dV(r′)forr∈Ω(R3). | (8) |
Therefore, given the displacement u, the harmonic potential A is given by (8). Recapitulating the above reasoning, we calculate the displacement field u, generated by the harmonic function A through the Papkovich representation (3) and then we face the inverse problem of determining this harmonic function A, which generates the displacement field u, via the integro-differential formula (8). This procedure is also invertible, in the sense that we can start with solutions of the form (8) and recover the displacement field via the Papkovich form (3). Both ways lead to the same result, independently of the coordinate system.
Aiming to approach applications that require an easily amenable expression for the displacement, we focus our attention to the spherical geometry [19], so in terms of the spherical coordinate system
x1=rsinθcosφ,x2=rsinθsinφandx3=rcosθ | (9) |
for 0⩽r<+∞, 0⩽θ⩽π and 0⩽φ<2π, the differential operators appearing into the Papkovich representation assume the forms
∇=ˆr∂∂r+ˆθˆθr∂∂θ+ˆφˆφrsinθ∂∂φandΔ=1r2∂∂r(r2∂∂r)+1r2sinθ∂∂θ(sinθ∂∂θ)+1r2sin2θ∂2∂φ2, | (10) |
where ˆr, ˆθˆθ, ˆφˆφ (see [19]) denote the coordinate vectors of the spherical system with position vector r=rˆr and unit dyadic ˜I=ˆr⊗ˆr+ˆθˆθ⊗ˆθˆθ+ˆφˆφ⊗ˆφˆφ. For every value of the nature n⩾0, there exist (2n+1) linearly independent complex spherical surface harmonics Ymn of degree n and of order |m|⩽n [20], given in the orthonormalized form by
Ymn(ˆr)=√2n+14π(n−|m|)!(n+|m|)!P|m|n(cosθ)eimφforˆr=(θ,φ)∈[0,π]×[0,2π) | (11) |
and with inner product
∫∫◯S2Ymn(ˆr)ˉYmn(ˆr)ds(ˆr)=1forˆr=(θ,φ)∈[0,π]×[0,2π), | (12) |
wherein P|m|n are the associated Legendre functions of the first kind [20], ˉYmn denote the complex conjugate surface spherical harmonics and S2 is the unit sphere in R3. Thereafter, we introduce the interior umn,in (regular as r→0+) and the exterior umn,ex (regular as r→+∞), solid spherical harmonic eigenfunctions for every n=0,1,2,... and m=−n,...,−1,0,1,...,n, which are given by the expression
umn,in(r)=rnYmn(ˆr)andumn,ex(r)=r−(n+1)Ymn(ˆr)forr∈Ω(R3), | (13) |
respectively. Relationships (13) comprise a complete set of eigenfunctions for harmonic functions and belong to the kernel space of the Laplace's operator from (10), i.e. Δumn,in=0 and Δumn,ex=0 for n⩾0 and |m|⩽n, while they are obtained once the classical method of separation of variables [21,22] is applied.
Adopting the above mathematical analysis, the harmonic function A in differential representation (3) admits series expansion in terms of functions (13), i.e.
A(r)=+∞∑n=0n∑m=−n[cmn,inumn,in(r)+cmn,exumn,ex(r)]forr∈Ω(R3), | (14) |
where
cmn,in=cm,1n,inˆx1+cm,2n,inˆx2+cm,3n,inˆx3andcmn,ex=cm,1n,exˆx1+cm,2n,exˆx2+cm,3n,exˆx3withn≥0and|m|≤n | (15) |
are arbitrary constant coefficients. Expansion (14) expresses the completeness of the interior and the exterior solid spherical harmonics. Consequently, substituting the potential (14) into the general solution (3), we obtain
u(r)=12(λ+2μ)+∞∑n=0n∑m=−n[(λ+3μ)cmn,inumn,in(r)−(λ+μ)(cmn,in⋅r)∇umn,in(r)+(λ+3μ)cmn,exumn,ex(r)−(λ+μ)(cmn,ex⋅r)∇umn,ex(r)]forr∈Ω(R3), | (16) |
where we used the trivial differential identities ∇(r⋅A)=A+(∇⊗A)⋅r (note that ∇⊗r=˜I) and ∇(cmn,yumn,y)=∇umn,y⊗cmn,y for y=in,ex, where n⩾0 and |m|⩽n. Nevertheless, expression (16) needs further elaboration in order for its processing to be feasible. To this direction, we analytically work as follows. Since the vector character of the harmonic function A is reflected upon the constant coefficients, which are written in Cartesian coordinates, we are obliged to work in the Cartesian system. This is attainable but requires the expression of the displacement field u in Cartesian coordinates involving constants and spherical surface harmonics. Therefore, we are able to transfer the connection between A and u to the corresponding connection via the constant coefficients. In order to do that, it is necessary to express the terms ∇umn,in and ∇umn,ex for n⩾0 and |m|⩽n as a function of spherical surface harmonics in Cartesian coordinates. This is possible, since these terms belong to the subspace that is produced by the spherical surface harmonics, while this task requires certain steps.
Hence, in the interest of making this work complete and independent, we provide recurrence relations for the associated Legendre functions [20], which, by definition of the conveniently chosen variable x=cosθ∈[−1,1], since θ∈[0,π], they are furnished by the Rodrigues formula
P|m|n(x)=(1−x2)|m|/212nn!d|m|+ndx|m|+n(x2−1)nforn≥0and|m|≤n, | (17) |
where obviously P|m|n(x)=0 if |m|>n. Thus, the associated Legendre functions of the first kind satisfy
(2n+1)xP|m|n(x)=(n+|m|)P|m|n−1(x)+(n−|m|+1)P|m|n+1(x), | (18) |
(2n+1)√1−x2P|m|n(x)=P|m|+1n+1(x)−P|m|+1n−1(x)=(n+|m|)(n+|m|−1)P|m|−1n−1(x)−(n−|m|+1)(n−|m|+2)P|m|−1n+1(x), | (19) |
2mx√1−x2P|m|n(x)=P|m|+1n(x)+(n+|m|)(n−|m|+1)P|m|−1n(x) | (20) |
and the first derivative relation
ddxP|m|n(x)=(1−x2)−1/2[−|m|x(1−x2)−1/2P|m|n(x)+P|m|+1n(x)], | (21) |
all relationships (18)–(21) being provided for every value of n⩾0 and |m|⩽n. On the other hand, the classical trigonometric functions imply the trivial identities
sinφsin|m|φ=12[cos(|m|−1)φ−cos(|m|+1)φ], | (22) |
cosφcos|m|φ=12[cos(|m|−1)φ+cos(|m|+1)φ], | (23) |
cosφsin|m|φ=12[sin(|m|+1)φ+sin(|m|−1)φ], | (24) |
sinφcos|m|φ=12[sin(|m|+1)φ−sin(|m|−1)φ] | (25) |
for any |m|⩽n (n⩾0). The derivation process to obtain handy expressions for ∇umn,in and ∇umn,ex with n⩾0 and |m|⩽n involves long and tedious calculations, however it is useful to provide the basic steps that gives us the necessary tools to accomplish this task. Under this aim, the procedure we follow is based on the utilization of the gradient operator within (10) in terms of the decomposition of the orthonormal spherical basis to the Cartesian one via the formulae
ˆr=sinθcosφˆx1+sinθsinφˆx2+cosθˆx3, | (26) |
ˆθˆθ=cosθcosφˆx1+cosθsinφˆx2−sinθˆx3, | (27) |
ˆφˆφ=−sinφˆx1+cosφˆx2. | (28) |
Thus, substituting above (26)–(28) into the spherical gradient operator from (10), we initially invoke the surface spherical harmonics from (11) into the corresponding solid spherical harmonics (13), next we write eimφ in the form e±i|m|φ=cos|m|φ±isin|m|φ for every |m|⩽n (n⩾0) and, finally, by virtue of cosθ=x and sinθ=√1−x2 with x∈[−1,1], we use recurrence relations (17)–(21) for the associated Legendre functions of the first kind, as well as relationships (22)–(25) for the trigonometric functions, so as to perform our analysis. Consequently, combining properly the relative terms to reproduce surface spherical harmonics, we derive for the interior solid spherical harmonic eigenfunctions umn,in in Ω(R3) the expressions
∇umn,in(r)=12κmn[(n+|m|)(n+|m|−1)Ym−1n−1(ˆr)κm−1n−1−Ym+1n−1(ˆr)κm+1n−1]rn−1ˆx1+i2κmn[(n+|m|)(n+|m|−1)Ym−1n−1(ˆr)κm−1n−1+Ym+1n−1(ˆr)κm+1n−1]rn−1ˆx2+κmn[(n+|m|)Ymn−1(ˆr)κmn−1+Ym+1n−1(ˆr)κm+1n−1]rn−1ˆx3forn≥0andm=−n,...,−1,1,...,n | (29) |
and for the case m=0
∇u0n,in(r)=−12κ0n[Y1n−1(ˆr)κ1n−1−Y−1n−1(ˆr)κ−1n−1]rn−1ˆx1+i2κ0n[Y1n−1(ˆr)κ1n−1+Y−1n−1(ˆr)κ−1n−1]rn−1ˆx2+κ0n[nY0n−1(ˆr)κ0n−1−Y−1n−1(ˆr)κ−1n−1]rn−1ˆx3forn≥0, | (30) |
while for the exterior solid spherical harmonic eigenfunctions umn,ex in Ω(R3), we similarly receive the expressions
∇umn,ex(r)=12κmn[(n−|m|+1)(n−|m|+2)Ym−1n+1(ˆr)κm−1n+1−Ym+1n+1(ˆr)κm+1n+1]r−(n+2)ˆx1+i2κmn[(n−|m|+1)(n−|m|+2)Ym−1n+1(ˆr)κm−1n+1+Ym+1n+1(ˆr)κm+1n+1]r−(n+2)ˆx2−κmn[(n−|m|+1)Ymn+1(ˆr)κmn+1+Ym+1n+1(ˆr)κm+1n+1]r−(n+2)ˆx3forn≥0andm=−n,...,−1,1,...,n | (31) |
and for the case m=0
∇u0n,ex(r)=−12κ0n[Y1n+1(ˆr)κ1n+1+Y−1n+1(ˆr)κ−1n+1]r−(n+2)ˆx1+i2κ0n[Y1n+1(ˆr)κ1n+1−Y−1n+1(ˆr)κ−1n+1]r−(n+2)ˆx2−κ0n[(n+1)Y0n+1(ˆr)κ0n+1]r−(n+2)ˆx3forn≥0, | (32) |
where
κmn=√2n+14π(n−|m|)!(n+|m|)!forn≥0and|m|≤n | (33) |
are the normalizing constants of the spherical surface harmonics. It is obvious, from the definition of the associated Legendre functions, that
Ym−n(ˆr)≡0forn≥0and|m|≤nwithˆr=(θ,φ)∈[0,π]×[0,2π), | (34) |
while
Ymn(ˆr)≡0forn≥0and|m|>nwithˆr=(θ,φ)∈[0,π]×[0,2π), | (35) |
rendering formulae (29)–(32) with (33) applicable for every value of n and m.
In order to validate the correctness and demonstrate the effectiveness of the produced formulae, we provide a simple example of evaluating ∇u11,in, utilizing two different approaches, one with direct calculation and the other with the aid of expression (29) for n=1 and m=1, showing that both the two results coincide. Towards this direction, from (11) and (33) we obtain
Y11(ˆr)=κ11P11(cosθ)eiφ=κ11sinθeiφwithκ11=√38πforˆr=(θ,φ)∈[0,π]×[0,2π), | (36) |
since P11(cosθ)=sinθ, therefore (13) yields
u11,in(r)=r1Y11(ˆr)=κ11rP11(cosθ)eiφ=κ11rsinθeiφforr∈Ω(R3). | (37) |
Thus, by direct action of the gradient operator from (10) on the interior harmonic (37) and in view of equations (26)–(28), it holds
∇u11,in(r)=[ˆr∂∂r+ˆθˆθr∂∂θ+ˆφˆφrsinθ∂∂φ](κ11rsinθeiφ)=κ11ˆrsinθeiφ+κ11ˆθˆθcosθeiφ+iκ11ˆφˆφeiφ=κ11(sinθcosφˆx1+sinθsinφˆx2+cosθˆx3)sinθeiφ+κ11(cosθcosφˆx1+cosθsinφˆx2−sinθˆx3)cosθeiφ+iκ11(−sinφˆx1+cosφˆx2)eiφ=κ11(cosφ−isinφ)ˆx1eiφ+κ11(sinφ+icosφ)ˆx2eiφforr∈Ω(R3) | (38) |
or
∇u11,in(r)=κ11ˆx1+iκ11ˆx2=κ11(ˆx1+iˆx2)forr∈Ω(R3) | (39) |
because eiφ=cosφ+isinφ. Otherwise, if we utilize the derived formula (29) for n=m=1 and we take profit from the property (35), yielding Y10≡0 and Y20≡0, we have
∇u11,in(r)=12κ11[2Y00(ˆr)κ00]ˆx1+i2κ11[2Y00(ˆr)κ00]ˆx2forr∈Ω(R3) | (40) |
or
∇u11,in(r)=κ11ˆx1+iκ11ˆx2=κ11(ˆx1+iˆx2)forr∈Ω(R3), | (41) |
since κ00=1 and Y00=1. Obviously, (39) and (41) are identical, not only validating the obtained expression (29), but also confirming that our approach is much more efficient and faster than the direct calculation of the gradient on solid spherical harmonics. It is evident that following a similar way, we can reproduce the subspace of any space of harmonic functions in the Cartesian basis, while a general procedure can be readily established.
Concluding, the result (16) for the displacement field, via the Papkovich differential representation, can be now rewritten in a handy form, by substituting the outcomes (29)–(32) with (33)–(35) into (16), using the Cartesian expressions for the implicated constant coefficients (15) and the interior, as well as the exterior harmonic eigenfunctions (13). Even though the current theory is applied on spherical boundaries, it can be extended to any surface. For instance, we indicate as an example that in a scattering problem by an arbitrary compact body we can always consider a sphere surrounding the scatterer. In this case, utilizing Green's second identity [22] we can transfer any information from the boundary of the scatterer to the sphere including the scatterer. This way, we reduce the exterior problem to the spherical geometry, where the boundary conditions are globally introduced. After this, it is obvious that the displacement field is now easily accessible to be applied to any kind of boundary value problem that concerns the wide area of linear isotropic elasticity.
The purpose of this section is to demonstrate the usefulness and the efficiency of the proposed analytical methodology, by invoking some special case problems in spherical geometry, which can be found in simple but quite important physical problems in linear isotropic elasticity. The specific choices of the examples come from specific boundary value problems appearing in the theory of low-frequency scattering [23]. In fact, the assumed elastic fields describe the leading low-frequency approximations of the incident excitation field. Consequently, the forthcoming examples come from real physical problems and they are not artificial.
We consider a spherical body of radius a, whose center coincides with the center of the Cartesian coordinate system. Therefore, the surrounding boundary ∂Ω≡S corresponds to the spherical variable at r=a for any θ∈[0,π] and φ∈[0,2π). This spherical surface separates the domain of interest Ω(R3)≡Ω into two subdomains, one interior Ω− for every r<a and one exterior Ω+ for every r>a, such as Ω=Ω−∪Ω+∪S. In order to facilitate our calculations, we define the ratio of the main phase velocities in an elastic medium [14] via the formula
τ=√μλ+2μ, | (42) |
hence, the differential representation of the displacement field (3) becomes
u−(r)=τ2+12A−(r)+τ2−12[∇⊗A−(r)]⋅r,whereΔA−(r)=0forr∈Ω− | (43) |
inside the sphere and
u+(r)=τ2+12A+(r)+τ2−12[∇⊗A+(r)]⋅r,whereΔA+(r)=0forr∈Ω+ | (44) |
outside the sphere, where we have once more utilized the identity ∇(r⋅A∓)=A∓+(∇⊗A∓)⋅r. The harmonic potentials A− and A+ admit expansions similar to (14) with (15), taking profit of their interior and exterior character, respectively, thus, using (13), we have
A−(r)=+∞∑n=0n∑m=−ncmn,inumn,in(r)=+∞∑n=0n∑m=−ncmn,inrnYmn(ˆr)forr∈Ω−, | (45) |
in order for the potential A− to be regular at the origin and
A+(r)=+∞∑n=0n∑m=−ncmn,exumn,ex(r)=+∞∑n=0n∑m=−ncmn,exr−(n+1)Ymn(ˆr)forr∈Ω+, | (46) |
in order for the potential A+ to remain bounded as we move towards infinity. The unknown vector constant coefficients
cmn,y=cm,1n,yˆx1+cm,2n,yˆx2+cm,3n,yˆx3fory=in,ex,wheren≥0and|m|≤n | (47) |
within (45) and (46) are calculated, when a particular set of boundary conditions is applied on the surface boundary S at r=a. In the sequel, we focus ourselves in providing general solutions of particular interest in Ω− and Ω+ for specific values of the degree n⩾0, without getting involved with the boundary value problem itself.
– 1st Example: If c is a constant vector, then it is not difficult to prove, performing direct substitution, that the Papkovich potentials
A−(r)=c,r<aandA+(r)=cr,r>a | (48) |
generate the displacement fields
u−(r)=τ2+12cforr∈Ω− | (49) |
and
u+(r)=[τ2+12˜I−τ2−12ˆr⊗ˆr]⋅crforr∈Ω+, | (50) |
respectively.
– 2nd Example: If ˜C is a constant tensor, let us introduce the Papkovich potentials
A−(r)=˜C⋅r,r<aandA+(r)=˜C⋅rr3,r>a, | (51) |
which, by virtue of the identities
∇⊗A−(r)=(∇⊗˜C)⋅r+(∇⊗r)⋅˜C⊺=˜C⊺ | (52) |
and
∇⊗A+(r)=(∇⊗˜C)⋅rr3+(∇⊗rr3)⋅˜C⊺=˜I−3ˆr⊗ˆrr3⋅˜C⊺ | (53) |
generate the displacement fields
u−(r)=(τ2+12˜C+τ2−12˜C⊺)⋅rforr∈Ω− | (54) |
and
u+(r)=[τ2+12˜C+τ2−12(˜I−3ˆr⊗ˆr)⋅˜C⊺]⋅rr3forr∈Ω+, | (55) |
respectively. Here, we must remark that, if we consider the potentials A−=r⋅˜C=˜C⊺⋅r and A+=(r/r3)⋅˜C=˜C⊺⋅(r/r3), we obtain the displacements generated by ˜C⊺ instead of ˜C. On the other hand, if ˜C is symmetric, i.e. if ˜C=˜C⊺, then the interior and exterior displacements (54) and (55) are rewritten as
u−(r)=τ2(˜C⋅r)forr∈Ω−andu+(r)=[τ2˜C−32(τ2−1)(˜C:ˆr⊗ˆr)˜I]⋅rr3forr∈Ω+, | (56) |
respectively. In particular, if ˜C=˜I, then (56) simplify to
u−(r)=τ2rforr∈Ω−andu+(r)=3−τ22rr3forr∈Ω+, | (57) |
respectively. Finally, if ˜C is antisymmetric, i.e. if ˜C⊺=−˜C, then obviously
˜C⊺=12(˜C⊺+˜C⊺)=12(˜C⊺−˜C) | (58) |
and since
˜C⊺:ˆr⊗ˆr=˜C:ˆr⊗ˆr, | (59) |
we obtain
˜C⊺:ˆr⊗ˆr=12˜C⊺:ˆr⊗ˆr−12˜C:ˆr⊗ˆr=0. | (60) |
Taking into account (58)–(60), the interior and exterior fields (54) and (55) are transformed to
u−(r)=˜C⋅rforr∈Ω−andu+(r)=˜C⋅rr3forr∈Ω+, | (61) |
respectively, which coincide with the harmonic potentials (51).
– 3rd Example: If ˜C is a constant dyadic, then we want to find the displacement field that the exterior harmonic potential
A+(r)=∇[˜C:∇⊗∇1r],r>awith˜C=3∑i=1ai⊗bi, | (62) |
produces. To prove it we first show that
∇1r=−rr3and∇⊗∇1r=∇⊗(−rr3)=3ˆr⊗ˆr−˜Ir3. | (63) |
Then, combining (62) and (63), we have
A+(r)=∇[˜C:3ˆr⊗ˆr−˜Ir3]=3∑i=1∇[ai⊗bi:(3r5r⊗r−1r3˜I)]=3∑i=1∇[3r5(ai⋅r)(bi⋅r)−1r3(ai⋅bi)]=3∑i=1[−15r7ai⊗bi:r⊗r⊗r+3r5ai⊗bi⋅r+3r5bi⊗ai⋅r+3r5(ai⋅bi)r]=3r53∑i=1[−5ai⊗bi:ˆr⊗ˆr⊗r+ai⊗bi⋅r+bi⊗ai⋅r+(ai⋅bi)r] | (64) |
or
A+(r)=3r5[(˜C:˜I)˜I+(˜C+˜C⊤)−5(˜C:ˆr⊗ˆr)˜I]⋅r,r>a. | (65) |
Thereafter, the gradient of the Papkovich potential (65), i.e. ∇⊗A+, is then calculated, using classical identities and algebra as follows,
∇⊗A+(r)=3∇[((˜C:˜I)˜I+(˜C+˜C⊤))⋅rr5]−153∑i=1∇[1r7(ai⋅r)(bi⋅r)r]=3(∇⊗rr5)⋅[(˜C:˜I)˜I+(˜C+˜C⊤)]⊤−153∑i=1[−7r5ai⊗bi:ˆr⊗ˆr⊗ˆr+1r5ai⊗bi⋅ˆr⊗ˆr+1r5bi⊗ai⋅ˆr⊗ˆr+1r5(ai⊗bi:ˆr⊗ˆr)˜I]=3[−5r5ˆr⊗ˆr+1r5˜I]⋅[(˜C:˜I)˜I+(˜C+˜C⊤)]+15[7r5(˜C:ˆr⊗ˆr)ˆr⊗ˆr−1r5(˜C+˜C⊤)⋅ˆr⊗ˆr−1r5(˜C:ˆr⊗ˆr)˜I],r>a. | (66) |
Contracting ∇⊗A+ of (66) with r from the right, we obtain
[∇A+(r)]⋅r=[−15r5ˆr⊗ˆr+3r5˜I]⋅[(˜C:˜I)r+(˜C+˜C⊤)⋅r]+15[7r5(˜C:ˆr⊗ˆr)r−1r5(˜C+˜C⊤)⋅r−1r5(˜C:ˆr⊗ˆr)r]=−15r5(˜C:˜I)r+3r5(˜C:˜I)r−15r5(˜C+˜C⊤):ˆr⊗ˆr⊗r+3r5(˜C+˜C⊤)⋅r+90r5(˜C:ˆr⊗ˆr)r−15r5(˜C+˜C⊤)⋅r=−12r5(˜C:˜I)r+60r5(˜C:ˆr⊗ˆr)r−12r5(˜C+˜C⊤)⋅r=−43r5[(˜C:˜I)˜I+(˜C+˜C⊤)−5(˜C:ˆr⊗ˆr)˜I]⋅r=−4A+(r),r>a. | (67) |
Substituting (67) into the Papkovich representation (44) for the exterior field, we obtain
u+(r)=τ2+12A+(r)+τ2−12[−4A+(r)]=5−3τ22A+(r)forr∈Ω+, | (68) |
which is the requested. Herein, we remark that, since A+ is harmonic and u+ solves the equation of elastostatics (3), we obtain that ∇∇⋅A+=0. Indeed,
∇∇⋅A+(r)=∇∇⋅[∇(˜C:∇⊗∇1r)]=∇Δ[˜C:∇⊗∇1r]=∇[˜C:∇⊗∇(Δ1r)]=0,r>a. | (69) |
Another interesting remark is that if ˜C is an antisymmetric dyadic, i.e. if ˜C⊺=−˜C, then it is not hard to prove
A+(r)=∇[˜C:∇∇1r]=0,r>aandthereforeu+(r)=0forr∈Ω+. | (70) |
This means that the potential A+ is symmetric and it generates a symmetric displacement. In particular, any ˜C generates a symmetric A+ and thus, a symmetric u+.
– 4th Example: In this case, we consider the equation of elastostatics (1) for f=0, which in terms of (42) is rewritten as
τ2Δu−(r)+(1−τ2)∇[∇⋅u−(r)]=0forr∈Ω− | (71) |
and we wish to find the coefficients α,β,γ for which the linear combination of the displacement field
u−(r)=α˜S:r⊗r⊗r+βr2˜S⋅r+γSr2rforr∈Ω− | (72) |
belongs to the ker[τ2Δ+(1−τ2)∇∇⋅], where ˜S=a⊗b=b⊗a=˜S⊺ is a constant symmetric dyadic and S=a⋅b is its trace. This example makes use of the previous three examples, as well as the linearity of the Navier equation, to expand the physical important interior displacement field in a three-dimensional solution subspace and to calculate the actual solution by calculating the three scalar coefficients. We primarily have to evaluate the following expressions,
Δ(˜S:r⊗r⊗r)=Δ[(a⋅r)(b⋅r)r]=Δ[(a⋅r)(b⋅r)]r+(a⋅r)(b⋅r)Δr+2∇[(a⋅r)(b⋅r)]r⋅∇⊗r=2(a⋅b)r+2(a⊗b+b⊗a)⋅r=2Sr+4˜S⋅r, | (73) |
Δ(r2˜S⋅r)=(Δr2)˜S⋅r+r2Δ(˜S⋅r)+2(∇r2)⋅∇(˜S⋅r)=6˜S⋅r+4r⋅˜S=10˜S⋅r, | (74) |
Δ(r2r)=6r+2⋅2r⋅˜I=10r, | (75) |
∇∇⋅(˜S:r⊗r⊗r)=∇∇⋅[(a⋅r)(b⋅r)r]=∇[(a⊗b+b⊗a):r⊗r+(a⋅r)(b⋅r)3]=5∇(˜S:r⊗r)=5(a⊗b+b⊗a)⋅r=10˜S⋅r, | (76) |
∇∇⋅(r2˜S×r)=∇[2r⋅˜S⋅r+r2∇⋅(˜S⋅r)]=∇[2(a⋅r)(b⋅r)+r2˜S:˜I]=2(a⊗b+b⊗a)⋅r+2Sr=4˜S⋅r+2Sr | (77) |
and
∇∇⋅(r2r)=∇[2r⋅r+3r2]=5∇r2=10r, | (78) |
wherein we have assumed that ˜S=(a⊗b+b⊗a)/2, so that ˜S:r⊗r=(a⋅r)(b⋅r) and S=a⋅b. In view of the above formulae (73)–(78) and with respect to (71) and (72), we now obtain
[τ2Δ+(1−τ2)∇∇⋅]u−(r)=ατ2[2Sr+4˜S⋅r]+βτ2[10˜S⋅r]+γτ2[10r]S+α(1−τ2)[10˜S⋅r]+β(1−τ2)[2Sr+4˜S⋅r]+γ(1−τ2)[10r]S=[2ατ2S+10γτ2S+2β(1−τ2)S+10γ(1−τ2)S]r+[4ατ2S+10βτ2+10α(1−τ2)+4β(1−τ2)]˜S⋅r=[2ατ2S+2β(1−τ2)S+10γS]r+[α(10−6τ2)+β(4+6τ2)]˜S⋅r=0forr∈Ω−, | (79) |
which implies
τ2Sα+(1−τ2)Sβ=−5γSand(5−3τ2)α+(2+3τ2)β=0, | (80) |
leading to the constants
α=3τ2+2,β=3τ2−5andγ=1−2τ2. | (81) |
Hence, the solution is provided via (72) by
u−(r)=(3τ2+2)˜S:r⊗r⊗r+(3τ2−5)r2˜S⋅r+(1−2τ2)Sr2rforr∈Ω−, | (82) |
which is the interior displacement field of the specific static elasticity problem.
Recapitulating, we offered the analytical solution for the displacement fields of four interesting problems in linear isotropic elastostatics in the absence of body forces, inside and outside a spherical boundary in the absence of boundary conditions. These solutions are solid and can be directly derived from the theory of differential representation, which we developed and analyzed in the previous section.
In this study, we presented a mathematical method for recovering the main spherical components of the well-known Navier equation in linear isotropic elastostatics, under the circumstance of no external forces present. To this end, we primarily combined the Hooke's and Newton's law, via the correlation of the displacement field with the strain, the stress and the stiffness tensors, reproducing the linearized equation of dynamic isotropic elasticity. In the sequel, we introduced the Papkovich differential representation, which offered solutions in terms of scalar and vector harmonic functions. Then, connection formulae were obtained, by which we transformed any solution of the Navier system from the Papkovich to the potential-type eigenform and vice versa. In order to enhance this procedure, we implied the commonly used spherical geometry and we calculated the time-independent displacement field, generated by the well-known spherical harmonic eigenfunctions. In order to demonstrate the effectiveness of our method, we presented some important degenerate cases for the evaluation of the interior and the exterior displacement field on either side of a spherical boundary.
Work under progress involves research directed towards the extension of the current analysis to more complicated geometries, e.g. spheroidal and ellipsoidal, producing ready-to-use functions and their Navier counterparts.
The authors have declared no conflict of interest.
[1] | A. A. Kilbas, H M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Amsterdam: Elsevier, 2006. |
[2] | D. Baleanu, Z. B. Guvenc, J. A. Tenreiro Machado, New trends in nanotechnology and fractional calculus applications, Springer, 2010. |
[3] |
M. Turkyilmazoglu, An efficient computational method for differential equations of fractional type, CMES-Comput. Model. Eng. Sci., 133 (2022), 1–9. http://dx.doi.org/10.32604/cmes.2022.020781 doi: 10.32604/cmes.2022.020781
![]() |
[4] |
A. Akgul, M. Inc, E. Karatas, D. Baleanu, Numerical solutions of fractional differential equations of Lane-Emden type by an accurate technique, Adv. Difference Equ., 2015 (2015), 220. http://doi.org/10.1186/s13662-015-0558-8 doi: 10.1186/s13662-015-0558-8
![]() |
[5] | B. Riemann, Versuch einer allgemeinen Auffassung der Integration und Differentiation, Cambridge: Cambridge University Press, 1847. |
[6] | M. Caputo, Elasticita e Dissipazione, Bologna: Zanichelli, 1969. |
[7] | I. Podlubny, Fractional differential equations, New York: Academic Press, 1998. |
[8] | K. S. Miller, B. Ross, An introduction to fractional calculus and fractional differential equations, Wiley, 1993. |
[9] | J. Liouville, Memoire sur quelques questions de geometrie et de mecanique, et sur un nouveaugenre de calcul pour resoudre ces questions, J. Ecole Polytech, 1832, 1–69. |
[10] |
W. Lin, Global existence theory and chaos control of fractional differential equations, J. Math. Anal. Appl., 332 (2007), 709–726. https://doi.org/10.1016/j.jmaa.2006.10.040 doi: 10.1016/j.jmaa.2006.10.040
![]() |
[11] |
N. A. Shah, H. A. Alyousef, S. A. El-Tantawy, J. D. Chung, Analytical investigation of fractional-order Korteweg-de-Vries-type equations under Atangana-Baleanu-Caputo operator: modeling nonlinear waves in a plasma and fluid, Symmetry, 14 (2022), 739. https://doi.org/10.3390/sym14040739 doi: 10.3390/sym14040739
![]() |
[12] |
M. M. Al-Sawalha, N. Amir, R. Shah, M. Yar, Novel analysis of fuzzy fractional Emden-Fowler equations within new iterative transform method, J. Funct. Spaces, 2022 (2022), 7731135. https://doi.org/10.1155/2022/7731135 doi: 10.1155/2022/7731135
![]() |
[13] |
P. Veeresha, D. G. Prakasha, H. M. Baskonus, Solving smoking epidemic model of fractional order using a modified homotopy analysis transform method, Math. Sci., 13 (2019), 115–128. https://doi.org/10.1007/S40096-019-0284-6 doi: 10.1007/S40096-019-0284-6
![]() |
[14] |
D. G. Prakasha, P. Veeresha, H. M. Baskonus, Analysis of the dynamics of hepatitis E virus using the Atangana-Baleanu fractional derivative, Eur. Phys. J. Plus, 134 (2019), 241. https://doi.org/10.1140/epjp/i2019-12590-5 doi: 10.1140/epjp/i2019-12590-5
![]() |
[15] |
M. M. Al-Sawalha, O. Y. Ababneh, R. Shah, K. Nonlaopon, Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators, AIMS Mathematics, 8 (2023), 2308–2336. https://doi.org/10.3934/math.2023120 doi: 10.3934/math.2023120
![]() |
[16] |
M. M. Al-Sawalha, R. P. Agarwal, O. Y. Ababneh, W. Weera, A reliable way to deal with fractional-order equations that describe the unsteady flow of a polytropic gas, Mathematics, 10 (2022), 2293. http://doi.org/10.3390/math10132293 doi: 10.3390/math10132293
![]() |
[17] |
D. Baleanu, G. C. Wu, S. D. Zeng, Chaos analysis and asymptotic stability of generalized Caputo fractional differential equations, Chaos Solitons Fractals, 102 (2017), 99–105. https://doi.org/10.1016/j.chaos.2017.02.007 doi: 10.1016/j.chaos.2017.02.007
![]() |
[18] | C. S. Drapaka, S. Sivaloganathan, A fractional model of continuum mechanics, J. Elasticity, 107 (2012), 105–123. |
[19] |
M. M. Al-Sawalha, A. S. Alshehry, K. Nonlaopon, R. Shah, O. Y. Ababneh, Fractional view analysis of delay differential equations via numerical method, AIMS Mathematics, 7 (2022), 20510–20523. https://doi.org/10.3934/math.20221123 doi: 10.3934/math.20221123
![]() |
[20] |
J. M. Cruz-Duarte, J. Rosales-Garcia, C. RodrigoCorrea-Cely, A. Garcia-Perez, J. GabrielAvina-Cervantes, A closed form expression for the Gaussian-based Caputo-Fabrizio fractional derivative for signal processing applications, Commun. Nonlinear Sci. Numer. Simul., 61 (2018), 138–148. https://doi.org/10.1016/j.cnsns.2018.01.020 doi: 10.1016/j.cnsns.2018.01.020
![]() |
[21] |
B. K. Singh, A novel approach for numeric study of 2D biological population model, Cogent Math., 3 (2016), 1261527. https://doi.org/10.1080/23311835.2016.1261527 doi: 10.1080/23311835.2016.1261527
![]() |
[22] |
E. Scalas, R. Gorenflo, F. Mainardi, Fractional calculus and continuous-time finance, Phys. A, 284 (2000), 376–384. https://doi.org/10.1016/S0378-4371(00)00255-7 doi: 10.1016/S0378-4371(00)00255-7
![]() |
[23] |
X. Gong, L. Wang, Y. Mou, H. Wang, X. Wei, W. Zheng, et al., Improved four-channel PBTDPA control strategy using force feedback bilateral teleoperation system, Int. J. Control, 20 (2022), 1002–1017. https://doi.org/10.1007/s12555-021-0096-y doi: 10.1007/s12555-021-0096-y
![]() |
[24] |
F. W. Meng, A. P. Pang, X. F. Dong, C. Han, X. P. Sha, H-infinity optimal performance design of an unstable plant under Bode integral constraint, Complexity, 2018 (2018), 4942906. https://doi.org/ 10.1155/2018/4942906 doi: 10.1155/2018/4942906
![]() |
[25] |
H. Zheng, S. F. Jin, A multi-source fluid queue based stochastic model of the probabilistic offloading strategy in a MEC system with multiple mobile devices and a single MEC server, Int. J. Appl. Math. Comput. Sci., 32 (2022), 125–138. https://doi.org/10.34768/amcs-2022-0010 doi: 10.34768/amcs-2022-0010
![]() |
[26] |
S. Momani, Z. Odibat, Analytical solution of a time-fractional Navier-Stokes equation by Adomian decomposition method, Appl. Math. Comput., 177 (2006), 488–494. https://doi.org/10.1016/j.amc.2005.11.025 doi: 10.1016/j.amc.2005.11.025
![]() |
[27] |
M. Alqhtani, K. M. Saad, R. Shah, W. Weera, W. M. Hamanah, Analysis of the fractional-order local poisson equation in fractal porous media, Symmetry, 14 (2022), 1323. https://doi.org/10.3390/sym14071323 doi: 10.3390/sym14071323
![]() |
[28] |
M. M. Al-Sawalha, A. S. Alshehry, K. Nonlaopon, R. Shah, O. Y. Ababneh, Approximate analytical solution of time-fractional vibration equation via reliable numerical algorithm, AIMS Mathematics, 7 (2022), 19739–19757. https://doi.org/10.3934/math.20221082 doi: 10.3934/math.20221082
![]() |
[29] |
G. C. Wu, A fractional variational iteration method for solving fractional nonlinear differential equations, Comput. Math. Appl., 61 (2011), 2186–2190. https://doi.org/10.1016/j.camwa.2010.09.010 doi: 10.1016/j.camwa.2010.09.010
![]() |
[30] |
D. G. Prakasha, P. Veeresha, M. S. Rawashdeh, Numerical solution for (2+1)-dimensional time-fractional coupled Burger equations using fractional natural decomposition method, Math. Methods Appl. Sci., 42 (2019), 3409–3427. https://doi.org/10.1002/mma.5533 doi: 10.1002/mma.5533
![]() |
[31] |
S. Mukhtar, R. Shah, S. Noor, The numerical investigation of a fractional-order multi-dimensional model of Navier-Stokes equation via novel techniques, Symmetry, 14 (2022), 1102. https://doi.org/10.3390/sym14061102 doi: 10.3390/sym14061102
![]() |
[32] | M. Areshi, A. Khan, R. Shah, K. Nonlaopon, Analytical investigation of fractional-order Newell-Whitehead-Segel equations via a novel transform. AIMS Mathematics, 7 (2022), 6936–6958. https://doi.org/10.3934/math.2022385 |
[33] |
E. M. Elsayed, R. Shah, K. Nonlaopon, The analysis of the fractional-order Navier-Stokes equations by a novel approach, J. Funct. Spaces, 2022 (2022), 8979447. https://doi.org/10.1155/2022/8979447 doi: 10.1155/2022/8979447
![]() |
[34] |
Z. H. Xie, X. A. Feng, X. J. Chen, Partial least trimmed squares regression, Chemometr. Intell. Lab. Syst., 221 (2022), 104486. https://doi.org/10.1016/j.chemolab.2021.104486 doi: 10.1016/j.chemolab.2021.104486
![]() |
[35] | V. N. Kovalnogov, M. I. Kornilova, Y. A. Khakhalev, D. A. Generalov, T. E. Simos, C. Tsitouras, Fitted modifications of Runge-Kutta-Nystrom pairs of orders 7(5) for addressing oscillatory problems, Math. Methods Appl. Sci., 2022. https://doi.org/10.1002/mma.8510 |
[36] | K. Liu, Z. X. Yang, W. F. Wei, B. Gao, D. L. Xin, C. M. Sun, et al., Novel detection approach for thermal defects: Study on its feasibility and application to vehicle cables, High Volt., 2022. https://doi.org/10.1049/hve2.12258 |
[37] |
B. A. Kupershmidt, A super Korteweg-de-Vries equations: an integrable system, Phys. Lett. A, 102 (1984), 213–215. https://doi.org/10.1016/0375-9601(84)90693-5 doi: 10.1016/0375-9601(84)90693-5
![]() |
[38] |
E. Fan, Uniformly constructing a series of explicit exact solutions to nonlinear equations in mathematical physics, Chaos Solitons Fractals, 16 (2003), 819–839. https://doi.org/10.1016/S0960-0779(02)00472-1 doi: 10.1016/S0960-0779(02)00472-1
![]() |
[39] |
M. Inc, On numerical soliton solution of the Kaup-Kupershmidt equation and convergence analysis of the decomposition method, Appl. Math. Comput., 172 (2006), 72–85. https://doi.org/10.1016/j.amc.2005.01.120 doi: 10.1016/j.amc.2005.01.120
![]() |
[40] | M. A. Ablowitz, P. A. Clarkson, Solitons, nonlinear evolution equations and inverse scattering, Cambridge: Cambridge University Press, 1991. |
[41] |
H. W. Tam, X. B. Hu, Two integrable differential-difference equations exhibiting soliton solutions of the Kaup-Kupershimdt equation type, Phys. Lett. A, 272 (2000), 174–183. https://doi.org/10.1016/S0375-9601(00)00422-9 doi: 10.1016/S0375-9601(00)00422-9
![]() |
[42] |
X. J. Yang, A new integral transform method for solving steady heat-transfer problem, Thermal Sci., 20 (2016), s639–s642. https://doi.org/10.2298/TSCI16S3639Y doi: 10.2298/TSCI16S3639Y
![]() |
[43] |
M. Caputo, M. Fabrizio, On the singular kernels for fractional derivatives. Some applications to partial differential equations, Prog. Fract. Differ. Appl., 7 (2021), 79–82. http://doi.org/10.18576/pfda/070201 doi: 10.18576/pfda/070201
![]() |
[44] |
S. Ahmad, A. Ullah, A. Akgul, M. De la Sen, A novel homotopy perturbation method with applications to nonlinear fractional order KdV and Burger equation with exponential-decay kernel, J. Funct. Spaces, 2021 (2021), 8770488. https://doi.org/10.1155/2021/8770488 doi: 10.1155/2021/8770488
![]() |
[45] |
D. G. Prakasha, N. S. Malagi, P. Veeresha, B. C. Prasannakumara, An efficient computational technique for time-fractional Kaup-Kupershmidt equation, Numer. Methods Partial Differential Equations, 37 (2021), 1299–1316. https://doi.org/10.1002/num.22580 doi: 10.1002/num.22580
![]() |
1. | D. Labropoulou, P. Vafeas, D. M. Manias, G. Dassios, Generalized Solutions in Isotropic and Anisotropic Elastostatics, 2025, 157, 0374-3535, 10.1007/s10659-025-10123-x |