Studying particle-laden oscillatory channel flow constitutes an important step towards understanding practical application. This study aims to take a step forward in our understanding of the role of turbulence on fine-particle transport in an oscillatory channel and the back effect of fine particles on turbulence modulation using an Eulerian–Eulerian framework. In particular, simulations presented in this study are selected to investigate wave-induced fine sediment transport processes in a typical coastal setting. Our modelling framework is based on a simplified two-way coupled formulation that is accurate for particles of small Stokes number (St). As a first step, the instantaneous particle velocity is calculated as the superposition of the local fluid velocity and the particle settling velocity while the higher-order particle inertia effect neglected. Correspondingly, only the modulation of carrier flow is due to particle-induced density stratification quantified by the bulk Richardson number, Ri. In this paper, we fixed the Reynolds number to be ReΔ = 1000 and varied the bulk Richardson number over a range (Ri = 0, 1 × 10−4, 3 × 10−4 and 6 × 10−4). The simulation results reveal critical processes due to different degrees of the particle–turbulence interaction. Essentially, four different regimes of particle transport for the given ReΔ are observed: (i) the regime where virtually no turbulence modulation in the case of very dilute condition, i.e. Ri ~ 0; (ii) slightly modified regime where slight turbulence attenuation is observed near the top of the oscillatory boundary layer. However, in this regime a significant change can be observed in the concentration profile with the formation of a lutocline; (iii) regime where flow laminarization occurs during the peak flow, followed by shear instability during the flow reversal. A significant reduction in the oscillatory boundary layer thickness is also observed; (iv) complete laminarization due to strong particle-induced stable density stratification.

The spatiotemporal evolution of a viscoelastic jet depends on the relative magnitude of capillary, viscous, inertial and elastic stresses. The interplay of capillary and elastic stresses leads to the formation of very thin and stable filaments between drops, or to ‘beads-on-a-string’ structure. In this paper, we show that by understanding the physical processes that control different stages of the jet evolution it is possible to extract transient extensional viscosity information even for very low viscosity and weakly elastic liquids, which is a particular challenge in using traditional rheometers. The parameter space at which a forced jet can be used as an extensional rheometer is numerically investigated by using a one-dimensional nonlinear free-surface theory for Oldroyd-B and Giesekus fluids. The results show that even when the ratio of viscous to inertio-capillary time scales (or Ohnesorge number) is as low as Oh ~ 0.02, the temporal evolution of the jet can be used to obtain elongational properties of the liquid.

The flow through a compressor passage without and with incoming free-stream grid turbulence is simulated. At moderate Reynolds number, laminar-to-turbulence transition can take place on both sides of the aerofoil, but proceeds in distinctly different manners. The direct numerical simulations (DNS) of this flow reveal the mechanics of breakdown to turbulence on both surfaces of the blade. The pressure surface boundary layer undergoes laminar separation in the absence of free-stream disturbances. When exposed to free-stream forcing, the boundary layer remains attached due to transition to turbulence upstream of the laminar separation point. Three types of breakdowns are observed; they combine characteristics of natural and bypass transition. In particular, instability waves, which trace back to discrete modes of the base flow, can be observed, but their development is not independent of the Klebanoff distortions that are caused by free-stream turbulent forcing. At a higher turbulence intensity, the transition mechanism shifts to a purely bypass scenario. Unlike the pressure side, the suction surface boundary layer separates independent of the free-stream condition, be it laminar or a moderate free-stream turbulence of intensity Tu ~ 3%. Upstream of the separation, the amplification of the Klebanoff distortions is suppressed in the favourable pressure gradient (FPG) region. This suppression is in agreement with simulations of constant pressure gradient boundary layers. FPG is normally stabilizing with respect to bypass transition to turbulence, but is, thereby, unfavourable with respect to separation. Downstream of the FPG section, a strong adverse pressure gradient (APG) on the suction surface of the blade causes the laminar boundary layer to separate. The separation surface is modulated in the instantaneous fields of the Klebanoff distortion inside the shear layer, which consists of forward and backward jet-like perturbations. Separation is followed by breakdown to turbulence and reattachment. As the free-stream turbulence intensity is increased, Tu ~ 6.5%, transitional turbulent patches are initiated, and interact with the downstream separated flow, causing local attachment. The calming effect, or delayed re-establishment of the boundary layer separation, is observed in the wake of the turbulent events.

