No description available.

Three-dimensional (3D) turbulent boundary layers (TBLs) play a crucial role in determining the aerodynamic properties of most aero-mechanical devices. However, accurately predicting these flows remains a challenge due to the complex nonlinear and nonlocal physics involved, which makes it difficult to develop universally applicable models. This limitation is particularly significant as the industry increasingly relies on simulations to make decisions in high-consequence environments, such as the certification or aircraft, and high-fidelity simulation methods that don't rely on modeling are prohibitively expensive. To address this challenge, it is essential to gain a better understanding of the physics underlying 3D TBLs. This research aims to improve the predictive accuracy of turbulence models in 3D TBLs by examining the impact of model assumptions underpinning turbulent constitutive relations, which are fundamental building blocks of every turbulence model. Specifically, the study focuses on the relevance and necessity of nonlinear and nonlocal model assumptions for accurately predicting 3D TBLs. The study considers the attached 3D boundary layer flow over the textbf{Be}nchmark textbf{V}alidation textbf{E}xperiment for textbf{R}ANS/textbf{L}ES textbf{I}nvestiagtions (BeVERLI) Hill as a test case and corresponding particle image velocimetry data for the investigation. In a first step, the BeVERLI Hill experiment is comprehensively described, and the important characteristics of the flow over the BeVERLI Hill are elucidated, including complex symmetry breaking characteristics of this flow. Reynolds-averaged Navier-Stokes simulations of the case using standard eddy viscosity models are then presented to establish the baseline behavior of local and linear constitutive relations, i.e., the standard Boussinesq approximation. The tested eddy viscosity models fail in the highly accelerated hill top region of the BeVERLI hill and near separation. In a further step, several nonlinear and nonlocal turbulent constitutive relations, including the QCR model, the model by Gatski and Speziale, and the difference-quotient model by Egolf are used as metrics to gauge the impact of nonlinearities and nonlocalities for the modeling of 3D TBLs. It is shown that nonlinear and nonlocal approaches are essential for effective 3D TBL modeling. However, simplified reduced-order models could accurately predict 3D TBLs without high computational costs. A constitutive relation with local second-order nonlinear mean strain relations and simplified nonlocal terms may provide such a minimal model. In a final step, the structure and response of non-equilibrium turbulence to continuous straining are studied to reveal new scaling laws and structural models. / Doctor of Philosophy / Airplanes and other flying objects rely on the way air flows around them to generate lift and stay in the sky. This airflow can be very complex, especially close to the surface of the object, where it is affected by friction with the object. This friction generates a layer of air called a boundary layer, which can become turbulent and lead to complex patterns of airflow. The boundary layer is generated by the friction between the air and the surface of the object, which causes the air molecules to "stick" to the surface. This sticking creates a layer of slow-moving air that slows down the flow of air around the object. This loss of momentum creates drag, which is one of the main factors that resist the motion of objects in the air. The slowing of the air flow in the boundary layer is due to the viscosity of the air, which is a measure of how resistant the air is to deformation. The molecules in the air have a tendency to stick together, making it difficult for them to move past each other. This resistance causes the momentum of the air to be lost as it flows over the surface of the object because air molecules close to the surface "pull" on the ones farther away. Understanding how turbulent boundary layers (TBLs) work is essential to accurately predict the airflow around these objects using computer simulations. However, it's challenging because TBLs involve complex physics that are difficult to model accurately. This research focuses on a specific type of TBL called a three-dimensional (3D) TBL. This study looks at how different assumptions affect the accuracy of computer simulations that predict this type of airflow. It is found that using more complex models that take into account nonlinear and nonlocal physics can help predict 3D TBLs more accurately. However, these models are computationally expensive, and it is also found that simpler models can work well enough and are cheaper. This research further establishes important physical relations of the mechanisms pertaining 3D TBLs that could support the advancement of current models.

