Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters

Chi Zhai

^1,2,

Chunxi Yang

and

Jing Na

^1,*

Faculty of Mechanical and Electrical Engineering, Kunming University of Science and Technology, Kunming 650093, China

Faculty of Chemical Engineering, Kunming University of Science and Technology, Kunming 650500, China

Author to whom correspondence should be addressed.

Sensors2022,22(23), 9334;https://doi.org/10.3390/s22239334

Submission received: 12 November 2022 /Revised: 25 November 2022 /Accepted: 27 November 2022 /Published: 30 November 2022

(This article belongs to the CollectionMulti-Sensor Information Fusion)

Downloadkeyboard_arrow_down

Browse Figures

Review Reports Versions Notes

Abstract

Information fusion integrates aspects of data and knowledge mostly on the basis that system information is accumulative/distributive, but a subtle case emerges for a system with bifurcations, which is always un-linearizable and exacerbates information acquisition and presents a control problem. In this paper, the problem of an un-linearizable system related to system observation and control is addressed, and Andronov–Hopf bifurcation is taken as a typical example of an un-linearizable system and detailed. Firstly, the properties of a linear/linearized system is upon commented. Then, nonlinear degeneracy for the normal form of Andronov–Hopf bifurcation is analyzed, and it is deduced that the cubic terms are an integral part of the system. Afterwards, the theoretical study on feedback stabilization is conducted between the normal-form Andronov–Hopf bifurcation and its linearized counterpart, where stabilization using washout-filter-aided feedback is compared, and it is found that by synergistic controller design, the dual-conjugate-unstable eigenvalues can be stabilized by single stable washout filter. Finally, the high-dimensional ethanol fermentation model is taken as a case study to verify the proposed bifurcation control method.

Keywords:

Andronov–Hopf bifurcation;washout filter;normal-form degeneracy;feedback stabilization

1. Introduction

In nonlinear systems covering ecology, economics, sociology and many other fields, complex dynamic evolutionary behaviors, such as multiplicity, self-oscillation, and even chaos arise in the vicinity of bifurcation points [1], and various methods are adopted to enhance sensory ability for accurate manipulation and control: the bifurcation oscillator is applied as a quantum state amplifier [2], the oscillating clearance of the rolling-element bear is taken as the bifurcating parameter with repetitive errors fused into the control system [3], and deep learning is implemented for the control of the Van del Pol system [4]. Mathematically, a bifurcation takes place at a parameter value where the system loses structural stability with respect to parameter variations, i.e., a phase portrait around the equilibrium is not locally topologically conjugate to the phase portraits around the equilibrium at nearby parameter values [5,6], which leads the solution manifold to change significantly, and exacerbates system information acquisition and control implementation [7,8].

Bifurcation control [9] refers to the task of designing a controller to suppress or reduce existing bifurcating dynamics, thereby achieving desirable dynamical behaviors, i.e., to obtain an asymptotically stable phase portrait. Feedback is often applied on a system with bifurcations, where the methods of linear/nonlinear feedback, washout-filter-aided feedback, impulse control, vibrational control, delayed feedback [10], etc., are frequently adopted. Nonlinear analysis and control based on harmonic balance approximations or nonlinear frequency response is also studied [11]. Recently, the idea of the anti-control of bifurcations [12] is proposed to segregate unwanted dynamics by deliberately introducing simple bifurcations into the system.

However, bifurcation control sometimes exhibits contradictory outcomes: improper control actions or process–model mismatch might cause the closed-loop system to move across the stable manifold, leading to more complicated dynamic behaviors [13], and one cause is that many control laws are developed on the basis of linear system theory (gain-scheduling, feedback linearization, etc.), but most dynamic systems with bifurcation are un-linearizable, i.e., truncating out the nonlinear terms alters the structural stability of the original system, and the higher-ordered-terms might not shrink to zero as the time approaches infinity. On the other hand, normal form degeneracy provides an effective way to evaluate those nonlinear terms, where types of bifurcations are expressed with quadratic or cubic degeneracy, and state-feedback based on normal forms and invariant sets are developed [14,15].

The main goal of this paper was to address the problem of an un-linearizable system in bifurcation control, where a comparative study is conducted on both the original and linearized models using washout-filter-aided state feedback. Through washout filters, a feedback signal is expected to be well geared with the open-loop dynamics and generate asymptotic outcomes. Abed and co-workers [16,17,18,19] implemented washout filters for bifurcation control, and feedback stabilization theory is proposed through auxiliary state variables that are inter-coupled. For systems with varying bifurcation parameters, the adaptive control scheme is also developed [12]. Since washout-aided feedback is introduced to migrate out of unstable eigenvalues, Marquardt [20] suggested manipulating directly on the unstable eigenspace, and introducing a pair of complex conjugate eigenvalues instead. Zhang [21] further employed eigenvalue assignment to achieve bifurcation control with guaranteed transient dynamics. Chen [22] proposed a systematic way of choosing constants and control parameters for the control of a high-dimensional bifurcating system. However, for a complicated system, i.e., both single- and dual-unstable eigenvalues exist, and controller design could be rather sophisticated.

The aforementioned washout-filter-aided control was developed on the linearized system, and the linearization procedure might change the structural stability of the original system and cause redundant control actions. In this work, the elementary properties of a linearizable system are offered, which provides a comparison with the un-linearizable system when these properties fail. Then, the system under Andronov–Hopf bifurcation is taken as a typical case of an un-linearizable one and the related analysis and bifurcation control strategy are developed. Since washout-filter-aided control can migrate the unstable eigenvalues, numerical bifurcation analysis on the closed-loop system reveals that, by inserting only one washout filter to any of the dual-unstable eigenspaces, the system might be stabilized. Compared with previous developed control strategies [23], our method could downsize the numbers of washout filters, facilitating the application for more practical and general cases, especially for high-dimensional systems.

The inconsistency of feedback stabilization and numerical bifurcation analysis stems from the fact that many nonlinear systems with bifurcation are un-linearizable [24]. This work is organized as follows: Through the analysis of the normal form of Andronov–Hopf bifurcation, the properties of an un-linearizable system are detailed in this paper; then, washout-filter-aided feedback on the un-linearizable system is analyzed in contrast with previously developed washout-filter-aided control strategies; and the continuous ethanol fermentation model is taken as an example of the un-linearizable system and is discussed in detail, and finally, concluding remarks are provided last.

2. The Un-Linearizable System

2.1. Properties of the Linear System

Previously to the discussion on the nonlinear dynamic system, basic properties concerning thelinear system are provided. Consider a system exhibiting outcomes (characters)M, and which is detected with finite numbers of sensible causes

x_{1}, \dots, x_{n}

, related by functionf; then, a linear/linearizable system starts as a weighted contribution of individual variables

M = f (x_{1}, x_{2}, \dots, x_{n}) \approx \frac{\partial f}{\partial x_{1}} x_{1} + \frac{\partial f}{\partial x_{2}} x_{2} \dots + \frac{\partial f}{\partial x_{n}} x_{n} .

(1)

