Model equations that govern the evolution of internal gravity waves at the interface of two immiscible inviscid fluids are derived. These models follow from the original Euler equations under the sole assumption that the waves are long compared to the undisturbed thickness of one of the fluid layers. No smallness assumption on the wave amplitude is made. Both shallow and deep water configurations are considered, depending on whether the waves are assumed to be long with respect to the total undisturbed thickness of the fluids or long with respect to just one of the two layers, respectively. The removal of the traditional weak nonlinearity assumption is aimed at improving the agreement with the dynamics of Euler equations for large-amplitude waves. This is obtained without compromising much of the simplicity of the previously known weakly nonlinear models. Compared to these, the fully nonlinear models' most prominent feature is the presence of additional nonlinear dispersive terms, which coexist with the typical linear dispersive terms of the weakly nonlinear models. The fully nonlinear models contain the Korteweg–de Vries (KdV) equation and the Intermediate Long Wave (ILW) equation, for shallow and deep water configurations respectively, as special cases in the limit of weak nonlinearity and unidirectional wave propagation. In particular, for a solitary wave of given amplitude, the new models show that the characteristic wavelength is larger and the wave speed is smaller than their counterparts for solitary wave solutions of the weakly nonlinear equations. These features are compared and found in overall good agreement with available experimental data for solitary waves of large amplitude in two-fluid systems.

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.

The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.

The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.

A direct numerical simulation is presented of an elliptical instability observed in the laboratory within an elliptically distorted, rapidly rotating, fluid-filled cylinder (Malkus 1989). Generically, the instability manifests itself as the pairwise resonance of two different inertial modes with the underlying elliptical flow. We study in detail the simplest ‘subharmonic’ form of the instability where the waves are a complex conjugate pair and which at weakly supercritical elliptical distortion should ultimately saturate at some finite amplitude (Waleffe 1989; Kerswell 1992). Such states have yet to be experimentally identified since the flow invariably breaks down to small-scale disorder. Evidence is presented here to support the argument that such weakly nonlinear states are never seen because they are either unstable to secondary instabilities at observable amplitudes or neighbouring competitor elliptical instabilities grow to ultimately disrupt them. The former scenario confirms earlier work (Kerswell 1999) which highlights the generic instability of inertial waves even at very small amplitudes. The latter represents a first numerical demonstration of two competing elliptical instabilities co-existing in a bounded system.

Numerical studies of two-dimensional, transonic flows of dense gases of retrograde type, known as BZT gases, around thin airfoils are presented. The computations are guided by a recent asymptotic theory of Rusak & Wang (1997). It provides a uniformly valid solution of the flow around the entire airfoil surface which is composed of outer and inner solutions. A new transonic small-disturbance (TSD) equation solver is developed to compute the nonlinear BZT gas flow in the outer region around most of the airfoil. The flow in the inner region near the nose of the airfoil is computed by solving the problem of a sonic flow around a parabola. Numerical results of the composite solutions calculated from the asymptotic formula are compared with the solutions of the Euler equations. The comparison demonstrates that, in the leading order, the TSD solutions of BZT gas flows represent the essence of the flow character around the airfoil as computed from the Euler equations. Furthermore, guided by the asymptotic formula, the computational results demonstrate the similarity rules for transonic flows of BZT gases. There are differences between the self-similar cases that may be related to the error associated with the accuracy of the asymptotic solution. A discussion on the flow patterns around an airfoil at transonic speeds and at various upstream thermodynamic conditions is also presented. The paper provides important guidelines for future studies on this subject.