Streamwise and quasi-streamwise elongated structures have been shown to play a significant role in turbulent shear flows. We model the mean behaviour of fully turbulent plane Couette flow using a streamwise constant projection of the Navier–Stokes equations. This results in a two-dimensional three-velocity-component (2D/3C) model. We first use a steady-state version of the model to demonstrate that its nonlinear coupling provides the mathematical mechanism that shapes the turbulent velocity profile. Simulations of the 2D/3C model under small-amplitude Gaussian forcing of the cross-stream components are compared to direct numerical simulation (DNS) data. The results indicate that a streamwise constant projection of the Navier–Stokes equations captures salient features of fully turbulent plane Couette flow at low Reynolds numbers. A systems-theoretic approach is used to demonstrate the presence of large input–output amplification through the forced 2D/3C model. It is this amplification coupled with the appropriate nonlinearity that enables the 2D/3C model to generate turbulent behaviour under the small-amplitude forcing employed in this study.

Flow control induced by a single wire that is attached on the outer surface and parallel to the span of a stationary circular cylinder is investigated experimentally. The Reynolds number has a value of 10 000 and the wire diameter is nearly two orders of magnitude smaller than the cylinder diameter, while being larger than the thickness of the unperturbed boundary layer forming around the cylinder. A technique of high-image-density particle image velocimetry is used to characterize mean and unsteady structures of the separating shear layer and the near wake. Only one of the shear layers is directly perturbed by the surface wire. This disturbance, however, has important global consequences over the entire near wake, provided that the wire is located within a certain range of angular positions with respect to the approach flow. Over this range, there are two angles that can be defined as critical on the basis of the streamwise extent of the near-wake structure. In a simplified sense, these critical angles are associated with significant extension and contraction of the near wake, relative to the wake in the absence of the effect of a surface disturbance. The critical angle of the wire that yields the most significant extension of the near wake is also found to lead to bistable oscillations of the separating shear layer at irregular time intervals, much longer than the time scale associated with the classical Kármán vortex shedding. The foregoing two critical states of extension and contraction of the near wake are, respectively, linked to attenuation or enhancement of the Kármán instability. Moreover, the onset of the shear-layer instability, Reynolds stress, Strouhal number and the transverse extent of shear-layer flapping are all shown to depend on the angular position of the wire within the defined range of angles.

We investigate the influence of the thermal properties of the boundaries in turbulent Rayleigh–Bénard convection on analytical upper bounds on convective heat transport. We model imperfectly conducting bounding plates in two ways: using idealized mixed thermal boundary conditions (BCs) of constant Biot number η, continuously interpolating between the previously studied fixed temperature (η = 0) and fixed flux (η = ∞) cases; and by explicitly coupling the evolution equations in the fluid in the Boussinesq approximation through temperature and flux continuity to identical upper and lower conducting plates. In both cases, we systematically formulate a bounding principle and obtain explicit upper bounds on the Nusselt number Nu in terms of the usual Rayleigh number Ra measuring the average temperature drop across the fluid layer, using the ‘background method’ developed by Doering and Constantin. In the presence of plates, we find that the bounds depend on σ = d/λ, where d is the ratio of plate to fluid thickness and λ is the conductivity ratio, and that the bounding problem may be mapped onto that for Biot number η = σ. In particular, for each σ > 0, for sufficiently large Ra (depending on σ) we show that Nu ≤ c(σ) R1/3 ≤ CRa1/2, where C is a σ-independent constant, and where the control parameter R is a Rayleigh number defined in terms of the full temperature drop across the entire plate–fluid–plate system. In the Ra → ∞ limit, the usual fixed temperature assumption is a singular limit of the general bounding problem, while fixed flux conditions appear to be most relevant to the asymptotic Nu–Ra scaling even for highly conducting plates.

