Large-eddy simulation (LES) of wind and wave forced oceanic turbulent boundary layers is performed using the residual-based variational multiscale method (RBVMS) and near-wall modeling. Wind and surface gravity wave forcing generates Langmuir turbulence characterized by Langmuir circulation (LC) with largest scales consisting of streamwise vortices aligned in the direction of the wind, acting as a secondary flow structure to the primary wind-driven component of the flow. The LES here is representative of a shallow water continental shelf flow (10 to 30 meters in depth) far from lateral boundaries in which LC engulfs the full depth of the water column and disrupts the bottom log layer. Field measurements indicate that occurrence of full-depth LC is typical during the passage of storms. The RBVMS method with quadratic NURBS (Non-Uniform Rational B-splines) with near-wall resolution is shown to possess good convergence characteristics for this flow. The use of near-wall modeling facilitates simulations with expanded domains over horizontal directions. Thus, these simulations are able to resolve multiple Langmuir cells permitting analysis of the interaction between the cells. Results in terms of velocity statistics are presented from simulations performed with various domain sizes and distinct near-wall treatments: (1) the classical treatment based on prescription of the wall shear stress assuming a law of the wall and (2) a recent treatment based on weak imposition of the no-slip bottom boundary condition.

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.

In this thesis study, Detached Eddy Simulation turbulence model is studied in two dimension mainly for flow over single element airfoils in high Reynolds numbers to gain experience with model before applying it to a three dimensional simulations. For this aim, Spalart-Allmaras and standard DES ,DES97, turbulence models are implemented to parallel, viscous, hybrid grid flow solver. The flow solver ,Set2d, is written in FORTRAN language. The Navier-Stokes equations are discretized by first order accurately cell centered finite volume method and solved explicitly by using Runge-Kutta dual time integration technique. Inviscid fluxes are computed using Roe flux difference splitting method. The numerical simulations are performed in parallel environment using domain decomposition and PVM library routines for inter-process communications. To take into account the effect of unsteadyness after the convergence is ensured by local time stepping technique for four order magnitude drop in density residual, global time stepping is applied for 20000 iterations. The solution algorithm is validated aganist the numerical and experimental studies for single element airfoils in subsonic and transonic flows. It is seen that Spalart-Allmaras and DES97 turbulence models give the same results in the non-seperated flows. Grey area is investigated by changing $C_{DES}$ coefficient. Modeled Stress Depletion which cause reduction of eddy viscosity is observed.

This thesis deals with Large Eddy Simulation (LES) of impinging air jets. The impinging jet configuration features heated circular jets impinging onto a flat plate. The problem addressed here is of generic nature, with applications in many engineering devices, such as cooling of components in gas turbines, in cars and electronic devices. The flow is inherently unsteady and contains relatively slowly varying coherent structures. Therefore, LES is the method of choice when the Reynolds number is large enough to exclude Direct Numerical Simulations (DNS). The present LES model is a basic model without explicit Sub-Grid-Scale (SGS) modeling and without explicit filtering. Instead, the numerical scheme is used to account for the necessary amount of dissipation. By using the computational grid as a filter the cutoff wavenumber depends explicitly on the grid spacing. The underlying computational grid is staggered and constructed in a Cartesian coordinate system. Heat transfer is modeled by the transport equation for a passive scalar. This is possible due to the negligible influence of buoyancy which implies constant density throughout the flow field. The present method provides accurate results for simple geometries in an efficient manner. A great variety of inlet conditions have been considered in order to elucidate how the dynamics of the flow and heat transfer are affected. The considered studies include top-hat and mollified mean velocity profiles subjected to random and sinusoidal perturbations and top-hat profiles superimposed with solid body rotation. It has been found that the shape of the mean inlet velocity profile has a decisive influence on the development of the flow and scalar fields, whereas the characteristics of the imposed artificial disturbances (under consideration) have somewhat weaker effect. In order to obtain results unequivocally comparable to experimental data on turbulent impinging jets both space and time correlations of the inflow data must be considered, so also the spectral content. This is particularly important if the region of interest is close to the velocity inlet, i.e. for small nozzle-to-plate spacings. Within this work mainly small nozzle-toplate spacings are considered (within the range of 0.25 and 4 nozzle diameters), which emphasizes the importance of the inflow conditions. Thus, additional to the basic methods also turbulent inflow conditions, acquired from a precursor pipe simulation, have been examined. Both for swirling and non-swirling flows. This method emulates fully developed turbulent pipe flow conditions and is the best in the sense of being well defined, but it demands a great deal of computing power and is also rather inflexibility. In case of the basic randomly perturbed methods the top-hat approach has been found to produce results in closest agreement with those originating from turbulent inlet conditions. In the present simulations the growth of individual instability modes is clearly detected. The character of the instability is strongly influenced by the imposed boundary conditions. Due to the lack of correlation random superimposed fluctuations have only a weak influence on the developing flow field. The shape of the mean profile, on the other hand, influences both the growth rate and the frequency of the dominant modes. The top-hat profile yields a higher natural frequency than the mollified. Furthermore, for the top-hat profile coalescence of pairs of vortices takes place within the shear-layer of the axial jet, whereas for the mollified profile (for the considered degree of mollification) it takes place within the wall jet. This indicates that the transition process is delayed for smoother profiles. The amount of wall heat transfer is directly influenced by the character of the convective vortical structures. For the mollified cases wall heat transfer originates predominantly from the dynamics of discrete coherent structures. The influence from eddy structures is low and hence Reynolds analogy is applicable, at least in regions of attached flow. The top-hat and the turbulent inflow conditions yield a higher rate of incoherent small scale structures. This strongly affects the character of wall heat transfer. Also the applied level of swirl at the velocity inlet has significant influence on the rate of heat transfer. The turbulence level increases with swirl, which is positive for heat transfer, and so also the spreading of the jet. The latter effect has a negative influence on wall heat transfer, particularly in the center most regions. This however depends also on the details of the inflow data. / QC 20100831

<p>Large-scale coherent structures are systematically investigated in terms of their geometric attributes, importance toward describing turbulent exchange of energy, momentum and mass as well as their relationship to landscape features in the context of land-atmosphere interaction. In the first chapter, we present the motivation of this work as well as a background review of large-scale coherent structures in land-atmosphere interaction. In the second chapter, the methodology of large-eddy simulation (LES) and the proper orthogonal decomposition (POD) is introduced. LES was used to serve as a virtual laboratory to simulate typical scenarios in land-atmosphere interaction and the POD was used as the major technique to educe the coherent structures from turbulent flows in land-atmosphere interaction. In the third chapter, we justify the use of the LES to simulate the realistic coherent structures in the atmospheric boundary layer (ABL) by comparing results obtained from LES of the ABL and direct numerical simulation (DNS) of channel flow. In the fourth chapter, we investigate the effects of a wide range of vegetation density on the coherent structures within the air space within and just above the canopy (the so-called canopy sublayer, CSL). The fifth chapter presents an analysis of the coherent structures across a periodic forest-clearing-forest transition in the steamwise direction. The sixth chapter focuses on the role of coherent structures in explaining scalar dissimilarity in the CSL. The seventh chapter summarizes this dissertation and provides suggestions for future study.</p> / Dissertation

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.

Two new models to calculate the species instantaneous and filtered reaction rates for multi-step, multi-species chemical kinetics mechanisms are developed based on the artificial neural networks (ANN) approach. The proposed methodologies depend on training the ANNs off-line on a thermo-chemical database representative of the actual composition and turbulence level of interest. The thermo-chemical database is constructed by stand-alone linear eddy mixing (LEM) model simulations under both premixed and non-premixed conditions, where the unsteady interaction of turbulence with chemical kinetics is included as a part of the training database. In this approach, the information regarding the actual geometry of interest is not needed within the LEM computations. The developed models are validated extensively on the large eddy simulations (LES) of (i) premixed laminar-flame-vortex-turbulence interaction, (ii) temporally mixing non-premixed flame with extinction-reignition characteristics, and (iii) stagnation point reverse flow combustor, which utilizes exhaust gas re-circulation technique. Results in general are satisfactory, and it is shown that the ANN provides considerable amount of memory saving and speed-up with reasonable and reliable accuracy. The speed-up is strongly affected by the stiffness of the reduced mechanism used for the computations, whereas the memory saving is considerable regardless.

A cause de sa température élevée, le jet propulsif d'un avion de combat est une source de rayonnement infrarouge qui rend l'appareil très vulnérable. Ce rayonnement peut toutefois être réduit en accélérant le mélange du jet avec l'atmosphère. <br />Cette thèse est consacrée à la simulation numérique d'un jet propulsif réaliste contrôlé par des injections radiales et à l'analyse physique des mécanismes d'augmentation de son mélange. <br />Deux types de simulations, RANS et ZDES, ont été réalisés sur la base du modèle de Spalart-Allmaras. Dans le modèle ZDES, une nouvelle longueur caractéristique de maille est formulée et améliore sensiblement la prévision de la région initiale du jet. Globalement, les simulations ZDES restituent fidèlement le champ moyen du jet supersonique sans et avec contrôle, aussi bien les cellules de détente/compression que la diffusion turbulente. <br />L'analyse physique est dédiée à la compréhension d'une part des mécanismes compressibles concentrées au coeur du jet et d'autre part des mécanismes tourbillonnaires périphériques ainsi qu'à l'évaluation de leurs rôles respectifs dans l'augmentation du mélange. Il en ressort que l'augmentation du mélange est exclusivement due aux mécanismes tourbillonnaires. Une étude paramétrique fournit des indications pour concevoir un mélangeur efficace. L'analyse des tourbillons focalisée sur le régime lointain quasi-bidimensionnel souligne leur dynamique moyenne et fait apparaître l'action des fluctuations turbulentes sur leur taux de dégénérescence. Enfin, deux régimes de contrôle sont identifiés et associés aux pénétrations respectivement quasi-stationnaire et intermittente des jets secondaires.

The major source of aircraft noise at take-off is jet noise. If jet noise is not adequately addressed environmental impact concerns will constrain the planned growth of the air transport system. A considerable amount of research worldwide has therefore been aimed at identifying ways to reduce jet noise including development of a predictive tool that can estimate the noise generated by new nozzle designs. Current noise prediction techniques, however, still require the input of empirically calibrated noise source models and their performance is still inadequate. In addition, development of detailed noise source identification measurements and the associated understanding of how to control (and reduce) the noise at the source has been limited. The fundamental turbulence property which acts as the source of propagating noise in shear layers is the two-point space-time velocity correlation (Rijkl). Very few measurements exist for this property to guide model development. It is therefore the aim of the work reported in this thesis to provide new experimental data that helps identify the turbulence sources located within the shear layer of jets. The technique of Partical Imaging Velocimetry (PIV) is used to capture directly the flowfield and all relevant turbulent statistics.