Convergent-Divergent (C-D) nozzles are used in rocket engines to produce thrust as a reaction to the acceleration of hot combustion chamber gases in the opposite direction. To maximize the engine performance at high altitudes, large area ratio, bell-shaped or contoured nozzles are used. At lower altitudes, the exit pressure of these nozzles is lower than the ambient pressure. During this over-expanded condition, the nozzle-internal flow adapts to the ambient pressure through an oblique shock. But the boundary layer inside the divergent portion of the nozzle is unable to withstand the pressure rise associated with the shock, and consequently flow separation is induced. Numerical simulation of separated flows in rocket nozzles is challenging because the existing turbulence models are unable to correctly predict shock-induced flow separation. The present thesis addresses this problem. Axisymmetric, steady-state, Reynolds-Averaged Navier-Stokes (RANS) simulations of a conical nozzle and three sub-scale contoured nozzles were carried out to numerically predict flow separation in over-expanded rocket nozzles at different nozzle pressure ratios (NPR). The conical nozzle is the JPL 45◦-15◦ and the contoured nozzles are the VAC-S1, the DLR-PAR and the VAC-S6-short. The commercial CFD code ANSYS FLUENT 13 was first validated for simulation of separated cold gas flows in the VAC-S1 nozzle. Some modeling issues in the numerical simulations of flow separation in rocket nozzles were determined. It is recognized that compressibility correction, nozzle-lip thickness and upstream-extension of the external domain are the sources of uncertainty, besides turbulence modeling. In high-speed turbulent flows, compressibility is known to affect dissipation rate of turbulence kinetic energy. As a consequence, a reduction in the spreading rate of supersonic mixing layers occurs. Whereas, the standard turbulence models are developed and calibrated for incompressible flows and hence, do not account for this effect. ANSYS FLUENT uses the compressibility correction proposed by Wilcox [1] which modifies the turbulence dissipation terms based on turbulent Mach number. This, as shown in this thesis, may not be appropriate to the prediction of flow separation in rocket nozzles. Simulation results of the standard SST model, with and without the compressibility correction, are compared with the experimental data at NPR=22 for the DLR-PAR nozzle. Compressibility correction is found to cause under-prediction of separation location and hence its use in the prediction of flow separation is not recommended. In the literature, computational domains for the simulation of DLR subscale nozzles have thick nozzle-lips whereas for the VAC subscale nozzles they have no nozzle-lip. Effect of nozzle-lip thickness on flow separation is studied in the DLR-PAR nozzle by varying its nozzle-lip thickness. It is found that nozzle-lip thickness significantly influences both separation location and post-separation pressure recovery by means of the recirculation bubbles formed at the nozzle-lip. Usually, experimental values of free stream turbulence are unknown. So conventionally, to minimize solution dependence on the boundary conditions specified for the ambient flow, the computational domain external to the nozzle is extended in the upstream direction. Its effect on flow separation is studied in the DLR-PAR nozzle through simulations conducted with and without this domain extension. No considerable effect on separation location and pressure recovery is found. The two eddy-viscosity based turbulence models, Spalart-Allmaras (SA) model and Shear Stress Transport (SST) model, are well known to predict separation location better than other eddy-viscosity models, but with moderate success. Their performances, in terms of predicting separation location and post-separation wall pressure distribution, were compared with each other and evaluated against experimental data for the conical and two contoured nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the range of NPRs and to the type of nozzle. Depending on NPR, the SST model either under-predicts or over-predicts Free Shock Separation (FSS). Moreover, it also fails to capture Restricted Shock Separation (RSS). With compressibility correction, it under-predicts separation at all NPRs to a greater extent. Even though RSS is captured by using compressibility correction, the transition from FSS to RSS is over-predicted [2]. Early efforts by few researchers to improve predictions of nozzle flow separation by realizability corrections to turbulence models have not been successful, especially in terms of capturing both the separation types. Therefore, causes of turbulence modeling failure in predicting nozzle flow separation correctly were further investigated. It is learnt that limiting of the shear stress inside boundary layer, due to Bradshaw’s assumption, and over-prediction of jet spreading rate are the causes of SST model’s failure in predicting nozzle flow separation correctly. Based on this physical reasoning, values of the a 1 parameter and the two diffusion coefficients σk,2 and σω,2 were empirically modified to match the predicted wall pressure distributions with experimental data of the DLR-PAR and the VAC-S6-short nozzles. The results confirm that accurate prediction of flow separation in rocket nozzles indeed depends on the correct prediction of spreading rate of the supersonic separation-jet. It is demonstrated that accurate RANS simulation of flow separation in rocket nozzles over a wide range of NPRs is feasible by modified values of the diffusion coefficients in turbulence model.