Surface waves superimposed upon a larger-scale flow are blocked and reflected at the points where the group velocities balance the convection by the larger-scale flow. In this study, we first extended the theory of Shyu & Phillips (1990) to the situation when short deep-water gravity waves propagate obliquely upon a steady unidirectional irrotational current and are reflected by it. In this case, the uniformly valid solution and the WKBJ solution of the short waves were derived from the Laplace equation and the kinematical and dynamical boundary conditions. These solutions in terms of some parameters (the expressions for which have also been deduced in this case) take the same forms as those derived by Shyu & Phillips, which by referring to Smith's (1975) theory can even be proved to be valid for gravity waves in an intermediate-depth region and near a curved moving caustic induced by an unsteady multidirectional irrotational current. In this general case, the expressions for certain parameters in these solutions cannot be obtained so that their values must be estimated in a numerical calculation. The algorithm for estimates of some of these parameters that are responsible for the amplitude of the reflected wave not being equal to that of the incident wave in the vicinity of the caustic and therefore are crucial for the computer calculation of the ray solution to be continued after reflection, was illustrated through numerical tests. This algorithm can avoid the error magnification phenomenon that occurred in the previous estimates of the reflected wave in the vicinity of the caustic using the action conservation principle directly. The forms of the solutions have also been utilized to clarify the wave profiles near caustics in a general situation, which indicate that in storm conditions freak waves characterized by a steeper forward face preceded by a deep trough will probably occur in the caustic regions.

The breaking of internal waves propagating in a stratified fluid of constant buoyancy frequency on a sloping boundary was investigated numerically. It was found that at the boundary, nonlinear non-resonant interactions between the incident and reflected waves produced higher-mode waves. These modes had frequencies greater than the local buoyancy frequency and so could not radiate from the interaction region. The energy level of trapped waves increased with time and subsequently led to overturning of the density field. At the critical frequency, when the reflected wave propagated in a direction parallel to the slope, wave overturning occurred near the wall, but the point of overturning moved off the bottom as the propagation angle changed away from that of the bottom slope as the waves became increasingly supercritical. The internal wave reflection coefficient generally increased as the effects of nonlinearity and viscosity decreased, but depended strongly on the forcing frequency and the angle of the sloping boundary.

This paper concerns a relaxation of the assumption of uniform mixture composition in the interior of sonoluminescence bubbles. Intense temperature and pressure gradients within the bubble drive relative mass diffusion which overwhelms diffusion driven by concentration gradients. This thermal and pressure diffusion results in a robust compositional inhomogeneity in the bubble which lasts several orders of magnitude longer than the temperature peak or light pulse at the main collapse of the bubble. This effect has important consequences for control of sonoluminescence, gas dynamics, sonochemistry, and the physics of light production.

Unsteady flow in collapsible tubes has been widely studied for a number of different physiological applications; the principal motivation for the work of this paper is the study of blood flow in the jugular vein of an upright, long-necked subject (a giraffe). The one-dimensional equations governing gravity- or pressure-driven flow in collapsible tubes have been solved in the past using finite-difference (MacCormack) methods. Such schemes, however, produce numerical artifacts near discontinuities such as elastic jumps. This paper describes a numerical scheme developed to solve the one-dimensional equations using a more accurate upwind finite volume (Godunov) scheme that has been used successfully in gas dynamics and shallow water wave problems. The adapatation of the Godunov method to the present application is non-trivial due to the highly nonlinear nature of the pressure–area relation for collapsible tubes.

The code is tested by comparing both unsteady and converged solutions with analytical solutions where available. Further tests include comparison with solutions obtained from MacCormack methods which illustrate the accuracy of the present method.

Finally the possibility of roll waves occurring in collapsible tubes is also considered, both as a test case for the scheme and as an interesting phenomenon in its own right, arising out of the similarity of the collapsible tube equations to those governing shallow water flow.

The intermittency of pressure and pressure gradient in stationary isotropic turbulence at low to moderate Reynolds numbers is studied by direct numerical simulation (DNS) and theoretically. The energy spectra scale in Kolmogorov units as required by the universal-equilibrium hypothesis, but the pressure spectra do not. It is found that the variances of the pressure and pressure gradient are larger than those computed using the Gaussian approximation for the fourth-order moments of velocity, and that the variance of the pressure gradient, normalized by Kolmogorov units, increases roughly as [Rscr ]1/2λ, where [Rscr ]λ is the Taylor microscale Reynolds number. A theoretical explanation of the Reynolds number dependence is presented which assumes that the small-scale pressure field is driven by coherent small-scale vorticity–strain domains. The variance of the pressure gradient given by the model is the product of the variance of ui,juj,i, the source term of the Poisson equation for pressure, and the square of an effective length of the small-scale coherent vorticity–strain structures. This length can be expressed in terms of the Taylor and Kolmogorov microscales, and the ratio between them gives the observed Reynolds number dependence. Formal asymptotic matching of the spectral scaling observed at small scales in the DNS with the classical scaling at large scales suggests that at high Reynolds numbers the pressure spectrum in these forced flows consists of three scaling ranges which are joined by two inertial ranges, the classical k−7/3 range and a k−5/3 range at smaller scale. It is not possible, within the classical Kolmogorov theory, to determine the length scale at which the inertial range transition occurs because information beyond the energy dissipation rate is required.

The effect of solute concentration on hindered diffusion of sphere-like colloidal solutes in stiff polymer hydrogels is examined theoretically and experimentally. In the theoretical development, it is shown that the presence of the gel fibres enhances the effect of concentration on the thermodynamic driving force for gradient diffusion, while simultaneously reducing the effect of concentration on the hydrodynamic drag. The result is that gradient diffusion depends more strongly on solute concentration in gels than it does in pure solution, by an amount that depends on the partition coefficient and hydraulic permeability of the gel–solute system. Quantitative calculations are made to determine the concentration-dependent diffusivity correct to first order in solute concentration. In order to compare the theoretical predictions with experimental data, rates of diffusion have been measured for nonionic micelles and globular proteins in solution and agarose hydrogels at two gel concentrations. The measurements were performed by using holographic interferometry, through which one monitors changes in refractive index as gradient diffusion takes place within a transparent gel. If the solutes are modelled as spheres with short-range repulsive interactions, then the experimentally measured concentration dependence of the diffusivities of both the protein and micelles is in good agreement with the theoretical predictions.

In this paper we consider the causal response of the inviscid shear-layer flow over an elastic surface to excitation by a time-harmonic line force. In the case of uniform flow, Brazier-Smith & Scott (1984) and Crighton & Oswell (1991) have analysed the long-time limit of the response. They find that the system is absolutely unstable for sufficiently high flow speeds, and that at lower speeds there exist certain anomalous neutral modes with group velocity directed towards the driver (in contradiction of the usual radiation condition of out-going disturbances). Our aim in this paper is to repeat their analysis for more realistic shear profiles, and in particular to determine whether or not the uniform-flow results can be regained in the limit in which the shear-layer thickness on a length scale based on the fluid loading, denoted ε, becomes small. For a simple broken-line linear shear profile we find that the results are qualitatively similar to those for uniform flow. However, for the more realistic Blasius profile very significant differences arise, essentially due to the presence of the critical layer. In particular, we find that as ε → 0 the minimum flow speed required for absolute instability is pushed to considerably higher values than was found for uniform flow, leading us to conclude that the uniform-flow problem is an unattainable singular limit of our more general problem. In contrast, we find that the uniform-flow anomalous modes (written as exp (ikx − iωt), say) do persist for non-zero shear over a wide range of ε, although now becoming non-neutral. Unlike the case of uniform flow, however, the k-loci of these modes can now change direction more than once as the imaginary part of ω is increased, and we describe the connection between this behaviour and local properties of the dispersion function. Finally, in order to investigate whether or not these anomalous modes might be realizable at a finite time after the driver is switched on, we evaluate the double Fourier inversion integrals for the unsteady flow numerically. We find that the anomalous mode is indeed present at finite time, once initial transients have propagated away, not only for impulsive start-up but also when the forcing amplitude is allowed to grow slowly from a small value at some initial instant. This behaviour has significant implications for the application of standard radiation conditions in wave problems with mean flow.