JDFP - Programme 2018

Nage phototactique de microalgues

F. Moisy, H. de Maleprade, K. Leptos, R. Goldstein : Université Paris-Sud. DAMTP, Université de Cambridge

Les algues unicellulaires flagellées de la famille des Chlamydomonas sont capables de nager en direction de la lumière (phototaxis positive), une stratégie bénéfique pour la photosynthèse. Si les mécanismes de réponse phototactique sont assez bien compris à l'échelle de la cellule unique, leur extension au cas de petites colonies multicellulaires reste mystérieuse : en l'absence de système nerveux centralisé, comment la réponse phototactique de chaque cellule se coordonne-t-elle pour conduire à une nage efficace de la colonie en direction de la lumière ? Des expériences réalisées sur des petites colonies de Gonium (assemblées de 16 cellules cousines de Chlamydomonas) permettent d'éclaircir les mécanismes de cette coordination.

Anneaux tourbillonnaires allongés

Ivan Delbende, Guillermo Artana, Pierre Audier, Denisse Sciamarella : Sorbonne Université, LIMSI ; FIUBA-CONICET, LFD, Buenos Aires.

Nous présentons des résultats expérimentaux et numériques sur le développement de jets synthétiques obtenus en sortie d'un orifice rectangulaire allongé modélisant la glotte. Des mesures par PIV synchronisée ont permis de reconstruire la dynamique tourbillonnaire 3D, et révèlent la séparation du jet en deux parties. Des expériences numériques sur des anneaux tourbillonnaires allongés sont présentées pour vérifier si le mécanisme de l'instabilité de Crow, originellement identifié dans les sillages d'avion, est à l'oeuvre dans ce phénomène glottique.

Metal pad rolling in small cylindrical Al-reduction cells : perturbation theory vs. multiphase MHD simulations

W. Herreman, C. Nore (1). J.-L. Guermond (2). L. Cappanera (3), N. Weber (4). : (1) LIMSI-CNRS, Université Paris-Sud, Université Paris-Saclay. (2)Texas A&M University. (3) Rice University (4) HZDR Dresden.

Dans des cellules d'electrolyse d'Aluminium, on observe l'apparition spontanée d'ondes de gravité tournantes à l'interface entre le cryolite (sel fondu dissolvant l'Alumina) et l'Aluminum liquide. Nous formulons un modèle perturbatif pour la déstabilisation de la cellule, qui inclue la dissipation visqueuse et magnétique habituellement ignorée. Le modèle théorique précis est comparé quantitativement avec des simulations numériques directes effectuées avec plusieurs codes numériques indépendants.

Convection in soft planets: when a plume can break a plate...

Anne Davaille (1), Suzanne Smrekar (2) : (1) FAST, CNRS / Univ. Paris-Sud / Univ. Paris-Saclay, Orsay, France (2) Jet Propulsion Laboratory, NASA/Caltech, Pasadena, USA

The cooling and evolution of a planet strongly depends on the convective regime prevailing in its « mantle » (e.g. the 2900 km-thick enveloppe of solid rocks under our feet on Earth). However, the conditions necessary for convection to generate quasi-continuous resurfacing on Earth, episodic and/or local resurfacing on Venus, or no resurfacing on Mars, remain strongly debated. The difficulty comes from the complexity of rocks rheology : viscous at high temperature and on long time-scale, brittle at low temperature and short time-scale. This « soft matter » behaviour can be recovered at much smaller scale in colloidal dispersions. We therefore have developed experiments of thermal and solutal convection in those systems. We found a new convective regime where hot rising plumes can break cold plates on the surface. Scaling analysis shows that it could operate on Venus today, but only on a hotter and younger Earth. It therefore could explain a number of observations on Venus and Earth.

Numerical modelling of turbulent convection in an asymmetric rough Rayleigh-Bénard cell

Mebarek Belkadi (1,3), Anne Sergent (1,2), Bérengère Podvin : (1) LIMSI, CNRS, Université Paris-Saclay. (2) Sorbonne Université

