La simulation des grandes échelles (SGE ou LES pour large eddy simulation) commence à être très utilisée dans l’industrie. Par résolution directe des structures turbulents de grande tailles, le calcul LES est capable de calculer le bruit générée par la voilure ou de prédire avec précision le décollement de coin dans une configuration très simplifiée du compresseur. L’un des problèmes les plus importants pour effectuer un calcul LES est de fournir des conditions d’entrée avec des champs turbulents.Pour une approche hybride RANS/LES (RANS pour Reynolds Averaged Navier-Stokes), les conditions d’entrée turbulentes pour un calcul LES sont générées à l’aide des solutions fournies par le calcul RANS en amont. Il existe plusieurs méthodes pour générer les conditions d’entrée pour LES. Elles peuvent principalement être classées en deux catégories : 1) simulation avec pré-calcul ; 2) la méthode de turbulence synthétique. La simulation avec pré-calcul consiste à effectuer un calcul LES indépendant pour générer un champ turbulent comme conditions d’entrée pour alimenter le calcul principal. Cette méthode peut obtenir des turbulences de haute qualité, mais elle augmente considérablement le temps de calcul et le stockage des données. Le champ turbulent généré par la méthode de turbulence synthétique exige une « distance de adaptation », pendante laquelle le champ turbulent devient pleinement développé. L’objectif principal pour améliorer ce genre de méthodes est donc de diminuer cette distance nécessaire.Dans cette thèse, la méthode de vortex, qui est une approche de turbulence synthétique, est présentée et améliorée. A travers des expériences numériques, les paramètres de la méthode de vortex améliorée sont systématiquement optimisés. L’application à l’écoulement en canal plan et à couche limite en plaque plane, montrent que la méthode de vortex améliorée génère de manière efficace pour fournir des conditions d’entrée pour LES. Dans le cas de l’écoulement en canal plan, la distance d’adaptation nécessaire pour le rétablissement de la turbulence est d’environ 6 fois la demi-hauteur du canal. Pour le cas de l’écoulement en plaque plane, cette distance est environ 21 fois l’épaisseur de la couche limite. Enfin, dans le but de qualifier la turbulence obtenue par des calculs LES, nous utilisons les coefficients de dissymétrie des dérivées des fluctuations de vitesse, et, nous les introduisons comme un nouveau critère pour la qualité de LES. / Large eddy simulation is becoming an important numerical tool in industry recently. Resolving large scale turbulent motions directly, LES is capable to compute the aeroacoustic noise generated by the airfoil or to precisely capture the corner separation in a linear compressor cascade. The main challenge to perform a LES calculation is to prescribe a realistic unsteady inflow field. For hybrid RANS/LES approaches, inflow conditions for downstream LES region must be generated from the upstream RANS solutions. There exist several methods to generate inflow conditions for LES. They can mainly be divided into two categories: 1) Precursor simulation; 2) Synthetic turbulence methods. Precursor simulation requires to run a separate calculation to generate a turbulent ow or a database to feed the main computation. This kind of methods can generate high quality turbulence. However, it requires heavy extra computing load. Synthetic turbulence methods consist in generating a fluctuating velocity field, and within a short “adaptation distance”, the field get fully developed. So main goal of synthetic turbulence methods is to decrease the required adaptation distance. The vortex method which is a synthetic turbulence method is presented and improved here. Parameters of the improved vortex method are optimized systematically with a series of calculations in this thesis. Applications on channel and flat-plate flows show that the improved vortex method is effective in generating the LES inflow conditions. The adaptation distance required for turbulence recovery is about 6 times the half channel height for channel flow, and 21 times the boundary-layer thickness (at the inlet of vortex) for at-plate ow. The velocity-derivative skewness is used to qualify the generated turbulence, and is introduced as a new criterion of LES calculation.

