Past exam of the topics in applied mathematics course of the University of Oxford / 2008 / A / 3

An unbiased random walk of bacteria leads, on the continuum scale, to diffusion with coefficient $D$, represented by $D{B}_{XX}$. The cell proliferation rate $\nu GB$ is proportional both to the available biological growth factor concentration and to the bacterial density. Consumption of biological growth factor at rate $\mu GB$ gives its negative reaction term.

There is no diffusion term for $G$ because the biological growth factor is attached to the substrate rather than freely transported. There is also no chemotactic drift in this model: the specified motion is random, while the biological growth factor changes cell proliferation rather than directional movement. Both reaction terms vanish if either bacteria or biological growth factor are absent.

Choose a positive reference biological growth factor concentration $G_{ref}$, for example its concentration ahead of the travelling wave. SetwhereThe chain rule then transforms the equations intoEach displayed coefficient is one with the chosen scales. Hence the dimensionless system isThe free choice of $G_{ref}$ expresses the freedom to normalize the as-yet unspecified concentration scale.

Take profiles $b=b(y)$ and $g=g(y)$ with $y=x-ct$, and use a prime for differentiation with respect to $y$. The chain rule gives $b_t=-c{b}^{\prime }$, $b_{xx}=b^{\prime \prime }$ and $g_t=-c{g}^{\prime }$. ThereforeThese are coupled ordinary differential equations for a right-moving travelling wave with speed $c>0$. A first-order formulation, using $v=b^{\prime }$, is $b^{\prime }=v$, $v^{\prime }=-cv-gb$, $g^{\prime }=gb/c$. In particular a positive profile has $g^{\prime }>0$, consistent with depleted substrate behind the travelling front and fresh biological growth factor ahead.

Substitute $gb=c{g}^{\prime }$ into the bacterial profile equation. This giveswhich integrates once toAt the asymptotic constant states the profile derivative tends to zero. The state ahead of the travelling wave, where $b\to 0$ and $g\to 1$, therefore gives $C=c$. The state behind it, where $g\to 0$ and $b\to b_0$, gives $C=c{b}_0$. Since $c>0$,This first integral expresses the matching between biomass produced and biological growth factor consumed after the chosen nondimensionalization.

At the leading edge $b\to 0$ and $g\to 1$. The linearized bacterial equation iswith characteristic exponentsFor $0<c<2$ they form a nonreal complex conjugate pair. A nonzero profile approaching the travelling front then has an exponentially decaying oscillatory leading term, and repeatedly changes sign. This is incompatible with the assumed positive bacterial density. The leading edge must therefore have real characteristic exponents, requiringAt $c=2$ the repeated characteristic exponent is $-1$ and polynomial factors in the asymptotic profile may occur. The argument supplies a necessary lower bound; it is not an independent existence proof for every speed satisfying it.

The first integral and $b_0=1$ yield the linear identity $b^{\prime }+cb=c(1-g)$. Multiply by the integrating factor $e^{cy}$ and integrate from $-\infty$. Since the bounded limit $b(-\infty )=1$ makes $e^{cy}b(y)\to 0$ there,Positivity of $g$ makes the last integral nonnegative, and strictly positive at any finite point of a strictly positive travelling wave. Thuswith equality approached only at the far-behind boundary. This bound uses the first integral rather than an unsupported assumption that the bacterial profile is monotone.

For the modified strain the travelling wave equations areIntegrating the second across the entire travelling wave gives ${\int }_{\mathbb{R}}gbdy=c$, since the biological growth factor profile still goes from $0$ to $1$. Integrating the first, using the constant asymptotic states and vanishing derivative at infinity, givesConsequentlyThe subtracted integral is strictly positive because the profiles are positive and $\eta >0$. It is finite: $b$ is bounded by continuity and its finite endpoint limits, and $\int gb=c$. The additional density-dependent inhibition therefore reduces the bacterial density left behind; positivity of $b_1$ is part of the assumed positive-wave solution.