Next Article in Journal
Motion along a Space Curve with a Quasi-Frame in Euclidean 3-Space: Acceleration and Jerk
Next Article in Special Issue
A Comparative Study of the Genetic Deep Learning Image Segmentation Algorithms
Previous Article in Journal
A DVA-Beam Element for Dynamic Simulation of DVA-Beam System: Modeling, Validation and Application
Previous Article in Special Issue
Analysis of Tempered Fractional Calculus in Hölder and Orlicz Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Analysis and Control for a Bioreactor in Fractional Order

by
Angelo M. Tusset
1,*,
Danilo Inacio
2,
Maria E. K. Fuziki
3,
Priscilla M. L. Z. Costa
1 and
Giane G. Lenzi
4
1
Department of Production Engineering, Federal University of Technology-Paraná, Paraná-Doutor Washington Subtil Chueire St. 330, Ponta Grossa 84017-220, Brazil
2
Department of Electrical Engineering, Federal University of Technology-Paraná, Paraná-Doutor Washington Subtil Chueire St. 330, Ponta Grossa 84017-220, Brazil
3
Department of Chemical Engineering, State University of Maringá, Colombo Ave. 5790, Maringá 87020-900, Brazil
4
Department of Chemical Engineering, Federal University of Technology-Paraná, Paraná-Doutor Washington Subtil Chueire St. 330, Ponta Grossa 84017-220, Brazil
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(8), 1609; https://doi.org/10.3390/sym14081609
Submission received: 21 June 2022 / Revised: 26 July 2022 / Accepted: 2 August 2022 / Published: 4 August 2022

Abstract

:
In this paper, a mathematical model was developed to describe the dynamic behavior of a bioreactor in which a fermentation process takes place. The analysis took into account the bioreactor temperature controlled by the refrigerant fluid flow through the reactor jacket. An optimal LQR control acting in the water flow through a jacket was used in order to maintain the reactor temperature during the process. For the control design, a reduced-order model of the system was considered. Given the heat transfer asymmetry observed in reactors, a model considering the fractional order heat exchange between the reactor and the jacket using the Riemann–Liouville differential operators was proposed. The numerical simulation demonstrated that the proposed control was efficient in maintaining the temperature at the desired levels and was robust for disturbances in the inlet temperature reactor. Additionally, the proposed control proved to be easy to apply in real life, bypassing the singularity problem and the difficulty of initial conditions for real applications that can be observed when considering Riemann–Liouville differential operators.

1. Introduction