Cette thèse propose une modélisation numérique 2D des écoulements en zone de jet de rive avec un code Volume-Of-Fluid. Dans un premier temps, le détail de la structure interne de l’écoulement de jet de rive généré par l’effondrement d’un front d’onde turbulent sur une plage imperméable lisse est étudié. Le modèle numérique est ensuite étendu à la simulation des écoulements en milieu poreux internes à la plage. L’utilisation d’une unique équation de quantité de mouvement (VARANS) et de la méthode 1-fluide, permet de résoudre simultanément les écoulements de l’eau et de l’air à la surface et dans une plage perméable. Ce modèle a été confronté à une série de cas tests analytiques et à de récentes mesures expérimentales. Les résultats numériques montrent l’aptitude du modèle VOF-VARANS à reproduire les écoulements en zone de jet de rive sur une plage imperméable fixe. / A 2D numerical modeling of flows in the swash zone is proposed using a Volume-Of- Fluid code. The detailed flow structure of a bore-driven swash event over an impermeable beach is first studied. The numerical model is then developed to account for porous media flow within the beach. The unique VARANS momentum equation and 1-fluid method used allow to solve simultaneously both surface and subsurface flows of air and water phases in the swash zone. This model is validated against a series of analytical tests cases and confronted to recent experimental measurements. The numerical results highlight the ability of the VOF-VARANS model to reproduce swash flows over and within a permeable beach.

The primary requirements for a modern industrial gas turbine consist of a continuous trend of an increasing efficiency combined with very low emissions in a robust, cost-effective manner. To fulfil these tasks a high turbine inlet temperature together with advanced dry low NOX combustion chambers are employed. These dry low NOX combustion chambers generate a rather flat temperature profile compared to previous generation gas turbines, which have a rather parabolic temperature profile before the nozzle guide vane. This means that the nozzle guide vane endwall heat load for modern gas turbines is much higher compared to previous generation gas turbines. Therefore the prediction of the nozzle guide vane flow field and endwall heat transfer is crucial for the engineering task of the design layout of the vane endwall cooling system. The present study is directed towards establishing new in-depth aerodynamic and endwall heat transfer knowledge for an advanced nozzle guide vane of a modern industrial gas turbine. To reach this objective the physical processes and effects which cause the different flow fields and the endwall heat transfer pattern in a baseline configuration, a combustion chamber variant, a heat shield variant without and with additional cooling air and a cavity variant without and with additional cooling air have been investigated. The variants, which differ from the simplified baseline configuration, apply design elements which are commonly used in real modern gas turbines. This research area is crucial for the nozzle guide vane endwall heat transfer, especially for the advanced design of the nozzle guide vane of a modern industrial gas turbine and has so far hardly been investigated in the open literature. For the experimental aerodynamic and endwall heat transfer research of the baseline configuration of the advanced nozzle guide vane geometry a new low pressure, low temperature test facility has been developed, designed and constructed, since no experimental heat transfer data exist in the open literature for this type of vane configuration. The new test rig consists of a linear cascade with the baseline configuration of the advanced nozzle guide vane geometry with four upscaled airfoils and three flow passages. For the aerodynamic tests the two middle airfoils and the hub and the tip endwall are instrumented with pressure taps to monitor the Mach number distribution. For the heat transfer tests the temperature distribution on the hub endwall is measured via thermography. The analysis of these measurements, including comparisons to research in the open literature shows that the new test rig generates accurate and reproducible results which give confidence that it is a reliable tool for the experimental aerodynamic and heat transfer research on the advanced nozzle guide vane of a modern industrial gas turbine. Previous own research work together with the numerical analysis performed in another part of the project as well as conclusions from a detailed literature study lead to the conclusion that advanced Navier-Stokes CFD tools with the v2-f turbulence model are most suitable for the calculation of the flow field and the endwall heat transfer of turbine vanes and blades. Therefore this numerical tool, validated against different vane and blade geometries and for different flow conditions, has been chosen for the numerical aerodynamic and endwall heat transfer research of the advanced nozzle guide vane of a modern industrial gas turbine. The evaluation of the numerical and experimental investigations of the baseline configuration of the advanced design of a nozzle guide vane shows the flow field of an advanced mid-loaded airfoil design with the features to reduce total airfoil losses. For the hub endwall of the baseline configuration of the advanced design of a nozzle guide vane the flow characteristics and heat transfer features of the classical vane endwall secondary flow model can be detected with a very weak intensity and geometric extension compared to the studies of less advanced vane geometries in the open literature. A detailed analysis of the numerical simulations and the experimental data showed very good qualitative and quantitative agreement for the three-dimensional flow field and the endwall heat transfer. These findings, together with the evaluations obtained from the open literature, lead to the conclusions that selected CFD software Fluent together with the applied v2-f turbulence model exhibits a high level of general applicability and is not tuned to a special vane or blade geometry. Therefore the CFD code Fluent with the v2-f turbulence model has been selected for the research of the influence of the several geometric variants of the baseline configuration on the flow field and the hub endwall heat transfer of the advanced nozzle guide vane of a modern industrial gas turbine. Most of the vane endwall heat transfer research in the open literature has been carried out only for baseline configurations of the flow path between combustion chamber and nozzle guide vane. Such a simplified geometry consists of a long, planar undisturbed approach length upstream of the nozzle guide vane. The design of real modern industrial gas turbines however requires often significant variations from this baseline configuration consisting of air-cooled heat shields and purged cavities between the combustion chamber and the nozzle guide vane. A detailed evaluation of the flow field and the endwall heat transfer shows major differences between the baseline and the heat shield configuration. The heat shield in front of the airfoil of the nozzle guide vane influences the secondary flow field and the endwall heat transfer pattern strongly. Additional cooling air, released under the heat shield has a distinctive influence as well. Also the cavity between the combustion chamber and the nozzle guide vane affects the secondary flow field and the endwall heat transfer pattern. Here the influence of additional cavity cooling air is more decisive. The results of the detailed studies of the geometric variants are applied to formulate guidelines for an optimized design of the flow path between the combustion chamber and the nozzle guide vane and the nozzle guide vane endwall cooling configuration of next-generation industrial gas turbines. / QC 20100917

Ce travail est une contribution à la simulation numérique d'écoulements turbulents autour de corps non profilés. Après avoir précisé les ingrédients numériques et les modèles de turbulence utilisés dans nos simulations, nous présentons une étude sur l'évaluation des effets de la procédure dynamique des modèles de sous-maille dans un modèle VMS-LES et une approche hybride RANS/VMS-LES. Des problèmes d'écoulements autour d'un cylindre seul et en tandem sont considérés. Nous étudions ensuite le comportement de modèles de turbulence hybrides pour la simulation d'écoulements en régime sous-critique autour d'un cylindre circulaire. Le calcul de l'écoulement autour d'un cylindre de section rectangulaire par l'approche VMS-LES est aussi présenté. Enfin, dans une dernière partie, après avoir fait une revue des travaux importants sur les schémas d'avancement en temps multirate, nous proposons une nouvelle approche explicite multirate par agglomération de volumes finis que nous appliquons à des calculs d'écoulements turbulents complexes en utilisant un modèle de turbulence hybride. / This work is a contribution to the numerical simulation of turbulent flows around bluff bodies. After specifying the numerical ingredients and the turbulence models used in our simulations, we present a study on the impact of the dynamic sub-grid scale modeling in VMS-LES model and a RANS/VMS-LES hybrid turbulence approach. Simulations of flows around a cylinder and a tandem are performed. Next, we assess the behaviour of some hybrid turbulence models for the simulation of flows around a circular cylinder in the subcritical regime. The computation of the flow around a rectangular cylinder with the VMS-LES approach is also presented. At last, after a review of some important works on multirate time advancing schemes, we propose a new volume-agglomeration explicit multirate approach that is applied to the computation of complex turbulent flows by a hybrid turbulence model.

