--- title: Doubly Degenerate Nutrient Taxis Model url: https://www.emergentmind.com/topics/doubly-degenerate-nutrient-taxis-model type: topic --- # Doubly Degenerate Nutrient Taxis Model The doubly degenerate nutrient taxis model denotes a class of chemotaxis-consumption systems in which the transport of a motile population \(u\) depends simultaneously on population density and nutrient or signal concentration \(v\), so that diffusion and, in many formulations, the taxis flux itself vanish when either \(u\) or \(v\) becomes small. Representative equations include the one-dimensional system \[ u_t=(uvu_x)_x-(u^2vv_x)_x+uv,\qquad v_t=v_{xx}-uv, \] its multidimensional variants \[ u_t=\nabla\cdot(uv\nabla u)-\chi\nabla\cdot(u^\alpha v\nabla v)+\ell uv,\qquad v_t=\Delta v-uv, \] and parabolic-elliptic consumption systems of the form \[ u_t=\nabla\cdot(D(u)\nabla u)-\nabla\cdot(u\nabla v),\qquad 0=\Delta v-uv. \] The class arose in the modeling of bacterial aggregation, especially for \emph{Bacillus subtilis} in nutrient-poor environments, and has developed into a distinct branch of degenerate cross-diffusion theory with results on derivation, weak and classical solvability, uniqueness, boundedness, sharp thresholds, and long-time stabilization [1301.5058], [1711.03015], [2508.03268]. ## 1. Defining structure and principal model families The defining feature is the simultaneous degeneration of bacterial transport with respect to more than one state variable. In the prototypical cross-diffusion form, \[ u_t=\nabla\cdot(uv\nabla u)-\nabla\cdot(u^2v\nabla v)+\text{reaction terms}, \] the diffusion coefficient \(uv\) vanishes if either \(u=0\) or \(v=0\), while the taxis coefficient \(u^2v\) vanishes even more strongly in low-density or low-nutrient regions. This is the sense in which the model is “doubly degenerate” in the literature on bacterial nutrient taxis [2501.06821], [2602.21164]. Several closely related formulations appear in current research. One family replaces \(uv\nabla u\) by \(u^{m-1}v\nabla u\) and \(u^2v\nabla v\) by \(f(u)v\nabla v\), leading to \[ u_t=\nabla\cdot\!\left(u^{m-1}v\nabla u\right)-\nabla\cdot\!\left(f(u)v\nabla v\right)+\ell uv,\qquad v_t=\Delta v-uv, \] with growth conditions such as \(f(s)\le C_fs^\alpha\) [2409.02741]. Another family uses the exponent \(l\) in \[ u_t=\nabla\cdot(u^{l-1}v\nabla u)-\nabla\cdot(u^lv\nabla v)+uv,\qquad v_t=\Delta v-uv, \] or its logistic counterpart with \(u-u^2\) במקום \(uv\) [2411.01133], [2410.19833], [2601.04845]. A further variant places the degeneracy in the signal dependence of the diffusion, \[ u_t=\Delta(u\phi(v)),\qquad v_t=\Delta v-uv, \] with \(\phi(0)=0\) and \(\phi'(0)>0\), so that the \(u\)-equation loses ellipticity as nutrient approaches vacuum [2209.12724]. The parabolic-elliptic consumption system \[ u_t=\nabla\cdot(D(u)\nabla u)-\nabla\cdot(u\nabla v),\qquad 0=\Delta v-uv \] is not doubly degenerate in the same explicit \(uv\)-sense, but it belongs to the same nutrient-taxis/consumption line and exhibits a sharp diffusion threshold under radial symmetry [2601.05338]. Across these formulations, \(u\) is interpreted as bacterial or organism density and \(v\) as nutrient or signal concentration. The term \(-uv\) in the \(v\)-equation models consumption, while \(uv\) or \(\ell uv\) in the \(u\)-equation models nutrient-supported proliferation. No-flux boundary conditions are standard on bounded domains; the whole-line Cauchy problem and parabolic-elliptic problems with prescribed nutrient boundary data also occur [2508.07503], [2601.05338]. ## 2. Biological motivation and kinetic derivation The biological origin of the class lies in models for bacterial colony formation on thin agar plates. A foundational formulation extends the Kawasaki et al. aggregation model by taking the bacterial diffusivity to be nonlinear and degenerate, \[ D_b=\sigma nb, \] and by introducing a nutrient-chemotactic flux \[ J_c=\sigma n b^2 \chi(n)\nabla n,\qquad \chi(n)=\frac{\chi_0K_d}{(K_d+n)^2}, \] which is explicitly constructed to be compatible with the same nonlinear mobility law [1301.5058]. After nondimensionalization, this yields \[ n_t=\Delta n-nb,\qquad b_t=\nabla\cdot(\sigma nb\nabla b)+nb-\chi_0\nabla\cdot\!\left(\frac{\sigma nb^2}{(1+n)^2}\nabla n\right). \] Here the transport vanishes when either bacteria or nutrient are scarce, reflecting reduced motility in nutrient-poor or low-density regions. A formal kinetic justification was later given through the parabolic limit of a velocity-jump process of Othmer-Dunbar-Alt type [1711.03015]. In that derivation, individual-cell motion is described by a Boltzmann-type transport equation with an isotropic leading turning rate \[ \lambda_0(\xi,\tau)=\frac{\mu_0}{\rho(\xi,\tau)S(\xi,\tau)} \] and a Schnitzer-type first-order perturbation \[ \lambda_1(v,\hat S)=\kappa(S)\,(v\cdot\nabla_\xi S). \] Under parabolic scaling and a Hilbert expansion, the effective diffusion tensor becomes proportional to \(\rho S\), and the chemotactic drift becomes proportional to \(\rho^2S\kappa(S)\nabla S\). The resulting macroscopic bacterial equation has the form \[ \partial_\tau \rho =\nabla_\xi\cdot(\sigma \rho S\nabla_\xi\rho) -\nabla_\xi\cdot\!\left(\sigma \rho^2S\frac{\chi_0K_d}{(K_d+S)^2}\nabla_\xi S\right), \] which recovers the doubly degenerate diffusion and nutrient-taxis flux of the phenomenological model [1711.03015]. This derivation is significant because it identifies the degeneracy as a consequence of microscopic movement rules rather than a purely ad hoc constitutive choice. It also explains why diffusion and taxis carry the same nutrient-density prefactors: both emerge from the same turning-rate mechanism. ## 3. Solvability theory across dimensions and geometries Existence theory depends strongly on dimension, whether the nutrient equation is parabolic or elliptic, and whether damping terms are present. The current state of the theory covers whole-space, bounded-domain, parabolic-parabolic, and parabolic-elliptic settings. | Setting | Representative feature | Main conclusion | |---|---|---| | \(\mathbb R\), Cauchy problem [2508.07503] | doubly degenerate diffusion and taxis | global weak solution under \(u_0,v_0\in W^{1,\infty}\cap L^1\) and extra entropy/gradient conditions | | 1D bounded interval [2602.21164] | parabolic-elliptic model with external supply \(f\not\equiv 0\) | global weak solution for any nonnegative \(u_0\in W^{1,\infty}\) | | 2D bounded convex domain [2405.20637] | \(u_t=\nabla\cdot(uv\nabla u)-\chi\nabla\cdot(u^2v\nabla v)+\ell uv\) | global bounded weak solution, continuous in \(u\), essentially smooth in \(v\) | | 2D smooth domain [2601.04845] | quadratic logistic source \(u-u^2\) | global weak solution uniformly bounded in time | | 3D bounded convex domain [2508.03268] | bacterial response parameter \(0\le \alpha<2\) | global weak solvability for all \(0\le \alpha<2\) | | Ball in arbitrary dimension [2601.05338] | parabolic-elliptic consumption system with radial data | global bounded classical solution for \(\alpha<1\) | In one space dimension, the bounded-interval system \[ u_t=(uvu_x)_x-(u^2vv_x)_x+uv,\qquad v_t=v_{xx}-uv \] already has a global weak-solution theory, and uniqueness is available under the additional condition \(1/u^0\in L^1(\Omega)\) [2501.06821]. On the whole line, global weak solutions were obtained by regularizing on expanding intervals \(B_{1/\varepsilon}=(-1/\varepsilon,1/\varepsilon)\), deriving \(\varepsilon\)-independent entropy bounds, and passing to the limit via Aubin-Lions [2508.07503]. For the parabolic-elliptic problem \[ u_t=(uvu_x)_x-(u^2vv_x)_x+uv,\qquad 0=v_{xx}-uv+f(x,t), \] global weak solvability holds for nonnegative \(u_0\in W^{1,\infty}(\Omega)\) and \(f\in C^1(\bar\Omega\times[0,\infty))\), \(f\not\equiv0\) [2602.21164]. In two dimensions, the model \[ u_t=\nabla\cdot(uv\nabla u)-\chi\nabla\cdot(u^2v\nabla v)+\ell uv,\qquad v_t=\Delta v-uv \] admits global bounded weak solutions in smooth bounded convex domains for reasonably regular initial data [2405.20637]. The logistic variant \[ u_t=\nabla\cdot(uv\nabla u)-\nabla\cdot(u^2v\nabla v)+u-u^2,\qquad v_t=\Delta v-uv \] admits global weak solutions uniformly bounded in time on general smooth planar domains; the paper attributes the strengthening primarily to the quadratic degradation term \(-u^2\) [2601.04845]. More general two-dimensional large-data results are available for \[ u_t=\nabla\cdot(u^{m-1}v\nabla u)-\nabla\cdot(f(u)v\nabla v)+\ell uv,\qquad v_t=\Delta v-uv, \] with weak solvability for \(1\le m<3\) and classical solvability for \(3\le m<4\) under explicit parameter ranges [2409.02741]. For the \(u^{l-1}v\nabla u-u^lv\nabla v\) family, global weak solutions exist in 2D for \(l\in[1,3]\), and uniform boundedness is established for \(1\max\{2,n\}\), generates a unique global bounded classical solution in arbitrary dimension \(n\ge1\). Earlier work gives nonexistence of global solutions for some radial data when \(D(\xi)\le K_D(\xi+1)^{-\alpha}\) with \(\alpha>1\). The critical exponent \[ \alpha_c^{(\mathrm{cons})}=1 \] is therefore sharp and dimension-independent for this nutrient-taxis/consumption system [2601.05338]. For explicitly doubly degenerate parabolic-parabolic systems, the thresholds are more geometry- and structure-sensitive. In the family \[ u_t=\nabla\cdot(u^{l-1}v\nabla u)-\nabla\cdot(u^lv\nabla v)+uv, \] the one-dimensional theory yields global boundedness for all \(l\ge1\), whereas the two-dimensional theory yields global weak solvability for \(l\in[1,3]\) and uniform boundedness for \(1