No abstract

Abstract. Pattern formation in systems with a conserved quantity is considered by studying the appropriate amplitude equations. The conservation law leads to a large-scale neutral mode that must be included in the asymptotic analysis for pattern formation near onset. Near a stationary bifurcation, the usual Ginzburg-Landau equation for the amplitude of the pattern is then coupled to an equation for the large-scale mode. These amplitude equations show that for certain parameters all roll-type solutions are unstable. This new instability differs from the Eckhaus instability in that it is amplitude-driven and is supercritical. Beyond the stability boundary, there exist stable stationary solutions in the form of strongly modulated patterns. The envelope of these modulations is calculated in terms of Jacobi elliptic functions and, away from the onset of modulation, is closely approximated by a sech profile. Numerical simulations indicate that as the modulation becomes more pronounced, the envelope broadens. A number of applications are considered, including convection with fixed-flux boundaries and convection in a magnetic field, resulting in new instabilities for these systems.

We consider the flow of a viscous incompressible fluid in a parallel-walled channel, driven by steady uniform suction through the porous channel walls. A similarity transformation reduces the Navier-Stokes equations to a single partial differential equation (PDE) for the stream function, with two-point boundary conditions. We discuss the bifurcations of the steady solutions first, and show how a pitchfork bifurcation is unfolded when a symmetry of the problem is broken.Then we describe time-dependent solutions of the governing PDE, which we calculate numerically. We analyse these unsteady solutions when there is a high rate of suction through one wall, and the other wall is impermeable: there is a limit cycle composed of an explosive phase of inviscid growth, and a slow viscous decay. The inviscid phase ‘almost’ has a finite-time singularity. We discuss whether solutions of the governing PDE, which are exact solutions of the Navier-Stokes equations, may develop mathematical singularities in a finite time.When the rates of suction at the two walls are equal so that the problem is symmetrical, there is an abrupt transition to chaos, a ‘homoclinic explosion’, in the time-dependent solutions as the Reynolds number is increased. We unfold this transition by perturbing the symmetry, and compare direct numerical integrations of the governing PDE with a recent theory for ‘Lorenz-like’ dynamical systems. The chaos is found to be very sensitive to symmetry breaking.

Topological chaos may be used to generate highly effective laminar mixing in a simple batch stirring device. Boyland, Aref & Stremler (2000) have computed a material stretch rate that holds in a chaotic flow, provided it has appropriate topological properties, irrespective of the details of the flow. Their theoretical approach, while widely applicable, cannot predict the size of the region in which this stretch rate is achieved. Here, we present numerical simulations to support the observation of Boyland et al. that the region of high stretch is comparable with that through which the stirring elements move during operation of the device. We describe a fast technique for computing the velocity field for either inviscid, irrotational or highly viscous flow, which enables accurate numerical simulation of dye advection. We calculate material stretch rates, and find close agreement with those of Boyland et al., irrespective of whether the fluid is modelled as inviscid or viscous, even though there are significant differences between the flow fields generated in the two cases. IntroductionStatic and dynamic mixing devices are important in many industries, e. ) have demonstrated, in an unusual blend of ad hoc experimentation and abstract mathematics, that flows with the topology of certain braids achieve a material stretch rate which can be determined quantitatively, given only the topology of the flow. However, a key feature not predicted by their theoretical considerations is the size of the domain in which this stretch rate is attained. Indeed, according to the theory, this domain may have measure zero, and if this were the case then the theory would have little practical impact. Here we provide numerical results that support the observations of Boyland et al., that the chaotic region is in fact commensurate with the region of fluid through which the stirring elements move during operation of the device. We should make clear at the outset that we use the terminology 'topological chaos' in the same sense as Boyland et al. (2000), to

Motivated by devices such as the atomic force microscope, we compute the drag experienced by a cylindrical body of circular or rectangular cross-section oscillating at small amplitude near a plane wall. The body lies parallel to the wall and oscillates normally to it; the body is assumed to be long enough for the dominant flow to be two-dimensional. The flow is parameterized by a frequency parameter γ 2 (a Strouhal number) and the wall-body separation ∆ (scaled on body radius). Numerical solutions of the unsteady Stokes equations obtained using finite-difference computations in bipolar coordinates (for circular cross-sections) and boundary-element computations (for rectangular cross-sections) are used to determine the drag on the body. Numerical results are validated and extended using asymptotic predictions (for circular cylinders) obtained at all extremes of (γ, ∆)-parameter space. Regions in parameter space for which the wall has a significant effect on drag are identified.

The long-time behaviour of many dynamical systems may be effectively predicted by a low-dimensional model that describes the evolution of a reduced set of variables. We consider the question of how to equip such a low-dimensional model with appropriate initial conditions, so that it faithfully reproduces the long-term behaviour of the original high-dimensional dynamical system. Our method involves putting the dynamical system into normal form, which not only generates the low-dimensional model, but also provides the correct initial conditions for the model. We illustrate the method with several examples.

No abstract

We study homoclinic snaking in the cubic-quintic Swift-Hohenberg equation (SHE) close to the onset of a subcritical pattern-forming instability. Application of the usual multiple-scales method produces a leading-order stationary front solution, connecting the trivial solution to the patterned state. A localized pattern may therefore be constructed by matching between two distant fronts placed back-to-back. However, the asymptotic expansion of the front is divergent, and hence should be truncated. By truncating optimally, such that the resultant remainder is exponentially small, an exponentially small parameter range is derived within which stationary fronts exist. This is shown to be a direct result of the 'locking' between the phase of the underlying pattern and its slowly varying envelope. The locking mechanism remains unobservable at any algebraic order, and can only be derived by explicitly considering beyond-allorders effects in the tail of the asymptotic expansion, following the method of Kozyreff and Chapman as applied to the quadratic-cubic SHE Kozyreff 2009 Physica D 238 319-54, Kozyreff and Chapman 2006 Phys. Rev. Lett. 97 44502). Exponentially small, but exponentially growing, contributions appear in the tail of the expansion, which must be included when constructing localized patterns in order to reproduce the full snaking diagram. Implicit within the bifurcation equations is an analytical formula for the width of the snaking region. Due to the linear nature of the beyond-all-orders calculation, the bifurcation equations contain an analytically indeterminable constant, estimated in the previous work by Chapman and Kozyreff using a best fit approximation. A more accurate estimate of the equivalent constant in the cubic-quintic case is calculated from the iteration of a recurrence relation, and the subsequent analytical bifurcation diagram compared with numerical simulations, with good agreement.