For the Brusselator model with parameters

PHYSICS · Multi-part · Fermi University Physics · Thermal, Statistical & Complex Systems · Complex and collective systems · Complex and collective systems

Many Parts, One Behaviour, Ch Pattern formation: instabilities that choose a wavelength

Problem

For the Brusselator model with parameters AA, BB and diffusion coefficients DuD_u, DvD_v (dimensionless units), the Turing threshold is BT=(1+ADu/Dv)2B_{\rm T} = (1+A\sqrt{D_u/D_v})^2, the Hopf (oscillatory) threshold is BH=1+A2B_{\rm H} = 1+A^2, and the critical wavenumber obeys qc2=A/DuDvq_{\rm c}^2 = A/\sqrt{D_uD_v}. Take A=1.5A=1.5, Du=1D_u=1, Dv=10D_v=10. ((a)) Find BTB_{\rm T} and BHB_{\rm H}, and confirm that the Turing instability sets in before the Hopf instability as BB is raised from below. ((b)) Find the critical wavenumber qcq_{\rm c} and the pattern wavelength Λc\Lambda_{\rm c}. ((c)) Assume that, as for the Swift-Hohenberg equation, the Brusselator's stationary bifurcation here is supercritical, so that past threshold the pattern amplitude grows as the square root of the distance above it, A∗(B)∝B−BTA^*(B)\propto\sqrt{B-B_{\rm T}} (confirming this needs a weakly nonlinear calculation not done here). Using this scaling, find the ratio A∗(3.0)/A∗(2.5)A^*(3.0)/A^*(2.5) of the pattern amplitudes at B=3.0B=3.0 and B=2.5B=2.5 (both values lie between BTB_{\rm T} and BHB_{\rm H}, so the system is in the stationary pattern regime at each).

Browse the problem library