Cette thèse présente la Simulation des Grandes Echelles (LES) de l’injection, de la pulvérisation et de la cavitation dans un injecteur pour les applications liées aux moteurs à combustion interne. Pour la modélisation de l’atomisation, on utilise le modèle ELSA (Eulerian Lagrangian Spray Atomization). Le modèle résout la fraction volumique du combustible liquide ainsi que la densité de surface d’interface liquide-gaz pour décrire le processus complet d’atomisation. Dans cette thèse, l’écoulement à l’intérieur de l’injecteur est également pris en compte pour une étude ultérieure de l’atomisation. L’étude présente l’application du modèle ELSA à un injecteur Diesel typique, à la fois dans le contexte de RANS et de LES.Le modèle est validé à l’aide de données expérimentales disponibles dans Engine Combustion Network (ECN). Le modèle ELSA, qui est normalement conçu pour les interfaces diffuses (non résolues), lorsque l’emplacement exact de l’interface liquide-gaz n’est pas pris en compte, est étendu pour fonctionner avec une formulation de type Volume of Fluid (VOF) de flux à deux phases, où l’interface est explicitement résolu. Le couplage est réalisé à l’aide de critères IRQ (Interface Resolution Quality), qui prennent en compte à la fois la courbure de l’interface et la quantité modélisée de la surface de l’interface. Le modèle ELSA est développé en premier lieu en considérant les deux phases comme incompressibles. L’extension à la phase compressible est également brièvement étudiée dans cette thèse. Il en résulte une formulation ELSA compressible qui prend en compte la densité variable de chaque phase. En collaboration avec l’Imperial College de Londres, la formulation de la fonction de densité de probabilité (PDF) avec les champs stochastiques est également explorée afin d’étudier l’atomisation. Dans les systèmes d’injection de carburant modernes, la pression locale à l’intérieur de l’injecteur tombe souvent en dessous de la pression de saturation en vapeur du carburant, ce qui entraîne une cavitation. La cavitation affecte le flux externe et la formulation du spray. Ainsi, une procédure est nécessaire pour étudier le changement de phase ainsi que la formulation du jet en utilisant une configuration numérique unique et cohérente. Une méthode qui couple le changement de phase à l’intérieur de l’injecteur à la pulvérisation externe du jet est développée dans cette thèse. Ceci est réalisé en utilisant le volume de formulation de fluide où l’interface est considérée entre le liquide et le gaz; le gaz est composé à la fois de vapeur et d’airambiant non condensable. / This thesis presents Large Eddy Simulation (LES) of fuel injection, atomization and cavitation inside the fuel injector for applications related to internal combustion engines. For atomization modeling, Eulerian Lagrangian Spray Atomization (ELSA) model is used. The model solves for volume fraction of liquid fuel as well as liquid-gas interface surface density to describe the complete atomization process. In this thesis, flow inside the injector is also considered for subsequent study of atomization. The study presents the application of ELSA model to a typical diesel injector, both in the context of RANS and LES. The model is validated with the help of experimental data available from Engine Combustion Network (ECN). The ELSA model which is normally designed for diffused (unresolved) interfaces, where the exact location of the liquid-gas interface is not considered, is extended to work with Volume of Fluid (VOF) type formulation of two phase flow, where interface is explicitly resolved. The coupling is achieved with the help of Interface Resolution Quality (IRQ) criteria, that takes into account both the interface curvature and modeled amount of interface surface. ELSA model is developed first considering both phases as incompressible, the extension to compressible phase is also briefly studied in this thesis, resulting in compressible ELSA formulation that takes into account varying density in each phase. In collaboration with Imperial College London, the Probability Density Function (PDF) formulation with Stochastic Fields is also explored to study atomization. In modern fuel injection systems, quite oftenthe local pressure inside the injector falls below the vapor saturation pressure of the fuel, resulting in cavitation. Cavitation effects the external flow and spray formulation. Thus, a procedure is required to study the phase change as well as jet formulation using a single and consistent numerical setup. A method is developed in this thesis that couples the phase change inside the injector to the external jet atomization. This is achieved using the volume of fluid formulation where the interface is considered between liquid and gas; gas consists of both the vapor and non condensible ambient air.

Une revue bibliographique présente le problème de l'érosion et les études antérieures. Une campagne d'essais est menée à l'aide de deux érodimètres (JET/HET) sur un ensemble de sols naturels, compactés suivant le protocole Proctor. L'interprétation des essais est effectuée avec les méthodes de la littérature. Les paramètres obtenus dépendent de l'appareil utilisé avec un mécanisme identique. Une méthode d'interprétation basée sur l'énergie hydraulique dissipée est développée pour les deux érodimètres. Pour une masse érodée donnée, les énergies sont similaires. L'indice de résistance à l'érosion proposé permet de regrouper sur un meme plan les essais. Une analyse statistique aboutit à l'expression de l'indice de résistance à l'érosion en quelques paramètres au voisinage de l'optimum Proctor. Une étude de la variabilité de la sensibilité à l'érosion met en évidence l'influence du degré de saturation et la corrélation de la perméabilité avec l'indice de résistance à l'érosion.

L'objectif de cette thèse est de modéliser et de développer des outils numériques adaptés à l'étude de l'écoulement des eaux souterraines ainsi que la propagation des polluants en milieux poreux. La motivation de ce travail est un benchmark du GDR Momas et de l'Andra pour la simulation de la propagations 3-D des radionucléides autour d'un stockage profond de déchets nucléaires. Premièrement on a construit une nouvelle méthode d'éléments finis mixtes sur un maillage formé d'hexaèdres généraux. La convergence de la méthode est prouvée et confirmée par des tests numériques. Deuxièment, nous présentons une méthode de discrétisation en temps pour une équation d'advection telle que des pas de temps différents sont utilisés dans différents sous-domaines afin de prendre en compte les hétérogèneités.<br />Enfin une méthode numérique pour le calcul de transport de contaminants est proposée. Les techniques précédentes sont implémentées en 3-D et des résultats numériques sont présentés sur le benchmark 3-D champ lointain du GDR Momas et de l'Andra.

