= Solution
The <bacteria> disperse by <diffusion> with coefficient $D_B$. Their additional <chemotaxis> <flux> is
$$
J_{B,\mathrm{chem}}=\frac{\chi_B B}{(N_0+N)^2}N_X.
$$
It points up the <nutrient> <concentration> <gradient>, and its negative spatial <derivative> enters the balance equation. The positive $N_0$ prevents singular sensitivity at zero <nutrient> <concentration>; increasing the background <nutrient> <concentration> reduces the drift response to a fixed <gradient>. The <logistic growth equation> term $rB(1-B/K)$ represents net reproduction with low-population growth rate $r$ and <carrying capacity> $K$.
The <nutrient> undergoes <diffusion> with coefficient $D_N$. Its production term $aB^2/(B_0^2+B^2)$ is an increasing, bounded response to the <bacteria>, reaching half its maximum $a$ at $B=B_0$. It has the algebraic form of a <Hill equation> with exponent two; the equation alone does not specify a molecular mechanism. The term $dBN$ removes <nutrient> at a rate proportional to both <bacteria> abundance and <nutrient> <concentration>.
The full <fluxes> are
$$
\boxed{J_B=-D_B B_X+\frac{\chi_BB}{(N_0+N)^2}N_X,\qquad J_N=-D_NN_X.}
$$
The homogeneous <Neumann boundary conditions> make both <fluxes> vanish at the endpoints, so neither <bacteria> nor <nutrient> crosses the boundaries. The <initial conditions> specify their starting spatial distributions. This model includes <nutrient> production by <bacteria>, rather than a separate external supply term.
Back to article page
