We study the dynamics of simple reactions where the chemical species are confined on a general, time-modulated surface, and subjected to externally-imposed stirring. The study of these inhomogeneous effects requires a model based on a reaction-advect
ion-diffusion equation, which we derive. We use homogenization methods to show that up to second order in a small scaling parameter, the modulation effects on the concentration field are asymptotically equivalent for systems with or without stirring. This justifies our consideration of the simpler reaction-diffusion model, where we find that by modulating the substrate, we can modify the reaction rate, the total yield from the reaction, and the speed of front propagation. These observations are confirmed in three numerical case studies involving the autocatalytic and bistable reactions on the torus and a sinusoidally-modulated substrate
The Camassa-Holm equation (CH) is a well known integrable equation describing the velocity dynamics of shallow water waves. This equation exhibits spontaneous emergence of singular solutions (peakons) from smooth initial conditions. The CH equation h
as been recently extended to a two-component integrable system (CH2), which includes both velocity and density variables in the dynamics. Although possessing peakon solutions in the velocity, the CH2 equation does not admit singular solutions in the density profile. We modify the CH2 system to allow dependence on average density as well as pointwise density. The modified CH2 system (MCH2) does admit peakon solutions in velocity and average density. We analytically identify the steepening mechanism that allows the singular solutions to emerge from smooth spatially-confined initial data. Numerical results for MCH2 are given and compared with the pure CH2 case. These numerics show that the modification in MCH2 to introduce average density has little short-time effect on the emergent dynamical properties. However, an analytical and numerical study of pairwise peakon interactions for MCH2 shows a new asymptotic feature. Namely, besides the expected soliton scattering behavior seen in overtaking and head-on peakon collisions, MCH2 also allows the phase shift of the peakon collision to diverge in certain parameter regimes.
We investigate the emergence of singular solutions in a non-local model for a magnetic system. We study a modified Gilbert-type equation for the magnetization vector and find that the evolution depends strongly on the length scales of the non-local e
ffects. We pass to a coupled density-magnetization model and perform a linear stability analysis, noting the effect of the length scales of non-locality on the systems stability properties. We carry out numerical simulations of the coupled system and find that singular solutions emerge from smooth initial data. The singular solutions represent a collection of interacting particles (clumpons). By restricting ourselves to the two-clumpon case, we are reduced to a two-dimensional dynamical system that is readily analyzed, and thus we classify the different clumpon interactions possible.
The Cahn-Hilliard equation describes phase separation in binary liquids. Here we study this equation with spatially-varying sources and stirring, or advection. We specialize to symmetric mixtures and time-independent sources and discuss stirring stra
tegies that homogenize the binary fluid. By measuring fluctuations of the composition away from its mean value, we quantify the amount of homogenization achievable. We find upper and lower bounds on our measure of homogenization using only the Cahn-Hilliard equation and the incompressibility of the advecting flow. We compare these theoretical bounds with numerical simulations for two model flows: the constant flow, and the random-phase sine flow. Using the sine flow as an example, we show how our bounds on composition fluctuations provide a measure of the effectiveness of a given stirring protocol in homogenizing a phase-separating binary fluid.
