In a process involving electromagnetic shaping, a high-frequency electromagnetic field is used to deform a liquid conductor into a required shape. This is particularly relevant to processes such as levitation melting. In this paper the stability of such configurations are investigated. The second variation of an appropriate energy functional is derived whose minimum states correspond to stable configurations, thus providing a stability criterion. As an example, this is applied to the shaping of a levitated cylinder of circular cross-section and to an almost spherical axisymmetric shape. In both cases we find that these shapes are unstable. We then consider enclosing the entire shaping device in a metal shield, thus preventing the escape of the magnetic field. It is then shown that in general the shield has a stabilizing effect, whose exact nature depends on the topology of the liquid shape and on the field structure on its surface. This differing behaviour is discussed for two-dimensional spherical and toroidal shapes.

The stages of evolution of a quast-steady breaker from the onest of a capillary pattern to a fully evolved breaking wave are cgaracterized using high-image-density particle image velocimetry, which provides instrantaneous representations of the free surface and the patterns of vorticity beneath it. The initial stage, which sets in at a low value of Froude number, involves a capillary pattern along each trough-crest surface of a quasi-stationary wave. The successive crests of the capillary pattern exhihit increasing scale and culminate in a single largest-scale crest of the free surface. Immediately upstream of the large-scale crest, the capillary pattern shows counterclockwise concentrations of vorticity at its troughs and regions of clockwise vorticity beneath its crests. The onset of the final, largest-scale crest exhihits two forms: one involving no flow sparation; and the other exhibiting a small-scale separaed mixing layer. At an intermediate value of Froude number, a breaker occurs and the acpillary pattern is replaced by large-scale distortions of the free surface. The onset of separation, which involves flow deceleration along a region of the free surface having a large radius of curvature, leads to formation of a long mixing layeer, which has substantial levels of vorticity. Downstream of this breaker, the long-wavelength wave pattern is suppressed. At the largest value of Froude number, the onset of flow sparation rapidly occurs in conjunction with an abrupt change in slope of the surface, giving rise to vorticity concentrationa in the mixing layer.

A fluid droplet subjected to shear flow deforms and rotates in the flow. In the presence of a wall the droplet migrates with respect to a material element in the undisturbed flow field. Neglecting fluid inertia, the Stakes problem for the droplet is solved using a boundary integral technique. It is shown how the time-dependent deformation, orientation, circulation and droplet viscosity. The migration velocities are calculated in the directions parallel and perpendicular to the wall, and compared with theoretical models and expeeriments. The results reveal some of the shortcomings of existiong models although not all diserepancies between our calculations and known experiments could be clarified.

The form and evolution of three-dimensional standing waves in deep water are calculated analytically from Zakharov's equation and computationally from the full nonlinear bounddary value problem. The water is contained in a basin with a square cross-cection, when three-dimensional properties to pairs of sides are the same. It is found that non-periodic standing waves commonly follow forms of cyclic recurrence over times. The two-dimensional Stokes type of periodic standing waves (dominated by the fundamental harmonic) are shown to be unstable to three dimensional disturbances, but over long times the waves return cyclically close to their initial state. In contrast, the three-dimensional Stokes type of periodic standing waves are found to be stabel to small disturbances. New two-dimensional periodic standing waves with amplitude maxima at other than the fundamental harmonic have been investigated recently (Bryant & Stiassnie 1994). The equivalent three-dimensional standing waves are described here. The new two-dimensional periodic standing waves, like the two-dimensional Stokes standing waves, are found to be unstable to three-dimensional disturbances, and to exhibit cyclic recurrence over long times. Only some of the new three-dimensional periodic standing waves are found to be stable to small disturbances.

We have examined the linear stability of the fully developed mixed-convection flow in a differentially heated tall vertical channel under non-Boussinesq conditions. The Three-dimensional analysis of the stability problem was reduced to an equivalent two-dimensional one by the use of Squire's transformation. The resulting eigenvalue problem was solved using an integral Chebyshev pseudo-spectral method. Although Squire's theorem cannot be proved analytically, two-dimensional disturbances are found to be the most unstable in all cases. The influence of the non-Boussinesq effects on the stability was studied. We have investigated the dependence of the critical Grashof and Reynolds numbers on the temperature difference. The results show that four different modes of instability are possible, two of which are new and due entirely to non-Boussinesq effects.

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees.

There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases ∝ (H − 1)−2 for (H − 1) [Lt ] 1 (thickness of wall much less than the tube radius), and decreases ∝ (H−4 for H [Gt ] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

Using the assumptions of an incompressible and viscous flow at large Reynolds number, we derive the evolution equation for the wave action density of an instability wave travelling on top of a laminar free-shear flow. The instability is considered to be viscous; the purpose of the present work is to include the cumulative effect of the (locally) small viscous correction to the wave, over length and time scales on which the underlying base flow appears inhomogeneous owing to its viscous diffusion. As such, we generalize our previous work for inviscid waves. This generalization appears as an additional (but usually non-negligible) term in the equation for the wave action. The basic structure of the equation remains unaltered.

The transport of heavy, polydispersed particles and the inter-phase transfer of kinetic energy due to the viscous drag forces is measured experimentally in a turbulent shear layer. To study the effect of the large-scale vortex pairing event, the shear layer is forced simultaneously with a fundarmental and subharmonic perturbation. It is shown that vortex pairing plays a homogenizing role on the particulate field, but hte amount of homogenization is strongly dependent upon the particle's viscous relaxtion time, the eddy turnover time, as well as the time the particles interact with each scale prior to a pairing event. Thus, even though the smaller size particles become well-mixed across the large eddies, the larger sizes are still dispersed in an inhormogeneous fashion. It is also found that the kinetic energy transfer between the phases occurs inhomogeneously with energy being exchanged predominantly in a sublayer just outside the region of maximum turbulence intensity. The kinetic energy transfer is shown to exhibit notable positive and negative peaks located beneath the cores and stagnation points of the large-scale eddy field, and these peaks are shown to result from the irrotational velocity perturbations created by the vortices. This energy exchange mechanism remains a prominent process as long as the Stokes number of the particles relative to the vortices is of order unity.

It is shown that the expressions for the correlation tensors of homogeneous axisymmetric turbulence can be considerably simplifed compared to previous analyses of Batchelor (1946) and Chandrasekhar (1950). Representations of the axisymmetric two-point correlations tensors are found, such that each measurable correlation corresponds to a single scalar function, and moreover such that the equations of continuity relating different tensor components to each other take the most simple form. Reflectional symmetry in planes normal to but not in planes through the axis of symmetry is demanded, which allows a full description of states with rotation about the axis of symmetry. The second and third-order velocity correlation tensors and the first-order pressure velocity correlation tensor are analysed with the new method. Small separation expansions of the correlation functions yield the quantities which have to be measured to determine various terms in the governing equations for the Reynolds stresses and the dissipation tensor. A scalar Poisson equation for the pressure-strain is derived. and the single-pint solution is written as a sum of integrals over measurable correlation functions. The simplified analysis can be of great experimental importance. It reveals in a simple way how a full experimental picture of homogeneous axisymmetric turbulence can be obtained by measuring components of the velocity at two points at variable distance from each other on a line perpendicular to the mean flow in a wind tunnel. By using the Fouricr Besscl transform it is also shown that the three-demensional energy, transfer, and pressure stain spectra can be extracted from such measurements.

The diffraction of regular waves by a vertical circular cylinder in finite depth water is considered, within the frame of potential theory. The wave slope kA is assumed to be small so that successive boundary value problems at orders kA, k2A2, and k3A3 can be formulated. Here we focus on the third-order (k3A3) problem but restrict ourselves to the triple-frequency component of the diffraction potential. The method of resolution is based on eigenfunction expansions and on the integral equation technique with the classical Green function expressed in cylindrical coordinates. Third-order (triple-frequency) loads are calculated and compared with experimental measurements and approximate methods based on long-wave theories.

A numerical technique, based on the boundary integral method, is developed to allow the modelling of unsteady free-surface flows at large Reynolds numbers in cases where the surface is contaminated by some surface-active compound. This requires the method to take account of the tangential stress condition at the interface and is achieved through a boundary-layer analysis. The constitutive relation that forms the surface stress condition is assumed to be of the Boussinesq type and allows the incorporation of surface shear and dilatational viscous forces as well as Marangoni effects due to gradients in surface tension. Sorption kinetics can be included in the model, allowing calculations for both soluble and insolube surfactants. Application of the numerical model to the problem of bursting gas bubbles at a free surface shows the greatest effect to be due to surface dilatational viscosity which drastically reduces the amount of surface compression and can slow and even prevent the information of a liquid jet. Surface tension gradients give dilatational elasticity to the surface and thus also significantly prevent surface compression. Surface shear viscosity has a smaller effect on the interface motion but results in initially increased surface concentrations due to the sweeping up of surface particles ahead of the inward-moving surface wave.

This is a study of the interactions of solitary waves climbing up a circular island. A series of large-scale laboratory experiments with waves of different incident height-to-depth ratios and different crest lengths is described. Detailed two-dimensional run-up height measurements and time histories of surface elevations around the island are presented. A numerical model based on the two-dimensional shallow-water wave equations including runup calculations was developed. Numerical model predictions agreed very well with the laboratory data and the model was used to study wave trapping and the effect of slope. Under certain conditions, enhanced runup and wave trapping on the lee side of the island were observed, suggesting a possible explanation for the devastation reported by field surveys in Babi Island off Flores, Indonesia, and in Okushiri Island, Japan.

The bifurcation of steady periodic waves from irrotational inviscid streamflows is considered. Normalizing the flux Q to unity leaves two other natural quantities R (pressure head) and S (flowforce) to parameterize the wavetrain. In a well-known paper, Benjamin & Lighthill (1954) presented calculations within a cnoidal-wave theory which suggested that the corresponding values of R and S lie inside the cusped locus traced by the sub- and supercritical streamflows. This rule has been applied since to many other flow scenarios. In this paper, regular expansions for the streamfunction and profile are constructed for a wave forming on a subcritical stream and thence values for R and S are calculated. These describe, locally, how wave brances in (R, S) parameter space point inside the streamflow cusp. Accurate numerics using a boundry-integral solver show how these constant-period branches extend globally and map out parameter space. The main result is to show that the large-amplitude branches for all steady Stokes’ waves lie surprisingly close to the subcritical stream branch, This has important consequences for the feasibility of undular bores (as opposed to hydraulic jumps) in obstructed flow. Moreover, the transition from the ‘long-wave region’ towards the ‘deep-water limit’ is char-acterized by an extreme geometry, bith of the wave branches and how they sit inside each other. It is also shown that a single (Q, R, S) trriple may represent more than one wave since the global branches can overlap in (R, S) parameter space. This non-uniqueness is not that associated with the known premature maxima of wave propertties as functions of wave amplitude near waves of greatest height.

We consider the solidification of a binary alloy in a mushy layer and analyse the system near the onset of buoyancy-driven convection in the layer. We employ a neareutectic approximation and consider the limit of large far-field temperature. These asymptotic limits allow us to examine the rich dynamics of the mushy layer in the form of small deviations from the classical case of convection in a horizontal porous layer of uniform permeability. Of particular interest are the effects of the asymmetries in the basic state and the non-uniform permeability in the mushy layer, which lead to transcritically bifurcating convection with hexagonal planform. We obtain a set of three coupled amplitude equations describing the evolution of small-amplitude convecting states in the mushy layer. These equations are analysed to determine the stability of and competition between two-dimensional roll and hexagonal convection patterns. We find that either rolls or hexagons can be stable. Furthermore, hexagons with either upflow or downflow at the centres can be stable, depending on the relative strengths of different physical mechanisms. We determine how to adjust the control parameters to minimize the degree of subcriticality of the bifurcation and hence render the system globally more stable. Finally, the amplitude equations reveal the presence of a new oscillatory instability.

An immersed boundary technique is used to model a riblet covered surface on one wall of a channel bounding fully developed turbulent flow. The conjecture that the beneficial drag reduction effect of riblets is a result of the damping of cross-flow velocity fluctuations is then examined. This possibility has been discussed by others but is unverified. The damping effect is explicitly modelled by applying a cross-flow damping force field in elongated streamwise zones with a height and spacing corresponding to the riblet crests. The same trends are observed in the turbulence profiles above both riblet and damped surfaces, thus supporting cross-flow damping as a beneficial mechanism. It is found in the examples presented that the effect of the riblets on the mean flow field quantities (mean velocity profile, velocity fluctuations, Reynolds shear stress, and low–speed sreak spacing) is small. The riblests cause a relatively small drag reduction of about 4%, a figure that is in rough agreement with experiments and other computations. The simulations also suggest a mechanism for the observed displacement of the turbulence quantities away from the wall.

The immersed boundary technique used to model the riblets consists of creating an externally imposed spatially localized body force which opposes the flow velocity and creates a riblet-like surface. For unstead viscous flow the calculation of the force is done with a feedback scheme in which the velocity is used to iteratively determine the desired value. In particular, the surface body force is determined by the relation f(xs, t) = α ∫ t0U(xs,t′)dt′ + βU(xs, t) for surface points xs, velocity U time t and negative constants α and β. All simulations are done with a spectral code in a single computational domain without any mapping of the mesh. The combination of the immersed boundary and spectral techniques can potentially be used to solve other problems having complex geometry and flow physics.