---
title: Fractional-Order Chemostat Model (FOCM)
url: https://www.emergentmind.com/topics/fractional-order-chemostat-model-focm
type: topic
---

# Fractional-Order Chemostat Model (FOCM)

The Fractional-Order Chemostat Model (FOCM) is a chemostat formulation in which substrate and biomass dynamics are governed by fractional differential operators rather than ordinary derivatives, with the explicit purpose of representing memory effects, periodic operation, and time-dependent microbial adaptation in biological water treatment. In the recent control-oriented formulation, the model uses the Caputo fractional derivative with sliding memory (CFDS), periodic boundary conditions (PBCs), and Contois growth kinetics to describe pollutant removal under periodic dilution-rate modulation. Within that framework, the primary optimization target is the minimization of the average substrate concentration over one cycle, and the reported numerical results include reductions of average substrate concentration by up to \(40\%\) relative to steady-state operation [2508.01022]. A closely related mathematical framework establishes existence, uniqueness, positivity, boundedness, and stability of periodic solutions for the sliding-memory fractional chemostat itself, thereby supplying the well-posedness theory on which the control problem rests [2506.04420].

## 1. Fractional formulation and periodic-memory structure

The state variables are the substrate concentration \(s(t)\) or \(S(t)\), the biomass concentration \(x(t)\) or \(X(t)\), and, in the control setting, the dilution rate \(D(t)\), with bounds \(D_{\min}\le D(t)\le D_{\max}\). The inflow substrate concentration is \(s_{in}\) or \(S_{in}\), the yield coefficient is \(Y>0\), the fractional order satisfies \(\alpha\in(0,1)\), and the sliding memory length satisfies \(L>0\). The control-oriented model also introduces a dynamic scaling parameter \(\vartheta>0\), appearing as \(\vartheta^{1-\alpha}\), which modulates the amplitude of the fractional dynamics. The specific growth law is the Contois model,
\[
\mu(s,x)=\frac{\mu_{\max}s}{Kx+s},
\]
with \(\mu_{\max}>0\) and \(K>0\) [2508.01022].

The CFDS replaces the fixed-origin history integral of the classical Caputo derivative by a finite sliding window. The left-sided operator is
\[
{}^{\mathrm{MC}_{L}}D_t^\alpha f(t)=\frac{1}{\Gamma(1-\alpha)}\int_{t-L}^{t}(t-\tau)^{-\alpha}f'(\tau)\,d\tau,
\]
and the right-sided operator used in the adjoint system is
\[
{}^{\mathrm{MC}_{L+}}D_t^\alpha f(t)=-\frac{1}{\Gamma(1-\alpha)}\int_{t}^{t+L}(\tau-t)^{-\alpha}f'(\tau)\,d\tau.
\]
With PBCs, the full fractional-order chemostat system is
\[
{}^{\mathrm{MC}_{L}}D_t^\alpha s(t)=\vartheta^{1-\alpha}\left[-\frac{1}{Y}\mu(s(t),x(t))x(t)+D(t)(s_{in}-s(t))\right],
\]
\[
{}^{\mathrm{MC}_{L}}D_t^\alpha x(t)=\vartheta^{1-\alpha}[\mu(s(t),x(t))-D(t)]x(t),
\]
together with
\[
s(t)=s(t+T),\qquad x(t)=x(t+T),\qquad D(t)=D(t+T).
\]
A key periodicity-preserving identity is
\[
\int_{0}^{T}{}_L^{MC}D_t^\alpha s(t)\,dt=0.
\]
This finite-window construction is central: unlike the standard Caputo derivative integrating from a fixed initial time, CFDS is designed to preserve periodicity of \(T\)-periodic absolutely continuous functions, and the parameters \(\alpha\), \(L\), and \(\vartheta\) jointly shape memory effects, responsiveness, and the magnitude of the dynamics [2508.01022; 2506.04420].

## 2. Relation to the classical chemostat and reduction to a one-dimensional FDE

The baseline integer-order chemostat model is written as
\[
\frac{dX}{dt}=(\mu(S)-u)X,\qquad 
\frac{dS}{dt}=u(S_{in}-S)-\frac{1}{Y}\mu(S)X.
\]
The fractional formulation modifies these balances by replacing \(d/dt\) with a sliding-memory fractional operator and, in the control formulation, multiplying the right-hand side by \(\vartheta^{1-\alpha}\). This introduces non-local dependence into the growth and washout terms and changes the transient regime from memoryless ODE dynamics to hereditary FDE dynamics [2508.01022].

A defining structural result is the reduction of the two-dimensional system to a one-dimensional fractional differential equation. The transformation
\[
z(t)=Y(s_{in}-s(t))-x(t)
\]
yields
\[
{}_L^{\mathrm{MC}}D_t^\alpha z(t)=-\vartheta^{1-\alpha}D(t)z(t)
\]
in the control setting, or \( {}^{\mathrm{MC}_{L}}D_t^\alpha z(t)=-D(t)z(t)\) in the foundational formulation. Under the periodicity assumptions and the energy argument given in the analysis, there are no non-trivial periodic solutions for \(z\), hence
\[
x(t)=Y(s_{in}-s(t)).
\]
The reduced substrate equation becomes
\[
{}_L^{\mathrm{MC}}D_t^\alpha s(t)=\vartheta^{1-\alpha}[D(t)-\nu(s(t))](s_{in}-s(t)),
\]
with
\[
\nu(s)=\mu\big(s,Y(s_{in}-s)\big)=\frac{\mu_{\max}s}{KY(s_{in}-s)+s}.
\]
In the corresponding non-scaled framework the same reduction gives
\[
{}^{\mathrm{MC}_{L}}D_t^\alpha S(t)=\big(D(t)-\nu(S(t))\big)\big(S_{in}-S(t)\big).
\]
This reduction is computationally decisive because it preserves the essential substrate dynamics needed for optimization while eliminating the second state through the exact identity \(X(t)=Y(S_{in}-S(t))\) [2508.01022; 2506.04420].

For constant dilution \(D<\mu_{\max}\), the reduced equation has a unique interior equilibrium satisfying \(D=\nu(\bar S)\), namely
\[
\bar S=\frac{D K Y S_{in}}{\mu_{\max}+D(KY-1)},\qquad 
\bar X=Y(S_{in}-\bar S).
\]
The alternative equilibrium is the washout state \((S_{in},0)\) [2506.04420].

## 3. Existence, uniqueness, positivity, and stability

The periodic theory is developed in the space \(AC_T\) of absolutely continuous \(T\)-periodic functions. A standard feasible state set is
\[
X=\{\,S\in AC_T:\ 0\le S(t)\le S_{in}\,\},
\]
and admissible controls are \(T\)-periodic \(L^\infty\) functions satisfying the average constraint and pointwise bounds. In the foundational analysis, the right-hand side satisfies Carathéodory regularity, and the reduced problem is rewritten as a Volterra integral equation over the sliding window. The induced fixed-point operator is continuous, maps the compact convex set into itself, and has relatively compact image, so Schauder’s fixed-point theorem yields existence of at least one \(T\)-periodic Carathéodory solution [2506.04420].

The control-oriented analysis proves that, under \(s(0)\in(0,s_{in})\), \(x(0)>0\), and \(D(t)<\mu_{\max}\) for all \(t\ge 0\), the reduced fractional optimal periodic control problem has at least one optimal periodic solution \((s^*,D^*)\). The argument uses compactness and convexity of the feasible state set, weak-\(^*\) compactness of the admissible control set via Banach–Alaoglu, continuity of the control-to-state map through the Volterra formulation, and continuity of the objective functional \(J(D)\) [2508.01022].

Positivity and boundedness are built into the reduced dynamics. At \(S=0\), the drift satisfies \(f(t,0)=D(t)S_{in}>0\), and at \(S=S_{in}\), \(f(t,S_{in})=0\), which enforces \(0\le S(t)\le S_{in}\). Because \(X(t)=Y(S_{in}-S(t))\), biomass remains nonnegative, and if \(S(0)<S_{in}\) then \(X(t)>0\) for all \(t>0\). The biologically feasible invariant set is
\[
\Omega=\{(S,X)\in\mathbb{R}^2:\ 0\le S\le S_{in},\ 0\le X\le YS_{in}\}.
\]
For stability, the reduced one-dimensional dynamics \(g(S)=(D-\nu(S))(S_{in}-S)\) satisfy
\[
g'(\bar S)=-\nu'(\bar S)(S_{in}-\bar S)<0,
\]
so linearization gives algebraic decay governed by the Mittag–Leffler function. In the full linear fractional system, asymptotic stability requires eigenvalues \(\lambda\) of the Jacobian to satisfy
\[
|\arg(\lambda)|>\frac{\alpha\pi}{2}.
\]
At washout, stability holds for \(D>\mu_{\max}\) and fails for \(D<\mu_{\max}\) [2506.04420].

Uniqueness is conditional rather than universal. In the foundational paper, uniqueness of the \(T\)-periodic Carathéodory solution is established under any of three stated conditions, including \(L\ge T\), \(KY>1\), and \(D(t)>\mu_{\max}\) for all \(t\), or conditions involving
\[
s^*=S_{in}\frac{\sqrt{KY}}{\sqrt{KY}+1}
\]
and upper bounds on \(D(t)\) or \(\bar D\) by \(\nu(s^*)\). In the control problem, uniqueness of the optimal periodic control is proved under \(KY\ne 1\) and additional conditions such as
\[
s(0)\le s=\frac{\sqrt{KY}\,s_{in}}{\sqrt{KY}+1},
\]
together with either \(D(t)\le \nu(s)\) for all \(t\) or \(\bar D\le \nu(s)\) [2506.04420; 2508.01022].

## 4. Optimal periodic control and bang-bang structure

The optimal periodic control problem minimizes the cycle-averaged substrate concentration,
\[
\min_D J(D),\qquad J(D)=\frac{1}{T}\int_{0}^{T}s(t)\,dt,
\]
subject to
\[
\frac{1}{T}\int_{0}^{T}D(t)\,dt=\bar D,
\]
the state and control bounds
\[
0\le s(t)\le s_{in},\qquad x(t)>0,\qquad D_{\min}\le D(t)\le D_{\max},
\]
and the periodicity conditions on \(s\), \(x\), and \(D\). The model allows non-constant optimal controls, and the paper states that if \(KY\ne 1\) and \(\alpha\in(0,1)\), then non-constant optimal periodic solutions may exist with improved average substrate concentration compared to steady-state operation. The perturbation analysis distinguishes three regimes: for \(KY<1\), concavity of \(\nu\) implies \([\nu(\cdot)]_{av}<\bar D\), enabling \(s^*<\bar s\); for \(KY>1\), improvements are possible but not guaranteed; for \(KY=1\), averages match steady-state [2508.01022].

