Showing papers in "Journal of Fluid Mechanics in 2010"

Dynamic mode decomposition of numerical and experimental data

Abstract: The description of coherent features of fluid flow is essential to our understanding of fluid-dynamical and transport processes. A method is introduced that is able to extract dynamic information from flow fields that are either generated by a (direct) numerical simulation or visualized/measured in a physical experiment. The extracted dynamic modes, which can be interpreted as a generalization of global stability modes, can be used to describe the underlying physical mechanisms captured in the data sequence or to project large-scale problems onto a dynamical system of significantly fewer degrees of freedom. The concentration on subdomains of the flow field where relevant dynamics is expected allows the dissection of a complex flow into regions of localized instability phenomena and further illustrates the flexibility of the method, as does the description of the dynamics within a spatial framework. Demonstrations of the method are presented consisting of a plane channel flow, flow over a two-dimensional cavity, wake flow behind a flexible membrane and a jet passing between two cylinders.

Assessment of direct numerical simulation data of turbulent boundary layers

Abstract: Statistics obtained from seven different direct numerical simulations (DNSs) pertaining to a canonical turbulent boundary layer (TBL) under zero pressure gradient are compiled and compared. The considered data sets include a recent DNS of a TBL with the extended range of Reynolds numbers Reθ = 500–4300. Although all the simulations relate to the same physical flow case, the approaches differ in the applied numerical method, grid resolution and distribution, inflow generation method, boundary conditions and box dimensions. The resulting comparison shows surprisingly large differences not only in both basic integral quantities such as the friction coefficient cf or the shape factor H12, but also in their predictions of mean and fluctuation profiles far into the sublayer. It is thus shown that the numerical simulation of TBLs is, mainly due to the spatial development of the flow, very sensitive to, e.g. proper inflow condition, sufficient settling length and appropriate box dimensions. Thus, a DNS has to be considered as a numerical experiment and should be the subject of the same scrutiny as experimental data. However, if a DNS is set up with the necessary care, it can provide a faithful tool to predict even such notoriously difficult flow cases with great accuracy.

A critical-layer framework for turbulent pipe flow

Abstract: A model-based description of the scaling and radial location of turbulent fluctuations in turbulent pipe flow is presented and used to illuminate the scaling behaviour of the very large scale motions. The model is derived by treating the nonlinearity in the perturbation equation (involving the Reynolds stress) as an unknown forcing, yielding a linear relationship between the velocity field response and this nonlinearity. We do not assume small perturbations. We examine propagating helical velocity response modes that are harmonic in the wall-parallel directions and in time, permitting comparison of our results to experimental data. The steady component of the velocity field that varies only in the wall-normal direction is identified as the turbulent mean profile. A singular value decomposition of the resolvent identifies the forcing shape that will lead to the largest velocity response at a given wavenumber–frequency combination. The hypothesis that these forcing shapes lead to response modes that will be dominant in turbulent pipe flow is tested by using physical arguments to constrain the range of wavenumbers and frequencies to those actually observed in experiments. An investigation of the most amplified velocity response at a given wavenumber–frequency combination reveals critical-layer-like behaviour reminiscent of the neutrally stable solutions of the Orr–Sommerfeld equation in linearly unstable flow. Two distinct regions in the flow where the influence of viscosity becomes important can be identified, namely wall layers that scale with R+1/2 and critical layers where the propagation velocity is equal to the local mean velocity, one of which scales with R+2/3 in pipe flow. This framework appears to be consistent with several scaling results in wall turbulence and reveals a mechanism by which the effects of viscosity can extend well beyond the immediate vicinity of the wall. The model reproduces inner scaling of the small scales near the wall and an approach to outer scaling in the flow interior. We use our analysis to make a first prediction that the appropriate scaling velocity for the very large scale motions is the centreline velocity, and show that this is in agreement with experimental results. Lastly, we interpret the wall modes as the motion required to meet the wall boundary condition, identifying the interaction between the critical and wall modes as a potential origin for an interaction between the large and small scales that has been observed in recent literature as an amplitude modulation of the near-wall turbulence by the very large scales.

Discrete particle simulation of particle–fluid flow: model formulations and their applicability

Abstract: The approach of combining computational fluid dynamics (CFD) for continuum fluid and the discrete element method (DEM) for discrete particles has been increasingly used to study the fundamentals of coupled particle–fluid flows. Different CFD–DEM models have been used. However, the origin and the applicability of these models are not clearly understood. In this paper, the origin of different model formulations is discussed first. It shows that, in connection with the continuum approach, three sets of formulations exist in the CFD–DEM approach: an original format set I, and subsequent derivations of set II and set III, respectively, corresponding to the so-called model A and model B in the literature. A comparison and the applicability of the three models are assessed theoretically and then verified from the study of three representative particle–fluid flow systems: fluidization, pneumatic conveying and hydrocyclones. It is demonstrated that sets I and II are essentially the same, with small differences resulting from different mathematical or numerical treatments of a few terms in the original equation. Set III is however a simplified version of set I. The testing cases show that all the three models are applicable to gas fluidization and, to a large extent, pneumatic conveying. However, the application of set III is conditional, as demonstrated in the case of hydrocyclones. Strictly speaking, set III is only valid when fluid flow is steady and uniform. Set II and, in particular, set I, which is somehow forgotten in the literature, are recommended for the future CFD–DEM modelling of complex particle–fluid flow.