An exhaustive parametric study of the transition scenario in the wake of oblate spheroids and flat cylinders placed with their rotation axis parallel to the flow is presented. The flatness of the investigated objects is classified by the aspect ratio χ defined as χ = d/a for spheroids (with d the diameter and a the length of the polar axis) and as χ = d/h) for cylinders (with h the cylinder height). We find a significant qualitative similarity between both configurations. At large aspect ratios (χ > 2.3 for spheroids and χ ≥ 4 for cylinders), the secondary bifurcation giving rise to a periodic state without planar symmetry is subcritical with a hysteresis interval of about two Reynolds number units. For spheroids, the sphere-like scenario is recovered only at aspect ratios very close to one (χ ≥ 1 are considered), while for cylindrical bodies the same holds for χ ≤ 1.7. For intermediate aspect ratios, a domain of states with non-zero net helicity separates states typical for the sphere wake from those of an infinitely flat disk.

This paper is the second part of the work on linear and nonlinear stability of buoyancy-driven coastal currents. Part 1, concerning a passive lower layer, was presented in the companion paper Gula & Zeitlin (J. Fluid Mech., vol. 659, 2010, p. 69). In this part, we use a fully baroclinic two-layer model, with active lower layer. We revisit the linear stability problem for coastal currents and study the nonlinear evolution of the instabilities with the help of high-resolution direct numerical simulations. We show how nonlinear saturation of the ageostrophic instabilities leads to reorganization of the mean flow and emergence of coherent vortices. We follow the same lines as in Part 1 and, first, perform a complete linear stability analysis of the baroclinic coastal currents for various depths and density ratios. We then study the nonlinear evolution of the unstable modes with the help of the recent efficient two-layer generalization of the one-layer well-balanced finite-volume scheme for rotating shallow water equations, which allows the treatment of outcropping and loss of hyperbolicity associated with shear, Kelvin–Helmholtz type, instabilities. The previous single-layer results are recovered in the limit of large depth ratios. For depth ratios of order one, new baroclinic long-wave instabilities come into play due to the resonances among Rossby and frontal- or coastal-trapped waves. These instabilities saturate by forming coherent baroclinic vortices, and lead to a complete reorganization of the initial current. As in Part 1, Kelvin fronts play an important role in this process. For even smaller depth ratios, short-wave shear instabilities with large growth rates rapidly develop. We show that at the nonlinear stage they produce short-wave meanders with enhanced dissipation. However, they do not change, globally, the structure of the mean flow which undergoes secondary large-scale instabilities leading to coherent vortex formation and cutoff.

Numerical investigation of the compressible flow past a wavy cylinder was carried out using large-eddy simulation for a free-stream Mach number M∞ = 0.75 and a Reynolds number based on the mean diameter Re = 2 × 105. The flow past a corresponding circular cylinder was also calculated for comparison and validation against experimental data. Various fundamental mechanisms dictating the intricate flow phenomena, including drag reduction and fluctuating force suppression, shock and shocklet elimination, and three-dimensional separation and separated shear-layer instability, have been studied systematically. Because of the passive control of the flow over a wavy cylinder, the mean drag coefficient of the wavy cylinder is less than that of the circular cylinder with a drag reduction up to 26%, and the fluctuating force coefficients are significantly suppressed to be nearly zero. The vortical structures near the base region of the wavy cylinder are much less vigorous than those of the circular cylinder. The three-dimensional shear-layer shed from the wavy cylinder is more stable than that from the circular cylinder. The vortex roll up of the shear layer from the wavy cylinder is delayed to a further downstream location, leading to a higher-base-pressure distribution. The spanwise pressure gradient and the baroclinic effect play an important role in generating an oblique vortical perturbation at the separated shear layer, which may moderate the increase of the fluctuations at the shear layer and reduce the growth rate of the shear layer. The analysis of the convective Mach number indicates that the instability processes in the shear-layer evolution are derived from oblique modes and bi-dimensional instability modes and their competition. The two-layer structures of the shear layer are captured using the instantaneous Lamb vector divergence, and the underlying dynamical processes associated with the drag reduction are clarified. Moreover, some phenomena relevant to the compressible effect, such as shock waves, shocklets and shock/turbulence interaction, are analysed. It is found that the shocks and shocklets which exist in the circular cylinder flow are eliminated for the wavy cylinder flow and the wavy surface provides an effective way of shock control. As the shock/turbulence interaction is avoided, a significant drop of the turbulent fluctuations around the wavy cylinder occurs. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the passive control of the compressible flow past a wavy surface.