We present direct numerical simulations of turbulent buoyant convection over a rough heated plate placed in a Rayleigh Bénard cell. The roughness is introduced by a set of cubic obstacles regularly spaced, modelled by using an immersed boundary method. This study aims at clarifying interactions between the large scale circulation filling the box, plume emission and the enhancement of the heat transfer. The simulations are performed in a box-shaped cell at fixed Prandtl number (Pr = 4.38) with the Rayleigh number Ra ranging from 5x10e5 to 5x10e9. As expected, results show an enhancement of heat transfer measured by the Nusselt number Nu, depending on the relative sizes of the mean boundary-layer thicknesses and the obstacle height in agreement with previous experimental and numerical studies from the literature. By varying Ra, we investigate the successive regimes of the turbulent heat transport, from inactive roughness to a regime where the heat transfer relative increase is larger than the relative surface increase induced by roughness addition.

Stripes instability of an oscillating non Brownian iso-dense suspension of spheres

Y. L. Roht{1,2}, I. Ippolito{2}, J.P. Hulin{1}, D. Salin{1} et G. Gauthier{1} : 1 : Laboratoire FAST, Univ. Paris Sud, CNRS, Université Paris-Saclay, Orsay 2 : Universidad de Buenos-Aires, Facultad de Ingenieria, Grupo de Medios Porosos, Paseo Colon 850, 1063, Buenos Aires, Argentina

We analyze experimentally the behavior of a non-Brownian, iso-dense suspension of spheres submitted to periodic square wave oscillations of the flow in a Hele-Shaw cell of gap H. We do observe an instability of the initially homogeneous concentration in form of concentration variation stripes transverse to the flow. The wavelength of these regular spatial structures scales roughly as the gap of the cell and is independent of the particle concentration and of the period of oscillation. This instability requires large enough particle volume fractions ϕ≥0.25, a gap large enough compared to the spheres diameter (H/d≥8) and is observed in a rather broad range of periods of the square waves and amplitudes of the fluid displacement. We map the domain of existence of this instability in the space of the control parameters. The analysis of the concentration distribution across the gap shows that the instability is likely to be associated to the migration of the particles from the center towards the cell walls. In order to account for the main features of this stripes instability, we use the theory of longitudinal instability due to normal stresses difference and recent measurements of the first normal stresses difference dependence with particles concentration.

Etude expérimentale d’un film liquide tombant un plan incliné dans un canal fortement confiné

Yiqin Li, Sophie Mergui, Georg Dietze : FAST, CNRS, UPMC

An experimental setup has been developed to study the falling water film interacting with a counter-current gas flow under strong confinement. It consists in an inclined glass plate (150cm *27cm) placed on a massive framework with adjustable inclination angle in the range 0-20°. A 70cm long strongly confined rectangular channel is created at upstream with a gap height of 5.1mm. Chromatic Confocal Imaging is used to measure the liquid film thickness. A linear stability study without counter-current gas flow has been performed. It has showed that the confinement tend to stabilize the interfacial waves.

Stabilization of thin liquid films by confinement

G. Lavalle, Y. Li, S. Mergui, N. Grenier, G. Dietze : LIMSI, CNRS, Université Paris-Sud, Université Paris-Saclay FAST, CNRS, Université Paris-Sud, Université Paris-Saclay

We revisit the linear stability of a falling liquid film interacting with a counter-current gas flow. In this study, we describe the effect of confinement on the interfacial instabilities and show that the confinement can strongly stabilize the liquid film. For instance, we find that the critical Reynolds number can be increased by up to 30% at moderate relative confinement. The stabilization mechanism is due to the perturbed interfacial shear stress by reducing the inertial effects. Direct numerical simulations of gas-liquid flows within confined channels are also performed.

Simulation numérique de l'interaction onde de choc couche limite

Ismaïl BEN HASSAN SAÏDI et Christian TENAUD, Guillaume FOURNIER : LIMSI-CNRS, Université Paris-Sud, LMEE, Université d'Évry-Val-d'Essonne

Under certain circumstances (High Mach number, large shock angle...), the interaction between an incident shock wave impinging a boundary layer may create an unsteady separation bubble. This bubble, as well as the subsequent reflected shock, are known to oscillate in a low-frequency streamwise motion that can spread over several tenth of the boundary layer thickness.The origin of those oscillations, however still unclear, has been related either to the shedding of vortices in the mixing layer downstream of the separation, or to the turbulent structures in the incoming boundary layer. In this presentation, we aim to introduce the study of the shock wave boundary layer interaction currently performed at LIMSI-CNRS using the in-house parallel (MPI) Finite-Volume based DNS/LES CHORUS.