Sharp-interface limit of the Cahn-Hilliard model for moving contact lines

Abstract: Diffuse-interface models may be used to compute moving contact lines because the Cahn–Hilliard diffusion regularizes the singularity at the contact line. This paper investigates the basic questions underlying this approach. Through scaling arguments and numerical computations, we demonstrate that the Cahn–Hilliard model approaches a sharp-interface limit when the interfacial thickness is reduced below a threshold while other parameters are fixed. In this limit, the contact line has a diffusion length that is related to the slip length in sharp-interface models. Based on the numerical results, we propose a criterion for attaining the sharp-interface limit in computing moving contact lines.

Turbulent boundary layers and channels at moderate Reynolds numbers

Abstract: The behaviour of the velocity and pressure fluctuations in the outer layers of wall-bounded turbulent flows is analysed by comparing a new simulation of the zero-pressure-gradient boundary layer with older simulations of channels. The 99 % boundary-layer thickness is used as a reasonable analogue of the channel half-width, but the two flows are found to be too different for the analogy to be complete. In agreement with previous results, it is found that the fluctuations of the transverse velocities and of the pressure are stronger in the boundary layer, and this is traced to the pressure fluctuations induced in the outer intermittent layer by the differences between the potential and rotational flow regions. The same effect is also shown to be responsible for the stronger wake component of the mean velocity profile in external flows, whose increased energy production is the ultimate reason for the stronger fluctuations. Contrary to some previous results by our group, and by others, the streamwise velocity fluctuations are also found to be higher in boundary layers, although the effect is weaker. Within the limitations of the non-parallel nature of the boundary layer, the wall-parallel scales of all the fluctuations are similar in both the flows, suggesting that the scale-selection mechanism resides just below the intermittent region, y/δ = 0.3–0.5. This is also the location of the largest differences in the intensities, although the limited Reynolds number of the boundary-layer simulation (Reθ ≈ 2000) prevents firm conclusions on the scaling of this location. The statistics of the new boundary layer are available from http://torroja.dmt.upm.es/ftp/blayers/.

Three-dimensional direct numerical simulation of a turbulent lifted hydrogen jet flame in heated coflow: a chemical explosive mode analysis

Abstract: A chemical explosive mode analysis (CEMA) was developed as a new diagnostic to identify flame and ignition structure in complex flows. CEMA was then used to analyse the near-field structure of the stabilization region of a turbulent lifted hydrogen–air slot jet flame in a heated air coflow computed with three-dimensional direct numerical simulation. The simulation was performed with a detailed hydrogen–air mechanism and mixture-averaged transport properties at a jet Reynolds number of 11000 with over 900 million grid points. Explosive chemical modes and their characteristic time scales, as well as the species involved, were identified from the Jacobian matrix of the chemical source terms for species and temperature. An explosion index was defined for explosive modes, indicating the contribution of species and temperature in the explosion process. Radical and thermal runaway can consequently be distinguished. CEMA of the lifted flame shows the existence of two premixed flame fronts, which are difficult to detect with conventional methods. The upstream fork preceding the two flame fronts thereby identifies the stabilization point. A Damkohler number was defined based on the time scale of the chemical explosive mode and the local instantaneous scalar dissipation rate to highlight the role of auto-ignition in affecting the stabilization points in the lifted jet flame.

Radial boundary layer structure and Nusselt number in Rayleigh–Bénard convection

Abstract: Results from direct numerical simulation (DNS) for three-dimensional Rayleigh–Benard convection in a cylindrical cell of aspect ratio 1/2 and Prandtl number Pr=0.7 are presented. They span five decades of Rayleigh number Ra from 2 × 106 to 2 × 1011. The results are in good agreement with the experimental data of Niemela et al. (Nature, vol. 404, 2000, p. 837). Previous DNS results from Amati et al. (Phys. Fluids, vol. 17, 2005, paper no. 121701) showed a heat transfer that was up to 30% higher than the experimental values. The simulations presented in this paper are performed with a much higher resolution to properly resolve the plume dynamics. We find that in under-resolved simulations the hot (cold) plumes travel further from the bottom (top) plate than in the better-resolved ones, because of insufficient thermal dissipation mainly close to the sidewall (where the grid cells are largest), and therefore the Nusselt number in under-resolved simulations is overestimated. Furthermore, we compare the best resolved thermal boundary layer profile with the Prandtl–Blasius profile. We find that the boundary layer profile is closer to the Prandtl–Blasius profile at the cylinder axis than close to the sidewall, because of rising plumes close to the sidewall

On the wake-induced vibration of tandem circular cylinders: the vortex interaction excitation mechanism

Abstract: The mechanism of wake-induced vibrations (WIV) of a pair of cylinders in a tandem arrangement is investigated by experiments. A typical WIV response is characterized by a build-up of amplitude persisting to high reduced velocities; this is different from a typical vortex-induced vibration (VIV) response, which occurs in a limited resonance range. We suggest that WIV of the downstream cylinder is excited by the unsteady vortex–structure interactions between the body and the upstream wake. Coherent vortices interfering with the downstream cylinder induce fluctuations in the fluid force that are not synchronized with the motion. A favourable phase lag between the displacement and the fluid force guarantees that a positive energy transfer from the flow to the structure sustains the oscillations. If the unsteady vortices are removed from the wake of the upstream body then WIV will not be excited. An experiment performed in a steady shear flow turned out to be central to the understanding of the origin of the fluid forces acting on the downstream cylinder.

Modulation of isotropic turbulence by particles of Taylor length-scale size

Abstract: This study investigates the two-way coupling effects of finite-size solid spherical particles on decaying isotropic turbulence using direct numerical simulation with an immersed boundary method. We fully resolve all the relevant scales of turbulence around freely moving particles of the Taylor length-scale size, 1.2≤d/λ≤2.6. The particle diameter and Stokes number in terms of Kolmogorov length- and time scales are 16≤d/η≤35 and 38≤τp/τk≤178, respectively, at the time the particles are released in the flow. The particles mass fraction range is 0.026≤φm≤1.0, corresponding to a volume fraction of 0.01≤φv≤0.1 and density ratio of 2.56≤ρp/ρf≤10. The maximum number of dispersed particles is 6400 for φv=0.1. The typical particle Reynolds number is of O(10). The effects of the particles on the temporal development of turbulence kinetic energy E(t), its dissipation rate (t), its two-way coupling rate of change Ψp(t) and frequency spectra E(ω) are discussed.In contrast to particles with d η, is that E(t) is always smaller than that of the single-phase flow. In addition, Ψp(t) is always positive for particles with d > η, whereas it can be positive or negative for particles with d < η.

Non-equilibrium effects in capillarity and interfacial area in two-phase flow: dynamic pore-network modelling

Abstract: Current macroscopic theories of two-phase flow in porous media are based on the extended Darcy's law and an algebraic relationship between capillary pressure and saturation. Both of these equations have been challenged in recent years, primarily based on theoretical works using a thermodynamic approach, which have led to new governing equations for two-phase flow in porous media. In these equations, new terms appear related to the fluid–fluid interfacial area and non-equilibrium capillarity effects. Although there has been a growing number of experimental works aimed at investigating the new equations, a full study of their significance has been difficult as some quantities are hard to measure and experiments are costly and time-consuming. In this regard, pore-scale computational tools can play a valuable role. In this paper, we develop a new dynamic pore-network simulator for two-phase flow in porous media, called DYPOSIT. Using this tool, we investigate macroscopic relationships among average capillary pressure, average phase pressures, saturation and specific interfacial area. We provide evidence that at macroscale, average capillary pressure–saturation–interfacial area points fall on a single surface regardless of flow conditions and fluid properties. We demonstrate that the traditional capillary pressure–saturation relationship is not valid under dynamic conditions, as predicted by the theory. Instead, one has to employ the non-equilibrium capillary theory, according to which the fluids pressure difference is a function of the time rate of saturation change. We study the behaviour of non-equilibrium capillarity coefficient, specific interfacial area, and its production rate versus saturation and viscosity ratio.A major feature of our pore-network model is a new computational algorithm, which considers capillary diffusion. Pressure field is calculated for each fluid separately, and saturation is computed in a semi-implicit way. This provides more numerical stability, compared with previous models, especially for unfavourable viscosity ratios and small capillary number values.

Influence of nozzle-exit boundary-layer conditions on the flow and acoustic fields of initially laminar jets

Abstract: Round jets originating from a pipe nozzle are computed by large-eddy simulations (LES) to investigate the effects of the nozzle-exit conditions on the flow and sound fields of initially laminar jets. The jets are at Mach number 0.9 and Reynolds number 105, and exhibit exit boundary layers characterized by Blasius velocity profiles, maximum root-mean-square (r.m.s.) axial velocity fluctuations between 0.2 and 1.9% of the jet velocity, and momentum thicknesses varying from 0.003 to 0.023 times the jet radius. The far-field noise is determined from the LES data on a cylindrical surface by solving the acoustic equations. Jets with a thinner boundary layer develop earlier but at a slower rate, yielding longer potential cores and lower centreline turbulent intensities. Adding random pressure disturbances of low magnitude in the nozzle also increases the potential core length and reduces peak r.m.s. radial velocity fluctuations in the shear layer. In all the jets, the shear-layer transition is dominated by vortex rolling-ups and pairings, which generate strong additional acoustic components, but also amplify the downstream-dominant low-frequency noise component when the exit boundary layer is thick. The introduction of inlet noise however results in weaker pairings, thus spectacularly reducing their contributions to the sound field. This high sensitivity to the initial conditions is in good agreement with experimental observations.

Direct numerical simulation of hypersonic turbulent boundary layers. Part 2. Effect of wall temperature

Abstract: In this paper, we perform direct numerical simulation (DNS) of turbulent boundary layers at Mach 5 with the ratio of wall-to-edge temperature Tw/Tδ from 1.0 to 5.4 (Cases M5T1 to M5T5). The influence of wall cooling on Morkovin's scaling, Walz's equation, the standard and modified strong Reynolds analogies, turbulent kinetic energy budgets, compressibility effects and near-wall coherent structures is assessed. We find that many of the scaling relations used to express adiabatic compressible boundary-layer statistics in terms of incompressible boundary layers also hold for non-adiabatic cases. Compressibility effects are enhanced by wall cooling but remain insignificant, and the turbulence dissipation remains primarily solenoidal. Moreover, the variation of near-wall streaks, iso-surface of the swirl strength and hairpin packets with wall temperature demonstrates that cooling the wall increases the coherency of turbulent structures. We present the mechanism by which wall cooling enhances the coherence of turbulence structures, and we provide an explanation of why this mechanism does not represent an exception to the weakly compressible hypothesis.

Aeromechanics of passive rotation in flapping flight

Abstract: Flying insects and robots that mimic them flap and rotate (or ‘pitch’) their wings with large angular amplitudes. The reciprocating nature of flapping requires rotation of the wing at the end of each stroke. Insects or flapping-wing robots could achieve this by directly exerting moments about the axis of rotation using auxiliary muscles or actuators. However, completely passive rotational dynamics might be preferred for efficiency purposes, or, in the case of a robot, decreased mechanical complexity and reduced system mass. Herein, the detailed equations of motion are derived for wing rotational dynamics, and a blade-element model is used to supply aerodynamic force and moment estimates. Passive-rotation flapping experiments with insect-scale mechanically driven artificial wings are conducted to simultaneously measure aerodynamic forces and three-degree-of-freedom kinematics (flapping, rotation and out-of-plane deviation), allowing a detailed evaluation of the blade-element model and the derived equations of motion. Variations in flapping kinematics, wing-beat frequency, stroke amplitude and torsional compliance are made to test the generality of the model. All experiments showed strong agreement with predicted forces and kinematics, without variation or fitting of model parameters.

Formation of turbulent patterns near the onset of transition in plane Couette flow

Abstract: The formation of turbulent patterns in plane Couette flow is investigated near the onset of transition, using numerical simulation in a very large domain of size 800 h x 2h x 356 h. Based on a maxi ...

Channel flow over large cube roughness: a direct numerical simulation study

Abstract: Computations of channel flow with rough walls comprising staggered arrays of cubes having various plan area densities are presented and discussed. The cube height h is 12.5% of the channel half-depth and Reynolds numbers (u? h/?) are typically around 700 – well into the fully rough regime. A direct numerical simulation technique, using an immersed boundary method for the obstacles, was employed with typically 35 million cells. It is shown that the surface drag is predominantly form drag, which is greatest at an area coverage around 15%. The height variation of the axial pressure force across the obstacles weakens significantly as the area coverage decreases, but is always largest near the top of the obstacles. Mean flow velocity and pressure data allow precise determination of the zero-plane displacement (defined as the height at which the axial surface drag force acts) and this leads to noticeably better fits to the log-law region than can be obtained by using the zero-plane displacement merely as a fitting parameter. There are consequent implications for the value of von K´arm´ an’s constant. As the effective roughness of the surface increases, it is also shown that there are significant changes to the structure of the turbulence field around the bottom boundary of the inertial sublayer. In distinct contrast to twodimensional roughness (longitudinal or transverse bars), increasing the area density of this three-dimensional roughness leads to a monotonic decrease in normalized vertical stress around the top of the roughness elements. Normalized turbulence stresses in the outer part of the flows are nonetheless very similar to those in smooth-wall flows.

Linear non-normal energy amplification of harmonic and stochastic forcing in the turbulent channel flow

Abstract: The linear response to stochastic and optimal harmonic forcing of small coherent perturbations to the turbulent channel mean flow is computed for Reynolds numbers ranging from Re_tau=500 to Re_tau=20000. Even though the turbulent mean flow is linearly stable, it is nevertheless able to sustain large amplifications by the forcing. The most amplified structures consist of streamwise elongated streaks that are optimally forced by streamwise elongated vortices. For streamwise elongated structures, the mean energy amplification of the stochastic forcing is found to be, to a first approximation, inversely proportional to the forced spanwise wavenumber while it is inversely proportional to its square for optimal harmonic forcing in an intermediate spanwise wavenumber range. This scaling can be explicitly derived from the linearised equations under the assumptions of geometric similarity of the coherent perturbations and of logarithmic base flow. Deviations from this approximate power-law regime are apparent in the premultiplied energy amplification curves that reveal a strong influence of two different peaks. The dominant peak scales in outer units with the most amplified spanwise wavelength of $\lambda_z \approx 3.5 h$ while the secondary peak scales in wall units with the most amplified $\lambda_z^+\approx 80$. The associated optimal perturbations are almost independent of the Reynolds number when respectively scaled in outer and inner units. In the intermediate wavenumber range the optimal perturbations are approximatively geometrically similar. Furthermore, the shape of the optimal perturbations issued from the initial value, the harmonic forcing and the stochastic forcing analyses are almost indistinguishable. The optimal streaks corresponding to the large-scale peak strongly penetrate into the inner layer, where their amplitude is proportional to the mean-flow profile. At the wavenumbers corresponding to the large-scale peak, the optimal amplifications of harmonic forcing are at least two orders of magnitude larger than the amplifications of the variance of stochastic forcing and both increase with the Reynolds number. This confirms the potential of the artificial forcing of optimal large-scale streaks for the flow control of wall-bounded turbulent flows.

Balanced and Unbalanced Routes to Dissipation in an Equilibrated, Eady Flow

Abstract: The oceanic general circulation is forced at large scales and is unstable to mesoscale eddies. Large-scale currents and eddy flows are approximately in geostrophic balance. Geostrophic dynamics is characterized by an inverse energy cascade except for dissipation near the boundaries. In this paper, we confront the dilemma of how the general circulation may achieve dynamical equilibrium in the presence of continuous large-scale forcing and the absence of boundary dissipation. We do this with a forced horizontal flow with spatially uniform rotation, vertical stratification and vertical shear in a horizontally periodic domain, i.e. a version of Eady's flow carried to turbulent equilibrium. A direct route to interior dissipation is presented that is essentially non-geostrophic in its dynamics, with significant submesoscale frontogenesis, frontal instability and breakdown, and forward kinetic energy cascade to dissipation. To support this conclusion, a series of simulations is made with both quasigeostrophic and Boussinesq models. The quasigeostrophic model is shown as increasingly inefficient in achieving equilibration through viscous dissipation at increasingly higher numerical resolution (hence Reynolds number), whereas the non-geostrophic Boussinesq model equilibrates with only weak dependence on resolution and Rossby number.

Spatial structure of a turbulent boundary layer with irregular surface roughness

Abstract: Particle image velocimetry experiments were performed to study the impact of realistic roughness on the spatial structure of wall turbulence at moderate Reynolds number. This roughness was replicated from an actual turbine blade damaged by deposition of foreign materials and its features are quite distinct from most roughness characterizations previously considered as it is highly irregular and embodies a broad range of topographical scales. The spatial structure of flow over this rough surface near the outer edge of the roughness sublayer is contrasted with that of smooth-wall flow to identify any structural modifications due to roughness. Hairpin vortex packets are observed in the outer layer of the rough-wall flow and are found to contribute heavily to the Reynolds shear stress, consistent with smooth-wall flow. While similar qualitative consistency is observed in comparisons of smooth- and rough-wall two-point correlations, some quantitative differences are also apparent. In particular, a reduction in the streamwise extent of two-point correlations of streamwise velocity is noted which could be indicative of a roughness-induced modification of outer-layer vortex organization. Proper orthogonal decomposition analysis reveals the streamwise coherence of the larger scales to be most sensitive to roughness while the spatial characteristics of the smaller scales appear relatively insensitive to such effects.

Turbulent flame–wall interaction: a direct numerical simulation study

Abstract: A turbulent flame–wall interaction (FWI) configuration is studied using three-dimensional direct numerical simulation (DNS) and detailed chemical kinetics. The simulations are used to investigate the effects of the wall turbulent boundary layer (i) on the structure of a hydrogen–air premixed flame, (ii) on its near-wall propagation characteristics and (iii) on the spatial and temporal patterns of the convective wall heat flux. Results show that the local flame thickness and propagation speed vary between the core flow and the boundary layer, resulting in a regime change from flamelet near the channel centreline to a thickened flame at the wall. This finding has strong implications for the modelling of turbulent combustion using Reynolds-averaged Navier–Stokes or large-eddy simulation techniques. Moreover, the DNS results suggest that the near-wall coherent turbulent structures play an important role on the convective wall heat transfer by pushing the hot reactive zone towards the cold solid surface. At the wall, exothermic radical recombination reactions become important, and are responsible for approximately 70% of the overall heat release rate at the wall. Spectral analysis of the convective wall heat flux provides an unambiguous picture of its spatial and temporal patterns, previously unobserved, that is directly related to the spatial and temporal characteristic scalings of the coherent near-wall turbulent structures.

Streamwise vortices in shear flows: harbingers of transition and the skeleton of coherent structures

Abstract: The relationship between asymptotic descriptions of vortex–wave interactions and more recent work on ‘exact coherent structures’ is investigated. In recent years immense interest has been focused on so-called self-sustained processes in turbulent shear flows where the importance of waves interacting with streamwise vortex flows has been elucidated in a number of papers. In this paper, it is shown that the so-called ‘lower branch’ state which has been shown to play a crucial role in these self-sustained processes is a finite Reynolds number analogue of a Rayleigh vortex–wave interaction with scales appropriately modified from those for external flows to Couette flow, the flow of interest here. Remarkable agreement between the asymptotic theory and numerical solutions of the Navier–Stokes equations is found even down to relatively small Reynolds numbers, thereby suggesting the possible importance of vortex–wave interaction theory in turbulent shear flows. The relevance of the work to more general shear flows is also discussed.

Drag reduction of flexible plates by reconfiguration

Abstract: Through an extensive and systematic experimental investigation of two geometries of flexible plates in air, it is shown that a properly defined scaled Cauchy number allows collapsing all drag measurements of the reconfiguration number. In the asymptotic regime of large deformation, it is shown that the Vogel exponents that scale the drag with the flow velocity for different geometries of plates can be predicted with a simple dimensional analysis reasoning. These predicted Vogel exponents are in agreement with previously published models of reconfiguration. The mechanisms responsible for reconfiguration, namely area reduction and streamlining, are studied with the help of a simple model for flexible plates based on an empirical drag formulation. The model predicts well the reconfiguration observed in the experiments and shows that for a rectangular plate, the effect of streamlining is prominent at the onset of reconfiguration, but area reduction dominates in the regime of large deformation. Additionally, the model demonstrates for both geometries of plates that the reconfiguration cannot be described by a single value of the Vogel exponent. The Vogel exponent asymptotically approaches constant values for small and for very large scaled Cauchy numbers, but in between both extremes it varies significantly over a large range of scaled Cauchy number.

The influence of initial conditions on turbulent mixing due to Richtmyer–Meshkov instability†

Abstract: This paper investigates the influence of different three-dimensional multi-mode initial conditions on the rate of growth of a mixing layer initiated via a Richtmyer–Meshkov instability through a series of well-controlled numerical experiments. Results are presented for large-eddy simulation of narrowband and broadband perturbations at grid resolutions up to 3 × 109 points using two completely different numerical methods, and comparisons are made with theory and experiment. It is shown that the mixing-layer growth is strongly dependent on initial conditions, the narrowband case giving a power-law exponent θ ≈ 0.26 at low Atwood and θ ≈ 0.3 at high Atwood numbers. The broadband case uses a perturbation power spectrum of the form P(k) ∝ k−2 with a proposed theoretical growth rate of θ = 2/3. The numerical results confirm this; however, they highlight the necessity of a very fine grid to capture an appropriately broad range of initial scales. In addition, an analysis of the kinetic energy decay rates, fluctuating kinetic energy spectra, plane-averaged volume fraction profiles and mixing parameters is presented for each case.

Experimental characterization of the agitation generated by bubbles rising at high Reynolds number

Abstract: An experimental investigation of the flow generated by a homogeneous population of bubbles rising in water is reported for three different bubble diameters (d = 1.6, 2.1 and 2.5 mm) and moderate gas volume fractions (0.005≤α≤0.1). The Reynolds numbers, Re = V 0 d/ν, based on the rise velocity V 0 of a single bubble range between 500 and 800. Velocity statistics of both the bubbles and the liquid phase are determined within the homogeneous bubble swarm by means of optical probes and laser Doppler anemometry. Also, the decaying agitation that takes place in the liquid just after the passage of the bubble swarm is investigated from high-speed particle image velocimetry measurements. Concerning the bubbles, the average velocity is found to evolve as V 0 α -0.1 whereas the velocity fluctuations are observed to be almost independent of α. Concerning the liquid fluctuations, the probability density functions adopt a self-similar behaviour when the gas volume fraction is varied, the characteristic velocity scaling as V 0 α 0.4 . The spectra of horizontal and vertical liquid velocity fluctuations are obtained with a resolution of 0.6 mm. The integral length scale Λ is found to be proportional to V 2 0 /g or equivalently to d/C d0 , where g is the gravity acceleration and C d0 the drag coefficient of a single rising bubble. Normalized by using Λ, the spectra are independent on both the bubble diameter and the volume fraction. At large scales, the spectral energy density evolves as the power -3 of the wavenumber. This range starts approximately from Λ and is followed for scales smaller than Λ/4 by a classic -5/3 power law. Although the Kolmogorov microscale is smaller than the measurement resolution, the dissipation rate is however obtained from the decay of the kinetic energy after the passage of the bubbles. It is found to scale as α 0.9 V 3 0 /Λ. The major characteristics of the agitation are thus expressed as functions of the characteristics of a single rising bubble. Altogether, these results provide a rather complete description of the bubble-induced turbulence.

Stability and receptivity characteristics of a laminar separation bubble on an aerofoil

Abstract: Stability characteristics of aerofoil flows are investigated by linear stability analysis of time-averaged velocity profiles and by direct numerical simulations with time-dependent forcing terms. First the wake behind an aerofoil is investigated, illustrating the feasibility of detecting absolute instability using these methods. The time-averaged flow around an NACA-0012 aerofoil at incidence is then investigated in terms of its response to very low-amplitude hydrodynamic and acoustic perturbations. Flow fields obtained from both two- and three-dimensional simulations are investigated, for which the aerofoil flow exhibits a laminar separation bubble. Convective stability characteristics are documented, and the separation bubble is found to exhibit no absolute instability in the classical sense; i.e. no growing disturbances with zero group velocity are observed. The flow is however found to be globally unstable via an acoustic-feedback loop involving the aerofoil trailing edge as a source of acoustic excitation and the aerofoil leading-edge region as a site of receptivity. Evidence suggests that the feedback loop may play an important role in frequency selection of the vortex shedding that occurs in two dimensions. Further simulations are presented to investigate the receptivity process by which acoustic waves generate hydrodynamic instabilities within the aerofoil boundary layer. The dependency of the receptivity process to both frequency and source location is quantified. It is found that the amplitude of trailing-edge noise in the fully developed simulation is sufficient to promote transition via leading-edge receptivity.

Tuning turbines in a tidal channel

Abstract: As tidal turbine farms grow they interact with the larger scale flow along a channel by increasing the channel's drag coefficient. This interaction limits a channel's potential to produce power. A 1D model for a tidal channel is combined with a theory for turbines in a channel to show that the tuning of the flow through the turbines and the density of turbines in a channel's cross-section also interact with the larger scale flow, via the drag coefficient, to determine the power available for production. To maximise turbine efficiency, i.e. the power available per turbine, farms must occupy the largest fraction of a channel's cross-section permitted by navigational and environmental constraints. Maximising of power available with these necessarily densely packed farms requires turbines to be tuned for a particular channel and turbine density. The optimal through-flow tuning fraction varies from near 1/3 for small farms occupying a small fraction of the cross-section, to near 1 for large farms occupying most of the cross-section. Consequently, tunings are higher than the optimal through-flow tuning of 1/3 for an isolated turbine from the classic turbine theory. Large optimally tuned farms can realise most of a channel's potential. Optimal tunings are dependent on the number of turbines per row, the number of rows, as well as the channel geometry, the background bottom friction coefficient and the tidal forcing.

Deformation of a hydrophobic ferrofluid droplet suspended in a viscous medium under uniform magnetic fields

Abstract: The effect of applied magnetic fields on the deformation of a biocompatible hydrophobic ferrofluid drop suspended in a viscous medium is investigated numerically and compared with experimental data. A numerical formulation for the time-dependent simulation of magnetohydrodynamics of two immiscible non-conducting fluids is used with a volume-of-fluid scheme for fully deformable interfaces. Analytical formulae for ellipsoidal drops and near-spheroidal drops are reviewed and developed for code validation. At low magnetic fields, both the experimental and numerical results follow the asymptotic small deformation theory. The value of interfacial tension is deduced from an optimal fit of a numerically simulated shape with the experimentally obtained drop shape, and appears to be a constant for low applied magnetic fields. At high magnetic fields, on the other hand, experimental measurements deviate from numerical results if a constant interfacial tension is implemented. The difference can be represented as a dependence of apparent interfacial tension on the magnetic field. This idea is investigated computationally by varying the interfacial tension as a function of the applied magnetic field and by comparing the drop shapes with experimental data until a perfect match is found. This estimation method provides a consistent correlation for the variation in interfacial tension at high magnetic fields. A conclusion section provides a discussion of physical effects which may influence the microstructure and contribute to the reported observations.

Dual resonance in vortex-induced vibrations at subcritical and supercritical Reynolds numbers

Abstract: An experimental study is performed on the vortex induced vibrations of a rigid flexibly mounted circular cylinder placed in a crossflow. The cylinder is allowed to oscillate in combined crossflow and in-line motions, and the ratio of the nominal in-line and transverse natural frequencies is varied systematically. Experiments were conducted on a smooth cylinder at subcritical Reynolds numbers between 15 000 and 60 000 and on a roughened cylinder at supercritical Reynolds numbers between 320 000 and 710 000, with a surface roughness equal to 0.23 % of the cylinder diameter. Strong qualitative and quantitative similarities between the subcritical and supercritical experiments are found, especially when the in-line natural frequency is close to twice the value of the crossflow natural frequency. In both Reynolds number regimes, the test cylinder may exhibit a 'dual-resonant' response, resulting in resonant crossflow motion at a frequency f υ , near the Strouhal frequency, and resonant in-line motion at 2 f υ . This dual resonance is shown to occur over a relatively wide frequency region around the Strouhal frequency, accompanied by stable, highly repeatable figure-eight cylinder orbits, as well as large third-harmonic components of the lift force. Under dual-resonance conditions, both the subcritical and the supercritical response is shown to collapse into a narrow parametric region in which the effective natural-frequency ratio is near the value 2, regardless of the nominal natural-frequency ratio. Some differences are noted in the magnitudes of forces and the cylinder response between the two different Reynolds number regimes, but the dual-resonant response and the resulting force trends are preserved despite the large Reynolds number difference.

Three-dimensional vortex organization in a high-Reynolds-number supersonic turbulent boundary layer

Abstract: Tomographic particle image velocimetry was used to quantitatively visualize the three-dimensional coherent structures in a supersonic (Mach 2) turbulent boundary layer in the region between y/? = 0.15 and 0.89. The Reynolds number based on momentum thickness Re? = 34000. The instantaneous velocity fields give evidence of hairpin vortices aligned in the streamwise direction forming very long zones of low-speed fluid, consistent with Tomkins & Adrian (J. Fluid Mech., vol. 490, 2003, p. 37). The observed hairpin structure is also a statistically relevant structure as is shown by the conditional average flow field associated to spanwise swirling motion. Spatial low-pass filtering of the velocity field reveals streamwise vortices and signatures of large-scale hairpins (height > 0.5?), which are weaker than the smaller scale hairpins in the unfiltered velocity field. The large-scale hairpin structures in the instantaneous velocity fields are observed to be aligned in the streamwise direction and spanwise organized along diagonal lines. Additionally the autocorrelation function of the wall-normal swirling motion representing the large-scale hairpin structure returns positive correlation peaks in the streamwise direction (at 1.5? distance from the DC peak) and along the 45° diagonals, which also suggest a periodic arrangement in those directions. This is evidence for the existence of a spanwise–streamwise organization of the coherent structures in a fully turbulent boundary layer.

The effect of Reynolds number on the dynamics and wakes of freely rising and falling spheres

Abstract: In this paper, we study the effect of the Reynolds number (Re) on the dynamics and vortex formation modes of spheres rising or falling freely through a fluid (where Re = 100–15000). Since the oscillation of freely falling spheres was first reported by Newton (University of California Press, 3rd edn, 1726, translated in 1999), the fundamental question of whether a sphere will vibrate, as it rises or falls, has been the subject of a number of investigations, and it is clear that the mass ratio m* (defined as the relative density of the sphere compared to the fluid) is an important parameter to define when vibration occurs. Although all rising spheres (m* 1) always move without vibration. However, in contrast with previous studies, we discover that a whole regime of buoyant spheres rise through a fluid without vibration. It is only when one passes below a critical value of the mass ratio, that the sphere suddenly begins to vibrate periodically and vigorously in a zigzag trajectory within a vertical plane. The critical mass is nearly constant over two ranges of Reynolds number (m*crit ≈ 0.4 for Re = 260–1550 and m*crit ≈ 0.6 for Re > 1550). We do not observe helical or spiral trajectories, or indeed chaotic types of trajectory, unless the experiments are conducted in disturbed background fluid. The wakes for spheres moving rectilinearly are comparable with wakes of non-vibrating spheres. We find that these wakes comprise single-sided and double-sided periodic sequences of vortex rings, which we define as the ‘R’ and ‘2R’ modes. However, in the zigzag regime, we discover a new ‘4R’ mode, in which four vortex rings are created per cycle of oscillation. We find a number of changes to occur at a Reynolds number of 1550, and we suggest the possibility of a resonance between the shear layer instability and the vortex shedding (loop) instability. From this study, ensuring minimal background disturbances, we have been able to present a new regime map of dynamics and vortex wake modes as a function of the mass ratio and Reynolds number {m*, Re}, as well as a reasonable collapse of the drag measurements, as a function of Re, onto principally two curves, one for the vibrating regime and one for the rectilinear trajectories.