The fractional Pontryagin maximum principle yields the Hamiltonian
\[
H(s,p,D)=\frac{1}{T}s(t)+p(t)\,\vartheta^{1-\alpha}[D(t)-\nu(s(t))](s_{in}-s(t)),
\]
the state equation
\[
{}_L^{\mathrm{MC}}D_t^\alpha s(t)=\vartheta^{1-\alpha}[D(t)-\nu(s(t))](s_{in}-s(t)),
\]
and the right-sided adjoint equation
\[
{}^{\mathrm{MC}_{L+}}D_t^\alpha p(t)=-\frac{1}{T}+p(t)\,\vartheta^{1-\alpha}\left[\nu'(s(t))(s_{in}-s(t))+D(t)-\nu(s(t))\right],
\]
with periodic transversality \(p(0)=p(T)\). The switching function is
\[
\phi(t)=p(t)\,\vartheta^{1-\alpha}(s_{in}-s(t)).
\]
Since \(s_{in}-s(t)>0\) and \(\vartheta^{1-\alpha}>0\), the optimal control satisfies the bang-bang law
\[
D^*(t)=
\begin{cases}
D_{\max}, & \phi(t)<0,\\
D_{\min}, & \phi(t)>0.
\end{cases}
\]
Singular arcs are ruled out because \(\phi(t)\equiv 0\) would imply \(p(t)=0\), contradicting the adjoint equation \( {}^{\mathrm{MC}_{L+}}D_t^\alpha p=-1/T\). Accordingly, the optimal control is bang-bang, and the number of switches per period is even due to periodicity [2508.01022].

## 5. Numerical realization and reported performance

The reported implementation uses a Fourier–Gegenbauer pseudospectral (FG-PS) method on \(N\) equispaced collocation points. The CFDS term is approximated through precomputed integration matrices tailored to periodic problems, and the discretized reduced control problem is solved as a constrained NLP by MATLAB’s `fmincon` with the `sqp` algorithm. Because the pseudospectral approximation of a bang-bang profile exhibits Gibbs oscillations, the workflow includes an edge-detection correction: Fourier coefficients of the predicted \(D(t)\) are computed, discontinuities are detected on the spectral profile, the bang-bang control is reconstructed from the detected switches, and the state \(s(t)\) is corrected with `fsolve` using the reconstructed \(D(t)\) in a predictor-corrector approach [2508.01022].

Typical parameters used in the examples are \(N=300\)–\(400\) collocation points, \(M=400\)–\(500\) interpolation points, \(\alpha=0.85\), \(L=5\), \(\vartheta=0.25\ \mathrm{h}\), \(T=15\ \mathrm{h}\), \(D_{\min}=0.02\ \mathrm{h}^{-1}\), \(D_{\max}=1.95\ \mathrm{h}^{-1}\), \(\mu_{\max}=2\ \mathrm{h}^{-1}\), \(K=5\), \(Y=1\), \(s_{in}=8\ \mathrm{mg/L}\), and \(\bar D=0.5\ \mathrm{h}^{-1}\). In the base case, the average substrate concentration is reduced from the steady-state value \(\bar s=5\ \mathrm{mg/L}\) to \(s^*\approx 3.622\ \mathrm{mg/L}\), corresponding to an approximately \(27.56\%\) improvement. With \(\vartheta=32\), the paper reports \(s^*\approx 3.001\ \mathrm{mg/L}\) versus steady-state \(5\ \mathrm{mg/L}\), a \(39.98\%\) reduction, described as nearly \(40\%\) [2508.01022].

The sensitivity analysis attributes distinct operational roles to \(\alpha\), \(L\), \(\vartheta\), and \(T\). Lower \(\alpha\) produces stronger memory and more switches per period; higher \(\alpha\) yields fewer switches, while preserving the bang-bang structure. The function \(s^*(\alpha)\) is non-monotonic for small \(\alpha\), with a peak around \(\alpha=0.3\), and decreases monotonically for \(\alpha\ge 0.4\). As \(L\) increases, \(s^*\) first rises slightly from \(L=0.5\) to \(1.5\), then decreases up to approximately \(10\), with diminishing returns beyond that range. Increasing \(\vartheta\) monotonically decreases \(s^*\), and increasing the period \(T\) over the tested range \(1\)–\(20\ \mathrm{h}\) also decreases \(s^*\). Because the optimal control is discontinuous, convergence of the state approximation is algebraic rather than spectral, although the switching times are reported as accurate across resolutions to machine precision [2508.01022].

## 6. Interpretation, limitations, and adjacent research directions

The FOCM interprets \(\alpha\) and \(L\) as descriptors of microbial adaptation latency and non-local response, while \(\vartheta\) tunes dynamic responsiveness. The periodic boundary conditions and the periodicity-preserving structure of CFDS align the mathematics with cyclic operational strategies such as scheduled dilution or diurnal loading. The foundational analysis states the principal assumptions explicitly: perfect mixing, constant \(S_{in}\), stationary parameters \((\mu_{\max},K,Y)\), Contois kinetics, periodic operation, and a finite sliding memory window representing localized historical dependence [2508.01022; 2506.04420].

The principal limitations are equally explicit. The model is single-species and single-substrate; parameters are assumed constant and identifiable; fractional operators and non-smooth controls introduce computational overhead; Gibbs phenomena require numerical mitigation; and uniqueness is only guaranteed under specific conditions, so multiple local optima may arise in the general case. The stated practical guidelines are correspondingly operational: choose \(\alpha\) close to \(1\) for faster-responding communities and fewer control switches; use a moderately large \(L\) to capture useful memory without excessive computational burden; and increase \(\vartheta\) to enhance responsiveness and reduce \(s^*\), while balancing actuator limits and system stability. Proposed extensions include multi-species and multi-substrate models, inhibition kinetics such as Haldane kinetics, stochastic forcing through uncertain influent \(s_{in}(t)\) or random shocks, adaptive or robust fractional control, and distributed systems based on fractional PDEs [2508.01022; 2506.04420].

A recurrent misconception in the broader fractional-compartment literature is that fractional modeling is exhausted by replacing an integer-order derivative with a fractional one. A related stochastic study argues instead for a physically grounded construction in which fractional behavior arises from semi-Markov dynamics and Mittag–Leffler waiting times, yielding Riemann–Liouville operators at the mean-field level and exact SSA realizations for small populations. That work also emphasizes that deterministic fractional equations can diverge substantially from stochastic dynamics when the population is small. This suggests that future FOCM research may need to distinguish more sharply between phenomenological sliding-memory formulations for periodic bioprocess control and stochastic semi-Markov derivations for chemostat systems operating far from the continuum limit [2312.05268].

Source: https://www.emergentmind.com/topics/fractional-order-chemostat-model-focm