Simulation aux grandes échelles d'écoulements diphasiques turbulents à phase liquide dispersée

Les écoulements diphasiques turbulents sont présents dans de nombreux systèmes industriels (moteur à piston, turbines à gaz, moteurs fusée...). La compréhension fine de telles configurations s'avèrent de nos jours nécessaire pour limiter notamment les émissions de polluants et de gaz à effet de serre, et la consommation des énergies fossiles. Nous nous intéressons ici à la simulation aux grandes échelles des écoulements diphasiques turbulents, permettant de capturer une large partie du spectre de la turbulence, et ainsi être capable de prédire des phénomènes instables ou transitoires. La phase dispersée est ici modélisée par une approche eulérienne, en raison de ses avantages dans le contexte du calcul haute performance. Le travail de cette thèse a consisté à étendre le formalisme eulérien existant dans le code AVBP à la simulation de sprays polydisperses dans des écoulements turbulents. Pour cela, le Formalisme Eulérien Mésoscopique (FEM) a été couplé à une approche Multi-fluide. Cette nouvelle approche, intitulée Formalisme Eulérien Mésoscopique Multi-fluide (FEMM), a été évaluée sur des cas simples canoniques, permettant de bien caractériser le comportement autant en terme de dynamique turbulente que d'effets polydisperses. Les stratégies numériques disponibles dans le code de calcul AVBP sont aussi analysées, afin d'en cerner les limites pour la simulation eulérienne d'une phase liquide. Ce nouveau formalisme est finalement appliqué à la configuration aéronautique MERCATO, pour laquelle on dispose de résultats numériques obtenus avec d'autres approches (FEM et approche lagrangienne), et de résultats expérimentaux. Un accord satisfaisant avec l'expérience est montré pour toutes les approches, même si le FEM, monodisperse, obtient de moins bon résultats en terme de fluctuations. D'autres résultats expérimentaux s'avèrent nécessaires pour évaluer les approches et déterminer quelle est la plus prédictive pour cette configuration, notamment concernant la fraction massique de kerosene, autant en phase liquide qu'en phase gazeuse. ABSTRACT : Turbulent two-phase flows are encountered in several industrial devices (piston engine, gas turbine, rocket engine...). A fine understanding of such configurations is mandatory to face problems of pollutant emissions, greenhouse gas, and fossil fuel rarefaction. The Large Eddy Simulation seems to be a good candidate. This kind of simulation captures a wide part of turbulence spectrum, and thus allows to predict instabilities and transient phenomena. The dispersed phase is simulated using an Eulerian approach, which seems to be more suitable than lagrangian methods for High Performance Computing. The present work consists in the extension to polydisperse flows of the existing eulerian formalism in the AVBP code. The Mesoscopic Eulerian Formalism (MEF) is coupled with the Multifluid approach. This new formalism, called Multifluid Mesoscopic Eulerian Formalism, is evaluated on simple test cases, showing the ability of such approach to capture turbulent and polydisperse effects. Numerical strategies available in AVBP are also evaluated, in order to emphasize on their limiting aspects for the eulerian simulation of a dispersed phase. The new formalism is finally applied to the simulation of the aeronautical configuration called MERCATO. Several experimental results are available, as well as numerical results using FEM and lagrangian approach. Results show a good agreement between experiments and numerical results, even if FEM results are worse concerning the fluctuations. New experimental results are necessary to determine which is the best approach, especially in terms of liquid and gas kerosene mass fraction

On the Direct Numerical Simulation of moderate-Stokes-number turbulent particulate flows using Algebraic-Closure-Based and Kinetic-Based Moment Methods

In turbulent particulate flows, the occurrence of particle trajectory crossings (PTC) is the main constraint on classical monokinetic Eulerian methods. To handle such PTC, accounting for high-order moments of the particle velocity distribution is mandatory. In the simplest case, second-order moments are needed. To retrieve these moments, two solutions are proposed in the literature: the Algebraic-Closure-Based Moment Method (ACBMM) and the Kinetic-Based Moment Method (KBMM). The ACBMM provides constitutive relations for the random-uncorrelated-motion (RUM) particle kinetic stress tensor as algebraic closures based on physical arguments (Simonin et al. 2002; Kaufmann et al. 2008; Masi 2010; Masi & Simonin 2012). These closures rely on the internal energy, namely the RUM particle kinetic energy, which is obtained using an additional transport equation. Alternatively, it is possible to directly solve for the second-order moment by providing a closure for the third-order correlation. The KBMM proposes a kinetic description, that is, the number density function (NDF) is reconstructed based on the resolved moments and on a presumed shape. In the present work, an isotropic Gaussian and the anisotropic Gaussian closure of Vié et al. (2012) are used. The goal of the present study is to provide a first comparison between ACBMM and KBMM, using the same robust numerical methods, in order to highlight differences and common points. The test case is a 2D turbulent flow with a mean shear

