
In this paper, a new modified hybrid explicit group (MHEG) iterative method is presented for the efficient and accurate numerical solution of a time-fractional diffusion equation in two space dimensions. The time fractional derivative is defined in the Caputo sense. In the proposed method, a Laplace transformation is used in the temporal domain, and, for the spatial discretization, a new finite difference scheme based on grouping strategy is considered. The unique solvability, unconditional stability and convergence are thoroughly proved by the matrix analysis method. Comparison of numerical results with analytical and other approximate solutions indicates the viability and efficiency of the proposed algorithm.
Citation: Fouad Mohammad Salama, Nur Nadiah Abd Hamid, Norhashidah Hj. Mohd Ali, Umair Ali. An efficient modified hybrid explicit group iterative method for the time-fractional diffusion equation in two space dimensions[J]. AIMS Mathematics, 2022, 7(2): 2370-2392. doi: 10.3934/math.2022134
[1] | Kaile Chen, Yunyun Liang, Nengqiu Zhang . Global existence of strong solutions to compressible Navier-Stokes-Korteweg equations with external potential force. AIMS Mathematics, 2023, 8(11): 27712-27724. doi: 10.3934/math.20231418 |
[2] | Hui Fang, Yihan Fan, Yanping Zhou . Energy equality for the compressible Navier-Stokes-Korteweg equations. AIMS Mathematics, 2022, 7(4): 5808-5820. doi: 10.3934/math.2022321 |
[3] | Qingkun Xiao, Jianzhu Sun, Tong Tang . Uniform regularity of the isentropic Navier-Stokes-Maxwell system. AIMS Mathematics, 2022, 7(4): 6694-6701. doi: 10.3934/math.2022373 |
[4] | Abdelkader Moumen, Ramsha Shafqat, Azmat Ullah Khan Niazi, Nuttapol Pakkaranang, Mdi Begum Jeelani, Kiran Saleem . A study of the time fractional Navier-Stokes equations for vertical flow. AIMS Mathematics, 2023, 8(4): 8702-8730. doi: 10.3934/math.2023437 |
[5] | Jianxia He, Ming Li . Existence of global solution to 3D density-dependent incompressible Navier-Stokes equations. AIMS Mathematics, 2024, 9(3): 7728-7750. doi: 10.3934/math.2024375 |
[6] | Xiaoxia Wang, Jinping Jiang . The pullback attractor for the 2D g-Navier-Stokes equation with nonlinear damping and time delay. AIMS Mathematics, 2023, 8(11): 26650-26664. doi: 10.3934/math.20231363 |
[7] | Jae-Myoung Kim . Time decay rates for the coupled modified Navier-Stokes and Maxwell equations on a half space. AIMS Mathematics, 2021, 6(12): 13423-13431. doi: 10.3934/math.2021777 |
[8] | Linlin Tan, Meiying Cui, Bianru Cheng . An approach to the global well-posedness of a coupled 3-dimensional Navier-Stokes-Darcy model with Beavers-Joseph-Saffman-Jones interface boundary condition. AIMS Mathematics, 2024, 9(3): 6993-7016. doi: 10.3934/math.2024341 |
[9] | Yina Lin, Qian Zhang, Meng Zhou . Global existence of classical solutions for the 2D chemotaxis-fluid system with logistic source. AIMS Mathematics, 2022, 7(4): 7212-7233. doi: 10.3934/math.2022403 |
[10] | Peng Lv . Structured backward errors analysis for generalized saddle point problems arising from the incompressible Navier-Stokes equations. AIMS Mathematics, 2023, 8(12): 30501-30510. doi: 10.3934/math.20231558 |
In this paper, a new modified hybrid explicit group (MHEG) iterative method is presented for the efficient and accurate numerical solution of a time-fractional diffusion equation in two space dimensions. The time fractional derivative is defined in the Caputo sense. In the proposed method, a Laplace transformation is used in the temporal domain, and, for the spatial discretization, a new finite difference scheme based on grouping strategy is considered. The unique solvability, unconditional stability and convergence are thoroughly proved by the matrix analysis method. Comparison of numerical results with analytical and other approximate solutions indicates the viability and efficiency of the proposed algorithm.
Ecology is the scientific study of how living organisms interact with one another and with their surroundings. Interactions between species and between organisms and their physical and chemical environments involve all individuals of the same type. Aquatic ecology is the study of these interactions in all aquatic environments, such as rivers, lakes, streams, ponds, wetlands, and oceans. The physical properties of aquatic organisms have an impact on all types of species found there. The organisms in a specific ecosystem are directly impacted by environmental elements like nutrient content, water flow, temperature, and shelter. Only those species that can adapt to and thrive in a specific habitat will flourish because of the resources at hand. Predation and competition for resources (such as food and habitat) have an impact on species abundance and variety in aquatic ecosystems, as allow interactions between living organisms. In turn, living creatures in a given environment can have an impact on certain characteristics of that environment [1,2,3,4].
We are all aware of this, algal blooms are triggered by a combination of external and internal forces, algal population aggregation and migration are examples of the latter [5]. In the last 10 years, numerous scholars have examined how the growth of algae affects the reproduction of aquatic organisms in environmental studies and aquatic biology [6,7,8,9,10,11,12]. The tight gel layer of microcystis aggregation in particular, filter feeding fish have a hard time digesting it. Even if the filter-feeding fish eat them, a substantial percentage of undamaged active algal cells will remain in the fish feces. Xie and Liu [13] have demonstrated that aggregations of microcystis can generate more algal toxins than solitary cells, reducing the danger of predation. On other hand, Sommer et.al. [14] have shown the food web in the aquatic environment that there are fish that feed on algae and predatory fish that feed on small fish (prey).
However, on the mentioned area, there has been little scientific investigation into the population dynamics and mathematical modeling of the dynamic impacts of algal aggregation. Furthermore, population dynamics and the ecological mathematical model a complete research technique combines coordination, organization, and purpose to create three foundations for studying biological systems. Analysis of dietary intake, the view from the model, and dynamic presentation are the three kinds (see [15,16]). The basic goal of population dynamics modeling is to gain a better knowledge of how populations interact and how they are affected by internal and external factors [17,18]. However, a thorough investigation of ecological mathematical models and population dynamics has been made in recent decade. In addition, a slew of additional articles has yielded key research findings in population dynamics and ecological mathematical models (we refer [19,20,21,22,23]).
Here we remark that authors [24] represented the following system
{dNdt=rN(1−NK1)−α1(N−m1)Pc1+N−m1,dPdt=d1P(1−NK1)+(β1α1(N−m1)Pc1+N−m1)−γ1P, | (1) |
which has created a fresh aquatic ecological model to explain the microcytic aeruginosa aggregation impact and filter-feeding fish. Further, the model describes the connection between microcytic aeruginosa and filter feeding fish.
By memorizing the foregoing, and depending on the investigation carried out by Sommer et.al. [14], we divided the aquatic environment into the following:
Our aquatic community=microcystis aeruginosa(N)+filter feeding fish(P)+predatory fish(F)Our food chain:N→P→F. |
According to recent studies, predators have an impact on prey populations not only by killing them directly, but also through instilling fear in prey animals. This significantly affects their pace of reproduction (see [25,26]). Predator fear may affect the behavioral patterns and physiology of some prey species. This is more terrible than killing someone directly [27,28]. Wang et al. [26], which has suggested that raising the cost of fear can help to stabilize the system by incorporating the fear effect into a predator-prey model. Pandey et al. [29] investigated a tri-topic food chain model with Holing type-Ⅱ functional response by incorporating fear impact. Mentioned authors have discovered that a low amount of fear causes chaotic behavior, whereas a high level of fear allows the chaotic oscillations to be controlled. Predator fear may affect the behavioral patterns and physiology of some prey species; this is more terrible than killing someone directly.
In this study, we have expanded the first system from the second dimension to the third dimension, where we added a very important component in the aquatic ecological model. The said target is the large predatory fish that has an effect on each of the first two components (the microcytic aeruginosa aggregation and filter feeding fish). In addition, we studied the effect of fear of predatory fish (big fish) on the filter feeding fish because the main source of growth for filters feeding from microcytic aeruginosa fish, when they were fed on microcytic aeruginosa. On other hand, the harvesting effect on the predatory fish is considered.
Keeping the above discussion in mind, we add the function of fear for food filter feeding fish from predatory fish. Therefore, the mathematical formulation of aquatic environment models takes the form
{dNdt=rN(1−Nk)−α1(N−m1)Pc1+N−m1dPdt=(11+δF)(β1α1(N−m1)Pc1+N−m1)dFdt=β2ϵ1PF−cEF,−ϵ1PF | (2) |
where
● r is the rate of intrinsic growth of microcystis aeruginosa.
● k is the maximum environmental capacity.
● α1 is the grazing coefficient.
● m1 is the microcystis aeruginosa aggregation parameter.
● c1 is the half-saturation constant.
● β1 is the absorption coefficient.
● c is the harvesting's capacity coefficient.
● E is the subject of big fish population to a combined harvesting effort.
● δ is the level of fear which divides anti- predatory fish behaviour of the filter feeding fish.
● β2 is efficiency of conversion of consumed filter-feeding fish into big fish.
In the process of dynamic systems analysis, we try as much as possible to facilitate the calculations. Therefore, since model (2) has 11 parameters, hence, we are going to do the Non-dimensionalization for all variables and parameters as:
N=c1x,P=rc1α1y,F=1δz,T=1rt,k1=kc1,k2=m1c1,k3=β1α1r,k4=ϵ1rδ,k5=β2ϵ1c1α1,k6=cEr |
Hence, the dimensional system (2) corresponding the following non-dimensional system:
{dxdT=x(1−xk1)−(x−k2)yx+1−k2dydT=k31+z((x−k2)yx+1−k2)−k4yzdzdT=k5yz−k6z. | (3) |
All parameters in model (3) are expected to be positive. The following are the initial conditions for system (3)
x(0)≥0,y(0)≥0,z(0)≥0. |
In this part, we first demonstrate the positivity and boundedness of solutions for the system (3). These are quite significant in terms of the model's validity.
Lemma 1. For every T>0, all solutions (x(T),y(T),z(T)) of the system (3) with initial values (x(0),y(0),z(0))∈R3+, are positive.
Proof. It is simple to show that the solutions of the system (3) are nonnegative under the initial conditions: x(0)≥0,y(0)≥0, and z(0)≥0.
In the first, second and third instances, we have gotten respectively as
x′(T)≥xyx+1−k2≥0y′(T)=y(T)[k31+z(T)((x(T)−k2)x(T)+1−k2)−k4z]z′(T)=z(T)[k5y(T)−k6]. |
So
x(T)≥x(0)eT∫0y(s)x(s)+1−k2ds≥0,y(T)=y(0)eT∫0[k31+z(s)((x(s)−k2)y(s)x(s)+1−k2)−k4z(s)]ds≥0z(T)=z(0)eT∫0(k5y(s)−k6)ds≥0. |
Therefore, all solutions of the model (3) are positive.
The system (3) exhibits continuous interaction functions and continuous partial derivatives, indicating that the solution exists and is unique. Furthermore, to guarantee convergence of the solution to an attractor, the solution of the system (3) is proven uniformly bounded, as shown in the following theorem.
Lemma 2. The feasible region of system (3), where all solutions can be found is given by
Σ={(x,y,z)∈R3+;x≤k1;L1≤ξ2W}. | (4) |
Proof. We have the following result from the first equation: dxdT≤x(1−xk1), therefore, when T→∞, then x≤k1. Now let us consider L1=x+1k3y+k4k3k5z, the derivative of L1 with regard to time can thus be expressed as
dL1dT=x(1−xk1)−(x−k2)yx+1−k2+11+z((x−k2)yx+1−k2)−k4k3yz+k4k3yz−k4k6k3k5z≤2x+1k3y−x−1k3y−k4k6k3k5z≤(2k1+ξk3)−WL1, |
where W =minT→∞{1,k6},ξ1=supT→∞y(T). Then, using direct computation, one can see that for T→∞, we obtain
L1≤ξ2W, | (5) |
where ξ2=2k1+ξk3.
Before starting with the existence of equilibrium points of (3) and its discussion, we present the Jacobian matrix of (3) as:
J(x,y,z)=(1−2k1x−y(x+1−k2)2−x−k2x+1−k20k31+z(y(x+1−k2)2)k31+z(x−k2x+1−k2)−k4z−k3(1+z)2((x−k2)yx+1−k2)−k4y0k5zk5y−k6). | (6) |
The equilibrium of the system (3) is as follows:
● The trivial equilibrium point u0(0,0,0), which always exists. To study its local stability, the characteristic equation is computed from the given Jacobian matrix
J0(0,0,0)=(1k21−k200−k2k31−k2000−k6), | (7) |
as
(1−λ)(−k2k31−k2−λ)(−k6−λ)=0. | (8) |
Clearly, (8) has three distinct roots λ1=1,λ2=−k2k31−k2 and λ3=−k6. The Routh–Hurwitz criterion states that the trivial equilibrium point is locally asymptotically stable, when its all eigenvalue has negative real parts. Notice that only λ3<0, and λ1 is purely positive, so the trivial equilibrium point u0(0,0,0) is saddle node point.
● The x-axis equilibrium point u1(k1,0,0), which is always exists. The characteristic equation of the following Jacobin matric around u1
J1(k1,0,0)=(−1−k1−k2k1−k2+100k3(k1−k2)k1−k2+1000−k6), | (9) |
is given by
(−1−λ)(k3(k1−k2)k1−k2+1−λ)(−k6−λ)=0. | (10) |
Clearly (10) has three distinct roots λ1=−1,λ2=k3(k1−k2)k1−k2+1 and λ3=−k6. The Routh–Hurwitz criterion establishes that the x-axis equilibrium point is locally asymptotically stable when its all eigenvalue has negative real parts. Notice that λ1,λ3<0, and λ2 is negative if k1<k2<k1+1. Hence, under this condition, the x-axis equilibrium point u1(k1,0,0) is locally asymptotically stable.
Theorem 1. The trivial equilibrium point u0(0,0,0) and the x-axis equilibrium point u1(k1,0,0) exist always and exhibit the following local behaviors:
1. the trivial equilibrium point u0(0,0,0) is saddle point.
2. the x-axis equilibrium point u1(k1,0,0) is locally asymptotically stable provided that:
−1<k1−k2<0. |
In the next part, we have to discussed the existence of internal equilibrium point. If the internal equilibrium does not exist, the model (3) will not be durable, implying that the filter-feeding fish and predatory fish would become extinct. But this is not good for ecological balance. As a result, it is critical that we investigate the model's complicated dynamic having properties of the Model (3). We know the existence condition of equilibrium from the foregoing analysis. Therefore, the next step is to investigate stability. (Note that we only discuss equilibrium stability when it exists.)
● The interior equilibrium point is denoted by u∗(x∗,y∗,z∗), where
y∗=k6k5z∗=−1+√1−4(k3(x∗−k2)k4(k2−1−x∗))2 |
Hence, z∗ is positive provided that:
k2<min{x∗,1,} | (11) |
and x∗ is a root of the following third degree polynomial
x3+R1x2+R2x+R3=0, | (12) |
where,
R1=1−k1−k2R2=k1(k2+y∗−1)R3=−k1k2y∗. |
Therefore, according the Descartes role, the following Lemma summarizes the existence of the interior equilibrium point u∗(x∗,y∗,z∗).
Lemma 3: Under the condition (11), the model (2) has:
1) a unique interior equilibrium point u∗(x∗,y∗,z∗) where one of the following is satisfied
ⅰ- 1>k1+k2,k6>k5(1−k2);
ⅱ- 1<k1+k2,k6<k5(1−k2);
ⅲ- 1>k1+k2,k6<k5(1−k2).
2) A triple interior equilibrium points u∗1(x∗1,y∗,z∗),u∗2(x∗2,y∗,z∗) and u∗3(x∗3,y∗,z∗) if the following is satisfied
ⅰ- 1<k1+k2,k6>k5(1−k2).
To simplify the calculations, let's defined the following two quantities
k∗03=11+z∗,k∗13=1x∗+1−k2. |
The Jacobin matrix (6) around u∗(x∗,y∗,z∗) is given by
J3(x∗,y∗,z∗)=(1−2k1x∗−k∗213y∗−k∗13(x∗−k2)0k∗03k∗213k3y∗k∗03k∗213k3(x∗−k2)−k4z∗−[k3k∗032k∗13(x∗−k2)−k4]y∗0k5z∗0). | (13) |
The characteristic equations of (13) is,
λ3+A∗λ2+B∗λ+C∗=0, | (14) |
where,
A∗=−1+2k1x∗+k∗213y∗+k4(z∗(1+k∗13)),B∗=k4(2k1x∗−1)(2z∗)+k5[k3k∗032k∗13(x∗−k2)+k4]y∗z∗+k∗132[k3k∗032k∗13(x∗−k2)+k4]y∗z∗+k∗03k3k∗133(1−k2)(1−k∗13)C∗=(2k1x∗+k∗213y∗−1)(k∗03k3k∗13(x∗−k2)y∗+k4y∗). |
Given the coefficient A∗, which according to Routh-Hwartze Criteria must be both it and coefficient C∗ with a positive sign, then A∗ is positive using the condition
x∗>k12. | (15) |
It is quite easy to verify that the same condition (15), leads to C∗ being positive as well. According to Ruth-Hwartz, the characteristic equation (14) represents eigenvalues with negative real parts in order for us to achieve the stability property of the positive point, if A∗,C∗ and A∗B∗−C∗>0. Therefore, by forward calculations its easy to verify that A∗B∗−C∗>0 provided the conditions (11) and (15). Hence, the positive equilibrium point u∗(x∗,y∗,z∗) is locally asymptotically stable.
Theorem 3. Suppose that the conditions of existence for the positive equilibrium point u∗(x∗,y∗,z∗), the point u∗ is locally asymptotically stable if the condition (15) is met.
Remark 1. The situation of extinction of all species is continually in flux. Other equilibrium states, on the other hand, are conditionally stable. Even if there is no microcystis aeruginosa in the environment, this might be taken as implying that both filter-feeding and predatory fish cannot go extinct at the same time. Because of maximum environmental capacity and microcystis aeruginosa aggregation, this occurrence may be explained by the model (3).
In this part, the system (3) is rebuilt in a vector form to explore the local bifurcation that may occur around each equilibrium point.
dPdt=g(P),P=(x,y,z)T,g=(g1,g2,g3)T, |
where gi,i=1,2,3 are functions which given in the right hand side of dimensionless the system (3). The general second derivative of the vector g can now be represented as follows using simple computations:
D2g(P,β)(H,H)=([−2k1+2y(x+1−k2)3]h21−[2(x+1−k2)2]h1h2[−k31+z2y(x+1−k2)3]h21+2[k31+z2(x+1−k2)2]h1h2−2[k3(1+z)2y(x+1−k2)2]h1h3−2[k3(1+z)2x−k2x+1−k2+k4]h2h3−[k3(1+z)2y{k2−x}(x+1−k2)]h232k5h3h2), | (16) |
here U=(h1,h2,h3)T be any eigenvector and k1 be a bifurcation parameter.
Theorem 4: Assume that the condition (15) is not satisfied. Then, system (3) near the positive equilibrium point u∗(x∗,y∗,z∗) has a saddle-node bifurcation.
Proof. According to (3) and condition (15), a simple calculations reveals that the Jacobian matrix of system (1), stated by (7), at E0 with k1≡k1∗=2x∗1−k∗213y∗ has zero eigenvalue λ02(k1∗)=0 and can be written as follows:
J3(x∗,y∗,z∗)=(1−2k1x∗−k∗213y∗−k∗13(x∗−k2)0k∗03k∗213k3y∗k∗03k∗213k3(x∗−k2)−k4z∗−[k3k∗032k∗13(x∗−k2)−k4]y∗0k5z∗0). | (17) |
Assume that V=(v1,v2,v3)T is the eigenvector to λ∗=0, then J3∗V=0 yields:
V=[[k3k∗032k∗13(x∗−k2)−k4]k∗03k∗213k3v30v3], | (18) |
here v3 be any real, non-zero number.
Similarly, if W=(w1,w2,w3)T is the eigenvector to, λ∗=0 of the matrix (J3∗)T, then (J3∗)TW=0 results in
W=[−k5z∗k∗13(x∗−k2)w30w3], | (19) |
where w3 be any real number, not zero. Moreover, it is easy to verify that
dgdk1=gk1∗=(2xk21,0,0)Tandgk1∗(u∗,k1∗)=((x∗k1∗)2,0,0)T. |
Then
WTgk1(u∗,k1∗)=−(1−k∗213y∗)2k5z∗k∗13(x∗−k2)≠0. | (20) |
System (3) near u* features saddle – node bifurcation, according to Sotomayer's theorem.
Now by substitute u∗ and k1∗ in (16) with V instead of H it's observed that
D2g(u∗,k1∗)(V,V)=([−2k1∗+2y∗(x∗+1−k2)3]v21[−k31+z∗2y∗(x∗+1−k2)3]v21−2[k3(1+z∗)22y∗(x∗+1−k2)2]v1v3+2[k3(1+z∗)2y∗{k2−x∗}(x∗+1−k2)3]v230) | (21) |
Accordingly, we obtain that:
WT(D2f(u∗,k1∗)(V,V))=[−2k1∗+2y∗(x∗+1−k2)3]v21w1≠0. |
Hence due to Sotomayer's theorem a saddle–node bifurcation is occur.
Some numerical simulations are offered in this part to back up our analytical and mathematical findings. These simulations also reveal the fascinating complicated of the system's behavior. In the following section, we will show several examples to back up our earlier findings. For a hypothetical collection of data, numerical simulations are run using the matlab 7.0.1 software suite.
Suppose that the following hypothetical collection of data;
k1=0.5,k2=0.35,k3=0.6,k4=0.55,k5=1.1,k6=0.7. | (22) |
Looking at the scalar hypothetical set (22), we notice that, it's given the coefficients of the polynomial R1=0.15>0 and R2=−0.0068<0. By solving model (3) we have, u∗=(0.4402,0.6364,0.0833). Since, k2<min{0.4402,1}, so the positive condition k2<min{x∗,1} for z∗ is satisfied. Therefore, according the condition (11) and point (1-ⅰ) of lemma 3 all the conditions of existence for the positive point are fulfilled. On other hand, the local stability condition (15) is satisfied. Hence, Figure 1(a) and Figure 1(b) show us the positive equilibrium point u∗(x∗,y∗,z∗), which appears to be locally asymptotically stable.
Suppose that the following hypothetical collection of data;
k1=7,k2=1,k3=0.17,k4=0.72,k5=0.35,k6=0.52. | (23) |
Looking at the scalar hypothetical values given in (23), we notice that at the coefficients of the polynomial R1=−7<0 and R2=10.4>0. By solving model (3), we have, u∗=(5.4396,1.4857,0.1654). Since, k2≤min{5.4396,1}, so the positive condition k2<min{x∗,1}, for z∗ is satisfied. Therefore, according to the condition (11) and point (1-ⅱ) of Lemma 3, all the conditions of existence for the positive point are fulfilled. On other hand, the local stability condition (15), at k22=3.5<5.4396 is satisfied. Hence, Figure 2 (a) and (b) shows us the positive equilibrium point u∗(x∗,y∗,z∗), which appears to be locally asymptotically stable.
Suppose that the following hypothetical collection of data
k1=15,k2=15.5,k3=0.501,k4=0.55,k5=0.31,k6=0.17. | (24) |
The reader can notice that the imposed set (24) that contains k1<k2<k1+1 makes the second root of the polynomial (10) with a negative sign. Therefore, the point u1, where the solution went is a stable point (see Figure 3).
Suppose that the following hypothetical collection of data;
k1=15,k2=17,k3=0.501,k4=0.55,k5=0.31,k6=0.17. | (25) |
When we chose the hypothesis data set (25), one can notice that k1<k2 but k2−k1>1 which makes the sign of the second eigenvalue of the characteristic equation (10) negative. So, according to Roth-Hwartize criteria, the point u1 will be completely unstable equilibrium point (see Figure 4 (a). On other hand, This data set gives values for the parameters R1=−31.5 and R2=17.0484 which are represent the coefficients of a polynomial (12) whose number of roots determines the existence of the positive equilibrium point u∗. As R1=−31.5<0 and R2=17.0484>0, so, according to Lemma 3, the model (3) admit a unique positive equilibrium point. The local stability of positive point is investigate since the condition (15) which gives a sufficient condition for routh howartize theorem is satisfied as 14.1606 = x∗>k12=7.5 (see Figure 4 (b)).
Suppose that the following hypothetical collection of data;
k1=0.5,k2=0.45,k3=0.501,k4=0.55,k5=0.31,k6=0.17. | (26) |
This data set gives values for the parameters R1=−31.5 and R2=17.0484 which are represent the coefficients of a polynomial (12) whose number of roots determines the existence of the positive equilibrium point u∗. As R1=0.1000>0 and R2=−0.0516<0, so, according to Lemma 3-ii of the model (3) admit a unique positive equilibrium point. The local stability of positive point is investigate since the condition (15) which gives a sufficient condition for routh howartize theorem is satisfied as 0.4640=x∗>k12=0.25 (see Figure 5).
In this article, we investigated the fear effect on filter feeding fish from predatory fish and harvesting effect of predatory fish of an aquatic ecological model of kind microcytic aeruginosa-filter feeding fish-predatory fish model with aggregation effect of microcytic aeruginosa. In the begging, to reduce the number of parameters of the model (2), we use non dimensionlization. The positivity and feasible region for the solution of the model (3) have been derived. In addition, we have proved some qualitative analysis by using Routh-Hurwitz stability criterion. In the numerical section, we have presented the approximate solution and their graphical presentation in various Figures including Figure 6 (a), Figure 6 (b), Figure 7 (a), Figure 7 (b) and Figure 8 to Figure 13 respectively. Some graphs of time series and phase portrait are given in this section also. As, we explained the outcomes one at a time by figuring out the solution and displaying the figures, this section demonstrated the validity of the theoretical conclusions obtained by this study. In future, the proposed study will be updated by using the concept of fractional calculus.
There does not exist any conflict of interest.
Author Kamal Shah and Thabet Abdeljawad would like to thank Prince Sultan University for paying the APC and support through TAS research lab.
[1] |
H. R. Ghehsareh, A. Zaghian, S. M. Zabetzadeh, The use of local radial point interpolation method for solving two-dimensional linear fractional cable equation, Neural Comput. Appl., 29 (2018), 745–754. doi: 10.1007/s00521-016-2595-y. doi: 10.1007/s00521-016-2595-y
![]() |
[2] |
Y. L. Zhao, T. Z. Huang, X. M. Gu, W. H. Luo, A fast second-order implicit difference method for time-space fractional advection-diffusion equation, Numer. Func. Anal. Opt., 41 (2020), 257–293. doi: 10.1080/01630563.2019.1627369. doi: 10.1080/01630563.2019.1627369
![]() |
[3] |
M. Hussain, S. Haq, Weighted meshless spectral method for the solutions of multi-term time fractional advection-diffusion problems arising in heat and mass transfer, Int. J. Heat Mass Tran., 129 (2019), 1305–1316. doi: 10.1016/j.ijheatmasstransfer.2018.10.039. doi: 10.1016/j.ijheatmasstransfer.2018.10.039
![]() |
[4] |
M. Abbaszadeh, A. Mohebbi, A fourth-order compact solution of the two dimensional modified anomalous fractional sub-diffusion equation with a nonlinear source term, Comput. Math. Appl., 66 (2013), 1345–1359. doi: 10.1016/j.camwa.2013.08.010. doi: 10.1016/j.camwa.2013.08.010
![]() |
[5] |
G. M. Zaslavsky, Chaos, fractional kinetics, and anomalous transport, Phys. Rep., 371 (2002), 461–580. doi: 10.1016/S0370-1573(02)00331-9. doi: 10.1016/S0370-1573(02)00331-9
![]() |
[6] | S. G. Samko, A. A. Kilbas, O. I. Marichev, Fractional integrals and derivatives: Theory and applications, Yverdon: Gordon and Breach, 1993. |
[7] | I. Podlubny, Fractional differential equations, mathematics in science and engineering, New York: Academic Press, 1999. |
[8] | B. Guo, X. Pu, F. Huang, Fractional partial differential equations and their numerical solutions, Singapore: World Scientific, 2015. |
[9] | A. Kochubei, Y. Luchko, V. E. Tarasov, I. Petra, Handbook of fractional calculus with applications, Berlin: De Gruyter Grand Forks, 2019. |
[10] |
J. Y. Shen, Z. Z. Sun, R. Du, Fast finite difference schemes for time fractional diffusion equations with a weak singularity at initial time, E. Asian J. Appl. Math., 8 (2018), 834–858. doi: 10.4208/eajam.010418.020718. doi: 10.4208/eajam.010418.020718
![]() |
[11] |
A. Chen, C. Li, A novel compact adi scheme for the time-fractional subdiffusion equation in two space dimensions, Int. J. Comput. Math., 93 (2016), 889–914. doi: 10.1080/00207160.2015.1009905. doi: 10.1080/00207160.2015.1009905
![]() |
[12] |
G. H. Gao, Z. Z. Sun, Y. N. Zhang, A finite difference scheme for fractional sub-diffusion equations on an unbounded domain using artificial boundary conditions, J. Comput. Phys., 231 (2012), 2865–2879. doi: 10.1016/j.jcp.2011.12.028. doi: 10.1016/j.jcp.2011.12.028
![]() |
[13] |
M. Tamsir, N. Dhiman, D. Nigam, A. Chauhan, Approximation of Caputo time-fractional diffusion equation using redefined cubic exponential B-spline collocation technique, AIMS Mathematics, 6 (2021), 3805–3820. doi: 10.3934/math.2021226. doi: 10.3934/math.2021226
![]() |
[14] | J. Shen, X. M. Gu, Two finite difference methods based on an H2N2 interpolation for two-dimensional time fractional mixed diffusion and diffusion-wave equations, Discrete Cont. Dyn. B, 2021. doi: 10.3934/dcdsb.2021086. |
[15] |
X. M. Gu, Y. L. Zhao, X. L. Zhao, B. Carpentieri, Y. Y. Huang, A note on parallel preconditioning for the all-at-once solution of Riesz fractional diffusion equations, Numer. Math. Theor. Meth. Appl., 14 (2021), 893–919. doi: 10.4208/nmtma.OA-2020-0020. doi: 10.4208/nmtma.OA-2020-0020
![]() |
[16] |
X. M. Gu, H. W. Sun, Y. L. Zhao, X. Zheng, An implicit difference scheme for time-fractional diffusion equations with a time-invariant type variable order, Appl. Math. Lett., 120 (2021), 107270. doi: 10.1016/j.aml.2021.107270. doi: 10.1016/j.aml.2021.107270
![]() |
[17] |
Y. Xu, Y. Zhang, J. Zhao, Backward difference formulae and spectral galerkin methods for the riesz space fractional diffusion equation, Math. Comput. Simulat., 166 (2019), 494–507. doi: 10.1016/j.matcom.2019.07.007. doi: 10.1016/j.matcom.2019.07.007
![]() |
[18] |
X. Gao, B. Yin, H. Li, Y. Liu, Tt-m FE method for a 2D nonlinear time distributed-order and space fractional diffusion equation, Math. Comput. Simulat., 181 (2021), 117–137. doi: 10.1016/j.matcom.2020.09.021. doi: 10.1016/j.matcom.2020.09.021
![]() |
[19] |
X. Zheng, H. Wang, An error estimate of a numerical approximation to a hidden-memory variable-order space-time fractional diffusion equation, SIAM J. Numer. Anal., 58 (2020), 2492–2514. doi: 10.1137/20M132420X. doi: 10.1137/20M132420X
![]() |
[20] |
X. Zheng, H. Wang, A Hidden-memory variable-order time-fractional optimal control model: Analysis and approximation, SIAM J. Control Optim., 59 (2021), 1851–1880. doi: 10.1137/20M1344962. doi: 10.1137/20M1344962
![]() |
[21] |
X. Zheng, H. Wang, Optimal-order error estimates of finite element approximations to variable-order time-fractional diffusion equations without regularity assumptions of the true solutions, IMA J. Numer. Anal., 41 (2021), 1522–1545. doi: 10.1093/imanum/draa013. doi: 10.1093/imanum/draa013
![]() |
[22] |
H. Wang, X. Zheng, Wellposedness and regularity of the variable-order time-fractional diffusion equations, J. Math. Anal. Appl., 475 (2019), 1778–1802. doi: 10.1016/j.jmaa.2019.03.052. doi: 10.1016/j.jmaa.2019.03.052
![]() |
[23] |
Y. Luchko, Initial-boundary-value problems for the one-dimensional time-fractional diffusion equation, Fract. Calc. Appl. Anal., 15 (2012), 141–160. doi: 10.2478/s13540-012-0010-7. doi: 10.2478/s13540-012-0010-7
![]() |
[24] |
B. Jin, B. Li, Z. Zhou, Numerical analysis of nonlinear subdiffusion equations, SIAM J. Numer. Anal., 56 (2018), 1–23. doi: 10.1137/16M1089320. doi: 10.1137/16M1089320
![]() |
[25] |
X. Zheng, H. Wang, Wellposedness and regularity of a variable-order space-time fractional diffusion equation, Anal. Appl., 18 (2020), 615–638. doi: 10.1142/S0219530520500013. doi: 10.1142/S0219530520500013
![]() |
[26] |
H. Fu, M. K. Ng, H. Wang, A divide-and-conquer fast finite difference method for space-time fractional partial differential equation, Comput. Math. Appl., 73 (2017), 1233–1242. doi: 10.1016/j.camwa.2016.11.023. doi: 10.1016/j.camwa.2016.11.023
![]() |
[27] | F. M. Salama, N. H. M. Ali, Fast O(N) hybrid method for the solution of two dimensional time fractional cable equation, Compusoft, 8 (2019), 3453–3461. |
[28] |
F. M. Salama, N. H. M. Ali, Computationally efficient hybrid method for the numerical solution of the 2D time fractional advection-diffusion equation, Int. J. Math. Eng. Manag., 5 (2020), 432–446. doi: 10.33889/IJMEMS.2020.5.3.036. doi: 10.33889/IJMEMS.2020.5.3.036
![]() |
[29] |
X. M. Gu, S. L. Wu, A parallel-in-time iterative algorithm for Volterra partial integro-differential problems with weakly singular kernel, J. Comput. Phys., 417 (2020), 109576. doi: 10.1016/j.jcp.2020.109576. doi: 10.1016/j.jcp.2020.109576
![]() |
[30] |
X. L. Lin, M. K. Ng, H. W. Sun, A multigrid method for linear systems arising from time-dependent two-dimensional space-fractional diffusion equations, J. Comput. Phys., 336 (2017), 69–86. doi: 10.1016/j.jcp.2017.02.008. doi: 10.1016/j.jcp.2017.02.008
![]() |
[31] |
J. Ren, Z. Z. Sun, W. Dai, New approximations for solving the caputo-type fractional partial differential equations, Appl. Math. Model., 40 (2016), 2625–2636. doi: 10.1016/j.apm.2015.10.011. doi: 10.1016/j.apm.2015.10.011
![]() |
[32] |
N. A. Khan, S. Ahmed, Finite difference method with metaheuristic orientation for exploration of time fractional partial differential equations, Int. J. Appl. Comput. Math., 7 (2021), 1–22. doi: 10.1007/s40819-021-01061-y. doi: 10.1007/s40819-021-01061-y
![]() |
[33] |
A. Ahmadian, S. Salahshour, M. Ali-Akbari, F. Ismail, D. Baleanu, A novel approach to approximate fractional derivative with uncertain conditions, Chaos Soliton. Fract., 104 (2017), 68–76. doi: 10.1016/j.chaos.2017.07.026. doi: 10.1016/j.chaos.2017.07.026
![]() |
[34] |
F. M. Salama, N. H. M. Ali, N. N. Abd Hamid, Fast O(N) hybrid Laplace transform-finite difference method in solving 2D time fractional diffusion equation, J. Math. Comput. Sci., 23 (2021), 110–123. doi: 10.22436/jmcs.023.02.04. doi: 10.22436/jmcs.023.02.04
![]() |
[35] |
N. H. M. Ali, L. M. Kew, New explicit group iterative methods in the solution of two dimensional hyperbolic equations, J. Comput. Phys., 231 (2012), 6953–6968. doi: 10.1016/j.jcp.2012.06.025. doi: 10.1016/j.jcp.2012.06.025
![]() |
[36] |
N. H. Mohd Ali, A. Mohammed Saeed, Convergence analysis of the preconditioned group splitting methods in boundary value problems, Abstr. Appl. Anal., 2012 (2012), 867598. doi: 10.1155/2012/867598. doi: 10.1155/2012/867598
![]() |
[37] |
L. M. Kew, N. H. M. Ali, New explicit group iterative methods in the solution of three dimensional hyperbolic telegraph equations, J. Comput. Phys., 294 (2015), 382–404. doi: 10.1016/j.jcp.2015.03.052. doi: 10.1016/j.jcp.2015.03.052
![]() |
[38] |
A. Saudi, J. Sulaiman, Robot path planning using four point-explicit group via nine-point laplacian (4EG9L) iterative method, Procedia Engineering, 41 (2012), 182–188. doi: 10.1016/j.proeng.2012.07.160. doi: 10.1016/j.proeng.2012.07.160
![]() |
[39] |
N. H. M. Ali, A. M. Saeed, Preconditioned modified explicit decoupled group for the solution of steady state navier-stokes equation, Appl. Math. Inform. Sci., 7 (2013), 1837–1844. doi: 10.12785/amis/070522. doi: 10.12785/amis/070522
![]() |
[40] |
M. A. Khan, N. H. M. Ali, N. N. Abd Hamid, A new fourth-order explicit 270 group method in the solution of two-dimensional fractional rayleigh–stokes problem for a heated generalized second-grade fluid, Adv. Diff. Equ., 2020 (2020), 598. doi: 10.1186/s13662-020-03061-6. doi: 10.1186/s13662-020-03061-6
![]() |
[41] |
N. Abdi, H. Aminikhah, A. H. Sheikhani, J. Alavi, M. Taghipour, An efficient explicit decoupled group method for solving two–dimensional fractional Burgers' equation and its convergence analysis, Adv. Math. Phys., 2021 (2021), 6669287. doi: 10.1155/2021/6669287. doi: 10.1155/2021/6669287
![]() |
[42] |
N. Abdi, H. Aminikhah, A. R. Sheikhani, High-order rotated grid point iterative method for solving 2D time fractional telegraph equation and its convergence analysis, Comp. Appl. Math., 40 (2021), 54. doi: 10.1007/s40314-021-01451-4. doi: 10.1007/s40314-021-01451-4
![]() |
[43] |
F. M. Salama, N. H. M. Ali, N. N. Abd Hamid, Efficient hybrid group iterative methods in the solution of two-dimensional time fractional cable equation, Adv. Differ. Equ., 2020 (2020), 257. doi: 10.1186/s13662-020-02717-7. doi: 10.1186/s13662-020-02717-7
![]() |
[44] |
N. Moraca, Bounds for norms of the matrix inverse and the smallest singular value, Linear Algebra Appl., 429 (2008), 2589–2601. doi: 10.1016/j.laa.2007.12.026. doi: 10.1016/j.laa.2007.12.026
![]() |
[45] |
A. Ali, N. H. M. Ali, On skewed grid point iterative method for solving 2d hyperbolic telegraph fractional differential equation, Adv. Differ. Equ., 2019 (2019), 303. doi: 10.1186/s13662-019-2238-6. doi: 10.1186/s13662-019-2238-6
![]() |
1. | Alaa Khadim Mohammed, Salam Jasim Majeed, N. Aldahan, A.J. Ramadhan, Fear induce bistability in an ecoepidemiological model involving prey refuge and hunting cooperation, 2024, 97, 2117-4458, 00150, 10.1051/bioconf/20249700150 | |
2. | Zafer Öztürk, Ali Yousef, Halis Bilgil, Sezer Sorgun, A Fractional-order mathematical model to analyze the stability and develop a sterilization strategy for the habitat of stray dogs, 2024, 14, 2146-5703, 134, 10.11121/ijocta.1418 | |
3. | Asmaa Abdulhussein Aziz, Azhar Abbas Majeed, 2024, 3097, 0094-243X, 080031, 10.1063/5.0209494 | |
4. | Shegaye L. Cheru, Kiros G. Kebedow, Tesfaye T. Ega, Prey-Predator Model of Holling Type II Functional Response with Disease on Both Species, 2024, 0971-3514, 10.1007/s12591-024-00677-y | |
5. | Isaac K. Adu, Fredrick A. Wireko, Mojeeb Al-R. El-N. Osman, Joshua Kiddy K. Asamoah, A fractional order Ebola transmission model for dogs and humans, 2024, 24, 24682276, e02230, 10.1016/j.sciaf.2024.e02230 | |
6. | Ashraf Adnan Thirthar, Prabir Panja, Salam Jasim Majeed, Kottakkaran Sooppy Nisar, Dynamic interactions in a two-species model of the mammalian predator–prey system: The influence of Allee effects, prey refuge, water resources, and moonlights, 2024, 11, 26668181, 100865, 10.1016/j.padiff.2024.100865 | |
7. | Tianxu Wang, Hao Wang, Stoichiometry-Dependent Fear Effect in a Food Chain Model, 2024, 84, 0036-1399, 783, 10.1137/23M1581613 | |
8. | Deniu Yang, The soliton and periodic solutions for the perturbed sine–cosine–Gordon equation, 2024, 47, 0170-4214, 13367, 10.1002/mma.10195 | |
9. | Hanna Rasool Hadi, Azher Abbas Majeed, 2024, 3097, 0094-243X, 080032, 10.1063/5.0209445 | |
10. | Chenyu Liang, Hangjun Zhang, Yancong Xu, Libin Rong, Bifurcation Analysis of a New Aquatic Ecological Model with Aggregation Effect and Harvesting, 2023, 33, 0218-1274, 10.1142/S0218127423501808 | |
11. | Ashraf Adnan Thirthar, Shireen Jawad, Salam J. Majeed, Kottakkaran Sooppy Nisar, Impact of wind flow and global warming in the dynamics of prey–predator model, 2024, 15, 26667207, 100424, 10.1016/j.rico.2024.100424 | |
12. | Ashraf Adnan Thirthar, Zahraa Albatool Mahdi, Prabir Panja, Santanu Biswas, Thabet Abdeljawad, Mutualistic behaviour in an interaction model of small fish, remora and large fish, 2024, 0228-6203, 1, 10.1080/02286203.2024.2392218 | |
13. | Ankur Jyoti Kashyap, Hemanta Kumar Sarmah, Complex Dynamics in a Predator–Prey Model with Fear Affected Transmission, 2024, 0971-3514, 10.1007/s12591-024-00698-7 | |
14. | Mohammed Cherif Bahi, Salma Bahramand, Rashid Jan, Salah Boulaaras, Hassan Ahmad, Rafik Guefaifia, Fractional view analysis of sexual transmitted human papilloma virus infection for public health, 2024, 14, 2045-2322, 10.1038/s41598-024-53696-8 |