It is important for practical and legal reasons that water and sediments in disturbed areas around the mining operation should be controlled. The construction of a well-designed drainage system that controls erosion and thus restores the proper hydraulic function of the surface is one of the most important post-disturbance features which should be done as part of the mining activities. However, even with the best planning and design, unless proper construction practices are adapted; both the disturbed and reclaimed areas are very much likely to be susceptible to erosion, sedimentation and stability problems. In order to tackle the problem, guidelines on how to design and construct the drainage system should be well prepared. The main objective of this study was to prepare guidelines for the proper design, construction and monitoring of the water drainage management system in the study areas (Kankberg, Maurliden and Renström-Petiknäs). This report has analysed the results from the outcome of HEC-RAS software for the case study of the new ditch around the Maurliden mine site and integrated with different guidelines. Based on the results of the HEC-RAS, the most common problems in the drainage system have been identified. Moreover the thesis project identified important physical parameters such as cross-sections and slopes of the representative ditch which affect the function of the ditch in the study areas. Hydraulic parameters such as velocity which is very important for designing the type of lining and also Froude number which is very important in identifying the type of flow whether it is super-critical, critical or sub-critical were identified. The latter helps to select the type of guideline to be used between steep slope and mild slope.

This work presents the development and analysis of a system that combines the concepts of light field microscopy and particle image velocimetry (PIV) to measure three-dimensional velocities within a microvolume. Rectanglar microchannels were fabricated with dimensions on the order of 350-950 micrometers using a photolithographic process and polydimethylsiloxane (PDMS). The flow was seeded with fluorescent particles and pumped through microchannels at Reynolds numbers ranging from 0.016 to 0.028. Flow at Reynolds numbers in the range of 0.02 to 0.03 was seeded with fluorescent particles and pumped through microchannels. A light field microscope with a lateral resolution of 6.25 micrometers and an axial resolution of 15.5 micrometers was designed and built based on the concepts described by Levoy et al. Light field images were captured continuously at a frame rate of 3.9 frames per second using a Canon 5D Mark II DSLR camera. Each image was post processed to render a stack of two-dimensional images. The focal stacks were further post processed using various methods including bandpass filtering, 3D deconvolution, and intensity-based thresholding, to remove effects of diffraction and blurring. Subsequently, a multi-pass, three-dimensional PIV algorithm was used to measure channel velocities. Results from PIV analysis were compared with an analytical solution for fully-developed cases, and with CFD simulations for developing flows. Relative errors for fully-developed flow measurements, within the light field microscope refocusing range, were approximately 5% or less. Overall, the main limitations are the reduction in lateral resolution, and the somewhat low axial resolution. Advantages include the relatively low cost, ease of incorporation into existing micro-PIV systems, simple self-calibration process, and potential for resolving instantaneous three-dimensional velocities in a microvolume.

[pt] Modelos RANS (Reynolds Average Navier-Stokes) estão entre os modelos mais empregados para resolver escoamentos turbulentos, devido a seu baixo custo computacional. A maioria dos modelos RANS usa a aproximação de Boussinesq, baseada em uma relação linear entre a parte deviatórica do tensor de Reynolds e o tensor taxa de deformação, com a viscosidade turbulenta sendo o parâmetro positivo de proporcionalidade. Contudo, esses modelos falham em várias situações, e um grande esforço tem sido feito pela comunidade científica com intuito de melhorar a previsibilidade do modelo desenvolvendo modelos não lineares. Análises de modelos de ordem superior empregando tensores ortogonais objetivos têm mostrado que estes são muito promissores para melhorar a previsão dos componentes normais do tensor de Reynolds. No presente trabalho, modelos não lineares baseados no quadrado do tensor taxa de deformação e no tensor não persistência de deformação foram avaliados para uma faixa de número de Reynolds baseados na velocidade de atrito, variando de 395 até 5200. Novas funções de parede foram desenvolvidas, utilizando energia cinética turbulenta e o módulo do tensor taxa de deformação para determinar a velocidade e comprimento característicos. Além disso, um novo modelo turbulento de uma-equação baseado somente na equação de transporte da energia cinética turbulenta foi proposto juntamente com uma equação de fechamento algébrica para modelar a dissipação da energia cinética turbulenta. Os resultados dos modelos para escoamento em canal foram comparados com os dados DNS, apresentando uma melhor aderência aos dados DNS em comparação com os resultados de outros modelos RANS encontrados na literatura. / [en] Reynolds Average Navier Stokes (RANS) models are among the most employed models to solve turbulent flows, due to their low computational cost. The majority of RANS models use the Boussinesq approximation, based on a linear relation between the deviatoric part of Reynolds stress tensor and the rate of strain tensor, with the turbulent viscosity as the positive proportionality parameter. However, these models fail in several situations, and a great deal of effort has been made by the scientific community aiming to improve model prediction through the development of non-linear models. Analysis of higher-order models employing objective orthogonal tensors has shown that these are very promising to improve the prediction of the normal components of the Reynolds stress. In this work, non-linear models based on the square of the rate-strain tensor and non-persistence tensor were examined for a range of friction Reynolds number from 395 to 5200. New wall damping functions were developed, employing the turbulent kinetic energy and intensity of the rate of strain tensor to determine the turbulent characteristic velocity and length. Further, a new one-equation turbulent model based only on the turbulent kinetic energy transport equation was proposed coupled with an algebraic closure equation to model the turbulent kinetic energy dissipation. The models prediction for a channel flow were compared with DNS data and presented a better adherence to the DNS data, than the results of other RANS models available in the literature.

