Current climate models often predict fractional cloud cover on the basis of a diagnostic probability density function (PDF) describing the subgrid-scale variability of the total water specific humidity, qt, favouring schemes with limited complexity. Standard shapes are uniform or triangular PDFs the width of which is assumed to scale with the gridbox mean qt or the grid-box mean saturation specific humidity, qs. In this study, the qt variability is analysed from large-eddy simulations for two stratocumulus, two shallow cumulus, and one deep convective cases. We find that in most cases, triangles are a better approximation to the simulated PDFs than uniform distributions. In two of the 24 slices examined, the actual distributions were so strongly skewed that the simple symmetric shapes could not capture the PDF at all. The distribution width for either shape scales acceptably well with both the mean value of qt and qs, the former being a slightly better choice. The qt variance is underestimated by the fitted PDFs, but overestimated by the existing parameterisations. While the cloud fraction is in general relatively well diagnosed from fitted or parameterised uniform or triangular PDFs, it fails to capture cases with small partial cloudiness, and in 10 – 30% of the cases misdiagnoses clouds in clear skies or vice-versa. The results suggest choosing a parameterisation with a triangular shape, where the distribution width would scale with the grid-box mean qt using a scaling factor of 0.076. This, however, is subject to the caveat that the reference simulations examined here were partly for rather small domains and driven by idealised boundary conditions.

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.

Hot wire anemometers have been used in several wind velocity sensors deployed in Mars. They are based in keeping the temperature of a surface at a constant value, above the ambient. This is done by means of a heater controlled with an electronic system. The cooling rate of each point at the sensor surface can be used to calculate the wind velocity and direction. However, due to turbulent fluctuations, the cooling rate is not constant even in the case of constant velocity. Moreover, RANS simulations cannot estimate such fluctuations as they only provide an estimation of the averaged flow field. The goal of this work has been to estimate such fluctuations and the e↵ect they might have on the sensor readings. To do so, the turbulent cooling rate (Nusselt number) of a sensor with a generic shape, under the typical conditions to be find in Mars, has been simulated using high performance LES (Large Eddy Simulation) simulations and compared with RANS and URANS simulations.

Les réseaux de tuyauterie industriels sont le siège de niveaux importants de vibrations induites par l’écoulement qui peuvent mener à la rupture par fatigue des installations. La présente étude se concentre sur l’analyse et modélisation simplifiée de la source vibratoire associée au passage d’un écoulement liquide turbulent par un coude à 90°. Une approche combinant expériences et simulation a été conduite. Une boucle de circulation d’eau munie d’un coude transparent a été conçue afin de permettre des me-sures de vitesse à l’intérieur du coude. Pour cette finalité, les techniques de Vélocimétrie par Imagerie de Particules (PIV) plane et stéréoscopique ont été utilisées. La pression pariétale et les vibrations du banc d’essais ont été mesurées simultanément. Plusieurs configurations d’écoulement ont été testées afin d’obtenir une riche base de données couplées reliant l’écoulement fluide à l’excitation dynamique des parois et, finalement, à la réponse vibratoire de la structure. En parallèle, l’écoulement instationnaire d’eau dans le coude a été simulé au moyen d’une approche du type Simulation des Grandes Echelles (LES). La simulation fluide a permis d’étudier en détails la topologie de l’écoulement turbulent au passage du coude ainsi que le champ instationnaire de pression fluctuante induit sur la paroi. Finale-ment, un ensemble d’outils statistiques a été appliqué aux données expérimentales et numériques afin de proposer un modèle simplifié des transferts qui relient l’écoulement turbulent à la sollicitation dynamique de la structure contenant le coude. / The flow of dense fluids within thin-walled piping systems may lead to significant levels of Flow-Induced Vibration, mainly in the vicinities of singularities such as obstacles inserted into the flow, sudden changes of cross-sectional area or flow direction. This study focuses on the analysis and reduced-order modelling of the vibrational source associated with the turbulent flow of liquids through a 90° elbow.A mixed experimental-computational approach is undertaken. A closed water loop containing a transparent elbow was designed in order to allow for fluid velocity measurements inside the singularity. To this purpose, planar and stereoscopic Particle Image Velocimetry (PIV) were employed. Wall pressure fluctuations and structural vibrations were measured simultaneously. Several flow configurations were tested in order to obtain a large coupled database linking the flow to the dynamic excitation, and then to the vibration response of the structure.In parallel, the unsteady water flow through the elbow was computed using Large-Eddy Simulation (LES). The fluid simulation allowed for a detailed study of the turbulent flow through the singularity and the unsteady pressure field induced on the piping walls. Finally, a set of statistical tools was applied to both experimental and computational data in order to propose a reduced-order model of the transfer function that links the tur-bulent flow to the dynamic excitation of the elbowed piping structure.

Interactions between flames and thermal phenomena are the guiding thread of this work. Flamesproduce heat indeed, but can also be affected by it. Large Eddy Simulations (LES) are used hereto investigate these interactions, with a focus on two main topics: wall heat transfer andcombustion control. In a first part, wall heat transfer in a rocket engine sub-scale CH4/O2 burner isstudied. In the context of launchers re-usability and cost reduction, which are major challenges,new propellant combinations are considered and wall heat fluxes have to be precisely predicted.The aim of this work is to evaluate LES needs and performances to simulate this kind ofconfiguration and provide a computational methodology permitting to simulate variousconfigurations. Numerical results are compared to experimental data provided by the TechnischeUniversität München (Germany). In a second part, combustion control by means of NanosecondRepetitively Pulsed (NRP) plasma discharges is studied. Modern gas turbine systems use indeedlean combustion with the aim of reducing fuel consumption and pollutant emissions. Lean flamesare however known to be prone to instabilities and combustion control can play a major role in thisdomain. A phenomenological model which considers the plasma discharges as a heat source isdeveloped and applied to a swirl-stabilized CH4/Air premixed lean burner. LES are performed inorder to evaluate the effects of the NRP discharges on the flame. Numerical results are comparedwith experimental observations made at the King Abdulla University of Science and Technology(Saudi Arabia).

Pour des raisons de sécurité, les moteurs aéronautiques doivent pouvoir redémarrer en vol sur toute leur plage d'opération. Mais les contraintes sur les émissions polluantes nécessitent le développement de nouvelles chambres de combustion dont la conception peut détériorer les capacités d'allumage du moteur. Afin d'améliorer la compréhension du processus d'allumage et d'aider à l'optimisation de la conception, les recherches actuelles combinent les études expérimentales de plus en plus complexes et les simulation numériques hautes fidélités. Dans ce travail, l'étude numérique du processus d'allumage des chambres de combustion aéronautiques, de l'étincelle à la propagation azimutale de la flamme, est conduite avec plusieurs objectifs: améliorer la robustesse et la confiance de l'outil LES pour l'étude de l'allumage, étudier les mécanismes qui affectent l'allumage dans des conditions représentatives des conditions réelles et enfin améliorer les méthodes bas-ordre pour la prédiction des performances d'allumage. Dans une première partie, la SGE d'un monobruleur installé au CORIA permet de mettre en évidence les bons résultats de la LES et de construire une base de données pour l'analyses des mécanismes d'extinction. Ces données sont aussi utilisées pour développer une méthodologie permettant de prédire les performances d'allumage à bas coût en utilisant les résultats d'une SGE non-réactive. Dans une seconde partie, la propagation inter-secteur est investiguée par l'étude de deux cas expérimentaux et la SGE est capable de reproduire les modes de propagation mais aussi les temps d'allumage avec précision. Sur la bases de ces bons résultats, une analyse plus fine de la simulation permet d'identifier les mécanismes qui contrôlent la propagation de la flamme. / For safety reasons, in-flight relight of the engine must be guaranteed over a wide range of operating conditions but the increasing stringency of pollutant emission constraints requires the development of new aero-engine combustor whose design might be detrimental to the ignition capability. To improve the knowledge of the ignition process in aeronautical gas turbines and better combine conflicting technological solutions, current research relies on both complex experimental investigation and high fidelity numerical simulations. In this work, numerical study of the ignition process in gas turbines from the energy deposit to the light-around is performed with several objectives: increase the level of confidence of Large Eddy Simulations tool for the analysis of the ignition process, investigate the mechanisms controlling ignition in conditions representative of realistic aeronautical gas turbine flows and improve the low-order methodologies for the prediction of ignition performance. In a first part, LES of the single burner installed at CORIA (France) is carried out and allows to highlight the LES accuracy and to build a database on which the main mechanisms controlling the ignition success are identified. Based on these results, a methodology is developed to predict the ignition performance at a low computational cost using the non-reacting flow statistics only. In a second part, the light-around process is studied on two experimental set-ups and the very good agreement of the LES results with experiments is the starting point from an analysis of the mechanisms driving the flame propagation process.

Investigating flow and turbulence structure around a barrier mounted on the ground or placed in its vicinity is a fundamental problem in wind engineering because of many practical applications related to protection against adverse effects induced by major wind storms (e.g., hurricanes) and snow events (e.g., snow fences used to reduce adverse effects of snow drifting on the roads). In this work the focus is on the case when the obstacle/barrier is porous and the shape of the obstacle is close to a high-aspect-ratio rectangular cylinder situated in the vicinity of the ground. The study employs a range of numerical and experimental techniques to achieve this goal that include 3D LES and 2D RANS numerical simulations, and RTK survey and 3D photogrammetry techniques to measure ground elevations and snow deposits in the field. In the first part of the study, high-resolution large eddy simulations are used to understand the fundamental flow physics of flow past 2D solid and porous vertical plates with a special focus on describing the unsteady wind loads on the obstacle, vortical structure of the turbulent wake, spectral content of the wake, the separated shear layers and of the characteristics of the large-scale vortex shedding behind the plate, if present. Results show that LES can accurately predict mean flow and turbulence statistics around solid/porous cylinders. Then, a detailed parametric study of flow past vertical solid and porous plates situated in the vicinity of a horizontal bed is performed for the purpose of understanding changes in the mean flow structure, turbulence statistics and dynamics of large scale coherent structures as a function of the main nondimensional geometrical parameters (bottom gap for solid and porous plates, and porosity and average hole size of porous plates) and flow variables (e.g., bed roughness) that affect the wake flow. In particular, the LES flow fields allowed clarifying how the interactions between the bottom and the top separated shear layers change with increasing bottom gap and what is the effect of the bleeding flow on the interactions between the separated shear layers that determine the coherence of the large-scale eddies at large distances from the wake. In the second part of the thesis, a novel methodology based on field monitoring of the snow deposits and RANS numerical simulations is proposed to improve the design of snow fences and in particular the design of lightweight plastic snow fences that are commonly used to protect roads in the US Midwest against the snow drifting. The goal of the design optimization procedure is to propose a snow fence design that can retain a considerable amount of snow within a shorter downwind distance compared to fences of standard design. A major contribution of the present thesis was the development of a novel non-intrusive image-based technique that can be used to quantitatively estimate the temporal evolution of the volume of snow trapped by a fence over long periods of time. This technique is based on 3-D close range photogrammetry. Results showed that this technique can produce estimations of the snow deposits of comparable accuracy to that given by commonly used methods. This is the first application of this type of techniques to measurements of the snow deposits.

Nous cherchons à modéliser numériquement un jet flottant air-hélium dans une cavité avec deux ouvertures à partir de simulations aux grandes échelles (LES) et de simulations numériques directes (DNS). La configuration considérée est basée sur une étude expérimentale menée au CEA de Saclay reproduisant une fuite d'hydrogène en environnement confiné. La dimension de la cavité a été choisie pour permettre une transition laminaire-turbulent intervenant environ à la mi-hauteur de la cavité. Cette étude porte principalement sur trois points majeurs : l'influence des conditions aux limites sur le développement du jet et son interaction avec l'environnement extérieur, la validité du modèle numérique qui est analysée en comparant la distribution de vitesse obtenue numériquement aux mesures expérimentales (PIV) et, enfin, la compréhension de la distribution air-hélium et le phénomène de stratification qui s'établit à l'intérieur de la cavité. Nous observons dans un premier temps que des conditions limites de pression constante appliquées directement au ras des évents conduisent à une sous-estimation du débit volumique d'air entrant dans la cavité et donc à une surestimation de la masse de l'hélium à l'intérieur de la cavité, ce qui n'est pas acceptable dans un contexte d'évaluation du risque hydrogène. En revanche, la prise en compte, dans le domaine de calcul, d'une région extérieure à la cavité prédit correctement le flux d'air entrant. Les résultats numériques sont alors en bon accord avec les données PIV. Il a été montré que les prédictions de la DNS, par rapport à la LES, concordent mieux avec les mesures de vitesse par PIV. Le champ de concentration prédit numériquement présente une couche homogène en haut de la cavité, dont la concentration est en accord avec le modèle théorique de Linden et al. 1990. Cependant, sa position et son épaisseur ne correspondent pas au modèle. Ceci est principalement dû aux interactions directes entre le jet flottant et, d'une part, avec les limites solides de la cavité et d'autre part, avec l'environnement extérieur. L'analyse statistique concernant la production de la flottabilité de l'énergie cinétique turbulente (TKE) a permis d'identifier les limites du jet flottant. / We present numerical results from large eddy simulations (LES) and coarse direct numerical simulations (DNS) of an air-helium buoyant jet rising in a two vented cavity. The geometrical configuration mimics the helium release experimental set-up studied at CEA Saclay in the framework of security assessment of hydrogen-based systems with an indoor usage. The dimension of the enclosure was chosen to ensure a laminar-turbulent transition occurring at about the middle height of the cavity. This study focuses mainly on three key points : the influence of the boundary conditions on the jet development and its interaction with the exterior environment, the validity of the numerical model which is analyzed by comparing the numerical velocity distribution versus the measured particle image velocimetry (PIV) ones, and finally understanding the distribution of air-helium and the stratification phenomenon that takes place inside the cavity. We observe at first that applying constant pressure outlet boundary conditions directly at the vent surfaces underestimates the volumetric flow rate of air entering the enclosure and thus overestimate the helium mass inside the cavity. On the contrary, modelling an exterior region in the computational domain better predicts the air flow-rate entrance and numerical results matches better with the experimental PIV data. It has been figured out that the coarse DNS predictions match better with the velocities PIV measurements, compared to the LES. Numerical prediction of the helium field depicts a homogeneous layer formed at the top of the cavity, with a concentration in good agreement with the theoretical model of Linden et al. 1990. However, the position and the thickness of the layer do not correspond to the theory. This is mainly due to the direct interactions between the buoyant jet and both the solid boundaries of the cavity and the exterior environment. Statistical analysis regarding the buoyancy production of the turbulent kinetic energy (TKE) served to identify the limits of the buoyant jet.

Among the different types of flows encountered in practical applications, the physics of turbulent separated flows and turbulence induced secondary motion are not fully understood despite the large amount of previous experimental and numerical work. The objectives of this work are to study theoretically and computationally the conditions at the separation and reattachment point, the numerical simulation of turbulence induced secondary motion in non-circular ducts, and to provide a comprehensive test of different RANS models of these types of flow. In a theoretical study of flow separation, a Lagrangian approach was first used to derive an Eulerian criterion, which associates separation and reattachment points to a critical point in the eigenvalues of the Cauchy-Green tensor. A turbulent separated boundary layer under the influence of an adverse pressure gradient was simulated using DNS and LES techniques. A bootstrapping method was used to obtain high fidelity results at a relatively high Reynolds number with which the performance of some of the most commonly used eddy-viscosity turbulence models was evaluated. The DNS and LES results were used to assess the consistency of the different terms in the k−e , ζ −f , k −ω and Spalart-Allmaras models. Different wall-modelling techniques were employed for the calculation of separated boundary layers. The exact values of the modelled terms were calculated using the reference DNS and LES dataset. These results were used for both a priori and a posteriori tests. It was determined that the eddy-viscosity assumption works well, and that anisotropic effects are not significant in separated boundary layer. For the secondary flow calculation in non-circular ducts, direct numerical simulations of turbulent flow in square and skewed ducts were carried out to determine the effect of the duct (rhombus) included angle on both the mean and turbulence energy budgets. Two skewed ducts, with included angles of 30 and 60 degrees, were simulated. The capability of different turbulence models to predict the secondary velocity field was investigated. Results obtained from a non-linear stress-strain constitutive relation was found to be fairly accurate for the flows at the range of Reynolds number considered in this study. / Thesis (Ph.D, Mechanical and Materials Engineering) -- Queen's University, 2010-11-26 13:52:18.361

The increased fatigue loads and decreased power output of a wind turbine placed in the wake of another turbine is a well-known problem when building new wind-power farms and a subject of intensive research. These problems are caused by the velocity gradients and high turbulence levels present in the wake of a turbine. In order to better estimate the total power output and life time of a wind-power farm, knowledge about the development and stability of wind-turbine wakes is crucial. In the present thesis, the flow field around small-scale model wind turbines has been investigated experimentally in two wind tunnels. The flow velocity was measured with both hot-wire anemometry and particle image velocimetry. To monitor the turbine performance, the rotational frequency, the power output and the total drag force on the turbine were also measured. The power and thrust coefficients for different tip-speed ratios were calculated and compared to the blade element momentum method, with a reasonable agreement. The same method was also used to design and manufacture new turbine blades, which gave an estimate of the distribution of the lift and drag forces along the blades. The influence of the inlet conditions on the turbine and the wake properties was studied by subjecting the turbine to both uniform in flow and different types of boundary layer in flows. In order to study the stability and development of the tip vortices shed from the turbine blades, a new experimental setup for phase-locked measurements was constructed. The setup made it possible to introduce perturbations of different frequencies and amplitudes, located in the rear part of the nacelle. With a newly developed method, it was possible to characterize the vortices and follow their development downstream, using only the streamwise velocity component. Measurements were also performed on porous discs placed in different configurations. The results highlighted the importance of turbine spacings. Both the measurements on the turbine and the discs were also used to compare with large eddy simulations using the actuator disc method. The simulations managed to predict the mean velocity fairly well in both cases, while larger discrepancies were seen in the turbulence intensity. / <p>QC 20140424</p>