Furthermore, the coupling relation between causes and outcomes is decoupled through variable substitution and combination, e.g., a 2 × 2 input–output system with a linearized form is given as follows,

M = [\begin{matrix} M_{1} \\ M_{2} \end{matrix}] = [\begin{matrix} f_{1} (x_{1}, x_{2}) \\ f_{1} (x_{1}, x_{2}) \end{matrix}] \approx [\begin{matrix} \frac{\partial f_{1}}{\partial x_{1}} & \frac{\partial f_{1}}{\partial x_{2}} \\ \frac{\partial f_{2}}{\partial x_{1}} & \frac{\partial f_{2}}{\partial x_{2}} \end{matrix}] [\begin{matrix} x_{1} \\ x_{2} \end{matrix}] .

(2)

Assume that the Jacobian matrix is non-singular, and a coordinate change

y = A x

exists that could diagonalize Equation (2)

\begin{matrix} [\begin{matrix} M_{1} \\ M_{2} \end{matrix}] \approx [\begin{matrix} \frac{\partial f_{1}}{\partial x_{1}} & \frac{\partial f_{1}}{\partial x_{2}} \\ \frac{\partial f_{2}}{\partial x_{1}} & \frac{\partial f_{2}}{\partial x_{2}} \end{matrix}] [\begin{matrix} x_{1} \\ x_{2} \end{matrix}] = A^{- 1} [\begin{matrix} J_{11} & 0 \\ 0 & J_{22} \end{matrix}] A [\begin{matrix} x_{1} \\ x_{2} \end{matrix}] = A^{- 1} [\begin{matrix} J_{11} y_{1} \\ J_{22} y_{2} \end{matrix}], \\ w h e r e : \\ [\begin{matrix} \frac{\partial f_{1}}{\partial x_{1}} & \frac{\partial f_{1}}{\partial x_{2}} \\ \frac{\partial f_{2}}{\partial x_{1}} & \frac{\partial f_{2}}{\partial x_{2}} \end{matrix}] = A^{- 1} [\begin{matrix} J_{11} & 0 \\ 0 & J_{22} \end{matrix}] A, y = [\begin{matrix} y_{1} \\ y_{2} \end{matrix}] = A [\begin{matrix} x_{1} \\ x_{2} \end{matrix}] . \end{matrix}

(3)

Equation (3) indicates that for a linear/linearizable system, outcomes can be departed as a sum of the weighted causes separately, and sets of causes can be found to correlate every outcome exclusively.

The above

l i n e a r s y s t e m

properties are held for the linear-time-invariant (LTI) dynamic system. Given in state-space form, a multi-input–multi-output (MIMO) system undergoes Laplace transform and obtains similar distributive and decoupling properties