The evolution and dynamics of a shallow-water vortex system with high initial Reynolds numbers are investigated experimentally without background rotation. A single vortex is generated by rotating a water mass at the centre of an experimental tank using a bottomless cylinder with internal sectors. The surface velocity field is observed via particle image velocimetry. The experimentally observed vorticity fields indicate that strong shallowness (the ratio of the cylinder diameter to the water depth) and high Reynolds number contribute to the formation of large-scale coherent structures in the form of a tripolar vortex system. The shallow-water vortices with high initial Reynolds numbers experience the transition from turbulent to laminar regimes in their decay process. The proposed first-order vortex decay model predicts that a shallow-water vortex decays as t−1 in the initial turbulent stage and as e−t in the later laminar stage due to horizontal diffusion and bottom friction. The estimated transition time scale from the turbulent to laminar stage increases with initial vortex Reynolds number and with shallowness. By taking the vortex expansion into consideration, the second-order vortex decay model is also presented. The azimuthally ensemble-averaged data elucidate effects of the vortex instabilities and of turbulent energy transfer on the formation of large-scale coherent flow structures. Normal mode analysis of the vortex systems is conducted to study the effect of shallowness and Reynolds number on the generation of two-dimensional large-scale coherent structures. The results show that the perturbation wavenumber of mode 2 is the fastest-growing instability in shallow-water conditions, and its effect depends on initial Reynolds number and shallowness.

Measurements of the Nusselt number Nu and of properties of the large-scale circulation (LSC) for turbulent Rayleigh–Bénard convection are presented in the presence of rotation about a vertical axis at angular speeds 0 ≤ Ω ≲ 2 rad s−1. The sample chamber was cylindrical with a height equal to the diameter, and the fluid contained in it was water. The LSC was studied by measuring sidewall temperatures as a function of azimuthal position. The measurements covered the Rayleigh-number range 3 × 108 ≲ Ra ≲ 2 × 1010, the Prandtl-number range 3.0 ≲ Pr ≲ 6.4 and the Rossby-number range 0 ≤ (1/Ro ∝ Ω) ≲ 20. At modest 1/Ro, we found an enhancement of Nu due to Ekman-vortex pumping by as much as 20%. As 1/Ro increased from zero, this enhancement set in discontinuously at and grew above 1/Roc. The value of 1/Roc varied from about 0.48 at Pr = 3 to about 0.35 at Pr = 6.2. At sufficiently large 1/Ro (large rotation rates), Nu decreased again, due to the Taylor–Proudman (TP) effect, and reached values well below its value without rotation. The maximum enhancement increased with increasing Pr and decreasing Ra and, we believe, was determined by a competition between the Ekman enhancement and the TP depression. The temperature signature along the sidewall of the LSC was detectable by our method up to 1/Ro ≃ 1. The frequency of cessations α of the LSC grew dramatically with increasing 1/Ro, from about 10−5 s−1 at 1/Ro = 0 to about 2 × 10−4 s−1 at 1/Ro = 0.25. A discontinuous further increase of α, by about a factor of 2.5, occurred at 1/Roc. With increasing 1/Ro, the time-averaged and azimuthally averaged vertical thermal gradient along the sidewall first decreased and then increased again, with a minimum somewhat below 1/Roc. The Reynolds number of the LSC, determined from oscillations of the time correlation functions of the sidewall temperatures, was constant within our resolution for 1/Ro ≲ 0.3 and then decreased with increasing 1/Ro. The retrograde rotation rate of the LSC circulation plane exhibited complex behaviour as a function of 1/Ro even at small rotation rates corresponding to 1/Ro < 1/Roc.

In order to investigate the turbulent drag reduction phenomenon and understand its mechanism, direct numerical simulation (DNS) was carried out on decaying homogeneous isotropic turbulence (DHIT) with and without polymer additives. We explored the polymer effect on DHIT from the energetic viewpoint, i.e. the decay of the total turbulent kinetic energy and energy distribution at each scale in Fourier space and from the phenomenological viewpoint, i.e. the alterations of vortex structures, the enstrophy and the strain. It was obtained that in DHIT with polymer additives the decay of the turbulent kinetic energy is faster than that in the Newtonian fluid case and a modification of the turbulent kinetic energy transfer process for the Newtonian fluid flow is observed due to the release of the polymer elastic energy into flow structures at certain small scales. Besides, we deduced the transport equations of the enstrophy and the strain, respectively, for DHIT with polymer additives. Based on the analyses of these transport equations, it was found that polymer additives depress both the enstrophy and the strain in DHIT as compared to the Newtonian fluid case, indicating the inhibition effect on small-scale vortex structures and turbulence intensity by polymers.

Smooth flat-plate turbulent boundary layers (TBLs) have been studied for nearly a century. However, there is a relative dearth of measurements at Reynolds numbers typical of full-scale marine and aerospace transportation systems (Reθ = Ueθ/ν > 105, where Ue = free-stream speed, θ = TBL momentum thickness and ν = kinematic viscosity). This paper presents new experimental results for the TBL that forms on a smooth flat plate at nominal Reθ values of 0.5 × 105, 1.0 × 105 and 1.5 × 105. Nominal boundary layer thicknesses (δ) were 80–90mm, and Karman numbers (δ+) were 17000, 32000 and 47000, respectively. The experiments were conducted in the William B. Morgan Large Cavitation Channel on a polished (k+ < 0.2) flat-plate test model 12.9m long and 3.05m wide at water flow speeds up to 20ms−1. Direct measurements of static pressure and mean wall shear stress were obtained with pressure taps and floating-plate skin friction force balances. The TBL developed a mild favourable pressure gradient that led to a streamwise flow speed increase of ~2.5% over the 11m long test surface, and was consistent with test section sidewall and model surface boundary-layer growth. At each Reθ, mean streamwise velocity profile pairs, separated by 24cm, were measured more than 10m from the model's leading edge using conventional laser Doppler velocimetry. Between these profile pairs, a unique near-wall implementation of particle tracking velocimetry was used to measure the near-wall velocity profile. The composite profile measurements span the wall-normal coordinate range from y+ < 1 to y > 2δ. To within experimental uncertainty, the measured mean velocity profiles can be fit using traditional zero-pressure-gradient (ZPG) TBL asymptotics with some modifications for the mild favourable pressure gradient. The fitted profile pairs satisfy the von-Kármán momentum integral equation to within 1%. However, the profiles reported here show distinct differences from equivalent ZPG profiles. The near-wall indicator function has more prominent extrema, the log-law constants differ slightly, and the profiles' wake component is less pronounced.

The effects of rotation on turbulent flow with separation and reattachment are investigated by means of direct numerical simulations. The backward-facing step configuration is rotated about a spanwise axis such that the sudden expansion of the channel is on the pressure side. The upstream flow is a fully developed plane Poiseuille flow subjected to orthogonal-mode rotation, which subsequently detaches from the step corner and eventually reattaches further downstream. The size of the resulting separation bubble with recirculating flow diminishes monotonically with increasing rotation rates and the reattachment distance is reduced from about 7 to 3 step heights. This is ascribed to the augmentation of the cross-stream turbulence intensity in the anti-cyclonic shear layer formed between the bulk flow and the recirculating eddy due to the destabilizing influence of the Coriolis force. The spanwise-oriented vortex cells or roller eddies found in non-rotating shear layers were disrupted by the enhanced turbulence. The flow along the planar wall is subjected to an adverse pressure gradient induced by the sudden expansion. The stabilizing influence of the system rotation in this cyclonic shear layer tends to damp the turbulence, the flow becomes susceptible to flow separation, and a substantial cyclonic recirculation bubble is observed at the highest rotation rates. The resulting meandering of the bulk flow is associated with interactions between the anti-cyclonic shear layer at the stepped side and the cyclonic shear flow along the planar surface. These give rise to enhanced turbulence levels at the cyclonic side in spite of the otherwise stabilizing influence of the Coriolis force. Exceptionally high velocity fluctuations in the spanwise direction are observed in the vicinity of flow reattachment behind the step and ascribed to longitudinal Taylor–Görtler-like roll cells which extend into the backflow region. These roll cells arise from a centrifugal instability mechanism associated with the convex streamline curvature in the reattachment zone.

The two-dimensional problem of gravity-driven laminar flow of a thin layer of fluid down a heated wavy inclined surface is discussed. The coupled effect of bottom topography, variable surface tension and heating has been investigated both analytically and numerically. A stability analysis is conducted while nonlinear simulations are used to validate the stability predictions and also to study thermocapillary effects. The governing equations are based on the Navier–Stokes equations for a thin fluid layer with the cross-stream dependence eliminated by means of a weighted residual technique. Comparisons with experimental data and direct numerical simulations have been carried out and the agreement is good. New interesting results regarding the combined role of surface tension and sinusoidal topography on the stability of the flow are presented. The influence of heating and the Marangoni effect are also deduced.

Two-dimensional forced liquid sloshing in a circular tank is studied by the multimodal method which uses an expansion in terms of the natural modes of free oscillations in the unforced tank. Incompressible inviscid liquid, irrotational flow and linear free-surface conditions are assumed. Accurate natural sloshing modes are constructed in an analytical form. Based on these modes, the ‘multimodal’ velocity potential of both steady-state and transient forced liquid motions exactly satisfies the body-boundary condition, captures the corner-point behaviour between the mean free surface and the tank wall and accurately approximates the free-surface conditions. The constructed multimodal solution provides an accurate description of the linear forced liquid sloshing. Surface wave elevations and hydrodynamic loads are compared with known experimental and nonlinear computational fluid dynamics results. The linear multimodal sloshing solution demonstrates good agreement in transient conditions of small duration, but fails in steady-state nearly-resonant conditions. Importance of the free-surface nonlinearity with increasing tank filling is explained.

A field experiment – the Snow Horizontal Array Turbulence Study (SnoHATS) – has been performed over an extensive glacier in Switzerland in order to study small-scale turbulence in the stable atmospheric surface layer, and to investigate the role, dynamics and modelling of the subgrid scales (SGSs) in the context of large-eddy simulations. The a priori data analysis aims at comparing the role and behaviour of the SGSs under stable conditions with previous studies under neutral or unstable conditions. It is found that the SGSs in a stable surface layer remain an important sink of temperature variance and turbulent kinetic energy from the resolved scales and carry a significant portion of the fluxes when the filter scale is larger than the distance to the wall. The fraction of SGS fluxes (out of the total fluxes) is found to be independent of stability. In addition, the stress–strain alignment is similar to the alignment under neutral and unstable conditions. The model coefficients vary considerably with stability but in a manner consistent with previous findings, which also showed that scale-dependent dynamic models can capture this variation. Furthermore, the variation of the coefficients for both momentum and heat SGS fluxes can be shown to be better explained by stability parameters based on vertical gradients, rather than vertical fluxes. These findings suggest that small-scale turbulence dynamics and SGS modelling under stable conditions share many important properties with neutral and convective conditions, and that a unified approach is thus possible. This paper concludes with a discussion of some other challenges for stable boundary-layer simulations that are not encountered in the neutral or unstable cases.

The analysis of flow over a slowly perturbed bottom (when perturbations have a typical length scale much larger than channel height) is often based on the shallow-water (or Saint-Venant) equations with the addition of a wall-friction term which is a local function of the mean velocity. By this choice, small sinusoidal disturbances of wall stress and mean velocity are bound to be in phase with each other. In contrast, studies of shorter-scale disturbances have long established that a phase lead develops between wall stress and mean velocity, with a crucial destabilizing effect on sediment transport along an erodible bed. The purpose of this paper is to calculate the wall-shear stress under large length-scale conditions and provide corrections to the Saint-Venant model.