An improved multi-strategy beluga whale optimization for global optimization problems

Hongmin Chen; Zhuo Wang; Di Wu; Heming Jia; Changsheng Wen; Honghua Rao; Laith Abualigah; Hongmin Chen; Zhuo Wang; Di Wu; Heming Jia; Changsheng Wen; Honghua Rao; Laith Abualigah

doi:10.3934/mbe.2023592

Mathematical Biosciences and Engineering

2023, Volume 20, Issue 7: 13267-13317. doi: 10.3934/mbe.2023592

Previous Article Next Article

Research article

An improved multi-strategy beluga whale optimization for global optimization problems

1.
School of Information Engineering, Sanming University, Sanming 365004, China
2.
School of Education and Music, Sanming University, Sanming 365004, China
3.
Prince Hussein Bin Abdullah College for Information Technology, Al Al-Bayt University, Mafraq 130040, Jordan
4.
College of Engineering, Yuan Ze University, Taoyuan, Taiwan
5.
Hourani Center for Applied Scientific Research, Al-Ahliyya Amman University, Amman 19328, Jordan
6.
MEU Research Unit, Middle East University, Amman 11831, Jordan
7.
Applied Science Research Center, Applied Science Private University, Amman 11931, Jordan
8.
School of Computer Sciences, Universiti Sains Malaysia, Pulau Pinang 11800, Malaysia

Academic Editor: Weining Wu

Received: 07 April 2023 Revised: 22 May 2023 Accepted: 23 May 2023 Published: 09 June 2023

This paper presents an improved beluga whale optimization (IBWO) algorithm, which is mainly used to solve global optimization problems and engineering problems. This improvement is proposed to solve the imbalance between exploration and exploitation and to solve the problem of insufficient convergence accuracy and speed of beluga whale optimization (BWO). In IBWO, we use a new group action strategy (GAS), which replaces the exploration phase in BWO. It was inspired by the group hunting behavior of beluga whales in nature. The GAS keeps individual belugas whales together, allowing them to hide together from the threat posed by their natural enemy, the tiger shark. It also enables the exchange of location information between individual belugas whales to enhance the balance between local and global lookups. On this basis, the dynamic pinhole imaging strategy (DPIS) and quadratic interpolation strategy (QIS) are added to improve the global optimization ability and search rate of IBWO and maintain diversity. In a comparison experiment, the performance of the optimization algorithm (IBWO) was tested by using CEC2017 and CEC2020 benchmark functions of different dimensions. Performance was analyzed by observing experimental data, convergence curves, and box graphs, and the results were tested using the Wilcoxon rank sum test. The results show that IBWO has good optimization performance and robustness. Finally, the applicability of IBWO to practical engineering problems is verified by five engineering problems.

Keywords:

Citation: Hongmin Chen, Zhuo Wang, Di Wu, Heming Jia, Changsheng Wen, Honghua Rao, Laith Abualigah. An improved multi-strategy beluga whale optimization for global optimization problems[J]. Mathematical Biosciences and Engineering, 2023, 20(7): 13267-13317. doi: 10.3934/mbe.2023592

Related Papers:

[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

Abstract

1. Introduction

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.

2. Materials and methods

2.1. Definition of Caputo derivative

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

$D_{t}^{α} f (t) = \frac{1}{Γ (n - α)} \int_{0}^{t} \frac{f^{(n)} τ}{{(t - τ)}^{α - n + 1}} d τ,$ (2.1)

where Г denotes the function of gamma, n is an integer, where $n - 1 < α < n$ [18].

2.2. Adams-type predictor corrector method

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

$\begin{array}{l} y_{h}^{P} (t_{n + 1}) = \sum_{k = 0}^{α - 1} \frac{t_{n + 1}^{k}}{k!} y_{0}^{(k)} + \frac{1}{Γ (α)} \sum_{j = 0}^{n} b_{j, n + 1} f (t_{j}, y_{h} (t_{j})), \end{array}$ (2.2)

meanwhile the corrector formula can be determined by

$y_{h} (t_{n + 1}) = \sum_{k = 0}^{α - 1} \frac{t_{n + 1}^{k}}{k!} + \frac{h^{α}}{Γ (α + 2)} f (t_{n + 1}, y_{h}^{P} (t_{n + 1})) + \frac{h^{α}}{Γ (α + 2)} \sum_{j = 0}^{n} f (t_{j}, y_{h} (t_{j})) .$ (2.3)

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, $y_{h}^{P} (t_{n + 1})$ in Eq. (2.2).

(ii) Evaluate $f (t_{n + 1}, y_{h}^{P} (t_{n + 1}))$ .

(iii) Calculate the corrector step, $y_{h} (t_{n + 1})$ in Eq. (2.3).

(iv) Evaluate $f (t_{n + 1}, y_{h} (t_{n + 1}))$ .

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.

2.3. Stability analysis of fractional-order system

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 [5], [6], [26]. The commensurate system of fractional-order where $x \in ℝ$ and $0 < α < 1$ is locally asymptotically stable if the eigenvalues of the Jacobian matrix evaluated at the steady-state point is satisfied by

$| a r g (λ) | > \frac{α π}{2} .$ (2.4)

Proposition 1 [5], [6], [26]. Suppose the characteristic polynomial is $P (λ) = λ^{2} + b λ + c$ of the Jacobian matrix which evaluated at the steady-state. The eigenvalues of the Jacobian matrix will satisfy Eq. (2.4) in Theorem 1 if

$b > 0, c > 0,$ (2.5)

$b < 0, 4 c > b^{2}, | t a n^{- 1} (\frac{\sqrt{4 c - b^{2}}}{b}) | > \frac{α π}{2} .$ (2.6)

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.

2.4. Hopf bifurcation analysis of fractional-order system

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.

2.4.1. Hopf bifurcation analysis of fractional order α

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 [28]. Assume α* as the critical value of the fractional-order. When bifurcation parameter α passes over critical value α*, which is $α^{*} \in (0,1)$ , Hopf bifurcation occurs at the steady-state point if the following conditions are satisfied

(i) 1. The characteristic equation of chemostat system has a pair of complex conjugate roots, $λ_{1,2} = p \pm i q$ , while the other eigenvalues are negative real roots.

(ii) Critical value $m (α^{*}) = \frac{α^{*} π}{2} - m i n | a r g (λ) | = 0$ .

(iii) $\frac{d m (α)}{d α} |_{α = α^{*}} \neq 0$ (condition of transversality).

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 $a r g (λ)$ is equivalent to $a r c t a n (\frac{q}{p})$ . Thus, the solution of critical value $m (α^{*})$ can be written as

$α^{*} = \frac{α^{*} π}{2} - a r c t a n (\frac{q}{p}) = 0, α^{*} \in (0,1) .$ (2.7)

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 $m (α^{*}) = \frac{α^{*} π}{2} - m i n | a r g (x) | = 0$ . For condition (iii), the condition of $m (α)$ changes when bifurcation parameter $α$ passes over critical value α*. For example, the steady-state point is asymptotically stable for $0 < α < α^{*}$ and unstable when $α < α^{*} < 1$ . Thus, Hopf bifurcation exists at $α = α^{*}$ .

2.4.2. Hopf bifurcation analysis of parameters value

In studying the dynamic process of chemostat, the parameters such as $Q, S_{0}, µ, k, γ$ and β are usually used as the bifurcation parameter since these parameters have significant effects on the dynamic process of the system of fractional-order and integer-order. Fractional-order α is considered fixed and the initial substrate concentration S₀ is studied as the control parameter. The existence of the Hopf bifurcation can be stated as in Theorem 3.

Theorem 3 [28]. Assume $S_{0}^{*}$ as the critical value of the fractional-order. When bifurcation parameter S₀ passes over critical value $S_{0}^{*}$ , Hopf bifurcation occurs at the steady-state point if the following conditions are satisfied

(i) The characteristic equation of chemostat system has a pair of complex conjugate roots, $λ_{1,2} = p (S_{0}) \pm i q (S_{0})$ , while the other eigenvalues are negative real roots.

(ii) Critical value $m (S_{0}^{*}) = \frac{α π}{2} - m i n | a r g (λ (S_{0}^{*})) | = 0$ .

(iii) $\frac{d m (S_{0}^{*})}{d (S_{0}^{*})} |_{S_{0} = S_{0}^{*}} \neq 0$ (condition of transversality).

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 $S_{0}^{*}$ and when $a r g [λ (S_{0}^{*})]$ is equivalent to $a r c t a n [\frac{q (S_{0}^{*})}{p (S_{0}^{*})}]$ . Thus, the solution of critical value $m (S_{0}^{*})$ can be written as

$S_{0}^{*} = \frac{α π}{2} - a r c t a n (\frac{q (S_{0}^{*})}{p (S_{0}^{*})}) = 0.$ (2.8)

For condition (iii), the condition of $m (S_{0}^{*})$ changes when bifurcation parameter S₀ passes over critical value $S_{0}^{*}$ . For example, the steady-state point is asymptotically stable when $0 < S_{0} < S_{0}^{*}$ and unstable when $S_{0} < S_{0}^{*} < 1$ . Thus, Hopf bifurcation exists at $S_{0} = S_{0}^{*}$ [28].

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.

Figure 1. Mathematical analysis of fractional-order chemostat model.

DownLoad: Full-Size Img PowerPoint

3. Results and Discussion

3.1. Fractional-order chemostat model

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

$\begin{array}{l} \frac{d S}{d t} & = Q (S_{0} - S) - \frac{µ S X}{(k + S) (γ + β S)}, \\ \frac{d X}{d t} & = Q (- X) + \frac{µ S X}{k + S}, \end{array}$ (3.1)

with the initial value of $X_{0} = 0$ , where the sterile feed case was assumed. The integer-order chemostat system of Eq. (3.1) is extended to the fractional-order differential equation

$\begin{array}{l} \frac{d^{α} S}{d t^{α}} & = Q (S_{0} - S) - \frac{µ S X}{(k + S) (γ + β S)}, \\ \frac{d^{α} X}{d t^{α}} & = Q (- X) + \frac{µ S X}{k + S} . \end{array}$ (3.2)

Let Eq. (3.2) equal to zero in order to find the steady-state solutions

$Q (S_{0} - S) - \frac{µ S X}{(k + S) (γ + β S)} = 0,$ (3.3)

$Q (- X) + \frac{µ S X}{k + S} = 0.$ (3.4)

By solving Eq. (3.4), the following solutions are obtained

$S^{*} = \frac{k Q}{µ - Q},$ (3.5)

$X^{*} = 0.$ (3.6)

From Eq. (3.3), if $X^{*} = 0$ , then $S^{*} = S_{0}$ . If $S^{*} = \frac{k Q}{µ - Q}$ , then

$X^{*} = \frac{(k Q + Q S_{0} - S_{0} µ) (- Q γ + k Q β + γ µ)}{{(µ - Q)}^{2}} .$ (3.7)

Hence, the solutions of steady-state for the chemostat model are

(i) Washout:

$(S_{0}^{*}, X_{0}^{*}) = (S_{0},0) .$ (3.8)

(ii) No Washout:

$(S_{1}^{*}, X_{1}^{*}) = (ρ, (S_{0} - ρ) (γ + β ρ)),$ (3.9)

where

$\begin{array}{l} ρ = \frac{k Q}{µ - Q} . \end{array}$

The steady-state solutions are physically meaningful if their components are positive. Therefore, $S_{0} > 0$ for the washout steady-state solution exists by biological meaning. The no-washout steady-state solution will only exist when $0 < ρ < S_{0}$ . The Jacobian matrix of no washout steady-state as in Eq. (3.10) can be used to investigate the stability properties of the fractional-order chemostat model.

$J = [\begin{matrix} - Q + \frac{X (- k γ + S^{2} β) µ}{{(k + S)}^{2} {(γ + S β)}^{2}} & \frac{S µ}{(k + S) (γ + S β)} \\ \frac{k X µ}{{(k + S)}^{2}} & - Q + \frac{S µ}{k + S} \end{matrix}] .$ (3.10)

3.2. Stability of washout and no-washout steady-state solution

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

$J = [\begin{matrix} - Q & - \frac{S_{0} µ}{(k + S_{0}) (γ + S_{0} β)} \\ 0 & - Q + \frac{S_{0} µ}{k + S_{0}} \end{matrix}] .$ (3.11)

The eigenvalues of this matrix are

$λ_{1} = - Q,$ (3.12)

$λ_{2} = \frac{- k Q - Q S_{0} + S_{0} µ}{k + S_{0}} .$ (3.13)

The eigenvalues in Eq. (3.12) and Eq. (3.13) are real. The washout steady-state solution is stable if Q>0 and $ρ > S_{0}$ where $ρ = \frac{k Q}{µ - Q}$ . The steady-state solution in Eq. (3.9) represents the no-washout situation. No-washout situation is where the cell mass is not removed and stay growth in the chemostat.This state is important. The steady-state solution is substituted into the Jacobian matrix in Eq. (3.10) and can be written as

$J = [\begin{matrix} - Q + \frac{µ (S_{0} - ρ) (β ρ^{2} - k γ)}{{(k + ρ)}^{2} (β ρ - γ)} & - \frac{µ ρ}{(k + ρ) (β ρ + γ)} \\ \frac{k µ (S_{0} - ρ) (β ρ - γ)}{{(k + ρ)}^{2}} & - Q + \frac{µ ρ}{k + ρ} \end{matrix}] .$ (3.14)

The eigenvalues of the Jacobian matrix in terms of the characteristic polynomial are

$P (λ) = λ^{2} + b λ + c,$ (3.15)

where

$b = 2 Q + \frac{µ (S_{0} - 2 ρ)}{(k + ρ)} + \frac{µ ρ (ρ - S_{0})}{{(k + ρ)}^{2}} + \frac{β µ ρ (ρ - S_{0})}{(k + ρ) (γ + β ρ)},$ (3.16)

and

$\begin{array}{l} c = & Q^{2} + \frac{Q µ (S_{0} - 2 ρ)}{(k + ρ)} + \frac{(µ ρ - S_{0} µ) (µ ρ - Q ρ)}{{(k + ρ)}^{2}} + \frac{µ^{2} ρ^{2} (S_{0} - ρ)}{{(k + ρ)}^{3}} \\ - \frac{Q β µ ρ (S_{0} - ρ)}{(k + ρ) (γ + β ρ)} + \frac{(S_{0} µ - µ ρ) (µ ρ γ + 2 µ β ρ^{2})}{{(k + ρ)}^{2} (γ + β ρ)} \\ + \frac{(µ ρ - S_{0} µ) (µ ρ^{2} γ + µ β ρ^{3})}{{(k + ρ)}^{3} (γ + β ρ)} . \end{array}$ (3.17)

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 $\frac{γ}{β} > P_{1}$ , the no-washout steady-state solution of the system in Eq. (3.2) is asymptotically stable. P₁ can be written as

$P_{1} = - \frac{µ ρ (ρ - S_{0})}{2 Q (k + ρ)} - \frac{ρ (ρ - S_{0})}{(S_{0} - ρ)} - \frac{(k + ρ)}{µ ρ (ρ - S_{0})} - ρ,$ (3.18)

(ii) If b<0 or equivalent to $\frac{γ}{β} < P_{1}$ and $t a n^{- 1} (\frac{\sqrt{4 c - b^{2}}}{b}) > \frac{α π}{2}$ , the no-washout steady-state solution of the system in Eq. (3.2) is asymptotically stable. The condition of $t a n^{- 1} (\frac{\sqrt{4 c - b^{2}}}{b}) > \frac{α π}{2}$ is also equivalent to $4 c o s^{2} (\frac{α π}{2}) c > b^{2}$ , which can be simplified as $\frac{γ}{β} > P_{2}$ . Then, this case can be concluded and written as $P_{2} < \frac{γ}{β} < P_{1}$ where

$P_{2} = \frac{µ ρ (ρ - S_{0})}{2 c o s (\frac{α π}{2}) \sqrt{c} (k + ρ)} - \frac{µ ρ (ρ - S_{0})}{2 Q (k + ρ)} - \frac{ρ (ρ - S_{0})}{(S_{0} - ρ)} - k - 2 ρ .$ (3.19)

Then, if $P_{2} < \frac{γ}{β} < P_{1}$ , the no-washout steady-state solution of the system in Eq. (3.2) is asymptotically stable.

The parameter values of the fractional-order chemostat model are provided in Table 1. The initial substrate concentration, S₀ 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.

Table 1. Parameter values.

Parameters	Description	Values	Units
k	Saturation constant	1.75	gl⁻¹
Q	Dilution rate	0.02	l²gr⁻¹
µ	Maximum growth rate	0.3	h⁻¹
γ	Constant in yield coefficient	0.01	–
β	Constant in yield coefficient	5.25	lg⁻¹
S₀	Input concentration of substrate	1	gl⁻¹

| Show Table

DownLoad: CSV

Figure 2. Stability diagram of the steady-state solutions when $α = 1$ .

DownLoad: Full-Size Img PowerPoint

The washout steady-state solution is stable if Q > 0 and $ρ > S_{0}$ . From Figure 2, it shows that the unstable solution of washout steady-state, as the eigenvalues did not fulfil the condition of $ρ > S_{0}$ . By choosing the appropriate parameter values, the unstable washout steady-state solution could ensure that the washout condition does not occur in the chemostat. Meanwhile, the solution of no-washout steady-state is stable.

3.3. Hopf bifurcation analysis of the order of fractional-order system

The steady-state solutions of the fractional-order chemostat model for the parameter values given in Table 1 are

(i) Washout:

$(S_{0}^{*}, X_{0}^{*}) = (1,0),$ (3.20)

(ii) No Washout:

$(S_{1}^{*}, X_{1}^{*}) = (\frac{1}{8}, \frac{3731}{6400}) .$ (3.21)

The eigenvalues obtained from the washout steady-state solution are

$λ_{1} = - \frac{1}{50},$ (3.22)

$λ_{2} = - \frac{49}{550},$ (3.23)

and the eigenvalues from the no-washout steady-state solution are

$λ_{1} = \frac{- 2552 + \sqrt{411098126 i}}{399750},$ (3.24)

$λ_{2} = \frac{- 2552 - \sqrt{411098126 i}}{399750},$ (3.25)

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

$m (α^{*}) = \frac{α^{*} π}{2} - m i n | a r g (λ) | = 0,$ (3.26)

$α^{*} = \frac{2}{π} m i n | a r g (λ) |,$ (3.27)

where

$a r g (λ) = a r c t a n (\frac{q}{p}),$ (3.28)

$α^{*} = \frac{2}{π} a r c t a n (\frac{q}{p}) = \frac{2}{π} a r c t a n (\frac{\frac{\sqrt{411098126}}{399750}}{\frac{2552}{399750}}) = 0.9202904711 \approx 0.9.$ (3.29)

Value of p and q are obtained from Eq. (3.24) and Eq. (3.25) by assuming parameter value in Table 1. Hence, when $α^{*} = 0.9$ , the chemostat system in Eq. (3.2) shows Hopf bifurcation, at which the system stability would be altered.

Figure 3 is plotted to determine the dynamic behaviour at the Hopf bifurcation point. The phase portrait diagrams of cell mass concentration against substrate concentration are plotted for values of order of the fractional is $α = 0.9$ .

Figure 3. Phase portrait plot of fractional-order chemostat system with $α = 0.9$ .

DownLoad: Full-Size Img PowerPoint

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 $α < 0.9$ or otherwise. This shows that increasing or decreasing the value of α may destabilise the stable state of the chemostat system. Therefore, these results show that the running state of the fractional-order chemostat system is affected by the value of α.

3.4. Hopf bifurcation analysis of initial concentration of substrate

The initial concentration of the substrate, S₀, was chosen as the control parameter, while fractional order α was fixed. The solutions of steady-state of the fractional-order chemostat model with S₀ as the control parameter are

(i) Washout:

$(S_{0}^{*}, X_{0}^{*}) = (S_{0},0),$ (3.30)

(ii) No Washout:

$(S_{1}^{*}, X_{1}^{*}) = (\frac{1}{8}, \frac{533 (- 1 + 8 S_{0})}{6400}) .$ (3.31)

The eigenvalues obtained from the washout steady-state solution are

$λ_{1} = - \frac{1}{50},$ (3.32)

$λ_{2} = - \frac{7 (1 - 8 S_{0})}{50 (- 7 - 4 S_{0}},$ (3.33)

and the eigenvalues from the no-washout steady-state solution are

$λ_{1} = - 4204 + 1652 S_{0} + \frac{\sqrt{2} \sqrt{1364552 S_{0}^{2} - 245579768 S_{0} + 38666153}}{399750} i,$ (3.34)

$λ_{2} = - 4204 + 1652 S_{0} - \frac{\sqrt{2} \sqrt{1364552 S_{0}^{2} - 245579768 S_{0} + 38666153}}{399750} i .$ (3.35)

These satisfied the first condition of Hopf bifurcation in Theorem 3. There exist a pair of complex conjugate roots in terms of S₀, and the other eigenvalues were negative real roots in terms of S₀. 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

$m (S_{0}^{*}) = \frac{α^{*} π}{2} - m i n | a r g (λ) | = 0.$ (3.36)

By referring to the study by [18], Eq. (3.37) can also be calculated as

$\frac{q^{'} (S_{0}^{*}) p (S_{0}^{*}) - q (S_{0}^{*}) p^{'} (S_{0}^{*})}{q^{2} (S_{0}^{*}) + p^{2} (S_{0}^{*})} \neq 0.$ (3.37)

From the calculations, the critical value of the initial concentration of the substrate is $S_{0} = 2.54$ . When $S_{0} = 2.54$ , the chemostat system shows Hopf bifurcation, at which the stability of the system would be altered.

Figure 4 presents the phase portrait diagrams of concentration of cell mass against concentration of substrate when $α = 1$ for different values of the initial concentration of the substrate, which are $S_{0} = 2$ , $S_{0} = 2.54$ and $S_{0} = 3.5$ .

Figure 4. Phase portrait plot of chemostat system with $α = 1$ (a) $S_{0} = 2$ (b) $S_{0} = 2.54$ and (c) $S_{0} = 3.5$ .

DownLoad: Full-Size Img PowerPoint

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 Figure 3(a), the fractional-order chemostat system is in a stable state when the initial substrate concentration value is less than the Hopf bifurcation point, or $S_{0} < 2.54$ . Meanwhile, when the value of the initial substrate concentration passes through the Hopf bifurcation point, or $S_{0} \geq 2.54$ , the fractional-order chemostat system lost its stability. This shows that increasing the value of the initial substrate concentration may destabilise the stable state of the chemostat system. These results show that the running state of the fractional-order chemostat system is affected by the value of the initial substrate concentration. In real-life application, the value of the initial substrate should remain at $S_{0} \geq 2.54$ to ensure that the chemostat system is at the unstable state. This is because the unstable state is suitable for the production of cell mass [1]. Unstable state means the system always move away after small disturbance, so the system must be at the unstable state because there will be a change in amount of cell mass production.

Figure 5 depicts the phase portrait diagrams of cell mass concentration against substrate concentration when $α = 0.9$ for different values of the initial concentration of the substrate, which are $S_{0} = 2$ , $S_{0} = 2.54$ and $S_{0} = 3.5$ .

Figure 5. Phase portrait plot of chemostat system with $α = 0.9$ (a) $S_{0} = 2$ (b) $S_{0} = 2.54$ and (c) $S_{0} = 3.5$ .

DownLoad: Full-Size Img PowerPoint

The Hopf bifurcation points of system of fractional-order chemostat and system of integer-order chemostat are different. Figure 4 shows the fractional-order chemostat system at a stable state for all values of the initial substrate concentration when $α = 0.9$ . The chemostat system destabilised the stable state when the initial substrate concentration value is $S_{0} \geq 2.54$ , as shown in Figure 3(b) and Figure 3(c). This shows that the dynamic behaviour of the fractional-order chemostat system is different compared with the integer-order chemostat system. In actual application, the value of the initial substrate should remain at $S_{0} \geq 2.54$ to ensure that the chemostat system can be well controlled in order to be suitable for cell mass production.

4. Conclusions

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, S₀, 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 $S_{0} \geq 2.54$ to ensure that the chemostat system is at the unstable state since the unstable state is suitable for the production of cell mass. These dynamical analyses are important to provide suitable values of the fractional-order and the parameters in order to ensure the controllability and stability of the chemostat to suit the actual chemostat environment.

References

[1]	E. G. Talbi, Metaheuristics: from Design to Implementation, John Wiley & Sons, 2009. https://doi.org/10.1002/9780470496916
[2]	X. S. Yang, Nature-inspired optimization algorithms: Challenges and open problems, J. Comput. Sci., 46 (2020), 101104. https://doi.org/10.1016/j.jocs.2020.101104 doi: 10.1016/j.jocs.2020.101104
[3]	M. H. Nadimi-Shahraki, S. Taghian, S. Mirjalili, An improved grey wolf optimizer for solving engineering problems, Expert Syst. Appl., 166 (2021), 113917. https://doi.org/10.1016/j.eswa.2020.113917 doi: 10.1016/j.eswa.2020.113917
[4]	M. H. Nadimi-Shahraki, S. Taghian, S. Mirjalili, H. Faris, MTDE: An effective multi-trial vector-based differential evolution algorithm and its applications for engineering design problems, Appl. Soft Comput., 97 (2020), 106761. https://doi.org/10.1016/j.asoc.2020.106761 doi: 10.1016/j.asoc.2020.106761
[5]	S. R. Zhao, Y. L. Wu, S. Tan, J. R. Wu, Z. S. Cui, Y. G. Wang, QQLMPA: A quasi-opposition learning and Q-learning based marine predators algorithm, Expert Syst. Appl., 213 (2023), 119246. https://doi.org/10.1016/j.eswa.2022.119246 doi: 10.1016/j.eswa.2022.119246
[6]	C. T. Zhong, G. Li, Z. Zeng, Beluga whale optimization: A novel nature-inspired metaheuristic algorithm, Knowledge-Based Syst., 251 (2022), 109215. https://doi.org/10.1016/j.knosys.2022.109215 doi: 10.1016/j.knosys.2022.109215
[7]	J. Kennedy, R. Eberhart, Particle swarm optimization, in Proceedings of ICNN'95-International Conference on Neural Networks, IEEE, 4 (1995), 1942–1948. https://doi.org/10.1109/ICNN.1995.488968
[8]	S. Mirjalili, S. M Mirjalili, A. Lewis, Grey wolf optimizer, Adv. Eng. Software, 69 (2014), 46–61. https://doi.org/10.1016/j.advengsoft.2013.12.007 doi: 10.1016/j.advengsoft.2013.12.007
[9]	S. Mirjalili, A. Lewis, The whale optimization algorithm, Adv. Eng. Software, 95 (2016), 51–67. https://doi.org/10.1016/j.advengsoft.2016.01.008 doi: 10.1016/j.advengsoft.2016.01.008
[10]	M. Dorigo, M. Birattari, T. Stutzle, Ant colony optimization, IEEE Comput. Intell. Mag., 1 (2006), 28–39. https://doi.org/10.1109/MCI.2006.329691 doi: 10.1109/MCI.2006.329691
[11]	H. Jia, X. Peng, C. Lang, Remora optimization algorithm, Expert Syst. Appl. 185 (2021), 115665. https://doi.org/10.1016/j.eswa.2021.115665
[12]	S. Mirjalili, Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm, Knowledge-Based Syst., 89 (2015), 228–249. https://doi.org/10.1016/j.knosys.2015.07.006 doi: 10.1016/j.knosys.2015.07.006
[13]	A. E. Ezugwu, J. O. Agushaka, L. Abualigah, S. Mirjalili, A. H. Gandomi, Prairie dog optimization algorithm, Neural Comput. Appl., 34 (2022), 20017–20065. https://doi.org/10.1007/s00521-022-07530-9 doi: 10.1007/s00521-022-07530-9
[14]	A. Seyyedabbasi, F. Kiani, Sand Cat swarm optimization: A nature-inspired algorithm to solve global optimization problems, Eng. Comput., 2022 (2022), 1–25. https://doi.org/10.1007/s00366-022-01604-x doi: 10.1007/s00366-022-01604-x
[15]	R. V. Rao, V. J. Savsani, D. P. Vakharia, Teaching-learning-based optimization: a novel method for constrained mechanical design optimization problems, Comput.-Aided Des., 43 (2011), 303–315. https://doi.org/10.1016/j.cad.2010.12.015 doi: 10.1016/j.cad.2010.12.015
[16]	A. Kaveh, V. R. Mahdavi, Colliding bodies optimization: a novel meta-heuristic method, Comput. Struct., 139 (2014), 18–27. https://doi.org/10.1016/j.compstruc.2014.04.005 doi: 10.1016/j.compstruc.2014.04.005
[17]	T. T. Huan, A. J. Kulkarni, J. Kanesan, C. J. Huang, A. Abraham, Ideology algorithm: a socio-inspired optimization methodology, Neural Comput. Appl., 28 (2017), 845–876. https://doi.org/10.1007/s00521-016-2379-4 doi: 10.1007/s00521-016-2379-4
[18]	Y. Shi, Brain storm optimization algorithm, in Advances in Swarm Intelligence: Second International Conference, Springer, (2011), 303–309. https://doi.org/10.1007/978-3-642-21515-5_36
[19]	E. Atashpaz-Gargari, C. Lucas, Imperialist competitive algorithm: an algorithm for optimization inspired by imperialistic competition, in 2007 IEEE Congress on Evolutionary Computation, IEEE, (2007), 4661–4667. https://doi.org/10.1109/CEC.2007.4425083
[20]	Z. W. Geem, J. H. Kim, G. V. Loganathan, A new heuristic optimization algorithm: harmony search, Simulation, 76 (2001), 60–68. https://doi.org/10.1177/003754970107600201 doi: 10.1177/003754970107600201
[21]	S. He, Q. H. Wu, J. R. Saunders, Group search optimizer: an optimization algorithm inspired by animal searching behavior, IEEE Trans. Evol. Comput., 13 (2009), 973–990. https://doi.org/10.1109/TEVC.2009.2011992 doi: 10.1109/TEVC.2009.2011992
[22]	Y. Zhang, Z. Jin, Group teaching optimization algorithm: A novel metaheuristic method for solving global optimization problems, Expert Syst. Appl., 148 (2020), 113246. https://doi.org/10.1016/j.eswa.2020.113246 doi: 10.1016/j.eswa.2020.113246
[23]	S. Mirjalili, SCA: a sine cosine algorithm for solving optimization problems, Knowledge-Based Syst., 96 (2016), 120–133. https://doi.org/10.1016/j.knosys.2015.12.022 doi: 10.1016/j.knosys.2015.12.022
[24]	D. Bertsimas, J. Tsitsiklis, Simulated annealing, Stat. Sci., 8 (1993), 10–15. https://doi.org/10.1214/ss/1177011077 doi: 10.1214/ss/1177011077
[25]	R. A. Formato, Central force optimization, Prog. Electromagn. Res., 77 (2007), 425–491. http://doi.org/10.2528/PIER07082403 doi: 10.2528/PIER07082403
[26]	S. Mirjalili, S. M. Mirjalili, A. Hatamlou, Multi-verse optimizer: a nature-inspired algorithm for global optimization, Neural Comput. Appl., 27 (2016), 495–513. https://doi.org/10.1007/s00521-015-1870-7 doi: 10.1007/s00521-015-1870-7
[27]	L. Abualigah, A. Diabat, S. Mirjalili, M. Abd Elaziz, A. H. Gandomi, The arithmetic optimization algorithm, Comput. Methods Appl. Mech. Eng., 376 (2021), 113609. https://doi.org/10.1016/j.cma.2020.113609 doi: 10.1016/j.cma.2020.113609
[28]	A. Hatamlou, Black hole: A new heuristic optimization approach for data clustering, Inf. Sci., 222 (2013), 175–184. https://doi.org/10.1016/j.ins.2012.08.023 doi: 10.1016/j.ins.2012.08.023
[29]	E. Rashedi, H. Nezamabadi-Pour, S. Saryazdi, GSA: a gravitational search algorithm, Inf. Sci., 179 (2009), 2232–2248. https://doi.org/10.1016/j.ins.2009.03.004 doi: 10.1016/j.ins.2009.03.004
[30]	H. Du, X. Wu, J. Zhuang, Small-world optimization algorithm for function optimization//advances in natural computation: Second international conference, in Advances in Natural Computation: Second International Conference, ICNC 2006, Springer, (2006), 264–273. https://doi.org/10.1007/11881223_33
[31]	W, Banzhaf, J. R. Koza, C. Ryan, L. Spector, C. Jacob, Genetic programming, IEEE Intell. Syst. Appl., 15 (2000), 74–84. https://doi.org/10.1109/5254.846288 doi: 10.1109/5254.846288
[32]	K. V. Price, Differential evolution, Handb. Optim.: Classical Mod. Approach, 2013 (2013), 187–214. https://doi.org/10.1007/978-3-642-30504-7_8 doi: 10.1007/978-3-642-30504-7_8
[33]	X. Yao, Y. Liu, G. Lin, Evolutionary programming made faster, IEEE Trans. Evol. Comput., 3 (1999), 82–102. https://doi.org/10.1109/4235.771163 doi: 10.1109/4235.771163
[34]	D. Simon, Biogeography-based optimization, IEEE Trans. Evol. Comput., 12 (2008), 702–713. https://doi.org/10.1109/TEVC.2008.919004 doi: 10.1109/TEVC.2008.919004
[35]	J. H. Holland, Genetic algorithms, Sci. Am., 267 (1992), 66–73. https://doi.org/10.1038/scientificamerican0792-66 doi: 10.1038/scientificamerican0792-66
[36]	H. G. FBeyer, H. P. Schwefel, Evolution strategies-a comprehensive introduction, Nat. Comput., 1 (2002), 3–52. https://doi.org/10.1023/A:1015059928466 doi: 10.1023/A:1015059928466
[37]	M. Jaderyan, H. Khotanlou, Virulence optimization algorithm, Appl. Soft. Comput., 43 (2016), 596–618. https://doi.org/10.1016/j.asoc.2016.02.038 doi: 10.1016/j.asoc.2016.02.038
[38]	D. H. Wolpert, W. G. Macready, No free lunch theorems for optimization, IEEE Trans. Evol. Comput., 1 (1997), 67–82. https://doi.org/10.1109/4235.585893 doi: 10.1109/4235.585893
[39]	S. Wang, A. G. Hussien, H. Jia, L. Aualigah, R. Zheng, Enhanced remora optimization algorithm for solving constrained engineering optimization problems, Mathematics, 10 (2022), 1696. https://doi.org/10.3390/math10101696 doi: 10.3390/math10101696
[40]	Z. Cui, X. Hou, H. Zhou, W. Lian, J. Wu, Modified slime mould algorithm via levy flight, in 2020 13th International Congress on Image and Signal Processing, Biomedical Engineering and Informatics (CISP-BMEI), IEEE, (2020), 1109–1113. https://doi.org/10.1109/CISP-BMEI51763.2020.9263669
[41]	D. Wu, H. Rao, C. Wen, H. Jia, Q. Liu, L. Abualigah, Modified sand cat swarm optimization algorithm for solving constrained engineering optimization problems, Mathematics, 10 (2022), 4350. https://doi.org/10.3390/math10224350 doi: 10.3390/math10224350
[42]	M. H. Nadimi-Shahraki, H, Mohammad, H. D. Zamani, S. Mirjalili, M. A. Elaziz, MMKE: Multi-trial vector-based monkey king evolution algorithm and its applications for engineering optimization problems, Plos One, 18 (2023), e0280006. https://doi.org/10.1371/journal.pone.0280006 doi: 10.1371/journal.pone.0280006
[43]	M. H. Nadimi-Shahraki, H, Mohammad, A. A. Ewees, L. Abualigah, Mtv-mfo: Multi-trial vector-based moth-flame optimization algorithm, Symmetry, 13 (2021), 2388. https://doi.org/10.3390/sym13122388 doi: 10.3390/sym13122388
[44]	M. H. Nadimi-Shahraki, E. Moeini, S. Taghian, S. Mirjalili, DMFO-CD: a discrete moth-flame optimization algorithm for community detection, Algorithms, 14 (2021), 314. https://doi.org/10.3390/a14110314 doi: 10.3390/a14110314
[45]	Y. Yang, Y. Gao, S. Tan, S. Zhao, J. Wu, S. Gao, et al., An opposition learning and spiral modelling based arithmetic optimization algorithm for global continuous optimization problems, Eng. Appl. Artif. Intell., 113 (2022), 104981. https://doi.org/10.1016/j.engappai.2022.104981 doi: 10.1016/j.engappai.2022.104981
[46]	Y. Yang, C. Qian, H. Li, Y. Gao, J. Wu, C. J. Liu, et al., An efficient DBSCAN optimized by arithmetic optimization algorithm with opposition-based learning, J. Supercomput., 78 (2022), 19566–19604. https://doi.org/10.1007/s11227-022-04634-w doi: 10.1007/s11227-022-04634-w
[47]	M. H. Nadimi-Shahraki, H. Mohammad, S. Mirjalili, L. Abualigah, Binary aquila optimizer for selecting effective features from medical data: a covid-19 case study, Mathematics, 10 (2022), 1929. https://doi.org/10.3390/math10111929 doi: 10.3390/math10111929
[48]	S. Sharma, A. K. Saha, G. Lohar, Optimization of weight and cost of cantilever retaining wall by a hybrid metaheuristic algorithm, Eng. Comput., 2021 (2021), 1–27. https://doi.org/10.1007/s00366-021-01294-x doi: 10.1007/s00366-021-01294-x
[49]	M. Masdari, S. Barshandeh, Discrete teaching-learning-based optimization algorithm for clustering in wireless sensor networks, Intell. Humaniz. Comput., 11 (2020), 5459–5476. https://doi.org/10.1007/s12652-020-01902-6 doi: 10.1007/s12652-020-01902-6
[50]	H. R. Tizhoosh, Opposition-based learning: a new scheme for machine intelligence, in International Conference on Computational Intelligence for Modelling, Control and Automation and International Conference on Intelligent Agents, Web Technologies and Internet Commerce (CIMCA-IAWTIC'06), IEEE, (2005), 695–701. https://doi.org/10.1109/CIMCA.2005.1631345
[51]	M. Li, G. Xu, B. Fu, X. Zhao, Whale optimization algorithm based on dynamic pinhole imaging and adaptive strategy, J. Supercomput., 2022 (2022), 1–31. https://doi.org/10.1007/s11227-021-04116-5 doi: 10.1007/s11227-021-04116-5
[52]	N. A. Dodgson, Quadratic interpolation for image resampling, IEEE Trans. Image Process., 6 (1997), 1322–1326. https://doi.org/10.1109/83.623195 doi: 10.1109/83.623195
[53]	A. A. Heidari, S. Mirjalili, H. Faris, I. Aljarah, M. Mafarja, H. Chen, Harris hawks optimization: Algorithm and applications, Future Gener. Comput. Syst., 97 (2019), 849–872. https://doi.org/10.1016/j.future.2019.02.028 doi: 10.1016/j.future.2019.02.028
[54]	H. A. Alsattar, A. A. Zaidan, B. B. Zaidan, Novel meta-heuristic bald eagle search optimisation algorithm, Artif. Intell. Rev., 53 (2020), 2237–2264. https://doi.org/10.1007/s10462-019-09732-5 doi: 10.1007/s10462-019-09732-5
[55]	G. Wu, M. Rammohan, P. N. Suganthan, Problem definitions and evaluation criteria for the CEC 2017 competition on constrained real-parameter optimization, Natl. Univ. Def. Technol., 2017 (2017).
[56]	C. Wen, H. Jia, D. Wu, H. Rao, S. Li, Q. Liu, et al., Modified remora optimization algorithm with multistrategies for global optimization problem, Mathematics, 10 (2022), 3604. https://doi.org/10.3390/math10193604 doi: 10.3390/math10193604
[57]	E. Cuevas, P. Diaz, O. Camarena, E. Cuevas, P. Diaz, O. Camarena, Experimental analysis between exploration and exploitation, in Metaheuristic Computation: A Performance Perspective, Springer, (2021), 249–269. https://doi.org/10.1007/978-3-030-58100-8_10
[58]	J. O. Agushaka, A. E. Ezugwu, L. Abualigah, Dwarf mongoose optimization algorithm, Comput. Methods Appl. Mech. Eng., 391 (2022), 114570. https://doi.org/10.1016/j.cma.2022.114570 doi: 10.1016/j.cma.2022.114570
[59]	E. H. Houssein, N. Neggaz, M. E. Hosney, W. M. Mohamed, M. Hassaballah, Enhanced Harris hawks optimization with genetic operators for selection chemical descriptors and compounds activities, Neural Comput. Appl., 33 (2021), 13601–13618. https://doi.org/10.1007/s00521-021-05991-y doi: 10.1007/s00521-021-05991-y
[60]	W. Long, J. Jiao, X. Liang, S. Cai, A random opposition-based learning grey wolf optimizer, IEEE Access, 7 (2019), 113810–113825. https://doi.org/10.1109/ACCESS.2019.2934994 doi: 10.1109/ACCESS.2019.2934994
[61]	E. H. Houssein, N. Neggaz, M. E. Hosney, M. E. Hosney, W. M. Mohamed, M. Hassaballah, Enhanced Harris hawks optimization with genetic operators for selection chemical descriptors and compounds activities, Neural Comput. Appl., 33 (2021), 13601–13618. https://doi.org/10.1007/s00521-021-05991-y doi: 10.1007/s00521-021-05991-y
[62]	A. G. Hussien, An enhanced opposition-based salp swarm algorithm for global optimization and engineering problems, J. Ambient. Intell. Humaniz. Comput., 13 (2022), 129–150. https://doi.org/10.1007/s12652-021-02892-9 doi: 10.1007/s12652-021-02892-9
[63]	G. Sayed, A. Darwish, A. E. Hassanien, A new chaotic multi-verse optimization algorithm for solving engineering optimization problems, J. Exp. Theor. Artif. Intell., 30 (2018), 293–317. https://doi.org/10.1080/0952813X.2018.1430858 doi: 10.1080/0952813X.2018.1430858
[64]	L. Abualigah, M. A. Elaziz, P. Sumari, Z. W. Geem, A. H. Gandomi, Reptile Search Algorithm (RSA): A nature-inspired meta-heuristic optimizer, Expert Syst. Appl., 191 (2021), 116158. https://doi.org/10.1016/j.eswa.2021.116158 doi: 10.1016/j.eswa.2021.116158
[65]	H. Eskandar, A. Sadollah, A. Bahreininejad, M. Hamdi, Water cycle algorithm-a novel metaheuristic optimization method for solving constrained engineering optimization problems, Comput. Struct., 110 (2012), 151–166. https://doi.org/10.1016/j.compstruc.2012.07.010 doi: 10.1016/j.compstruc.2012.07.010
[66]	A. H. Gandomi, X. S. Yang, A. H. Alavi, Cuckoo search algorithm: A metaheuristic approach to solve structural optimization problems, Eng. Comput., 29 (2013), 17–35. https://doi.org/10.1007/s00366-011-0241-y doi: 10.1007/s00366-011-0241-y
[67]	A. Baykasoglu, S. Akpinar, Weighted superposition attraction (WSA): A swarm intelligence algorithm for optimization problems-part2: Constrained optimization, Appl. Soft Comput., 37 (2015), 396–415. https://doi.org/10.1016/j.asoc.2015.08.052 doi: 10.1016/j.asoc.2015.08.052
[68]	J. M. Czerniak, H. Zarzycki, D. Ewald, Aao as a new strategy in modeling and simulation of constructional problems optimization, Simul. Modell. Pract. Theory, 76 (2017), 22–33. https://doi.org/10.1016/j.simpat.2017.04.001 doi: 10.1016/j.simpat.2017.04.001
[69]	A. Baykasoglu, F. B. Ozsoydan, Adaptive firefly algorithm with chaos for mechanical design optimization problems, Appl. Soft Comput., 36 (2015), 152–164. https://doi.org/10.1016/j.asoc.2015.06.056 doi: 10.1016/j.asoc.2015.06.056

This article has been cited by:

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

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)