Bioreactors are commonly used in the industry for biochemical oxygen demand (BOD) and chemical oxygen demand (COD) reduction and ethanol production from glucose fermentation [1,2].
In alcoholic fermentation, glucose is degraded into ethanol and carbon dioxide through an anaerobic process. In the oxidative phase, some aromas are formed, while in the fermentation phase, the yeast metabolism converts sugars into ethyl alcohol and carbon dioxide and Saccharomyces cerevisiae, the yeast that is commonly used, as it is easy to use, and its biomass is generally reused [2,3]. In this process, temperature is an important variable to consider, since it acts directly in the fermentation efficiency and denaturation process of microorganism proteins in bioreactors [4].
Regarding the optimization of ethanol production in bioreactors, in addition to temperature control, it is necessary to monitor operating parameters such as the pH of the medium, the concentration of oxygen, the concentration of glucose, etc. [2]. Taking into account these variables, it is possible to obtain mathematical models for the fermentation process, and thus to design a control system that optimizes the production of ethanol [4,5,6,7].
A hybrid mathematical model for biological reactors was proposed by [8], to optimize and control microbial processes, by considering cellular intelligence. In [9], a fractional order model was proposed for a bioreactor system, considering a fractional differential equation for ‘biomass’ and an integer-order differential equation for the ‘substrate’ dynamics. For numerical analysis of gas and heat production by bioreactors, a mathematical model was proposed in [10].
Mathematical models that describe the behavior of bioreactors have helped in process control projects, enabling the optimization of ethanol production. In [7], a fractional order PID control was proposed for the reactor temperature control. In [11], non-linear control systems were considered seeking to preserve the influences of the system’s non-linearities. For control of chemical reactors in [12], a cascade control was proposed. Optimal control design was used by [13] in the controlling of the reactor jacket temperature. A PI-fuzzy controller was designed by [14] to regulate the temperature inside a fermentation tank. In [4], a control system for the fermentation process was proposed considering the use of an adaptive linear neural network (ADALINE) to control the process. A control system for a continuous bioreactor with an unknown reaction rate term subjected to input saturation is proposed in [15]. In [16], a PID controller for integrating systems with time delay was proposed using the direct synthesis method and multiple dominant pole-placement approach. The proposed control was applied in a nonlinear continuous stirred tank reactor model.
In [17], the non-linear control SDRE (state-dependent Riccati equation) was used to control the temperature of a bioreactor, by controlling the flow of liquid in the reactor jacket. A fractional order control system was proposed in [18] for a nonlinear mathematical model of a bioreactor modeled in fractional order.
Fractional order dynamic systems have shown a huge application perspective in several areas, such as control techniques, artificial intelligence, network science, mechanics, physical sciences, economics, biological medical treatment, physical systems, and so on [19,20,21,22,23,24,25,26]. In [27,28,29], the authors considered the application of an adaptive control to the chaotic behavior of neural network systems in fractional order. In [30], the sliding mode technique was applied to investigate the adaptive fuzzy finite-time backtracking control for a nonlinear fractional model. In [31], a delayed feedback controller was used to suppress the chaotic behavior of a fractional order hybrid optical model. Numerical simulations demonstrated the effectiveness of the proposed control. In [32], we proposed a structure of distributed interval observers for multi-agent systems of fractional order considering non-linearities in the presence of nonlinearity. Numerical simulations were presented to demonstrate the effectiveness of the proposed distributed range observer. In [33], a fractional fuzzy PID control was applied in the temperature maintenance of industrial systems. Numerical results showed that the proposed control system was robust for environmental changes caused by internal and external disturbances. In [34], a hybrid controller called controller (FO-F-PID) was proposed, composed of fractional order, fuzzy and proportional–integral–derivative components. The proposed control was applied in a power generation system, involving different independent power-producing units. In [35], the dynamic analysis and control of a generalized prey-predator delay stage structure model with a fear and prey refuge effect was performed. Numerical simulations demonstrated the influence of the order of the fractional derivative on the dynamics of the system and that the proposed control was efficient in controlling the emergence of Hopf bifurcation. In [36], a numerical investigation of a mathematical model of fractional order of diabetes and its complications was presented. In [37], the authors presented a proposal for the solutions of non-homogeneous fractional integral equations. In [38], it was possible to establish the existence and uniqueness of the boundary value problem solution for the non-integer variant of the classical Atkinson theorem in the oscillation of the Emden–Fowler equations, considering the Riemann–Liouville order derivative.
In this context, the objective of this work was to design a linear control system for the cooling fluid flow that passes through the fermenter jacket, in order to maintain the ideal temperature in the bioreactor, considering the asymmetric heat transfer effects observed in reactors [39]. Derivatives of fractional order were included in the model to represent the asymmetric heat transfer and hysteresis effect of temperature variation. The Riemann-Liouville differential operators were considered in this paper in the numerical simulations, as they are among the most popular and well-known operators, facilitating the reproduction of results by other researchers [40].
We considered the technique of linear control by the LQR (linear quadratic regulator) control and included fractional derivatives in the reactor and jacket temperature equations, thus including non-linearity in the heat exchange between the reactor and jacket.
The LQR control aims to design a practical control system that provides the desired operational performance and minimizes the operating cost. Moreover, according to [41,42], the LQR control has proven to be efficient in reactor control applications. The representation of heat transfer by fractional order models has received a lot of attention from researchers, due to its ability to more adequately represent the system’s non-linearity and memory effects [43,44,45]. Fractional order models have already been used in various fields of science and technology, such as heat conduction [46,47,48,49]. This paper contributes to the understanding of the reduced order control system, considering, in the LQR control design, only the equations that represent the temperature variation in the reactor and in the jacket. Such a strategy allows for easier control in real implementation and makes it possible to control the system efficiently, even with the inclusion of fractional derivatives in the reactor and jacket temperature equations, bypassing possible problems of singularities and determination of initial conditions observed in Riemann–Liouville differential operator applications. These problems are solved since the proposed control does not depend on fractional derivatives, making the system completely controllable.
This paper is organized as follows: Section 2 presents the mathematical model used and the system dynamics in the case in which reactor temperature is not controlled. In Section 3, the reactor temperature control is proposed, considering the application of the LQR control in the cooling liquid flow through the jacket, presenting the system dynamics with control, and the analysis of robustness for disturbances in the reactor inlet temperature. Section 4 presents the mathematical model considering both reactor and jacket temperature equations with fractional order derivatives. Numerical results of the system dynamics in fractional order, the efficiency of the proposed control for the system in fractional order, and its robustness for disturbances in the reactor inlet temperature are presented. Finally, the article is completed in Section 5.

2. Mathematical Model

Figure 1 shows the continuous fermentation reactor representation, with a stirring system and cooling jacket, and the variables considered in this paper.
Where CO2, Cs, Cx and Cp are the concentration of oxygen, glucose, biomass and ethanol, respectively. Tag, Tin and Tr, are the temperatures in the jacket, reactor inlet, and inside the reactor. Fag is the coolant flow, Fi is the fermentation vessel inlet flow, and Fe is the fermenter outlet flow [4,17].

2.1. Mass Balance Equations

