\documentclass{beamer}
%\usepackage[latin1]{inputenc}
\usepackage{graphicx}
\usepackage{hyperref}
\usepackage{amsmath,color}
\usepackage{amsthm}
\usepackage{amstext}
\usepackage{amsopn}
\usepackage{pst-plot,pst-node,pst-slpe}
\usepackage{hyperref}
\usepackage{epsfig}
\usepackage{subfigure}

\usetheme{Frankfurt}
\usecolortheme{seahorse}
%\useoutertheme{infolines}


\usepackage{epstopdf}
\DeclareGraphicsExtensions{.pdf,.eps,.png,.jpg,.mps}

\newcommand{\ds}{\displaystyle}
\newcommand{\g}{\gamma}
\newcommand{\ep}{\varepsilon}
\newcommand{\la}{\lambda}
\newcommand{\R}{{\mathbb R}}
\newcommand{\Om}{\Omega}
\newcommand{\noi}{\noindent}
\newcommand{\f}{\frac}


\title[Chemostat]{Further studies of a reaction-diffusion system for an unstirred chemostat with internal storage}
\author[Junping Shi]{Junping Shi}


\institute{
College of William and Mary, USA\\

\vspace{0.2in}

{\blue Joint work with Sze-Bi Hsu and Feng-Bin Wang}}


\date{{\green
}\\
{\blue Mathbiology Seminar}\\
{\red May 3, 2013} }

\begin{document}

\begin{frame}
\titlepage
\begin{center}
\includegraphics[height=0.5in]{wm1.eps}\hspace{0.3in}
\end{center}
\end{frame}




\section{Chemostat}

\begin{frame}\frametitle{Chemostat}
{\scriptsize

{\cyan Chemostat}: an apparatus for growing bacterial cultures at a constant rate by controlling the supply of nutrient medium

\begin{center}
\includegraphics[width=0.4\textwidth]{chemostats.eps}
\hspace{0.05in}
\includegraphics[width=0.4\textwidth]{chemostats2.eps}
\end{center}

}
\end{frame}

\begin{frame}\frametitle{Minimal chemostat model}
{\scriptsize

\vspace{0.1in}
{\magenta Simplest setup}: a vessel with volume $V$ with an inflow $F$ of the culture medium and nutrient (with a constant concentration $S^0$), and an outflow $F$ of the mixture of medium, nutrient and bacteria. Since the inflow and outflow are same, then the volume of the culture remains a constant. We assume that the culture is {\red well-stirred}.

\vspace{0.1in}\pause
{\cyan Variables}: $S(t)$ (concentration of nutrient), $N(t)$ (concentration of bacteria)

\vspace{0.1in}
{\magenta Model}:

\begin{equation*}
    \begin{cases}
    & \ds S_t= d(S^0-S)-\frac{1}{\gamma} \frac{mS}{a+S}N, \\
    & \ds N_t= \frac{mS}{a+S}N-dN.
    \end{cases}
\end{equation*}

$d=F/V$: dilution rate of the chemostat, \; $\gamma$: the growth yield or the conversion rate\\
$\ds  \frac{mS}{a+S}$: Monod growth function, similar to  Michaelis-Menten  enzyme kinetics function

\vspace{0.1in}\pause
{\cyan Convergence}: $\ds \left(S+\gamma^{-1}N\right)_t=d[S^0-(S+\gamma^{-1}N)]$ so $\ds \lim_{t\to\infty}[S(t)+\gamma^{-1}N(t)]=S^0$.

\vspace{0.1in}
{\magenta Limiting system}: $\ds N_t=\frac{m(\gamma S^0-N)N}{\gamma a+\gamma S^0-N}$, and $\ds\lim_{t\to\infty}N(t)=\gamma S^0$.

}
\end{frame}

\begin{frame}\frametitle{Competition model and competition exclusion}
{\scriptsize

Suppose there are two kinds of bacteria. Then

$ \ds S_t= d(S^0-S)-\frac{1}{\gamma_1} \frac{m_1S}{a_1+S}N_1-\frac{1}{\gamma_2} \frac{m_2S}{a_2+S}N_2, \;\;
    \ds (N_i)_t= \frac{m_iS}{a_i+S}N_i-dN_i.$

Rescaling: $N_1 \mapsto N_1/\gamma_1$, and $N_2 \mapsto N_2/\gamma_2$
\begin{equation*}
    \begin{cases}
    & \ds S_t= d(S^0-S)- \frac{m_1S}{a_1+S}N_1- \frac{m_2S}{a_2+S}N_2, \\
    & \ds (N_i)_t= \frac{m_iS}{a_i+S}N_i-dN_i,
    %& \ds (N_2)_t= \frac{m_2S}{a_2+S}N_2-dN_2.
    \end{cases}
\end{equation*}
{\cyan Conservation}: $\ds \lim_{t\to\infty} [S(t)+N_1(t)+N_2(t)]=S^0$

{\magenta Limiting system}:
\begin{equation*}
    \begin{cases}
    & \ds (N_1)_t= \frac{m_1(S^0-N_1-N_2)}{a_1+S^0-N_1-N_2}N_1-dN_1, \\
    & \ds (N_2)_t= \frac{m_2(S^0-N_1-N_2)}{a_2+S^0-N_1-N_2}N_2-dN_2.
    \end{cases}
\end{equation*}
{\blue [Hsu-Hubble-Waltman, 1978, SIAM-AM], [Hsu, 1978, SIAM-AM]} ($n$-species)\\
Let $b_i=m_i/d$ and $\la_i=a_i/(b_i-1)$ for $i=1,2$. \\
(i) If $b_1\le 1$ or $\la_i>S^0$, then $\ds\lim_{t\to\infty} N_i(t)=0$.\\
(ii) If $\la_1<\la_2$ and $\la_1<S^0$, then
$\ds\lim_{t\to\infty} N_1(t)=\frac{d(S^0-\la_1)(a_1+\la_1)}{\la_1}$ and $\ds\lim_{t\to\infty} N_2(t)=0$.


}
\end{frame}

\section{Diffusion}


\begin{frame}\frametitle{Diffusion: unstirred model}
{\scriptsize

{\blue [Hsu-Waltman, 1993, SIAM-AM]}
\begin{equation} \label{1.1}
 \begin{cases}
 \ds S_{t}= dS_{xx}-m_1 f_1\left(S\right)u-m_2 f_2\left(S\right)v,&\;\; x\in (0,1), \\
\ds u_{t}=du_{xx}+m_1 f_1\left(S\right)u, &\;\; x\in (0,1), \\
\ds v_{t}=dv_{xx}+m_2 f_2\left(S\right)v, &\;\; x\in (0,1),\\
S_x(0,t)=-S^{0}, \; S_{x}(1,t)+ \gamma S(1,t)=0, &\\
 w_{x}(0,t)=0, \; w_{x}(1,t)+ \gamma w(1,t)=0, &\;\; w=u,v,
 \end{cases}
\end{equation}
where
\begin{equation}\label{1.2}
    f_i(S)=\frac{S}{a_i+S}, \;\;
\end{equation}

Limiting system:
\begin{equation}\label{2.1}
\begin{cases}
 \ds u_{t}=du_{xx}+ m_1 f_1\left(\phi-u-v\right)u, & \;\; x\in (0,1),\; t>0,\\
 \ds  v_{t}=dv_{xx}+ m_2 f_2\left(\phi-u-v\right)v, & \;\; x\in (0,1),\; t>0,\\
  w_{x}(0,t)=0, \; w_{x}(1,t)+ \gamma w(1,t)=0, &\;\; w=u,v,
  \end{cases}
\end{equation}
where $\ds \phi(x)=S^{0}\left(\frac{1+\gamma}{\gamma}-x\right)$.



}
\end{frame}




\begin{frame}\frametitle{Coexistence}
{\scriptsize

 {\blue [Hsu-Waltman, 1993, SIAM-AM] [Hsu-Smith-Waltman, 1995, TAMS]}

 \vspace{0.1in}

 When $m_1> d\la_1(f_1(\phi))$ and $m_2> d\la_1(f_2(\phi))$, there exist two semi-trivial steady state solutions $(u_*,0)=(w_*(\cdot; m_1,a_1),0)$, and $(0,v_*)=(0,w_*(\cdot; m_2,a_2))$. For the coexistence steady states:

 \begin{enumerate}
   \item If $m_1> d\la_1(f_1(\phi-w_*(m_2,a_2)))$ and $m_2> d\la_1(f_2(\phi-w_*(m_1,a_1)))$, then both of $(u_*,0)$ and $(0,v_*)$ are unstable, and there exists at least one positive coexistence steady state solution which is stable.
   \item If $m_1< d\la_1(f_1(\phi-w_*(m_2,a_2)))$ and $m_2< d\la_1(f_2(\phi-w_*(m_1,a_1)))$, then both of $(u_*,0)$ and $(0,v_*)$ are locally asymptotically stable, and there exists at least one positive coexistence steady state solution.
 \end{enumerate}

If  $a_1=a_2$, a coexistence steady state solution only exists for a unique  $m_1=m_2=d\la_1(f_1(\phi-w_*(m_1,a_2)))$. A competition exclusion occurs when $m_1<m_2$ or $m_1>m_2$.

 \vspace{0.1in}
 {\red Question}: when $a_1\ne a_2$, how do you characterize the set of coexistence parameter ranges\\
 $\{(m_1,m_2): m_1> d\la_1(f_1(\phi-w_*(m_2,a_2))), \;m_2> d\la_1(f_2(\phi-w_*(m_1,a_1)))$ and \\
 $\{(m_1,m_2): m_1< d\la_1(f_1(\phi-w_*(m_2,a_2))), \;m_2< d\la_1(f_2(\phi-w_*(m_1,a_1)))$.

}
\end{frame}

\begin{frame}\frametitle{Numerical simulations}
{\scriptsize

\begin{center}
\includegraphics[width=0.4\textwidth]{chemo-1.eps}
\hspace{0.05in}
\includegraphics[width=0.4\textwidth]{chemo-2.eps}
\end{center}

Here $m_1=m_2=1.5$, $\gamma=1$, $d=0.05$, $a_2=1$.\\
Left (competition exclusion): $a_1=0.8$; Right (coexistence): $a_1=0.98$.

}
\end{frame}

\section{Internal storage}

\begin{frame}\frametitle{Droop model}
{\scriptsize

For the variable yield case, {\blue [Droop, 1973, J. Phycology]} formulated the following model with an internal storage of the nutrient:
\begin{equation}\label{eqn1.2}
\begin{cases}
  \ds S'=(S^{(0)}-S)D-f_1(S,Q_1)u-f_2(S,Q_2)v, \;\;&\\
 \ds u'=(\mu_1(Q_1)-D)u, \;\;&\\
  \ds Q_1'=f_1(S,Q_1)-\mu_1(Q_1)Q_1, \;\;&\\
  \ds v'=(\mu_2(Q_2)-D)v, \;\;&\\
  \ds Q_2'=f_2(S,Q_2)-\mu_2(Q_2)Q_2, \;\;&\\
  S(0)\ge 0,\; u(0)\ge 0,\; v(0)\ge 0,\; Q_1(0)\ge Q_{\min,1},\; Q_2(0) \ge Q_{\min,2}.
  \end{cases}
\end{equation}
Here $Q_i(t)$ ($i=1,2$) is the average amount of stored nutrient per cell of the $i$-th population.
$\mu_{i}(Q_{i})$
is the growth rate of the $i$-th population as a function of cell quota
$Q_{i}$, $f_{i}(S,Q_{i})$ is the per capital nutrient uptake rate per cell of the $i$-th population as a function of nutrient concentration $S$ and cell quota $Q_{i}$, and $Q_{\min,i}$ is the threshold cell
quota below which no growth of the $i$-th population occurs. Typically the functions $\mu_i(Q_i)$ and $f_{i}(S,Q_{i})$ satisfy  the following assumptions: ($i=1,2$)
\begin{enumerate}
  \item[$(H1)$] $\mu_i(Q_i)$ is continuously differentiable for $Q_i\ge Q_{\min,i}$, $\mu_i(Q_{\min,i})=0$,  $\mu_i(Q_i)>0$ and $\mu_i'(Q_i)>0$ for $Q\ge Q_{\min,i}$.
  \item[$(H2)$] $f_i(S,Q_i)$ is continuously differentiable for $S>0$ and $Q\ge Q_{\min,i}$, $f_i(0,Q_i)=0$, $f_i(S,Q_i)\ge 0$, $\ds \frac{\partial f_i}{\partial S}(S,Q_i)>0$ and  $\ds \frac{\partial f_i}{\partial Q_i}(S,Q_i)\leq 0$ for $S>0$ and $Q_i\ge Q_{\min,i}$.
\end{enumerate}

}
\end{frame}

\begin{frame}\frametitle{Functions}
{\scriptsize

Some specific functions $\mu_i(Q_i)$ and $f_{i}(S,Q_{i})$  are
\begin{equation}\label{eqn1.3}
\begin{split}
  \ds \mu_i(Q_i)&=\mu_{i\infty}\left(1-\frac{Q_{\min,i}}{Q_i}\right),\\
  \ds \mu_i(Q_i)&=\mu_{i\infty} \frac{Q_i-Q_{\min,i}}{K_i+Q_i-Q_{\min,i}},
\end{split}
\end{equation}
where $\mu_{i\infty}$ is the maximal growth rate of the $i$-th population,
and
\begin{equation}\label{eqn1.4}
\begin{split}
  f_i(S,Q_i)&=\ds\rho_i(Q_i)\frac{S}{k_i+S},\ \mbox{where}\  \rho_i(Q_i)\  \mbox{is defined by}\\
  \rho_i(Q_i)&=\ds\rho_{\max,i}^{\rm high}-(\rho_{\max,i}^{\rm high}-\rho_{\max,i}^{\rm low})
  \frac{Q_i-Q_{\min,i}}{Q_{\max,i}-Q_{\min,i}}, \;\; \mbox{or}\\
  \rho_i(Q_i)&=\ds\rho_{\max,i}\frac{Q_{\max,i}-Q_{i}}{Q_{\max,i}-Q_{\min,i}},\;\; \mbox{or}\\
  \rho_i(Q_i)&=\rho_i^c, \;\; \mbox{(constant)},
\end{split}
\end{equation}
where $Q_{\min,i}\leq Q_i\leq Q_{\max,i}$ (a maximal possible quota).

}
\end{frame}

\begin{frame}\frametitle{Michael Droop (1918-2011)}
{\scriptsize

Michael Droop is a well-known British marine biologist who is famous for his Cell-Quota model of phytoplankton growth.

\vspace{0.1in}

{\blue [1968]} M. Droop,  Vitamin B12 and marine ecology. IV. The kinetics of uptake, growth and inhibition in Monochrysis lutheri. J. Mar. Biol. Assoc. UK 48 (1968) 689--733. (cited 542 times on Google scholar)

{\blue [1973]} M. Droop, Some thoughts on nutrient limitation in algae, J. Phycol., 9 (1973), 264--272. (cited 523 times on Google scholar)



\begin{center}
\includegraphics[width=0.6\textwidth]{droop.eps}
\end{center}


}
\end{frame}

\begin{frame}\frametitle{Result for the ODE Droop model}
{\scriptsize

{\blue [Smith-Waltman, 1994, SIAM-AM]}

1. Trivial steady state solution $E_0=(\tilde{S}, 0, \tilde{Q}_1, 0, \tilde{Q}_2)$.\\
2. Semi-trivial steady state solution $E_1=(S^1, u^1, Q^1_1, 0, Q^1_2)$.\\
3. Semi-trivial steady state solution $E_2=(S^2, 0, Q^2_1, v^2, Q^2_2)$.\\

\vspace{0.1in}
$E_1$ exists if there exists $Q^1_1>Q_{min,1}$ such that $\mu_1(Q^1_1)=D$ and $DQ^1_1<f_1(S^{(0)}, Q^1_1)$.\\
$E_2$ exists if there exists $Q^2_1>Q_{min,2}$ such that $\mu_2(Q^2_2)=D$ and $DQ^2_2<f_2(S^{(0)}, Q^2_2)$.

\vspace{0.1in}
{\green Theorem}. Assume that the steady state solutions are all non-degenerate.\\
1. If $E_1$ and $E_2$ do not exist, then $E_0$ is globally asymptotically stable.\\
2. If $E_1$ exists but $E_2$ does not exist, then $E_1$ is globally asymptotically stable.\\
3. If $E_2$ exists but $E_1$ does not exist, then $E_2$ is globally asymptotically stable.\\
4. If  $E_1$ and $E_2$ both exist, and $S^1>S^2$, then $E_2$ is globally asymptotically stable.\\
5. If  $E_1$ and $E_2$ both exist, and $S^1<S^2$, then $E_1$ is globally asymptotically stable.\\
6. If  $E_1$ and $E_2$ both exist, and $S^1=S^2$, then there exists a family of steady states connecting $E_1$ and $E_2$.

\vspace{0.1in} {\red No coexistence steady state solutions except the degenerate case.}

}
\end{frame}







\section{Conclusions}
\begin{frame}\frametitle{Conclusions}
{\scriptsize

\begin{enumerate}
  \item The dynamics of the reaction-diffusion competitive system for an unstirred chemostat with internal storage is completed. A sharp threshold diffusion coefficient $d_{0,1}$ exists, and when $d$ is above $d_{0,1}$ both species become extinct. A second  sharp threshold diffusion coefficient $d_{0,2}$ so that species 1 always survives for $d_{0,2}<d<d_{0,1}$. For $0<d<d_{0,2}$, a coexistence steady state is possible.
  \item For one-species model, it is proved that the limiting steady state is larger for smaller $d$. The limiting quota $Q(x)$ appears to be smaller for smaller $d$, and $Q(x)$ also appears to be decreasing in $x$. This indicates that when the bacteria and nutrient mobility is low, bacteria tends to store less nutrient. And the nutrient storage rate is higher when it is closer to the nutrient source point ($x=0$ here). ({\red Prove the monotonicity of $Q_d(x)$ in $d$ and $x$.})
  \item Similar to the unstirred Monod chemostat model, it is hard to determine the parameter ranges where a coexistence steady state is achieved. The coexistence steady state is observed in numerical simulations for various parameters values. It appears that a coexistence steady state is more robust for smaller $d$. Note that a coexistence steady state in general does not exist for well-stirred chemostat models. ({\red Obtain explicit parameter range for the coexistence for both Monod and Droop unstirred chemostat models})
\end{enumerate}


}
\end{frame}

\begin{frame}\frametitle{}
\begin{center}
\includegraphics[height=1.4in]{thanks.eps}
\end{center}
\end{frame}




\end{document}