Turbulence is an open and challenging problem for mathematical approaches, physical modeling and numerical simulations. Numerical solutions contribute significantly to the understand of the nature and effects of turbulence. The focus of this thesis is the development of appropriate numerical methods for the computer simulation of turbulent flows. Many of the existing approaches to turbulence utilize analogies from kinetic theory. Degond & Lemou (J. Math. Fluid Mech., 4, 257-284, 2002) derived a k-✏ type turbulence model completely from kinetic theoretic framework. In the first part of this thesis, a numerical method is developed for the computer simulation based on this model. The Boltzmann equation used in the model has an isotropic, relaxation collision operator. The relaxation time in the collision operator depends on the microscopic turbulent energy, making it difficult to construct an efficient numerical scheme. In order to achieve the desired numerical efficiency, an appropriate change of frame is applied. This introduces a stiff relaxation source term in the equations and the concept of asymptotic preserving schemes is then applied to tackle the stiffness. Some simple numerical tests are introduced to validate the new scheme. In the second part of this thesis, alternative approaches are sought for more efficient numerical techniques. The Lattice Boltzmann Relaxation Scheme (LBRS) is a novel method developed recently by Rohan Deshmukh and S.V. Raghuram Rao for simulating compressible flows. Two different approaches for the construction of implicit sub grid scale -like models as Implicit Large Eddy Simulation (ILES) methods, based on LBRS, are proposed and are tested for Burgers turbulence, or Burgulence. The test cases are solved over a largely varying Reynolds number, demonstrating the efficiency of this new ILES-LBRS approach. In the third part of the thesis, as an approach towards the extension of ILES-LBRS to incompressible flows, an artificial compressibility model of LBRS is proposed. The modified framework, LBRS-ACM is then tested for standard viscous incompressible flow test cases.

The purpose of this diploma thesis was to find optimal operational parameters of Ranque–Hilsch vortex tube that would give the best results of temperature separation, and to create a numerical model of the device. Firstly, extensive research of current literature was done and analytical model of the tube was created. Then, the numerical model of the vortex tube was designed using Star-CCM+ software. Afterward, best fitting turbulence model was chosen to do the calculation and optimal geometrical parameters of the tube were obtained as a result of CFD simulation. Then, inlet nozzles for the tube were designed and manufactured. Experimental evaluation of the vortex tube and description of its optimal settings form a substantial part of the project. Finally, the data obtained from the experiment were compared to the results of numerical analysis and conclusions were deduced.

This dissertation investigates the effect of electromagnetic braking and gas injection on the fluid flow in a continuous casting slab mold numerically and makes verifications on basis of a small Liquid Metal Model for Continuous Casting of steel (mini-LIMMCAST). Numerical calculations were performed by means of the software package CFX with an implemented RANS-SST turbulence model. The non-isotropic nature of the MHD turbulence was taken into account by specific modifications of the turbulence model. The numerical results were validated by flow measurements at the mini-LIMMCAST facility. Numerical simulations disclose the damping effect on the flow closely depending on the wall conductance ratio. In addition, specific modifications of the turbulence model play a crucial role in reconstructing the peculiar phenomenon of an excitation of nonsteady, nonisotropic, large-scale flow perturbations caused by the application of the DC magnetic field.