Considering Figure 1, it is possible to determine the mass and energy balance equations for the bioreactor system applied to ethanol production, assessing the influence of temperature, oxygen and the yeast concentrations [17].
For the mass balance equations, it was assumed that the variation in the microorganism growth in the bioreactor depends on the yeast growth rate and the microorganism concentration in the outlet.
d ( C x ) d t = μ x C x C s K s + C s e K p C p F e V C x
where µ x is the specific growth rate (h−1) defined as follows [17]:
µ x = A 1 e E a 1 R T r + 273 A 2 e E a 2 R T r + 273
The variation of glucose concentration is given by:
d ( C S ) d t = 1 R S X µ x C x C S K s + C S e K p C p 1 R S P μ P C x C S K S 1 + C S e K p 1 C p + F i V C S , i n F e V C S
The variation in the concentration of ethanol inside the reactor can be related to the amount of ethanol produced and the fermenter output, and is defined as:
d ( C P ) d t = µ P C x C S K S 1 + C S e K P 1 C P F e V C P
According to [7], the maximum ethanol production depends on the ideal level of dissolved oxygen in the fermentation process; therefore, the rate of oxygen variation is given by the following equation:
d ( C O 2 ) d t = k l a C O 2 * C O 2 r O 2
where k l a = k l a 0 1.024 T r 20 .
The following relationship gives the oxygen concentration in the reactor:
C O 2 * = 14.6 0.3943 T r + 0.00714 T r 2 0.0000646 T r 3 10 H i I i
The equilibrium oxygen level in the fermenter depends on the temperature and the coefficient of ionic forces, given by:
H i I i = 0.5 H N a m N a C l M N a C l M N a V + 2 H C a m C a C O 3 M C a C O 3 M C a V + 2 H M g m M g C l 2 M M g C l 2 M M g V + 0.5 H C l m N a C l M N a C l + 2 m M g C l 2 M M g C l 2 M C l V + 2 H C O 3 m C a C O 3 M C a C O 3 M C O 3 V + 0.5 H H 10 p H + 0.5 H O H 10 14 p H
However, during the fermentation process in the reactor, the yeast will consume oxygen, the consumed oxygen being defined as:
r O 2 = μ O 2 C x Y O 2 C O 2 K O 2 + C O 2
The reactor temperature is given by equation:
d ( T r ) d t = F i V T i n + 273 F e V T r + 273 + r O 2 Δ H r 32 ρ r C h e a t , r K T A T T r T a g V ρ r C h e a t , r
The temperature variation of the cooling jacket:
d ( T a g ) d t = F a g V j T i n , a g T a g + K T A T T r T a g V j ρ a g C h e a t , a g
where V j is the jacket volume (L).
Considering Equations (1), (3)–(5), (9) and (10), we obtain the following system of differential equations:
d ( C x ) d t = μ x C x C s K s + C s e K p C p F e V C x d ( C S ) d t = 1 R S X μ x C x C S K s + C S e K p C p 1 R S P μ P C x C S K S 1 + C S e K p 1 C p + F i V C S , i n F e V C S d ( C P ) d t = μ P C x C S K S 1 + C S e K P 1 C P F e V C P d ( C O 2 ) d t = k l a C O 2 * C O 2 r O 2 d ( T r ) d t = F i V T i n + 273 F e V T r + 273 + r O 2 Δ H r 32 ρ r C h e a t , r K T A T T r T a g V ρ r C h e a t , r d ( T a g ) d t = F a g V j T i n , a g T a g + K T A T T r T a g V j ρ a g C h e a t , a g

2.2. Numerical Simulations

In numerical simulations, we considered the system of differential equations previously described in Equation (11), fourth order Runge-Kutta method, with integration step (h = 0.01), and parameters: A1 = 9.5 × 108, A2 = 2.55 × 1033, AT = 1, Cheat,ag = 4.18, Cheat,r = 4.18, Ea1 = 55,000, Ea2 = 220,000, kla0 = 38, Ks = 1.03, Ks1 = 1.68, KT = 3.6 × 105, R = 8.31, RSP = 0.435, Rsx = 0.607, HCa = −0.303, HCl = 0.844, HCO3 = 0.485, HH = −0.774, HMg = −0.314, HNa = −0.550, HOH = 0.941, mCaCO3 = 100, mMgCl2 = 100, mNaCl = 500, MCa = 40, MCaCO3 = 90, MCl = 35.5, MCO3 = 60, YO2 = 0.97, ΔHr = 518, µO2 = 0.5, µP = 1.79, ρag = 1000, ρr = 1080 [4,15].
In Figure 2, the system (11) behavior without control can be seen, considering Tin and Faq.
Figure 2 shows that the fermenter took approximately 100 h to stabilize and maintain its constant output. The results also show that in order to increase ethanol production, it is necessary to increase the yeast concentration, with more consumption of oxygen and glucose.

3. Temperature Control Strategy by LQR Control