No description available.

The interest towrds machine learning applied to turbulence has experienced a fast-paced growth in the last years. Thanks to deep-learning algorithms, flow-control stratigies have been designed, as well as tools to model and reproduce the most relevant turbulent features. In particular, the success of recurrent neural networks (RNNs) has been demonstrated in many recent studies and applications. The main objective of this project is to assess the capability of these networks to reproduce the temporal evolution of a minimal turbulent channel flow. We first obtain a data-driven model based on a modal decomposition in the Fourier domain (FFT-POD) on the time series sampled from the flow. This particular case of turbulent flow allows us to accurately simulate the most relevant coherent structures close to the wall. Long-short-term-memory (LSTM) networks and a Koopman-based framework (KNF) are trained to predict the temporal dynamics of the minimal channel flow modes. Tests with different configurations highlight the limits of the KNF method compared to the LSTM, given the complexity of the data-driven model. Long-term prediction for LSTM show excellent agreement from the statistical point of view, with errors below 2% for the best models. Furthermore, the analysis of the chaotic behaviour thorugh the use of the Lyapunov exponent and of the dynamic behaviour through Pointcaré maps emphasizes the ability of LSTM to reproduce the nature of turbulence. Alternative reduced-order models (ROMS), based on the identification of different turbulent structures, are explored and they continue to show a good potential in predicting the temporal dynamics of the minimal channel.

Les matériaux absorbants acoustiques, qui sont d’un intérêt stratégique en aéronautique pour la diminution passive du bruit des réacteurs d’avion, conduisent à une physique complexe où l’écoulement turbulent, des ondes acoustiques, et l’absorbant interagissent. Cette thèse porte sur la simulation de cette interaction dans le problème modèle d’un écoulement de canal turbulent avec des parois impédantes, par le biais de simulations numériques aux grandes échelles implicites, dans un contexte de calcul haute performance.Une étude est d’abord faite des grandes échelles dans un canal turbulent avec des parois rigides, en s’intéressant plus particulièrement à l’effet d’une faible compressibilité (Mach <3) sur les caractéristiques de ces échelles.Un canal turbulent avec une paroi de type impédance est ensuite simulé, avec une condition habituelle de périodicité dans le sens de l’écoulement. On observe que pour des faibles valeurs de la résistance et des fréquences de résonance basses, l’écoulement est instable, ce qui engendre une onde le long de l’absorbant, qui modifie la turbulence et augmente la trainée.Enfin, on se tourne vers une simulation de canal spatial en levant la condition de périodicité dans la direction de l’écoulement, ce qui permet d’introduire une onde acoustique en entrée de domaine. L’atténuation de l’onde dans l’écoulement turbulent est étudiée avec des parois rigides, puis un absorbant acoustique est introduit. Dans cette configuration plus réaliste, il est confirmé que l’écoulement peut devenir instable au bord amont de l’absorbant, ce qui empêche l’atténuation de l’onde acoustique incidente. / Acoustic liners are a key technology in aeronautics for the passive reduction of the noise generated by aircraft engines. They are employed in a complex flow scenario in which the acoustic waves, the turbulent flow, and the acoustic liner are interacting.During this thesis, in a context of high performance computing, a compressible Navier-Stokes solver has been developed to perform implicit large eddy simulations of a model problem of this interaction: a turbulent plane channel flow with one wall modeled as an impedance condition.As a preliminary step the wall-turbulence in rigid channel flows and associated large-scale motions are investigated. A straightforward algorithm to detect these flow features is developed and the effect of compressibility on the flow structures and their contribution to the drag are studied. Then, the interaction between the acoustic liner and turbulent flow is investigated assuming periodicity in the streamwise direction. It is shown that low resistance and low resonance frequency tend to trigger flow instability, which modifies the conventional wall-turbulence and also results in drag increase.Finally, the simulation of a spatial channel flow was addressed. In this case no periodicity is assumed and an acoustic wave can be injected at the inlet of the domain. The effect of turbulence on sound attenuation is studied without liner, before a liner is introduced on a part of the channel bottom wall. In this more realistic case, it is confirmed that low resistance acoustic liners trigger an instability at the leading edge of the liner, resulting in drag increase and excess noise generation.

The context of this thesis is the numerical simulation of turbulent flows at moderate Reynolds numbers and the improvement of the capabilities of an in-house 3D unsteady and incompressible flow solver called SFELES to simulate such flows. In addition to this abstract, this thesis includes five other chapters. The second chapter of this thesis presents the numerical methods implemented in the two CFD solvers used as part of this work, namely SFELES and PHASTA. The third chapter concentrates on the implementation of a new library called FlexMG. This library allows the use of various types of iterative solvers preconditioned by algebraic multigrid methods, which require much less memory to solve linear systems than a direct sparse LU solver available in SFELES. Multigrid is an iterative procedure that relies on a series of increasingly coarser approximations of the original 'fine' problem. The underlying concept is the following: low wavenumber errors on fine grids become high wavenumber errors on coarser levels, which can be effectively removed by applying fixed-point methods on coarser levels. Two families of algebraic multigrid preconditioners have been implemented in FlexMG, namely smooth aggregation-type and non-nested finite element-type. Unlike pure gridless multigrid, both of these families use the information contained in the initial fine mesh. A hierarchy of coarse meshes is also needed for the non-nested finite element-type multigrid so that our approaches can be considered as hybrid. Our aggregation-type multigrid is smoothed with either a constant or a linear least square fitting function, whereas the non-nested finite element-type multigrid is already smooth by construction. All these multigrid preconditioners are tested as stand-alone solvers or coupled with a GMRES (Generalized Minimal RESidual) method. After analyzing the accuracy of the solutions obtained with our solvers on a typical test case in fluid mechanics (unsteady flow past a circular cylinder at low Reynolds number), their performance in terms of convergence rate, computational speed and memory consumption is compared with the performance of a direct sparse LU solver as a reference. Finally, the importance of using smooth interpolation operators is also underlined in this work. The fourth chapter is devoted to the study of subgrid scale models for the large eddy simulation (LES) of turbulent flows. It is well known that turbulence features a cascade process by which kinetic energy is transferred from the large turbulent scales to the smaller ones. Below a certain size, the smallest structures are dissipated into heat because of the effect of the viscous term in the Navier-Stokes equations. In the classical formulation of LES models, all the resolved scales are used to model the contribution of the unresolved scales. However, most of the energy exchanges between scales are local, which means that the energy of the unresolved scales derives mainly from the energy of the small resolved scales. In this fourth chapter, constant-coefficient-based Smagorinsky and WALE models are considered under different formulations. This includes a classical version of both the Smagorinsky and WALE models and several scale-separation formulations, where the resolved velocity field is filtered in order to separate the small turbulent scales from the large ones. From this separation of turbulent scales, the strain rate tensor and/or the eddy viscosity of the subgrid scale model is computed from the small resolved scales only. One important advantage of these scale-separation models is that the dissipation they introduce through their subgrid scale stress tensor is better controlled compared to their classical version, where all the scales are taken into account without any filtering. More precisely, the filtering operator (based on a top hat filter in this work) allows the decomposition u' = u - ubar, where u is the resolved velocity field (large and small resolved scales), ubar is the filtered velocity field (large resolved scales) and u' is the small resolved scales field. At last, two variational multiscale (VMS) methods are also considered. The philosophy of the variational multiscale methods differs significantly from the philosophy of the scale-separation models. Concretely, the discrete Navier-Stokes equations have to be projected into two disjoint spaces so that a set of equations characterizes the evolution of the large resolved scales of the flow, whereas another set governs the small resolved scales. Once the Navier-Stokes equations have been projected into these two spaces associated with the large and small scales respectively, the variational multiscale method consists in adding an eddy viscosity model to the small scales equations only, leaving the large scales equations unchanged. This projection is obvious in the case of a full spectral discretization of the Navier-Stokes equations, where the evolution of the large and small scales is governed by the equations associated with the low and high wavenumber modes respectively. This projection is more complex to achieve in the context of a finite element discretization. For that purpose, two variational multiscale concepts are examined in this work. The first projector is based on the construction of aggregates, whereas the second projector relies on the implementation of hierarchical linear basis functions. In order to gain some experience in the field of LES modeling, some of the above-mentioned models were implemented first in another code called PHASTA and presented along with SFELES in the second chapter. Finally, the relevance of our models is assessed with the large eddy simulation of a fully developed turbulent channel flow at a low Reynolds number under statistical equilibrium. In addition to the analysis of the mean eddy viscosity computed for all our LES models, comparisons in terms of shear stress, root mean square velocity fluctuation and mean velocity are performed with a fully resolved direct numerical simulation as a reference. The fifth chapter of the thesis focuses on the numerical simulation of the 3D turbulent flow over a re-entry Apollo-type capsule at low speed with SFELES. The Reynolds number based on the heat shield is set to Re=10^4 and the angle of attack is set to 180º, that is the heat shield facing the free stream. Only the final stage of the flight is considered in this work, before the splashdown or the landing, so that the incompressibility hypothesis in SFELES is still valid. Two LES models are considered in this chapter, namely a classical and a scale-separation version of the WALE model. Although the capsule geometry is axisymmetric, the flow field in its wake is not and induces unsteady forces and moments acting on the capsule. The characterization of the phenomena occurring in the wake of the capsule and the determination of their main frequencies are essential to ensure the static and dynamic stability during the final stage of the flight. Visualizations by means of 3D isosurfaces and 2D slices of the Q-criterion and the vorticity field confirm the presence of a large meandering recirculation zone characterized by a low Strouhal number, that is St≈0.15. Due to the detachment of the flow at the shoulder of the capsule, a resulting annular shear layer appears. This shear layer is then affected by some Kelvin-Helmholtz instabilities and ends up rolling up, leading to the formation of vortex rings characterized by a high frequency. This vortex shedding depends on the Reynolds number so that a Strouhal number St≈3 is detected at Re=10^4. Finally, the analysis of the force and moment coefficients reveals the existence of a lateral force perpendicular to the streamwise direction in the case of the scale-separation WALE model, which suggests that the wake of the capsule may have some preferential orientations during the vortex shedding. In the case of the classical version of the WALE model, no lateral force has been observed so far so that the mean flow is thought to be still axisymmetric after 100 units of non-dimensional physical time. Finally, the last chapter of this work recalls the main conclusions drawn from the previous chapters.

Large-eddy simulation (LES) and unsteady Reynolds-averaged Navier-Stokes (URANS) calculations have been performed to investigate the effects of different mathematical models for scalar variance and its dissipation rate as applied to both a non-reacting bluff-body turbulent flow and an extension to a reacting case. In the conserved scalar formalism, the mean value of a thermo-chemical variable is obtained through the PDF-weighted integration of the local description over the conserved scalar, the mixture fraction. The scalar variance, one of the key parameters for the determination of a presumed β-function PDF, is obtained by solving its own transport equation with the unclosed scalar dissipation rate modelled using either an algebraic expression or a transport equation. The proposed approach is first applied to URANS and then extended to LES. Velocity, length and time scales associated with the URANS modelling are determined using the standard two-equation k-ε transport model. In contrast, all three scales required by the LES modelling are based on the Smagorinsky subgrid scale (SGS) algebraic model. The present study proposes a new algebraic and a new transport LES model for the scalar dissipation rate required by the transport equation for scalar variance, with a time scale consistent with the Smagorinsky SGS model.

From the Proceedings of the 1975 Meetings of the Arizona Section - American Water Resources Assn. and the Hydrology Section - Arizona Academy of Science - April 11-12, 1975, Tempe, Arizona / The importance of water resources and an increasing interest on improvement of out environment have identified the urgent need for methods to predict river response due to various changes resulting from proposed water resource planning. Fluvial geomorphology and hydraulic elements that are related to the interpretation and modeling of response to the problem are presented. Interpretation of alluvial rivers should be preceded by a qualitative analysis and information is presented which should be adequate to carry this out in most cases. This should be followed by a quantitative evaluation of channel response and water sediment routing using theory supplemented by physical and mathematical model studies of the system.