Inversion-Based and Optimal Feedforward Control for Population Dynamics with Input Constraints and Self-Competition in Chemostat Reactor Applications

Author(s):  
Anna-Carina Kurth ◽  
Kevin Schmidt ◽  
Oliver Sawodny

Abstract Through chemostat reactors, organisms can be observed under laboratory conditions. Hereby, the reactor contains the biomass, whose growth can be controlled via the dilution rate respectively the speed of a pump. Due to physical limitations, input constraints need to be considered. The population density in the reactor can be described by a hyperbolic nonlinear integro partial differential equation of first order. The steady-states and generalized eigenvalues and -modes of these integro partial differential equation are determined. In order to track a desired reference trajectory an optimal and an inversion-based feedforward control are designed. For the optimal feedforward control, the singular arc of the control is calculated and a switching strategy is stated, which explicitly considers the input constraints. For the inversion-based feedforward control, the integro partial differential equation is first linearized around the steady-state. To comply with the input constraints a control system simulator is designed. For the simulation model, the integro partial differential equation is approximated using Galerkin's method. Simulations show the functionality of the designed controls and provide the basis for comparison. The inversion-based feedforward control operates well near the steady-state, whereas the performance of the optimal feedforward control is not bounded to the proximity to the steady-state.

2021 ◽  
pp. 1-20
Author(s):  
STEPHEN TAYLOR ◽  
XUESHAN YANG

Abstract The functional partial differential equation (FPDE) for cell division, $$ \begin{align*} &\frac{\partial}{\partial t}n(x,t) +\frac{\partial}{\partial x}(g(x,t)n(x,t))\\ &\quad = -(b(x,t)+\mu(x,t))n(x,t)+b(\alpha x,t)\alpha n(\alpha x,t)+b(\beta x,t)\beta n(\beta x,t), \end{align*} $$ is not amenable to analytical solution techniques, despite being closely related to the first-order partial differential equation (PDE) $$ \begin{align*} \frac{\partial}{\partial t}n(x,t) +\frac{\partial}{\partial x}(g(x,t)n(x,t)) = -(b(x,t)+\mu(x,t))n(x,t)+F(x,t), \end{align*} $$ which, with known $F(x,t)$ , can be solved by the method of characteristics. The difficulty is due to the advanced functional terms $n(\alpha x,t)$ and $n(\beta x,t)$ , where $\beta \ge 2 \ge \alpha \ge 1$ , which arise because cells of size x are created when cells of size $\alpha x$ and $\beta x$ divide. The nonnegative function, $n(x,t)$ , denotes the density of cells at time t with respect to cell size x. The functions $g(x,t)$ , $b(x,t)$ and $\mu (x,t)$ are, respectively, the growth rate, splitting rate and death rate of cells of size x. The total number of cells, $\int _{0}^{\infty }n(x,t)\,dx$ , coincides with the $L^1$ norm of n. The goal of this paper is to find estimates in $L^1$ (and, with some restrictions, $L^p$ for $p>1$ ) for a sequence of approximate solutions to the FPDE that are generated by solving the first-order PDE. Our goal is to provide a framework for the analysis and computation of such FPDEs, and we give examples of such computations at the end of the paper.


2015 ◽  
Vol 47 (1) ◽  
pp. 89-94
Author(s):  
C.L. Yu ◽  
D.P. Gao ◽  
S.M. Chai ◽  
Q. Liu ◽  
H. Shi ◽  
...  

Frenkel's liquid-phase sintering mechanism has essential influence on the sintering of materials, however, by which only the initial 10% during isothermal sintering can be well explained. To overcome this shortage, Nikolic et al. introduced a mathematical model of shrinkage vs. sintering time concerning the activated volume evolution. This article compares the model established by Nikolic et al. with that of the Frenkel's liquid-phase sintering mechanism. The model is verified reliable via training the height and diameter data of cordierite glass by Giess et al. and the first-order partial differential equation. It is verified that the higher the temperature, the more quickly the value of the first-order partial differential equation with time and the relative initial effective activated volume to that in the final equibrium state increases to zero, and the more reliable the model is.


Sign in / Sign up

Export Citation Format

Share Document