Die vorliegende Arbeit beschreibt die Entwicklung des hochauflösenden Atmosphärenmodells ASAMgpu. Dabei handelt es sich um ein sogenanntes Grobstrukturmodell bei dem gröbere Strukturen mit typischen Skalen von Deka- bis Kilometern in der atmosphärischen Grenzschicht explizit aufgelöst werden. Hochfrequentere Anteile und deren Dissipation müssen dabei entweder explizit mit einem Turbulenzmodell oder, wie im Falle des beschriebenen Modells, implizit behandelt werden. Dazu wurde der Advektionsoperator mit einem dissipativen Upwind-Verfahren dritter Ordnung diskretisiert. Das Modell beinhaltet ein Zwei-Momenten-Schema zur Beschreibung mikrophysikalischer Prozesse. Ein weiterer wichtiger Aspekt ist die verwendete thermodynamische Variable, die einige Vorteile herkömmlicher Ansätze vereint. Im Falle adiabatischer Prozesse stellt sie eine Erhaltungsgröße dar und die Quellen und Senken im Falle von Phasenumwandlungen sind leicht ableitbar. Außerdem können die benötigten Größen Temperatur und Druck explizit berechnet werden. Das gesamte Modell wurde in C++ implementiert und verwendet OpenGL und die OpenGL Shader Language (GLSL) um die nötigen Berechnungen auf Grafikkarten durchzuführen. Durch diesen Ansatz können genannte Simulationen, für die bisher Supercomputer nötig waren, sehr preisgünstig und energieeffizient durchgeführt werden. Neben der Modellbeschreibung werden die Ergebnisse einiger erfolgreicher Test-Simulationen, darunter drei Fälle mit mariner bewölkter Grenzschicht mit flacher Cumulusbewölkung, vorgestellt.

Fourier-based stochastic simulation of wind fields commonly used in wind turbine loads computations is unable to account for contrasting states of atmospheric stability. Flow fields in the stable boundary layer (SBL), for instance, have characteristics such as enhanced wind shear and veering wind direction profiles; the influence of such characteristics on utility-scale wind turbine loads has not been studied. To investigate these influences, we use large-eddy simulation (LES) to generate inflow wind fields and to estimate load statistics for a 5-MW wind turbine model. In the first part of this thesis, we describe a procedure employing LES to generate SBL wind fields for wind turbine load computations. In addition, we study how large-scale atmospheric conditions affect the characteristics of wind fields and turbine loads. Next, in the second part, we study the contrasting characteristics of LES-SBL and stochastic NBL flow fields and their influences on wind turbine load statistics by isolating effects of the mean wind (shear) profile and of variation in wind direction and turbulence levels over the rotor sept area. Among large-scale atmospheric conditions, the geostrophic wind speed and surface cooling rate have the greatest influence on flow field characteristics and, thus, on wind turbine loads. Higher geostrophic winds lead to increased mean and standard deviation values of the longitudinal wind speed at hub height. Increased surface cooling rates lead to steeper shear profiles and appear to also increase fatigue damage associated with out-of-plane blade root moments. In summary, our studies suggest that LES may be effectively used to model wind fields in the SBL, to study characteristics of turbine-scale wind fields, and to assess turbine loads for conditions that are not typically examined in design standards. / text

