Laboratoire des Écoulements Géophysiques et Industriels

International audience

The identification issue of coherent vortices is investigated on the basis of direct numerical simulation (DNS) and large-eddy simulations (LES) of turbulent flows. It is first shown that the pressure Laplacian is positive within a low-pressure tube of small cross section enclosed by convex isobaric surfaces in a uniform-density flow. Since this quantity is related to the second invariant Q of ∇ u , the Q criterion (region where Q is positive) is a necessary condition for the existence of such tubes. This eduction scheme is compared to other classical methods in incompressible simulations of isotropic turbulence: a mixing layer, a channel flow and a backward-facing step. Q-isosurfaces turn out to display very nice coherent vortices. This criterion is also used in combination with a conditional sampling method to discuss the characteristics of quasi-longitudinal vortices in a manipulated channel flow. The contribution of near-wall vortical structures to velocity and vorticity fluctuations is clearly isolated. Finally, the Q criterion is applied to the eduction of Λ vortices in separated flows. The dynamics of these structures is discussed for an incompressible backward-facing step and a transonic flow past a rectangular cavity.

The advances in our understanding of extratropical atmosphere-ocean interaction over the past decade and a half are examined, focusing on the atmospheric response to sea surface temperature anomalies.The main goal of the paper is to assess what was learned from general circulation model (GCM) experiments over the recent two decades or so.Observational evidence regarding the nature of the interaction and dynamical theory of atmospheric anomalies forced by surface thermal anomalies is reviewed.Three types of GCM experiments used to address this problem are then examined: models with fixed climatological conditions and idealized, stationary SST anomalies; models with seasonally evolving climatology forced with realistic, time-varying SST anomalies; and models coupled to an interactive ocean.From representative recent studies, it is argued that the extratropical atmosphere does respond to changes in underlying SST although the response is small compared to internal (unforced) variability.Two types of interactions govern the response.One is an eddy-mediated process, in which a baroclinic response to thermal forcing induces and combines with changes in the position or strength of the storm tracks.This process can lead to an equivalent barotropic response that feeds back positively on the ocean mixed layer temperature.The other is a linear, thermodynamic interaction in which an equivalent-barotropic low-frequency atmospheric anomaly forces a change in SST and then experiences reduced surface thermal damping due to the SST adjustment.Both processes contribute to an increase in variance and persistence of low-frequency atmospheric anomalies and, in fact, may act together in the natural system.

The objective of this work was to observe and quantify the onset and evolution of localised deformation processes in sand with grain-scale resolution. The key element of the proposed approach is combining state-of-the-art X-ray micro tomography imaging with three-dimensional volumetric digital image correlation techniques. This allows not only the grain-scale details of a deforming sand specimen to be viewed, but also, and more importantly, the evolving three-dimensional displacement and strain fields throughout loading to be assessed. X-ray imaging and digital image correlation have been in the past applied individually to study sand deformation, but the combination of these two methods to study the kinematics of shear band formation at the grain scale is the first novel aspect of this work. Moreover, the authors have developed a completely original grain-scale volumetric digital image correlation method that permits the characterisation of the full kinematics (i.e. three-dimensional displacements and rotations) of all the individual sand grains in a specimen. The results obtained using the discrete volumetric digital image correlation confirm the importance of grain rotations associated with strain localisation.

▪ Abstract An overview of the near and far-field breakup and atomization of a liquid jet by a high speed annular gas jet is presented. The various regimes of liquid jet breakup are discussed in the parameter space of the liquid Reynolds number, the aerodynamic Weber number, and the ratio of the momentum fluxes between the gas and the liquid streams. Recent measurements of the gas-liquid interfacial instabilities are reviewed and used to analyze the underlying physical mechanisms involved in the primary breakup of the liquid jet. This process is shown to consist of the periodic stripping of liquid sheets, or ligaments, which subsequently break up into smaller lumps or drops. Models to predict the liquid shedding frequency, as well as the global parameters of the spray such as the liquid core length and spray spreading angle are discussed and compared with the experiments. The role of the secondary liquid breakup on the far-field atomization of the liquid jet is also considered, and an attempt is made to apply the classical turbulent breakup concepts to explain qualitatively the measurement of the far-field droplet size distribution and its dependence on the liquid to gas mass and momentum flux ratios. Models for the droplet breakup frequency in the far-field region of the jet, and for the daughter-size probability density function, which account for the effect of the liquid loading on the local turbulent dissipation rate in the gas, are discussed in the context of the statistical description of the spray in the far field. The striking effect of the addition of swirl in the gas stream is also examined.