Particle-laden flows forced by the disperse phase: Comparison between Lagrangian and Eulerian simulations

The goal of the present work is to assess the ability of Eulerian moment methods to reproduce the physics of two-way coupled particle-laden turbulent flow systems. Previous investigations have been focused on effects such as preferential concentration, and turbulence modulation, but in regimes in which turbulence is sustained by an imposed external forcing. We show that in such regimes, Eulerian methods need resolutions finer than nominal Kolmogorov scale in order to capture statistics of particle segregation, but gas and disperse phase velocity variances can be captured with resolutions comparable to the Kolmogorov length. The work is then extended to address the question whether Eulerian methods are suitable in scenarios in which the continuum field of interest (temperature or momentum) is itself primarily driven by particles. To this end we have extended our analysis to the problem of turbulence driven by heated particles (Zamansky et al. PoF 2014) and have assessed capabilities of Eulerian methods in capturing particle segregation, as well as statistics of the temperature and velocity fields. Separate investigations are developed for cases with and without buoyancy driven turbulence. For each case corresponding Lagrangian calculations are developed and convergence of statistics with respect to the number of particles is established. Then the statistically- converged Lagrangian and Eulerian results are compared. Results show that accurate capture of segregation by the Eulerian methods always requires resolutions much higher than the nominal Kolmogorov scale. In scenarios for which a continuum phase is forced by particles, results from Eulerian methods show some sensitivity of predicted continuum statistics to the mesh resolution. This sensitivity was found to be largest for the case of a temperature field forced by hot particles, but without presence of buoyancy. In this case a Eulerian method with nominal Kolmogorov resolution was found to be insufficient for capture of temperature statistics. When additional coupling between particles and continuum phase is introduced by including the buoyancy effects, this sensitivity is suppressed in the temperature field, but some sensitivity to the Eulerian mesh resolution were detected in the momentum fields

Reexamining the framework for intermittency in Lagrangian stochastic models for turbulent flows: a way to an original and versatile numerical approach

The characterization of intermittency in turbulence has its roots in the refined similarity hypotheses of Kolmogorov, and if no proper definition is to be found in the literature, statistical properties of intermittency were studied and models were developed in an attempt to reproduce it. The first contribution of this work is to propose a requirement list to be satisfied by models designed within the Lagrangian framework. Multifractal stochastic processes are a natural choice to retrieve multifractal properties of the dissipation. Among them, we investigate the Gaussian multiplicative chaos formalism, which requires the construction of a log-correlated stochastic process Xt. The fractional Gaussian noise of Hurst parameter H=0 is of great interest because it leads to a log correlation for the logarithm of the process. Inspired by the approximation of fractional Brownian motion by an infinite weighted sum of correlated Ornstein-Uhlenbeck processes, our second contribution is to propose a stochastic model: Xt=∫∞0Yxtk(x)dx, where Yxt is an Ornstein-Uhlenbeck process with speed of mean reversion x and k is a kernel. A regularization of k(x) is required to ensure stationarity, finite variance, and logarithmic autocorrelation. A variety of regularizations are conceivable, and we show that they lead to the aforementioned multifractal models. To simulate the process, we eventually design a new approach relying on a limited number of modes for approximating the integral through a quadrature XNt=∑Ni=1ωiYxit, using a conventional quadrature method. This method can retrieve the expected behavior with only one mode per decade, making this strategy versatile and computationally attractive for simulating such processes, while remaining within the proposed framework for a proper description of intermittency

Large eddy simulation of turbulent gas-dispersed liquid two-phase flows

Les écoulements diphasiques turbulents sont présents dans de nombreux systèmes industriels (moteur à piston, turbines à gaz, moteurs fusée...). La compréhension fine de telles configurations s'avèrent de nos jours nécessaire pour limiter notamment les émissions de polluants et de gaz à effet de serre, et la consommation des énergies fossiles. Nous nous intéressons ici à la simulation aux grandes échelles des écoulements diphasiques turbulents, permettant de capturer une large partie du spectre de la turbulence, et ainsi être capable de prédire des phénomènes instables ou transitoires. La phase dispersée est ici modélisée par une approche eulérienne, en raison de ses avantages dans le contexte du calcul haute performance. Le travail de cette thèse a consisté à étendre le formalisme eulérien existant dans le code AVBP à la simulation de sprays polydisperses dans des écoulements turbulents. Pour cela, le Formalisme Eulérien Mésoscopique (FEM) a été couplé à une approche Multi-fluide. Cette nouvelle approche, intitulée Formalisme Eulérien Mésoscopique Multi-fluide (FEMM), a été évaluée sur des cas simples canoniques, permettant de bien caractériser le comportement autant en terme de dynamique turbulente que d'effets polydisperses. Les stratégies numériques disponibles dans le code de calcul AVBP sont aussi analysées, afin d'en cerner les limites pour la simulation eulérienne d'une phase liquide. Ce nouveau formalisme est finalement appliqué à la configuration aéronautique MERCATO, pour laquelle on dispose de résultats numériques obtenus avec d'autres approches (FEM et approche lagrangienne), et de résultats expérimentaux. Un accord satisfaisant avec l'expérience est montré pour toutes les approches, même si le FEM, monodisperse, obtient de moins bon résultats en terme de fluctuations. D'autres résultats expérimentaux s'avèrent nécessaires pour évaluer les approches et déterminer quelle est la plus prédictive pour cette configuration, notamment concernant la fraction massique de kerosene, autant en phase liquide qu'en phase gazeuse.Turbulent two-phase flows are encountered in several industrial devices (piston engine, gas turbine, rocket engine...). A fine understanding of such configurations is mandatory to face problems of pollutant emissions, greenhouse gas, and fossil fuel rarefaction. The Large Eddy Simulation seems to be a good candidate. This kind of simulation captures a wide part of turbulence spectrum, and thus allows to predict instabilities and transient phenomena. The dispersed phase is simulated using an Eulerian approach, which seems to be more suitable than lagrangian methods for High Performance Computing. The present work consists in the extension to polydisperse flows of the existing eulerian formalism in the AVBP code. The Mesoscopic Eulerian Formalism (MEF) is coupled with the Multifluid approach. This new formalism, called Multifluid Mesoscopic Eulerian Formalism, is evaluated on simple test cases, showing the ability of such approach to capture turbulent and polydisperse effects. Numerical strategies available in AVBP are also evaluated, in order to emphasize on their limiting aspects for the eulerian simulation of a dispersed phase. The new formalism is finally applied to the simulation of the aeronautical configuration called MERCATO. Several experimental results are available, as well as numerical results using FEM and lagrangian approach. Results show a good agreement between experiments and numerical results, even if FEM results are worse concerning the fluctuations. New experimental results are necessary to determine which is the best approach, especially in terms of liquid and gas kerosene mass fraction

How Connected is Too Connected? Impact of Network Topology on Systemic Risk and Collapse of Complex Economic Systems

Abstract Economic interdependencies have become increasingly present in globalized production, financial and trade systems. While establishing interdependencies among economic agents is crucial for the production of complex products, they may also increase systemic risk due to failure propagation. It is crucial to identify how network connectivity impacts both the emergent production and risk of collapse of economic systems. In this paper we propose a model to study the effects of network structure on the behavior of economic systems by varying the density and centralization of connections among agents. The complexity of production increases with connectivity given the combinatorial explosion of parts and products. Emergent systemic risks arise when interconnections increase vulnerabilities. Our results suggest a universal description of economic collapse given in the emergence of tipping points and phase transitions in the relationship between network structure and risk of individual failure. This relationship seems to follow a sigmoidal form in the case of increasingly denser or centralized networks. The model sheds new light on the relevance of policies for the growth of economic complexity, and highlights the trade-off between increasing the potential production of the system and its robustness to collapse. We discuss the policy implications of intervening in the organization of interconnections and system features, and stress how different network structures and node characteristics suggest different directions in order to promote complex and robust economic systems

Large eddy simulations of a liquid fuel swirl burner: flame characterization for pilot and multipoint injection strategies

International audienceA laboratory-scale swirling burner fueled with dodecane is studied numerically using large eddy simulations. The burner is composed of two stages, allowing the use of two different injector types (namely pilot and multipoint) in order to study the influence of droplet size and initial position with fixed geometry and delivered power. For a chosen lean operating point, the two liquid injection types are tested, highlighting a dramatic influence on the flame stabilization process. When fuel is injected through the multipoint stage, evaporation and mixing are enhanced and a partially premixed mixture enters the combustion chamber. The flame then takes an ‘M’ shape, mainly controlled by the large inner and outer recirculation zones associated with this highly swirling flow and in which trapped burnt gases guarantee permanent ignition of the fresh mixture entering the chamber. The situation is much more complex when fuel is solely injected through the pilot nozzle. Due to the large amount of liquid fuel present in the pilot zone, premixing is not achieved and the flame must stabilize itself mainly in a hybrid combustion regime. This is only possible thanks to a very complex situation in this region, where hot evaporated fuel is trapped in front of the nozzle and oxygen is mainly coming from the large central recirculation zone. In that case, the flame takes a ‘tulip’ shape, with a stabilization point inside the injection device. Both flame shapes are compared using scatterplots and flame dynamics are analyzed