Stochastic Chaos in a Turbulent Swirling Flow

Davide Faranda, Yuzuru Sato, Brice Saint-Michel, Cecile Wiertel, Vincent Padilla, Berengere Dubrulle, Francois Daviaud : LSCE - CNRS - CEA Saclay - Université Paris Saclay

We report the experimental evidence of the existence of a random attractor in a fully developed turbulent swirling flow. By defining a global observable which tracks the asymmetry in the flux of angular momentum imparted to the flow, we can first reconstruct the associated turbulent attractor and then follow its route towards chaos. We further show that the experimental attractor can be modeled by stochastic Duffing equations, that match the quantitative properties of the experimental flow, namely the number of quasi-stationary states and transition rates among them, the effective dimensions, and the continuity of the first Lyapunov exponents. Such properties can neither be recovered using deterministic models nor using stochastic differential equations based on effective potentials obtained by inverting the probability distributions of the experimental global observables. Our findings open the way to low dimensional modeling of systems featuring a large number of degrees of freedom and multiple quasi-stationary states.

Footprint of air Turbulence on a water surface

S. Perrard 1,2, F. Moisy 1, M. Benzaquen 2, M. Rabaud 2 : 1: Université Paris-Sud, Laboratoire, 2: Ecole Polytechnique, LadHyX

Drainage de film de suspension

Emilie Dressaire, Alban Sauret : FAST, SVI/FAST

Lorsqu’un film mince de suspension s’écoule sur une surface solide, certaines particules peuvent rester piégées immobilisées sur le substrat. En effet, si l’épaisseur du film liquide devient comparable à la taille des particules, celles-ci déforment l’interface air/liquide, ce qui conduit à des interactions capillaires. Ces effets interfaciaux modifient le transport des particules et la dynamique du film liquide. Nous considérons ici le drainage gravitaire d’un film de suspension sur un plan incliné. En fonction de la concentration des particules, nous constatons que la dynamique de drainage ne peut pas être décrite par un modèle rhéologique continu, en raison du dépôt de particules sur le substrat.

Input-Output models for time-series analysis: an application to turbulent jets

Onofrio Semeraro (1), Francois Lusseyran (1), Luc Pastur (2) : (1) LIMSI, (2) ENSTA

The temporal dynamics associated with the coherent structures in a turbulent jet has been the subject of a recent publication from our group. These coherent structures -- often called "wavepackets" -- exhibit a remarkable spatiotemporal organization, despite turbulence. The characterization of their temporal dynamics was attempted by means of nonlinear statistical tools and classic tools from system identification. We found that linear models can reproduce much of the convective input-ouput behavior, while nonlinear, output-only models can reproduce only a small portion of the dynamical features of the system. This apparent contradiction is mainly due to the presence of an input in the model. In that sense, it is questionable whether output-only models are appropriated or not for such an analysis. To answer this question, we are currently reviewing our results by means of a recent algorithm proposed by Brunton et al. (Nature Comm. 8, 2017). The algorithm combines delay embedding and Koopman theory to decompose nonlinear dynamics into a linear model driven by a low-energy forcing extracted directly from the time-series. Preliminary results and open questions will be proposed during the presentation.

Influence de perturbations numériques sur la formation et le comportement d'un bulbe de décollement laminaire par simulations des grandes échelles.

Tommy Rigall, Benjamin Cotté, Philippe Lafon : IMSIA - UMR EDF/CNRS/CEA/ENSTA 9219, 91120, Palaiseau

Le contexte actuel de démocratisation des machines tournantes de petite taille telles que les drones ou les éoliennes domestiques pose la question des nuisances sonores qu’elles engendrent. Dans ce cadre, étudier et comprendre la structure de l’écoulement autour de ces structures est essentiel afin de prévoir le champ acoustique. La taille ré- duite de ces systèmes implique de faibles nombres de Reynolds, de l’ordre de 10^4 à 10^5 . A ces régimes, un bulbe de décollement laminaire est susceptible d’apparaître sur les pales. Au vu de l’influence de ce bulbe sur les caractéristiques aérodynamiques des ailes, capturer précisément les positions de séparation, de transition laminaire-turbulent et de recollement de la couche limite est important. Dans le cadre du présent travail, des calculs directs de bruit d’un profil NACA0012 à bas nombres de Reynolds et à différents angles d’incidence ont été effectués afin d’obtenir simultanément et précisément le bulbe de décollement laminaire ainsi que le champ acoustique associé. Pour cela, le code de calcul Code_Safari développé par EDF R&D a été utilisé. Il permet de résoudre les équations de Navier-Stokes compressibles en effectuant des simula- tions des grandes échelles par filtre de relaxation. Des schémas de différences finies d’ordre élevé ainsi que la méthode de maillage Chimère sont employés. La naissance, l’intermittence et l’éclatement du bulbe de décollement laminaire ont été observés et ses effets sur les performances aérodynamiques du profil seront présentés. Il sera également montré que la présence ou non du bulbe est très dépendante de l’amplification des perturbations numériques et notamment de l’ordre du schéma numérique et des zones d’interpolation employés. En particulier, la quantification de ces perturbations numériques, leur influence sur la formation du bulbe de décollement ainsi que sur les performances aérodynamiques du profil seront présentées. Il sera mis en avant que ces perturbations numériques font office de turbulence en amont du profil, dont l’intensité n’est pas contrôlée. La nature de l’instabilité provoquant la transition laminaire-turbulent de la couche limite autour du profil sera également abordée.

Influence of the thermo-elastic stresses on the onset of convective motions in carbopol.

Nicolo Sgreva, Anne Davaille : FAST, CNRS / Univ. Paris-Sud

In non-newtonian, yield-stress fluid such as Carbopol, it is assumed that deformations and flow can occur only if local stresses are larger than a critical value, the yield stress. However, this model is not sufficient to describe how to get there and build up the stresses, i.e. the transition from the solid to the fluid state, in natural convection experiments where a fluid initially at rest is heated from below. We explore here if and how thermo-elastic stresses can help reaching the yield stress. We developed a numerical simulation in order to investigate the development of thermo-elastic stresses in two different experimental situations: (1) the classical Rayleigh-Benard geometry where a steady-state linear temperature gradient is established between two horizontal boundaries heated form below and cooled form above, and (2) a localized source heating with a prescribed temperature history that heats the fluid from below. Both situations were studied in 2D cartesian and axisymmetric geometry with no-slip boundary conditions and with different combinations of no-slip and free-slip boundary conditions in order to replicate the experimental setups. Our results show that thermo-elastic stresses can indeed overcome the yield stress. It is therefore important to consider the first deformation stage (i.e. the elastic response) for an elasto-viscoplastic fluid subject to a temperature difference.

Formation des vermiculations dans les cavités souterraines. Application à la préservation des peintures pariétales de la grotte de Lascaux.

P. Freydier, F. Doumenc, B. Guerrier, J.P. Hulin, J. Martin, S. Mergui : FAST, UMR 7608

Climatologie en grottes peu profondes

F.Doumenc, P.Freydier, B.Guerrier, J.P.Hulin, J.Martin, S.Mergui,P.Y.Jeannin : FAST (UPSud-CNRS), ISSKA, Suisse

Statistiques conditionnelles des structures cohérentes en convection turbulente de Rayleigh-Bénard

Andrés Castillo-Castellanos, Anne Sergent, Bérengère Podvin & Maurice Rossi : LIMSI, CNRS, Université Paris-Saclay, F-91405 Orsay, Sorbonne Universités, UPMC Univ Paris 06, UMR 7190, Inst Jean Le Rond dAlembert, F-75005, Paris, Sorbonne Universités, UPMC Univ Paris 06, UFR d’Ingénierie, F-75005, Paris, CNRS, UMR 7190, Inst Jean Le Rond dAlembert, F-75005, Paris,

A partir de simulations numériques directes des équations de Boussinesq en régime faiblement turbulent et réalisées sur des temps physiques longs (plusieurs milliers d’unités de temps convectives), il est possible d’identifier sur une certaine plage de nombres de Rayleigh, deux régimes intermittents d’écoulement. Le premier régime est caractérisé par des renversements spontanés et souvent consécutifs de la circulation à grande échelle, tandis que le deuxième régime, plus rapide et moins prédominant, est caractérisé par la cessation prolongée de cette circulation. Nous adoptons une approche statistique conditionnelle pour étudier les deux régimes séparément, sur plusieurs centaines de renversements observés. Cette approche a permis préalablement de mettre en évidence un cycle générique décrivant énergétiquement les renversements standards en terme d’énergie cinétique et énergie potentielle disponible grâce à une remise à l’échelle temporelle des renversements.

Écoulement tournant avec suface libre: Quel écoulement de base pour les polygones tournants ?

Wen Yang, Ivan Delbende, Yann Fraigneau, Laurent Martin Witkowski : Sorbonne Université, LIMSI-CNRS, FAST

Des motifs polygonaux tournants sont souvent observés au niveau de la surface libre fortement déformée d'un liquide au sein d'une cuve cylindrique à fond tournant. Un modèle simplifié proposé dans la littérature comme écoulement de base pour étudier cette instabilité est le vortex de Rankine. Nous voudrions trouver un écoulement de base axisymétrique plus réaliste par des calculs DNS et des mesures expérimentales. L’écoulement se compose d’une rotation solide au centre et d'une circulation méridienne à la périphérie, souvent négligée. Les résultats numériques montrent que cette circulation est entourée par quatre couches limites dans le plan méridien. En particulier, l'écoulement radial rentrant de la couche de surface peut causer une survitesse de la vitesse azimutale, plus ou moins marquée selon que les simulations sont effectuées en 2D ou en 3D. Les mesures expérimentales confirment globalement les solutions numériques, avec cependant des petits écarts locaux sur l’écoulement méridien. En particulier, il semble que le modèle de Rankine rende bien compte de la distribution de vitesse azimutale, ce qui pose la question du rôle de la modélisation numérique de l'interface sur la structure de la couche de surface.

Surface libre dans un écoulement tournant : de l'idéalité numérique à la réalité expérimentale

Antoine Faugaret, Laurent Martin Witkowski, Yohann Duguet, Yann Fraigneau : Sorbonne Université, LIMSI-CNRS

Dans cette brève présentation, nous parlerons des instabilités survenant dans un écoulement tournant en présence d'une surface libre non deformée. En particulier, nous nous intéresserons à la condition limite à la surface et l'influence de celle-ci sur les écarts entre prévisions numériques et mesures expérimentales des Reynolds critiques et des modes d'instabilités.

Du contrôle d'écoulement par apprentissage automatique

Guy Cornejo Maceda, Bernd R. Noack, François Lusseyran, Marek Morzynski, Luc Pastur, Nan Deng : LIMSI/CNRS

Nous explorons les capacités des algorithmes génétiques à générer un modèle de contrôle en boucle fermée efficace pour diminuer la traîner dans une configuration d'écoulement modèle: le PinBall Fluidique. L'analyse du modèle, appris automatiquement, à vocation à nous révéler les propriétés de la dynamique de l'écoulement sous contrôle.

Etats métastables dans l'écoulement de Couette plan

David de Souza, Romain Monchaux : IMSIA

Nous étudions l'écoulement de Couette plan à des nombres de Reynolds proches de la transition à la turbulence. En réalisant des abaissements lents du nombre de Reynolds, nous observons une dynamique de retrait de la turbulence conduisant à la relaminarisation de l'écoulement. Cette relaminarisation fait apparaître des plateaux métastables qui peuvent être assez longuement visités par l'écoulement.

Bouncing dynamics of a sphere on a textured wall

Thibault Chastel, Angela Phor, Philippe Gondret, Anne Mongruel : Laboratoire FAST (UPSud/CNRS Orsay) & PMMH (UPMC/CNRS/ESPCI Paris)

In the present work, we investigate how the collisional process of a solid sphere onto a wall in a viscous liquid is influenced by a wall texture, with a special attention given to the bouncing transition. The sphere motion is measured by an interferometer with a high spatial and time resolution. The detailed bouncing dynamics is modeled with different force terms with special attention dedicated to fluid dissipation.

Mesure des interactions gouttes - plasma par mesures 4D.

Junghwa Yi, Antoine Renaud et Laurent Zimmer : CNRS, CentraleSupélec, Université Paris-Saclay

Ignition of sprays remains an important issue. The final faith of an event strongly depends on local conditions (positions and sizes of droplets as well as velocities). In the present paper, measurements performed in a volume (cylinder of 10 mm diameter) around a spark induced by focussing a Nd:YAG laser (8ns - 532nm) are discussed. Droplets of dodecane (Sauter Mean Diameter 20 microns) are convected at atmospheric pressure and ambiant temperature through a laminar counter-flow burner with velocities around 1 m/s. Defocussed images of droplets are taken using a continuous Nd:YAG laser and a 10kHz camera. Position and velocities of each droplets are measured by a 3D-optical flow approach whereas diameters are obtained using an IInterferometric Particle Imaging technique. Results show different behaviors of laser - droplets - plasma interactions, depending on initial position of droplets with respect to the focal point and laser induced plasma. In some cases, plasma may be created few mm before the focal point. On the other hand, droplets may be vaporized, split into smaller droplets or their trajectories be substantially changed due to the shock-wave and the rarefaction wave.

Adjoint-based mean flow reconstruction of laminar separated bubbles based on PIV data

Nicolò Fabbiane, Olivier Marquet and Benjamin Leclaire : ONERA - DAAA (MAPE)

In this work, experimental techniques and numerical simulations cooperate to produce reliable measurements of the flow. An adjoint-based data-assimilation procedure is used to match experimental data and RANS simulations. Boundary conditions and turbulent viscosity production -- i.e. the assimilated variables -- are tailored to minimise a quadratic cost function based on the difference between RANS simulations and a mean flow field obtained by PIV. The gradient of the cost function with respect to turbulent viscosity and boundary condition is determined by solving a steady adjoint problem and it is used to obtain a descent direction in a BFGS minimisation algorithm.

Miroirs et cristaux temporels pour les ondes (de surface)

V. Bacot 1, M. Labousse 2, E. Eddi 3, G. Durey 2, M. Fink 4, E. Fort 4 : 1) Laboratoire d'hydrodynamique de l'École polytechnique (LadHyX), 2) Laboratoire Gullivers, ESPCI Paris, 3) Laboratoire PMMH , ESPCI Paris, 4) Institut Langevin, ESPCI Paris

L’universalité des phénomènes ondulatoires classiques est largement décrite par la structure de l’équation de d’Alembert. Dans cette dernière, les variables spatiales et temporelles jouent des rôles analogues. Cette analogie bien connue suggère la possibilité de réaliser dans le domaine temporel des équivalents de phénomènes ondulatoires spatiaux standard, tels que la réflexion d’une onde par un miroir, en remplaçant la variation spatiale des propriétés du milieu qui y est mise en jeu par une variation temporelle équivalente. Nous avons démontré expérimentalement, en utilisant les ondes à la surface d’un liquide comme système ondulatoire modèle, la pertinence de cette approche et montré théoriquement sa généralisation à tout type d’onde. L’équivalent temporel d’un miroir conduit à la génération instantanée, depuis tout l'espace, d’une onde retournée temporellement c’est-à-dire dont l’évolution temporelle est inversée par rapport à l’onde initiale. Une version monofréquentielle de ce miroir a été réalisée et permet notamment de revisiter l’instabilité de Faraday à la surface d’un bain de liquide vibré comme l’équivalent temporel d’un miroir de Bragg. Ce dernier est également associé à un retournement temporel monochromatique, étroitement lié au phénomène appelé conjugaison de phase en optique. Dans ces deux équivalents temporels, les limites de l’analogie entre temps et espace induites par la causalité révèlent des différences très intéressantes entre les concepts présentés et leurs équivalents spatiaux.

Granular compaction by fluidization

G. Gauthier and P. Gondret : Laboratoire FAST, Université Paris Sud, CNRS, Université Paris Saclay

We presented experimental results on the compaction of a granular pile immersed in a viscous fluid submit to a continuous or bursting upward flow. An initial fluidized bed leads to well reproduced initial loose packing by the settling of grains when the high enough continuous upward flow is turned off. When the upward flow is then turned on again, we record the dynamical evolution of the bed packing. For a low enough continuous upward flow, below the critical velocity of fluidization, a slow compaction dynamics is observed. Strikingly, a slow compaction can be also observed in case of “fluidization taps” with bursts of fluid velocity higher than the critical fluidization velocity. The different compaction dynamics will be discussed when varying the different control parameters of these “fluidization taps”.

Internal dynamics in the front of free-surface viscoplastic surges : Development of a consistent thin-layer model and experimental confrontation

P. Freydier (1), G. Chambon (1), M. Naaim (1), J.P Vila (2 ) : (1) Irstea/Université Grenoble Alpes, (2) IMFT