\{\begin{matrix} x^{'} = A x + B u + ξ (x, u) \\ y = C x + D u \end{matrix} \overset{\leftrightarrow}{L a p l a c e} \frac{y}{u} = G (s) \approx C \frac{B}{s I - A} + D,

(4)

where

ξ

is the nonlinear terms that are truncated out in the linearizable dynamic system, whilsty andu represent outcomes and causes, respectively.

Mention that many nonlinear controls are developed on the linearized system. In practical applications, the presence of system nonlinearities causes the dynamic behavior to be qualitatively different from one operating regime to another; hence, one might obtain linear (at a specific set-point) models valid within a “small” region about the linearization point, and perform local designs for a set of operating conditions and then construct a gain-scheduling scheme that interpolates controller gains as the process traverses the operating region.

Accordingly, un-linearizable dynamic systems are those truncating out

ξ

as it alters the structural stability of the original system, and many dynamic systems with bifurcation are un-linearizable. Since a system under Andronov–Hopf bifurcation has dual- unstable eigenvalues with planar phase portrait, which is a typical multivariate, un-linearizable case, the following analysis and control is focused on the systems with Andronov–Hopf bifurcation.

2.2. Normal Form Andronov–Hopf Bifurcation: A Typical Un-Linearizable System

A standard approach to study the behavior of ordinary differential equations (ODE) around a bifurcation point is through normal form analysis, which starts by center manifold degeneracy. For a parameterized dynamical system,

x^{'} = f (x, a),

(5)

where

x \to f (x, α) \in R^{n}

and

α

is the bifurcation parameter, with equilibrium

x_{0}

satisfying

f (x_{0}, 0) = 0

. Assume eigenvalues of the Jacobian about

x_{0}

are

λ_{1, 2} = \pm i ω, λ_{k} < 0 (k = 3, 4, \dots, n)

, then, there exists n-2 invariant/center manifold that is tangent at the bifurcation point to the eigenspace of the neutrally stable eigenvalues

λ_{1, 2}

. Specifically, for a small neighborhood of

α \in (0, δ)

and within the invariant manifold exists a 2-dimensional attracting manifold and the coordinate change of the first two eigenspace; by introducing complex conjugatesz and

z^{*}

, Equation (5) around

x_{0}

can be given as

z^{'} = i ω z + g (z, z^{*}) z \in C^{1},

(6)

where

g (z, z^{*})

is the degenerated terms at

α = 0

and can be given explicitly through rules of implicit function theorem. Assume

S (z, z^{*}) \in R^{(n - 2)}

represents the stable manifold, the n-2 stable equivalence is provided as follows

ω (z^{*} \frac{\partial S}{\partial z} (z, z^{*}) - z \frac{\partial S}{\partial z^{*}} (z, z^{*})) = f_{3, \dots, n} (S (z, z^{*}), p (z, z^{*})),

(7)

and the n-2 solutions are exponentially convergent to

S (t)

x_{3, \dots, n} (t) = S (t) + O (e^{- β t}), a s t \to \infty,

(8)

where

β

is the positive convergent rate. Considering that

x_{3, . . ., n}

are locally asymptotic and therefore can be neglected in a local stability analysis around the bifurcation point, in the following discussion, about the Andronov–Hopf bifurcation, a two-dimensional topologically equivalent eigenspace containing dual-unstable (

λ_{1, 2} = \pm i ω

) is discussed instead of Equation (5).

For a general two-dimensional parameterized system, i.e.,

x = (x_{1}, x_{2}) \to f (x, α) \in R^{2}, α \in R

, linearization about the equilibrium gives

\begin{matrix} x^{'} = A (α) x + h (x, α), \\ w h e r e : \\ A (α) = (\begin{matrix} a & b \\ c & d \end{matrix}), \end{matrix}

(9)

where

α

,b,c andd are elements of the Jacobian

A (α)

, which is simplified as

A; h (x, α)

represents the higher-ordered term. The characteristic equation corresponding toA gives

\begin{matrix} λ^{2} - t r (A) λ + d e t (A) = 0 \\ \Rightarrow \{\begin{matrix} λ = \frac{1}{2} (t r (A) + \sqrt{t r {(A)}^{2} - 4 d e t (A)}) \\ λ^{*} = \frac{1}{2} (t r (A) - \sqrt{t r {(A)}^{2} - 4 d e t (A)}) \end{matrix}, \end{matrix}

(10)

where

t r (A)

and

d e t (A)

are the trace and determinant, respectively. Andronov–Hopf bifurcation takes place at

α = 0

satisfying

t r (A) = 0

and

d e t (A) = ω_{0}^{2}

. Assume in a small neighborhood

α \in (0, δ)

that the eigenvalues are written as

λ_{1} (α) = μ + i ω

and

λ_{2} (α) = μ - i ω

, where

μ = t r (A)

and

ω^{2} = t r {(A)}^{2} - d e t (A)

; then,

T (α)

exists to transform

A (α)

into canonical real form as follows

J = T (α) A (α) T^{- 1} (α) = (\begin{matrix} μ (α) & - ω (α) \\ ω (α) & μ (α) \end{matrix}),

(11)

where

J^{- 1} = J^{T} / (μ^{2} + ω^{2})

. Let

y = T (α) x

and Equation (9) turns to

\{\begin{matrix} y_{1}^{'} = μ (α) y_{1} - ω (α) y_{2} + {o (| y |}^{2}) \\ y_{2}^{'} = ω (α) y_{1} + μ (α) y_{2} + {o (| y |}^{2}) \end{matrix} .

(12)

Introducing

z = y_{1} + i y_{2}

and

z^{*} = y_{1} - i y_{2}

and Equation (12) presents one-dimensional equivalence as follows

z^{'} = λ (α) z + g (z, z^{*}),

(13)

with

g = o (z^{2})

.Let

q (α)

be the eigenvector of

λ_{1}

about

A (α)

and one can give an eigenvector

p (α)

corresponding to

λ_{2}

. Since

λ_{1}

and

λ_{2}

are complex conjugates,p andq can be chosen so that the inner product satisfies

〈 q, p 〉 = 1

, and

\begin{matrix} 〈 p, q^{*} 〉 = 〈 p, \frac{1}{λ_{2}} A q^{*} 〉 = \frac{1}{λ_{2}} 〈 A^{T} p, q^{*} 〉 = \frac{λ_{1}}{λ_{2}} 〈 p, q^{*} 〉 \\ \Rightarrow 〈 p, q^{*} 〉 = 0 . \end{matrix}

(14)

Equation (11) after transposition has

T^{T} (α) J^{T} = A^{T} (α) T^{T} (α)

, and if

p (α)

has the following form

p (α) = (\begin{matrix} p_{1} \\ p_{2} \end{matrix}) = (\begin{matrix} T_{11} - i T_{21} \\ T_{12} - i T_{22} \end{matrix}),

(15)

where

T_{i j}

isi-th line,j-th row element, thenz has the explicit form as follows

z = 〈 p (α), x 〉 .

(16)

Considering

〈 q, p 〉 = 1

and Equation(13), the original variablex has the only equivalent transform

x = (\begin{matrix} q \\ q^{*} \end{matrix}) (z z^{*}) .

(17)

Substituting Equation (17) into (16), one obtains

z^{'} = λ (α) z + 〈 p (α), F (z q + z^{*} q^{*}, α) 〉 .

(18)

Here, the nonlinear term

F (x, α)

in Equation (9) after Taylor expansion gives

\begin{matrix} F (x, 0) = \frac{1}{2} B (x, x) + \frac{1}{6} C (x, x, x) + {O (| | x | |}^{4}), \\ w i t h \{\begin{matrix} B_{1, 2} (x, y) = \sum_{j, k = 1}^{2} \frac{\partial^{2} F_{1, 2} (ξ, 0)}{\partial ξ_{j} \partial ξ_{k}} |_{ξ = 0} x_{j} y_{k} \\ C_{1, 2} (x, y, m) = \sum_{j, k = 1, l = 1}^{2} \frac{\partial^{3} F_{1, 2} (ξ, 0)}{\partial ξ_{j} \partial ξ_{k} \partial ξ_{m}} |_{ξ = 0} x_{j} y_{k} m_{l} \end{matrix} . \end{matrix}

(19)

Substituting Equation (19) into (18) and the second-order term B gives

z^{'} = λ z + \frac{〈 p, B (q, q) 〉}{2} z^{2} + 〈 p, B (q, q^{*}) 〉 z z^{*} + \frac{〈 p, B (q^{*}, q^{*}) 〉}{2} z^{* 2} + {O (| z |}^{3}) .

(20)

Through a coordinate change in Equation (21) the second-order term vanishes and gives Equation (22)

\begin{matrix} z = ω + \frac{〈 p, B (q, q) 〉}{2 λ} ω^{2} + & \frac{〈 p, B (q, q^{*}) 〉}{λ^{*}} ω ω^{*} + \frac{〈 p, B (q^{*}, q^{*}) 〉}{2 (2 λ^{*} - λ)} ω^{* 2} \end{matrix}

(21)

\begin{matrix} ω^{'} & = λ ω + O (| ω |^{3}) \end{matrix}

(22)

Similarly, the third-order term after coordinate change obtains

\begin{matrix} ω^{'} = λ ω + σ ω^{2} ω^{*} + {O (| ω |}^{4}), \\ w h e r e : \\ σ = \frac{〈 p, C (q, q, q^{*}) 〉}{2} . \end{matrix}

(23)

Next, let

ω = ρ e^{i ϕ}

and Equation (23) turns to

\begin{matrix} \{\begin{matrix} ρ^{'} = ρ (α - ρ^{2}) + Φ (ρ, φ) \\ φ^{'} = 1 + Ψ (ρ, φ) \end{matrix}, \\ w h e r e : \\ Φ = {O (| ρ |}^{4} {); Ψ = O (| ρ |}^{3}) . \end{matrix}

(24)

Suppose that

φ

performs constant rate rotation and Equation (24) mapped to the rotating plane gives

\begin{matrix} \frac{d ρ}{d φ} = \frac{ρ (α - ρ^{2}) + Φ}{1 + Ψ} = ρ (α - ρ^{2}) + (ρ, φ), \\ w h e r e : \\ R = O (| ρ |^{4}) . \end{matrix}

(25)

Assume that

ρ_{0} = 0

is the initial condition, solving Equation (24) that provides

ρ = e^{α φ} ρ_{0} + e^{α φ} \frac{1 - e^{2 α φ}}{2 α} ρ_{0}^{3} + O_{1} (| ρ_{0} |^{4}) .

(26)

For a fixed

α > 0

, there exists a limit cycle, and when

φ = 2 π

, the state

ρ_{1}

would turn back to the initial point

ρ_{1} = e^{2 π α} ρ_{0} - e^{2 π α} (2 π + O (α)) ρ_{0}^{3} + O_{2} (| ρ_{0} |^{4}) .

(27)

Comparing Equations (26) and (27), it is found that the 4th- or higher-order terms are unaffected by the formation of limit cycles, and hence can be truncated out. Furthermore, an Andronov–Hopf bifurcation is structurally equivalent to the following standard form

(\binom{{\hat{x}}_{1}^{'}}{{\hat{x}}_{2}^{'}}) = (\begin{matrix} α & - 1 \\ 1 & α \end{matrix}) (\begin{matrix} {\hat{x}}_{1} \\ {\hat{x}}_{2} \end{matrix}) \pm ({\hat{x}}_{1}^{2} + {\hat{x}}_{2}^{2}) (\begin{matrix} {\hat{x}}_{1} \\ {\hat{x}}_{2} \end{matrix}) .

(28)

As shown in the following schematic diagram, if Equation (28) is linearized by truncating out the nonlinear terms, the phase portrait shows that the outcome is asymptotically stable to the origin (Figure 1a); if the nonlinear terms at Equation (28) are preserved, the phase portraits of the Andronov–Hopf bifurcation would be super-critical (Figure 1b), where the periodic solution is stable, or sub-critical (Figure 1c), where the periodic solution is unstable, and the criterion for super- or sub-criticality is judged by the sign of the third-order term of Equation (28).

It is an obvious standard form of the Andronov–Hopf bifurcation system that differs from the linearized system as it either attracts or repels an oscillatory trajectory, while the linearized one exhibits a diverging trend in the whole phase portrait. Hence, the nonlinear terms in Equation (28) are indispensable, and a system with Andronov–Hopf bifurcation is un-linearizable.

Address that there exists two variables

{\hat{x}}_{1}

and

{\hat{x}}_{2}

interacting with one another to form a limit cycle around the Andronov–Hopf bifurcation point (

α > 0

). Furthermore, Equation (28) can be represented one-dimensionally using complex variables, which means manipulating either

{\hat{x}}_{1}

{\hat{x}}_{2}

to affect the dynamic behavior of the system directly (or equally). Hence, introducing one state feedback suffices stabilizing a system with Andronov–Hopf bifurcation. Notion, the sign of the third-order term in Equation (28), decides whether the bifurcation is super-critical (−) or sub-critical (+), and one can refer to [21] for more details.

3. Stabilizing Normal Form Andronov–Hopf Bifurcation

3.1. Washout Filter as a State Manipulator

Since washout filters reject steady-state signals and passing transient ones, it could be used as a state manipulator and by fine-tuned feedback, eliminates the instabilities originating from open-loop bifurcations.Figure 2 is the washout-filter-aided controller design scheme: by inserting washout filters between states and outputs, transient dynamic behaviors are separated out that can be used to cancel out the perturbations of the system through feedback; while outputs preserve the steady state information, this is important for a control system, meaning washout filters are inoperative until vibrational signals are detected. Assume that

x_{k}

is thek-th element of the state variables, andd is the filter time constant. The control law used in this paper is quite simple, andv is the reference input used to determine the actual operating point andk is a controller parameter.

3.2. Feedback Stabilization Based on the Linearized System

To stabilize the normal-form Andronov–Hopf bifurcation, washout-filter-aided control is constructed. Without loss of generality, the linearized system is extended as follows

x^{'} = A x + B u + h (x, u),

(29)

where the control

u = {[u_{1}, u_{2}]}^{T}

takes

u = v + k y

, andk is the feedback gain;B is an arbitrary input-state transfer matrix and

h (x, u)

represents the higher-ordered-terms.

Adding the auxiliary state satisfies

d z / d t = y

, and the closed-loop system is provided as follows

\begin{matrix} \{\begin{matrix} x^{'} = A x + B u + h (x, u) \\ z^{'} = P (x - z) \\ u = K (x - z) \end{matrix} . \end{matrix}

(30)

Stabilizing Equation (30) requires the following augmented matrix

A_{C}

being Hurwitz

\begin{matrix} (\begin{matrix} x^{'} \\ z^{'} \end{matrix}) & = (\begin{matrix} A + B K & - B K \\ P & - P \end{matrix}) (\begin{matrix} x \\ z \end{matrix}) \\ = A_{C} (\begin{matrix} x \\ z \end{matrix}) . \end{matrix}

(31)

Theorem 1.

If A is non-singular and

(A, B)

is stabilizable, there exists

(P, k)

, which makes

A_{C}

Hurwitz.

Proof.

A_{C}

after similarity transforms of

T_{1}

and

T_{2}

gives

\begin{matrix} T_{1} A_{C} T_{1}^{- 1} & = (\begin{matrix} I & 0 \\ 0 & P^{- 1} \end{matrix}) (\begin{matrix} A + B K & - B K \\ P & - P \end{matrix}) {(\begin{matrix} I & 0 \\ 0 & P^{- 1} \end{matrix})}^{- 1} \\ = (\begin{matrix} A + B K & - K B P \\ I & - P \end{matrix}), \\ T_{2} T_{1} A_{C} T_{1}^{- 1} T_{2}^{- 1} & = (\begin{matrix} I & M \\ 0 & I \end{matrix}) (\begin{matrix} A + B K & - B K \\ P & - P \end{matrix}) {(\begin{matrix} I & M \\ 0 & I \end{matrix})}^{- 1} \\ = (\begin{matrix} A + B K + M & - A M - K B M - M^{2} - K B P - M P \\ I & - M - P \end{matrix}) . \end{matrix}

(32)

Through the following small gain analysis, the control (

P, K

) is obtained

\begin{matrix} A M + K B M + M^{2} + K B P + M P = 0, \\ s . t . \{\begin{matrix} P = ε P_{1} + O (ε^{2}) \\ M = M_{0} + ε M_{1} + O (ε^{2}) \end{matrix} . \end{matrix}

(33)

Canceling the first- and second-order perturbational terms provides

\begin{matrix} \{\begin{matrix} (A + B K + M_{0}) M_{0} = 0 \\ (A + B K) M_{1} + A P_{1} = 0 \end{matrix}, \\ w h e r e : \\ M_{0} = - A - B K, \end{matrix}

(34)

which dictates,

M_{1} = - A P_{1} {(A + B K)}^{- 1} .

(35)

Then, around

M_{0}

in a small region, the onlyM satisfying Equation (33) is given

\begin{matrix} M & = M_{0} + ε M_{1} + O (ε^{2}) \\ = - A - B K - ε A P_{1} {(A + B K)}^{- 1} + O (ε^{2}) \end{matrix}

(36)

with

A_{C}

satisfying

\begin{matrix} T_{2} T_{1} A_{C} T_{1}^{- 1} T_{2}^{- 1} = (\begin{matrix} - ε A P_{1} {(A + B K)}^{- 1} + O (ε^{2}) & 0 \\ I & A_{C 2} (2, 2) \end{matrix}), \\ w h e r e : \\ A_{C 2} (2, 2) = A + B K + ε (A P_{1} {(A + B K)}^{- 1} - P_{1}) + O (ε^{2}) . \end{matrix}

(37)

When (

A, B

) is stabilizable,

A + B K

can be designed Hurwitz. SinceA is non-singular, setting

P_{1} = A^{- 1} (A + B K)

, then

- A P_{1} {(A + B K)}^{- 1}

is Hurwitz, so as

A_{C 2}

. Plus, proper feedback gaink can promise the good transient dynamics of the closed-loop system.

One can designP based on above perturbation theory, but how many washout filters needed to stabilize the system is still a problem, the following lemma provides some insights. □

Lemma 1.

The parity of the unstable eigenvalues of

- P

needs to be consistent with that of A to make

A_{C}

in Equation (31) be Hurwitz.

Proof.

From Equation (31), one has det(

A_{C}

) = det(A)det(

- P

); hence, sign(det(

A C

)) = sign(det(A)det(

- P

)). If all eigenvaluses of

A_{C}

have a negative real part, then sign(det(

A_{C}

)) =

{(- 1)}^{2 n}

, which leads to sign(det(A)det(

- P

)) =

{(- 1)}^{n - l} {(- 1)}^{n - r},

wherel represents the number of unstable eigenvalues ofA, andr represents the corresponding numbers of

- P

, i.e.,

l + r

is an even number.

Therefore, to stabilize the normal form Andronov–Hopf bifurcation based on information of Equation (29),

- P

needs to be a 2 × 2 matrix. The Laplace transform between state and control gives

\frac{u (s)}{x (s)} = K \frac{s I}{s I + P} .

(38)

Obviously, when the state

x (s)

is one-dimensional, the washout filter has

d = P

, but for a multi-dimensional system, e.g.,P is 2 × 2, and the parameters of the washout filters are difficult to provide

\begin{matrix} \frac{s I}{s I + ε p^{*}} = \frac{s}{Δ} (\begin{matrix} s + ε p_{22}^{*} & - ε p_{12}^{*} \\ - ε p_{21}^{*} & s + ε p_{11}^{*} \end{matrix}), \\ w h e r e : \\ Δ = s^{2} + ε (p_{11}^{*} + p_{22}^{*}) s + ε^{2} (p_{11}^{*} p_{22}^{*} - p_{12}^{*} p_{21}^{*}) . \end{matrix}

(39)

In accordance with the form of washout filters, the higher-ordered terms in Equation (39) need to be simplified. Considering only the transient dynamics are concerned with washout filters, and based on the

L a p l a c e e x t r e m e t h e o r e m

,s in Equation (39) takes maximal value, which leads to the idea of the form of washout filters

\frac{s I}{s I + ε P^{*}} = (\begin{matrix} \frac{s}{s + ε (P_{11}^{*} + P_{22}^{*})} & 0 \\ 0 & \frac{s}{s + ε (P_{11}^{*} + P_{22}^{*})} \end{matrix}) .

(40)

Meaning two washout filters with

d = t r (P)

aid the controller design, and Equation (40) is decoupled. □

3.3. Bifurcation Analysis on the Closed-Loop System

Consider only one washout filter is introduced to

x_{i}

with feedback

u = K_{i} y

, and (

i, j

) pairing generates four possible cases for a two-input–two-output system. The normal form Andronov–Hopf bifurcation is augmented and

h (x, u)

explicitly presented

\{\begin{matrix} x^{'} = A x + B u + h (x, u) \\ z^{'} = x_{i} - d z = y \end{matrix} \Rightarrow \{\begin{matrix} x^{'} = (\begin{matrix} a & - 1 \\ 1 & a \end{matrix}) (\begin{matrix} x_{1} \\ x_{2} \end{matrix}) - (x_{1}^{2} + x_{2}^{2}) (\begin{matrix} x_{1} \\ x_{2} \end{matrix}) + K_{i} y \\ y^{'} = f_{i} (x_{i}, K_{i} y) - d y \end{matrix},

(41)

wherey is the washout-filter-output corresponding to

x_{i}

. Based on Routh–Hurwitz stability criterion [25], one can deduce the control gain

K_{i}

of the washout filter that stabilizes Equation (41), and the Jacobean of the augmented system gives

\tilde{A} = [\begin{matrix} (\begin{matrix} a & - 1 \\ 1 & a \end{matrix}) & {(B K)}_{i} \\ \frac{\partial f_{i}}{\partial x_{1}} \frac{\partial f_{i}}{\partial x_{2}} & {(B K)}_{i} - d \end{matrix}] .

(42)

Expending Equation (42) along the last column gives the characteristic polynomial

P_{c}

\begin{matrix} P_{c} & = d e t (\tilde{A} - λ I_{2 + 1}) \\ = ({(B K)}_{i} - d - λ) (λ_{- i} - λ) + {(\begin{matrix} \frac{\partial f_{i}}{\partial x_{i}} \end{matrix})}_{n} \cdot {(B K)}_{i} (λ_{- i} - λ) \\ = [({(B K)}_{i} - d - λ) (λ - λ_{n}) - {(\frac{\partial f_{i}}{\partial x_{1}})}_{n} \cdot {(B K)}_{i}] (λ_{- i} - λ) . \end{matrix}

(43)

One can design pairs of (

K_{i}, d

) to make eigenvalues of [∼] in Equation (43) stable, but

λ_{- i}

is kept unchanged. Meaning that, for the linearized system, dual-unstable eigenvalues need two washout filters.

However, Equation (41) is a nonlinear control, and the dynamics of

f_{i}

could be compensated through feedback. The Jacobian from Equation (42) ignores nonlinear correlations between system and control feedback, which leads to redundant control loops. When numerical bifurcation analysis is implemented, it is found that only one washout filter can stabilize the closed-loop system represented by Equation (41).

One can design controllers through the identification of critical points for the emergence/elimination of unstable eigenvalues. For the case of Equation (41), the limit cycles when

α > 0

need to be eliminated (or stabilized); hence, with the aid of a washout filter framed inFigure 2, Andronov–Hopf bifurcation points in the (

K_{i}, d

) domain can belocated numerically (refer toAppendix A for more details).

Substantial configurations are provided inFigure 3, where the solid lines are Andronov–Hopf bifurcation curves. BT represents Bogdanov–Takens bifurcation, and ZH represents zero-Hopf bifurcation, both of which are codimension-2 bifurcations originating from Andronov–Hopf bifurcation curves. Since

x_{1}

and

x_{2}

in Equation (41) are somewhat symmetric, configurationsC1 andC4 are identical, and so areC2 andC3. For each stabilizable region,

d > 0

implies the analog signal after washout filters are stable, while

d < 0

means that the feedback signal is unstable.

As shown inFigure 3, the Andronov–Hopf bifurcation curves separate out parameter space

(K_{i}, d)

that stabilizes the closed-loop system. Mention

d = 0

in cases ofC1 andC4 indicates that traditional PID control stabilizes the system. It is observed

d < 0

is able to stabilize the system from a nonlinear analysis perspective, but the unstable washout filter generates diverging signals which are dangerous for real-time applications where interference is unavoidable and is not recommended. Exclusively for

d > 0

, it is foundC2 andC3 have broader stable margins thanC1 andC4, which implies manipulating variables and control feedback could be chosen such that stabilizing control is easy to achieve.

Based on an analysis inFigure 3, numerical simulations on the Simulink platform were conducted, where

α = 0.1

for Equation (41). As is shown inFigure 4, when the parameter pair (−5, 0.5) is adopted for caseC1, it is stabilizable. The perturbation output after washout filter damps down and so as the states

x_{1}

and

x_{2}

, and a positive phase shift between y and

x_{1}

implies that the oscillatory dynamics of the control feedback and forward signal are well geared to stabilize the closed-loop system. Since

x_{1}

and

x_{2}

are coupled in the self-oscillatory system, stabilizing

x_{1}

would cause

x_{2}

to damp down. Further study of the information acquisition by the washout filter shows that, the stable one obtains oscillatory signals, while the unstable one is exponentially divergent, where the overall stable character of the closed-loop system is obtained by the subtraction of the forward and backward signals. Hence, to stabilize the self-oscillatory dynamic system, stable washout filters are suggested.

4. Case Study

4.1. Bifurcation Analysis on the Continuous Ethanol Fermentation Model

One of the most important bio-industrial sectors is the production of ethanol, which can be used as a gasoline additive to improve vehicle performance and reduce carbon emission. However, complex microorganism growth kinetics and ethanol inhibition present challenges for continuous production at the industrial scale, which limits the contribution of bio-fuel to current energy supply, and hence an exquisite controller design is critical to stabilizing the process. The structured model for

Z y m o m o n a s m o b i l i s

fermentation [26] in a continuous stirred reactor is provided as follows, which considers the mass balance of the substrate (

C_{s}

, g/L), microorganism (

C_{x}

, g/L),K-compartment biomass (

C_{e}

, g/L), product (

C_{p}

, g/L) and the related parameters are given inTable 1

\begin{matrix} \frac{d C_{s}}{d_{t}} = \frac{μ \cdot C_{S} \cdot C_{e}}{Y_{s x} \cdot (K_{S} + C_{S})} - m_{S} \cdot C_{x} + D \cdot (C_{s f} - C_{s}) \\ \frac{d C_{X}}{d t} = \frac{μ_{m a x} \cdot C_{s} \cdot C_{e}}{K_{s} + C_{s}} - D C_{x} \\ \frac{d C_{e}}{d_{t}} = \frac{(k_{1} - k_{2} C_{p} + k_{3} C_{P}^{2}) C_{s} \cdot C_{e}}{(K_{s} + C_{s})} - D C_{e} \\ \frac{d C_{p}}{d_{t}} = \frac{μ_{m a x} \cdot C_{s} \cdot C_{e}}{Y_{p x} \cdot (K_{s} + C_{s})} + m_{p} \cdot C_{x} - D C_{p} \end{matrix},

(44)

where the dilution rateD and the substrate feed concentration

C_{s f}

are set as the design /operation parameters because of practical concerns.D is the reciprocal of the residence time for the reaction stream, which has a fundamental impact on the bioreactor behavior;

C_{s f}

always comes from upstream, which may be interfered with due to upstream uncertainties.

With this model, nonlinear analysis is conducted on the open-loop system. A two-parameter bifurcation diagram decomposes the parameter space (

D, C_{s f}

) into different sections, and each has distinct dynamic properties. Furthermore, as is shown inFigure 5, the Andronov–Hopf bifurcation curve divides the parameter space into

S 1

and

S 2

. Any designing point at

S 1

generates a self-sustained oscillator, while the point at

S 2

is in stable steady equilibrium. Particularly, when

C_{s f}

is set as 140 g/L,H is detected at

D_{H} = 0.05024

^{- 1}

.H is a supercritical Hopf bifurcation point due to its negative first Lyapunov coefficient (

l p

), and operating the process withD at the left hand side ofH leads to self-oscillations, while stable equilibriums are obtained on the right-hand side. It is obvious that the highest

C_{p}

is achieved whenD is infinitesimal, which means that the reactor size needs to be very large. Since a high

C_{p}

of the bioreactor outlet cuts down the separation cost of later processing, one should design the reactor by settingD as the small reasonable equipment cost.

We set the design point for Equation (44) as

P_{d} (0.022, 140)

in this work, which falls into the self-oscillatory section. The idea of designing a self-oscillatory process was inspired by [27,28]. Taking

P_{d}

point as example, time average

C_{p}

8.8 %

higher than the corresponding steady-state output of the same design condition. However, the presence of self-oscillatory dynamics brings a challenge for the control of the continuous system; therefore, bifurcation control with the aid of washout filters is demanded to stabilize this multi-dimensional system.

4.2. Washout-Filter-Aided Control of the Fermentation Model

The model in Equation (44) is four-dimensional, and dual-unstable eigenvalues exist when Andronov–Hopf bifurcations emerge. For

P_{d}

, the equilibrium for

C_{s}, C_{x}, C_{e}

and

C_{p}

are 0.8509 g/L, 21.6547 g/L, 0.01915 g/L and 58.71 g/L, respectively; and the corresponding eigenvalues are

- 2, - 0.07136, 0.01232 \pm 0.1411

. Dual-unstable eigenspace emerges between states

C_{e}

and

C_{p}

, indicating that the controlling variable would be chosen among them,

C_{p}

is set as the washout filter input, and it is assumed that

C_{p}

can be detected with online measuring devices. Then, the washout filter output is provided

\frac{y}{C_{p}} = \frac{s}{s + d} \Rightarrow y^{'} = \frac{μ_{m a x} \cdot C_{s} \cdot C_{e}}{Y_{p x} \cdot (K_{s} + C_{s})} + m_{p} \cdot C_{x} - D C_{p} - d y,

(45)

whered is the filter parameter, and to promise thaty is stable and convergent,d needs to be positive. The manipulating variable is set toD, which is related to the input flow rate. Compared to the nominal

D_{0}

= 0.022 h

^{- 1}

, the feedback byP control is provided

D = D_{0} + k y,

(46)

wherek is the feedback gain.

Combining Equations (44)–(46), one could implement bifurcation analysis on the closed-loop system. As is shown inFigure 6, when the washout-filter-aided P-control is facilitated on the continuous fermentation model, the controlling space (

k, d

) is separated into three sections by the Andronov–Hopf and neutral saddle curves, where only the upper-left one is stabilizable. Note that both Andronov–Hopf and neutral saddle bifurcations are codimension-1 bifurcations, and the bifurcations happening on these two curves are codimension-2 bifurcations. For the Andronov–Hopf bifurcation curve, three codimension-2 bifurcations, GH1 (0.1745, 0.0024), GH2 (0.0623, 0.0006) and ZH (0, 0.0005) are detected; and for the neutral saddle bifurcation curve, HH (0.02221, −0.03037) is detected. Mention for a given filter parameter, e.g.,

d = 2

, with the decrease ink, that the dynamics of the closed-loop system evolves from stable to saddle, and then diverges.

Further study is conducted by setting the control parameters to constant, e.g.,

d = 2

and

k = 1

, and observing the changes in the dynamics. As shown inFigure 7, between the two Andronov–Hopf bifurcation curves is the segment that can be stabilized by washout-filter-aided P-control. When the control is facilitated, stabilizing the designing point demands that

C_{p}

reaches the steady state, and then the washout filter outputy has a tendency to converge towards zero, and Equation (45) tends towards zero; hence, the closed-loop system has the same equilibrium as the open-loop one, and

C_{p} = 58.71

g/L is preserved at operational point

P_{d}

. As provided inFigure 8, the simulation results show converged dynamics after control

d = 2

and

k = 1

, where

C_{p}

andD reach

P_{d}

as time approaches infinity.

5. Conclusions

In the presented study, we have focused on a subtle case with bifurcations, which is always un-linearizable and exacerbates information acquisition and a control problem. Taking the system with Andronov–Hopf bifurcation as a typical example, normal form analysis reveals that the cubic terms are the integral part of the whole system, and with the consideration of the un-linearizable characteristics, the stabilization procedure adopting washout-filter-aided control could realize the migration of dual-unstable eigenvalues with only one washout filter, which is not possible by previous control strategies. Moreover, the continuous ethanol fermentation process is taken as a case study for stabilizing a high-dimensional Andronov–Hopf bifurcation system: (1) variables mapping to the dual-unstable eigenspace are chosen as the manipulated variable; and (2) a stable washout filter parameter is demanded to stabilize the dual-unstable eigenvalues. By washout-filter-aided P-control, the ethanol fermentation process is stabilized and the equilibrium is preserved.

Author Contributions

All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by C.Z., C.Y. and J.N. The first draft of the manuscript was written by C.Z. and both authors commented on previous versions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge the National Natural Science Foundation of China-Yunnan basic research project (Grant No. 202001AU070048) for their support.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Nomenclature

D	dilution rate, h $^{- 1}$
$K_{S}$	saturation constant, g/L
$k_{1}, k_{2}, k_{3}$	inhibition constant, L/g
$T_{0}$	inherent period, h
$P_{d}$	design point for unforced system
$r_{e}$	growth rate of the K-compartment, g/L h
$C_{s f}$	substrate feed concentration, g/L
t	time, h
x	biomass concentration, g/L
$Y_{s x}, Y_{p x}$	yield factor of biomass on substrate/product, g/g
$m_{s}$	maintenance item, g/g h
Greek letters
$δ$	Greek letters
$μ_{m a x}$	maximum specific growth rate, h $^{- 1}$
Subscripts
s	substrate
x	biomass
e	K-compartment
p	product

Appendix A

The following test function (Stephanos theorem) locates the critical points numerically. Stephanos theorem: If a

n \times n

matrixA has eigenvalues

μ_{1}, μ_{2} \dots μ_{n}

ordered in an arbitrary order, then

2 A Θ I_{n}

satisfies

{(μ_{i} + μ_{j})}_{n > i > j} > 1

, and the function

φ_{H} (x, α) = d e t (2 f_{x} (x, α) Θ I_{n})

vanishes at an Andronov–Hopf bifurcation point, where “det” represents the determinant, and

Θ

is the bi-alternate product.

Proof.

The eigenvalues of a matrix are preserved by similarity transforms as follows

(P Θ P) (2 A Θ I_{n}) {(P Θ P)}^{- 1} = 2 B Θ I_{n}

whereB is the transform of

A, B = P A P^{- 1}

. Hence, one can assume thatA is in upper triangular form (by reducing to the Jordan form if necessary). Then, the eigenvalues ofA become its diagonal elements. By applying the equation below,

{(2 A Θ I_{n})}_{(i, j) (k, l)} = α_{i k} \{\begin{matrix} - α_{i l} & i f k = j \\ α_{i k} & i f k \neq i a n d l = j \\ α_{i i} + α_{j j} & i f k = i a n d l = j \\ α_{j l} & i f k = i a n d l \neq j \\ - α_{j k} & i f l = i \\ 0 & e l s e \end{matrix}

2 A Θ I_{n}

is an upper matrix if its indices are ordered lexicographically, and its diagonal elements are given as Stephanos theorem claims. By defining the test equation

φ_{H} (x, α) = 0

, the critical point

α_{c}

can be located and inserted into the continuation sequence. □

References

Wang, X.; Zhai, C.; Zhao, Z.; Palazoglu, A.; Sun, W.; Ai, J. Analysis of the Onset of Chaos for the Belousov-Zhabotinsky Reaction.Comput. Aided Chem. Eng.2018,44, 1375–1380. [Google Scholar]
Pashin, D.; Bastrakova, M.; Satanin, A.; Klenov, N. Bifurcation Oscillator as an Advanced Sensor for Quantum State Control.Sensors2022,22, 6580. [Google Scholar] [CrossRef] [PubMed]
Ambrożkiewicz, B.; Syta, A.; Georgiadis, A.; Gassner, A.; Meier, N. Experimental Verification of the Impact of Radial Internal Clearance on a Bearing’s Dynamics.Sensors2022,22, 6366. [Google Scholar] [CrossRef] [PubMed]
Zhai, H.; Sands, T. Comparison of Deep Learning and Deterministic Algorithms for Control Modeling.Sensors2022,22, 6362. [Google Scholar] [CrossRef] [PubMed]
Xu, C.; Liu, Z.; Liao, M.; Li, P.; Xiao, Q.; Yuan, S. Fractional-order bidirectional associate memory (BAM) neural networks with multiple delays: The case of Hopf bifurcation.Math. Comput. Simul.2021,182, 471–494. [Google Scholar] [CrossRef]
Yang, R.; Nie, C.; Jin, D. Spatiotemporal dynamics induced by nonlocal competition in a diffusive predator-prey system with habitat complexity.Nonlinear Dyn.2022,110, 879–900. [Google Scholar] [CrossRef]
Krishan, R.; Verma, A. Assessment and Enhancement of Hopf Bifurcation Stability Margin in Uncertain Power Systems.Electr. Power Syst. Res.2022,206, 107783. [Google Scholar] [CrossRef]
Alfifi, H.Y. Stability and Hopf bifurcation analysis for the diffusive delay logistic population model with spatially heterogeneous environment.Appl. Math. Comput.2021,408, 126362. [Google Scholar] [CrossRef]
Shi, J.; He, K.; Fang, H. Chaos, Hopf bifurcation and control of a fractional-order delay financial system.Math. Comput. Simul.2022,194, 348–364. [Google Scholar] [CrossRef]
Fiedler, B.; Flunkert, V.; Georgi, M.; Hövel, P.; Schöll, E. Refuting the odd-number limitation of time-delayed feedback control.Phys. Rev. Lett.2007,98, 114101. [Google Scholar] [CrossRef] [Green Version]
Zhai, C.; Sun, W. Analytical approximation of a self-oscillatory reaction system using the Laplace-Borel transform.Chaos Solitons Fractals2021,142, 110508. [Google Scholar] [CrossRef]
Cheng, Z. Anti-control of Hopf bifurcation for Chen’s system through washout filters.Neurocomputing2010,73, 3139–3146. [Google Scholar] [CrossRef]
Hassouneh, M.A.; Lee, H.C.; Abed, E.H. Washout filters in feedback control: Benefits, limitations and extensions. In Proceedings of the 2004 American Control Conference, Boston, MA, USA, 30 June–2 July 2004; Volume 5, pp. 3950–3955. [Google Scholar]
Kang, W. Bifurcation control via state feedback for systems with a single uncontrollable mode.SIAM J. Control Optim.2000,38, 1428–1452. [Google Scholar] [CrossRef]
Krener, A.J.; Kang, W.; Chang, D.E. Control bifurcations.IEEE Trans. Autom. Control2004,49, 1231–1246. [Google Scholar] [CrossRef]
Abed, E.H. A simple proof of stability on the center manifold for Hopf bifurcation.SIAM Rev.1988,30, 487–491. [Google Scholar] [CrossRef] [Green Version]
Abed, E.H.; Fu, J.H. Local feedback stabilization and bifurcation control, I. Hopf bifurcation.Syst. Control Lett.1986,7, 11–17. [Google Scholar] [CrossRef] [Green Version]
Abed, E.H.; Fu, J.H. Local feedback stabilization and bifurcation control, II. Stationary bifurcation.Syst. Control Lett.1987,8, 467–473. [Google Scholar] [CrossRef] [Green Version]
Abed, E.H.; Wang, H.O.; Chen, R.C. Stabilization of period doubling bifurcations and implications for control of chaos.Phys. D Nonlinear Phenom.1994,8, 154–164. [Google Scholar] [CrossRef]
Laiou, M.C.; Mönnigmann, M.; Marquardt, W. Stabilization of nonlinear systems by bifurcation placement.IFAC Proc. Vol.2004,37, 675–680. [Google Scholar] [CrossRef]
Zhang, N.; Qiu, T.; Chen, B. Bifurcation control and eigenstructure assignment in continuous solution polymerization of vinyl acetate.Chin. J. Chem. Eng.2015,23, 1523–1529. [Google Scholar] [CrossRef]
Wang, H.; Zhang, N.; Qiu, T.; Zhao, J.; He, X.; Chen, B. Optimization of a continuous fermentation process producing 1,3-propane diol with Hopf singularity and unstable operating points as constraints.Chem. Eng. Sci.2014,116, 116. [Google Scholar] [CrossRef]
Ding, L.; Hou, C. Stabilizing control of Hopf bifurcation in the Hodgkin-Huxley model via washout filter with linear control term.Nonlinear Dyn.2010,60, 131–139. [Google Scholar] [CrossRef]
Suo, L.; Ren, J.; Zhao, Z.; Zhai, C. Study on the nonlinear dynamics of the continuous stirred tank reactors.Processes2020,8, 1436. [Google Scholar] [CrossRef]
Hurwitz, A. On the conditions under which an equation hasonly roots with negative real parts.Math. Ann1895,46, 273–284. [Google Scholar] [CrossRef]
Jobses, I.M.L.; Egberts, G.T.C.; Luyben, K.C.A.M.; Roels, J.A. Fermentation kinetics of Zymomonasmobilis at high ethanol concentrations: Oscillations in continuous cultures.Biotechnol. Bioeng.1986,28, 868–877. [Google Scholar] [CrossRef]
Douglas, J.M.; Rippin, D.W.T. Unsteady state processoperation.Chem. Eng. Sci1966,21, 305–315. [Google Scholar] [CrossRef]
Bailey, J.E. Periodic operation of chemical reactors: A review.Chem. Eng. Commun.1974,1, 111–124. [Google Scholar] [CrossRef]

Figure 1. Phase portrait of (a) linearized system, (b) supercritical; and (c) subcritical Andronov–Hopf bifurcations, where

α = 0.1

, and the arrows point out the time evolution direction.

Figure 1. Phase portrait of (a) linearized system, (b) supercritical; and (c) subcritical Andronov–Hopf bifurcations, where

α = 0.1

, and the arrows point out the time evolution direction.

Figure 2. Block diagram for bifurcation control with washout-filter-aided stabilization.

Figure 3. Analysis on the closed system of normal form Andronov–Hopf bifurcation.

Figure 4. Time series of the closed-loop system, where the washout-filter-aided feedback parameters are

d = 0.5, k =

−5.

Figure 4. Time series of the closed-loop system, where the washout-filter-aided feedback parameters are

d = 0.5, k =

−5.

Figure 5. Bifurcation analyses on the ethanol fermentation process: two-parameter continuation diagram (left); and one-parameter continuation diagram (right).

Figure 6. Two–parameter continuation diagram for the control variables, whereGH1 andGH1 are the general Hopf bifurcation,ZH is the zero Hopf bifurcation andHH is the Hopf–Hopf bifurcation.

Figure 7. Two-parameter continuation diagram for the designing parameters.

Figure 8. Stabilized outcomes for the closed-loop system, where

d = 2

and

k = 1

Figure 8. Stabilized outcomes for the closed-loop system, where

d = 2

and

k = 1

Table 1. Parameters used for the continuous bio-ethanol reactor model.

Parameter	Definition	Value
$μ_{m a x}$	Maximum specific biomass growth rate (h $^{- 1}$ )	1.0
$Y_{s x}$	Yield factor for substrate to biomass (g/g)	0.00244498
$Y_{p x}$	Yield factor for substrate to product (g/g)	0.00526315
$m_{s}$	Maintenance substrate consumption rate (g/g h $^{- 1}$ )	2.16
$K_{s}$	Monod constant (g/L)	0.5
$k_{1}$	Empirical constant (h $^{- 1}$ )	16
$k_{2}$	Empirical constant (L/g h $^{- 1}$ )	0.497
$k_{3}$	Empirical constant (L $^{2}$ /g $^{2}$ h $^{- 1}$ )	0.000383
$m_{p}$	Maintenance product consumption rate (g/g h $^{- 1}$ )	1.1
$C_{s f}$	Substrate feed concentration (g/L)	140

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zhai, C.; Yang, C.; Na, J. Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters.Sensors2022,22, 9334. https://doi.org/10.3390/s22239334

AMA Style

Zhai C, Yang C, Na J. Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters.Sensors. 2022; 22(23):9334. https://doi.org/10.3390/s22239334

Chicago/Turabian Style

Zhai, Chi, Chunxi Yang, and Jing Na. 2022. "Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters"Sensors 22, no. 23: 9334. https://doi.org/10.3390/s22239334

APA Style

Zhai, C., Yang, C., & Na, J. (2022). Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters.Sensors,22(23), 9334. https://doi.org/10.3390/s22239334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further detailshere.

Article Metrics

Article Access Statistics

For more information on the journal statistics, clickhere.

Multiple requests from the same IP address are counted as one view.

Movatterモバイル変換

Article Menu

Bifurcation Control on the Un-Linearizable Dynamic System via Washout Filters

Abstract

1. Introduction

2. The Un-Linearizable System

2.1. Properties of the Linear System

2.2. Normal Form Andronov–Hopf Bifurcation: A Typical Un-Linearizable System

3. Stabilizing Normal Form Andronov–Hopf Bifurcation

3.1. Washout Filter as a State Manipulator

3.2. Feedback Stabilization Based on the Linearized System

3.3. Bifurcation Analysis on the Closed-Loop System

4. Case Study

4.1. Bifurcation Analysis on the Continuous Ethanol Fermentation Model

4.2. Washout-Filter-Aided Control of the Fermentation Model

5. Conclusions

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

Nomenclature

Appendix A

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI