No description available.

Detached Eddy Simulations (DES) is a hybrid URANS-LES technique that was proposed to obtain computationally feasible solutions of high Reynolds number flows undergoing massive separation with reliable accuracy. Since its inception, DES has been applied to a wide variety of flow fields, but mostly limited to unbounded external aerodynamic flows. This is the first study to apply and validate DES to predict the internal flow and heat transfer in non-canonical flows of industrial relevance. The prediction capabilities of DES in capturing the effects of Coriolis forces, which are induced by rotation, and centrifugal buoyancy forces, which are induced by thermal gradients, are also authenticated. The accurate prediction of turbulent flows is sensitive to the level of turbulence predicted by the turbulence scheme. By treating the regions of interest in LES mode, DES allows the unsteadiness in these regions to develop and hence predicts the turbulence levels accurately. Additionally, this permits DES to capture the effects of system rotation and buoyancy. Computations on a rotating system (a sudden expansion duct) and a system subjected to thermal gradients (cavity with a heated wall) validate the prediction capability of DES. The application of DES is further extended to a non-canonical, internal flow which is of relevance in internal cooling of gas turbine blades. Computations of the fully developed flow and heat transfer shows that DES surpasses several shortcomings of the RANS model on which it is based. DES accurately predicts the primary and secondary flow features, the turbulence characteristics and the heat transfer in stationary ducts and in rotating ducts, where the effects of Coriolis forces and centrifugal buoyancy forces are dominant. DES computations are carried out at a computational cost that is almost an order of magnitude less than the LES with little compromise on the accuracy. However, the capabilities of DES in predicting the transition to turbulence are inadequate, as highlighted by the flow features and the heat transfer in the developing region of the duct. But once the flow becomes fully turbulent, DES predicts the flow physics and shows good quantitative agreement with the experiments and LES. / Ph. D.

The thesis deals with the flow in pipe bends and radial turbines geometries that are commonly found in an Internal Combustion Engine (ICE). The development phase of internal combustion engines relies more and more on simulations as an important complement to experiments. This is partly because of the reduction in development cost and the shortening of the development time. This is one of the reasons for the need of more accurate and predictive simulations. By using more complex computational methods the accuracy and predictive capabilities are increased. The disadvantage of using more sophisticated tools is that the computational time is increasing, making such tools less attractive for standard design purposes. Hence, one of the goals of the work has been to contribute to assess and improve the predictive capability of the simpler methods used by the industry. By comparing results from experiments, Reynolds Averaged Navier-Stokes (RANS) computations, and Large Eddy Simulations (LES) the accuracy of the different computational methods can be established. The advantages of using LES over RANS for the flows under consideration stems from the unsteadiness of the flow in the engine manifold. When such unsteadiness overlaps the natural turbulence the model lacks a rational foundation. The thesis considers the effect of the cyclic flow on the chosen numerical models. The LES calculations have proven to be able to predict the mean field and the fluctuations very well when compared to the experimental data. Also the effects of pulsatile exhaust flow on the performance of the turbine of a turbocharging system is assessed. Both steady and pulsating inlet conditions are considered for the turbine case, where the latter is a more realistic representation of the real flow situation inside the exhaust manifold and turbine. The results have been analysed using different methods: single point Fast Fourier Transforms (FFT), probe line means and statistics, area and volume based Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD). / Denna avhandling behandlar flödet i rörkrökar och radiella turbiner som vanligtvis återfinns i en förbränningsmotor. Utvecklingsfasen av förbränningsmotorer bygger mer och mer på att simuleringar är ett viktigt komplement till experiment. Detta beror delvis på minskade utvecklingskostnader men även på kortare utevklningstider. Detta är en av anledningarna till att man behöver mer exakta och prediktiva simuleringsmetoder. Genom att använda mer komplexa beräkningsmetoder så kan både nogrannheten och prediktiviteten öka. Nackdelen med att använda mer sofistikerade metoder är att beräkningstiden ökar, vilket medför att sådana verktyg är mindre attraktiva för standardiserade design ändamål. Härav, ett av målen med projektet har varit att bidra med att bedöma och förbättra de enklare metodernas prediktionsförmåga som används utav industrin. Genom att jämföra resultat från experiment, Reynolds Averaged Navier-Stokes (RANS) och Large Eddy Simulations (LES) så kan nogrannheten hos de olika simuleringsmetoderna fastställas. Fördelarna med att använda LES istället för RANS när det gäller de undersökta flödena kommer ifrån det instationära flödet i grenröret. När denna instationäritet överlappar den naturligt förekommande turbulensen så saknar modellen en rationell grund. Denna avhandling behandlar effekten av de cykliska flöderna på de valda numeriska modellerna. LES beräkningarna har bevisats kunna förutsäga medelfältet och fluktuationerna väldigt väl när man jämför med experimentell data. Effekterna som den pulserande avgasströmning har på turboladdarens turbin prestanda har också kunnat fastställas. Både konstant och pulserande inlopps randvillkor har används för turbinfallet, där det senare är ett mer realistiskt representation av den riktiga strömningsbilden innuti avgasgrenröret och turbinen. Resultaten har analyserats på flera olika sätt: snabba Fourier transformer (FFT) i enskilda punkter, medelvärden och statistik på problinjer, area och volumsbaserade metoder så som Proper Orthogonal Decomposition (POD) samt Dynamic Mode Decomposition (DMD). / <p>QC 20140919</p>

<p>Ground vehicles are subjected to crosswind from various origins such as weather, topography of the ambient environment (land, forest, tunnels, high bushes...) or surrounding traffic. The trend of lowering the weight of vehicles imposes a stronger need for understanding the coupling between crosswind stability, the vehicle external shape and the dynamic properties. Means for reducing fuel consumption of ground vehicles can also conflict with the handling and dynamic characteristics of the vehicle. Streamlined design of vehicle shapes to lower the drag can be a good example of this dilemma. If care is not taken, the streamlined shape can lead to an increase in yaw moment under crosswind conditions which results in a poor handling.</p><p>The development of numerical methods provides efficient tools to investigate these complex phenomena that are difficult to reproduce experimentally. Time accurate and scale resolving methods, like Detached-Eddy Simulations (DES), are particularly of interest, since they allow a better description of unsteady flows than standard Reynolds Average Navier-Stokes (RANS) models. Moreover, due to the constant increase in computational resources, this type of simulations complies more and more with industrial interests and design cycles.</p><p>In this thesis, the possibilities offered by DES to simulate unsteady crosswind aerodynamics of simple vehicle models in an industrial framework are explored. A large part of the work is devoted to the grid design, which is especially crucial for truthful results from DES. Additional concerns in simulations of unsteady crosswind aerodynamics are highlighted, especially for the resolution of the wind-gust boundary layer profiles. Finally, the transient behaviour of the aerodynamic loads and the flow structures are analyzed for several types of vehicles. The results simulated with DES are promising and the overall agreement with the experimental data available is good, which illustrates a certain reliability in the simulations. In addition, the simulations show that the force coefficients exhibit highly transient behaviour under gusty conditions.</p> / ECO2 Crosswind Stability and Unsteady Aerodynamics for Ground Vehicles

L'objectif de cette thèse est de développer un code de simulation complet pour le remplissage en fonderie de pièces de grandes dimensions (jusqu'à plusieurs mètres). Ce type de procédé fait entrer en jeu de nombreux phénomènes physiques couplés, nécessitant des méthodes numériques adaptées. La faible viscosité du métal liquide (de l'ordre de 10−6 m2/s) requiert l'emploi d'un modèle de turbulence basé sur un solveur Navier Stokes stabilisé et une méthode de suivi/capture d'interface. Nous avons pour cela choisi un approche stabilisée de type Variational Multi Scales (VMS), qui s'est révélée efficace pour simuler des nombres de Reynolds modérés, alliée à une méthode level-set permettant de déterminer de manière précise et à tout moment la position de l'interface liquide/air. La turbulence est quant à elle prise en compte grâce à un modèle dynamique de type Large Eddy Simulations (L.E.S.), ne faisant pas apparaître de paramètre empirique. Chacune de ces méthodes numériques a été confrontée à des résultats expérimentaux, numériques ou analytiques. Nous avons également conçu notre propre maquette expérimentale de remplissage d'eau, afin de valider le couplage des solveurs pour un cas représentatif. Une autre caractéristique de ces procédés à durée relativement longue (jusqu'à plusieurs dizaines de minutes) est l'importance des transferts thermiques, pouvant mener à la solidification du métal en cours de remplissage. Il convient donc de développer une méthode de résolution stabilisée de la thermique avec convection dominante. Cette méthode doit prendre en compte les variables turbulentes introduites précédemment. Enfin, nous proposons une méthode innovante pour simuler le changement de phase, basée sur une approche germination/ croissance avec fonction level-set. L'application de toutes ces méthodes au cas du remplissage avec glaçon mobile a enfin permis de valider la robustesse numérique de notre code et le bon couplage de ses différentes entités.

Three-dimensional large-eddy simulations (LES) of single and multiple jets in cross-flow (JICF) were conducted using the 19-bit Lattice Boltzmann Equation (LBE) method coupled with a conventional Navier-Stokes (NS) finite-volume scheme. In this coupled LBE-NS approach, the LBE-LES was employed to simulate the flow inside jet nozzles, while the NS-LES was used to simulate the cross-flow. The key application area was to study the micro-blowing technique (MBT) for drag control similar to recent experiments at NASA/GRC. A single jet in the cross-flow case was used for validation purposes, and results were compared with experimental data and full LBE-LES simulation. Good agreement with data was obtained. Transient analysis of flow structures was performed to investigate the contribution of flow structures to the counter-rotating vortex pair (CRVP) formation. It was found that both spanwise roller (at the lee side of the jet) and streamwise vortices (at the jet-side) contribute to the generation of the CRVP. Span-wise roller at the corner of the jet experiences high spanwise vortex compression as well as high streamwise vortex stretch. As a result, they get realigned, mix with the jet-side streamwise vortices, and eventually generate the CRVP. Furthermore, acoustic pulses were used to test the proper information exchange from the LBE domain to the NS domain, and vice-versa. Subsequently, MBT over a flat plate with porosity of 25 percent was simulated using nine jets in a compressible cross-flow at a Mach number of 0.4. Three cases with injection ratios of 0.003, 0.02 and 0.07 were conducted to investigate how the blowing rate impacts skin friction. It is shown that MBT suppressed the near-wall vortices and reduced the skin friction by up to 50 percent. This is in good agreement with experimental data.

The separation and reattachment of the suction surface boundary layer in a low pressure turbine is characterized using large-eddy simulation at Re=68,000 based on freestream velocity and suction surface length. A high pass filtered Smagorinsky model is used for modeling the sub-grid scales. The onset of time mean separation is at s=so = 0:61 and reattachment at s=so = 0:81, extending over 20% of the suction surface. The boundary layer is convectively unstable with a maximum reverse flow velocity of about 13% of freestream. The breakdown to turbulence occurs over a very short distance of suction surface which is followed by reattachment. Detailed investigations into the structure and kinematics of the bubble and turbulence statistics are presented. The vortex shed from the bubble, convects downstream and interacts with the trailing edge vortices increasing the turbulence intensity. On the suction side, dominant hairpin structures near the transitional and turbulent flow regime are observed. These hairpin vortices are carried by the freestream even downstream of the trailing edge of the blade with a possibility of reaching the next stage. Longitudinal streaks that evolve from the breakdown of hairpin vortices formed near the leading edge are observed on the pressure surface.

The main objective of this dissertation is to develop large-eddy simulation (LES) based computational tools for supersonic inlet and combustor design. In the recent past, LES methodology has emerged as a viable tool for modeling turbulent combustion. LES computes the large scale mixing process accurately, thereby providing a better starting point for small-scale models that describe the combustion process. In fact, combustion models developed in the context of Reynolds-averaged Navier Stokes (RANS) equations exhibit better predictive capability when used in the LES framework. The development of a predictive computational tool based on LES will provide a significant boost to the design of scramjet engines. Although LES has been used widely in the simulation of subsonic turbulent flows, its application to high-speed flows has been hampered by a variety of modeling and numerical issues. In this work, we develop a comprehensive LES methodology for supersonic flows, focusing on the simulation of scramjet engine components. This work is divided into three sections. First, a robust compressible flow solver for a generalized high-speed flow configuration is developed. By using carefully designed numerical schemes, dissipative errors associated with discretization methods for high-speed flows are minimized. Multiblock and immersed boundary method are used to handle scramjet-specific geometries. Second, a new combustion model for compressible reactive flows is developed. Subsonic combustion models are not directly applicable in high-speed flows due to the coupling between the energy and velocity fields. Here, a probability density function (PDF) approach is developed for high-speed combustion. This method requires solution to a high dimensional PDF transport equation, which is achieved through a novel direct quadrature method of moments (DQMOM). The combustion model is validated using experiments on supersonic reacting flows. Finally, the LES methodology is used to study the inlet-isolator component of a dual-mode scramjet. The isolator is a critical component that maintains the compression shock structures required for stable combustor operation in ramjet mode. We simulate unsteady dynamics inside an experimental isolator, including the propagation of an unstart event that leads to loss of compression. Using a suite of simulations, the sensitivity of the results to LES models and numerical implementation is studied. / text

La prédiction des flux aux parois joue un rôle déterminant dans le cycle de vie des chambres de combustion. Le transfert de chaleur de la flamme aux parois est entraîné, outre la convection, également par le rayonnement des gaz chauds au sein de la chambre. Afin d’intégrer les contributions convectives et radiatives au flux pariétal il est nécessaire de résoudre simultanément l’équation de transfert radiatif et les équations régissant l’écoulement réactif. Quand les méthodes de Monte Carlo sont couplées aux simulations aux grandes échelles (LES), de telles simulations deviennent très coûteuses. L’objectif de cette thèse est donc d’investiguer une technique pour améliorer l’efficacité de la méthode MC, basée sur un mécanisme alternatif d’échantillonnage appelée intégration Quasi-Monte Carlo (QMC). Au cours de cette thèse, la méthode QMC a été couplée à une simulation LES dans une configuration où le rayonnement joue un rôle très important : la flamme méthane-air de la chambre Oxytec. La comparaison entre les simulations couplées et non couplées avec les données expérimentales montre que le rayonnement thermique a un impact sur la topologie de l’écoulement et de la flamme. Enfin, un bon accord est trouvé entre le flux de chaleur pariétal prédit par la simulation et les données expérimentales. / The prediction of wall fluxes is a significant aspect in the life cycle of combustors, since it allows to prevent eventual wall damages. Heat transfer from flame to the walls is driven, apart from convection, also by radiation of burnt gases inside the chamber. In order to correctly account for both convective and radiative contributions to wall fluxes, the simultaneous solution of the radiative transfer equation (RTE) and the governing equations for reactive flows is required. However, multi-physics simulations where MC methods are coupled to Large Eddy Simulation (LES), remain very costly. The purpose of this study is then to investigate improvements of MC methods, by using an alternative sampling mechanism for numerical integration usually referred to as Quasi-Monte Carlo (QMC) integration. In this study, QMC method is coupled to Large Eddy Simulation (LES) of a configuration where the radiation plays an important role: the methane-air flame investigated during the experimental campaign Oxytec. Coupled and non-coupled simulations are compared and their comparison with experimental data shows that thermal radiation has an impact on both flow and flame topology. Finally a good agreement is found between numerical wall fluxes and experimental conductive fluxes.

This thesis deals with turbulent mixing processes occuring in internal combustion engines, when applying exhaust gas recirculation (EGR). EGR is a very efficient way to reduce emissions of nitrogen oxides (NOx) in internal combustion engines. Exhaust gases are recirculated and mixed with the intake air of the engine, thus reducing the oxygen concentration of the combustion gas and the maximum combustion tempera- ture. This temperature decrease results in a reduction of NOx emissions, since NOx is produced at high temperatures.The issue of NOx reduction is of high importance for current engine development (particularly for heavy-duty engines), since NOx is the main cause for smog formation and subject to increasingly stronger emission legislation. One of the practical problems when applying EGR is the non-uniformity of the mixture among and inside the cylinders deteriorating the engine and emission performance.The aim of this work is to develop and assess methods suited for the computation of turbulent mixing processes in engine conditions. More specifically, RANS and LES computations are considered. The flow structures responsible for the mixing are analyzed for two different T-junctions and a six-cylinder Scania engine-manifold. Shortcomings and advantages of the applied mixing models are explained.The main results are, that commonly applied scalar flux models for the RANS framework do not predict correct scalar flux directions. In stationary flow, the applied k-ε-model in combination with a gradient-diffusion-model gives too small mixing rates as compared to LES and experiments. Furthermore, the LES computations of the T-junctions show, that Dean vortices occuring due to the curvature of the flow are broken up and dissipated only a few diameters downstream of the junction. The RANS computations do not predict this break-up, giving fundamentally different flow structures and mixing distributions. In pulsating flow, a resonance between the natural stabilities and the pulsation frequency is found by LES results, which could not be predicted by RANS.Computations of the flow in a Scania intake manifold with generic boundary con- ditions indicate, that inlet pulsations are important for the mixing process and that the smoothing effect of URANS is not adequate for accurate mixing computations. LES, on the other hand, is more promising, since it is able to capture the physics of pulsating flows much better. / QC 20111117