
Despite the growing interest in studying the oscillatory behavior of delay differential equations of even-order, odd-order equations have received less attention. In this work, we are interested in studying the oscillatory behavior of two classes of odd-order equations with deviating arguments. We get more than one criterion to check the oscillation in different methods. Our results are an extension and complement to some results published in the literature.
Citation: A. Muhib, I. Dassios, D. Baleanu, S. S. Santra, O. Moaaz. Odd-order differential equations with deviating arguments: asymptomatic behavior and oscillation[J]. Mathematical Biosciences and Engineering, 2022, 19(2): 1411-1425. doi: 10.3934/mbe.2022065
[1] | Shaimaa A. M. Abdelmohsen, D. Sh. Mohamed, Haifa A. Alyousef, M. R. Gorji, Amr M. S. Mahdy . Mathematical modeling for solving fractional model cancer bosom malignant growth. AIMS Biophysics, 2023, 10(3): 263-280. doi: 10.3934/biophy.2023018 |
[2] | Mati ur Rahman, Mehmet Yavuz, Muhammad Arfan, Adnan Sami . Theoretical and numerical investigation of a modified ABC fractional operator for the spread of polio under the effect of vaccination. AIMS Biophysics, 2024, 11(1): 97-120. doi: 10.3934/biophy.2024007 |
[3] | Yasir Nadeem Anjam, Mehmet Yavuz, Mati ur Rahman, Amna Batool . Analysis of a fractional pollution model in a system of three interconnecting lakes. AIMS Biophysics, 2023, 10(2): 220-240. doi: 10.3934/biophy.2023014 |
[4] | Marco Menale, Bruno Carbonaro . The mathematical analysis towards the dependence on the initial data for a discrete thermostatted kinetic framework for biological systems composed of interacting entities. AIMS Biophysics, 2020, 7(3): 204-218. doi: 10.3934/biophy.2020016 |
[5] | Mehmet Yavuz, Kübra Akyüz, Naime Büşra Bayraktar, Feyza Nur Özdemir . Hepatitis-B disease modelling of fractional order and parameter calibration using real data from the USA. AIMS Biophysics, 2024, 11(3): 378-402. doi: 10.3934/biophy.2024021 |
[6] | Larisa A. Krasnobaeva, Ludmila V. Yakushevich . On the dimensionless model of the transcription bubble dynamics. AIMS Biophysics, 2023, 10(2): 205-219. doi: 10.3934/biophy.2023013 |
[7] | Carlo Bianca . Differential equations frameworks and models for the physics of biological systems. AIMS Biophysics, 2024, 11(2): 234-238. doi: 10.3934/biophy.2024013 |
[8] | Mohammed Alabedalhadi, Mohammed Shqair, Ibrahim Saleh . Analysis and analytical simulation for a biophysical fractional diffusive cancer model with virotherapy using the Caputo operator. AIMS Biophysics, 2023, 10(4): 503-522. doi: 10.3934/biophy.2023028 |
[9] | Bertrand R. Caré, Pierre-Emmanuel Emeriau, Ruggero Cortini, Jean-Marc Victor . Chromatin epigenomic domain folding: size matters. AIMS Biophysics, 2015, 2(4): 517-530. doi: 10.3934/biophy.2015.4.517 |
[10] | David Gosselin, Maxime Huet, Myriam Cubizolles, David Rabaud, Naceur Belgacem, Didier Chaussy, Jean Berthier . Viscoelastic capillary flow: the case of whole blood. AIMS Biophysics, 2016, 3(3): 340-357. doi: 10.3934/biophy.2016.3.340 |
Despite the growing interest in studying the oscillatory behavior of delay differential equations of even-order, odd-order equations have received less attention. In this work, we are interested in studying the oscillatory behavior of two classes of odd-order equations with deviating arguments. We get more than one criterion to check the oscillation in different methods. Our results are an extension and complement to some results published in the literature.
Numerous mathematical models have been established to predict and study the biological system. In the past four decades, there have been far-reaching research on improving cell mass production in chemical reactors [1]. The chemostat model is used to understand the mechanism of cell mass growth in a chemostat. A chemostat is an apparatus for continuous culture that contains bacterial populations. It can be used to investigate the cell mass production under controlled conditions. This reactor provides a dynamic system for population studies and is suitable to be used in a laboratory. A substrate is continuously added into the reactor containing the cell mass, which grows by consuming the substrate that enters through the inflow chamber. Meanwhile, the mixture of cell mass and substrate is continuously harvested from the reactor through the outflow chamber. The dynamics in the chemostat can be investigated by using the chemostat model [2].
Ordinary differential equations (ODEs) are commonly used for modelling biological systems. However, most biological systems behaviour has memory effects, and ODEs usually neglect such effects. The fractional-order differential equations (FDEs) are taken into account when describing the behaviour of the systems' equations. A FDEs is a generalisation of the ODEs to random nonlinear order [3]. This equation is more effective because of its good memory, among other advantages [4]–[8]. The errors occuring from the disregarded parameters when modelling of phenomena in real-life also can be reduced. FDEs are also used to efficiently replicate the real nature of various systems in the field of engineering and sciences [9]. In the past few decades, FDEs have been used in biological systems for various studies [5], [6], [10]–[17].
Since great strides in the study of FDEs have been developed, the dynamics in the chemostat can be investigated using the mathematical model of the chemostat in the form of FDEs. Moreover, there have been few studies on the expansion of the chemostat model with fractional-order theory. Thus, we deepen and complete the analysis on the integer-order chemostat model with fractional-order theory and discuss the stability of the equilibrium points of the fractional-order chemostat model. Next, the bifurcation analysis for the fractional-order chemostat model is conducted to identify the bifurcation point that can change the stability of the system. The analysis identifies the values of the fractional-order and the system parameters to ensure the operation of the chemostat is well-controlled.
Recently, there are many approaches to define fractional operators such as by Caputo, Riemann-Liouville, Hadamard, and Grunwald-Letnikov [3]. However, Caputo is often used due to its convenience in various applications [18]. Caputo is also useful to encounter an obstacle where the initial condition is done in the differential of integer-order [19]. In this paper, we applied Caputo derivative to define the system of fractional-order. The Caputo derivative for the left-hand side is defined as
where Г denotes the function of gamma, n is an integer, where
The Adams-type predictor-corrector method is one of the technique that have been proposed for fractional-order differential equations [19]–[21]. The Adams-type predictor-corrector method is a analysis of numerical algorithm that involves two basics steps: predictor and corrector. The predictor formula can be described as
meanwhile the corrector formula can be determined by
The predictor-corrector method is also called as the PECE (Predict, Evaluate, Correct, Evaluate) method [22]. The procedure of the predictor-corrector method can be explained as follows
(i) Calculate the predictor step,
(ii) Evaluate
(iii) Calculate the corrector step,
(iv) Evaluate
The procedure repeatedly predicts and corrects the value until the corrected value becomes a converged number [21]. This method is able to maintain the stability of the properties and has good accuracy. Moreover, this method also has lower computational cost than other methods [23]. The algorithm of the Adams-type predictor-corrector method proposed by [22] is shown in Appendix.
The conditions of stability for integer-order differential equations and fractional-order differential equations are different. Both systems could have the same steady-state points but different stability conditions [24], [25]. The stability condition for fractional-order differential equations can be stated by Theorem 1 and the Routh-Hurwitz stability condition as described by Proposition 1.
Theorem 1
Proposition 1
or
The stability theorem on the fractional-order systems and fractional Routh-Hurwitz stability conditions are introduced to analyze the stability of the model. The fractional Routh-Hurwitz stability conditions is specifically introduced for the eigenvalues of the Jacobian matrix that obtained in quadratic form. The proof of this proposition is shown in Appendix.
Bifurcation can be defined as any sudden change that occurs while a parameter value is varied in the differential equation system and it has a significant influence on the solution [27]. An unstable steady-state may becomes stable and vice versa. A slight changes in the parameter value may change the system's stability. Despite the steady-state point and the eigenvalues of the system of fractional-order are similar as the system of integer-order, the discriminant method used for the stability of the steady-state point is different. Accordingly, the Hopf bifurcation condition in the fractional-order system is slightly different as compared with the integer-order system.
Fractional order α can be selected as the bifurcation parameter in a fractional-order system, but this is not allowed in an integer-order system. The existence of Hopf bifurcation can be stated as in Theorem 2.
Theorem 2
(i) 1. The characteristic equation of chemostat system has a pair of complex conjugate roots,
(ii) Critical value
(iii)
Proof. Condition (i) is not easy to obtain due to the selected parameter's value. However, this condition can be managed under some confined conditions. In fact, the washout steady-state solution of the chemostat model has two negative real roots. The remaining two roots depend on the characteristic of the polynomial from the no-washout steady-state solution.
Condition (ii) can be satisfied with the existence of critical value α* and when
The integer system required p = 0 for the bifurcation's operating condition. For the fractional-order system, the operating condition of the system will change into
In studying the dynamic process of chemostat, the parameters such as
Theorem 3
(i) The characteristic equation of chemostat system has a pair of complex conjugate roots,
(ii) Critical value
(iii)
Proof. This theorem can be proved in the same way as Theorem 2. Therefore, condition (i) can be guaranteed. Condition (ii) can be satisfied with the existence of critical value
For condition (iii), the condition of
Firstly, determine the steady-states, Jacobian matrix and eigenvalues of the fractional-order chemostat model. The stability properties of the fractional-order chemostat model were estimated by using the stability and bifurcation analyses with FDEs by referring to Theorem 1 and Proposition 1. Then, determine the bifurcation point of fractional-order by referring to Theorem 2 and determine the bifurcation point of parameter values by referring to Theorem 3. Next, plot the phase portrait of fractional-order chemostat model by using Adam-types predictor-corrector method to study the dynamic behaviour of the system. Figure 1 depicts the flowchart of this research. This flowchart can be applied to all problems with suitable parameter values.
An integer-order chemostat model that considered a variable yield coefficient and the Monod growth model from [1] is studied in this section. The chemostat system can be written as
with the initial value of
Let Eq. (3.2) equal to zero in order to find the steady-state solutions
By solving Eq. (3.4), the following solutions are obtained
From
Hence, the solutions of steady-state for the chemostat model are
(i) Washout:
(ii) No Washout:
where
The steady-state solutions are physically meaningful if their components are positive. Therefore,
The solution of steady-state in Equation (3.8) represents the washout situation, where the cell mass is wholly removed from the reactor and where the substrate concentration is at the same stock as in the beginning. This state must always be unstable in order to ensure that the cell mass is able to grow in the chemostat. This is because the cell mass will be continuously removed from the chemostat if the washout steady-state is stable. The Jacobian matrix for the washout steady-state solution can be written as
The eigenvalues of this matrix are
The eigenvalues in
The eigenvalues of the Jacobian matrix in terms of the characteristic polynomial are
where
and
The eigenvalues of the no-washout steady-state solution were evaluated with Routh-Hurwitz condition in Proposition 1. Based on the eigenvalues in Eq. (3.15) and by referring to the study by [5], the eigenvalues' condition can be simplified as the following two cases
(i) If b>0 or equivalent to
(ii) If b<0 or equivalent to
Then, if
The parameter values of the fractional-order chemostat model are provided in Table 1. The initial substrate concentration, S0 and ρ were assumed as non-negative values to ensure that the steady-state solutions were physically meaningful. The stability diagram of the steady-state solutions is plotted in Figure 2.
Parameters | Description | Values | Units |
k | Saturation constant | 1.75 | gl−1 |
Q | Dilution rate | 0.02 | l2gr−1 |
µ | Maximum growth rate | 0.3 | h−1 |
γ | Constant in yield coefficient | 0.01 | – |
β | Constant in yield coefficient | 5.25 | lg−1 |
S0 | Input concentration of substrate | 1 | gl−1 |
The washout steady-state solution is stable if Q > 0 and
The steady-state solutions of the fractional-order chemostat model for the parameter values given in Table 1 are
(i) Washout:
(ii) No Washout:
The eigenvalues obtained from the washout steady-state solution are
and the eigenvalues from the no-washout steady-state solution are
Based on Eq. (3.22) to Eq. (3.25), these satisfied the first condition of Hopf bifurcation in Theorem 2. There exists a pair of complex conjugate roots and the other eigenvalues are negative real roots. The transversality condition as the third condition is also satisfied. The eigenvalues of the washout steady-state solution based on the chemostat system is not imaginary, and so there is no existence of Hopf bifurcation in the washout steady-state solution. According to Theorem 2, the critical value of the fractional-order as stated in the second condition can be obtained as
where
Value of p and q are obtained from
The running state of the fractional-order chemostat system when fractional order α at the Hopf bifurcation point is shown. The fractional-order chemostat system changed its stability once Hopf bifurcation occurred. Therefore, we conjecture that the system of fractional-order chemostat may be lost or gain its stability when the fractional order α is less than the Hopf bifurcation point, or
The initial concentration of the substrate, S0, was chosen as the control parameter, while fractional order α was fixed. The solutions of steady-state of the fractional-order chemostat model with S0 as the control parameter are
(i) Washout:
(ii) No Washout:
The eigenvalues obtained from the washout steady-state solution are
and the eigenvalues from the no-washout steady-state solution are
These satisfied the first condition of Hopf bifurcation in Theorem 3. There exist a pair of complex conjugate roots in terms of S0, and the other eigenvalues were negative real roots in terms of S0. The transversality condition as the third condition is also satisfied. According to Theorem 3, the critical value of the fractional order as stated in the second condition can be obtained as follows
By referring to the study by [18], Eq. (3.37) can also be calculated as
From the calculations, the critical value of the initial concentration of the substrate is
The change in the running state when the value of the initial substrate concentration passes through the Hopf bifurcation point is shown. The stability of the fractional-order chemostat system changed once Hopf bifurcation occurred. In
The Hopf bifurcation points of system of fractional-order chemostat and system of integer-order chemostat are different.
The stability analysis of the fractional-order chemostat model was conducted based on the stability theory of FDEs. The integer-order chemostat model was extended to the FDEs. There are two steady-state solutions obtained, which are washout and no-washout steady-state solutions. The Hopf bifurcation of the order of α occured at the solutions of steady-state when the Hopf bifurcation conditions is fulfilled. The results show that the increasing or decreasing the value of α may stabilise the unstable state of the chemostat system. Therefore, the running state of the fractional-order chemostat system is affected by the value of α. The Hopf bifurcation of the initial concentration of the substrate, S0, also occurred when the Hopf bifurcation condition is fulfilled. As the evidence from the phase portrait plots, increase the value of the initial substrate concentration may destabilise the stable state of the chemostat system. The value of the initial substrate should remain at
[1] |
R. P. Agarwal, C. Zhang, T. Li, Some remarks on oscillation of second order neutral differential equations, Appl. Math. Comput., 274 (2016), 178–181. doi: 10.1016/j.amc.2015.10.089. doi: 10.1016/j.amc.2015.10.089
![]() |
[2] |
M. Bohner, S. R. Grace, I. Jadlovska, Oscillation criteria for second-order neutral delay differential equations, Electron. J. Qual. Theo. Differ. Equ., 60 (2017), 1–12. doi: 10.14232/ejqtde.2017.1.60. doi: 10.14232/ejqtde.2017.1.60
![]() |
[3] |
A. Muhib, On oscillation of second-order noncanonical neutral differential equations, J. Inequal. Appl., 2021 (2021), 1–11. doi: 10.1186/s13660-021-02595-x. doi: 10.1186/s13660-021-02595-x
![]() |
[4] |
J. Dzurina, S. R. Grace, I. Jadlovska, T. Li, Oscillation criteria for second-order Emden–Fowler delay differential equations with a sublinear neutral term, Math. Nachr., 293 (2020), 910–922. doi: 10.1002/mana.201800196. doi: 10.1002/mana.201800196
![]() |
[5] |
G. E. Chatzarakis, O. Moaaz, T. Li, B. Qaraad, Some oscillation theorems for nonlinear second-order differential equations with an advanced argument, Adva. Differ. Eq., 2020 (2020), 1–17. doi: 10.1186/s13662-020-02626-9. doi: 10.1186/s13662-020-02626-9
![]() |
[6] |
O. Moaaz, E. M. Elabbasy, B. Qaraad, An improved approach for studying oscillation of generalized Emden–Fowler neutral differential equation, J. Ineq. Appl., 2020 (2020), 1–18. doi: 10.1186/s13660-020-02332-w. doi: 10.1186/s13660-020-02332-w
![]() |
[7] |
O. Moaaz, M. Anis, D. Baleanu, A. Muhib, More effective criteria for oscillation of second-order differential equations with neutral arguments, Mathematics, 8 (2020), 986. doi: 10.3390/math806098. doi: 10.3390/math806098
![]() |
[8] |
S. S. Santra, R. A. El-Nabulsi, Kh. M. Khedher, Oscillation of Second-Order Differential Equations With Multiple and Mixed Delays under a Canonical Operator, Mathematics, 9 (2021), 1323. doi: 10.3390/math9121323. doi: 10.3390/math9121323
![]() |
[9] |
S. S. Santra, Kh. M. Khedher, O. Moaaz, A. Muhib, S-W Yao, Second-order impulsive delay differential systems: necessary and sufficient conditions for oscillatory or asymptotic behavior, Symmetry, 13 (2021), 722. doi: 10.3390/sym13040722. doi: 10.3390/sym13040722
![]() |
[10] |
S. S. Santra, A. K. Sethi, O. Moaaz, Kh. M. Khedher, S-W. Yao, New oscillation theorems for second-order differential equations with canonical and non canonical operator via Riccati transformation, Mathematics, 9 (2021), 1111. doi: 10.3390/math9101111. doi: 10.3390/math9101111
![]() |
[11] |
A. Alghamdi, C. Cesarano, B. Almarri, O. Bazighifan, Symmetry and its importance in the oscillation of solutions of differential equations, Symmetry, 13 (2021), 650. doi: 10.3390/sym13040650. doi: 10.3390/sym13040650
![]() |
[12] |
O. Moaaz, Ch. Park, A. Muhib, O. Bazighifan, Oscillation criteria for a class of even-order neutral delay differential equations, J. Appl. Math. Comput., 63 (2020), 607–617. doi: 10.1007/s12190-020-01331-w. doi: 10.1007/s12190-020-01331-w
![]() |
[13] |
O. Moaaz, C. Cesarano, A. Muhib, Some new oscillation results for fourth-order neutral differential equations, Eur. J. Pure Appl. Math., 13 (2020), 185–199. doi: 10.29020/nybg.ejpam.v13i2.36. doi: 10.29020/nybg.ejpam.v13i2.36
![]() |
[14] |
O. Bazighifan, T. Abdeljawad, Q. M. Al-Mdallal, Differential equations of even-order with p-Laplacian like operators: qualitative properties of the solutions, Adv. Differ. Equ., 2021 (2021), 96. doi: 10.1186/s13662-021-03254-7. doi: 10.1186/s13662-021-03254-7
![]() |
[15] |
O. Bazighifan, F. Minhos, O. Moaaz, Sufficient conditions for oscillation of fourth-order neutral differential equations with distributed deviating arguments, Axioms, 9 (2020), 39. doi: 10.3390/axioms9020039. doi: 10.3390/axioms9020039
![]() |
[16] |
O. Moaaz, S. Furuichi, A. Muhib, New comparison theorems for the nth order neutral differential equations with delay inequalities, Mathematics, 8 (2020), 454. doi: 10.3390/math8030454. doi: 10.3390/math8030454
![]() |
[17] |
O. Moaaz, R. A. El-Nabulsi, O. Bazighifan, Oscillatory behavior of fourth-order differential equations with neutral delay, Symmetry, 12 (2020), 371. doi: 10.3390/sym12030371. doi: 10.3390/sym12030371
![]() |
[18] |
C. Park, O. Moaaz, O. Bazighifan, Oscillation results for higher order differential equations, Axioms, 9 (2020), 14. doi: 10.3390/axioms9010014. doi: 10.3390/axioms9010014
![]() |
[19] |
M. Aktas, A. Tiryaki, A. Zafer, Oscillation criteria for third-order nonlinear functional differential equations, Appl. Math. Lett., 23 (2010), 756–762. doi: 10.1016/j.aml.2010.03.003. doi: 10.1016/j.aml.2010.03.003
![]() |
[20] |
B. Baculikova, E. M. Elabbasy, S. H. Saker, J. Dzurina, Oscillation criteria for third-order nonlinear differential equations, Math. Slovaca, 58 (2008), 201–220. doi: 10.2478/s12175-008-0068-1. doi: 10.2478/s12175-008-0068-1
![]() |
[21] |
O. Moaaz, J. Awrejcewicz, A. Muhib, Establishing new criteria for oscillation of odd-order nonlinear differential equations, Mathematics, 8 (2020), 607–617. doi: 10.3390/math8060937. doi: 10.3390/math8060937
![]() |
[22] |
O. Moaaz, I. Dassios, W. Muhsin, A. Muhib, Oscillation Theory for Non-Linear Neutral Delay Differential Equations of Third Order, Appl. Sci., 10 (2020), 4855. doi: 10.3390/app10144855. doi: 10.3390/app10144855
![]() |
[23] |
O. Moaaz, E. M. Elabbasy, E. Shaaban, Oscillation criteria for a class of third order damped differential equations, Arab J. Math. Sci., 24 (2018), 16–30. doi: 10.1016/j.ajmsc.2017.07.001. doi: 10.1016/j.ajmsc.2017.07.001
![]() |
[24] |
O. Moaaz, D. Baleanu, A. Muhib, New aspects for non-existence of kneser solutions of neutral differential equations with odd-order, Mathematics, 8 (2020), 494. doi: 10.3390/math8040494. doi: 10.3390/math8040494
![]() |
[25] | S. H. Saker, J. Dzurina, On the oscillation of certain class of third-order nonlinear delay differential equations, Math. Bohem., 135 (2010), 225–237. |
[26] |
R. P. Agarwal, M. Bohner, T. Li, C. Zhang, Oscillation of third-order nonlinear delay differential equations, Taiwanese J. Math., 17 (2013), 545–558. doi: 10.11650/tjm.17.2013.2095. doi: 10.11650/tjm.17.2013.2095
![]() |
[27] |
T. Li, Yu. V. Rogovchenko, On asymptotic behavior of solutions to higher-order sublinear Emden-Fowler delay differential equations, Appl. Math. Lett., 67 (2017), 53–59. doi: 10.1016/j.aml.2016.11.007. doi: 10.1016/j.aml.2016.11.007
![]() |
[28] |
C. Zhang, R. P. Agarwal, T. Li, Oscillation and asymptotic behavior of higher-order delay differential equations with fanxiexian_myfhpfanxiexian_myfh-Laplacian like operators, J. Math. Anal. Appl., 409 (2014), 1093–1106. doi: 10.1016/j.jmaa.2013.07.066. doi: 10.1016/j.jmaa.2013.07.066
![]() |
[29] |
C. Zhang, T. Li, B. Sun, E. Thandapani, On the oscillation of higher-order half-linear delay differential equations, Appl. Math. Lett., 24 (2011), 1618–1621. doi: 10.1016/j.aml.2011.04.015. doi: 10.1016/j.aml.2011.04.015
![]() |
[30] |
B. Baculikova, J. Dzurina, Oscillation of third-order neutral differential equations, Math. Comput. Mode., 52 (2010), 215–226. doi: 10.1016/j.mcm.2010.02.011. doi: 10.1016/j.mcm.2010.02.011
![]() |
[31] |
T. Li, Y. V. Rogovchenko, On the asymptotic behavior of solutions to a class of third-order nonlinear neutral differential equations, Appl. Math. Lett., (2020), 106293. doi: 10.1016/j.aml.2020.106293. doi: 10.1016/j.aml.2020.106293
![]() |
[32] |
D. Lackova, The asymptotic properties of the solutions of the n-th order functional neutral differential equations, Comput. Appl. Math., 146 (2003), 385–392. doi: 10.1016/S0096-3003(02)00590-8. doi: 10.1016/S0096-3003(02)00590-8
![]() |
[33] |
M. Bohner, T. S. Hassan, T. Li, Fite-Hille-Wintner-type oscillation criteria for second-order half-linear dynamic equations with deviating arguments, Indag. Math. (N.S.), 29 (2018), 548–560. doi: 10.1016/j.indag.2017.10.006. doi: 10.1016/j.indag.2017.10.006
![]() |
[34] |
T. Li, N. Pintus, G. Viglialoro, Properties of solutions to porous medium problems with different sources and boundary conditions, Z. Angew. Math. Phys., 70 (2019), 1–18. doi: 10.1007/s00033-019-1130-2. doi: 10.1007/s00033-019-1130-2
![]() |
[35] | T. Li, G. Viglialoro, Boundedness for a nonlocal reaction chemotaxis model even in the attraction-dominated regime, Differ. Integral Equ., 34 (2021), 315–336. |
[36] |
B. Baculikova, J. Dzurina, On the oscillation of odd order advanced differential equations, Bound. Value Probl., 214 (2014), 214. doi: 10.1186/s13661-014-0214-3. doi: 10.1186/s13661-014-0214-3
![]() |
[37] | R. P. Agarwal, S. R. Grace, D. O'Regan, Oscillation Theory for difference and functional differential equations, Marcel Dekker, Kluwer Academic, Dordrecht, 2000. doi: 10.1007/978-94-015-9401-1. |
1. | Xiaomeng Ma, Zhanbing Bai, Sujing Sun, Stability and bifurcation control for a fractional-order chemostat model with time delays and incommensurate orders, 2022, 20, 1551-0018, 437, 10.3934/mbe.2023020 |
Parameters | Description | Values | Units |
k | Saturation constant | 1.75 | gl−1 |
Q | Dilution rate | 0.02 | l2gr−1 |
µ | Maximum growth rate | 0.3 | h−1 |
γ | Constant in yield coefficient | 0.01 | – |
β | Constant in yield coefficient | 5.25 | lg−1 |
S0 | Input concentration of substrate | 1 | gl−1 |