
This study focuses on developing efficient numerical techniques for solving the fractional Keller-Segel (KS) model, which is critical in explaining chemotaxis events. Within the Caputo operator framework, the study applied two unique methodologies: The Aboodh residual power series method (ARPSM) and the Aboodh transform iteration method (ATIM). These approaches were used to find precise solutions to the fractional KS equation, resulting in a better understanding of chemotactic behavior in biological systems. The comparative examination of the ARPSM and ATIM revealed their distinct strengths and applications in solving complicated fractional models. The work advances numerical approaches for fractional differential equations and improves our understanding of chemotaxis dynamics using a precise modeling approach.
Citation: Nader Al-Rashidi. Innovative approaches to fractional modeling: Aboodh transform for the Keller-Segel equation[J]. AIMS Mathematics, 2024, 9(6): 14949-14981. doi: 10.3934/math.2024724
[1] | Giacomo Ascione, Enrica Pirozzi . On a stochastic neuronal model integrating correlated inputs. Mathematical Biosciences and Engineering, 2019, 16(5): 5206-5225. doi: 10.3934/mbe.2019260 |
[2] | Virginia Giorno, Serena Spina . On the return process with refractoriness for a non-homogeneous Ornstein-Uhlenbeck neuronal model. Mathematical Biosciences and Engineering, 2014, 11(2): 285-302. doi: 10.3934/mbe.2014.11.285 |
[3] | Meng Gao, Xiaohui Ai . A stochastic Gilpin-Ayala mutualism model driven by mean-reverting OU process with Lévy jumps. Mathematical Biosciences and Engineering, 2024, 21(3): 4117-4141. doi: 10.3934/mbe.2024182 |
[4] | Buyu Wen, Bing Liu, Qianqian Cui . Analysis of a stochastic SIB cholera model with saturation recovery rate and Ornstein-Uhlenbeck process. Mathematical Biosciences and Engineering, 2023, 20(7): 11644-11655. doi: 10.3934/mbe.2023517 |
[5] | Zhihao Ke, Chaoqun Xu . Structure analysis of the attracting sets for plankton models driven by bounded noises. Mathematical Biosciences and Engineering, 2023, 20(4): 6400-6421. doi: 10.3934/mbe.2023277 |
[6] | Katrine O. Bangsgaard, Morten Andersen, James G. Heaf, Johnny T. Ottesen . Bayesian parameter estimation for phosphate dynamics during hemodialysis. Mathematical Biosciences and Engineering, 2023, 20(3): 4455-4492. doi: 10.3934/mbe.2023207 |
[7] | Yongxin Gao, Shuyuan Yao . Persistence and extinction of a modified Leslie-Gower Holling-type Ⅱ predator-prey stochastic model in polluted environments with impulsive toxicant input. Mathematical Biosciences and Engineering, 2021, 18(4): 4894-4918. doi: 10.3934/mbe.2021249 |
[8] | Xiaochen Sheng, Weili Xiong . Soft sensor design based on phase partition ensemble of LSSVR models for nonlinear batch processes. Mathematical Biosciences and Engineering, 2020, 17(2): 1901-1921. doi: 10.3934/mbe.2020100 |
[9] | M. Nagy, M. H. Abu-Moussa, Adel Fahad Alrasheedi, A. Rabie . Expected Bayesian estimation for exponential model based on simple step stress with Type-I hybrid censored data. Mathematical Biosciences and Engineering, 2022, 19(10): 9773-9791. doi: 10.3934/mbe.2022455 |
[10] | Aniello Buonocore, Luigia Caputo, Enrica Pirozzi, Maria Francesca Carfora . A simple algorithm to generate firing times for leaky integrate-and-fire neuronal model. Mathematical Biosciences and Engineering, 2014, 11(1): 1-10. doi: 10.3934/mbe.2014.11.1 |
This study focuses on developing efficient numerical techniques for solving the fractional Keller-Segel (KS) model, which is critical in explaining chemotaxis events. Within the Caputo operator framework, the study applied two unique methodologies: The Aboodh residual power series method (ARPSM) and the Aboodh transform iteration method (ATIM). These approaches were used to find precise solutions to the fractional KS equation, resulting in a better understanding of chemotactic behavior in biological systems. The comparative examination of the ARPSM and ATIM revealed their distinct strengths and applications in solving complicated fractional models. The work advances numerical approaches for fractional differential equations and improves our understanding of chemotaxis dynamics using a precise modeling approach.
The smoothed particle hydrodynamics (SPH) method is a Lagrangian representation that is based on a particle interpolation to calculate smooth field variables [1,2]. Because of its Lagrangian features, the SPH method has fine advantages. The SPH method is used to simulate several industrial applications by simulating highly dynamic and violent flows such as tank sloshing, oil flow in gearboxes, nozzles, and jet impingement. In addition, the SPH method works well with the free surface flows, surface tension in free surface flows, and breaking waves without a special treatment required in this area. Additionally, the simulations of multiphase flows are easily handled by the SPH method. In conventional grid-based approaches, both the Volume-of-Fluid (VOF) method and the Level set method are adept at accurately capturing multiphase interfaces. An additional strength of the SPH method is its ability to handle complex/moving geometries. He et al. [3] introduced a new, weakly-compressible, SPH model tailored for multi-phase flows with notable differences in the density, while also allowing for large Courant-Friedrichs-Lewy numbers. Their method involves adjusting the continuity equation to exclude contributions from particles of different phases, solely focusing on neighboring particles of the same phase. Subsequently, they utilize a Shepard interpolation to reset the pressure and density of the particles from other phases. Bao et al. [4] presented an upgraded method within the SPH aimed at simulating turbulent flows near walls with medium to high Reynolds numbers. Their approach combined the κ−ε turbulence model and a wall function to enhance the accuracy. To tackle numerical errors and challenges related to the boundary conditions, they split the second-order partial derivative term of the composite function containing the turbulent viscosity coefficient into two components. Ensuring uniform particles, smooth pressure fields, and efficient computation, they employed particle shifting techniques, the δ -SPH method, and graphics processing unit (GPU) parallel processing in their simulations. Huang et al. [5] devised an innovative approach that merged SPH with the Finite Difference Method (SPH-FDM) within an Eulerian framework to analyze the unsteady flow around a pitching airfoil. They introduced an improved Iterative Shifting Particle Technology (IPST) to accurately track the motion of pitching airfoils. The numerical results from the SPH-FDM approach within a Eulerian framework closely aligned with those from either the finite volume method (FVM) or the existing literature, affirming the effectiveness of their proposed method. Further applications of the SPH method and its incompressible scheme, entitled the ISPH method, can be found in the references [6,7,8,9,10,11,12,13,14,15,16,17]. In recent years, the ISPH-FVM coupling method refers to the integration of the ISPH method with the FVM. This hybrid approach combines the strengths of both techniques to accurately model the fluid flow phenomena, particularly in scenarios involving either complex geometries or multiphase flows. Xu et al. [18] introduced the ISPH-FVM coupling method to simulate a two-phase incompressible flow, and combined the advantages of ISPH in interface tracking and FVM in flow field calculations with enhancements in surface tension calculations inspired by VOF and Level Set methods, demonstrating an improved accuracy and stability in complex two-phase flow scenarios. Xu et al. [19] developed a 3D ISPH-FVM coupling technique to model a bubble ascent in a viscous stagnant liquid, and employed an exchange of information at overlapping regions between the ISPH and FVM schemes, with an evaluation of surface tension effects using the 3D continuum surface force model (CSF). They validated the accuracy of the interface capture through benchmark tests and compared simulations of a single bubble ascent against experimental and literature-based data; moreover, they explored the method's capabilities in various parameter settings and in studying the dynamics of a two-bubble coalescence, demonstrating the reliability and versatility of the approach in simulating the bubble ascent phenomena. Xu et al. [20] proposed an advanced mapping interpolation ISPH-FVM coupling technique tailored to simulate complex two-phase flows characterized by intricate interfaces. The nanofluids were a suspension of nanomaterials in the base fluids. The nanofluids were applied to enhance the thermal features of cooling devices and equipment in electronic cooling systems [21,22] and heat exchangers [23,24]. In 1995, Choi and Eastman [25] achieved an enhanced heat transfer by adding nanoparticles inside a base fluid. The double diffusion of nanofluids (Al2O3-H2O) in an enclosure under the magnetic field impacts was studied by the lattice Boltzmann method [26]. There are several uses of coupling magnetic field impacts on nanofluids in thermal engineering experiments [27,28,29,30,31,32]. The variable magnetic field on ferrofluid flow in a cavity was studied numerically by the ISPH method [33]. Aly [16] adopted the ISPH method to examine the effects of Soret-Dufour numbers and magnetic fields on the double diffusion from grooves within a nanofluid-filled cavity. To enhance the heat transfer in cavities, the wavy or extended surface or fins were applied. The fins were used in several mechanical engineering applications. Saeid [34] studied the existence of different fin shapes on a heating cavity. The numerical attempt based on the finite volume method was introduced for natural convection and thermal radiation in a cavity containing N thin fins [35]. Hussain et al. [36] examined the effects of fins and inclined magnetic field on mixed convection of a nanofluid inside single/double lid-driven cavities. The recent numerical attempts based on the ISPH method for double diffusion of a nanofluid in a finned cavity, which were introduced by Aly et al. [37,38]. This work examines the impacts of Cattaneo-Christov heat flux and magnetic impacts on the double diffusion of a nanofluid-filled cavity. The cavity contains a rotated wavy circular cylinder and four fins fixed on its borders. The ISPH method is utilized to handle the interactions between the rotational speed of an inner wavy cylinder with a nanofluid flow. The established ANN model can accurately predict the values of −Nu and −Sh, as demonstrated in this study, in which there is a high agreement between the prediction values generated by the ANN model and the target values.
The initial physical model and its particle generation are represented in Figure 1. The physical model contains a square cavity filled with a nanofluid and containing four fins and a wavy circular cylinder. The embedded wavy circular cylinder rotates around its center by a rotational velocity V=ω(r−r0) with adiabatic thermal/solutal conditions. The four fins are fixed at the center of the cavity borders with a low temperature/concentration Tc&Cc. The source is located on the left wall of a cavity with Th&Ch, the plane walls are adiabatic, and the right wall is maintained at Tc&Cc. It is assumed that the nanofluid flow is laminar, unsteady, and incompressible. Figure 2 depicts the architecture and fundamental configuration structure of the created multilayer perceptron (MLP) network. A model with 14 neurons in the hidden layer of the MLP network model was identified so that the number of neurons were included in the hidden layer [39,40]. To build ANN models, the data-collecting process must be idealistically structured. Three methods were used to arrange the 1, 020, 888 data points that were used to construct the ANN model. In the created MLP network model, the −Nu and −Sh values were estimated in the output layer, while the τ and Sr values were defined as input parameters in the input layer. A total of 70% of the data was set aside for model training, 15% for validation, and 15% for testing. The recommended training algorithm in the network model was the Levenberg-Marquardt (LM) algorithm. The LM algorithm effectively minimizes nonlinear least squares problems by iteratively adjusting parameters based on an objective function, utilizing the gradient descent and Gauss-Newton methods, while dynamically controlling the step size to prevent overshooting; this achieves convergence to the optimal parameter values, making it widely used in various fields for an efficient optimization. Below are the transfer functions that are employed in the MLP network's hidden and output layers:
f(x)=21+e(−2x)−1 | (1) |
purelin(x)=x. | (2) |
The mathematical equations utilized to compute the mean squared error (MSE), coefficient of determination (R), and margin of deviation (MoD) values, which were selected as the performance criteria, are outlined below [41]:
MSE=1N∑Ni=1(Xtarg(i)−Xpred(i))2 | (3) |
R=√1−∑Ni=1(Xtarg(i)−Xpred(i))2∑Ni=1(Xtarg(i))2 | (4) |
MoD(%)=[Xtarg−XpredXtarg]x100. | (5) |
The regulating equations in the Lagrangian description [38,42,43] are as follows:
∂u∂x+∂v∂y=0, | (6) |
dudt=−1ρnf∂p∂x+μnfρnf(∂2u∂x2+∂2u∂y2)−σnfB20ρnf(usin2γ−vsinγcosγ), | (7) |
dvdt=−1ρnf∂p∂y+μnfρnf(∂2v∂x2+∂2v∂y2)+(ρβT)nfρnfg(T−Tc)+(ρβC)nfρnfg(C−Cc)−σnfB20ρnf(vcos2γ−usinγcosγ), | (8) |
dTdt=∇·(αnf∇T)+1(ρCP)nf∇·(D1∇C)−δ1(u∂T∂x∂u∂x+v∂T∂y∂v∂y+u2∂2T∂x2+v2∂2T∂y2+2uv∂2T∂x∂y+u∂T∂y∂v∂x+v∂T∂x∂u∂y), | (9) |
dCdt=∇·(Dm∇C)+∇·(D2∇T). | (10) |
The dimensionless quantities [38,42,43] are as follows:
X=xL,Y=yL,τ=tαfL2,U=uLαf,V=vLαf,αf=kf(ρCP)f |
P=pL2ρnfα2f,θ=T−TcTh−Tc,Φ=C−CcCh−Cc. | (11) |
After substituting Eq (6) into the dimension form of governing Eqs (1)-(5), then the dimensionless equations [38,42,43] are as follows:
∂U∂X+∂V∂Y=0, | (12) |
dUdτ=−∂P∂X+μnfρnfαf(∂2U∂X2+∂2U∂Y2)−σnfρfσfρnfHa2Pr(Usin2γ−Vsinγcosγ), | (13) |
dVdτ=−∂P∂Y+μnfρnfαf(∂2V∂X2+∂2V∂Y2)+(ρβ)nfρnfβfRaPr(θ+NΦ)−σnfρfσfρnfHa2Pr(Vcos2γ−Usinγcosγ), | (14) |
dθdτ=1(ρCP)nfDu(∂2Φ∂X2+∂2Φ∂Y2)+αnfαf(∂2θ∂X2+∂2θ∂Y2)−δC(U∂θ∂X∂U∂X+V∂θ∂Y∂V∂Y+U2∂2θ∂X2+V2∂2θ∂Y2+2UV∂2θ∂X∂Y+U∂θ∂Y∂V∂X+V∂θ∂X∂U∂Y), | (15) |
dΦdτ=Sr(∂2θ∂X2+∂2θ∂Y2)+1Le(∂2Φ∂X2+∂2Φ∂Y2). | (16) |
The Soret number is Sr=D2αf((Th−Tc)(Ch−Cc)), δC=νfδ1L2 is the Cattaneo Christov heat flux parameter, and Du=D1αf((Ch−Cc)(Th−Tc)) is Dufour number. The Hartmann number is Ha=√σfμfB0L, the buoyancy ratio is N=βC(Ch−Cc)βT(Th−Tc), Rayleigh's number is Ra=gβT(Th−Tc)L3νfαf, the Prandtl number is Pr=νfαf, and Lewis's number is Le=αfDm.
Left wall, θ=Φ=1,U=0=V, (17)
Right wall, θ=Φ=0,U=0=V,
Plane walls, ∂Φ∂Y=∂θ∂Y=0,U=0=V,
Four fins, θ=Φ=0,U=0=V,
Wavy circular cylinder ∂Φ∂Y=∂θ∂Y=0,U=Urot,V=Vrot,
The mean Nusselt/Sherwood numbers are as follows:
−Nu=−∫10knfkf(∂θ∂X)dY, | (18) |
−Sh=−∫10(∂Φ∂X)dY. | (19) |
The nanofluid properties are as follows ([44,45,46]):
ρnf=ρf−ϕρf+ϕρnp | (20) |
(ρCp)nf=(ρCp)f−ϕ(ρCp)f+ϕ(ρCp)np | (21) |
(ρβ)nf=ϕ(ρβ)np−ϕ(ρβ)f+(ρβ)f | (22) |
σnf=σf(1+3ϕ(σnp/σf−1)(σnp/σf+2)−(σnp/σf−1)ϕ) | (23) |
αnf=knf(ρCp)nf | (24) |
knf=(knp+2kf)+ϕ((knp+2kf)kf−2ϕ(kf−knp)kf)(kf−knp)−1 | (25) |
μnf=(1−ϕ)−2.5μf. | (26) |
The solving steps in the ISPH method [42,43,47,48] are as follows:
First:
U∗=μnfΔτρnfαf(∂2U∂X2+∂2U∂Y2)n−ΔτσnfρfσfρnfHa2Pr(Unsin2γ−Vnsinγcosγ)+Un | (27) |
V∗=μnfΔτρnfαf(∂2V∂X2+∂2V∂Y2)n−ΔτσnfρfσfρnfHa2Pr(Vncos2γ−Unsinγcosγ)+Δτ(ρβ)nfρnfβfRaPr(θn+NΦn)+Vn. | (28) |
Second:
∇2Pn+1=1Δτ(∂U∗∂X+∂V∗∂Y). | (29) |
Third:
Un+1=U∗−Δτ(∂P∂X)n+1, | (30) |
Vn+1=V∗−Δτ(∂P∂Y)n+1. | (31) |
Fourth:
θn+1=θn+αnfΔταf(∂2θ∂X2+∂2θ∂Y2)n+Δτ(ρCP)nfDu(∂2Φ∂X2+∂2Φ∂Y2)n−ΔτδC(U∂θ∂X∂U∂X+V∂θ∂Y∂V∂Y+U2∂2θ∂X2+V2∂2θ∂Y2+2UV∂2θ∂X∂Y+U∂θ∂Y∂V∂X+V∂θ∂X∂U∂Y)n, | (32) |
Φn+1=Φn+ΔτLe(∂2Φ∂X2+∂2Φ∂Y2)n+ΔτSr(∂2θ∂X2+∂2θ∂Y2)n. | (33) |
Fifth:
rn+1=ΔτUn+1+rn. | (34) |
The shifting technique [49,50] is as follows:
ϕi'=ϕi+(∇ϕ)i·δrii'+O(δr2ii'). | (35) |
The ISPH method's flowchart encompasses initiating the computational domain and boundary conditions, advancing through the time steps to estimate particle densities, solving the pressure fields, computing the heat and mass transfer, updating the particle velocities and positions, and ultimately analyzing results, all to accurately simulate the fluid dynamics. The flowchart of the ISPH method is shown in Figure 3.
This section introduces several numerical tests of the ISPH approach, beginning with an examination of the natural convection around a circular cylinder within a square cavity. Figure 4 illustrates a comparison of isotherms during the natural convection (NC) at different Rayleigh numbers (Ra) between the ISPH approach and the numerical findings by Kim et al. [51]. The ISPH approach yielded results consistent with those obtained by Kim et al. [51].
Second, the ISPH method's validation for the natural convection in a cavity with a heated rectangle source is presented. Figure 5 compares the ISPH method with numerical simulations (Fluent software) and experimental data from Paroncini and Corvaro [51] for isotherms at various Rayleigh numbers and hot source heights. The ISPH method shows good agreement with both experimental and numerical results across different scenarios. Additionally, Figure 6 illustrates comparisons between the experimental, numerical, and ISPH method results for streamlines at different Rayleigh numbers and hot source lengths, confirming the accuracy of the ISPH method simulations without redundancy.
The performed numerical simulations are discussed in this section. The range of pertinent parameters is the Cattaneo Christov heat flux parameter (0≤δc≤0.001), the Dufour number (0≤Du≤2), the nanoparticle parameter (0≤ϕ≤0.1), the Soret number (0≤Sr≤2), the Hartmann number (0≤Ha≤80), the Rayleigh number (103≤Ra≤105), Fin's length (0.05≤LFin≤0.2), and the radius of a wavy circular cylinder (0.05≤RCyld≤0.3). Figure 7 shows the effects of δc on the temperature θ, concentration Φ, and velocity field V. Increasing δc slightly enhances the temperature, concentration, and velocity field within an annulus. Figure 8 shows the effects of Du on θ, Φ, and V. Physically, the Dufour number indicates the influence of concentration gradients on thermal energy flux. An increase in Du has a minor influence on θ and Φ. Increasing Du growths the velocity field. These results are relevant to the presence of a diabatic wavy cylinder inside a square cavity, which shrinks the contributions of Du. Figure 9 represents the effects of Ha on θ, Φ, and V. The Hartmann number (Ha) is a proportion of an electromagnetic force to a viscous force. Physically, the Lorentz forces represent the magnetic field that shrinks the flow velocity. Consequently, an increment in Ha shrinks the distributions of θ and Φ. The maximum of the velocity field declines by 87.09% as Ha powers from 0 to 80. The study of the inner geometric parameters including the length of inner fins LFin and the radius of a wavy circular cylinder RCyld on the features of temperature, concentration, and velocity field is represented in Figures 10 and 11. In Figure 10, an expanded of LFin supports the cooling area, which reduces the strength of θ and Φ within a cavity. As the inner fins represent blockages for the nanofluid flow in a square cavity, the maximum of the velocity field reduces by 48.65% as the LFin boosts from 0.05 to 0.2. In Figure 11, regardless of the adiabatic condition of an inner wavy cylinder, an increment in RCyld reduces the distributions of temperature and concentration across the cavity. Furthermore, due to the low rotating speed of the inner wavy cylinder, an expanded wavy cylinder reduces the velocity field. It is seen that the maximum of the velocity field shrinks by 55.42% according to an increase in RCyld from 0.05 to 0.3. Figure 12 represents the effects of ϕ on θ, Φ, and V. Here, due to the presence of an inner wavy cylinder, the contribution of ϕ is minor in enhancing the cooling area and reducing the strength of the concentration inside a cavity, whilst the physical contribution of ϕ in increasing the viscosity of a nanofluid appears well by majoring the velocity field. The velocity's maximum reduces by 36.52% according to an increasing concentration of nanoparticles from 0 to 10%. Figure 13 establishes the effects of Ra on θ, Φ, and V. As Ra represents the ratio of buoyancy to thermal diffusivity and it is considering a main factor in controlling the natural convection. Accordingly, there are beneficial enhancements in the strength of the temperature and concentration in a cavity according to an increase in Ra. Furthermore, the Rayleigh number effectively works to accelerate the velocity field in a cavity. Figure 14 indicates the effects of Sr on θ, Φ, and V. Physically, the Soret influence expresses diffusive movement that is created from a temperature gradient. Increasing Sr enhances the strength of the temperature; this increment strongly enhances the strength of the concentration. Increasing Sr from 0 to 2 supports the velocity's maximum by 134.3%. Figure 15 indicates the average −Nu and −Sh under the variations of Sr, RCyld, Ha, and δC. The average −Nu is enhanced according to an increment in Sr, whilst it decreases alongside an increase in RCyld and Ha. The average −Sh decreases alongside an increase of Sr, RCyld, and Ha. It is found that the variations of δC slightly changes the values of −Nu and −Sh.
Figure 16 displays the performance chart that was made to examine the produced ANN model's training phase. The graph illustrates how the MSE values, which were large at the start of the training phase, got smaller as each epoch progressed. By achieving the optimal validation value for each of the three data sets, the MLP model's training phase was finished. The error histogram derived from the training phase's data is displayed in Figure 17. Upon scrutinizing the error histogram, one observes that the error values are predominantly centered around the zero-error line. Additionally, relatively low numerical values for the errors were achieved. In Figure 18, the regression curve for the MLP model is displayed, indicating that R = 1 is consistently employed for training, validation, and testing. According to the results of the performance, regression, and error histograms, the built ANN model's training phase was successfully finished. Following the training phase's validation, the ANN model's prediction performance was examined. First, the goal values and the −Nu and −Sh values derived by the ANN model were compared in this context. The goal values for every data point are displayed in Figure 19, along with the −Nu and −Sh values that the ANN model produced. Upon a closer inspection of the graphs, it becomes evident that the −Nu and −Sh values derived by the ANN model perfectly align with the desired values.
The current work numerically examined the effects of a Cattaneo-Christov heat flux on the thermosolutal convection of a nanofluid-filled square cavity. The square cavity had four fins located on the cavity's boarders and a wavy circular cylinder, which was located in the cavity's center, and rotated in a circular form. This study checked the impacts of the Cattaneo Christov heat flux parameter, the Dufour number, a nanoparticle parameter, the Soret number, the Hartmann number, the Rayleigh number, Fin's length, and the radius of a wavy circular cylinder on the temperature, concentration, and velocity field. The results are highlighted within the following points:
• The parameters of the embedded blockages, including the fins length and radius of a wavy circular cylinder, worked well to improve the cooling area and reduce the strength of concentration as well as the velocity field.
• Adding an extra concentration of nanoparticles up to 10% enhanced the cooling area and reduced the maximum velocity field by 36.52%.
• The Rayleigh number beneficially contributed to power the double diffusion within an annulus and was effective in accelerating the nanofluid velocity.
• The Soret number effectively contributed to change the distributions of the temperature, the concentration, and the velocity field within an annulus as compared to the Dufour number.
• The excellent accord between the prediction values produced by the ANN model and the goal values indicated that the developed ANN model could forecast the −Nu and −Sh values with a remarkable precision.
Munirah Alotaibi: Conceptualization, Funding acquisition, Methodology, Writing – original draft; Abdelraheem M. Aly: Formal Analysis, Investigation, Validation, Visualization, Writing – review & editing. 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 extend their appreciation to the Deanship of Scientific Research at King Khalid University, Abha, Saudi Arabia, for funding this work through the Research Group Project under Grant Number (RGP. 2/38/45). This research was funded by the Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R522), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.
The authors declare that they have NO affiliations with or involvement in any organization or entity with any financial interest.
[1] | A. Carpinteri, F. Mainardi, Fractals and fractional calculus in continuum mechanics, Vienna: Springer, https://doi.org/10.1007/978-3-7091-2664-6 |
[2] |
L. Debnath, Recent applications of fractional calculus to science and engineering, Int. J. Math. Math. Sci., 2003 (2003), 753601. https://doi.org/10.1155/S0161171203301486 doi: 10.1155/S0161171203301486
![]() |
[3] |
J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci. Numer. Simul., 16 (2011), 1140–1153. https://doi.org/10.1016/j.cnsns.2010.05.027 doi: 10.1016/j.cnsns.2010.05.027
![]() |
[4] |
A. A. M. Arafa, S. Z. Rida, H. Mohamed, Approximate analytical solutions of Schnakenberg systems by homotopy analysis method, Appl. Math. Model., 36 (2012), 4789–4796. https://doi.org/10.1016/j.apm.2011.12.014 doi: 10.1016/j.apm.2011.12.014
![]() |
[5] |
P. Sunthrayuth, A. M. Zidan, S. W. Yao, R. Shah, M. Inc, The comparative study for solving fractional-order Fornberg Whitham equation via ρ-Laplace transform, Symmetry, 13 (2021), 784. https://doi.org/10.3390/sym13050784 doi: 10.3390/sym13050784
![]() |
[6] |
R. Shah, H. Khan, D. Baleanu, Fractional Whitham Broer Kaup equations within modified analytical approaches, Axioms, 8 (2019), 125. https://doi.org/10.3390/axioms8040125 doi: 10.3390/axioms8040125
![]() |
[7] |
H. M. Srivastava, R. Shah, H. Khan, M. Arif, Some analytical and numerical investigation of a family of fractional-order Helmholtz equations in two space dimensions, Math. Methods Appl. Sci., 43 (2020), 199–212. https://doi.org/10.1002/mma.5846 doi: 10.1002/mma.5846
![]() |
[8] |
H. Yasmin, N. H. Aljahdaly, A. M. Saeed, Investigating symmetric soliton solutions for the fractional coupled konno onno system using improved versions of a novel analytical technique, Mathematics, 11 (2023), 2686. https://doi.org/10.3390/math11122686 doi: 10.3390/math11122686
![]() |
[9] |
M. M. Al-Sawalha, R. Shah, A. Khan, O. Y. Ababneh, T. Botmart, Fractional view analysis of Kersten-Krasil'shchik coupled KdV-mKdV systems with non-singular kernel derivatives, AIMS Mathematics, 7 (2022), 18334-18359. https://doi.org/10.3934/math.20221010 doi: 10.3934/math.20221010
![]() |
[10] |
A. A. Alderremy, R. Shah, N. Iqbal, S. Aly, K. Nonlaopon, Fractional series solution construction for nonlinear fractional reaction-diffusion Brusselator model utilizing Laplace residual power series, Symmetry, 14 (2022), 1944. https://doi.org/10.3390/sym14091944 doi: 10.3390/sym14091944
![]() |
[11] |
S. Alshammari, M. M. Al-Sawalha, R. Shah, Approximate analytical methods for a fractional-order nonlinear system of Jaulent Miodek equation with energy-dependent Schrodinger potential, Fractal Fract., 7 (2023), 140. https://doi.org/10.3390/fractalfract7020140 doi: 10.3390/fractalfract7020140
![]() |
[12] |
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
![]() |
[13] |
M. Alqhtani, K. M. Saad, 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
![]() |
[14] |
H. Yasmin, A. S. Alshehry, A. H. Ganie, A. M. Mahnashi, Perturbed Gerdjikov Ivanov equation: Soliton solutions via Backlund transformation, Optik, 298 (2024), 171576. http://dx.doi.org/10.1016/j.ijleo.2023.171576 doi: 10.1016/j.ijleo.2023.171576
![]() |
[15] |
N. A. Pirim, F. Ayaz, A new technique for solving fractional order systems: Hermite collocation method, Appl. Math., 7 (2016), 2307–2323. http://dx.doi.org/10.4236/am.2016.718182 doi: 10.4236/am.2016.718182
![]() |
[16] |
V. Marinca, N. Herisanu, Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer, Int. Commun. Heat Mass Transfer, 35 (2008), 710–715. https://doi.org/10.1016/j.icheatmasstransfer.2008.02.010 doi: 10.1016/j.icheatmasstransfer.2008.02.010
![]() |
[17] | J. S. Duan, R. Rach, D. Baleanu, A. M. Wazwaz, A review of the Adomian decomposition method and its applications to fractional differential equations, Commun. Frac. Calc., 3 (2012), 73–99. |
[18] |
M. Khan, M. A. Gondal, I. Hussain, S. K. Vanani, A new comparative study between homotopy analysis transform method and homotopy perturbation transform method on a semi infinite domain, Math. Comput. Model., 55 (2012), 1143–1150. https://doi.org/10.1016/j.mcm.2011.09.038 doi: 10.1016/j.mcm.2011.09.038
![]() |
[19] |
A. Jabbari, H. Kheiri, A. Yildirim, Homotopy analysis and homotopy Pade methods for (1+1) and (2+1) dimensional dispersive long wave equations, Internat. J. Numer. Methods Heat Fluid Flow, 23 (2013), 692–706. http://dx.doi.org/10.1108/09615531311323818 doi: 10.1108/09615531311323818
![]() |
[20] |
R. K. Gazizov, A. A. Kasatkin, Construction of exact solutions for fractional order differential equations by the invariant subspace method, Comput. Math. Appl., 66 (2013), 576–584. https://doi.org/10.1016/j.camwa.2013.05.006 doi: 10.1016/j.camwa.2013.05.006
![]() |
[21] |
A. Prakash, P. Veeresha, D. G. Prakasha, M. Goyal, A new efficient technique for solving fractional coupled Navier-Stokes equations using q-homotopy analysis transform method, Pramana-J. Phys., 93 (2018), 6. https://doi.org/10.1007/s12043-019-1763-x doi: 10.1007/s12043-019-1763-x
![]() |
[22] |
R. K. Pandey, H. K. Mishra, Homotopy analysis Sumudu transform method for time-fractional third order dispersive partial differential equation, Adv. Comput. Math., 43 (2017), 365–383. https://doi.org/10.1007/s10444-016-9489-5 doi: 10.1007/s10444-016-9489-5
![]() |
[23] |
Z. H. Guo, O. Acan, S. Kumar, Sumudu transform series expansion method for solving the local fractional Laplace equation in fractal thermal problems, Thermal Sci., 20 (2016), 739–742. http://dx.doi.org/10.2298/TSCI16S3739G doi: 10.2298/TSCI16S3739G
![]() |
[24] |
K. K. Ali, M. Maneea, M. S. Mohamed, Solving nonlinear fractional models in superconductivity using the q-Homotopy analysis transform method, J. Math., 2023 (2023), 6647375. https://doi.org/10.1155/2023/6647375 doi: 10.1155/2023/6647375
![]() |
[25] |
Z. Y. Fan, K. K. Ali, M. Maneea, M. Inc, S. W. Yao, Solution of time fractional Fitzhugh-Nagumo equation using semi analytical techniques, Results Phys., 51 (2023), 106679. https://doi.org/10.1016/j.rinp.2023.106679 doi: 10.1016/j.rinp.2023.106679
![]() |
[26] |
K. K. Ali, F. E. A. Elbary, M. Maneea, Efficient techniques for nonlinear dynamics: A study of fractional generalized quintic Ginzburg-Landau equation, J. Taibah Univ. Sci., 18 (2024), 2333593. https://doi.org/10.1080/16583655.2024.2333593 doi: 10.1080/16583655.2024.2333593
![]() |
[27] | M. A. El-Tawil, S. N. Huseen, The q-homotopy analysis method (q-HAM), Int. J. Appl. Math. Mech., 8 (2012), 51–75. |
[28] | M. A. El-Tawil, S. N. Huseen, On convergence of the q-homotopy analysis method, Int. J. Contemp. Math. Sci., 8 (2013), 481–497. |
[29] |
Z. J. Liu, M. Y. Adamu, E. Suleiman, J. H. He, Hybridization of homotopy perturbation method and Laplace transformation for the partial differential equations, Thermal Sci., 21 (2017), 1843–1846. http://dx.doi.org/10.2298/TSCI160715078L doi: 10.2298/TSCI160715078L
![]() |
[30] | A. Prakash, H. Kaur, q-homotopy analysis transform method for space and time-fractional KdV-Burgers equation, Nonlinear Sci. Lett. A, 9 (2018), 44–61. |
[31] |
E. F. Keller, L. A. Segel, Initiation of slime mold aggregation viewed as an instability, J. Theoret. Biol., 26 (1970), 399–415. https://doi.org/10.1016/0022-5193(70)90092-5 doi: 10.1016/0022-5193(70)90092-5
![]() |
[32] |
A. Atangana, Extension of the Sumudu homotopy perturbation method to an attractor for one-dimensional Keller-Segel equations, Appl. Math. Model., 39 (2015), 2909–2916. https://doi.org/10.1016/j.apm.2014.09.029 doi: 10.1016/j.apm.2014.09.029
![]() |
[33] |
A. Atangana, B. S. T. Alkahtani, Analysis of the Keller-Segel model with a fractional derivative without singular kernel, Entropy, 17 (2015), 4439–4453. https://doi.org/10.3390/e17064439 doi: 10.3390/e17064439
![]() |
[34] |
A. Atangana, E. Alabaraoye, Solving a system of fractional partial differential equations arising in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller-Segel equations, Adv. Differ. Equ., 2013 (2013), 94. https://doi.org/10.1186/1687-1847-2013-94 doi: 10.1186/1687-1847-2013-94
![]() |
[35] |
M. Zayernouri, A. Matzavinos, Fractional Adams-Bashforth/Moulton methods: An application to the fractional Keller-Segel chemotaxis system, J. Comput. Phys., 317 (2016), 1–14. https://doi.org/10.1016/j.jcp.2016.04.041 doi: 10.1016/j.jcp.2016.04.041
![]() |
[36] |
S. Kumar, A. Kumar, I. K. Argyros, A new analysis for the Keller-Segel model of fractional order, Numer. Algorithms, 75 (2017), 213–228. https://doi.org/10.1007/s11075-016-0202-z doi: 10.1007/s11075-016-0202-z
![]() |
[37] |
M. A. Dokuyucu, D. Baleanu, E. Çelik, Analysis of Keller-Segel model with Atangana-Baleanu fractional derivative, Filomat, 32 (2018), 5633–5643. http://dx.doi.org/10.2298/FIL1816633D doi: 10.2298/FIL1816633D
![]() |
[38] |
X. Luo, M. Nadeem, M. Inc, S. Dawood, Fractional complex transform and homotopy perturbation method for the approximate solution of Keller-Segel model, J. Funct. Spaces, 2022 (2022), 9637098. https://doi.org/10.1155/2022/9637098 doi: 10.1155/2022/9637098
![]() |
[39] |
O. A. Arqub, Series solution of fuzzy differential equations under strongly generalized differentiability, J. Adv. Res. Appl. Math., 5 (2013), 31–52. http://dx.doi.org/10.5373/jaram.1447.051912 doi: 10.5373/jaram.1447.051912
![]() |
[40] |
O. A. Arqub, Z. Abo-Hammour, R. Al-Badarneh, S. Momani, A reliable analytical method for solving higher-order initial value problems, Discrete Dyn. Nat. Soc., 2013 (2013), 673829. http://dx.doi.org/10.1155/2013/673829 doi: 10.1155/2013/673829
![]() |
[41] |
O. A. Arqub, A. El-Ajou, Z. A. Zhour, S. Momani, Multiple solutions of nonlinear boundary value problems of fractional order: A new analytic iterative technique, Entropy, 16 (2014), 471–493. https://doi.org/10.3390/e16010471 doi: 10.3390/e16010471
![]() |
[42] |
A. El-Ajou, O. A. Arqub, S. Momani, Approximate analytical solution of the nonlinear fractional KdV-Burgers equation: A new iterative algorithm, J. Comput. Phys., 293 (2015), 81–95. https://doi.org/10.1016/j.jcp.2014.08.004 doi: 10.1016/j.jcp.2014.08.004
![]() |
[43] |
S. Rida, A. Arafa, A. Abedl-Rady, H. Abdl-Rahaim, Fractional physical differential equations via natural transform, Chinese J. Phys., 55 (2017), 1569–1575. https://doi.org/10.1016/j.cjph.2017.05.004 doi: 10.1016/j.cjph.2017.05.004
![]() |
[44] |
J. Zhang, Z. Wei, L. Li, C. Zhou, Least-squares residual power series method for the time-fractional differential equations, Complexity, 2019 (2019), 6159024. https://doi.org/10.1155/2019/6159024 doi: 10.1155/2019/6159024
![]() |
[45] |
Y. Xie, I. Ahmad, T. I. S. Ikpe, E. F. Sofia, H. Seno, What influence could the acceptance of visitors cause on the epidemic dynamics of a Reinfectious disease?: A mathematical model, Acta Biotheor., 72 (2024), 3. https://doi.org/10.1007/s10441-024-09478-w doi: 10.1007/s10441-024-09478-w
![]() |
[46] |
I. Jaradat, M. Alquran, K. Al-Khaled, An analytical study of physical models with inherited temporal and spatial memory, Eur. Phys. J. Plus, 133 (2018), 162. https://doi.org/10.1140/epjp/i2018-12007-1 doi: 10.1140/epjp/i2018-12007-1
![]() |
[47] | M. Alquran, K. Al-Khaled, S. Sivasundaram, H. M. Jaradat, Mathematical and numerical study of existence of bifurcations of the generalized fractional Burgers-Huxley equation, Nonlinear Stud., 24 (2017), 235–244. |
[48] |
I. Ahmad, H. Seno, An epidemic dynamics model with limited isolation capacity, Theory Biosci., 142 (2023), 259–273. https://doi.org/10.1007/s12064-023-00399-9 doi: 10.1007/s12064-023-00399-9
![]() |
[49] |
G. O. Ojo, N. I. Mahmudov, Aboodh transform iterative method for spatial diffusion of a biological population with fractional-order, Mathematics, 9 (2021), 155. https://doi.org/10.3390/math9020155 doi: 10.3390/math9020155
![]() |
[50] |
M. A. Awuya, G. O. Ojo, N. I. Mahmudov, Solution of space-time fractional differential equations using Aboodh transform iterative method, J. Math., 2022 (2022), 4861588. https://doi.org/10.1155/2022/4861588 doi: 10.1155/2022/4861588
![]() |
[51] |
M. A. Awuya, D. Subasi, Aboodh transform iterative method for solving fractional partial differential equation with Mittag-Leffler Kernel, Symmetry, 13 (2021), 2055. https://doi.org/10.3390/sym13112055 doi: 10.3390/sym13112055
![]() |
[52] |
M. I. Liaqat, S. Etemad, S. Rezapour, C. Park, A novel analytical Aboodh residual power series method for solving linear and nonlinear time-fractional partial differential equations with variable coefficients, AIMS Mathematics, 7 (2022), 16917–16948. http://dx.doi.org/10.3934/math.2022929 doi: 10.3934/math.2022929
![]() |
[53] |
M. I. Liaqat, A. Akgul, H. Abu-Zinadah, Analytical investigation of some time-fractional Black-Scholes models by the Aboodh residual power series method, Mathematics, 11 (2023), 276. https://doi.org/10.3390/math11020276 doi: 10.3390/math11020276
![]() |
[54] | K. S. Aboodh, The new integral transform'Aboodh transform, Glob. J. Pure Appl. Math., 9 (2013), 35–43. |
[55] | S. Aggarwal, R. Chauhan, A comparative study of Mohand and Aboodh transforms, Int. J. Res. Adv. Technol., 7 (2019), 520–529. |
[56] |
M. E. Benattia, K. Belghaba, Application of the Aboodh transform for solving fractional delay differential equations, Univ. J. Math. Appl., 3 (2020), 93–101. https://doi.org/10.32323/ujma.702033 doi: 10.32323/ujma.702033
![]() |
[57] |
B. B. Delgado, J. E. Macias-Diaz, On the general solutions of some non-homogeneous Div-curl systems with Riemann-Liouville and Caputo fractional derivatives, Fractal Fract., 5 (2021), 117. https://doi.org/10.3390/fractalfract5030117 doi: 10.3390/fractalfract5030117
![]() |
[58] |
S. Alshammari, M. Al-Smadi, I. Hashim, M. A. Alias, Residual power series technique for simulating fractional Bagley-Torvik problems emerging in applied physics, Appl. Sci., 9 (2019), 5029. https://doi.org/10.3390/app9235029 doi: 10.3390/app9235029
![]() |