NOWADAYS, a number of studies keep on demonstrating the existence of a strong relation between high concentrations of particulate matter (PM) and the prevalence of human morbidity and mortality. Large particles can be filtered in the nose or in the throat, while fine particles (about10 micrometer) can settle in the bronchi and lungs, leading to more serious consequences. According to Karagulian et al. the major sources of urban air pollution are traffic (25%), combustion and agriculture (22%), domestic fuel burning (20%), natural dust (18%) and industrial activities (15%).As a consequence, the detailed study of dispersion phenomena within the urban canopy becomes a target of great interest. To this end, Computational Fluid Dynamics (CFD) can be successfully employed to predict turbulence and dispersion patterns, accounting for a detailed characterization of the pollutant sources, complex obstacles and atmospheric stability classes.Despite being intrinsically different phenomena, turbulence and dispersion are closely related. It is universally accepted that, to reach accurate prediction of the concentration field, it is necessary to properly reproduce the turbulence one. For this reason, the present PhD thesis is split into two main Sections: one focused on turbulence modelling and the subsequent, centered on the dispersion modelling.Thanks to its good compromise between accuracy of results and calculation time, Reynolds-averaged Navier-Stokes (RANS) still represents a valid alternative to more resource-demanding methods. However, focusing on the models’ performance in urban studies, Large Eddy Simulation (LES) generally outperforms RANS results, even if the former is at least one order of magnitude more expensive. Stemming from this consideration, the aim of this work is to propose a variety of approaches meant to solve some of the major limitations linked to standard RANS simulation and to further improve its accuracy in disturbed flow fields, without renouncing to its intrinsic feasibility. The proposed models are suitable for the urban context, being capable of automatically switching from a formulation proper for undisturbed flow fields to one suitable for disturbed areas. For neutral homogeneous atmospheric boundary layer (ABL), a comprehensive approach is adopted, solving the issue of the erroneous stream-wise gradients affecting the turbulent profiles and able to correctly represent the various roughness elements. Around obstacles, more performing closures are employed. The transition between the two treatments is achieved through the definition of a Building Influence Area (BIA). The finalgoal is to offer more affordable alternatives to LES simulations without sacrificing a good grade of accuracy.Focusing on the dispersion modelling framework, there exists a number of parameters which have to be properly specified. In particular, the definition of the turbulent Schmidt number Sct, expressing the ratio of turbulent viscosity to turbulent mass diffusivity, is imperative. Despite its relevance, the literature does not report a clear guideline on the definition of this quantity. Nevertheless, the importance of Sct with respect to dispersion is undoubted and further demonstrated in the works of different authors. For atmospheric boundary layer flows, typical constant values range between 0.2 and 1.3. As a matter of fact, the local variability of Sct is supported by experimental evidence and by direct numerical simulations (DNS). These observations further suggest that the turbulent Schmidt number should be prescribed as a dynamic variable. Following these observations a variable turbulent Schmidt number formulation is proposed in this work. The latter stems from the same hypothesis of the variable formulation developed by Gorlé et al. Moreover, the relevant uncertain model parameters are optimized through uncertainty quantification (UQ). This formulation further increased the accuracy of the predictions, and was successfully verified by Di Bernardino et al. However, the turbulent Schmidt number resulting from this formulation is still intrinsically linked to the turbulence model employed, i.e. to the Cμ coefficient. To overcome this constraint, the nature and the dependencies of Sct were further analyzed through correlation studies and employing principal component analysis (PCA) on data obtained through the proposed ABL RANS model. Subsequently, the same data-driven technique was employed based on the high-fidelity outcomes of a delayed Detached Eddy Simulation (dDES) to derive a generalized turbulentSchmidt number formulation. The latter can be employed within a wide range of turbulence models, without limiting its variability. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished