The stability properties of 24 experimentally generated internal solitary waves (ISWs) of extremely large amplitude, all with minimum Richardson number less than 1/4, are investigated. The study is supplemented by fully nonlinear calculations in a three-layer fluid. The waves move along a linearly stratified pycnocline (depth h2) sandwiched between a thin upper layer (depth h1) and a deep lower layer (depth h3), both homogeneous. In particular, the wave-induced velocity profile through the pycnocline is measured by particle image velocimetry (PIV) and obtained in computation. Breaking ISWs were found to have amplitudes (a1) in the range , while stable waves were on or below this limit. Breaking ISWs were investigated for 0.27 < h2/h1 < 1 and 4.14 < h3/(h1 + h2) < 7.14 and stable waves for 0.36 < h2/h1 < 3.67 and 3.22 < h3/(h1 + h2) < 7.25. Kelvin–Helmholtz-like billows were observed in the breaking cases. They had a length of 7.9h2 and a propagation speed 0.09 times the wave speed. These measured values compared well with predicted values from a stability analysis, assuming steady shear flow with U(z) and ρ(z) taken at the wave maximum (U(z) horizontal velocity profile, ρ(z) density along the vertical z). Only unstable modes in waves of sufficient strength have the chance to grow sufficiently fast to develop breaking: the waves that broke had an estimated growth (of unstable modes) more than 3.3–3.7 times than in the strongest stable case. Evaluation of the minimum Richardson number (Rimin, in the pycnocline), the horizontal length of a pocket of possible instability, with wave-induced Ri < 14, (Lx) and the wavelength (λ), showed that all measurements fall within the range Rimin = −0.23Lx/λ + 0.298 ± 0.016 in the (Lx/λ, Rimin)-plane. Breaking ISWs were found for Lx/λ > 0.86 and stable waves for Lx/λ < 0.86. The results show a sort of threshold-like behaviour in terms of Lx/λ. The results demonstrate that the breaking threshold of Lx/λ = 0.86 was sharper than one based on a minimum Richardson number and reveal that the Richardson number was found to become almost antisymmetric across relatively thick pycnoclines, with the minimum occurring towards the top part of the pycnocline.

Direct numerical simulations (DNSs) are used to investigate the drag-reducing performance of superhydrophobic surfaces (SHSs) in turbulent channel flow. SHSs combine surface roughness with hydrophobicity and can, in some cases, support a shear-free air–water interface. Slip velocities, wall shear stresses and Reynolds stresses are considered for a variety of SHS microfeature geometry configurations at a friction Reynolds number of Reτ ≈ 180. For the largest microfeature spacing studied, an average slip velocity over 75% of the bulk velocity is obtained, and the wall shear stress reduction is found to be nearly 40%. The simulation results suggest that the mean velocity profile near the superhydrophobic wall continues to scale with the wall shear stress but is offset by a slip velocity that increases with increasing microfeature spacing.

In this paper, the regular reflection (RR) to Mach reflection (MR) transition of asymmetric shock waves is theoretically studied by employing the classical two- and three-shock theories. Computations are conducted to evaluate the effects of expansion fans, which are inherent flow structures in asymmetric reflection of shock waves, on the RR → MR transition. Comparison shows good agreement among the theoretical, numerical and experimental results. Some discrepancies between experiment and theory reported in previous studies are also explained based on the present theoretical analysis. The advanced RR → MR transition triggered by a transverse wave is also discussed for the interaction of a hypersonic flow and a double-wedge-like geometry.

We present an experimental investigation of the flow structure and vorticity field in the wake of a NACA-0012 airfoil pitching sinusoidally at small amplitude and high reduced frequencies. Molecular tagging velocimetry is used to quantify the characteristics of the vortex array (circulation, peak vorticity, core size, spatial arrangement) and its downstream evolution over the first chord length as a function of reduced frequency. The measured mean and fluctuating velocity fields are used to estimate the mean force on the airfoil and explore the connection between flow structure and thrust generation.

Results show that strong concentrated vortices form very rapidly within the first wavelength of oscillation and exhibit interesting dynamics that depend on oscillation frequency. With increasing reduced frequency the transverse alignment of the vortex array changes from an orientation corresponding to velocity deficit (wake profile) to one with velocity excess (reverse Kármán street with jet profile). It is found, however, that the switch in the vortex array orientation does not coincide with the condition for crossover from drag to thrust. The mean force is estimated from a more complete control volume analysis, which takes into account the streamwise velocity fluctuations and the pressure term. Results clearly show that neglecting these terms can lead to a large overestimation of the mean force in strongly fluctuating velocity fields that are characteristic of airfoils executing highly unsteady motions. Our measurements show a decrease in the peak vorticity, as the vortices convect downstream, by an amount that is more than can be attributed to viscous diffusion. It is found that the presence of small levels of axial velocity gradients within the vortex cores, levels that can be difficult to measure experimentally, can lead to a measurable decrease in the peak vorticity even at the centre of the flow facility in a flow that is expected to be primarily two-dimensional.

The steady two-dimensional laminar flow around a stationary circular cylinder has been investigated via a stabilized finite-element method. The Reynolds number Re is based on the cylinder diameter and free-stream speed. The results have been presented for 6 ≤ Re ≤ 40 and the blockages between 0.000125 and 0.80. The blockage B is the ratio of the cylinder diameter to the domain width. There is large scatter in the value of Res, reported in the literature, marking the onset of the flow separation. From the present study the Res is found to be 6.29, approximately for B = 0.005. The effect of the blockage on the characteristic flow parameters is found to be insignificant for B ≤ 0.01. The bubble length, separation angle and Res exhibit non-monotonic variation with the blockage. It is for the first time that such a behaviour for the separation angle and Res is being reported. Two types of boundary conditions at the lateral walls have been studied: the slip wall and towing tank. In general for high blockage, the results from the slip boundary condition are closer to the ones for the unbounded flow. In that sense, the use of the slip boundary condition as opposed to the towing tank boundary condition on the lateral walls is advocated. The bubble length, separation angle, base suction, total drag, pressure drag, viscous drag and maximum vorticity on the cylinder surface for the steady flow are found to vary as Re, Re−0.5, Re−1, Re−0.5, Re−0.64, Re−0.60 and Re0.5, respectively. The extrapolated results for the steady flow, for higher Re, are found to match quite well with the other results from the literature.

We perform a multi-scale non-local geometrical analysis of the structures extracted from the enstrophy and kinetic energy dissipation-rate, instantaneous fields of a numerical database of incompressible homogeneous isotropic turbulence decaying in time obtained by DNS in a periodic box. Three different resolutions are considered: 2563, 5123 and 10243 grid points, with kmax approximately 1, 2 and 4, respectively, the same initial conditions and Reλ ≈ 77. This allows a comparison of the geometry of the structures obtained for different resolutions. For the highest resolution, structures of enstrophy and dissipation evolve in a continuous distribution from blob-like and moderately stretched tube-like shapes at the large scales to highly stretched sheet-like structures at the small scales. The intermediate scales show a predominance of tube-like structures for both fields, much more pronounced for the enstrophy field. The dissipation field shows a tendency towards structures with lower curvedness than those of the enstrophy, for intermediate and small scales. The 2563 grid resolution case (kmax ≈ 1) was unable to detect the predominance of highly stretched sheet-like structures at the smaller scales in both fields. The same non-local methodology for the study of the geometry of structures, but without the multi-scale decomposition, is applied to two scalar fields used by existing local criteria for the eduction of tube- and sheet-like structures in turbulence, Q and [Aij]+, respectively, obtained from invariants of the velocity-gradient tensor and alike in the 10243 case. This adds the non-local geometrical characterization and classification to those local criteria, assessing their validity in educing particular geometries. Finally, we introduce a new methodology for the study of proximity issues among structures of different fields, based on geometrical considerations and non-local analysis, by taking into account the spatial extent of the structures. We apply it to the four fields previously studied. Tube-like structures of Q are predominantly surrounded by sheet-like structures of [Aij]+, which appear at closer distances. For the enstrophy, tube-like structures at an intermediate scale are primarily surrounded by sheets of smaller scales of the enstrophy and structures of dissipation at the same and smaller scales. A secondary contribution results from tubes of enstrophy at smaller scales appearing at farther distances. Different configurations of composite structures are presented.

Measurements of the bow waves generated by a rectangular flat plate, immersed at a draught D = 0.2 m, towed at constant speed U = 1.75 m s−1 in calm water and held at a heel angle 10° and a series of nine yaw angles α = 10°, 15°, 20°, 25°, 30°, 45°, 60°, 75° and 90° are reported. The measurements show that bow wave unsteadiness is significantly larger for the flat plate towed at yaw angles 30° ≤ α ≤ 90° than at 10° ≤ α ≤ 20°, which are associated with the unsteady and overturning bow wave regimes, respectively, separated by the boundary with g ≡ acceleration of gravity. These measurements of bow wave unsteadiness provide preliminary experimental validation of the foregoing simple theoretical relation for the boundary between the unsteady and overturning bow wave regimes for non-bulbous wedge-shaped ship bows with insignificant rake and flare. Extension of this relation to more complicated ship bows, notably bows with rake and flare, is also considered.

Along the lines of the restricted Euler equation (REE) for incompressible flows, we develop homogenized Euler equation (HEE) for describing turbulent velocity gradient dynamics of an isentropic compressible calorically perfect gas. Starting from energy and state equations, an evolution equation for pressure Hessian is derived invoking uniform (homogeneous) velocity gradient assumption. Behaviour of principal strain rates, vorticity vector alignment and invariants of the normalized velocity gradient tensor is investigated conditioned on dilatation level. The HEE results agree very well with the known behaviour in the incompressible limit. Indeed, at zero dilatation HEE reproduces the incompressible anisotropic pressure Hessian behaviour very closely. When compared against compressible direct numerical simulation results, the HEE accurately captures the strain rate behaviour at different dilatation levels. The model also recovers the fixed point behaviour of pressure-released (high-Mach-number limit) Burgers turbulence.

Three-dimensional numerical simulations of laminar flow around a circular cylinder with sinusoidal variation of cross-section along the spanwise direction, named ‘wavy cylinder’, are performed. A series of wavy cylinders with different combinations of dimensionless wavelength (λ/Dm) and wave amplitude (a/Dm) are studied in detail at a Reynolds number of Re = U∞Dm/ν = 100, where U∞ is the free-stream velocity and Dm is the mean diameter of a wavy cylinder. The results of variation of mean drag coefficient and root mean square (r.m.s.) lift coefficient with dimensionless wavelength show that significant reduction of mean and fluctuating force coefficients occurs at optimal dimensionless wavelengths λ/Dm of around 2.5 and 6 respectively for the different amplitudes studied. Based on the variation of flow structures and force characteristics, the dimensionless wavelength from λ/Dm = 1 to λ/Dm = 10 is classified into three wavelength regimes corresponding to three types of wake structures. The wake structures at the near wake of different wavy cylinders are captured. For all wavy cylinders, the flow separation line varies along the spanwise direction. This leads to the development of a three-dimensional free shear layer with periodic repetition along the spanwise direction. The three-dimensional free shear layer of the wavy cylinder is larger and more stable than that of the circular cylinder, and in some cases the free shear layer even does not roll up into a mature vortex street behind the cylinder. As a result, the mean drag coefficients of some of the typical wavy cylinders are less than that of a corresponding circular cylinder with a maximum drag coefficient reduction up to 18%. The r.m.s. lift coefficients are greatly reduced to practically zero at optimal wavelengths. In the laminar flow regime (60 ≤ Re ≤ 150), the values of optimal wavelength are Reynolds number dependent.

We present a second-order stochastic model of weakly nonlinear waves and develop theoretical expressions for the expected shape of large surface displacements. The model also leads to an exact theoretical expression for the statistical distribution of large wave crests in a form that generalizes the Tayfun distribution (Tayfun, J. Geophys. Res., vol. 85, 1980, p. 1548). The generalized distribution depends on a steepness parameter given by μ = λ3/3, where λ3 represents the skewness coefficient of surface displacements. It converges to the Tayfun distribution in narrowband waves, where both distributions describe the crests of all waves well. In broadband waves, the generalized distribution represents the crests of large waves just as well whereas the Tayfun distribution appears as an upper bound and tends to overestimate them. However, the theoretical nature of the generalized distribution presents practical difficulties in oceanic applications. We circumvent these by adopting an appropriate approximation for the steepness parameter. Comparisons with wind-wave measurements from the North Sea suggest that this approximation allows both distributions to assume an identical form with which we can describe the distribution of large wave crests fairly accurately. The same comparisons also show that third-order nonlinear effects do not appear to have any discernable effect on the statistics of large surface displacements or wave crests.

We have studied the temporal evolution of electro-kinetic flows in the vicinity of polarizable dielectric solids following the application of a ‘weak’ transient electric field. To obtain a macro-scale description in the limit of narrow electric double layers (EDLs), we have derived a pair of effective transient boundary conditions directly connecting the electric potentials across the EDL. Within the framework of the above assumptions, these conditions apply to a general transient electro-kinetic problem involving dielectric solids of arbitrary geometry and relative permittivity. Furthermore, the newly derived scheme is applicable to general transient and spatially non-uniform external fields. We examine the details of the physical mechanisms involved in the relaxation of the induced-charging process of the EDL adjacent to polarizable dielectric solids. It is thus established that the time scale characterizing the electrostatic relaxation increases with the dielectric constant of the solid from the Debye time (for the diffusion across the EDL) through the ‘intermediate’ scale (proportional to the product of the respective Debye- and geometric-length scales). Thus, the present rigorous analysis substantiates earlier results largely obtained by heuristic use of equivalent RC-circuit models. Furthermore, for typical values of ionic diffusivity and kinematic viscosity of the electrolyte solution, the latter time scale is comparable to the time scale of viscous relaxation in problems concerning microfluidic applications or micro-particle dynamics. The analysis is illustrated for spherical micro-particles. Explicit results are thus presented for the temporal evolution of electro-osmosis around a dielectric sphere immersed in unbounded electrolyte solution under the action of a suddenly applied uniform field, combining both induced charge and ‘equilibrium’ (fixed charge) contributions to the zeta potential. It is demonstrated that, owing to the time delay of the induced-EDL charging, the ‘equilibrium’ contribution to fluid motion (which is linear in the electric field) initially dominates the (quadratic) ‘induced’ contribution.

The dynamics and control of two-dimensional disturbances in the spatially evolving boundary layer on a flat plate are investigated from an input–output viewpoint. A set-up of spatially localized inputs (external disturbances and actuators) and outputs (objective functions and sensors) is introduced for the control design of convectively unstable flow configurations. From the linearized Navier–Stokes equations with the inputs and outputs, controllable, observable and balanced modes are extracted using the snapshot method. A balanced reduced-order model (ROM) is constructed and shown to capture the input–output behaviour of the linearized Navier–Stokes equations. This model is finally used to design a 2-feedback controller to suppress the growth of two-dimensional perturbations inside the boundary layer.

It is numerically known that normal modes for small disturbances to pipe Poiseuille flow have complex phase velocities c which form a Y-shaped set of discrete points in the fourth quadrant, when the Reynolds number R is large. In this paper, the eigenvalue problem of determining c for axisymmetric torsional disturbances is treated, and the Y-shaped distribution of these c is studied analytically (with a little aid from numerics) by using some asymptotic forms of a Whittaker function. As a result, a Y-shaped contour on which eigenvalues c are approximately located is obtained, independent of the wavenumber α and R, from simple equations which contain elementary functions only. Naturally, the location of each individual c depends on α and R. How it changes on the contour when large R is becoming still larger is explained. Several approximate values of c on the contour are compared with c computed by Schmid & Henningson (1994), and their agreement is seen to be good. Furthermore, the limit Re c → 2/3 as Im c → −∞ is rigorously proved when R is fixed at an arbitrary number, which is not required to be large. This limit is shown to be true also with eigenvalues c for axisymmetric meridional disturbances. The alternate distribution of c for torsional and meridional disturbances on a branch of the Y-shaped contour is explained.

We present a theoretical model for the wave-induced drift and horizontal deformation of an oil slick. The waves and the mean flow are coupled through the influence of the mean flow on the concentration of slick material, which in turn determines the damping rate of the waves and hence the transfer of momentum from the waves to the mean flow. We also briefly discuss a simplified version of the model that can be used when remote sensing data are available. With this simpler model the wave-induced forcing of the mean flow is obtained directly from observations of the wave field, hence knowledge of any specific slick properties is not required.

Drop coalescence was studied using numerical simulations. Liquid drops were made to coalesce with a body of the same liquid, either a reservoir or a drop of different size, each with negligible impact velocity. We considered either gas or liquid as a surrounding fluid, and experimental results are discussed for the gas–liquid set-up. Under certain conditions, a drop will not fully coalesce with the liquid reservoir, leaving behind a daughter drop. Partial coalescence is observed for systems of low viscosity, characterized by a small Ohnesorge number, where capillary waves remain sufficiently vigourous to distort the drop significantly. For drops coalescing with a flat interface, we determine the critical Ohnesorge number as a function of Bond number, as well as density and viscosity ratios of the fluids. Studying the coalescence of two drops of different sizes reveals that partial coalescence may occur in low-viscosity systems provided the size ratio of the drops exceeds a certain threshold. We also determine the extent to which the process of partial coalescence is self-similar and find that the viscosity of the drop has a large effect on the droplet's vertical velocity after pinch off. Finally, we report on the formation of satellite droplets generated after a first pinch off and on the ejection of a jet of tiny droplets during coalescence of a parent drop significantly deformed by gravity.

Three-dimensional and time-dependent simulations of viscoelastic Taylor–Couette flow of dilute polymer solutions are performed using a fully implicit parallel spectral time-splitting algorithm to discover flow patterns with various spatio-temporal symmetries, namely rotating standing waves (RSWs), disordered oscillations (DOs) and solitary vortex structures referred to as oscillatory strips (OSs) and diwhirls (DWs). A detailed account of the impact of flow transitions on molecular conformation and viscoelastic stress, velocity profiles, hydrodynamic drag force and energy spectra of time-dependent flow states is presented. Overall, predicted pattern selection and flow features compare very favourably with experimental observations. For elasticity number E, that signifies the ratio of elastic to viscous forces, >0.1, and when the shear rate (cylinder rotation speed) is increased above the linear stability threshold, the circular Couette flow (CCF) becomes unstable to RSWs which are characterized by a checkerboard-like pattern in the space–time plot of radial velocity, implying symmetry between inflow/outflow (I/O) regions. As the shear rate is further increased, perturbations that break the I/O symmetry are amplified leading to DOs and/or flame-like patterns with spectral mechanical energy transfer reminiscent of elastically induced low-Reynolds-number turbulence. However, when the shear rate is decreased from those at which such chaotic states are observed, the radially inward acting polymer body force created by flow-induced molecular stretching causes the development of narrow inflow regions surrounded by much broader weak outflow domains. This promotes the formation of solitary vortex structures, which can be stationary and axisymmetric (DWs) or time-dependent (OSs). The dynamics of the formation of these structures by merging and coalescence of vortex pairs and the implication of such events on instantaneous hydrodynamic force are studied. For O(1) values of E, OSs and DWs appear approximately at constant values of the We, defined as the ratio of polymer relaxation time to the inverse shear rate in the gap. As shear rate is decreased further, DWs decay to CCF although at We values less than the linear stability threshold. The flow transitions are hysteretic with respect to We, as evidenced by a plot of drag force versus We.

In the present study, the effects of the jet inflow conditions such as the initial momentum thickness (θ) and background disturbances on the downstream evolution of a circular jet are investigated using large eddy simulation (LES). We consider four different initial momentum thicknesses, D/θ = 50, 80, 120 and 180, and three different Reynolds numbers, ReD = UJD/ν = 3600, 104 and 105, where UJ is the jet inflow velocity and D is the jet diameter. The present study shows that the jet characteristics significantly depend on both the initial momentum thickness and the Reynolds number. For all the Reynolds numbers considered in this study, vortex rings are generated at an earlier position with decreasing initial momentum thickness. In case of a relatively low Reynolds number of ReD = 3600, at smaller initial momentum thickness, early growth of the shear layer due to the early generation of vortex ring leads to the occurrence of large-scale coherent structures in earlier downstream locations, which results in larger mixing enhancement and more rapid increase in turbulence intensity. However, at a high Reynolds number such as ReD = 105, with decreasing initial momentum thickness, rapid growth of the shear layer leads to the saturation of the shear layer and the generation of fine-scale turbulence structures, which reduces mixing and turbulence intensity. With increasing Reθ (= UJ θ/ν), the characteristic frequency corresponding to the shear layer mode (Stθ = fθ/UJ) gradually increases and reaches near 0.017 predicted from the inviscid instability theory. On the other hand, the frequency corresponding to the jet-preferred mode (StD = f D/UJ) varies depending on ReD and D/θ. From a mode analysis, we show that, in view of the energy of the axial velocity fluctuations integrated over the radial direction, the double-helix mode (mode 2) becomes dominant past the potential core, but the axisymmetric mode (mode 0) is dominant near the jet exit. In view of the local energy, the disturbances grow along the shear layer near the jet exit: for thick shear layer, mode 0 grows much faster than other modes, but modes 0–3 grow almost simultaneously for thin shear layer. However, past the potential core, the dominant mode changes from mode 0 near the centreline to mode 1 and then to mode 2 with increasing radial direction regardless of the initial shear layer thickness.