The unsteady Hele-Shaw problem is a model nonlinear system that, for a certain parameter ranger, exhibits the phenomenon known as viscous fingering. While not directly applicable to multiphase porous-media flow, it does prove to be an adequate mathematical model for unstable dieplacement in laboratory parallel-plate devices. We seek here to determine, by use of an accurate boundary-integral frount-tracking scheme, the extent to which the simplified system captures the canonical nonlinear behavior of displacement flows and, in particular, to ascertain the role of noise in such systems. We choose to study a particular pattern of injection and production “wells.” The pattern chosen is the isolated “five-spot,” that is a single source surrounded by four symmetrically-placed sinks in an infinite two-dimensional “reservoir.” In cases where the “pusher” fluid has negligible viscosity, sweep efficiency is calculated for a range of values of the single dimensionless parameter τ, an inverse capillary number. As this parameter is reduced, corresponding to increased flow rate or reduced interfacial tension, this efficiency decreases continuously. For small values of τ, these stable displacements change abruptly to a regime characterized by unstable competing fingers and a significant reduction in sweep efficiency. A simple stability argument appears to correctly predict the noise level required to transit from the stable to the competing-finger regimes. Published compilations of experimental results for sweep efficiency as a function of viscosity ratio showed an unexplained divergence when the pusher fluid is less viscous. Our simulations produce a similar divergence when, for a given viscosity ratio, the parameter τ is varied.

The problem of a source or sink submerged beneath a free surface is investigated in the infinite Froude number limit. Solutions are found for all cases in which the source is situated away from the bottom of the channel. Solutions are also found for the case where the source is situated above the asymptotic level of the free surface, giving fountain type free surface shapes.

The modulation of short gravity waves by long waves or currents is described for the situation when the flow is irrotational and when the short waves are described by linearised equations. Two cases are distinguished depending on whether the basic flow can be characterised as a deep-water current, or a shallow-water current. In both cases the basic flow has a current which has finite amplitude, while in the first case the free surface slope of the basic flow can be finite, but in the second case is small. The modulation equations are the local dispersion relation of the short waves, the kinematic equation for conservation of wave crests and the wave action equation. The results incorporate and extend the earlier work of Longuet-Higgins and Stewart [10, 11].

The leading-order interaction of short gravity waves with a dominant long-wave swell is calculated by means of Zakharov's [7] spectral formulation. Results are obtained both for a discrete train of short waves and for a localised wave packet comprising a spectrum of short waves.

The results for a discrete wavetrain agree with previous work of Longuet-Higgins & Stewart [5], and general agreement is found with parallel work of Grimshaw [4] which employed a very different wave-action approach.

A quasi-steady-state apprcncimation to the Navier-Stokes equation is the corresponding equation with nonhomogeneous forcing term f(x, t), but with the term Vt deleted. For solutions that are zero on the boundary, the difference z between the solution of the Navier-Stokes equation and the solution of this quasi-steady-state approximation is estimated in the L2 norm ║z║ with respect to the spatial variables. For sufficiently large viscosity or sufficiently small body force f, the inequality

holds for 0 < t ≤ T and certain real numbres C, β > 0.

The two-dimensional Helmholtz equation is studied for an infinite region with two semi-infinite plates extending to infinity in opposite directions and a finite duct in the overlapping region. The solution technique leads to coupled Wiener-Hopf equations, and subsequently to an infinite set of simultaneous linear equations. As an example, an asymptotic expansion is calculated and graphed for the case when the duct length divided by duct width is large enough to ensure damping of all but the zero mode wave in the duct.

In this paper, we study controllability and observability problems for the wave and heat equation in a spherical region in Rn, where the control enters in the mixed boundary condition. In the main result, we show that all “finite energy” initial states (i.e. (ω0, ν0) ∈ H1(Ω) × L2 (Ω)) can be steered to zero at time T, using a control f ∈ L2 (∂Ω × [0, T]), provided T > 2. On this basis, we use the duality principle to investigate initial observability for the wave equation. Applying the Fourier transform technique, we obtain controllability and observability results for the heat equation.

In many physical problems, the system tends quickly to a particular structure, which then evolves relatively slowly in space and time. Various methods exist to derive equations describing the slow evolution of the particular structure; for example, the method of multiple scales. However, the resulting equations are typically valid only for a limited range of the parameters. In order to extend the range of validity and to improve the accuracy, correction terms must be found for the equations. Here we describe a procedure, inspired by centre-manifold theory, which provides a systematic approach to calculating a sequence of successively more accurate approximations to the evolution of the principal structure in space and time.

The formal procedure described here raises a number of questions for future research. For example: what sort of error bounds can be obtained, do the approximations converge or are they strictly asymptotic, and what sort of boundary conditions are appropriate in a given problem?