Le contexte général de cette thèse s'inscrit dans le cadre de la modélisation eulérienne instationnaire des écoulements turbulents anisothermes gaz - particules. La modélisation de ces écoulements est cruciale pour de nombreuses applications industrielles et pour la prédiction de certains phénomènes naturels. Par exemple, la combustion diphasique dans les moteurs automobiles et aéronautiques est précédée par l'injection et la dispersion de carburant liquide dans la chambre de combustion. Les phénomènes mis en jeu exigent alors une prédiction locale tenant compte du caractère instationnaire de l'écoulement turbulent et de la présence de géométries complexes. De plus, de nombreuses études expérimentales et numériques récentes ont mis en évidence le rôle prépondérant de l'inertie des particules sur les mécanismes de dispersion et de concentration préférentielle en écoulement turbulent. Ceci rend donc indispensable la prise en compte de ces mécanismes dans la modélisation diphasique. Au cours de ce travail de thèse, une approche eulérienne locale et instantanée a été développée pour prédire les écoulements gaz-particules anisothermes et turbulents. Elle est basée sur l'approche statistique du Formalisme Eulérien Mésoscopique (MEF) introduite par Février et al. (JFM, 2005). Cette approche a été ici étendue aux variables thermiques pour la prise en compte du caractère anisotherme de l'écoulement. Cette approche a été ensuite utilisée dans le cadre de la méthode des moments (Kaufmann et al., JCP, 2008), et un système d'équations locales et instantanées pour la phase dispersée a été proposé. La modélisation au premier ordre exige la fermeture des moments de second ordre apparaissant dans les équations de la quantité de mouvement et de l'énergie. La proposition de telles relations constitutives fait l'objet d'une partie de la thèse. Afin de fournir une méthode capable de prédire le comportement local, instantané et anisotherme de la phase dispersée dans des configurations `a une échelle réaliste, les équations pour la phase dispersée ont été filtrées et une modélisation aux grandes échelles (LES) est effectuée. Cette modélisation étends, par la prise en compte des variables thermiques, le travail de Moreau et al. (FtaC, 2010) sur l'approche LES Euler-Euler en conditions isothermes. L'approche complète est enfin appliquée aux résultats de simulation numérique d'un jet plan turbulent gazeux froid, chargé en particules, dans une turbulence homogène isotrope chaude monophasique. / The aim of this thesis is to provide an Eulerian modeling for the dispersed phase interacting with unsteady non-isothermal turbulent flows. The modeling of these flows is crucial for several industrial applications and for predictions of natural events. Examples are the combustion chambers of areo engines where the combustion is preceded by the injection and dispersion of liquid fuel. The prediction of such phenomena involves a local modeling of the mixture for taking into account the unsteady behavior of the turbulent flow and the presence of complex geometries. Moreover, many experimental and numerical studies have recently highlighted the significant role of the particle inertia on the mechanisms of dispersion and preferential concentration. Accounting for such mechanisms is therefore essential for modeling the particle-laden turbulent flows. In this thesis, a local and instantaneous Eulerian approach able to describe and to predict the local behavior of inertial particles interacting with non-isothermal turbulent flows has been developed. It is based on the statistical approach known as Mesoscopic Eulerian formalism (MEF) introduced by Février et al. (JFM, 2005). The statistical approach has been extended to the thermal quantities in order to account for the non-isothermal conditions into the modeling. This formalism is then used in the framework of the moment approach (Kaufmann et al., JCP, 2008) and a system of local and instantaneous equations for the non-isothermal dispersed phase has been suggested. The first order modeling requires to close second-order moments appearing in momentum and energy equations. The proposal of such constitutive relations makes the object of a part of this study. In order to provide an Eulerian approach usable in real configurations at industrial scale, the equations of the dispersed phase are filtered and the approach developed in the framework of the Large-Eddy Simulations. From the work of Moreau et al. (FTaC, 2010), the Eulerian-Eulerian LES approach is then extended to non-isothermal conditions. The whole modeling is then a priori tested against numerical simulations of a cold planar turbulent particle-laden jet crossing a homogeneous isotropic decaying hot turbulence.

Ce travail présente des modèles d'écoulement en temps continu de flux de trafic routier. En premier lieu, il s'agit de modèles microscopiques de poursuite. Un modèle par systèmes d'équations différentielles couplées est proposé, basé sur le temps inter-véhiculaire. Ce modèle intègre un temps de réaction et des possibilités d'anticipation pour chaque véhicule. Les paramètres sont estimés par maximum de vraisemblance dans un modèle statistique à deux niveaux. Des simulations permettent de caractériser le comportement d'une file de véhicules. Dans une approche stochastique, un modèle d'évolution de la distance inter-véhiculaire est étudié à l'aide du processus Markovien de saut zero-range. L'introduction d'un temps de réaction tend à produire des ondes cinématiques. D'autre part, un modèle d'écoulement de trafic par le processus Markovien de saut des misanthropes est proposé. Il s'agit d'une modélisation au niveau mésoscopique, adaptée à la simulation de flux de trafic sur un réseau / This work presents different continuous-time traffic flow models. Microscopic models are considered first. A model by coupled differential equation system is proposed, based on the time gap. It incorporates a reaction time parameter and some anticipation possibilities, for each vehicle. The parameters are estimated by maximum likelihood over a two-level statistical model. Simulations allow to characterise the behaviour of a vehicles line. In a stochastic approach, a model of the distance gap evolution is studied with a zero-range process. The introduction of a reaction time parameter produces kinematics waves. On the other hand, traffic flow model by a misanthropes process is proposed. It is a mesoscopic approach, adapted to the simulation of traffic flow on a network

La rupture de la plaque d'athérome carotidienne est la première cause des infarctus cérébraux. Pour prévenir ces accidents, l'endartérectomie carotidienne est le traitement le plus utilisé. La vulnérabilité de la plaque est en relation avec les efforts que le sang applique sur la plaque. Ces actions sont différentes suivant les propriétés constitutives, mécaniques et géométriques de la plaque. Plusieurs auteurs ont développé des modèles numériques de la plaque d'athérome carotidienne à partir desquels une analyse mécanique a permis de caractériser les déformations et les contraintes en lien avec la rupture de la plaque. Néanmoins, les caractéristiques d'une plaque vulnérable sont encore mal connues. Dans ce manuscrit, un modèle numérique de plaque d'athérome carotidienne est développé en interaction fluide-structure dans le but mieux comprendre comment les actions mécaniques du sang sur la plaque sont affectées par les propriétés mécaniques et géométriques de la plaque. Plusieurs résultats sont en concordance avec la littérature : la vulnérabilité de la plaque est associée à la sévérité de sténose et à l'épaisseur de la chape fibreuse. De plus une analyse de l'écoulement du sang, de la déformation de la plaque et des contraintes dans la plaque révèle que les effets de l'écoulement du sang sont amplifiés si la plaque est courte, si la pente en amont de sténose est raide ou si la morphologie de la plaque est irrégulière et asymétrique. Ces résultats offrent de nouvelles perspectives dans la compréhension de la vulnérabilité de la plaque. / The rupture of carotid atheromatous plaques is the major cause of cerebrovascular thromboembolic events such as strokes and ischemic attacks. To prevent this issue, carotid endarterectomy is the preferred treatment. The vulnerability of the plaque is related the mechanical action of the blood onto the plaque. This action is different according to the plaque morphology, the plaque constitution and the mechanical properties of the constituents. Several authors developed computational models to perform mechanical analyses for carotid atherosclerotic plaques and to identify critical mechanical descriptors as stresses or strains related to plaque rupture. But the question of which plaque characteristics affect the plaque rupture is not closely elucidated. In this manuscript a fluid structure interaction model is developed, questioning how the mechanical action of the blood onto an atheromatous plaque is affected by the mechanical and geometrical properties of the plaque. Many results are in agreement with the literature: the vulnerability of atheromatous plaques is related to the degree of severity of the endoluminal stenosis and the thickness of the fibrous cap. Moreover the resulting flow patterns, wall shear stresses, plaque deformations and stresses in the fibrous cap reveal that the effects of the blood flow are amplified if the plaque is short, if the slope upstream stenosis is steep or if the plaque morphology is irregular and asymmetric. These results offer new perspectives for understanding the vulnerability of plaques.

Deux stratégies de contrôle ont été successivement mises en œuvre pour réduire la traînée et l’énergie cinétique turbulente d’un canal plan en régime turbulent (Re[tau]=180) par soufflage/aspiration aux parois. L’objectif principal était de prouver qu’une simulation aux grandes échelles (LES) pouvait être utilisée de manière pertinente comme modèle réduit des équations de Navier-Stokes et ainsi diminuer fortement les coûts numériques. Une approche heuristique dite de contrôle par opposition a d’abord été employée. Les résultats montrent que l'efficacité énergétique est maximale pour une position du plan de détection différente de celle qui correspond au maximum de réduction de traînée. Par ailleurs, nos résultats confirment que la réduction de traînée diminue avec l'augmentation du nombre de Reynolds. Par la suite, une procédure de contrôle optimal a été utilisée en considérant différentes fonctionnelles objectif (traînée, énergie cinétique au temps terminal, énergie cinétique moyen). Pour Re[tau]=100, le contrôle est parvenu à relaminariser complètement l’écoulement (réduction de traînée de l'ordre de 50 %) en prenant comme fonctionnelle coût l’énergie cinétique au temps terminal. Pour cette même fonctionnelle coût, une réduction importante de traînée de l'ordre de 55 % est encore obtenue à Re[tau] =180 mais sans atteindre la relaminarisation. Nos résultats confirment que pour minimiser la traînée de l’écoulement, il est plus efficace de considérer comme objectif l’énergie cinétique que directement la traînée. Enfin, il est essentiel pour la convergence de la minimisation que le système optimal soit résolu sur un horizon temporel suffisamment long / Two control strategies were successively implemented to reduce the drag and the turbulent kinetic energy of a plane channel flow in turbulent regime (Re[tau]=180). Wall transpiration (unsteady blowing/suction) with zero net mass flux is used as the control. The main objective was to prove that a large eddy simulation (LES) could be relevant as a reduced-order model of the Navier-Stokes equations and thus strongly reduce the numerical costs. A heuristic approach known as opposition control was initially employed. The results show that the energetic efficiency is maximum for a position of the detection plane different from that which corresponds to the maximum of drag reduction. In addition, our results confirm that the drag reduction decreases with the increase of the Reynolds number. Then, an optimal control procedure was used by considering different cost functional (drag, terminal turbulent kinetic energy, mean turbulent kinetic energy). At Re[tau] =100, control managed to fully relaminarize the flow (drag reduction of about 50%) by considering as cost functional the terminal kinetic energy. For this same cost functional, an important drag reduction of about 55% is still obtained at Re[tau] =180 but without reaching the relaminarization. Our results show that to minimize the flow drag, it is more effective to consider the kinetic energy as cost functional than directly the drag. Lastly, it is essential for the convergence of the minimization that the optimality system is solved on a sufficiently long time horizon

L’objet de cette thèse est de développer des méthodologies permettant d’identifier la distribution spatiale des valeurs de perméabilité dans des échantillons de roches. Nous avons tout d’abord développé en laboratoire des expériences d’injection de fluide miscible très visqueux dans des échantillons initialement saturés par une saumure peu visqueuse. Pendant l’injection, l’évolution au cours du temps de la pression différentielle entre les deux faces de l’échantillon a été enregistrée par des capteurs de pression. En outre, des mesures scanner ont fourni une carte 3D de la porosité ainsi que des cartes 3D décrivant la distribution spatiale des concentrations dans l’échantillon à différents temps. Nous avons mis en place une méthode d’interprétation donnant directement le profil 1D de la perméabilité le long de la direction d’écoulement à partir de la pression différentielle mesurée au cours du temps. Cette méthode a été validée numériquement et expérimentalement. Puis, afin d’affiner la description de l’agencement des valeurs de perméabilité dans l’échantillon, c’est à dire d’obtenir un modèle 3D de perméabilité représentatif de l’échantillon, nous avons développé une méthodologie itérative de calage des pressions et des concentrations. Cette méthode passe par deux étapes : une optimisation simple pour capturer l’hétérogénéité dans la direction de l’écoulement et une optimisation complexe pour capturer l’hétérogénéité transverse. Cette méthode a été validée à partir de tests numériques. La méthode a été appliquée à deux des expériences d’injection de fluide visqueux. Nous avons pu alors déterminer des modèles de perméabilité capables de reproduire assez bien les données de pression et de concentration acquises pendant l’injection / The objective of this study is to develop new methodologies to identify the spatial distribution of permeability values inside the heterogeneous core samples. We developed laboratory viscous miscible displacements by injecting high viscosity glycerin into the core samples initially saturated by low viscosity brine. The pressure drop across the samples was measured as a function of time until breakthrough. Meanwhile, CT scan measurements provided a 3D porosity map plus several 3D maps of concentration distribution inside the core samples at different times. A simple permeability mapping technique was developed deducing a one-dimensional permeability profile along the flow direction from the measured pressure drop data. The method was validated with both numerical and laboratory experiments. To go beyond one-dimensional characterization of permeability into cores, we developed an iterative process for matching pressure and concentration data. This method consisted of two steps: a simple optimization for capturing the permeability heterogeneity along the flow direction axis and a complex optimization for capturing transversal permeability heterogeneities. The methodology was validated by numerical data. It was also applied to the data collected from two laboratory viscous miscible displacements. We showed that the final 3D permeability models reproduce well the measured pressure drop and concentration data

On s'intéresse dans ce travail à la modélisation d'un écoulement diphasique gaz-liquide co-courant descendant dans les réacteurs à lit fixe de particules, procédé largement utilisé dans le domaine industriel. En raison de la complexité de l'écoulement, induite par les nombreuses configurations multiphasiques pouvant coexister au sein du lit, les modèles développés directement à l'échelle du réacteur sont généralement issus d'approches semi-empiriques, en considérant l'écoulement homogène. Or, il a été observé que des hétérogénéités locales, géométrique et hydrodynamique, telle qu'une mal-distribution de la phase liquide, entrainaient une diminution du taux de réaction et conduisait les modèles existants à surestimer la productivité d'un réacteur. La nécessité de prendre en compte les phénomènes microscopiques dans un modèle macroscopique à l'échelle du réacteur rend l'utilisation d'approches multi-échelles indispensable. L'écoulement étant cependant d'une nature complexe, le changement d'échelle ne peut se faire de façon directe et nécessite donc la mise en place d'outils de modélisation adaptés à une échelle intermédiaire. Dans une première étape, la méthode de simulation numérique directe ``Volume-Of-Fluid'' (VOF) est validée dans le cas d'un film ruisselant dans un tube capillaire. Cette méthode est ensuite utilisée, à l'échelle microscopique, afin de proposer et de valider des relations de fermeture pour un modèle de type ``réseau de pores'' pouvant être utilisé à une échelle intermédiaire, celle du Volume Elémentaire Représentatif. Ce changement d'échelle est tout d'abord effectué dans le cas d'un lit fixe en deux dimensions, c'est-à-dire un empilement de cylindres entre deux plaques. Cette configuration permet la mise en place d'un dispositif expérimental qui, couplé à des simulations VOF 2D à plus grande échelle, valide l'approche de type "réseau de pores" adoptée. Le modèle réseau est ensuite étendu au cas d'un lit fixe réel, c'est-à-dire en trois dimensions, dont la géométrie est obtenue par micro-tomographie. Les lois de comportement locales sont redéfinies à l'aide de simulations numériques directes à l'échelle microscopique. Les résultats provenant de simulations de type « réseaux de pore » sont ensuite confrontés, dans le cas d'une répartition homogène des phases, aux modèles 1D habituellement utilisés pour les écoulements diphasiques en lit fixe. Enfin, une campagne expérimentale est menée afin d'observer, par imagerie scanner, l'étalement d'un jet de liquide sur un empilement de grains. Une comparaison qualitative est ensuite effectuée entre les observations expérimentales et les simulations numériques réseaux dans le cas spécifique de l'étalement d'un jet de liquide / We study in this work the modelling of two-phase cocurrent downflows in fixed bed reactors, a process widely used in industry. Due to the flow complexity, i.e., the presence of different interface configurations and, therefore, different phase interactions, most models have been developed using empirical approaches, with the assumption of a homogeneous flow in the reactor. However, several studies showed that local heterogeneities, geometric and hydrodynamic, such as the liquid distribution, could have a great influence on the flow at the reactor-scale and, therefore, on the reactor performance. Consider the microscopic phenomena in a macroscopic model require the use of multi-scale approaches. However, due to the flow complexity, the upscaling cannot be done directly and requires the development of modelling tools suitable for an intermediate scale. In a first step, the direct numerical method \ Volume-Of-Fluid" (VOF) is validated in the case of a two-phase flow in a capillary tube with the presence of a thin film. Then, this method is used, at a microscopic level to propose and validate closure laws for a pore-network model which will be used to simulate the flow at the intermediate scale. This upscaling approach is first tested in a two-dimensional case,i.e., an array of cylinders between two walls. This configuration allows the set up of an experimental approach, coupled with 2D VOF simulations at the intermediate scale, in order to validate the pore-network approach. The pore-network approach is then extended to a real fixed bed, i.e. in three dimensions, whose geometry is obtained by micro-tomography. Local laws of the pore-network model are redefined using direct numerical simulations at a microscopic scale. Pore-network simulations are then compared, for a homogenous phase distribution, with 1D models typically used for two-phase flow in fixed beds. Finally, an experimental campaign was set up to observe, by imaging scanner, the spreading of a liquid jet on a fixed bed pilot. A qualitative comparison is then performed between experimental observations and pore-network simulations in the specific case of the spreading of a liquid jet

Le but de cette étude expérimentale est de caractériser la dynamique de particules solides, à flottabilité nulle, incluse dans un écoulement turbulent isotrope en décroissance libre. Les particules utilisées sont de forme sphérique et ont un diamètre de 4 à 5 fois plus grand que l'échelle spatial de Kolmogorov de l'écoulement. De part leur taille, les particules ont également un nombre de Stokes proche de l'unité. On s'attend alors à ce que ces particules aient une dynamique différente de celle du fluide environnant. Dans cette étude, ont se propose de quantifier les différences de vitesses entre les deux phases à l'aide d'une technique de vélocimétrie simultanée / In this experimental study, the focus is made on the characterization of the dynamics of solid neutrally buoyant particles embedded in a freely decaying, nearly isotropic turbulence, with a weak mean flow. The particles are spherical with diameters several times larger than the Kolmogorov scale. The study of this flow configuration is still challenging both theoretically and numerically. Due to large particle sizes, the local flow around particles can not be considered as uniform and due to fluid-particle density ratio of around unity, the history and Basset forces cannot be neglected in comparison with the viscous drag force. Particle equation of motion is then fully non-linear, in contrast to the equation for heavy particles with diameters smaller then the Kolmogorov scale, for which only the Stokes drag is considered. In several experimental and numerical studies, the effect of particle size on velocity and acceleration statistics has been investigated (Homann and Bec 2010 ; Qureshi et al. 2008 ; Ouellette et al. 2008 ; Xu and Bodenschatz 2008). In the case of isotropic turbulence, Homann and Bec (2010) show that while the PDF of the particle velocity normalized by the square root of its variance does not vary with particle size, the variance itself is size dependent. A scaling relation for particle velocity variance has been proposed by using the Faxen correction (Gatignol 1983) which takes into account the non uniformity of the fluid flow at the scale of the particle. The aim of our research is to further study the dependence of particle dynamics on particle size. To that purpose, a turbulence generator has been set-up and the resulting turbulence is characterized. Then the flow was seeded with millimeter sized, neutrally-buoyant particles and the velocity of the two phases have been measured simultaneously. Simultaneous measurements of particle and surrounding fluid velocities show that although the global velocity statistics of the two phases have comparable values, the particles may have different local velocity from the velocity of the neighboring fluid