Figure 3 shows the bioreactor representation considering the introduction of coolant flow control in the jacket.
Figure 3 shows that the control strategy of the reactor temperature acts in the control of the coolant flow control inside the jacket, thus allowing the cooling jacket to maintain the temperature inside the reactor by heat exchange.
Considering that, the main objective is to control the reactor temperature, which can simplify, for control project, the system from six differential equations to two differential equations, taking into account only the reactor and jacket temperature equations. The strategy used is similar to the one used successfully by [50].
The system can be considered as follows:
d ( T r ) d t = F e V + K T A T V ρ r C h e a t , r T r + K T A T V ρ r C h e a t , r T a g + F i V T i n + 273 273 F e V + r O 2 Δ H r 32 ρ r C h e a t , r d ( T a g ) d t = K T A T V j ρ a g C h e a t , a g T r K T A T V j ρ a g C h e a t , a g T a g F a g V j T a g + F a g V j T i n , a g
The system (12) can be presented by the following matrix form:
d ( X ) d t = A X + G + B U
where X = T r T a g , A = F e V + K T A T V ρ r C h e a t , r K T A T V ρ r C h e a t , r K T A T V j ρ a g C h e a t , a g K T A T V j ρ a g C h e a t , a g , B = 0 T a g V j + T i n , a g V j , U = F a q and G = F i V T i n + 273 273 F e V + r O 2 Δ H r 32 ρ r C h e a t , r 0 , where U is the feedback control, it is defined as follows:
U = R 1 B T P e = K e
where e   =   [ X X * ] , X represents the system states and X * the desired states. The matrix P is obtained by solving the Ricatti equation:
A T P + P A P B R 1 B T P + Q = 0
The functional cost for the control problem U is given by:
J = 1 2 t 0 ( e T Q e + U T R U ) d t
where Q and R are positive definite matrices.
The controllability is given by:
M = [ B n × m A n × n B n × m ]
If rank of the matrix M is 2, the system (12) is controllable. Considering the matrices A, B and M, it is possible to obtain: A = 0 . 1307 0 . 0797 1 . 7225 1 . 7225 , B = 0 0.1800 and M = B A B = 0 0.1800 0 . 0144 0 . 3100 .
The rank of the matrix is 2, which ensures that the system (13) is controllable. Defining: Q = 35000 655000 100 100 1 and R = 200 . Considering F a g   =   U , and
U = k 1 , 1 T r T r * k 1 , 2 T a g T a g *
where T r * and T a g * are the desired temperatures for the reactor and jacket, and k 1 , 1 and k 1 , 2 are the gains of the control, obtained according to Equation (14).
According to [51,52], the optimum temperature for the alcoholic fermentation of Saccharomyces cerevisiae is 32 °C. In order to optimize alcoholic fermentation, the desired temperatures considered were T r * = 32   ° C and T a g * = 24   ° C . Solving Equation (15) and replacing T r * = 32   ° C and T a g * = 24   ° C in Equation (18), the control signal is obtained:
U = 10557.7595 T r 32 + 88.5189 T a g 24
Figure 4 shows the states’ variations with the proposed control: microorganism concentrations, glucose, alcohol and O2 concentrations, as well as the reactor temperature and jacket temperature variation, taking into account the bioreactor mathematical model for the fermentation and ethanol production represented by Equation (11), and considering coolant flow rate variation (Fag).
Figure 4 shows that the control maintained the ideal temperature at approximately 32 °C in the reactor, with less than ±5% error, which is considered to be an acceptable variation, according to [47]. As can be seen in Figure 4, in order to reduce the temperature of the reactor, it was necessary to rapidly increase the coolant flow in the jacket, reducing the temperature of the cooling agent and thus reducing the temperature of the reactor to 32 °C. The results also showed that when the reactor temperature was lower than 32 °C, the liquid flow was zero. Table 1 presents a comparison of the process variables with and without the temperature control.
As can be seen in Table 1, with the proposed control, it was possible to maintain the reactor temperature at 32 °C, which provided an increase of 16.29% in the microorganisms’ productivity in generating ethanol (productivity coefficient ( C p C x ): with 32 °C = 12.6 and 36 °C = 10), consuming 20.9% and 25.4% less glucose and O2.
Proposed Control Sensitivity to Variations in the Reactor Inlet Temperature
To analyze the control robustness, a variation in the reactor inlet temperature was applied in the following form:
T i t = 25 + 5 s i n π 12 t
In Figure 5, it is possible to observe the reactor temperature variation considering the perturbation according to Equation (20).
Figure 5 shows that the proposed control is robust for periodic variation in reactor inlet temperature (Tin), with less than 0.125% error. It can also be observed that the jacket temperature variation balances the reactor temperature close to the desired value, with variation between 21.04 and 27.33 °C. In addition, it is possible to observe that the coolant flow rate necessary to keep the reactor temperature close to the desired 32 °C varied between 29 and 157 L/h in a permanent regime case.

4. Dynamic Analysis and Control for the Fractional-Order Case

Given what was discussed in the present work until here, and taking into account that the heat exchange between the reactor and the jacket can be represented in fractional order, in the present section, both the reactor and jacket temperature equations were considered to have a fractional order derivative.
According to [53], differential equations may involve Riemann–Liouville differential operators of fractional order q > 0 , and they are generally considered in the following form [54,55,56]:
D q x ( t ) = 1 Γ m q t 0 t x ( m ) ( x ) t x q m + 1 d x
where m = q , which is the first integer not smaller than q.
The system (11) in fractional order is described as follows:
d q 1 C x d τ q 1 = μ x C x C s K s + C s e K p C p F e V C x d q 2 C S d τ q 2 = 1 R S X μ x C x C S K s + C s e K p C p 1 R S P μ P C x C S K S 1 + C s e K p 1 C p + F i V C S , i n F e V C S d q 3 C p d τ q 3 = μ P C x C S K S 1 + C s e K p 1 C p F e V C P d q 4 C O 2 d τ q 4 = ( k l a ) ( C O 2 * C O 2 ) r O 2 d q 5 T r d τ q 5 = F i V ( T i n + 273 ) F e V ( T r + 273 ) + r O 2 Δ H r 32 ρ r C h e a t , r K T A T ( T r T a g ) V ρ r C h e a t , r d q 6 T a g d τ q 6 = F a g V j ( T i n , a g T a g ) + K T A T ( T r T a g ) V j ρ a g C h e a t , a g
where q 1 = q 2 = q 3 = q 4 = 1 , q 5 and q 6 represent the order of the derivatives, and 0.3 < q 5 , 6 < 1 .
The numerical simulations for behavior analysis of the system (22) in fractional order uses the algorithm proposed by [33].
Figure 6 show the state variations considering the bioreactor mathematical model for fermentation and ethanol production with fractional order represented by Equation (22).
Considering the order variation of the fractional derivative in the reactor and in the jacket, it can be observed that for ( q 5 , 6 < 0.85 ), the reactor temperature lower is than the ideal temperature of 32 °C, and for ( q 5 , 6 > 0.85 ), the reactor temperature is higher than the desired value. We can also observe that the lower the value of ( q 5 , 6 ), the lower the temperature inside the reactor, implying a reduction in alcohol production and an increase in glucose and oxygen concentrations.
The results are in accordance with what is expected for the system in fractional order, given that, according to [9], the dynamic responses of a bioreactor in fractional order are inherently slower than those obtained from the integer-order model, and their velocity increases with the order of the system.

4.1. Proposal of the Temperature Control Strategy by LQR Control in Fractional Order Case

