Paper 3, Section II, C

Mathematical Biology
Part II, 2012

Consider the two-variable reaction-diffusion system

ut=au+u2v+2uvt=bu2v+d2v\begin{gathered} \frac{\partial u}{\partial t}=a-u+u^{2} v+\nabla^{2} u \\ \frac{\partial v}{\partial t}=b-u^{2} v+d \nabla^{2} v \end{gathered}

where a,ba, b and dd are positive constants.

Show that there is one possible spatially homogeneous steady state with u>0u>0 and v>0v>0 and show that it is stable to small-amplitude spatially homogeneous disturbances provided that γ<β\gamma<\beta, where

γ=bab+a and β=(a+b)2\gamma=\frac{b-a}{b+a} \quad \text { and } \quad \beta=(a+b)^{2}

Now assuming that the condition γ<β\gamma<\beta is satisfied, investigate the stability of the homogeneous steady state to spatially varying perturbations by considering the timedependence of disturbances whose spatial form is such that 2u=k2u\nabla^{2} u=-k^{2} u and 2v=k2v\nabla^{2} v=-k^{2} v, with kk constant. Show that such disturbances vary as epte^{p t}, where pp is one of the roots of

p2+(βγ+dk2+k2)p+dk4+(βdγ)k2+βp^{2}+\left(\beta-\gamma+d k^{2}+k^{2}\right) p+d k^{4}+(\beta-d \gamma) k^{2}+\beta

By comparison with the stability condition for the homogeneous case above, give a simple argument as to why the system must be stable if d=1d=1.

Show that the boundary between stability and instability (as some combination of β,γ\beta, \gamma and dd is varied) must correspond to p=0p=0.

Deduce that dγ>βd \gamma>\beta is a necessary condition for instability and, furthermore, that instability will occur for some kk if

d>βγ{1+2γ+21γ+1γ2}.d>\frac{\beta}{\gamma}\left\{1+\frac{2}{\gamma}+2 \sqrt{\frac{1}{\gamma}+\frac{1}{\gamma^{2}}}\right\} .

Deduce that the value of k2k^{2} at which instability occurs as the stability boundary is crossed is given by

k2=βdk^{2}=\sqrt{\frac{\beta}{d}}