Hydration of the <1 μm size fraction of SWy-1 source clay (low-charge montmorillonite) was studied by modeling of X-ray diffraction (XRD) patterns recorded under controlled relative humidity (RH) conditions on Li-, Na-, K-, Mg-, Ca-, and Sr-saturated specimens. The quantitative description of smectite hydration, based on the relative proportions of different layer types derived from the fitting of experimental XRD patterns, was consistent with previous reports of smectite hydration. However, the coexistence of smectite layer types exhibiting contrasting hydration states was systematically observed, and heterogeneity rather than homogeneity seems to be the rule for smectite hydration. This heterogeneity can be characterized qualitatively using the standard deviation of the departure from rationality of the 00l reflection series (ξ), which is systematically larger than 0.4 Å when the prevailing layer type accounts for ~70% or less of the total layers (~25% of XRD patterns examined). In addition, hydration heterogeneities are not distributed randomly within smectite crystallites, and models describing these complex structures involve two distinct contributions, each containing different layer types that are interstratifed randomly. As a result, the different layer types are partially segregated in the sample. However, these two contributions do not imply the actual presence of two populations of particles in the sample. XRD profile modeling also has allowed the refinement of structural parameters, such as the location of interlayer species and the layer thickness corresponding to the different layer types, for all interlayer cations and RH values. From the observed dependence of the latter parameter on the cation ionic potential (v/r; v = cation valency and r = ionic radius) and on RH, the following equations were derived: Layer thickness (1W) = 12.556 + 0.3525 × (v/r . 0.241) × (v × RH . 0.979) Layer thickness (2W) = 15.592 + 0.6472 × (v/r . 0.839) × (v × RH . 1.412) which allow the quantification of the increase of layer thickness with increasing RH for both 1W (one water) and 2W (two water) layers. In addition, for 2W layers, interlayer H 2 O molecules are probably distributed as a unique plane on each side of the central interlayer cation. This plane of H 2 O molecules is located at ~1.20 Å from the central interlayer cation along the c* axis

Unsteady cavitation in a Venturi-type section was simulated by two-dimensional computations of viscous, compressible, and turbulent cavitating flows. The numerical model used an implicit finite volume scheme (based on the SIMPLE algorithm) to solve Reynolds-averaged Navier-Stokes equations, associated with a barotropic vapor/liquid state law that strongly links the density variations to the pressure evolution. To simulate turbulence effects on cavitating flows, four different models were implemented (standard k-ε RNG; modified k-ε RNG; k-ω with and without compressibility effects), and numerical results obtained were compared to experimental ones. The standard models k-ε RNG and k-ω without compressibility effects lead to a poor description of the self-oscillation behavior of the cavitating flow. To improve numerical simulations by taking into account the influence of the compressibility of the two-phase medium on turbulence, two other models were implemented in the numerical code: a modified k-ε model and the k-ω model including compressibility effects. Results obtained concerning void ratio, velocity fields, and cavitation unsteady behavior were found in good agreement with experimental ones. The role of the compressibility effects on turbulent two-phase flow modeling was analyzed, and it seemed to be of primary importance in numerical simulations.

The instability of a partial cavity induced by the development of a re-entrant jet is investigated on the basis of experiments conducted on a diverging step. Detailed visualizations of the cavity behaviour allowed us to identify the domain of the re-entrant jet instability which leads to classical cloud cavitation. The surrounding regimes are also investigated, in particular the special case of thin cavities which do not oscillate in length but surprisingly exhibit a re-entrant jet of periodical behaviour. The velocity of the re-entrant jet is measured from visualizations, in the case of both cloud cavitation and thin cavities. The limits of the domain of the re-entrant jet instability are corroborated by velocity fluctuation measurements. By varying the divergence and the confinement of the channel, it is shown that the extent of the auto-oscillation domain primarily depends upon the average adverse pressure gradient in the channel. This conclusion is corroborated by the determination of the pressure gradient on the basis of LDV measurements which shows a good correlation between the domain of the cloud cavitation instability and the region of high adverse pressure gradient. A simple phenomenological model of the development of the re-entrant jet in an adverse pressure gradient confirms the strong influence of the pressure gradient on the development of the re-entrant jet and particularly on its thickness. An ultrasonic technique is developed to measure the re-entrant jet thickness, which allowed us to compare it with the cavity thickness. By considering an estimate of the characteristic height of the perturbations developing on the interface of the cavity and of the re-entrant jet, it is shown that cloud cavitation requires negligible interaction between both interfaces, i.e. a thick enough cavity. In the case of thin cavities, this interaction becomes predominant; the cavity interface breaks at many points, giving birth to small-scale vapour structures unlike the large-scale clouds which are periodically shed in the case of cloud cavitation. The low-frequency content of the cloud cavitation instability is investigated using spectral analysis of wall pressure signals. It is shown that the characteristic frequency of cloud cavitation corresponds to a Strouhal number of about 0.2 whatever the operating conditions and the cavity length may be, provided the Strouhal number is computed on the basis of the maximum cavity length. For long enough cavities, another peak is observed in the spectra, at lower frequency, which is interpreted as a surge-type instability. The present investigations give insight into the instabilities that a partial cavity may undergo, and particularly the re-entrant jet instability. Two parameters are shown to be of most importance in the analysis of the re-entrant jet instability: the adverse pressure gradient and the cavity thickness compared to the re-entrant jet thickness. The present results allowed us to conduct a qualitative phenomenological analysis of the stability of partial cavities on cavitating hydrofoils. It is conjectured that cloud cavitation should occur for short enough cavities, of the order of half the chordlength, whereas the instability often observed at the limit between partial cavitation and super-cavitation is here interpreted as a cavitation surge-type instability.

Large-Eddy Simulations of Turbulence is an ideal introduction for people new to LES [large-eddy simulation], direct numerical simulation and Reynolds-averaged Navier-Stokes simulation, and as a reference for researchers. Of particular interest in the text are the detailed discussion, in Chapter 2, of vorticity, pressure, and the velocity gradient tensor, quantities useful for probing the results of a simulation, particularly when looking for coherent vortices and coherent structures. Chapters 4 and 5 feature an in-depth discussion of spectral subgrid-scale modeling. Although physical-space models are generally more readily applied, spectral models give insight into the requirements and limitations in subgrid-scale modeling and backscattering. A third special feature is the detailed discussion in Chapter 7, of large-eddy simulation of compressible flows previously only available in articles scattered throughout the literature. This will be of interest to those dealing with supersonic flows, combustion, astrophysics, and other related topics.

We develop, calibrate and test a dataset intended to drive global ocean hindcasts simulations of the last five decades. This dataset provides surface meteorological variables needed to estimate air-sea fluxes and is built from 6-hourly surface atmospheric state variables of ERA40. We first compare the raw fields of ERA40 to the CORE.v1 dataset of Large and Yeager (2004), used here as a reference, and discuss our choice to use daily radiative fluxes and monthly precipitation products extracted from satellite data rather than their ERA40 counterparts. Both datasets lead to excessively high global imbalances of heat and freshwater fluxes when tested with a prescribed climatological sea surface temperature. After identifying unrealistic time discontinuities (induced by changes in the nature of assimilated observations) and obvious global and regional biases in ERA40 fields (by comparison to high quality observations), we propose a set of corrections. Tropical surface air humidity is decreased from 1979 onward, representation of Arctic surface air temperature is improved using recent observations and the wind is globally increased. These corrections lead to a significant decrease of the excessive positive global imbalance of heat. Radiation and precipitation fields are then submitted to a small adjustment (in zonal mean) that yields a near-zero global imbalance of heat and freshwater. A set of 47-year-long simulations is carried out with the coarse-resolution (2° × 2°) version of the NEMO OGCM to assess the sensitivity of the model to the proposed corrections. Model results show that each of the proposed correction contributes to improve the representation of central features of the global ocean circulation.

It is well known that subgrid models such as Smagorinsky's cannot be used for the spatially growing simulation of the transition to turbulence of flat-plate boundary layers, unless large-amplitude perturbations are introduced at the upstream boundary: they are over-dissipative, and the flow simulated remains laminar. This is also the case for the structure-function model (SF) of Métais & Lesieur (1992). In the present paper we present a sequel to this model, the filtered-structure-function (FSF) model. It consists of removing the large-scale fluctuations of the field before computing its second-order structure function. Analytical arguments confirm the superiority of the FSF model over the SF model for large-eddy simulations of weakly unstable transitional flows. The FSF model is therefore used for the simulation of a quasi-incompressible (M∞ = 0.5) boundary layer developing spatially over an adiabatic flat plate, with a low level of upstream forcing. With the minimal resolution 650 × 32 × 20 grid points covering a range of streamwise Reynolds numbers Rex1 ε [3.4 × 105, 1.1 × 106], transition is obtained for 80 hours of time-processing on a CRAY 2 (whereas DNS of the whole transition takes about ten times longer). Statistics of the LES are found to be in acceptable agreement with experiments and empirical laws, in the laminar, transitional and turbulent parts of the domain. The dynamics of low-pressure and high-vorticity distributions is examined during transition, with particular emphasis on the neighbourhood of the critical layer (defined here as the height of the fluid travelling at a speed equal to the phase speed of the incoming Tollmien–Schlichting waves). Evidence is given that a subharmonic-type secondary instability grows, followed by a purely spanwise (i.e. time-independent) mode which yields peak-and-valley splitting and transition to turbulence. In the turbulent region, flow visualizations and local instantaneous profiles are provided. They confirm the presence of low- and high-speed streaks at the wall, weak hairpins stretched by the flow and bursting events. It is found that most of the vorticity is produced in the spanwise direction, at the wall, below the high-speed streaks. Isosurfaces of eddy viscosity confirm that the FSF model does not perturb transition much, and acts mostly in the vicinity of the hairpins.

Abstract The Cévennes–Vivarais Mediterranean Hydrometeorological Observatory (OHM-CV) is a research initiative aimed at improving the understanding and modeling of the Mediterranean intense rain events that frequently result in devastating flash floods in southern France. A primary objective is to bring together the skills of meteorologists and hydrologists, modelers and instrumentalists, researchers and practitioners, to cope with these rather unpredictable events. In line with previously published flash-flood monographs, the present paper aims at documenting the 8–9 September 2002 catastrophic event, which resulted in 24 casualties and an economic damage evaluated at 1.2 billion euros (i.e., about 1 billion U.S. dollars) in the Gard region, France. A description of the synoptic meteorological situation is first given and shows that no particular precursor indicated the imminence of such an extreme event. Then, radar and rain gauge analyses are used to assess the magnitude of the rain event, which was particularly remarkable for its spatial extent with rain amounts greater than 200 mm in 24 h over 5500 km2. The maximum values of 600–700 mm observed locally are among the highest daily records in the region. The preliminary results of the postevent hydrological investigation show that the hydrologic response of the upstream watersheds of the Gard and Vidourle Rivers is consistent with the marked space–time structure of the rain event. It is noteworthy that peak specific discharges were very high over most of the affected areas (5–10 m3 s−1 km−2) and reached locally extraordinary values of more than 20 m3 s−1 km−2. A preliminary analysis indicates contrasting hydrological behaviors that seem to be related to geomorphological factors, notably the influence of karst in part of the region. An overview of the ongoing meteorological and hydrological research projects devoted to this case study within the OHM-CV is finally presented.

The near- and far-field break-up and atomization of a water jet by a high-speed annular air jet are examined by means of high-speed flow visualizations and phase Doppler particle sizing techniques. Visualization of the jet's near field and measurements of the frequencies associated with the gas–liquid interfacial instabilities are used to study the underlying physical mechanisms involved in the primary break-up of the water jet. This process is shown to consist of the stripping of water sheets, or ligaments, which subsequently break into smaller lumps or drops. An entrainment model of the near-field stripping of the liquid is proposed, and shown to describe the measured liquid shedding frequencies. This simplified model explains qualitatively the dependence of the shedding frequency on the air/water momentum ratio in both initially laminar and turbulent water jets. The role of the secondary liquid break-up in the far-field atomization of the water jet is also investigated, and an attempt is made to apply the classical concepts of local isotropy to explain qualitatively the measurement of the far-field droplet size distribution and its dependence on the water to air mass and momentum ratios. Models accounting for the effect of the local turbulent dissipation rate in the gas on both the break-up and coalescence of the droplets are developed and compared with the measurements of the variation of the droplet size along the jet's centreline. The total flux of kinetic energy supplied by the gas per unit total mass of the spray jet was found to be the primary parameter determining the secondary break-up and coalescence of the droplets in the far field.

This paper summarizes some applications of ultrasonic vibrations regarding heat transfer enhancement techniques. Research literature is reviewed, with special attention to examples for which ultrasonic technology was used alongside a conventional heat transfer process in order to enhance it. In several industrial applications, the use of ultrasound is often a way to increase productivity in the process itself, but also to take advantage of various subsequent phenomena. The relevant example brought forward here concerns heat exchangers, where it was found that ultrasound not only increases heat transfer rates, but might also be a solution to fouling reduction.

▪ Abstract We review the mechanisms of steepening and breaking for internal gravity waves in a continuous density stratification. After discussing the instability of a plane wave of arbitrary amplitude in an infinite medium at rest, we consider the steepening effects of wave reflection on a sloping boundary and propagation in a shear flow. The final process of breaking into small-scale turbulence is then presented. The influence of those processes upon the fluid medium by mean flow changes is discussed. The specific properties of wave turbulence, induced by wave-wave interactions and breaking, are illustrated by comparative studies of oceanic and atmospheric observations, as well as laboratory and numerical experiments. We then review the different attempts at a statistical description of internal gravity wave fields, whether weakly or strongly interacting.

The chemical forms of zinc (Zn) in the Zn-tolerant and hyperaccumulator Arabidopsis halleri and in the non-tolerant and nonaccumulator Arabidopsis lyrata subsp. petraea were determined at the molecular level by combining chemical analyses, extended x-ray absorption spectroscopy (EXAFS), synchrotron-based x-ray microfluorescence, and muEXAFS. Plants were grown in hydroponics with various Zn concentrations, and A. halleri specimens growing naturally in a contaminated site were also collected. Zn speciation in A. halleri was independent of the origin of the plants (contaminated or non-contaminated) and Zn exposure. In aerial parts, Zn was predominantly octahedrally coordinated and complexed to malate. A secondary organic species was identified in the bases of the trichomes, which contained elevated Zn concentrations, and in which Zn was tetrahedrally coordinated and complexed to carboxyl and/or hydroxyl functional groups. This species was detected thanks to the good resolution and sensitivity of synchrotron-based x-ray microfluorescence and muEXAFS. In the roots of A. halleri grown in hydroponics, Zn phosphate was the only species detected, and is believed to result from chemical precipitation on the root surface. In the roots of A. halleri grown on the contaminated soil, Zn was distributed in Zn malate, Zn citrate, and Zn phosphate. Zn phosphate was present in both the roots and aerial part of A. lyrata subsp. petraea. This study illustrates the complementarity of bulk and spatially resolved techniques, allowing the identification of: (a) the predominant chemical forms of the metal, and (b) the minor forms present in particular cells, both types of information being essential for a better understanding of the bioaccumulation processes.

Over 20 global ocean tide models have been developed since 1994, primarily as a consequence of analysis of the precise altimetric measurements from TOPEX/POSEIDON and as a result of parallel developments in numerical tidal modeling and data assimilation. This paper provides an accuracy assessment of 10 such tide models and discusses their benefits in many fields including geodesy, oceanography, and geophysics. A variety of tests indicate that all these tide models agree within 2–3 cm in the deep ocean, and they represent a significant improvement over the classical Schwiderski 1980 model by approximately 5 cm rms. As a result, two tide models were selected for the reprocessing of TOPEX/POSEIDON Geophysical Data Records in late 1995. Current ocean tide models allow an improved observation of deep ocean surface dynamic topography using satellite altimetry. Other significant contributions include their applications in an improved orbit computation for TOPEX/POSEIDON and other geodetic satellites, to yield accurate predictions of Earth rotation excitations and improved estimates of ocean loading corrections for geodetic observatories, and to allow better separation of astronomical tides from phenomena with meteorological and geophysical origins. The largest differences between these tide models occur in shallow waters, indicating that the current models are still problematic in these areas. Future improvement of global tide models is anticipated with additional high‐quality altimeter data and with advances in numerical techniques to assimilate data into high‐resolution hydrodynamic models.

Abstract A 2D numerical model is proposed to simulate unsteady cavitating flows. The Reynolds‐averaged Navier–Stokes equations are solved for the mixture of liquid and vapour, which is considered as a single fluid with variable density. The vapourization and condensation processes are controlled by a barotropic state law that relates the fluid density to the pressure variations. The numerical resolution is a pressure‐correction method derived from the SIMPLE algorithm, with a finite volume discretization. The standard scheme is slightly modified to take into account the cavitation phenomenon. That numerical model is used to calculate unsteady cavitating flows in two Venturi type sections. The choice of the turbulence model is discussed, and the standard RNG k–εmodel is found to lead to non‐physical stable cavities. A modified k–εmodel is proposed to improve the simulation. The influence of numerical and physical parameters is presented, and the numerical results are compared to previous experimental observations and measurements. The proposed model seems to describe the unsteady cavitation behaviour in 2D geometries well. Copyright © 2003 John Wiley & Sons, Ltd.

In the tropical oceans the density mixed layer does influence the ocean‐atmosphere interactions, a feature that is fundamental to the development of energetic climatic events such as the El Niño‐Southern Oscillation (ENSO) phenomenon. The aim of this work is to take advantage of existing sea surface temperature (SST) and velocity together with future sea surface salinity (SSS) satellite‐derived data by assimilating these data in a primitive equation model in order to improve the modeled mixed layer. As satellite SSS data are not yet available, and to better analyze how assimilation works, we performed twin experiments in a context that renders assimilation conclusive with regard to real experiment. Two simulations using different forcing were used for the experiments. They have errors that are comparable to the ones between any simulation and real observations. The assimilation scheme is an adaptive version of the SEEK filter [ Pham et al., 1998 ] readily usable for assimilating real data. An assimilation experiment was conducted covering the 1997–1998 ENSO event to analyze the main characteristics of the actual mixed layer. Looking in particular at the relevance of SSS, SST, mixed layer depth, and barrier layer thickness, satellite‐derived data prove to be useful to better simulate the oceanic mixed layer. Velocity data are specially needed to control the zonal equatorial current. Interestingly, assimilation of surface‐only data still worked well below the mixed layer, and some improvements were detectable in terms of barrier layer thickness, even though the limit of the assimilation scheme was reached.

We present an experimental characterization of preferential concentration and clustering of inertial particles in a turbulent flow obtained from Voronoï diagram analysis. Several results formerly obtained from various data processing techniques are successfully recovered and further analyzed with Voronoï tesselations as the main single tool. We introduce a simple and nonambiguous way to identify particle clusters. We emphasize the maximum preferential concentration for particles with Stokes numbers around unity and the self-similar nature of clustering and we report new unpredicted results concerning clusters inner concentration dependence on Stokes number and global seeding density. Some of these experimental observations can be consistently interpreted in the context of the so-called sweep-stick mechanism. Finally, we stress the great potential of Voronoï analysis that offers important openings for new investigations of particle laden flows in terms, for instance, of simultaneous Lagrangian statistics of particle dynamics and local concentration field.