In Figure 7, it is possible to observe the stabilization point of the system (22) in the case in which the system is under the action of the proposed control of the coolant flow in the jacket, and the reactor and jacket temperature derivatives are represented in fractional order.
As can be seen in Figure 7, the proposed control is efficient in bringing the temperature of the reactor to 32 °C for the cases in which the temperature of the reactor is above the considered setpoint. For cases in which the temperature was lower than 32 °C, the control was null, since the liquid flow control in the jacket was only able to cool the reactor. In this case, it can be seen that the control acted only in cases in which the order was higher than ( q 5 = q 6 = 0.75 ).
In Figure 8, the control applied to the system (22) is shown considering ( q 5 = q 6 = 0.95 ), which are the values that bring the reactor temperature close to 36 °C, the temperature in which the yeast proteins begin to denature [52].
Analyzing the results presented in Figure 8, it can be concluded that the proposed control is efficient in bringing the temperature of the reactor to 32 °C, keeping the temperature of the reactor at the ideal value for alcohol production even if the system is represented by a fractional order model.

4.2. Proposed Control Sensitivity to Variations in the Reactor Inlet Temperature in Fractional Order Case

In Figure 9, it is possible to observe the reactor temperature variation considering the disturbance of the inlet temperature according to Equation (20) and reactor and jacket temperature in fractional order considering 0.3 q 5 , 6 1 .
As can be seen in the results presented in Figure 9, the control remained robust even when the system was represented in fractional order and suffered variations in the inlet temperature, as it was able to keep the reactor temperature at 32 °C, with and without disturbance. It could also be observed that the jacket temperature variation that balances the reactor temperature is similar to that already observed for the system with control in the integer-order model (Equation (11)).

5. Conclusions

Using numerical simulations, a dynamic analysis of a nonlinear mathematical model of an alcoholic fermentation reactor was presented. The use of the optimal LQR control made it possible to obtain a control system in a reduced order state space, allowing the maintenance of a constant reactor temperature by the variation of the water flow in the cooling jacket with only two fixed gains from the controller. In order to analyze the non-linearities of the heat exchange between the cooling jacket and the reactor, the use of a fractional order in the reactor temperature equation and in the coolant temperature equation was proposed.
The numerical results presented demonstrated that the proposed control strategy was efficient in keeping the interior temperature of the reactor at 32 °C and was robust for variations in the reactor inlet temperature, both in cases of fractional order or integer order. The numerical results are similar to those obtained experimentally as presented by [17], demonstrating the versatility of the proposed control to be used experimentally in future studies.
The results presented demonstrated that the proposed LQR control in reduced order contributes to the body of knowledge on temperature control of bioreactors, mainly considering the variation of the reactor temperature in fractional order. The presented control strategy can eliminate possible problems of singularities and problems of initial condition determination observed in the applications of the Riemann–Liouville differential operators for the cases in which control in fractional order is used in real systems. Such problems do not occur with the proposed control since it does not depend on fractional derivatives of the temperatures and given that systems in reduced order are completely controllable.

Author Contributions

Conceptualization, A.M.T. and D.I.; methodology, G.G.L.; software, A.M.T.; formal analysis, M.E.K.F.; investigation, P.M.L.Z.C.; resources, G.G.L.; writing—original draft preparation, P.M.L.Z.C.; writing—review and editing, M.E.K.F.; visualization, G.G.L.; supervision, A.M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the Capes, Fundação Araucária, and CNPq agency. The first author thanks CNPq for the financial support (process: 310562/2021-0).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AState matrix
BControl matrix
CsGlucose concentration in the reactor
CxBiomass concentration
CpEthanol concentration
CO2Oxygen concentration
Cs,inSubstrate concentration
FagFlow of the coolant
FeFermenter outlet flow
FeV−1Dilution rate of the reactor
FiFermentation vessel flow
KsGrowth substrate constant
KpGrowth inhibition constant by ethanol
KlaVolumetric mass transfer coefficient
MControllability matrix
PPressure
PRicatti Equation
q Derivative order
QPositive definite matrix
RPositive definite matrix
TTemperature
TagTemperature of the thermal jacket
T*agDesired temperature of the jacket
TrOutlet temperature
T*rDesired temperature of the reactor
TinReactor inlet temperature
UFeedback control
VReactor volume
VjJacket volume
VrReagent flow
VtTank volume
XStates vector
X*Desired states vector
µ x Specific growth rate
Abbreviations
ADALINEAdaptive linear neural network
BODBiochemical oxygen demand
CODChemical oxygen demand
CSTRContinuous stirred tank reactor
DOFDegree of freedom
FOPIDFractional order PID
IMCInternal model control
LDPELow density polyethylene
LQR Linear quadratic regulator
pHPotential of hydrogen
PIDProportional, integrative, derivative
SDREState-dependent Riccati equation

References

  1. Ahmed, F.N.; Lan, C.Q. Treatment of landfill leachate using membrane bioreactors: A review. Desalination 2012, 287, 41–54. [Google Scholar] [CrossRef]
  2. Basso, T.O.; de Kok, S.; Dario, M.; Espirito-Santo, J.C.A.; Müller, G.; Schlölg, P.S.; Silva, C.P.; Tonso, A.; Daran, J.M.; Gombert, A.K.; et al. Engineering topology and kinetics of sucrose metabolism in Saccharomyces cerevisiae for improved ethanol yield. Metab. Eng. 2011, 13, 694–703. [Google Scholar] [CrossRef] [PubMed]
  3. Muniz, J.B.; Marcelino, M.; Motta, M.D.; Schuler, A.; Alves da Motta, M. Influence of static magnetic fields on S. cerevisae biomass growth. Braz. Arch. Biol. Technol. 2007, l50, 515–520. [Google Scholar] [CrossRef]
  4. Pachauri, N.; Singh, V.; Rani, A. Two degree of freedom PID based inferential control of continuous bioreactor for ethanol production. ISA Trans. 2017, 68, 235–250. [Google Scholar] [CrossRef]
  5. Schubert, J.; Simutis, R.; Dors, M.; Havlik, I.; Lübbert, A. Bioprocess optimization and control: Application of hybrid modelling. J. Biotechnol. 1994, 35, 51–68. [Google Scholar] [CrossRef]
  6. Sultana, S.; Jamil, N.M.; Saleh, E.A.M.; Yousuf, A.; Faizal, C.K.M. A mathematical model for ethanol fermentation from oil palm trunk sap using Saccharomyces cerevisiae. J. Phys. Conf. Ser. 2017, 890, 012050. [Google Scholar] [CrossRef]
  7. Pachauri, N.; Rani, A.; Singh, V. Bioreactor temperature control using modified fractional order IMC-PID for ethanol production. Chem. Eng. Res. Des. 2017, 122, 97–112. [Google Scholar] [CrossRef]
  8. Patnaik, P. Hybrid Models for Biological Reactors: Performance and Possibilities. IREPHY 2018, 12, 1–7. [Google Scholar]
  9. Ahmad, W.; Abdel-Jabbar, N. Modeling and Simulation of a Fractional Order Bioreactor System. IFAC Proc. 2006, 39, 260–264. [Google Scholar] [CrossRef]
  10. Mohcine, A.; Gueraoui, K.; Mzred, Z.G.; Men-La-Yakhaf, S. Mathematic and Numerical Modeling of Biogas Production in the Bioreactive Plant for Valorizing Domestic Waste. Int. Rev. Mech. Eng. 2017, 11, 249–255. [Google Scholar] [CrossRef]
  11. Von Stosch, M.; Oliveira, R.; Peres, J.; Azevedo, S.F. A general hybrid semi-parametric process control framework. J. Process. Control 2012, 22, 1171–1181. [Google Scholar] [CrossRef]
  12. Romero-Bustamante, J.A.; Moguel-Castañeda, J.G.; Puebla, H.; Hernandez-Martinez, E. Robust Cascade Control for Chemical Reactors: An Approach based on Modelling Error Compensation. Int. J. Chem. React. Eng. 2017, 15, 1–15. [Google Scholar] [CrossRef]
  13. Yao, F.Z.; Lohi, A.; Upreti, S.R.; Dhib, R. Modeling, simulation and optimal control of ethylene polymerization in non-isothermal, high-pressure tubular reactors. Int. J. Chem. React. Eng. 2004, 2, 1–15. [Google Scholar] [CrossRef]
  14. Fonseca, R.R.; Schmitz, J.; Fileti, A.M.F.; Silva, F. A fuzzy–split range control system applied to a fermentation process. Bioresour. Technol. 2013, 142, 475–482. [Google Scholar] [CrossRef]
  15. Rincón, A.; Hoyos, F.E.; Candelo-Becerra, J.E. Conditioned Adaptive Control for an Uncertain Bioreactor with Input Saturation and Steep Settling Desired Output. Symmetry 2022, 14, 1232. [Google Scholar] [CrossRef]
  16. Zhang, W.; Cui, Y.; Ding, X. An Improved Analytical Tuning Rule of a Robust PID Controller for Integrating Systems with Time Delay Based on the Multiple Dominant Pole-Placement Method. Symmetry 2020, 12, 1449. [Google Scholar] [CrossRef]
  17. Bressan, D.C.; Ribeiro, M.A.; Lenzi, G.G.; Balthazar, J.M.; Tusset, A.M. A Note on SDRE Control Applied in the Fermentation Reactor. Int. Rev. Mech. Eng. 2019, 13, 576–586. [Google Scholar] [CrossRef]
  18. Estakhrouiyeh, M.R.; Vali, M.; Gharaveisi, A. Application of fractional order iterative learning controller for a type of batch bioreactor. IET Control. Theory Appl. 2016, 10, 1374–1383. [Google Scholar] [CrossRef]
  19. Tan, H.L.; Wu, J.W.; Bao, H.B. Event-triggered impulsive synchronization of fractional-order coupled neural networks. Appl. Math. Comput. 2022, 429, 127244. [Google Scholar] [CrossRef]
  20. Zhang, Z.; Wang, Y.N.; Zhang, J.; Ai, Z.Y.; Liu, F. Novel stability results of multivariable fractional-order system with time delay. Chaos Solitons Fractals 2022, 157, 111943. [Google Scholar] [CrossRef]
  21. Huang, C.D.; Wang, J.; Chen, X.P.; Cao, J.D. Bifurcations in a fractional-order BAM neural network with four different delays. Neural Netw. 2021, 141, 344–354. [Google Scholar] [CrossRef]
  22. Xu, C.J.; Liao, M.X.; Li, P.L.; Yuan, S. Impact of leakage delay on bifurcation in fractional-order complex-valued neural networks. Chaos Solitons Fractals 2021, 142, 110535. [Google Scholar] [CrossRef]
  23. Borah, M.; Das, D.; Gayan, A.; Fenton, F.; Cherry, E. Control and anticontrol of chaos in fractional-order models of Diabetes, HIV, Dengue, Migraine, Parkinson’s and Ebola virus diseases. Chaos Solitons Fractals 2021, 153, 111419. [Google Scholar] [CrossRef]
  24. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations to Methods of Their Solution and Some of Their Applications. In Mathematics in Science and Engineering; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  25. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, North-Holland Mathematical Studies; Elsevier (North-Holland) Science Publishers: Amsterdam, The Netherlands, 2006. [Google Scholar]
  26. Miller, K.S.; Ross, B. An Introduction to Fractional Calculus and Fractional Differential Equation; John Wiley and Sons: New York, NY, USA, 1993. [Google Scholar]
  27. Shi, J.P.; 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]
  28. Xu, C.J.; Zhang, W.; Aouiti, C.; Liu, Z.X.; Liao, M.X.; Li, P.L. Further investigation on bifurcation and their control of fractional order BAM neural networks involving four neurons and multiple delays. Math. Methods Appl. Sci. 2021, 2021, 1–24. [Google Scholar] [CrossRef]
  29. Xu, C.J.; Zhang, W.; Aouiti, C.; Liu, Z.X.; Yao, L.Y. Further analysis on dynamical properties of fractional-order bi-directional associative memory neural networks involving double delays. Math. Methods Appl. Sci. 2022, 2022, 1–10. [Google Scholar] [CrossRef]
  30. Xue, G.M.; Lin, F.N.; Li, S.G.; Liu, H. Adaptive fuzzy finite-time backstepping control of fractional-order nonlinear systems with actuator faults via command-filtering and sliding mode technique. Inf. Sci. 2022, 600, 189–208. [Google Scholar] [CrossRef]
  31. Li, P.; Gao, R.; Xu, C.; Li, Y. Chaos Suppression of a Fractional-Order Modificatory Hybrid Optical Model via Two Different Control Techniques. Fractal Fract. 2022, 6, 359. [Google Scholar] [CrossRef]
  32. Zhang, H.; Huang, J.; He, S. Fractional-Order Interval Observer for Multiagent Nonlinear Systems. Fractal Fract. 2022, 6, 355. [Google Scholar] [CrossRef]
  33. Liu, L.; Xue, D.; Zhang, S. General type industrial temperature system control based on fuzzy fractional-order PID controller. Complex Intell. Syst. 2021, 2021, 1–13. [Google Scholar] [CrossRef]
  34. Mahto, T.; Kumar, R.; Malik, H.; Khan, I.A.; Al Otaibi, S.; Albogamy, F.R. Design and Implementation of Frequency Controller forWind Energy-Based Hybrid Power System Using Quasi-Oppositional Harmonic Search Algorithm. Energies 2021, 14, 6459. [Google Scholar] [CrossRef]
  35. Lan, Y.; Shi, J.; Fang, H. Hopf Bifurcation and Control of a Fractional-Order Delay Stage Structure Prey-Predator Model with Two Fear Effects and Prey Refuge. Symmetry 2022, 14, 1408. [Google Scholar] [CrossRef]
  36. Srivastava, H.M.; Dubey, R.S.; Jain, M. A study of the fractional-order mathematical model of diabetes and its resulting complications. Math. Meth. Appl. Sci. 2019, 42, 4570–4583. [Google Scholar] [CrossRef]
  37. Kaewnimit, K.; Wannalookkhee, F.; Nonlaopon, K.; Orankitjaroen, S. The Solutions of Some Riemann–Liouville Fractional Integral Equations. Fractal Fract. 2021, 5, 154. [Google Scholar] [CrossRef]
  38. Baleanu, D.; Mustafa, O.G.; Agarwal, R.P. An existence result for a super linear fractional differential equation. Appl. Math. Lett. 2010, 23, 1129–1132. [Google Scholar] [CrossRef]
  39. Gosiewski, K.; Pawlaczyk-Kurek, A. Impact of Thermal Asymmetry on Efficiency of the Heat Recovery and Ways of Restoring Symmetry in the Flow Reversal Reactors. Int. J. Chem. React. Eng. 2019, 17, 20180021. [Google Scholar] [CrossRef]
  40. Atangana, A.; Secer, A. A Note on Fractional Order Derivatives and Table of Fractional Derivatives of Some Special Functions. Abstr. Appl. Anal. 2013, 2013, 279681. [Google Scholar] [CrossRef]
  41. Moghadam, A.A. LQ control of coupled hyperbolic PDEs and ODEs: Application to a CSTR-PFR system. IFAC Proc. 2010, 43, 721–726. [Google Scholar] [CrossRef]
  42. Sundari, S.; Nachiappan, A. Design of optimal Linear Quadratic Regulator for the stabilization of Continuous Stirred Tank Reactor (CSTR) Process. Int. J. Pure Appl. Math. 2018, 118, 1–19. [Google Scholar]
  43. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College: Bologna, Italy, 2010. [Google Scholar]
  44. Ortigueira, M.D. Fractional Calculus for Scientists and Engineers; Springer: Berlin/Heidelberg, Germany, 2011; Volume 84. [Google Scholar]
  45. Uchaikin, V. Fractional Derivatives for Physicists and Engineers, Background and Theory; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  46. Battaglia, J.; Lay, L.L.; Batsale, J.C.; Oustaloup, A.; Cois, O. Heat flux estimation through inverted non integer identification models. Int. J. Therm. Sci. 2000, 39, 374–389. [Google Scholar] [CrossRef]
  47. Ouhsaine, L.; Boukal, Y.; Ganaoui, E.M.; Darouachc, M.; Zasadzinskic, M.; Mimetb, A.; Radhyd, N.E. A general fractional-order heat transfer model for photovoltaic/termal hybrid systems and its observer design. Energy Procedia 2017, 139, 49–50. [Google Scholar] [CrossRef]
  48. Ezzat, M.A.; El-Bary, A.A. Effects of variable thermal conductivity and fractional order of heat transfer on a perfect conducting infinitely long hollow cylinder. Int. J. Therm. Sci. 2016, 108, 62–69. [Google Scholar] [CrossRef]
  49. Cecová, M.; Terpák, J. Fractional Heat Conduction Models and Thermal Diffusivity Determination. Math. Probl. Eng. 2015, 2015, 753936. [Google Scholar] [CrossRef]
  50. Tusset, A.M.; Piccirillo, V.; Balthazar, J.M.; Fonseca, R.M.L.R. On suppression of chaotic motions of a portal frame structure under non-ideal loading using a magneto-rheological damper. J. Theor. Appl. Mech. 2015, 53, 653–664. [Google Scholar] [CrossRef]
  51. Esfahanian, M.; Rad, A.S.; Khoshhal, S.; Najafpour, G. Mathematical modeling of continuous ethanol fermentation in a membrane bioreactor by pervaporation compared to conventional system: Genetic algorithm. Bioresour. Technol. 2016, 212, 62–71. [Google Scholar] [CrossRef]
  52. Amillastre, E.; Aceves-Lara, C.A.; Uribelarrea, J.L.; Alfenore, S.; Guillouet, S.E. Dynamic model of temperature impact on cell viability and major product formation during fed-batch and continuous ethanolic fermentation in Saccharomyces cerevisiae. Bioresour. Technol. 2012, 117, 242–250. [Google Scholar] [CrossRef]
  53. Tusset, A.M.; Balthazar, J.M.; Rocha, R.T.; Ribeiro, M.A.; Lenz, W.B. On suppression of chaotic motion of a nonlinear MEMS oscillator. Nonlinear Dyn. 2019, 99, 537–557. [Google Scholar] [CrossRef]
  54. Petras, I. Fractional-Order Nonlinear Systems: Modeling, Analysis and Simulation; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  55. Ribeiro, M.A.; Balthazar, J.M.; Lenz, W.B.; Rocha, R.T.; Tusset, A.M. Numerical Exploratory Analysis of Dynamics and Control of an Atomic Force Microscopy in Tapping Mode with Fractional Order. Shock Vib. 2020, 2020, 4048307. [Google Scholar] [CrossRef]
  56. Tusset, A.M.; Ribeiro, M.A.; Lenz, W.B.; Rocha, R.T.; Balthazar, J.M. Time Delayed Feedback Control Applied in an Atomic Force Microscopy (AFM) Model in Fractional-Order. J. Vib. Eng. Technol. 2020, 8, 327–335. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the bioreactor.
Figure 1. Schematic representation of the bioreactor.
Symmetry 14 01609 g001
Figure 2. Output concentration and temperature versus time curves (hours) without control.
Figure 2. Output concentration and temperature versus time curves (hours) without control.
Symmetry 14 01609 g002
Figure 3. Bioreactor schematic representation with temperature control in the reactor.
Figure 3. Bioreactor schematic representation with temperature control in the reactor.
Symmetry 14 01609 g003
Figure 4. Output concentration and temperature versus time curves (hours) with control.
Figure 4. Output concentration and temperature versus time curves (hours) with control.
Symmetry 14 01609 g004
Figure 5. Variation of temperature in the cooling jacket and coolant flow rate.
Figure 5. Variation of temperature in the cooling jacket and coolant flow rate.
Symmetry 14 01609 g005
Figure 6. State stabilization point, considering the systems without control and in fractional order representation, considering t = 5000 h, and ( 0.3 q 5 , 6 1 ).
Figure 6. State stabilization point, considering the systems without control and in fractional order representation, considering t = 5000 h, and ( 0.3 q 5 , 6 1 ).
Symmetry 14 01609 g006
Figure 7. State stabilization point, considering the systems with control and in fractional order representation, considering t = 5000 h, and ( 0.3 q 5 , 6 1 ).
Figure 7. State stabilization point, considering the systems with control and in fractional order representation, considering t = 5000 h, and ( 0.3 q 5 , 6 1 ).
Symmetry 14 01609 g007
Figure 8. Temperature and coolant flow rate for controlled systems in fractional order.
Figure 8. Temperature and coolant flow rate for controlled systems in fractional order.
Symmetry 14 01609 g008
Figure 9. State stabilization point, considering the systems with control and in fractional order representation, with perturbation in inlet temperature, considering t = 5000 h, and ( 0.3 q 5 , 6 1 ).
Figure 9. State stabilization point, considering the systems with control and in fractional order representation, with perturbation in inlet temperature, considering t = 5000 h, and ( 0.3 q 5 , 6 1 ).
Symmetry 14 01609 g009
Table 1. Comparison of the variable values in equilibrium, with and without temperature control.
Table 1. Comparison of the variable values in equilibrium, with and without temperature control.
C x C S C p C O 2 T r T a g
Without control1.44721.5215.70.348436.0130.01
With control1.08926.0213.740.437132.0124.18
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tusset, A.M.; Inacio, D.; Fuziki, M.E.K.; Costa, P.M.L.Z.; Lenzi, G.G. Dynamic Analysis and Control for a Bioreactor in Fractional Order. Symmetry 2022, 14, 1609. https://doi.org/10.3390/sym14081609

AMA Style

Tusset AM, Inacio D, Fuziki MEK, Costa PMLZ, Lenzi GG. Dynamic Analysis and Control for a Bioreactor in Fractional Order. Symmetry. 2022; 14(8):1609. https://doi.org/10.3390/sym14081609

Chicago/Turabian Style

Tusset, Angelo M., Danilo Inacio, Maria E. K. Fuziki, Priscilla M. L. Z. Costa, and Giane G. Lenzi. 2022. "Dynamic Analysis and Control for a Bioreactor in Fractional Order" Symmetry 14, no. 8: 1609. https://doi.org/10.3390/sym14081609

APA Style

Tusset, A. M., Inacio, D., Fuziki, M. E. K., Costa, P. M. L. Z., & Lenzi, G. G. (2022). Dynamic Analysis and Control for a Bioreactor in Fractional Order. Symmetry, 14(8), 1609. https://doi.org/10.3390/sym14081609

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

Article Metrics

Back to TopTop