The present thesis deals with the development and performance analysis ofanisotropy-resolving models for the small, unresolved scales (”sub-grid scales”,SGS) in large eddy simulation (LES). The models are characterised by a descriptionof anisotropy by use of explicit algebraic models for both the subgridscale(SGS) stress tensor (EASSM) and SGS scalar flux vector (EASSFM). Extensiveanalysis of the performance of the explicit algebraic SGS stress model(EASSM) has been performed and comparisons made with the conventionalisotropic dynamic eddy viscosity model (DEVM). The studies include LES ofplane channel flow at relatively high Reynolds numbers and a wide range ofresolutions and LES of separated flow in a channel with streamwise periodichill-shaped constrictions (periodic hill flow) at coarse resolutions. The formersimulations were carried out with a pseudo-spectral Navier–Stokes solver, whilethe latter simulations were computed with a second-order, finite-volume basedsolver for unstructured grids. The LESs of channel flow demonstrate that theEASSM gives a good description of the SGS anisotropy, which in turn gives ahigh degree of resolution independence, contrary to the behaviour of LES predictionsusing the DEVM. LESs of periodic hill flow showed that the EASSMalso for this case gives significantly better flow predictions than the DEVM.In particular, the reattachment point was much better predicted with the EASSMand reasonably well predicted even at very coarse resolutions, where theDEVM is unable to predict a proper flow separation.The explicit algebraic SGS scalar flux model (EASSFM) is developed toimprove LES predictions of complex anisotropic flows with turbulent heat ormass transfer, and can be described as a nonlinear tensor eddy diffusivity model.It was tested in combination with the EASSM for the SGS stresses, and itsperformance was compared to the conventional dynamic eddy diffusivity model(DEDM) in channel flow with and without system rotation in the wall-normaldirection. EASSM and EASSFM gave predictions of high accuracy for meanvelocity and mean scalar fields, as well as stresses and scalar flux components.An extension of the EASSM and EASSFM, based on stochastic differentialequations of Langevin type, gave further improvements. In contrast to conventionalmodels, these extended models are able to describe intermittent transferof energy from the small, unresolved scales, to the resolved large ones.The present study shows that the EASSM/EASSFM gives a clear improvementof LES of wall-bounded flows in simple, as well as in complex geometriesin comparison with simpler SGS models. This is also shown to hold for a widerange of resolutions and is particularly accentuated for coarse resolution. The advantages are also demonstrated both for high-order numerical schemes andfor solvers using low-order finite volume methods. The models therefore havea clear potential for more applied computational fluid mechanics. / <p>QC 20140304</p> / Explicit algebraic sub-grid scale modelling for large-eddy simulations

Large-eddy simulations are carried out to study the combined effects of roughness and favourable pressure gradient in boundary layer flows, where the high acceleration (on smooth walls) may cause flow reversion to the quasi-laminar state. A sand-grain roughness model is used, with the no-slip boundary condition modeled by an immersed boundary method. The properties and accuracies of the scheme are studied, the roughness model is validated, and the spatial-resolution requirements are determined. The roughness model is applied to boundary layers subject to mild or strong acceleration, with simulations carried out underlining the effects of three parameters: the acceleration parameter, the roughness height, and the inlet Reynolds number. The roughness effects are limited to the roughness sublayer; the outer layer is affected indirectly only, through the changes that roughness causes in the relaminarization and retransition processes. The roughness significantly affects the inner-layer quantities like the friction velocity and the friction coefficient, while the local Reynolds number, the outer-layer mean velocity, as well as the Reynolds stresses beyond the roughness sublayer, are not sensitive to the roughness. The acceleration decreases the Reynolds stresses in the overlap region and promotes a laminar-like velocity profile. The acceleration leads to stabilization of near-wall structures and causes one-dimensional turbulence. The roughness generates small-scale structures at the bottom wall, which disturb the larger structures originally stabilized by the pressure gradient, leading to a decrease in the Reynolds-stress anisotropy. Roughness increases the Reynolds stresses in the roughness sublayer and tends to restore the fully turbulence flow early. The inlet Reynolds number affects the flow stability by determining the viscous length scale compared to the roughness length scales, and by determining how far the roughness effect extents into the boundary layer. / Thesis (Master, Mechanical and Materials Engineering) -- Queen's University, 2011-10-17 11:19:08.063

In the present doctoral work we have developed and tested a model for dissolved oxygen (DO) transfer from water to underlying flat and cohesive sediment beds populated with DO-absorbing bacteria. The model couples Large-Eddy Simulation (LES) of turbulent transport in the water-column, a biogeochemical model for DO transport and consumption in the sediment, and Darcy’s Law for the pore water-driven solute dispersion and advection. The model’s predictions compare well against experimental data for low friction-Reynolds numbers (Re). The disagreement for higher Re is investigated by progressively increasing the complexity of the model. A sensitivity analysis shows that the sediment-oxygen uptake (or demand, SOD) is approximately proportional to the bacterial content of the sediment layer, and varies with respect to fluid dynamics conditions, in accordance to classic high-Schmidt-number mass-transfer laws. The non- linear transport dynamics responsible for sustaining a statistically steady SOD are investigated by temporal- and-spatial correlations and with the aid of instantaneous visualizations: the near-wall coherent structures modulate the diffusive sublayer, which exhibits complex spatial and temporal filtering behaviours; its slow and quasi-periodic regeneration cycle determines the streaky structure of the DO field at the sediment-water interface (SWI), retained in the deeper layers of the porous medium. Oxygen depletion dynamics are then simulated by preventing surface re-areation with turbulent mixing driven by an oscillating low-speed current — an idealization of hypolimnetic DO depletion in the presence of a non-equilibrium periodic forcing. The oxygen distribution exhibits a self-similar pattern of decay with, during the deceleration phase, oscillations modulated by the periodic ejection of peaks of high turbulent mass flux (pumping oxygen towards the SWI), generated at the edge of the diffusive sublayer at the end of the acceleration phase. These fronts of highly turbulent mixing propagate away from the SWI, at approximately constant speed, in layers of below-average oxygen concentration. Finally, the model has been tested in a real geophysical framework, reproducing published in-situ DO measurements of a transitional flow in the bottom boundary layer of lake Alpnach. A simple model for the SOD is then derived for eventual inclusion in RANSE biogeochemical management-type models for similar applications. / Thesis (Ph.D, Mechanical and Materials Engineering) -- Queen's University, 2012-07-04 11:13:24.936

The research focusses on computational analysis of vortex generation behind A-Pillar of simplified model derived from Jaguar XF that excludes air from the underside of vehicle. This vortex formation contributes in generating wall pressure fluctuations especially at speeds higher than 100km/hr. It is a collaborative work between Cranfield University and Jaguar Land Rover. Three dimensional pressure based incompressible flow using Large Eddy Simulation turbulence model is selected for computational analysis in FLUENT v14. This used high parallel computing systems available in Cranfield University. In the initial phase, three grid resolutions (coarse, medium and fine) were prepared in ICEM CFD with fine case consisting of 10 million cells. Qualitative analysis includes extraction of slices, 3-D and surface streamlines and pressure and velocity contours for capturing the unsteadiness due to the vortex formation over the front side glass surface. The iso-surface of Q captures the unsteadiness at the A-Pillar wake and side mirror wake over front side glass surface. It also reveals that the range of length scales captured were limited even at the finest grid resolution. Quantitative analysis compares the mean pressure (Cp) data with JLR results. Probes were located at 51 locations over the front side glass window that showed a good comparison; specifically for the fine grid; with maximum variation incurred at probes located in separation areas. For predicting the wall pressure fluctuations, a total of ten probes were located over the front side glass window surface. The surface pressure (static) data was recorded for 1 sec of flow-time and later imported in MATLAB for post-processing. The results obtained in 1/3rd octave band showed that the large scales were too energetic and small scales are not captured. However, comparing sound pressure levels with the Aero-acoustic Wind Tunnel (AWT); provided by JLR; it is concluded that either the grid is too coarse to resolve higher frequencies or the numerical modelling used is too dissipative to benefits the use of LES.

In this thesis, wall-resolved and wall-modeled large-eddy simulation (LES) have been employed to investigate turbulent flow and dispersion around a single and a group of wall-mounted bluff bodies which are partially and fully submerged in developing boundary layers, respectively. The dispersion is caused by a continuous release of a passive scalar from a ground-level point source located within the matrix of obstacles. The results have been validated through comparisons against the available experimental measurement data. Thorough physical analysis including investigation of the spatial evolution and temporal cascades of the kinetic and scalar energies, flow structures and their influences on dispersion of the concentration plume in the context of highly disturbed flows, and study of turbulence statistics for the flow and concentration fields have been performed to provide deeper insights into turbulent flow and dispersion in domains with complex geometries. An in-house code based on FORTRAN programming language, parallelized with MPI libraries has been developed, modified and optimized for conducting the simulations. The simulations have been conducted on public-domain supercomputers ofWest-Grid, specifically Orcinus and Grex, and also the local 256-core cluster system of the CFD LAB at the University of Manitoba.

Since its introduction in 1997, the use of Detached Eddy Simulation (DES) and similar hybrid turbulence techniques has become increasingly popular in the field of CFD. However, with increased use some of the limitations of the DES model have become apparent. One of these is the dependence of DES on grid construction, particularly regarding the point of transition between the Reynolds-Averaged Navier-Stokes and Large Eddy Simulation models. An additional issue that arises with unstructured grids is the definition of the grid spacing in the implementation of a DES length scale. To lay the ground work to study these effects the Spalart-Allmaras one-equation turbulence model, SA based DES hybrid turbulence model, and the Scale Adaptive Simulation hybrid turbulence model are implemented in an unstructured grid CFD code, UNCLE. The implemented SA based DES model is validated for flow over a three-dimensional circular cylinder for three different turbulent Reynolds numbers. Validation included studying the pressure, skin friction coefficient, centerline velocity distributions averaged in time and space. Tools to output the mean velocity profiles and Reynolds stresses were developed. A grid generation code was written to generate a two/three dimensional circular cylinder grid to simulate flow over the cylinder in UNCLE. The models implemented and validated, and the additional tools mentioned will be used in the future.

The normal injection of a liquid jet into a gaseous crossflow has many engineering applications. In this thesis, detailed numerical simulations based on the level set method are employed to understand the physical mechanism underlying the jet ``surface breakup''. The numerical observations reveal the existence of hydrodynamic instabilities on the jet periphery. The temporal growth of such azimuthal instabilities leads to the formation of interface corrugations, which are eventually sheared off of the jet surface as sheet-like structures. The sheets finally undergo disintegration into ligaments and drops during the surface breakup process. Temporal linear stability analyses are employed to understand the nature of these instabilities. To facilitate the analysis, analytical solutions for the flow fields of the jet and the crossflow are derived. We identify the ``shear instability'' as the primary destabilization mechanism in the flow. This inherently inviscid mechanism opposes the previously suggested mechanism of surface breakup (known as ``boundary layer stripping''), which is based on a viscous interpretation. The influence of the jet-to-crossflow density ratio on the flow stability are also studied. The findings show that a higher density jet leads to higher wavenumber instabilities on the jet surface and thereby subsequent smaller drops and ligaments. The stability characteristics of the most amplified modes (i.e., the wavenumber and corresponding growth rate) obtained from stability analyses and numerical simulations are in good agreement. The stability results of the jet also show that the density may have a non-monotonic stabilizing/destabilizing effect on the flow stability. To investigate such effect, the concept of wave resonance are employed to physically interpret the inviscid instability mechanism in two-phase flows with sharp interfaces and linear velocity profiles. We demonstrate that neutrally stable waves are formed due to the density jump in the flow, in addition to the well-known vorticity (Rayleigh) waves. Under certain conditions, such neutral waves are capable of resonating and generating unstable modes. The resonance of different pairs of neutral waves, therefore, results in either stabilizing or destabilizing effect of density variation. We predict similar reasoning behind the density behavior in the jet in crossflow configuration with smoothly varying velocity and density profiles.

A sub-grid scale closure for Large Eddy Simulation (LES) of turbulent combustion, based on physical space filtering of laminar flames is presented. The proposed formalism relies on a presumed probability density function (PDF) derived from the filtered laminar flames and flamelet tabulated chemistry. The combustion LES filter size is not fixed in this novel approach when sub-grid scale wrinkling occurs, but calibrated depending on the local level of unresolved scalar fluctuations. The model was validated by simulating 1D filtered laminar flames and 2D Bunsen flames. Subsequently, the model was tested on a 3D turbulent scenario by performing LES of the premixed and stratified configurations of the Cambridge swirl burner, experimentally studied by Sweeney and co-workers. Comparison of simulation and experiments for both the premixed and stratified configurations showed good agreement emphasizing the model characteristiscs. Instantaneous and time averaged LES data were analyzed to extract