Abstract
Elzaki decomposition method (EDM) is adopted to deal with the fractional-order relaxation and damped oscillation equation along with time-fractional spatial diffusion biological population model in different suitable habitat situations. In accordance with the graphs for the solutions obtained, the fractional relaxation exhibits super-slow phenomenon due to its extended descent, and fractional damped oscillation is an intermediate process that explains damped oscillation dynamic systems induced by some attenuated oscillations. The biological population model of time-fractional spatial diffusion depicts a rapid rise in population density in an ecosystem in a suitable habitat that is migrating from an unfavourable zone.
1. Introduction
Fractional calculus is used to analyse, acclimate and implement more efficient and capable systems in numerous real-world fields, like economics, physics, life sciences and engineering including financial economy, electromagnetics, viscoelasticity, epidemiological models, diffusion in cell membranes, fluid mechanics, image processing. Because of their memory and heredity concerns, fractional differential equations are best characterized to model physical and biological processes [Citation1,Citation2]. Various analytical and numerical approaches are taken in many research attempts to achieve exact and approximate solutions for fractional-order differential equations.
Elzaki transform introduced in Ref. [Citation3] is itself and with combinations of other techniques is being used for solving various classical and fractional, linear and non-linear, ordinary and partial differential equations describing various physical and biological process. In Ref. [Citation4], the Elzaki transform is used to solve partial differential equations, whereas in Ref. [Citation5] it is used to solve non-linear partial differential equations of fractional order. In Ref. [Citation6], the Homotopy perturbation method with Elzaki transform is used to solve a system of non-linear partial differential equations. In Ref. [Citation7], the Elzaki transformed method is used to provide an analytical solution of a linear fractionally damped oscillator. In Ref. [Citation8], the radial diffusivity and shock wave equations are solved using the Elzaki variational iteration method, and in Ref. [Citation9], a numerical method based on the Elzaki transform and He’s Homotopy perturbation method (HPM) is developed for solving non-linear partial differential equations arising in spatial flow that characterize the general biological population.
The Elzaki decomposition method (EDM), which is a combination of the Elzaki transform and the Adomian decomposition method first developed in Ref. [Citation10]. EDM is a straightforward and practical approach for solving both linear and non-linear differential equations. It is highly recommended to solve linear and non-linear differential equation with shorten calculation part as compared to other methods.
In this work, firstly the EDM is applied to fractional relaxation and damped oscillation equations with a Caputo fractional derivative of order with and , respectively, to find analytic solution involving Mittag–Leffler function. The plots for different values of are analysed for the solutions obtained. After that, application of EDM is adopted to solve and analyse time-fractional spatial diffusion biological population model for both Malthusian law and Verhulst law in suitable habitat situations from various unfavourable zones.
2. Preliminaries
The Elzaki transform of a function is defined by the integral [Citation3] (2.1) (2.1) where M being finite constant and being finite or infinite.
The inverse of Elzaki transform is denoted by .
ELzaki transform of the power series function is [Citation3] given by (2.2) (2.2) Also, for any (2.3) (2.3) The Elzaki transform of Caputo fractional derivative of a function is given by [Citation5] (2.4) (2.4)
3. Fractional relaxation oscillation equation
In physical, chemical and biological processes, an oscillator is something that exhibits a rhythmic periodic reaction. Many real mechanical, radio-technical, biological and other things have oscillatory processes in which a slow smooth transition of an object’s status over a finite length of time shifts to an irregular change of status over an incredibly short period. The behaviour of a physical system returning to equilibrium after being disrupted provides the basis for a relaxation oscillation. A damped oscillating system is an oscillator that fades away over time owing to energy loss, such as a swinging pendulum, a weight on a spring or a resistor-inductor-capacitor circuit. Relaxation and damped oscillations are described by ordinary differential equations of order one and two, respectively.
In fractional oscillators, the amplitude of the oscillation decays. Because the fractional oscillator is accompanied by a damping mechanism [Citation11], total energy in physical systems dissipates. A comprehensive survey of recent work in the subject of fractional calculus application to dynamic solid mechanics problems is provided in Ref. [Citation12]. We present a fractional differentiation-based model of relaxation and damped oscillation that is more relevant for the physical nature of diverse oscillatory processes with damping.
In order to investigate the mechanical behaviour of objects fractional-order derivatives are preferably applicable in place of classical integer-order. In the relaxation and oscillation models, fractional derivatives are used to for accurate modelling of super-slow relaxation and damped oscillation dynamic systems that require accurate modelling as [Citation13]: (3.1) (3.1) where is a constant that is not negative. For , this Equation (3.1) is known as the fractional relaxation oscillation equation. Which is generated by replacing the first and second-order time derivatives in the related classical relaxation and damped oscillation equations with a Caputo fractional derivative of order with and correspondingly.
3.1. Elzaki decomposition method on fractional relaxation oscillation equation
For a field variable function and a given continuous function with t ≥ 0, we consider the fractional relaxation oscillation differential equation of order as [Citation13]: (3.2) (3.2) (3.3) (3.3) where is a positive constant and denotes Caputo fractional derivative of order .
On taking the Elzaki transform of (3.2) by using (2.4), subject to the initial condition (3.3), we have (3.4) (3.4) On applying Elzaki inverse transform on (3.4) and decomposing as , we have (3.5) (3.5) Applying the Adomian decomposition method, the components of using (2.2), recursively be obtained by (3.6) (3.6) (3.7) (3.7) The components can be obtained when the Elzaki Transform and the Elzaki inverse transform are applied to (3.6) and (3.7) while taking some particular continuous function , resulting the general solution as In particular, if we take in (3.2) and (3.3), we arrive at
Problem 3.1: Consider the following fractional relaxation oscillation differential equation (3.8) (3.8) Solution: In view of (3.6) and (3.7), here we have (3.9) (3.9) which gives the solution for (3.8), on using (2.3) as (3.10) (3.10) where , denotes the Mittag–Leffler function [Citation14].
For the sake of convenience, we will set and to show graphs of the solution (3.10) of fractional relaxation oscillation differential Equation (3.8) in Figure , without losing generality.
The solution of the fractional relaxation Equation (3.10) decays significantly more slowly than the solution of the classical relaxation equation ().Because of its protracted fall, fractional relaxation is often characterized as a super-slow phenomenon.
Next, if we take in (3.2) and (3.3), we reach at
Problem 3.2: Consider the following damped oscillation fractional differential equation as follow: (3.11) (3.11) In view of (3.6) and (3.7), here we have (3.12) (3.12) The components are obtained as
Resulting in the general solution of (3.11) as (3.13) (3.13) where and denote the Mittag–Leffler functions [Citation14,Citation15]. This solution (3.13) is same as the solution obtained by Achar et al. [Citation16] by using Laplace transform method in solving integral equation of motion of harmonic oscillator.
Again, for convenience, we will set and to show graphs of the solution (3.13) of damped oscillation fractional differential Equation (3.11) in Figure , without losing generality.
The solution of damped oscillation fractional differential equation demonstrates asymptotic algebraic decay rather than permanent oscillations, in contrast to the solution of classical oscillation equation . Only some attenuated oscillations appear as a temporary condition in the fractional solution, with their number and initial amplitude increasing with . Because of the presence of intermediate components, fractional oscillations are usually referred to as an intermediate process.
4. Fractional biological population equation
Population growth, which is one of the world’s primary issues, is influenced by a variety of biological and non-biological variables. Globalization, which is a sort of geographical diffusion, is causing population to expand at an alarming rate. The spread of diseases has been accelerated by worldwide migration, but global advancements in public health, medical and family planning technology have played a critical role in global increases in life expectancy [Citation17–19]. To better understand how dispersal impacts species evolution, population dynamics and distribution, mathematical biologists typically use models in which movement is represented by a spatial diffusion process [Citation20,Citation21].
A proliferation of a phenomenon that has a source and diffuses outwards through space and time from constrained sources, such as the spread of a fire or pollution, is known as spatial diffusion. Reference [Citation22] demonstrated the conceptual applicability of the spatial diffusion process to represent redistribution of species in space due to random migration of individuals for any population density.
The spatial diffusion of the population is based on the premise of dispersal based on the irregular motion of each individual in a biological population [Citation23]. Mathematical model of spatial diffusion can be used in a variety of settings, including innovation acceptance, cultural development, epidemiology, settlement, migration and landscape evolution [Citation24,Citation25]. Diffusion systems-based mathematical modelling of biological populations, with applications in ecology, epidemiology and neural systems are focused in Ref. [Citation26].
In recent research, the concept of fractional generalization has been demonstrated to be a successful technique for modelling a variety of natural events. Various fractional diffusion models have been developed to predict the human brain’s response to external stimuli [Citation27], the role of structural heterogeneity in repolarization dispersion in cardiac electrical propagation [Citation28], charge transfer in dye-sensitized solar cells [Citation29], dust aerosols floating in Mars’ atmosphere that cause attenuation of solar radiation traversing the atmosphere [Citation30], and fluid transport through porous media [Citation31], impact on raise of environmental pollution and occurrence in biological populations by presenting a mathematical model relating to incomplete H-function [Citation32] and many others. In these models, fractional derivatives have been incorporated into partial differential equations, allowing for the description of anomalous diffusion within fractal spaces.
The non-linear partial differential Equation (4.1) describes the spatial diffusion model of biological population in a region based on the classical model given by Skellam [Citation20] (4.1) (4.1) where is the population density, that provides the overall population at a given period by integrating the number of inhabitants per unit volume throughout the region . The population supply is a non-linear function of and , expresses the rate at which inhabitants are added to the population (per unit volume) through births and deaths.
Now, we consider constitutive equations for as the Malthusian law and Verhulst law given by and , respectively. If , the entire space is eventually inhabited, regardless of the beginning data; if , and the population is originally contained within a finite interval, then the population is permanently contained inside a finite interval.
Also, focus on the case has a specific form; this allows us to apply the theory to the case in which several well-known results, such as showing the diffusion velocity of inhabitants at the wave front determines the finite speed of propagation for disturbances moving through unpopulated areas, can be demonstrated. (4.2) (4.2) The partial differential Equation (4.1) will then be reduced to (4.3) (4.3) Now, we suppose that without losing generality; this normalization is simply performed by changing spatial coordinates.
We investigate a Caputo fractional generalization of spatial diffusion model of biological population [Citation21] to simulate genuine movements arising from random scattering of species, which will be relevant to many natural biological circumstances. The initial-value problem under study entails locating a density distribution that satisfies certain conditions, with the underlying region assumed to be the real line R is given by (4.4) (4.4) the initial condition: (4.5) (4.5)
3.2. Elzaki decomposition method on time-fractional biological population equation
Now, by taking fractional generalization of biological population problem given by (4.4) and (4.5) with respect to time, we arrive at the following fractional-order biological population problem (4.6) (4.6) the initial condition: (4.7) (4.7) where denotes Caputo fractional derivative of order of population density and population supply .
On taking the Elzaki transform of (4.6), by using (2.4) subject to the initial condition (4.7), we have (4.8) (4.8) On applying Elzaki inverse transform on (4.8) and decomposing as , we have (4.9) (4.9) where represent the Adomian polynomials the non-linear term , which
for Malthusian law is given by (4.10) (4.10) (4.11) (4.11) (4.12) (4.12) and so on,
for Verhulst law are given by (4.13) (4.13) (4.14) (4.14) (4.15) (4.15) and so on.
Identifying the zero component by and the following recurrence relation will be used to determine the remaining components (4.16) (4.16) (4.17) (4.17) The components can be obtained when the Elzaki Transform and its inverse are applied to (4.17), resulting the general solution as .
Now by taking for Malthusian law and the effect of random dispersal [Citation20], in (4.6) and (4.7), and by considering following two cases.
Case (i): An unfavourable zone that extends endlessly in one direction borders a suitable ecosystem, i.e. .
Case (ii): A suitable habitat surrounds an unfavourable segment on both ends, i.e. .
Problem 4.1: We reach at the following fractional model of biological population (4.18) (4.18) with initial conditions for the two cases
when , .
when , .
Solution: Case (i): In view of (4.16) and (4.17), here we have (4.19) (4.19) which gives (4.20) (4.20) (4.21) (4.21) and so on.
Resulting in the solution of (4.18) in Case (i) as (4.22) (4.22) For the sake of convenience, we will set and to show graphs of the solution (4.22) in Figure , without losing generality
For , these graphs show a rapid increase in population density in an ecosystem that is migrating from an unfavourable zone that never ends, as alpha decreases by solution to the classical biological population differential equation in Case (i).
Case (ii): Similarly, we can find the solution of (4.18) and analyse graphically as (4.23) (4.23)
For the sake of convenience, we will set and to show graphs of the solution (4.23) in Figure , without losing generality
For , as declines by the solution to the conventional biological population differential equation in Case (ii), these graphs demonstrate a rapid increase in population density in a suitable habitat around an unfavourable segment on both ends.
Next, we illustrate the special case for Verhulst law, along with , which is based on probabilistic considerations of walk, in which individuals either stay put or migrate in a population-decreasing direction [Citation33], an unfavourable zone encroaches upon a suitable habitat and extends endlessly in one direction. That is, , .
Problem 4.2: We reach at the following fractional biological population equation (4.24) (4.24) with initial conditions , .
Solution: In view of (3.6) and (3.7), here we have (4.25) (4.25) which in view of (4.13), (4.14) and (4.15) gives us and so on.
Resulting in the solution of (4.24) with given initial condition as (4.26) (4.26) For the sake of convenience, we will set and to show graphs of the solution (4.23) in Figure , without losing generality
The population density has increased, as seen in these graphs. Individuals that migrate to or remain in an undesirable area that never ends in one direction encroach on a suitable habitat as lowers by the solution to the classical biological population differential equation for Verhulst law.
5. Conclusion
We have obtained the analytic solutions for fractional relaxation oscillation, fractional damped oscillation and time-fractional biological population differential equations, using the Elzaki decomposition method (EDM). The effectiveness of the method is demonstrated by the predicted accurate analytical solutions and visual observations obtained by varying the order of the fractional derivative to produce various outcomes. The fractional relaxation equation solution decays substantially slower than the classical relaxation equation solution. The fractional solution for damped oscillation fractional differential equation only exhibits a few reduced oscillations as a temporary condition, with their number and initial amplitude increasing with . In the time-fractional biological population equation, the population supply has been considered according to Malthusian law and Verhulst law. Two cases have been studied for Malthusian law, in which an unfavourable zone that extends endlessly in one direction and unfavourable territory is encircled on both ends by a suitable habitat. In all the cases, for , the graphs indicate a sharp rise in population density as falls and people shift from unfavourable areas towards a good environment.
Acknowledgements
The authors acknowledge the anonymous reviewers for their purported insights, their comments and suggestions for enhancing an earlier draft of the manuscript, as well as for helping to shape the paper as it is now.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Data availability statement
No data were used for this study.
References
- Etemad S, Avci I, Kumar P, et al. Some novel mathematical analysis on the fractal–fractional model of the AH1N1/09 virus and its generalized Caputo-type version. Chaos, Solitons Fractals. 2022;162:112511.
- Butt AI, Ahmad W, Rafiq M, et al. Numerical analysis of Atangana-Baleanu fractional model to understand the propagation of a novel corona virus pandemic. Alexandria Eng J. 2022;61(9):7007–7027.
- Elzaki TM. The new integral transform Elzaki transform. Global J Pure Appl Math. 2011;7(1):57–64.
- Elzaki TM. Application of new transform ‘Elzaki transform’ to partial differential equations. Glob J Pure Appl Math. 2011;7(1):65–70.
- Neamaty A, Agheli B, Darzi R. New integral transform for solving nonlinear partial Differential equations of fractional order. Theory Approx Appl. 2013;10(1):69–86.
- Elzaki TM, Biazar J. Homotopy perturbation method and Elzaki transform for solving system of nonlinear partial differential equations. World Appl Sci J. 2013;24(7):944–948.
- Suleman M, Lu D, Rahman JU, et al. Analytical solution of linear fractionally damped oscillator by Elzaki transformed method. DJ J Eng Appl Math. 2018;4(2):49–57.
- Elzaki TM, Kim H. The solution of radial diffusivity and shock wave equations by Elzaki variational iteration method. Int J Math Anal. 2015;9:1065–1071.
- Rahman J, Lu D, Suleman M, et al. He–Elzaki method for spatial diffusion of biological population. Fractals. 2019 Aug 13;27(05):1950069.
- Khalid M, Sultana M, Zaidi F, et al. An Elzaki transform decomposition algorithm applied to a class of non-linear differential equations. Natural Sci. Res. 2015;5:48–56.
- Achar BN, Hanneken JW, Clarke T. Damping characteristics of a fractional oscillator. Physica A. 2004;339(3-4):311–319.
- Rossikhin YA, Shitikova MV. Application of fractional calculus for dynamic problems of solid mechanics: novel trends and recent results. Appl Mech Rev. 2010;63(1):010801.
- Gorenflo R, Mainardi F. Fractional relaxation-oscillation phenomena. In Handbook of fractional calculus and applications. 2019. p. 45–74.
- Mittag-Leffler GM. Sur la nouvelle fonction. CR Acad Sci Paris. 1903;137(2):554–558.
- Wiman A. Über den Fundamentalsatz in der Teorie der Funktionen Ea (x). Acta Mathematica. 1905;29:191–201.
- Achar BN, Hanneken JW, Enck T, et al. Dynamics of the fractional oscillator. Physica A. 2001;297(3-4):361–367.
- Dayan F, Rafiq M, Ahmed N, et al. Design and numerical analysis of fuzzy nonstandard computational methods for the solution of rumor based fuzzy epidemic model. Physica A. 2022;127542.
- Raza A, Baleanu D, Yousaf M, et al. Modeling of anthrax disease via efficient computing techniques. Intell. Autom. Soft Comput. 2022;32(2):1109–1124.
- Naveed M, Baleanu D, Raza A, et al. Treatment of polio delayed epidemic model via computer simulations. CMC-Comput Mater Continua. 2022;70(2):3415–3431.
- Skellam JG. Random dispersal in theoretical populations. Biometrika. 1951;38(1/2):196–218.
- Gurtin ME, MacCamy RC. On the diffusion of biological populations. Math Biosci. 1977;33(1-2):35–49.
- Okubo A. Diffusion and ecological problems: mathematica1 models; 1980.
- Pielou EC. An introduction to mathematical ecology. New York: Wiley-Inter-science; 1969.
- Morrill R, Gaile GL, Thrall GI. Spatial diffusion; 2020.
- Kheirkhah F, Hajipour M, Baleanu D. The performance of a numerical scheme on the variable-order time-fractional advection-reaction-subdiffusion equations. Appl Numer Math. 2022;178:25–40.
- Upadhyay RK, Iyengar SR. Spatial dynamics and pattern formation in biological populations. New York: Chapman and Hall/CRC; 2021.
- Namazi H, Kulish VV. Fractional diffusion based modelling and prediction of human brain response to external stimuli. Comput Math Methods Med. 2015;2015:1–11.
- Bueno-Orovio A, Kay D, Grau V, et al. Fractional diffusion models of cardiac electrical propagation: role of structural heterogeneity in dispersion of repolarization. J R Soc Interface. 2014;11(97):20140352.
- Maldon B, Thamwattana N. A fractional diffusion model for dye-sensitized solar cells. Molecules. 2020;25(13):2966.
- Ali I, Malik N, Chanane B. Fractional diffusion model for transport through porous media; 2014.
- Jiménez S, Usero D, Vázquez L, et al. Fractional diffusion models for the atmosphere of Mars. Fractal Fract. 2017;2(1):1.
- Purohit SD, Khan AM, Suthar DL, et al. The impact on raise of environmental pollution and occurrence in biological populations pertaining to incomplete H-function. Natl Acad Sci Lett. 2021;44:263–266. DOI:10.1007/s40009-020-00996-y
- Gurney WS, Nisbet RM. The regulation of inhomogeneous populations. J Theor Biol. 1975;52(2):441–457.