CFD simulation of an unbaffled stirred tank reactor driven by a magnetic rod: assessment of turbulence models.

PubMed

Li, Jiajia; Deng, Baoqing; Zhang, Bing; Shen, Xiuzhong; Kim, Chang Nyung

2015-01-01

A simulation of an unbaffled stirred tank reactor driven by a magnetic stirring rod was carried out in a moving reference frame. The free surface of unbaffled stirred tank was captured by Euler-Euler model coupled with the volume of fluid (VOF) method. The re-normalization group (RNG) k-ɛ model, large eddy simulation (LES) model and detached eddy simulation (DES) model were evaluated for simulating the flow field in the stirred tank. All turbulence models can reproduce the tangential velocity in an unbaffled stirred tank with a rotational speed of 150 rpm, 250 rpm and 400 rpm, respectively. Radial velocity is underpredicted by the three models. LES model and RNG k-ɛ model predict the better tangential velocity and axial velocity, respectively. RNG k-ɛ model is recommended for the simulation of the flow in an unbaffled stirred tank with magnetic rod due to its computational effort.

Numerical investigation of a helicopter combustion chamber using LES and tabulated chemistry

NASA Astrophysics Data System (ADS)

Auzillon, Pierre; Riber, Eléonore; Gicquel, Laurent Y. M.; Gicquel, Olivier; Darabiha, Nasser; Veynante, Denis; Fiorina, Benoît

2013-01-01

This article presents Large Eddy Simulations (LES) of a realistic aeronautical combustor device: the chamber CTA1 designed by TURBOMECA. Under nominal operating conditions, experiments show hot spots observed on the combustor walls, in the vicinity of the injectors. These high temperature regions disappear when modifying the fuel stream equivalence ratio. In order to account for detailed chemistry effects within LES, the numerical simulation uses the recently developed turbulent combustion model F-TACLES (Filtered TAbulated Chemistry for LES). The principle of this model is first to generate a lookup table where thermochemical variables are computed from a set of filtered laminar unstrained premixed flamelets. To model the interactions between the flame and the turbulence at the subgrid scale, a flame wrinkling analytical model is introduced and the Filtered Density Function (FDF) of the mixture fraction is modeled by a β function. Filtered thermochemical quantities are stored as a function of three coordinates: the filtered progress variable, the filtered mixture fraction and the mixture fraction subgrid scale variance. The chemical lookup table is then coupled with the LES using a mathematical formalism that ensures an accurate prediction of the flame dynamics. The numerical simulation of the CTA1 chamber with the F-TACLES turbulent combustion model reproduces fairly the temperature fields observed in experiments. In particular the influence of the fuel stream equivalence ratio on the flame position is well captured.

LES of Swirling Reacting Flows via the Unstructured scalar-FDF Solver

NASA Astrophysics Data System (ADS)

Ansari, Naseem; Pisciuneri, Patrick; Strakey, Peter; Givi, Peyman

2011-11-01

Swirling flames pose a significant challenge for computational modeling due to the presence of recirculation regions and vortex shedding. In this work, results are presented of LES of two swirl stabilized non-premixed flames (SM1 and SM2) via the FDF methodology. These flames are part of the database for validation of turbulent-combustion models. The scalar-FDF is simulated on a domain discretized by unstructured meshes, and is coupled with a finite volume flow solver. In the SM1 flame (with a low swirl number) chemistry is described by the flamelet model based on the full GRI 2.11 mechanism. The SM2 flame (with a high swirl number) is simulated via a 46-step 17-species mechanism. The simulated results are assessed via comparison with experimental data.

Grid Effects on LES Thermo-Acoustic Limit-Cycle of a Full Annular Aeronautical Engine

NASA Astrophysics Data System (ADS)

Wolf, Pierre; Gicquel, Laurent Y. M.; Staffelbach, Gabriel; Poinsot, Thierry

Recent developments in large scale computer architectures allow Large Eddy Simulation (LES) to be considered for the prediction of turbulent reacting flows in geometries encountered in industry. To do so, various difficulties must be overcome and the first one is to ensure that proper meshes can be used for LES. Indeed, the quality of meshes is known to be a critical factor in LES of reacting flows. This issue becomes even more crucial when LES is used to compute large configurations such as full annular combustion chambers. Various analysis of mesh effects on LES results have been published before but all are limited to single-sector computational domains. However, real annular gas-turbine engines contain ten to twenty of such sectors and LES must also be used in such full chambers for the study of ignition or azimuthal thermo-acoustic interactions. Instabilities (mostly azimuthal modes involving the full annular geometry) remain a critical issue to aeronautical or power-generation industries and LES seems to be a promising path to properly apprehend such complex unsteady couplings. Based on these observations, mesh effects on LES in a full annular gas-turbine combustion chamber (including its casing) is studied here in the context of its azimuthal thermo-acoustic response. To do so, a fully compressible, multi-species reacting LES is used on two meshes yielding two fully unsteady turbulent reacting predictions of the same configuration. The two tetrahedra meshes contain respectively 38 and 93 millions cells. Limit-cycles as obtained by the two LES are gauged against each other for various flow quantities such as mean velocity profiles, flame position and temperature fields. The thermo-acoustic limit-cycles are observed to be relatively indepedent of the grid resolution which comforts the use of LES tools to provide insights and understanding of the mechanisms triggering the coupling between the system acoustic eigenmodes and combustion.

A fully coupled bolus-esophageal-gastric model for esophageal emptying based on the immersed boundary method

NASA Astrophysics Data System (ADS)

Kou, Wenjun; Pandolfino, John E.; Kahrilas, Peter J.; Patankar, Neelesh A.

2016-11-01

In this work, we develop a fully coupled bolus-esophageal-gastric model to study esophageal emptying based on the immersed boundary method. The model includes an esophageal segment, an ellipsoid-shaped stomach, and a bolus. It can easily handle the passive and active function of the lower esophageal sphincter (LES). Two groups of case studies are presented. The first group is about the influence from tissue anisotropy. Simulation shows that the weaker (or more compliant) part suffers from a higher wall shear stress and higher pressure load when the bolus is filled in and emptied from the LES segment. This implies a degradation cycle in which a weaker tissue becomes much weaker due to an increased load, a possible pathway to the esophageal lower diverticulum. The second group is about bulge formation resulting from asymmetric anatomy and a compliant LES. In particular, we find a right bulge tends to develop for a compliant LES. The bulge is most pronounced with a highest stiffness of the gastric wall. This implies that the competition between the LES stiffness and gastric wall stiffness might be another factor related to the esophageal lower diverticulum. The support of Grant R01 DK56033 and R01 DK079902 from NIH is gratefully acknowledged.

Large eddy simulation of reactive pollutants in a deep urban street canyon: Coupling dynamics with O3-NOx-VOC chemistry.

PubMed

Zhong, Jian; Cai, Xiao-Ming; Bloss, William James

2017-05-01

A large eddy simulation (LES) model coupled with O 3 -NO x -VOC chemistry is implemented to simulate the coupled effects of emissions, mixing and chemical pre-processing within an idealised deep (aspect ratio = 2) urban street canyon under a weak wind condition. Reactive pollutants exhibit significant spatial variations in the presence of two vertically aligned unsteady vortices formed in the canyon. Comparison of the LES results from two chemical schemes (simple NO x -O 3 chemistry and a more comprehensive Reduced Chemical Scheme (RCS) chemical mechanism) shows that the concentrations of NO 2 and O x inside the street canyon are enhanced by approximately 30-40% via OH/HO 2 chemistry. NO, NO x , O 3 , OH and HO 2 are chemically consumed, while NO 2 and O x (total oxidant) are chemically produced within the canyon environment. Within-canyon pre-processing increases oxidant fluxes from the canyon to the overlying boundary layer, and this effect is greater for deeper street canyons (as found in many traditional European urban centres) than shallower (lower aspect ratio) streets. There is clear evidence of distinct behaviours for emitted chemical species and entrained chemical species, and positive (or negative) values of intensities of segregations are found between pairs of species with similar (or opposite) behaviour. The simplified two-box model underestimated NO and O 3 levels, but overestimated NO 2 levels for both the lower and upper canyon compared with the more realistic LES-chemistry model. This suggests that the segregation effect due to incomplete mixing reduces the chemical conversion rate of NO to NO 2 . This study reveals the impacts of nonlinear O 3 -NO x -VOC photochemical processes in the incomplete mixing environment and provides a better understanding of the pre-processing of emissions within canyons, prior to their release to the urban boundary layer, through the coupling of street canyon dynamics and chemistry. Copyright © 2017 Elsevier Ltd. All rights reserved.

Self-similarity and flow characteristics of vertical-axis wind turbine wakes: an LES study

NASA Astrophysics Data System (ADS)

Abkar, Mahdi; Dabiri, John O.

2017-04-01

Large eddy simulation (LES) is coupled with a turbine model to study the structure of the wake behind a vertical-axis wind turbine (VAWT). In the simulations, a tuning-free anisotropic minimum dissipation model is used to parameterise the subfilter stress tensor, while the turbine-induced forces are modelled with an actuator line technique. The LES framework is first validated in the simulation of the wake behind a model straight-bladed VAWT placed in the water channel and then used to study the wake structure downwind of a full-scale VAWT sited in the atmospheric boundary layer. In particular, the self-similarity of the wake is examined, and it is found that the wake velocity deficit can be well characterised by a two-dimensional multivariate Gaussian distribution. By assuming a self-similar Gaussian distribution of the velocity deficit, and applying mass and momentum conservation, an analytical model is developed and tested to predict the maximum velocity deficit downwind of the turbine. Also, a simple parameterisation of VAWTs for LES with very coarse grid resolutions is proposed, in which the turbine is modelled as a rectangular porous plate with the same thrust coefficient. The simulation results show that, after some downwind distance (x/D ≈ 6), both actuator line and rectangular porous plate models have similar predictions for the mean velocity deficit. These results are of particular importance in simulations of large wind farms where, due to the coarse spatial resolution, the flow around individual VAWTs is not resolved.

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

NASA Technical Reports Server (NTRS)

Drozda, Tomasz G.; Quinlan, Jesse R.; Pisciuneri, Patrick H.; Yilmaz, S. Levent

2012-01-01

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

CFD study of mixing miscible liquid with high viscosity difference in a stirred tank

NASA Astrophysics Data System (ADS)

Madhania, S.; Cahyani, A. B.; Nurtono, T.; Muharam, Y.; Winardi, S.; Purwanto, W. W.

2018-03-01

The mixing process of miscible liquids with high viscosity difference is crucial role even though the solution mutually dissolved. This paper describes the mixing behaviour of the water-molasses system in a conical-bottomed cylindrical stirred tank (D = 0.28 m and H = 0.395 m) equipped with a side-entry Marine propeller (d = 0.036 m) under the turbulence regime using a three-dimensional and transient CFD-simulation. The objective of this work is to compare the solution strategies was applied in the computational analysis to capture the detail phenomena of mixing two miscible liquid with high viscosity difference. Four solution strategies that have been used are the RANS Standards k-ε (SKE) model as the turbulence model coupled with the Multiple Reference Frame (MRF) method for impeller motion, the RANS Realizable k-ε (RKE) combine with the MRF, the Large Eddy Simulation (LES) coupled with the Sliding Mesh (SM) method and the LES-MRF combination. The transient calculations were conducted with Ansys Fluent 17.1 version. The mixing behaviour and the propeller characteristic are to be compared and discussed in this work. The simulation results show the differences of flow pattern and the molasses distribution profile for every solution strategy. The variation of the flow pattern which happened in each solution strategy showing an instability of the mixing process in stirred tank. The LES-SM strategy shows the realistic direction of flow than another solution strategies.

Computation of unsteady turbomachinery flows: Part 2—LES and hybrids

NASA Astrophysics Data System (ADS)

Tucker, P. G.

2011-10-01

The choice of turbulence model can have a strong impact on results for many turbomachinery zones. Palliative corrections to them and also transition modeling can have a further profound solution impact. The spectral gaps necessary for theoretically valid URANS solutions are also lacking in certain turbomachinery zones. Large Eddy Simulation (LES) alleviates the serious area of turbulence modeling uncertainty but with an extreme increase in computational cost. However, there seems a lack of validation data to explore in depth the performance of LES and thus strategies to refine it. LES best practices are needed. Although LES is, obviously, much less model dependent than RANS, grids currently used for more practical simulations are clearly insufficiently fine for the LES model and numerical schemes not to be playing an excessively strong role. Very few turbomachinery simulations make use of properly constructed, correlated turbulence inflow. Even if this is attempted, most measurement sets are incomplete and lack an adequate basis for modeling this inflow. Gas turbines are highly complex coupled systems and hence inflow and outflow boundary condition specification needs to go beyond just synthesizing turbulent structures and preventing their reflection. Despite the strong limitations of the dissipative Smagorinsky model, it still sees the most wide spread use, generally, in excessively dissipative flow solvers. Monotone Integrated LES (MILES) related approaches, hybrid LES-RANS and more advanced LES models seem to have an equal but subservient frequency of use in turbomachinery applications. Clearly the introduction of a RANS layer can have a substantial accuracy penalty. However, it does allow LES to be rationally used, albeit in a diluted sense for industrial applications. The Reynolds numbers found in turbomachinery are substantial. However, in certain areas evidence suggests they will not be enough to ensure a long inertial subrange and hence the use of standard LES modeling practices. Despite the excessively coarse grids used in much of the LES work reviewed, with essentially RANS based codes, meaningful results are often gained. This can perhaps be attributed to the choice of cases, these being ones for which RANS modeling gives extremely poor performance. It is a concern that for practical turbomachinery LES studies grid densities used tend to have an Reynolds number scaling to a strong negative power.

Large Eddy Simulation of Spatially Developing Turbulent Reacting Shear Layers with the One-Dimensional Turbulence Model

NASA Astrophysics Data System (ADS)

Hoffie, Andreas Frank

Large eddy simulation (LES) combined with the one-dimensional turbulence (ODT) model is used to simulate spatially developing turbulent reacting shear layers with high heat release and high Reynolds numbers. The LES-ODT results are compared to results from direct numerical simulations (DNS), for model development and validation purposes. The LES-ODT approach is based on LES solutions for momentum and pressure on a coarse grid and solutions for momentum and reactive scalars on a fine, one-dimensional, but three-dimensionally coupled ODT subgrid, which is embedded into the LES computational domain. Although one-dimensional, all three velocity components are transported along the ODT domain. The low-dimensional spatial and temporal resolution of the subgrid scales describe a new modeling paradigm, referred to as autonomous microstructure evolution (AME) models, which resolve the multiscale nature of turbulence down to the Kolmogorv scales. While this new concept aims to mimic the turbulent cascade and to reduce the number of input parameters, AME enables also regime-independent combustion modeling, capable to simulate multiphysics problems simultaneously. The LES as well as the one-dimensional transport equations are solved using an incompressible, low Mach number approximation, however the effects of heat release are accounted for through variable density computed by the ideal gas equation of state, based on temperature variations. The computations are carried out on a three-dimensional structured mesh, which is stretched in the transverse direction. While the LES momentum equation is integrated with a third-order Runge-Kutta time-integration, the time integration at the ODT level is accomplished with an explicit Forward-Euler method. Spatial finite-difference schemes of third (LES) and first (ODT) order are utilized and a fully consistent fractional-step method at the LES level is used. Turbulence closure at the LES level is achieved by utilizing the Smagorinsky model. The chemical reaction is simulated with a global single-step, second-order equilibrium reaction with an Arrhenius reaction rate. The two benchmark cases of constant density reacting and variable density non-reacting shear layers used to determine ODT parameters yield perfect agreement with regards to first and second-order flow statistics as well as shear layer growth rate. The variable density non-reacting shear layer also serves as a testing case for the LES-ODT model to simulate passive scalar mixing. The variable density, reacting shear layer cases only agree reasonably well and indicate that more work is necessary to improve variable density coupling of ODT and LES. The disagreement is attributed to the fact that the ODT filtered density is kept constant across the Runge-Kutta steps. Furthermore, a more in-depth knowledge of large scale and subgrid turbulent kinetic energy (TKE) spectra at several downstream locations as well as TKE budgets need to be studied to obtain a better understanding about the model as well as about the flow under investigation. The local Reynolds number based on the one-percent thickness at the exit is Redelta ≈ 5300, for the constant density reacting and for the variable density non-reacting case. For the variable density reacting shear layer, the Reynolds number based on the 1% thickness is Redelta ≈ 2370. The variable density reacting shear layers show suppressed growth rates due to density variations caused by heat release. This has also been reported in literature. A Lewis number parameter study is performed to extract non-unity Lewis number effects. An increase in the Lewis number leads to a further suppression of the growth rate, however to an increase spread of second-order flow statistics. Major focus and challenge of this work is to improve and advance the three-dimensional coupling of the one-dimensional ODT domains while keeping the solution correct. This entails major restructuring of the model. The turbulent reacting shear layer poses a physical challenge to the model because of its nature being a statistically stationary, non-decaying inhomogeneous and anisotropic turbulent flow. This challenge also requires additions to the eddy sampling procedure. Besides physical advancements, the LES-ODT code is also improved regarding its ability to use general cuboid geometries, an array structure that allows to apply boundary conditions based on ghost-cells and non-uniform structured meshes. The use of transverse grid-stretching requires the implementation of the ODT triplet map on a stretched grid. Further, advancing subroutine structure handling with global variables that enable serial code speed-up and parallelization with OpenMP are undertaken. Porting the code to a higher-level language, object oriented, finite-volume based CFD platform, like OpenFoam for example that allows more advanced array and parallelization features with graphics processing units (GPUs) as well as parallelization with the message passing interface (MPI) to simulate complex geometries is recommended for future work.

A Study of Strain Rate Effects for Turbulent Premixed Flames with Application to LES of a Gas Turbine Combustor Model

DOE PAGES

Kemenov, Konstantin A.; Calhoon, William H.

2015-03-24

Large-scale strain rate field, a resolved quantity which is easily computable in large-eddy simulations (LES), could have profound effects on the premixed flame properties by altering the turbulent flame speed and inducing local extinction. The role of the resolved strain rate has been investigated in a posterior LES study of GE lean premixed dry low NOx emissions LM6000 gas turbine combustor model. A novel approach which is based on the coupling of the lineareddy model with a one-dimensional counter-flow solver has been applied to obtain the parameterizations of the resolved premixed flame properties in terms of the reactive progress variable,more » the local strain rate measure, and local Reynolds and Karlovitz numbers. The strain rate effects have been analyzed by comparing LES statistics for several models of the turbulent flame speed, i.e, with and without accounting for the local strain rate effects, with available experimental data. The sensitivity of the simulation results to the inflow velocity conditions as well as the grid resolution have been also studied. Overall, the results suggest the necessity to represent the strain rate effects accurately in order to improve LES modeling of the turbulent flame speed.« less

Des schémas équivalents pour les circuits couplés multi-enroulements

NASA Astrophysics Data System (ADS)

Keradec, J. P.; Cogitore, B.; Laveuve, E.; Bensoam, M.

1994-04-01

The aim of this paper is to represent the electrical behaviour of any number of magnetically coupled windings with couplers and inductors. Two methods, mathematicaly justified, are proposed. The second one introduces only positive inductances. As an exemple, it is applied to the representation of a three column three phase transformer. The obtained circuits supply the requisite guide to design more complete circuits which allow the high frequency behaviour of wound components to be taken into account, especialy in electronics simulation softwares. Le but de cet article est de traduire le comportement électrique d'un nombre quelconque d'enroulements magnétiquement couplés, par des coupleurs et des inductances. Deux méthodes, établies mathématiquement, sont proposées. La seconde n'introduit que des inductances positives. A titre d'exemple, elle est appliquée à la représentation d'un transformateur triphasé à trois colonnes. Les schémas obtenus fournissent l'indispensable ossature de schémas plus complets, aptes à représenter le comportement haute fréquence des composants bobinés, notamment dans un logiciel de simulation électronique.

High speed turbulent reacting flows: DNS and LES

NASA Technical Reports Server (NTRS)

Givi, Peyman

1990-01-01

Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.

Dust devil characteristics and associated dust entrainment based on large-eddy simulations

NASA Astrophysics Data System (ADS)

Klose, Martina; Kwidzinski, Nick; Shao, Yaping

2015-04-01

The characteristics of dust devils, such as occurrence frequency, lifetime, size, and intensity, are usually inferred from in situ field measurements and remote sensing. Numerical models, e.g. large-eddy simulation (LES) models, have also been established as a tool to investigate dust devils and their structures. However, most LES models do not contain a dust module. Here, we present results from simulations using the WRF-LES model coupled to the convective turbulent dust emission (CTDE) scheme of Klose et al. (2014). The scheme describes the stochastic process of aerodynamic dust entrainment in the absence of saltation. It therefore allows for dust emission even below the threshold friction velocity for saltation. Numerical experiments have been conducted for different atmospheric stability and background wind conditions at 10 m horizontal resolution. A dust devil tracking algorithm is used to identify dust devils in the simulation results. The detected dust devils are statistically analyzed with regard to e.g. radius, pressure drop, lifetime, and turbulent wind speeds. An additional simulation with higher horizontal resolution (2 m) is conducted for conditions, which are especially favorable for dust devil development, i.e. unstable atmospheric stratification and weak mean winds. The higher resolution enables the identification of smaller dust devils and a more detailed structure analysis. Dust emission fluxes, dust concentrations, and dust mass budgets are calculated from the simulations. The results are compared to field observations reported in literature.

Dynamic Turbulence Modelling in Large-eddy Simulations of the Cloud-topped Atmospheric Boundary Layer

NASA Technical Reports Server (NTRS)

Kirkpatrick, M. P.; Mansour, N. N.; Ackerman, A. S.; Stevens, D. E.

2003-01-01

The use of large eddy simulation, or LES, to study the atmospheric boundary layer dates back to the early 1970s when Deardor (1972) used a three-dimensional simulation to determine velocity and temperature scales in the convective boundary layer. In 1974 he applied LES to the problem of mixing layer entrainment (Deardor 1974) and in 1980 to the cloud-topped boundary layer (Deardor 1980b). Since that time the LES approach has been applied to atmospheric boundary layer problems by numerous authors. While LES has been shown to be relatively robust for simple cases such as a clear, convective boundary layer (Mason 1989), simulation of the cloud-topped boundary layer has proved more of a challenge. The combination of small length scales and anisotropic turbulence coupled with cloud microphysics and radiation effects places a heavy burden on the turbulence model, especially in the cloud-top region. Consequently, over the past few decades considerable effort has been devoted to developing turbulence models that are better able to parameterize these processes. Much of this work has involved taking parameterizations developed for neutral boundary layers and deriving corrections to account for buoyancy effects associated with the background stratification and local buoyancy sources due to radiative and latent heat transfer within the cloud (see Lilly 1962; Deardor 1980a; Mason 1989; MacVean & Mason 1990, for example). In this paper we hope to contribute to this effort by presenting a number of turbulence models in which the model coefficients are calculated dynamically during the simulation rather than being prescribed a priori.

Development of a sharp interface model that simulates coastal aquifer flow with the coupled use of GIS

NASA Astrophysics Data System (ADS)

Gemitzi, Alexandra; Tolikas, Demetrios

A simulation program, which works seamlessly with GIS and simulates flows in coastal aquifers, is presented in the present paper. The model is based on the Galerkin finite element discretization scheme and it simulates both steady and transient freshwater and saltwater flow, assuming that the two fluids are separated by a sharp interface. The model has been verified in simple cases where analytical solutions exist. The simulation program works as a tool of the GIS program, which is the main database that stores and manages all the necessary data. The combined use of the simulation and the GIS program forms an integrated management tool offering a simpler way of simulating and studying saline intrusion in coastal aquifers. Application of the model to the Yermasogia aquifer illustrates the coupled use of modeling and GIS techniques for the examination of regional coastal aquifer systems. Pour étudier un système aquifère côtier, nous avons développé un modèle aux éléments finis en quasi 3-D qui simule les écoulements d'eau douce et d'eau salée en régime aussi bien permanent que transitoire. Les équations qui les régissent sont discrétisées par un schéma de discrétisation de Garlekin aux éléments finis. Le modèle a été vérifié dans des cas simples où il existe des solutions analytiques. Toutes les données nécessaires sont introduites et gérées grâce à un logiciel de gestion de SIG. Le programme de simulation est utilisé comme un outil du logiciel de SIG, constituant ainsi un outil de gestion intégrée dont le but est de simuler et d'étudier l'intrusion saline dans les aquifères côtiers. L'application du modèle à l'aquifère de Yermasogia illustre l'utilisation couplée de la modélisation et des techniques de SIG pour l'étude des systèmes aquifères côtiers régionaux. Se ha desarrollado un modelo casi tridimensional de elementos finitos para simular el flujo de agua dulce y salada, tanto en régimen estacionario como en transitorio, en sistemas acuíferos costeros, bajo la hipótesis de separación por medio de una interfaz abrupta. Las ecuaciones del modelo han sido discretizadas mediante un esquema de Galerkin de discretización en elementos finitos. El modelo ha sido verificado en casos sencillos para los que existe solución analítica. Todos los datos necesarios se introducen y gestionan con un Sistema de Información Geográfica [SIG] por ordenador. El programa de simulación forma parte del programa de SIG, constituyendo una herramienta integrada de gestión para estudiar la intrusión salina en acuíferos costeros. La aplicación del modelo al acuífero de Yermasogia ilustra el uso acoplado de las técnicas de modelación y de SIG con el fin de examinar sistemas acuíferos costeros a escala regional.

LES FOR SIMULATING THE GAS EXCHANGE PROCESS IN A SPARK IGNITION ENGINE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ameen, Muhsin M; yang, xiaofeng; kuo, tang-wei

2015-01-01

The gas exchange process is known to be a significant source of cyclic variability in Internal Combustion Engines (ICE). Traditionally, Large Eddy Simulations (LES) are expected to capture these cycle-to-cycle variations. This paper reports a numerical effort to establish best practices for capturing cyclic variability with LES tools in a Transparent Combustion Chamber (TCC) spark ignition engine. The main intention is to examine the sensitivity of cycle averaged mean and Root Mean Square (RMS) flow fields and Proper Orthogonal Decomposition (POD) modes to different computational hardware, adaptive mesh refinement (AMR) and LES sub-grid scale (SGS) models, since these aspects havemore » received little attention in the past couple of decades. This study also examines the effect of near-wall resolution on the predicted wall shear stresses. LES is pursued with commercially available CONVERGE code. Two different SGS models are tested, a one-equation eddy viscosity model and dynamic structure model. The results seem to indicate that both mean and RMS fields without any SGS model are not much different than those with LES models, either one-equation eddy viscosity or dynamic structure model. Computational hardware results in subtle quantitative differences, especially in RMS distributions. The influence of AMR on both mean and RMS fields is negligible. The predicted shear stresses near the liner walls is also found to be relatively insensitive to near-wall resolution except in the valve curtain region.« less

A Large-eddy Simulation Study of Vertical Axis Wind Turbine Wakes in the Atmospheric Boundary Layer

NASA Astrophysics Data System (ADS)

Shamsoddin, Sina; Porté-Agel, Fernando

2016-04-01

Vertical axis wind turbines (VAWTs) offer some advantages over their horizontal axis counterparts, and are being considered as a viable alternative to conventional horizontal axis wind turbines (HAWTs). Nevertheless, a relative shortage of scientific, academic and technical investigations of VAWTs is observed in the wind energy community with respect to HAWTs. Having this in mind, in this work, we aim to study the wake of a single VAWT, placed in the atmospheric boundary layer, using large-eddy simulation (LES) coupled with actuator line model (ALM). It is noteworthy that this is the first time that such a study is being performed. To do this, for a typical 1 MW VAWT design, first, the variation of power coefficient with both the chord length of the blades and the tip-speed ratio is analyzed using LES-ALM, and an optimum combination of chord length and tip-speed ratio is obtained. Subsequently, the wake of a VAWT with these optimum specifications is thoroughly examined by showing different relevant mean and turbulent wake flow statistics. Keywords: vertical axis wind turbine (VAWT); VAWT wake; Atmospheric Boundary Layer (ABL); large eddy simulation (LES); actuator line model (ALM); turbulence.

On the influences of key modelling constants of large eddy simulations for large-scale compartment fires predictions

NASA Astrophysics Data System (ADS)

Yuen, Anthony C. Y.; Yeoh, Guan H.; Timchenko, Victoria; Cheung, Sherman C. P.; Chan, Qing N.; Chen, Timothy

2017-09-01

An in-house large eddy simulation (LES) based fire field model has been developed for large-scale compartment fire simulations. The model incorporates four major components, including subgrid-scale turbulence, combustion, soot and radiation models which are fully coupled. It is designed to simulate the temporal and fluid dynamical effects of turbulent reaction flow for non-premixed diffusion flame. Parametric studies were performed based on a large-scale fire experiment carried out in a 39-m long test hall facility. Several turbulent Prandtl and Schmidt numbers ranging from 0.2 to 0.5, and Smagorinsky constants ranging from 0.18 to 0.23 were investigated. It was found that the temperature and flow field predictions were most accurate with turbulent Prandtl and Schmidt numbers of 0.3, respectively, and a Smagorinsky constant of 0.2 applied. In addition, by utilising a set of numerically verified key modelling parameters, the smoke filling process was successfully captured by the present LES model.

An interior penalty stabilised incompressible discontinuous Galerkin-Fourier solver for implicit large eddy simulations

NASA Astrophysics Data System (ADS)

Ferrer, Esteban

2017-11-01

We present an implicit Large Eddy Simulation (iLES) h / p high order (≥2) unstructured Discontinuous Galerkin-Fourier solver with sliding meshes. The solver extends the laminar version of Ferrer and Willden, 2012 [34], to enable the simulation of turbulent flows at moderately high Reynolds numbers in the incompressible regime. This solver allows accurate flow solutions of the laminar and turbulent 3D incompressible Navier-Stokes equations on moving and static regions coupled through a high order sliding interface. The spatial discretisation is provided by the Symmetric Interior Penalty Discontinuous Galerkin (IP-DG) method in the x-y plane coupled with a purely spectral method that uses Fourier series and allows efficient computation of spanwise periodic three-dimensional flows. Since high order methods (e.g. discontinuous Galerkin and Fourier) are unable to provide enough numerical dissipation to enable under-resolved high Reynolds computations (i.e. as necessary in the iLES approach), we adapt the laminar version of the solver to increase (controllably) the dissipation and enhance the stability in under-resolved simulations. The novel stabilisation relies on increasing the penalty parameter included in the DG interior penalty (IP) formulation. The latter penalty term is included when discretising the linear viscous terms in the incompressible Navier-Stokes equations. These viscous penalty fluxes substitute the stabilising effect of non-linear fluxes, which has been the main trend in implicit LES discontinuous Galerkin approaches. The IP-DG penalty term provides energy dissipation, which is controlled by the numerical jumps at element interfaces (e.g. large in under-resolved regions) such as to stabilise under-resolved high Reynolds number flows. This dissipative term has minimal impact in well resolved regions and its implicit treatment does not restrict the use of large time steps, thus providing an efficient stabilization mechanism for iLES. The IP-DG stabilisation is complemented with a Spectral Vanishing Viscosity (SVV) method, in the z-direction, to enhance stability in the continuous Fourier space. The coupling between the numerical viscosity in the DG plane and the SVV damping, provides an efficient approach to stabilise high order methods at moderately high Reynolds numbers. We validate the formulation for three turbulent flow cases: a circular cylinder at Re = 3900, a static and pitch oscillating NACA 0012 airfoil at Re = 10000 and finally a rotating vertical-axis turbine at Re = 40000, with Reynolds based on the circular diameter, airfoil chord and turbine diameter, respectively. All our results compare favourably with published direct numerical simulations, large eddy simulations or experimental data. We conclude that the DG-Fourier high order solver, with IP-SVV stabilisation, proves to be a valuable tool to predict turbulent flows and associated statistics for both static and rotating machinery.

Hybrid LES RANS technique based on a one-equation near-wall model

NASA Astrophysics Data System (ADS)

Breuer, M.; Jaffrézic, B.; Arora, K.

2008-05-01

In order to reduce the high computational effort of wall-resolved large-eddy simulations (LES), the present paper suggests a hybrid LES RANS approach which splits up the simulation into a near-wall RANS part and an outer LES part. Generally, RANS is adequate for attached boundary layers requiring reasonable CPU-time and memory, where LES can also be applied but demands extremely large resources. Contrarily, RANS often fails in flows with massive separation or large-scale vortical structures. Here, LES is without a doubt the best choice. The basic concept of hybrid methods is to combine the advantages of both approaches yielding a prediction method, which, on the one hand, assures reliable results for complex turbulent flows, including large-scale flow phenomena and massive separation, but, on the other hand, consumes much fewer resources than LES, especially for high Reynolds number flows encountered in technical applications. In the present study, a non-zonal hybrid technique is considered (according to the signification retained by the authors concerning the terms zonal and non-zonal), which leads to an approach where the suitable simulation technique is chosen more or less automatically. For this purpose the hybrid approach proposed relies on a unique modeling concept. In the LES mode a subgrid-scale model based on a one-equation model for the subgrid-scale turbulent kinetic energy is applied, where the length scale is defined by the filter width. For the viscosity-affected near-wall RANS mode the one-equation model proposed by Rodi et al. (J Fluids Eng 115:196 205, 1993) is used, which is based on the wall-normal velocity fluctuations as the velocity scale and algebraic relations for the length scales. Although the idea of combined LES RANS methods is not new, a variety of open questions still has to be answered. This includes, in particular, the demand for appropriate coupling techniques between LES and RANS, adaptive control mechanisms, and proper subgrid-scale and RANS models. Here, in addition to the study on the behavior of the suggested hybrid LES RANS approach, special emphasis is put on the investigation of suitable interface criteria and the adjustment of the RANS model. To investigate these issues, two different test cases are considered. Besides the standard plane channel flow test case, the flow over a periodic arrangement of hills is studied in detail. This test case includes a pressure-induced flow separation and subsequent reattachment. In comparison with a wall-resolved LES prediction encouraging results are achieved.

PARALLEL PERTURBATION MODEL FOR CYCLE TO CYCLE VARIABILITY PPM4CCV

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ameen, Muhsin Mohammed; Som, Sibendu

This code consists of a Fortran 90 implementation of the parallel perturbation model to compute cyclic variability in spark ignition (SI) engines. Cycle-to-cycle variability (CCV) is known to be detrimental to SI engine operation resulting in partial burn and knock, and result in an overall reduction in the reliability of the engine. Numerical prediction of cycle-to-cycle variability (CCV) in SI engines is extremely challenging for two key reasons: (i) high-fidelity methods such as large eddy simulation (LES) are required to accurately capture the in-cylinder turbulent flow field, and (ii) CCV is experienced over long timescales and hence the simulations needmore » to be performed for hundreds of consecutive cycles. In the new technique, the strategy is to perform multiple parallel simulations, each of which encompasses 2-3 cycles, by effectively perturbing the simulation parameters such as the initial and boundary conditions. The PPM4CCV code is a pre-processing code and can be coupled with any engine CFD code. PPM4CCV was coupled with Converge CFD code and a 10-time speedup was demonstrated over the conventional multi-cycle LES in predicting the CCV for a motored engine. Recently, the model is also being applied to fired engines including port fuel injected (PFI) and direct injection spark ignition engines and the preliminary results are very encouraging.« less

Study of the flow field past dimpled aerodynamic surfaces: numerical simulation and experimental verification

NASA Astrophysics Data System (ADS)

Binci, L.; Clementi, G.; D'Alessandro, V.; Montelpare, S.; Ricci, R.

2017-11-01

This work presents the study of the flow field past of dimpled laminar airfoil. Fluid dynamic behaviour of these elements has been not still deeply studied in the scientific community. Therefore Computational Fluid-Dynamics (CFD) is here used to analyze the flow field induced by dimples on the NACA 64-014A laminar airfoil at Re = 1.75 · 105 at α = 0°. Reynolds Averaged Navier-Stokes (RANS) equations and Large-Eddy Simulations (LES) were compared with wind tunnel measurements in order to evaluate their effectiveness in the modeling this kind of flow field. LES equations were solved using a specifically developed OpenFOAM solver adopting an L-stable Singly Diagonally Implicit Runge-Kutta (SDIRK) technique with an iterated PISO-like procedure for handling pressure-velocity coupling within each RK stage. Dynamic Smagorinsky subgrid model was employed. LES results provided good agreement with experimental data, while RANS equations closed with \\[k-ω -γ -\\overset{}{\\mathop{{{\\operatorname{Re}}θ, \\text{t}}}} \\] approach overstimates laminar separation bubble (LSB) extension of dimpled and un-dimpled configurations. Moreover, through skin friction coefficient analysis, we found a different representation of the turbulent zone between the numerical models; indeed, with RANS model LSB seems to be divided in two different parts, meanwhile LES model shows a LSB global reduction.

Self-Employment among Same-Sex and Opposite-Sex Couples in Canada.

PubMed

Waite, Sean; Denier, Nicole

2016-05-01

This study presents novel evidence on the relationship between sexual orientation and self-employment. Using data from the 2001 and 2006 Census of Canada and the 2011 Canadian National Household Survey, we explore the propensity for self-employment among same- and opposite-sex couples. We examine the demographic, human capital, and family characteristics of coupled gay men and lesbians relative to their coupled heterosexual counterparts to offer potential mechanisms generating differences in rates of self-employment. Our analysis further considers occupational variability in the likelihood of self-employment. We find that gay men are less likely and lesbians more likely than heterosexuals to be self-employed; however, there is significant variation across occupations. Gay men are more likely to be self-employed in arts and culture, sales and service, and natural and applied sciences, but less likely in business, finance, and health-related occupations. Lesbians are much more likely to be self-employed in health-related occupations, natural and applied sciences, and arts and culture. Marriage and having children are significant predictors of self-employment for coupled heterosexual women but not lesbians. Cette étude présente des évidences empiriques concernant la relation entre l'orientation sexuelle et le travail indépendant. Utilisant des données provenant du Recensement du Canada de 2001 et de 2006, ainsi que l'Enquête nationale auprès des ménages (ENM) de 2011, nous explorons la tendance du travail indépendant parmi les couples de même sexe et ceux de sexe opposé. Ainsi, nous examinons les caractéristiques démographiques, du capital humain et familiales des couples gais et lesbiens par rapport à leurs homologues hétérosexuels, afin de démontrer une corrélation entre l'orientation sexuelle et la probabilité d'être travailleur indépendant - une causation qui nous semble évident et que nous analysons plus en profondeur. Nous concluons d'ailleurs que les hommes gais sont moins probables - et les lesbiennes plus probables - d'être travailleurs indépendants que les hétérosexuels. Cependant, il existe des variantes significatives entre les différentes occupations. Les hommes homosexuels tendent plutôt à être travailleurs indépendants dans les domaines des arts et de la culture, la vente et les services, ainsi que les sciences pures et les sciences naturelles. Par contre, ils le sont moins dans les domaines des affaires, de la finance et des domaines reliés à la santé. De leur côté, les femmes homosexuelles tendent plutôt à être travailleuses indépendantes dans les domaines reliés à la santé, aux sciences pures et aux sciences naturelles, ainsi qu'aux arts et la culture. Finalement, le fait d'être marié et d'avoir des enfants sont des indicateurs significatifs du travail indépendant pour les femmes en couples hétérosexuels, mais pas pour les lesbiennes. © 2016 Canadian Sociological Association/La Société canadienne de sociologie.

Large-eddy simulations of a solid-rocket booster jet

NASA Astrophysics Data System (ADS)

Paoli, Roberto; Poubeau, Adele; Cariolle, Daniel

2014-11-01

Emissions from solid-rocket boosters are responsible for a severe decrease in ozone concentration in the rocket plume during the first hours after a launch. The main source of ozone depletion is due to hydrogen chloride that is converted into chlorine in the high temperature regions of the jet (afterburning). The objective of this study is to evaluate the active chlorine concentration in the plume of a solid-rocket booster using large-eddy simulations. The gas is injected through the entire nozzle of the booster and a local time-stepping method based on coupling multi-instances of a fluid solver is used to extend the computational domain up to 600 nozzle exit diameters. The methodology is validated for a non-reactive case by analyzing the flow characteristics of supersonic co-flowing under expanded jets. Then, the chemistry of chlorine is studied offline using a complex chemistry solver and the LES data extracted from the mean trajectories of sample fluid particles. Finally, the online chemistry is analyzed by means of the multispecies version of the LES solver using a reduced chemistry scheme. The LES are able to capture the mixing of the exhaust with ambient air and the species concentrations, which is also useful to initialize atmospheric simulations on larger domains.

PEVC-FMDF for Large Eddy Simulation of Compressible Turbulent Flows

NASA Astrophysics Data System (ADS)

Nouri Gheimassi, Arash; Nik, Mehdi; Givi, Peyman; Livescu, Daniel; Pope, Stephen

2017-11-01

The filtered density function (FDF) closure is extended to a ``self-contained'' format to include the subgrid scale (SGS) statistics of all of the hydro-thermo-chemical variables in turbulent flows. These are the thermodynamic pressure, the specific internal energy, the velocity vector, and the composition field. In this format, the model is comprehensive and facilitates large eddy simulation (LES) of flows at both low and high compressibility levels. A transport equation is developed for the joint ``pressure-energy-velocity-composition filtered mass density function (PEVC-FMDF).'' In this equation, the effect of convection appears in closed form. The coupling of the hydrodynamics and thermochemistry is modeled via a set of stochastic differential equation (SDE) for each of the transport variables. This yields a self-contained SGS closure. For demonstration, LES is conducted of a turbulent shear flow with transport of a passive scalar. The consistency of the PEVC-FMDF formulation is established, and its overall predictive capability is appraised via comparison with direct numerical simulation (DNS) data.

LES-ODT Simulations of Turbulent Reacting Shear Layers

NASA Astrophysics Data System (ADS)

Hoffie, Andreas; Echekki, Tarek

2012-11-01

Large-eddy simulations (LES) combined with the one-dimensional turbulence (ODT) simulations of a spatially developing turbulent reacting shear layer with heat release and high Reynolds numbers were conducted and compared to results from direct numerical simulations (DNS) of the same configuration. The LES-ODT approach is based on LES solutions for momentum on a coarse grid and solutions for momentum and reactive scalars on a fine ODT grid, which is embedded in the LES computational domain. The shear layer is simulated with a single-step, second-order reaction with an Arrhenius reaction rate. The transport equations are solved using a low Mach number approximation. The LES-ODT simulations yield reasonably accurate predictions of turbulence and passive/reactive scalars' statistics compared to DNS results.

Using Reconstructed POD Modes as Turbulent Inflow for LES Wind Turbine Simulations

NASA Astrophysics Data System (ADS)

Nielson, Jordan; Bhaganagar, Kiran; Juttijudata, Vejapong; Sirisup, Sirod

2016-11-01

Currently, in order to get realistic atmospheric effects of turbulence, wind turbine LES simulations require computationally expensive precursor simulations. At times, the precursor simulation is more computationally expensive than the wind turbine simulation. The precursor simulations are important because they capture turbulence in the atmosphere and as stated above, turbulence impacts the power production estimation. On the other hand, POD analysis has been shown to be capable of capturing turbulent structures. The current study was performed to determine the plausibility of using lower dimension models from POD analysis of LES simulations as turbulent inflow to wind turbine LES simulations. The study will aid the wind energy community by lowering the computational cost of full scale wind turbine LES simulations, while maintaining a high level of turbulent information and being able to quickly apply the turbulent inflow to multi turbine wind farms. This will be done by comparing a pure LES precursor wind turbine simulation with simulations that use reduced POD mod inflow conditions. The study shows the feasibility of using lower dimension models as turbulent inflow of LES wind turbine simulations. Overall the power production estimation and velocity field of the wind turbine wake are well captured with small errors.

Comparative evaluation of polarimetric and bi-spectral cloud microphysics retrievals: Retrieval closure experiments and comparisons based on idealized and LES case studies

NASA Astrophysics Data System (ADS)

Miller, D. J.; Zhang, Z.; Ackerman, A. S.; Platnick, S. E.; Cornet, C.

2016-12-01

A remote sensing cloud retrieval simulator, created by coupling an LES cloud model with vector radiative transfer (RT) models is the ideal framework for assessing cloud remote sensing techniques. This simulator serves as a tool for understanding bi-spectral and polarimetric retrievals by comparing them directly to LES cloud properties (retrieval closure comparison) and for comparing the retrieval techniques to one another. Our simulator utilizes the DHARMA LES [Ackerman et al., 2004] with cloud properties based on marine boundary layer (MBL) clouds observed during the DYCOMS-II and ATEX field campaigns. The cloud reflectances are produced by the vectorized RT models based on polarized doubling adding and monte carlo techniques (PDA, MCPOL). Retrievals are performed utilizing techniques as similar as possible to those implemented on their corresponding well known instruments; polarimetric retrievals are based on techniques implemented for polarimeters (POLDER, AirMSPI, and RSP) and bi-spectral retrievals are performed using the Nakajima-King LUT method utilized on a number of spectral instruments (MODIS and VIIRS). Retrieval comparisons focus on cloud droplet effective radius (re), effective variance (ve), and cloud optical thickness (τ). This work explores the sensitivities of these two retrieval techniques to various observation limitations, such as spatial resolution/cloud inhomogeneity, impact of 3D radiative effects, and angular resolution requirements. With future remote sensing missions like NASA's Aerosols/Clouds/Ecosystems (ACE) planning to feature advanced polarimetric instruments it is important to understand how these retrieval techniques compare to one another. The cloud retrieval simulator we've developed allows us to probe these important questions in a realistically relevant test bed.

Validating LES for Jet Aeroacoustics

NASA Technical Reports Server (NTRS)

Bridges, James; Wernet, Mark P.

2011-01-01

Engineers charged with making jet aircraft quieter have long dreamed of being able to see exactly how turbulent eddies produce sound and this dream is now coming true with the advent of large eddy simulation (LES). Two obvious challenges remain: validating the LES codes at the resolution required to see the fluid-acoustic coupling, and the interpretation of the massive datasets that are produced. This paper addresses the former, the use of advanced experimental techniques such as particle image velocimetry (PIV) and Raman and Rayleigh scattering, to validate the computer codes and procedures used to create LES solutions. This paper argues that the issue of accuracy of the experimental measurements be addressed by cross-facility and cross-disciplinary examination of modern datasets along with increased reporting of internal quality checks in PIV analysis. Further, it argues that the appropriate validation metrics for aeroacoustic applications are increasingly complicated statistics that have been shown in aeroacoustic theory to be critical to flow-generated sound, such as two-point space-time velocity correlations. A brief review of data sources available is presented along with examples illustrating cross-facility and internal quality checks required of the data before it should be accepted for validation of LES.

Five-equation and robust three-equation methods for solution verification of large eddy simulation

NASA Astrophysics Data System (ADS)

Dutta, Rabijit; Xing, Tao

2018-02-01

This study evaluates the recently developed general framework for solution verification methods for large eddy simulation (LES) using implicitly filtered LES of periodic channel flows at friction Reynolds number of 395 on eight systematically refined grids. The seven-equation method shows that the coupling error based on Hypothesis I is much smaller as compared with the numerical and modeling errors and therefore can be neglected. The authors recommend five-equation method based on Hypothesis II, which shows a monotonic convergence behavior of the predicted numerical benchmark ( S C ), and provides realistic error estimates without the need of fixing the orders of accuracy for either numerical or modeling errors. Based on the results from seven-equation and five-equation methods, less expensive three and four-equation methods for practical LES applications were derived. It was found that the new three-equation method is robust as it can be applied to any convergence types and reasonably predict the error trends. It was also observed that the numerical and modeling errors usually have opposite signs, which suggests error cancellation play an essential role in LES. When Reynolds averaged Navier-Stokes (RANS) based error estimation method is applied, it shows significant error in the prediction of S C on coarse meshes. However, it predicts reasonable S C when the grids resolve at least 80% of the total turbulent kinetic energy.

Development of a Wind Plant Large-Eddy Simulation with Measurement-Driven Atmospheric Inflow: Preprint

DOE Office of Scientific and Technical Information (OSTI.GOV)

Quon, Eliot; Churchfield, Matthew; Cheung, Lawrence

This paper details the development of an aeroelastic wind plant model with large-eddy simulation (LES). The chosen LES solver is the Simulator for Wind Farm Applications (SOWFA) based on the OpenFOAM framework, coupled to NREL's comprehensive aeroelastic analysis tool, FAST. An atmospheric boundary layer (ABL) precursor simulation was constructed based on assessments of meteorological tower, lidar, and radar data over a 3-hour window. This precursor was tuned to the specific atmospheric conditions that occurred both prior to and during the measurement campaign, enabling capture of a night-to-day transition in the turbulent ABL. In the absence of height-varying temperature measurements, spatiallymore » averaged radar data were sufficient to characterize the atmospheric stability of the wind plant in terms of the shear profile, and near-ground temperature sensors provided a reasonable estimate of the ground heating rate describing the morning transition. A full aeroelastic simulation was then performed for a subset of turbines within the wind plant, driven by the precursor. Analysis of two turbines within the array, one directly waked by the other, demonstrated good agreement with measured time-averaged loads.« less

Development of a Wind Plant Large-Eddy Simulation with Measurement-Driven Atmospheric Inflow

DOE Office of Scientific and Technical Information (OSTI.GOV)

Quon, Eliot W.; Churchfield, Matthew J.; Cheung, Lawrence

This paper details the development of an aeroelastic wind plant model with large-eddy simulation (LES). The chosen LES solver is the Simulator for Wind Farm Applications (SOWFA) based on the OpenFOAM framework, coupled to NREL's comprehensive aeroelastic analysis tool, FAST. An atmospheric boundary layer (ABL) precursor simulation was constructed based on assessments of meteorological tower, lidar, and radar data over a 3-hour window. This precursor was tuned to the specific atmospheric conditions that occurred both prior to and during the measurement campaign, enabling capture of a night-to-day transition in the turbulent ABL. In the absence of height-varying temperature measurements, spatiallymore » averaged radar data were sufficient to characterize the atmospheric stability of the wind plant in terms of the shear profile, and near-ground temperature sensors provided a reasonable estimate of the ground heating rate describing the morning transition. A full aeroelastic simulation was then performed for a subset of turbines within the wind plant, driven by the precursor. Analysis of two turbines within the array, one directly waked by the other, demonstrated good agreement with measured time-averaged loads.« less

Effect of artificial length scales in large eddy simulation of a neutral atmospheric boundary layer flow: A simple solution to log-layer mismatch

NASA Astrophysics Data System (ADS)

Chatterjee, Tanmoy; Peet, Yulia T.

2017-07-01

A large eddy simulation (LES) methodology coupled with near-wall modeling has been implemented in the current study for high Re neutral atmospheric boundary layer flows using an exponentially accurate spectral element method in an open-source research code Nek 5000. The effect of artificial length scales due to subgrid scale (SGS) and near wall modeling (NWM) on the scaling laws and structure of the inner and outer layer eddies is studied using varying SGS and NWM parameters in the spectral element framework. The study provides an understanding of the various length scales and dynamics of the eddies affected by the LES model and also the fundamental physics behind the inner and outer layer eddies which are responsible for the correct behavior of the mean statistics in accordance with the definition of equilibrium layers by Townsend. An economical and accurate LES model based on capturing the near wall coherent eddies has been designed, which is successful in eliminating the artificial length scale effects like the log-layer mismatch or the secondary peak generation in the streamwise variance.

PNNL - WRF-LES - Convective - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a convective atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on July 4, 2012. The dataset was used to assess the LES models for simulation of canonical convective ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

ANL - WRF-LES - Convective - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a convective atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on July 4, 2012. The dataset was used to assess the LES models for simulation of canonical convective ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

LLNL - WRF-LES - Neutral - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a neutrally stratified atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on Aug. 17, 2012. The dataset was used to assess LES models for simulation of canonical neutral ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

ANL - WRF-LES - Neutral - TTU

DOE Data Explorer

Kosovic, Branko

2018-06-20

This dataset includes large-eddy simulation (LES) output from a neutrally stratified atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on Aug. 17, 2012. The dataset was used to assess LES models for simulation of canonical neutral ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

LANL - WRF-LES - Neutral - TTU

DOE Data Explorer

Kosovic, Branko

2018-06-20

This dataset includes large-eddy simulation (LES) output from a neutrally stratified atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on Aug. 17, 2012. The dataset was used to assess LES models for simulation of canonical neutral ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

LANL - WRF-LES - Convective - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a convective atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on July 4, 2012. The dataset was used to assess the LES models for simulation of canonical convective ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

Hydrodynamic Drag Reduction

DTIC Science & Technology

2015-04-01

Computational Engineering unstructured RANS/LES/DES solver , Tenasi, was used to predict drag and simulate the free surface flow around the ACV over a...using a second-order accurate Roe approximate Riemann scheme, while viscous fluxes are evaluated using a second-order directional derivative approach...Predictions of rigid body ship motions for the SI75 container ship in incident waves and methodology for a one-way coupling of the Tenasi flow solver

Sand Waves in Environmental Flows: Insights gained by LES

NASA Astrophysics Data System (ADS)

Sotiropoulos, Fotis

2014-11-01

In fluvial and coastal environments, sediment transport processes induced by near-bed coherent structures in the turbulent boundary layer developing over a mobile sediment bed result in the formation of dynamically rich sand waves, or bed forms, which grow and migrate continuously. Bed form migration alters streambed roughness and provides the primary mechanism for transporting large amounts of sediment through riverine systems impacting the morphology, streambank stability, and ecology of waterways. I will present recent computational advances, which have enabled coupled, hydro-morphodynamic large-eddy simulation (LES) of turbulent flow in mobile-bed open channels. Numerical simulations: 1) elucidate the role of near-bed sweeps in the turbulent boundary layer as the mechanism for initiating the instability of the initially flat sand bed; 2) show how near-bed processes give rise to aperiodic eruptions of suspended sediment at the free surface; and 3) clarify the mechanism via which sand waves migrate. Furthermore, in agreement with recent experimental observations, the computed spectra of the resolved velocity fluctuations above the bed exhibit a distinct spectral gap whose width increases with distance from the bed. The spectral gap delineates the spectrum of turbulence from that of slowly evolving coherent structures associated with sand wave migration. The talk will also present computational results demonstrating the feasibility of carrying out coupled, hydro-morphodynamic LES of large dunes migrating in meandering streams and rivers with embedded hydraulic structures and discuss future challenges and opportunities. This work was supported by NSF Grants EAR-0120914 and EAR-0738726, and National Cooperative Highway Research Program Grant NCHRP-HR 24-33.

Numerical simulation of sloshing with large deforming free surface by MPS-LES method

NASA Astrophysics Data System (ADS)

Pan, Xu-jie; Zhang, Huai-xin; Sun, Xue-yao

2012-12-01

Moving particle semi-implicit (MPS) method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation (LES) by coupling with a sub-particle-scale (SPS) turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is "neighbor particle". In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no "neighbor particles" in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

Nested mesoscale-to-LES modeling of the atmospheric boundary layer in the presence of under-resolved convective structures

DOE PAGES

Mazzaro, Laura J.; Munoz-Esparza, Domingo; Lundquist, Julie K.; ...

2017-07-06

Multiscale atmospheric simulations can be computationally prohibitive, as they require large domains and fine spatiotemporal resolutions. Grid-nesting can alleviate this by bridging mesoscales and microscales, but one turbulence scheme must run at resolutions within a range of scales known as the terra incognita (TI). TI grid-cell sizes can violate both mesoscale and microscale subgrid-scale parametrization assumptions, resulting in unrealistic flow structures. Herein we assess the impact of unrealistic lateral boundary conditions from parent mesoscale simulations at TI resolutions on nested large eddy simulations (LES), to determine whether parent domains bias the nested LES. We present a series of idealized nestedmore » mesoscale-to-LES runs of a dry convective boundary layer (CBL) with different parent resolutions in the TI. We compare the nested LES with a stand-alone LES with periodic boundary conditions. The nested LES domains develop ~20% smaller convective structures, while potential temperature profiles are nearly identical for both the mesoscales and LES simulations. The horizontal wind speed and surface wind shear in the nested simulations closely resemble the reference LES. Heat fluxes are overestimated by up to ~0.01 K m s –1 in the top half of the PBL for all nested simulations. Overestimates of turbulent kinetic energy (TKE) and Reynolds stress in the nested domains are proportional to the parent domain's grid-cell size, and are almost eliminated for the simulation with the finest parent grid-cell size. Furthermore, based on these results, we recommend that LES of the CBL be forced by mesoscale simulations with the finest practical resolution.« less

Nested mesoscale-to-LES modeling of the atmospheric boundary layer in the presence of under-resolved convective structures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mazzaro, Laura J.; Munoz-Esparza, Domingo; Lundquist, Julie K.

Multiscale atmospheric simulations can be computationally prohibitive, as they require large domains and fine spatiotemporal resolutions. Grid-nesting can alleviate this by bridging mesoscales and microscales, but one turbulence scheme must run at resolutions within a range of scales known as the terra incognita (TI). TI grid-cell sizes can violate both mesoscale and microscale subgrid-scale parametrization assumptions, resulting in unrealistic flow structures. Herein we assess the impact of unrealistic lateral boundary conditions from parent mesoscale simulations at TI resolutions on nested large eddy simulations (LES), to determine whether parent domains bias the nested LES. We present a series of idealized nestedmore » mesoscale-to-LES runs of a dry convective boundary layer (CBL) with different parent resolutions in the TI. We compare the nested LES with a stand-alone LES with periodic boundary conditions. The nested LES domains develop ~20% smaller convective structures, while potential temperature profiles are nearly identical for both the mesoscales and LES simulations. The horizontal wind speed and surface wind shear in the nested simulations closely resemble the reference LES. Heat fluxes are overestimated by up to ~0.01 K m s –1 in the top half of the PBL for all nested simulations. Overestimates of turbulent kinetic energy (TKE) and Reynolds stress in the nested domains are proportional to the parent domain's grid-cell size, and are almost eliminated for the simulation with the finest parent grid-cell size. Furthermore, based on these results, we recommend that LES of the CBL be forced by mesoscale simulations with the finest practical resolution.« less

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

NASA Astrophysics Data System (ADS)

Kim, Jeonglae; Pope, Stephen B.

2014-05-01

A turbulent lean-premixed propane-air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.

Investigating the Impact of Surface Heterogeneity on the Convective Boundary Layer Over Urban Areas Through Coupled Large-Eddy Simulation and Remote Sensing

NASA Technical Reports Server (NTRS)

Dominguez, Anthony; Kleissl, Jan P.; Luvall, Jeffrey C.

2011-01-01

Large-eddy Simulation (LES) was used to study convective boundary layer (CBL) flow through suburban regions with both large and small scale heterogeneities in surface temperature. Constant remotely sensed surface temperatures were applied at the surface boundary at resolutions of 10 m, 90 m, 200 m, and 1 km. Increasing the surface resolution from 1 km to 200 m had the most significant impact on the mean and turbulent flow characteristics as the larger scale heterogeneities became resolved. While previous studies concluded that scales of heterogeneity much smaller than the CBL inversion height have little impact on the CBL characteristics, we found that further increasing the surface resolution (resolving smaller scale heterogeneities) results in an increase in mean surface heat flux, thermal blending height, and potential temperature profile. The results of this study will help to better inform sub-grid parameterization for meso-scale meteorological models. The simulation tool developed through this study (combining LES and high resolution remotely sensed surface conditions) is a significant step towards future studies on the micro-scale meteorology in urban areas.

Performance assessment of Large Eddy Simulation (LES) for modeling dispersion in an urban street canyon with tree planting

NASA Astrophysics Data System (ADS)

Moonen, P.; Gromke, C.; Dorer, V.

2013-08-01

The potential of a Large Eddy Simulation (LES) model to reliably predict near-field pollutant dispersion is assessed. To that extent, detailed time-resolved numerical simulations of coupled flow and dispersion are conducted for a street canyon with tree planting. Different crown porosities are considered. The model performance is assessed in several steps, ranging from a qualitative comparison to measured concentrations, over statistical data analysis by means of scatter plots and box plots, up to the calculation of objective validation metrics. The extensive validation effort highlights and quantifies notable features and shortcomings of the model, which would otherwise remain unnoticed. The model performance is found to be spatially non-uniform. Closer agreement with measurement data is achieved near the canyon ends than for the central part of the canyon, and typical model acceptance criteria are satisfied more easily for the leeward than for the windward canyon wall. This demonstrates the need for rigorous model evaluation. Only quality-assured models can be used with confidence to support assessment, planning and implementation of pollutant mitigation strategies.

Multiscale Data Assimilation for Large-Eddy Simulations

NASA Astrophysics Data System (ADS)

Li, Z.; Cheng, X.; Gustafson, W. I., Jr.; Xiao, H.; Vogelmann, A. M.; Endo, S.; Toto, T.

2017-12-01

Large-eddy simulation (LES) is a powerful tool for understanding atmospheric turbulence, boundary layer physics and cloud development, and there is a great need for developing data assimilation methodologies that can constrain LES models. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility has been developing the capability to routinely generate ensembles of LES. The LES ARM Symbiotic Simulation and Observation (LASSO) project (https://www.arm.gov/capabilities/modeling/lasso) is generating simulations for shallow convection days at the ARM Southern Great Plains site in Oklahoma. One of major objectives of LASSO is to develop the capability to observationally constrain LES using a hierarchy of ARM observations. We have implemented a multiscale data assimilation (MSDA) scheme, which allows data assimilation to be implemented separately for distinct spatial scales, so that the localized observations can be effectively assimilated to constrain the mesoscale fields in the LES area of about 15 km in width. The MSDA analysis is used to produce forcing data that drive LES. With such LES workflow we have examined 13 days with shallow convection selected from the period May-August 2016. We will describe the implementation of MSDA, present LES results, and address challenges and opportunities for applying data assimilation to LES studies.

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

NASA Technical Reports Server (NTRS)

Bellan, Josette; Taskinoglu, Ezgi

2012-01-01

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Validating LES for Jet Aeroacoustics

NASA Technical Reports Server (NTRS)

Bridges, James

2011-01-01

Engineers charged with making jet aircraft quieter have long dreamed of being able to see exactly how turbulent eddies produce sound and this dream is now coming true with the advent of large eddy simulation (LES). Two obvious challenges remain: validating the LES codes at the resolution required to see the fluid-acoustic coupling, and the interpretation of the massive datasets that result in having dreams come true. This paper primarily addresses the former, the use of advanced experimental techniques such as particle image velocimetry (PIV) and Raman and Rayleigh scattering, to validate the computer codes and procedures used to create LES solutions. It also addresses the latter problem in discussing what are relevant measures critical for aeroacoustics that should be used in validating LES codes. These new diagnostic techniques deliver measurements and flow statistics of increasing sophistication and capability, but what of their accuracy? And what are the measures to be used in validation? This paper argues that the issue of accuracy be addressed by cross-facility and cross-disciplinary examination of modern datasets along with increased reporting of internal quality checks in PIV analysis. Further, it is argued that the appropriate validation metrics for aeroacoustic applications are increasingly complicated statistics that have been shown in aeroacoustic theory to be critical to flow-generated sound.

Evapotranspiration and cloud variability at regional sub-grid scales

NASA Astrophysics Data System (ADS)

Vila-Guerau de Arellano, Jordi; Sikma, Martin; Pedruzo-Bagazgoitia, Xabier; van Heerwaarden, Chiel; Hartogensis, Oscar; Ouwersloot, Huug

2017-04-01

In regional and global models uncertainties arise due to our incomplete understanding of the coupling between biochemical and physical processes. Representing their impact depends on our ability to calculate these processes using physically sound parameterizations, since they are unresolved at scales smaller than the grid size. More specifically over land, the coupling between evapotranspiration, turbulent transport of heat and moisture, and clouds lacks a combined representation to take these sub-grid scales interactions into account. Our approach is based on understanding how radiation, surface exchange, turbulent transport and moist convection are interacting from the leaf- to the cloud scale. We therefore place special emphasis on plant stomatal aperture as the main regulator of CO2-assimilation and water transpiration, a key source of moisture source to the atmosphere. Plant functionality is critically modulated by interactions with atmospheric conditions occurring at very short spatiotemporal scales such as cloud radiation perturbations or water vapour turbulent fluctuations. By explicitly resolving these processes, the LES (large-eddy simulation) technique is enabling us to characterize and better understand the interactions between canopies and the local atmosphere. This includes the adaption time of vegetation to rapid changes in atmospheric conditions driven by turbulence or the presence of cumulus clouds. Our LES experiments are based on explicitly coupling the diurnal atmospheric dynamics to a plant physiology model. Our general hypothesis is that different partitioning of direct and diffuse radiation leads to different responses of the vegetation. As a result there are changes in the water use efficiencies and shifts in the partitioning of sensible and latent heat fluxes under the presence of clouds. Our presentation is as follows. First, we discuss the ability of LES to reproduce the surface energy balance including photosynthesis and CO2 soil respiration coupled to the dynamics of a convective boundary layer. LES results are compared with a complete set of surface and upper-air meteorological and carbon-dioxide observations gathered during a representative day at the 213-meter meteorological tall tower at Cabauw. Second, we perform systematic numerical experiments under a wide range of background wind conditions and stomatal aperture response time. Our analysis unravel how thin clouds, characterized by lower values of the cloud optical depth, have a different impact on evapotranspiration compared to thick clouds due to differences in the partitioning between direct and diffuse radiation at canopy level. Related to this detailed simulation, we discuss how new instrumental techniques, e.g. scintillometery, enable us to obtain new observational insight of the coupling between clouds and vegetation. We will close the presentation with open questions regarding the need to include parameterizations for these interactions at short spatiotemporal scales in regional or climate models.

LES-based filter-matrix lattice Boltzmann model for simulating fully developed turbulent channel flow

NASA Astrophysics Data System (ADS)

Zhuo, Congshan; Zhong, Chengwen

2016-11-01

In this paper, a three-dimensional filter-matrix lattice Boltzmann (FMLB) model based on large eddy simulation (LES) was verified for simulating wall-bounded turbulent flows. The Vreman subgrid-scale model was employed in the present FMLB-LES framework, which had been proved to be capable of predicting turbulent near-wall region accurately. The fully developed turbulent channel flows were performed at a friction Reynolds number Reτ of 180. The turbulence statistics computed from the present FMLB-LES simulations, including mean stream velocity profile, Reynolds stress profile and root-mean-square velocity fluctuations greed well with the LES results of multiple-relaxation-time (MRT) LB model, and some discrepancies in comparison with those direct numerical simulation (DNS) data of Kim et al. was also observed due to the relatively low grid resolution. Moreover, to investigate the influence of grid resolution on the present LES simulation, a DNS simulation on a finer gird was also implemented by present FMLB-D3Q19 model. Comparisons of detailed computed various turbulence statistics with available benchmark data of DNS showed quite well agreement.

Coupled large-eddy simulation and morphodynamics of a large-scale river under extreme flood conditions

NASA Astrophysics Data System (ADS)

Khosronejad, Ali; Sotiropoulos, Fotis; Stony Brook University Team

2016-11-01

We present a coupled flow and morphodynamic simulations of extreme flooding in 3 km long and 300 m wide reach of the Mississippi River in Minnesota, which includes three islands and hydraulic structures. We employ the large-eddy simulation (LES) and bed-morphodynamic modules of the VFS-Geophysics model to investigate the flow and bed evolution of the river during a 500 year flood. The coupling of the two modules is carried out via a fluid-structure interaction approach using a nested domain approach to enhance the resolution of bridge scour predictions. The geometrical data of the river, islands and structures are obtained from LiDAR, sub-aqueous sonar and in-situ surveying to construct a digital map of the river bathymetry. Our simulation results for the bed evolution of the river reveal complex sediment dynamics near the hydraulic structures. The numerically captured scour depth near some of the structures reach a maximum of about 10 m. The data-driven simulation strategy we present in this work exemplifies a practical simulation-based-engineering-approach to investigate the resilience of infrastructures to extreme flood events in intricate field-scale riverine systems. This work was funded by a Grant from Minnesota Dept. of Transportation.

Self-contained filtered density function

DOE PAGES

Nouri, Arash G.; Nik, Mehdi B.; Givi, Pope; ...

2017-09-18

The filtered density function (FDF) closure is extended to a “self-contained” format to include the subgrid-scale (SGS) statistics of all of the hydro-thermo-chemical variables in turbulent flows. These are the thermodynamic pressure, the specific internal energy, the velocity vector, and the composition field. In this format, the model is comprehensive and facilitates large-eddy simulation (LES) of flows at both low and high compressibility levels. A transport equation is developed for the joint pressure-energy-velocity-composition filtered mass density function (PEVC-FMDF). In this equation, the effect of convection appears in closed form. The coupling of the hydrodynamics and thermochemistry is modeled via amore » set of stochastic differential equation for each of the transport variables. This yields a self-contained SGS closure. We demonstrated how LES is conducted of a turbulent shear flow with transport of a passive scalar. Finally, the consistency of the PEVC-FMDF formulation is established, and its overall predictive capability is appraised via comparison with direct numerical simulation (DNS) data.« less

Self-contained filtered density function

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nouri, Arash G.; Nik, Mehdi B.; Givi, Pope

The filtered density function (FDF) closure is extended to a “self-contained” format to include the subgrid-scale (SGS) statistics of all of the hydro-thermo-chemical variables in turbulent flows. These are the thermodynamic pressure, the specific internal energy, the velocity vector, and the composition field. In this format, the model is comprehensive and facilitates large-eddy simulation (LES) of flows at both low and high compressibility levels. A transport equation is developed for the joint pressure-energy-velocity-composition filtered mass density function (PEVC-FMDF). In this equation, the effect of convection appears in closed form. The coupling of the hydrodynamics and thermochemistry is modeled via amore » set of stochastic differential equation for each of the transport variables. This yields a self-contained SGS closure. We demonstrated how LES is conducted of a turbulent shear flow with transport of a passive scalar. Finally, the consistency of the PEVC-FMDF formulation is established, and its overall predictive capability is appraised via comparison with direct numerical simulation (DNS) data.« less

LES Modeling of Supersonic Combustion at SCRAMJET Conditions

NASA Astrophysics Data System (ADS)

Vane, Zachary; Lacaze, Guilhem; Oefelein, Joseph

2016-11-01

Results from a series of large-eddy simulations (LES) of the Hypersonic International Flight Research Experiment (HIFiRE) are examined with emphasis placed on the coupled performance of the wall and combustion models. The test case of interest corresponds to the geometry and conditions found in the ground based experiments performed in the HIFiRE Direct Connect Rig (HDCR) in dual-mode operation. In these calculations, the turbulence and mixing characteristics of the high Reynolds number turbulent boundary layer with multi-species fuel injection are analyzed using a simplified chemical model and combustion closure to predict the heat release measured experimentally. These simulations are then used to identify different flame regimes in the combustor section. Concurrently, the performance of an equilibrium wall-model is evaluated in the vicinity of the fuel injectors and in the flame-holding cavity where regions of boundary layer and thermochemical non-equilibrium are present. Support for this research was provided by the Defense Advanced Research Projects Agency (DARPA).

Self-contained filtered density function

NASA Astrophysics Data System (ADS)

Nouri, A. G.; Nik, M. B.; Givi, P.; Livescu, D.; Pope, S. B.

2017-09-01

The filtered density function (FDF) closure is extended to a "self-contained" format to include the subgrid-scale (SGS) statistics of all of the hydro-thermo-chemical variables in turbulent flows. These are the thermodynamic pressure, the specific internal energy, the velocity vector, and the composition field. In this format, the model is comprehensive and facilitates large-eddy simulation (LES) of flows at both low and high compressibility levels. A transport equation is developed for the joint pressure-energy-velocity-composition filtered mass density function (PEVC-FMDF). In this equation, the effect of convection appears in closed form. The coupling of the hydrodynamics and thermochemistry is modeled via a set of stochastic differential equation for each of the transport variables. This yields a self-contained SGS closure. For demonstration, LES is conducted of a turbulent shear flow with transport of a passive scalar. The consistency of the PEVC-FMDF formulation is established, and its overall predictive capability is appraised via comparison with direct numerical simulation (DNS) data.

A priori study of subgrid-scale flux of a passive scalar in isotropic homogeneous turbulence.

PubMed

Chumakov, Sergei G

2008-09-01

We perform a direct numerical simulation (DNS) of forced homogeneous isotropic turbulence with a passive scalar that is forced by mean gradient. The DNS data are used to study the properties of subgrid-scale flux of a passive scalar in the framework of large eddy simulation (LES), such as alignment trends between the flux, resolved, and subgrid-scale flow structures. It is shown that the direction of the flux is strongly coupled with the subgrid-scale stress axes rather than the resolved flow quantities such as strain, vorticity, or scalar gradient. We derive an approximate transport equation for the subgrid-scale flux of a scalar and look at the relative importance of the terms in the transport equation. A particular form of LES tensor-viscosity model for the scalar flux is investigated, which includes the subgrid-scale stress. Effect of different models for the subgrid-scale stress on the model for the subgrid-scale flux is studied.

A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing

NASA Technical Reports Server (NTRS)

Bellan, Josette; Taskinoglu, Ezgi

2012-01-01

Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.

High-Accuracy Near-Surface Large-Eddy Simulation with Planar Topography

DTIC Science & Technology

2015-08-03

Navier-Stokes equation, in effect randomizing the subfilter-scale (SFS) stress divergence. In the intervening years it has been discovered that this...surface stress models do introduce spurious effects that force deviations from LOTW at the first couple grid levels adjacent to the surface. Fig. 10 shows...SFS stress is sufficiently overwhelming to produce the overshoot. When the LES is moved into the HAZ so that the viscous effects causing the

Convective - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a convective atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on July 4, 2012. The dataset was used to assess the LES models for simulation of canonical convective ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

LANL - Convective - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a convective atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on July 4, 2012. The dataset was used to assess the LES models for simulation of canonical convective ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

LANL - Neutral - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a neutrally stratified atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on Aug. 17, 2012. The dataset was used to assess LES models for simulation of canonical neutral ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

Les violences conjugales à Dakar

PubMed Central

Soumah, Mohamed Maniboliot; Issa, Abdoul Wahab; Ndiaye, Mor; Ndoye, El Hadj Oumar; Sow, Mamadou Lamine

2015-01-01

L'objectif était d’évaluer les aspects épidémiologiques des violences conjugales, identifier les facteurs de risques et les différents types de violences conjugales, évaluer les conséquences des violences conjugales sur la santé des victimes, afin d'améliorer la prise en charge des victimes et la prévention du phénomène. Il s'est agit d'une étude transversale effectuée de décembre 2012 à janvier 2013 à Dakar. Les données ont été recueillies, après consentement, sur fiche d'enquête anonyme soumise à toute personne volontaire vivant en couple et résidant à Dakar. L'analyse statistique a été effectuée avec le logiciel SPSS 13.0. Le nombre de personnes victimes de violences conjugales était de 60 soit 37,30% dont 31 femmes (51,70%) et 29 hommes (48,30%). Le sex-ratio était de 0,93. Parmi les victimes, 53 étaient scolarisées soit 88,30%. Le régime matrimonial était de type monogame dans 39 cas (65%) et polygame dans 21 cas (35%). La vie en couple durait depuis moins de 11 ans dans 60% des cas et durait de 11 ans à 20 ans au plus dans 26,6% des cas. L’étude des types de violences montrait la fréquence des agressions physiques. Les armes utilisées étaient surtout les armes naturelles. Les principaux facteurs de risque de violence conjugale sont les facteurs sociodémographiques, culturels et économiques comme le jeune âge, l'inégalité du genre, les jeunes couples, la précarité, le niveau d'instruction élevé. La prise en charge des victimes et la prévention du phénomène restent insuffisantes dans nos pays. PMID:26918077

Detailed characteristics of drop-laden mixing layers: LES predictions compared to DNS

NASA Technical Reports Server (NTRS)

Okong'o, N.; Leboissetier, A.; Bellan, J.

2004-01-01

Results have been compared from Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) of a temporal mixing layer laden with evaporating drops, to assess the ability of LES to reproduce detailed characteristics of DNS.

NREL - SOWFA - Neutral - TTU

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kosovic, Branko

This dataset includes large-eddy simulation (LES) output from a neutrally stratified atmospheric boundary layer (ABL) simulation of observations at the SWIFT tower near Lubbock, Texas on Aug. 17, 2012. The dataset was used to assess LES models for simulation of canonical neutral ABL. The dataset can be used for comparison with other LES and computational fluid dynamics model outputs.

Comparison of Measured and Numerically Simulated Turbulence Statistics in a Convective Boundary Layer Over Complex Terrain

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rai, Raj K.; Berg, Larry K.; Kosović, Branko

High resolution numerical simulation can provide insight into important physical processes that occur within the planetary boundary layer (PBL). The present work employs large eddy simulation (LES) using the Weather Forecasting and Research (WRF) model, with the LES domain nested within mesoscale simulation, to simulate real conditions in the convective PBL over an area of complex terrain. A multiple nesting approach has been used to downsize the grid spacing from 12.15 km (mesoscale) to 0.03 km (LES). A careful selection of grid spacing in the WRF Meso domain has been conducted to minimize artifacts in the WRF-LES solutions. The WRF-LESmore » results have been evaluated with in situ and remote sensing observations collected during the US Department of Energy-supported Columbia BasinWind Energy Study (CBWES). Comparison of the first- and second-order moments, turbulence spectrum, and probability density function (PDF) of wind speed shows good agreement between the simulations and data. Furthermore, the WRF-LES variables show a great deal of variability in space and time caused by the complex topography in the LES domain. The WRF-LES results show that the flow structures, such as roll vortices and convective cells, vary depending on both the location and time of day. In addition to basic studies related to boundary-layer meteorology, results from these simulations can be used in other applications, such as studying wind energy resources, atmospheric dispersion, fire weather etc.« less

Impact of a highly detailed emission inventory on modeling accuracy

NASA Astrophysics Data System (ADS)

Taghavi, M.; Cautenet, S.; Arteta, J.

2005-03-01

During Expérience sur Site pour COntraindre les Modèles de Pollution atmosphérique et de Transport d'Emissions (ESCOMPTE) campaign (June 10 to July 14, 2001), two pollution events observed during an intensive measurement period (IOP2a and IOP2b) have been simulated. The comprehensive Regional Atmospheric Modeling Systems (RAMS) model, version 4.3, coupled online with a chemical module including 29 species is used to follow the chemistry of a polluted zone over Southern France. This online method takes advantage of a parallel code and use of the powerful computer SGI 3800. Runs are performed with two emission inventories: the Emission Pre Inventory (EPI) and the Main Emission Inventory (MEI). The latter is more recent and has a high resolution. The redistribution of simulated chemical species (ozone and nitrogen oxides) is compared with aircraft and surface station measurements for both runs at regional scale. We show that the MEI inventory is more efficient than the EPI in retrieving the redistribution of chemical species in space (three-dimensional) and time. In surface stations, MEI is superior especially for primary species, like nitrogen oxides. The ozone pollution peaks obtained from an inventory, such as EPI, have a large uncertainty. To understand the realistic geographical distribution of pollutants and to obtain a good order of magnitude in ozone concentration (in space and time), a high-resolution inventory like MEI is necessary. Coupling RAMS-Chemistry with MEI provides a very efficient tool able to simulate pollution plumes even in a region with complex circulations, such as the ESCOMPTE zone.

Mesure et retroaction sur un qubit multi-niveaux en electrodynamique quantique en circuit non lineair

NASA Astrophysics Data System (ADS)

Boissonneault, Maxime

L'electrodynamique quantique en circuit est une architecture prometteuse pour le calcul quantique ainsi que pour etudier l'optique quantique. Dans cette architecture, on couple un ou plusieurs qubits supraconducteurs jouant le role d'atomes a un ou plusieurs resonateurs jouant le role de cavites optiques. Dans cette these, j'etudie l'interaction entre un seul qubit supraconducteur et un seul resonateur, en permettant cependant au qubit d'avoir plus de deux niveaux et au resonateur d'avoir une non-linearite Kerr. Je m'interesse particulierement a la lecture de l'etat du qubit et a son amelioration, a la retroaction du processus de mesure sur le qubit de meme qu'a l'etude des proprietes quantiques du resonateur a l'aide du qubit. J'utilise pour ce faire un modele analytique reduit que je developpe a partir de la description complete du systeme en utilisant principalement des transfprmations unitaires et une elimination adiabatique. J'utilise aussi une librairie de calcul numerique maison permettant de simuler efficacement l'evolution du systeme complet. Je compare les predictions du modele analytique reduit et les resultats de simulations numeriques a des resultats experimentaux obtenus par l'equipe de quantronique du CEASaclay. Ces resultats sont ceux d'une spectroscopie d'un qubit supraconducteur couple a un resonateur non lineaire excite. Dans un regime de faible puissance de spectroscopie le modele reduit predit correctement la position et la largeur de la raie. La position de la raie subit les decalages de Lamb et de Stark, et sa largeur est dominee par un dephasage induit par le processus de mesure. Je montre que, pour les parametres typiques de l'electrodynamique quantique en circuit, un accord quantitatif requiert un modele en reponse non lineaire du champ intra-resonateur, tel que celui developpe. Dans un regime de forte puissance de spectroscopie, des bandes laterales apparaissent et sont causees par les fluctuations quantiques du champ electromagnetique intra-resonateur autour de sa valeur d'equilibre. Ces fluctuations sont causees par la compression du champ electromagnetique due a la non-linearite du resonateur, et l'observation de leur effet via la spectroscopie d'un qubit constitue une premiere. Suite aux succes quantitatifs du modele reduit, je montre que deux regimes de parametres ameliorent marginalement la mesure dispersive d'un qubit avec un resonateur lineaire, et significativement une mesure par bifurcation avec un resonateur non lineaire. J'explique le fonctionnement d'une mesure de qubit dans un resonateur lineaire developpee par une equipe experimentale de l'Universite de Yale. Cette mesure, qui utilise les non-linearites induites par le qubit, a une haute fidelite, mais utilise une tres haute puissance et est destructrice. Dans tous ces cas, la structure multi-niveaux du qubit s'avere cruciale pour la mesure. En suggerant des facons d'ameliorer la mesure de qubits supraconducteurs, et en decrivant quantitativement la physique d'un systeme a plusieurs niveaux couple a un resonateur non lineaire excite, les resultats presentes dans cette these sont pertinents autant pour l'utilisation de l'architecture d'electrodynamique quantique en circuit pour l'informatique quantique que pour l'optique quantique. Mots-cles: electrodynamique quantique en circuit, informatique quantique, mesure, qubit supraconducteur, transmon, non-linearite Kerr

A Novel Multi-scale Simulation Strategy for Turbulent Reacting Flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

James, Sutherland C.

In this project, a new methodology was proposed to bridge the gap between Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES). This novel methodology, titled Lattice-Based Multiscale Simulation (LBMS), creates a lattice structure of One-Dimensional Turbulence (ODT) models. This model has been shown to capture turbulent combustion with high fidelity by fully resolving interactions between turbulence and diffusion. By creating a lattice of ODT models, which are then coupled, LBMS overcomes the shortcomings of ODT, which are its inability to capture large scale three dimensional flow structures. However, by spacing these lattices significantly apart, LBMS can avoid the cursemore » of dimensionality that creates untenable computational costs associated with DNS. This project has shown that LBMS is capable of reproducing statistics of isotropic turbulent flows while coarsening the spacing between lines significantly. It also investigates and resolves issues that arise when coupling ODT lines, such as flux reconstruction perpendicular to a given ODT line, preservation of conserved quantities when eddies cross a course cell volume and boundary condition application. Robust parallelization is also investigated.« less

DOE Office of Scientific and Technical Information (OSTI.GOV)

Som, Sibendu; Wang, Zihan; Pei, Yuanjiang

A state-of-the-art spray modeling methodology, recently presented by Senecal et al. [ , , ], is applied to Large Eddy Simulations (LES) of vaporizing gasoline sprays. Simulations of non-combusting Spray G (gasoline fuel) from the Engine Combustion Network are performed. Adaptive mesh refinement (AMR) with cell sizes from 0.09 mm to 0.5 mm are utilized to further demonstrate grid convergence of the dynamic structure LES model for the gasoline sprays. Grid settings are recommended to optimize the accuracy/runtime tradeoff for LES-based spray simulations at different injection pressure conditions typically encountered in gasoline direct injection (GDI) applications. The influence of LESmore » sub-grid scale (SGS) models is explored by comparing the results from dynamic structure and Smagorinsky based models against simulations without any SGS model. Twenty different realizations are simulated by changing the random number seed used in the spray sub-models. It is shown that for global quantities such as spray penetration, comparing a single LES simulation to experimental data is reasonable. Through a detailed analysis using the relevance index (RI) criteria, recommendations are made regarding the minimum number of LES realizations required for accurate prediction of the gasoline sprays.« less

How trees uptake carbon, release water and cool themselves in air: a marriage between biophysics and turbulent fluid dynamics

NASA Astrophysics Data System (ADS)

Banerjee, Tirtha; Linn, Rodman

2017-11-01

Resolving the role of the biosphere as a terrestrial carbon sink and the nature of nonlinear couplings between carbon and water cycles across a very wide range of spatiotemporal scales constitute the scope of this work. To achieve this goal, plant physiology models are coupled with atmospheric turbulence simulations. The plant biophysics code is based on the following principles: (1) a model for photosynthesis; (2) a mass transfer model through the laminar boundary layer on leaves; (3) an optimal leaf water use strategy regulated by stomatal aperture variation; (4) a leaf-level energy balance to accommodate evaporative cooling. Leaf-level outputs are upscaled to plant, canopy and landscape scales using HIGRAD/FIRETEC, a high fidelity large eddy simulation (LES) framework developed at LANL. The coupled biophysics-CFD code can take inputs such as wind speed, light availability, ambient CO2 concentration, air temperature, site characteristics etc. and can deliver predictions for leaf temperature, transpiration, carbon assimilation, sensible and latent heat flux, which is used to illustrate the complex the complex interaction between trees and their surrounding environments. These simulation capabilities are being used to study climate feedbacks of forests and agroecosystems.

Tree-crown-resolving large-eddy simulation for evaluating greenery effects on urban heat environments

NASA Astrophysics Data System (ADS)

Matsuda, K.; Onishi, R.; Takahashi, K.

2017-12-01

Urban high temperatures due to the combined influence of global warming and urban heat islands increase the risk of heat stroke. Greenery is one of possible countermeasures for mitigating the heat environments since the transpiration and shading effect of trees can reduce the air temperature and the radiative heat flux. In order to formulate effective measures, it is important to estimate the influence of the greenery on the heat stroke risk. In this study, we have developed a tree-crown-resolving large-eddy simulation (LES) model that is coupled with three-dimensional radiative transfer (3DRT) model. The Multi-Scale Simulator for the Geoenvironment (MSSG) is used for performing building- and tree-crown-resolving LES. The 3DRT model is implemented in the MSSG so that the 3DRT is calculated repeatedly during the time integration of the LES. We have confirmed that the computational time for the 3DRT model is negligibly small compared with that for the LES and the accuracy of the 3DRT model is sufficiently high to evaluate the radiative heat flux at the pedestrian level. The present model is applied to the analysis of the heat environment in an actual urban area around the Tokyo Bay area, covering 8 km × 8 km with 5-m grid mesh, in order to confirm its feasibility. The results show that the wet-bulb globe temperature (WBGT), which is an indicator of the heat stroke risk, is predicted in a sufficiently high accuracy to evaluate the influence of tree crowns on the heat environment. In addition, by comparing with a case without the greenery in the Tokyo Bay area, we have confirmed that the greenery increases the low WBGT areas in major pedestrian spaces by a factor of 3.4. This indicates that the present model can predict the greenery effect on the urban heat environment quantitatively.

Implementation of an Online Chemistry Model to a Large Eddy Simulation Model (PALM-4U0

NASA Astrophysics Data System (ADS)

Mauder, M.; Khan, B.; Forkel, R.; Banzhaf, S.; Russo, E. E.; Sühring, M.; Kanani-Sühring, F.; Raasch, S.; Ketelsen, K.

2017-12-01

Large Eddy Simulation (LES) models permit to resolve relevant scales of turbulent motion, so that these models can capture the inherent unsteadiness of atmospheric turbulence. However, LES models are so far hardly applied for urban air quality studies, in particular chemical transformation of pollutants. In this context, BMBF (Bundesministerium für Bildung und Forschung) funded a joint project, MOSAIK (Modellbasierte Stadtplanung und Anwendung im Klimawandel / Model-based city planning and application in climate change) with the main goal to develop a new highly efficient urban climate model (UCM) that also includes atmospheric chemical processes. The state-of-the-art LES model PALM; Maronga et al, 2015, Geosci. Model Dev., 8, doi:10.5194/gmd-8-2515-2015), has been used as a core model for the new UCM named as PALM-4U. For the gas phase chemistry, a fully coupled 'online' chemistry model has been implemented into PALM. The latest version of the Kinetic PreProcessor (KPP) Version 2.3, has been utilized for the numerical integration of chemical species. Due to the high computational demands of the LES model, compromises in the description of chemical processes are required. Therefore, a reduced chemistry mechanism, which includes only major pollutants namely O3, NO, NO2, CO, a highly simplified VOC chemistry and a small number of products have been implemented. This work shows preliminary results of the advection, and chemical transformation of atmospheric pollutants. Non-cyclic boundaries have been used for inflow and outflow in east-west directions while periodic boundary conditions have been implemented to the south-north lateral boundaries. For practical applications, our approach is to go beyond the simulation of single street canyons to chemical transformation, advection and deposition of air pollutants in the larger urban canopy. Tests of chemistry schemes and initial studies of chemistry-turbulence, transport and transformations are presented.

Prediction of dynamic and mixing characteristics of drop-laden mixing layers using DNS and LES

NASA Technical Reports Server (NTRS)

Okong'o, N.; Leboissetier, A.; Bellan, J.

2004-01-01

Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) have been conducted of a temporal mixing layer laden with evaporating drops, in order to assess the ability of LES to reproduce dynamic and mixing aspects of the DNS which affect combustion, independently of combustion models.

A Model for Temperature Fluctuations in a Buoyant Plume

NASA Astrophysics Data System (ADS)

Bisignano, A.; Devenish, B. J.

2015-11-01

We present a hybrid Lagrangian stochastic model for buoyant plume rise from an isolated source that includes the effects of temperature fluctuations. The model is based on that of Webster and Thomson (Atmos Environ 36:5031-5042, 2002) in that it is a coupling of a classical plume model in a crossflow with stochastic differential equations for the vertical velocity and temperature (which are themselves coupled). The novelty lies in the addition of the latter stochastic differential equation. Parametrizations of the plume turbulence are presented that are used as inputs to the model. The root-mean-square temperature is assumed to be proportional to the difference between the centreline temperature of the plume and the ambient temperature. The constant of proportionality is tuned by comparison with equivalent statistics from large-eddy simulations (LES) of buoyant plumes in a uniform crossflow and linear stratification. We compare plume trajectories for a wide range of crossflow velocities and find that the model generally compares well with the equivalent LES results particularly when added mass is included in the model. The exception occurs when the crossflow velocity component becomes very small. Comparison of the scalar concentration, both in terms of the height of the maximum concentration and its vertical spread, shows similar behaviour. The model is extended to allow for realistic profiles of ambient wind and temperature and the results are compared with LES of the plume that emanated from the explosion and fire at the Buncefield oil depot in 2005.

Francis-99 turbine numerical flow simulation of steady state operation using RANS and RANS/LES turbulence model

NASA Astrophysics Data System (ADS)

Minakov, A.; Platonov, D.; Sentyabov, A.; Gavrilov, A.

2017-01-01

We performed numerical simulation of flow in a laboratory model of a Francis hydroturbine at three regimes, using two eddy-viscosity- (EVM) and a Reynolds stress (RSM) RANS models (realizable k-ɛ, k-ω SST, LRR) and detached-eddy-simulations (DES), as well as large-eddy simulations (LES). Comparison of calculation results with the experimental data was carried out. Unlike the linear EVMs, the RSM, DES, and LES reproduced well the mean velocity components, and pressure pulsations in the diffusor draft tube. Despite relatively coarse meshes and insufficient resolution of the near-wall region, LES, DES also reproduced well the intrinsic flow unsteadiness and the dominant flow structures and the associated pressure pulsations in the draft tube.

Etude théorique et expérimentale d'un processus thermique dans un limiteur supraconducteur de courant

NASA Astrophysics Data System (ADS)

Lévêque, J.; Netter, D.; Rezzoug, A.; Caron, J. P.

1995-12-01

The increasing of fault current level in electrical networks leads to a new interest for superconducting current limiters since 1983 when an ac wire with low losses has been developed. They are based on the natural transition from the superconducting state to the normal resistive state. In order to study the transition usual models are not sufficient. This paper deals with a numerical resolution of the electro-thermical coupled problem. Computed and experimental results are favorably compared. Some factors which affect the transition as well as the influence of the starting instant are studied. L'accroissement des courants de court-circuit dans les réseaux électriques a ravivé l'intérêt pour les limiteurs supraconducteurs de courant depuis la mise au point en 1983 de fil supraconducteur à faibles pertes en régime variable. Le principe de fonctionnement de ces limiteurs est fondé sur la transition vers l'état normal du matériau supraconducteur. Les modèles simplifiés s'avèrent insuffisants pour l'étude de la transition, une méthode de résolution numérique du problème électrothermique est donc faite. Les résultats des simulations sont comparés à ceux d'expérimentation puis deux études, l'une de facteurs influant sur la propagation de la transition, l'autre du type de court-circuit, sont menées.

From DNS to RANS: A Multi-model workflow to understand the Influence of Hurricanes on Generating Turbidity Currents in the Gulf of Mexico

NASA Astrophysics Data System (ADS)

Syvitski, J. P.; Arango, H.; Harris, C. K.; Meiburg, E. H.; Jenkins, C. J.; Auad, G.; Hutton, E.; Kniskern, T. A.; Radhakrishnan, S.

2016-12-01

A loosely coupled numerical workflow is developed to address land-sea pathways for sediment routing from terrestrial and coastal sources, across the continental shelf and ultimately down the continental slope canyon system of the northern Gulf of Mexico (GOM). Model simulations represent a range of environmental conditions that might lead to the generation of turbidity-currents. The workflow comprises: 1) A simulator for the water and sediment discharged from rivers into the GOM with WMBsedv2 with calibration using USGS and USACE gauged river data; 2) Domain grids and bathymetry (ETOPO2) for the ocean models and realistic seabed sediment texture grids (dbSEABED) for the sediment transport models; 3) A spectral wave action simulator (10 km resolution) (WaveWatch III) driven by GFDL - GFS winds; 4) A simulator for ocean dynamics (ROMS) forced with ECMWF ERA winds; 5) A simulator for seafloor resuspension and transport (CSTMS); 6) Simulators (HurriSlip) of seafloor failure and flow ignition locations for boundary input to a turbidity current model; and 7) A RANS turbidity current model (TURBINS) to route sediment flows down GOM canyons, providing estimates of bottom shear stresses. TURBINS was developed first as a DNS model and then converted to an LES model wherein a dynamic turbulence closure scheme was employed. Like most DNS to LES model comparisons (these being done by the UCSB team), turbulence scaling allowed for higher Re applications but were found still not capable of simulating field scale (GOM continental canyons) environments. The LES model was next converted to a non-hydrostatic RANS model capable of field scale applications but only with a daisy-chain approach to multiple model runs along the simulated canyon floor. These model adaptations allowed the workflow to be tested for the year 1-Oct-2007 to 30-Sep-2008 that included two domain Hurricanes (Ike and Gustav). The RANS-TURBINS employed further boundary simplifications on both sediment erosion and deposition in line with the ocean model ROMS-CSTMS.

Large-eddy simulation/Reynolds-averaged Navier-Stokes hybrid schemes for high speed flows

NASA Astrophysics Data System (ADS)

Xiao, Xudong

Three LES/RANS hybrid schemes have been proposed for the prediction of high speed separated flows. Each method couples the k-zeta (Enstrophy) BANS model with an LES subgrid scale one-equation model by using a blending function that is coordinate system independent. Two of these functions are based on turbulence dissipation length scale and grid size, while the third one has no explicit dependence on the grid. To implement the LES/RANS hybrid schemes, a new rescaling-reintroducing method is used to generate time-dependent turbulent inflow conditions. The hybrid schemes have been tested on a Mach 2.88 flow over 25 degree compression-expansion ramp and a Mach 2.79 flow over 20 degree compression ramp. A special computation procedure has been designed to prevent the separation zone from expanding upstream to the recycle-plane. The code is parallelized using Message Passing Interface (MPI) and is optimized for running on IBM-SP3 parallel machine. The scheme was validated first for a flat plate. It was shown that the blending function has to be monotonic to prevent the RANS region from appearing in the LES region. In the 25 deg ramp case, the hybrid schemes provided better agreement with experiment in the recovery region. Grid refinement studies demonstrated the importance of using a grid independent blend function and further improvement with experiment in the recovery region. In the 20 deg ramp case, with a relatively finer grid, the hybrid scheme characterized by grid independent blending function well predicted the flow field in both the separation region and the recovery region. Therefore, with "appropriately" fine grid, current hybrid schemes are promising for the simulation of shock wave/boundary layer interaction problems.

Étude des perturbations conduites et rayonnées dans une cellule de commutation

NASA Astrophysics Data System (ADS)

Costa, F.; Forest, F.; Puzo, A.; Rojat, G.

1993-12-01

The principles used in static conversion and the rise of the performances of the new switching devices contribue to increase the level of electromagnetic noises emitted by electronic converters. We have studied the way how these perturbations are created and coupled through their environment in conducted and radiated mode by a switching cell. This one can work in hard switching, zero current or voltage switching modes. We first outline the general problems of electromagnetic pollution and their metrology in converters. Then we describe the experimental environment. We analyse the mechanisms of generation of parasitic signals in a switching cell related to the electrical constraints and its switching mode. The simulated results, issued of the analytical models obtained, are confronted with the experimental ones. Then we show a method to calculate analytically the E and H near fields. It has been confirmed by experimental results. At last, we present, in a synthetic manner, the main results obtained, relative to the switching mode and the electrical constraints, using a new characterizing method. Theses results will allow the designer to incorporate the electromagnetic considerations in the conception of a converter. Les principes de commutation employés en conversion statique, l'évolution des performances statiques et dynamiques des composants, contribuent à faire des dispositifs de conversion statique de puissants générateurs de perturbations conduites et rayonnées. Nous nous sommes attachés à étudier les mécanismes de génération et de couplage des perturbations, tant en mode conduit que rayonné dans des structures à une seule cellule de commutation et fonctionnant selon les trois principaux modes de commutation : commutation forcée, à zéro de courant (ZCS), et à zéro de tension (ZVS). Après la mise en évidence de la problématique de pollution électromagnétique dans les structures et leur métrologie, nous décrivons l'environnement expérimental étudié. Nous analysons ensuite les principaux mécanismes produisant les perturbations au sein d'une cellule de commutation en introduisant un certain nombre de composants parasites. Les modèles sont simulés et confrontés aux résultats expérimentaux. Nous décrivons alors une méthode, validée expérimentalement et permettant de calculer les intensités des champs E et H proches émis. Enfin, nous présentons de façon synthétique les résultats observés selon les régimes de fonctionnement de la cellule de commutation et les contraintes électriques et dynamiques qu'elle subit. Nous avons, pour ce faire, développé une méthode originale de quantification des signaux perturbateurs. Les résultats obtenus doivent permettre d'intégrer les problèmes de pollution électromagnétique au stade de la conception d'un dispositif.

Large-Eddy Simulation: Current Capabilities, Recommended Practices, and Future Research

NASA Technical Reports Server (NTRS)

Georgiadis, Nicholas J.; Rizzetta, Donald P.; Fureby, Christer

2009-01-01

This paper presents the results of an activity by the Large Eddy Simulation (LES) Working Group of the AIAA Fluid Dynamics Technical Committee to (1) address the current capabilities of LES, (2) outline recommended practices and key considerations for using LES, and (3) identify future research needs to advance the capabilities and reliability of LES for analysis of turbulent flows. To address the current capabilities and future needs, a survey comprised of eleven questions was posed to LES Working Group members to assemble a broad range of perspectives on important topics related to LES. The responses to these survey questions are summarized with the intent not to be a comprehensive dictate on LES, but rather the perspective of one group on some important issues. A list of recommended practices is also provided, which does not treat all aspects of a LES, but provides guidance on some of the key areas that should be considered.

Large-Eddy Simulations of Dust Devils and Convective Vortices

NASA Astrophysics Data System (ADS)

Spiga, Aymeric; Barth, Erika; Gu, Zhaolin; Hoffmann, Fabian; Ito, Junshi; Jemmett-Smith, Bradley; Klose, Martina; Nishizawa, Seiya; Raasch, Siegfried; Rafkin, Scot; Takemi, Tetsuya; Tyler, Daniel; Wei, Wei

2016-11-01

In this review, we address the use of numerical computations called Large-Eddy Simulations (LES) to study dust devils, and the more general class of atmospheric phenomena they belong to (convective vortices). We describe the main elements of the LES methodology. We review the properties, statistics, and variability of dust devils and convective vortices resolved by LES in both terrestrial and Martian environments. The current challenges faced by modelers using LES for dust devils are also discussed in detail.

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

NASA Technical Reports Server (NTRS)

Balakumar, P.

2015-01-01

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Nesting large-eddy simulations within mesoscale simulations for wind energy applications

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lundquist, J K; Mirocha, J D; Chow, F K

2008-09-08

With increasing demand for more accurate atmospheric simulations for wind turbine micrositing, for operational wind power forecasting, and for more reliable turbine design, simulations of atmospheric flow with resolution of tens of meters or higher are required. These time-dependent large-eddy simulations (LES), which resolve individual atmospheric eddies on length scales smaller than turbine blades and account for complex terrain, are possible with a range of commercial and open-source software, including the Weather Research and Forecasting (WRF) model. In addition to 'local' sources of turbulence within an LES domain, changing weather conditions outside the domain can also affect flow, suggesting thatmore » a mesoscale model provide boundary conditions to the large-eddy simulations. Nesting a large-eddy simulation within a mesoscale model requires nuanced representations of turbulence. Our group has improved the Weather and Research Forecasting model's (WRF) LES capability by implementing the Nonlinear Backscatter and Anisotropy (NBA) subfilter stress model following Kosovic (1997) and an explicit filtering and reconstruction technique to compute the Resolvable Subfilter-Scale (RSFS) stresses (following Chow et al, 2005). We have also implemented an immersed boundary method (IBM) in WRF to accommodate complex terrain. These new models improve WRF's LES capabilities over complex terrain and in stable atmospheric conditions. We demonstrate approaches to nesting LES within a mesoscale simulation for farms of wind turbines in hilly regions. Results are sensitive to the nesting method, indicating that care must be taken to provide appropriate boundary conditions, and to allow adequate spin-up of turbulence in the LES domain.« less

Large eddy simulation of shock train in a convergent-divergent nozzle

NASA Astrophysics Data System (ADS)

Mousavi, Seyed Mahmood; Roohi, Ehsan

2014-12-01

This paper discusses the suitability of the Large Eddy Simulation (LES) turbulence modeling for the accurate simulation of the shock train phenomena in a convergent-divergent nozzle. To this aim, we selected an experimentally tested geometry and performed LES simulation for the same geometry. The structure and pressure recovery inside the shock train in the nozzle captured by LES model are compared with the experimental data, analytical expressions and numerical solutions obtained using various alternative turbulence models, including k-ɛ RNG, k-ω SST, and Reynolds stress model (RSM). Comparing with the experimental data, we observed that the LES solution not only predicts the "locations of the first shock" precisely, but also its results are quite accurate before and after the shock train. After validating the LES solution, we investigate the effects of the inlet total pressure on the shock train starting point and length. The effects of changes in the back pressure, nozzle inlet angle (NIA) and wall temperature on the behavior of the shock train are investigated by details.

An Immersed Boundary - Adaptive Mesh Refinement solver (IB-AMR) for high fidelity fully resolved wind turbine simulations

NASA Astrophysics Data System (ADS)

Angelidis, Dionysios; Sotiropoulos, Fotis

2015-11-01

The geometrical details of wind turbines determine the structure of the turbulence in the near and far wake and should be taken in account when performing high fidelity calculations. Multi-resolution simulations coupled with an immersed boundary method constitutes a powerful framework for high-fidelity calculations past wind farms located over complex terrains. We develop a 3D Immersed-Boundary Adaptive Mesh Refinement flow solver (IB-AMR) which enables turbine-resolving LES of wind turbines. The idea of using a hybrid staggered/non-staggered grid layout adopted in the Curvilinear Immersed Boundary Method (CURVIB) has been successfully incorporated on unstructured meshes and the fractional step method has been employed. The overall performance and robustness of the second order accurate, parallel, unstructured solver is evaluated by comparing the numerical simulations against conforming grid calculations and experimental measurements of laminar and turbulent flows over complex geometries. We also present turbine-resolving multi-scale LES considering all the details affecting the induced flow field; including the geometry of the tower, the nacelle and especially the rotor blades of a wind tunnel scale turbine. This material is based upon work supported by the Department of Energy under Award Number DE-EE0005482 and the Sandia National Laboratories.

Towards the computation of time-periodic inertial range dynamics

NASA Astrophysics Data System (ADS)

van Veen, L.; Vela-Martín, A.; Kawahara, G.

2018-04-01

We explore the possibility of computing simple invariant solutions, like travelling waves or periodic orbits, in Large Eddy Simulation (LES) on a periodic domain with constant external forcing. The absence of material boundaries and the simple forcing mechanism make this system a comparatively simple target for the study of turbulent dynamics through invariant solutions. We show, that in spite of the application of eddy viscosity the computations are still rather challenging and must be performed on GPU cards rather than conventional coupled CPUs. We investigate the onset of turbulence in this system by means of bifurcation analysis, and present a long-period, large-amplitude unstable periodic orbit that is filtered from a turbulent time series. Although this orbit is computed on a coarse grid, with only a small separation between the integral scale and the LES filter length, the periodic dynamics seem to capture a regeneration process of the large-scale vortices.

LES, DNS, and RANS for the Analysis of High-Speed Turbulent Reacting Flows

NASA Technical Reports Server (NTRS)

Colucci, P. J.; Jaberi, F. A.; Givi, P.

1996-01-01

A filtered density function (FDF) method suitable for chemically reactive flows is developed in the context of large eddy simulation. The advantage of the FDF methodology is its inherent ability to resolve subgrid scales (SGS) scalar correlations that otherwise have to be modeled. Because of the lack of robust models to accurately predict these correlations in turbulent reactive flows, simulations involving turbulent combustion are often met with a degree of skepticism. The FDF methodology avoids the closure problem associated with these terms and treats the reaction in an exact manner. The scalar FDF approach is particularly attractive since it can be coupled with existing hydrodynamic computational fluid dynamics (CFD) codes.

Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows

NASA Technical Reports Server (NTRS)

Givi, Peyman; Madnia, Cyrus K.; Steinberger, C. J.; Frankel, S. H.

1992-01-01

The principal objective is to extend the boundaries within which large eddy simulations (LES) and direct numerical simulations (DNS) can be applied in computational analyses of high speed reacting flows. A summary of work accomplished during the last six months is presented.

Development of a High Level Architecture Federation of Ship Replenishment at Sea

DTIC Science & Technology

2011-10-01

utiliser une infrastructure de simulation appelée architecture de haut niveau (HLA) afin de fournir des environne - ments de simulation interarmées...fournir un environnement de simulation qui modélise l’interactions entre les divers composants afin de simuler les conditions qui mènent aux

Radiation-Spray Coupling for Realistic Flow Configurations

NASA Technical Reports Server (NTRS)

El-Asrag, Hossam; Iannetti, Anthony C.

2011-01-01

Three Large Eddy Simulations (LES) for a lean-direct injection (LDI) combustor are performed and compared. In addition to the cold flow simulation, the effect of radiation coupling with the multi-physics reactive flow is analyzed. The flame let progress variable approach is used as a subgrid combustion model combined with a stochastic subgrid model for spray atomization and an optically thin radiation model. For accurate chemistry modeling, a detailed Jet-A surrogate mechanism is utilized. To achieve realistic inflow, a simple recycling technique is performed at the inflow section upstream of the swirler. Good comparison is shown with the experimental data mean and root mean square profiles. The effect of combustion is found to change the shape and size of the central recirculation zone. Radiation is found to change the spray dynamics and atomization by changing the heat release distribution and the local temperature values impacting the evaporation process. The simulation with radiation modeling shows wider range of droplet size distribution by altering the evaporation rate. The current study proves the importance of radiation modeling for accurate prediction in realistic spray combustion configurations, even for low pressure systems.

Development of a Hybrid RANS/LES Method for Compressible Mixing Layer Simulations

NASA Technical Reports Server (NTRS)

Georgiadis, Nicholas J.; Alexander, J. Iwan D.; Reshotko, Eli

2001-01-01

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modem day aircraft and also those of hypersonic vehicles currently under development. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS/LES method on stretched, non-Cartesian grids. The hybrid RANS/LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two-dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Actual LES calculations, performed in three spatial directions, indicated an initial vortex shedding followed by rapid transition to turbulence, which is in agreement with experimental observations.

LES models for incompressible magnetohydrodynamics derived from the variational multiscale formulation

NASA Astrophysics Data System (ADS)

Sondak, David; Oberai, Assad

2012-10-01

Novel large eddy simulation (LES) models are developed for incompressible magnetohydrodynamics (MHD). These models include the application of the variational multiscale formulation (VMS) of LES to the equations of incompressible MHD, a new residual-based eddy viscosity model (RBEVM,) and a mixed LES model that combines the strengths of both of these models. The new models result in a consistent numerical method that is relatively simple to implement. A dynamic procedure for determining model coefficients is no longer required. The new LES models are tested on a decaying Taylor-Green vortex generalized to MHD and benchmarked against classical and state-of-the art LES turbulence models as well as direct numerical simulations (DNS). These new models are able to account for the essential MHD physics which is demonstrated via comparisons of energy spectra. We also compare the performance of our models to a DNS simulation by A. Pouquet et al., for which the ratio of DNS modes to LES modes is 262,144. Additionally, we extend these models to a finite element setting in which boundary conditions play a role. A classic problem on which we test these models is turbulent channel flow, which in the case of MHD, is called Hartmann flow.

Simulations of turbulent rotating flows using a subfilter scale stress model derived from the partially integrated transport modeling method

NASA Astrophysics Data System (ADS)

Chaouat, Bruno

2012-04-01

The partially integrated transport modeling (PITM) method [B. Chaouat and R. Schiestel, "A new partially integrated transport model for subgrid-scale stresses and dissipation rate for turbulent developing flows," Phys. Fluids 17, 065106 (2005), 10.1063/1.1928607; R. Schiestel and A. Dejoan, "Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations," Theor. Comput. Fluid Dyn. 18, 443 (2005), 10.1007/s00162-004-0155-z; B. Chaouat and R. Schiestel, "From single-scale turbulence models to multiple-scale and subgridscale models by Fourier transform," Theor. Comput. Fluid Dyn. 21, 201 (2007), 10.1007/s00162-007-0044-3; B. Chaouat and R. Schiestel, "Progress in subgrid-scale transport modelling for continuous hybrid non-zonal RANS/LES simulations," Int. J. Heat Fluid Flow 30, 602 (2009), 10.1016/j.ijheatfluidflow.2009.02.021] viewed as a continuous approach for hybrid RANS/LES (Reynolds averaged Navier-Stoke equations/large eddy simulations) simulations with seamless coupling between RANS and LES regions is used to derive a subfilter scale stress model in the framework of second-moment closure applicable in a rotating frame of reference. This present subfilter scale model is based on the transport equations for the subfilter stresses and the dissipation rate and appears well appropriate for simulating unsteady flows on relatively coarse grids or flows with strong departure from spectral equilibrium because the cutoff wave number can be located almost anywhere inside the spectrum energy. According to the spectral theory developed in the wave number space [B. Chaouat and R. Schiestel, "From single-scale turbulence models to multiple-scale and subgrid-scale models by Fourier transform," Theor. Comput. Fluid Dyn. 21, 201 (2007), 10.1007/s00162-007-0044-3], the coefficients used in this model are no longer constants but they are some analytical functions of a dimensionless parameter controlling the spectral distribution of turbulence. The pressure-strain correlation term encompassed in this model is inspired from the nonlinear SSG model [C. G. Speziale, S. Sarkar, and T. B. Gatski, "Modelling the pressure-strain correlation of turbulence: an invariant dynamical systems approach," J. Fluid Mech. 227, 245 (1991), 10.1017/S0022112091000101] developed initially for homogeneous rotating flows in RANS methodology. It is modeled in system rotation using the principle of objectivity. Its modeling is especially extended in a low Reynolds number version for handling non-homogeneous wall flows. The present subfilter scale stress model is then used for simulating large scales of rotating turbulent flows on coarse and medium grids at moderate, medium, and high rotation rates. It is also applied to perform a simulation on a refined grid at the highest rotation rate. As a result, it is found that the PITM simulations reproduce fairly well the mean features of rotating channel flows allowing a drastic reduction of the computational cost in comparison with the one required for performing highly resolved LES. Overall, the mean velocities and turbulent stresses are found to be in good agreement with the data of highly resolved LES [E. Lamballais, O. Metais, and M. Lesieur, "Spectral-dynamic model for large-eddy simulations of turbulent rotating flow," Theor. Comput. Fluid Dyn. 12, 149 (1998)]. The anisotropy character of the flow resulting from the rotation effects is also well reproduced in accordance with the reference data. Moreover, the PITM2 simulations performed on the medium grid predict qualitatively well the three-dimensional flow structures as well as the longitudinal roll cells which appear in the anticyclonic wall-region of the rotating flows. As expected, the PITM3 simulation performed on the refined grid reverts to highly resolved LES. The present model based on a rational formulation appears to be an interesting candidate for tackling a large variety of engineering flows subjected to rotation.

Multi-Sensor Systems and Data Fusion for Telecommunications, Remote Sensing and Radar (les Systemes multi-senseurs et le fusionnement des donnees pour les telecommunications, la teledetection et les radars)

DTIC Science & Technology

1998-04-01

The result of the project is a demonstration of the fusion process, the sensors management and the real-time capabilities using simulated sensors...demonstrator (TAD) is a system that demonstrates the core ele- ment of a battlefield ground surveillance system by simulation in near real-time. The core...Management and Sensor/Platform simulation . The surveillance system observes the real world through a non-collocated heterogene- ous multisensory system

Evaluation of probabilistic flow in two unsaturated soils

NASA Astrophysics Data System (ADS)

Boateng, Samuel

2001-11-01

A variably saturated flow model is coupled to a first-order reliability algorithm to simulate unsaturated flow in two soils. The unsaturated soil properties are considered as uncertain variables with means, standard deviations, and marginal probability distributions. Thus, each simulation constitutes an unsaturated probability flow event. Sensitivities of the uncertain variables are estimated for each event. The unsaturated hydraulic properties of a fine-textured soil and a coarse-textured soil are used. The properties are based on the van Genuchten model. The flow domain has a recharge surface, a seepage boundary along the bottom, and a no-flow boundary along the sides. The uncertain variables are saturated water content, residual water content, van Genuchten model parameters alpha (α) and n, and saturated hydraulic conductivity. The objective is to evaluate the significance of each uncertain variable to the probabilistic flow. Under wet conditions, saturated water content and residual water content are the most significant uncertain variables in the sand. For dry conditions in the sand, however, the van Genuchten model parameters α and n are the most significant. Model parameter n and saturated hydraulic conductivity are the most significant for the wet clay loam. Saturated water content is most significant for the dry clay loam. Résumé. Un modèle d'écoulement variable en milieu saturé est couplé à un algorithme d'exactitude de premier ordre pour simuler les écoulements en milieu non saturé dans deux sols. Les propriétés des sols non saturés sont considérés comme des variables incertaines avec des moyennes, des écarts-types et des distributions de probabilité marginale. Ainsi chaque simulation constitue un événement d'écoulement non saturé probable. La sensibilité des variables incertaines est estimée pour chaque événement. Les propriétés hydrauliques non saturées d'un sol à texture fine et d'un sol à texture grossière sont utilisées. Les propriétés sont basées sur le modèle de van Genuchten. Le domaine d'écoulement possède une surface de recharge, une limite de fuite à sa base et des limites sans écoulement sur les côtés. Les variables incertaines sont la teneur en eau à saturation, la teneur en eau résiduelle, les paramètres alpha (α) et n du modèle de van Genuchten et la conductivité hydraulique à saturation. L'objectif est d'évaluer la signification de chacune des variables incertaines dans l'écoulement probabiliste. Dans des conditions humides, la teneur en eau à saturation et la teneur en eau résiduelle sont les variables incertaines les plus significatives dans le sable. Toutefois, dans des conditions sèches dans le sable, les paramètres α et n du modèle de van Genuchten sont les plus significatifs. Le paramètre n du modèle et la conductivité hydraulique à saturation sont les plus significatifs pour un sol argileux humide. La teneur en eau à saturation est très significative pour le sol argileux sec. Resumen. Se ha acoplado un modelo de flujo de saturación variable con un algoritmo de fiabilidad de primer orden con el fin de simular el flujo no saturado en dos tipos de suelos. Se ha tratado las propiedades del suelo no saturado como variables inciertas, a las que se asigna las medias, desviaciones estándar y distribuciones de probabilidad marginal correspondientes. Así, cada simulación constituye un evento probabilístico de flujo no saturado y la sensibilidad de las variables inciertas es estimada para cada evento. Se ha utilizado las propiedades de la conductividad hidráulica no saturada de dos suelos con dos tipos de textura - fina y gruesa - mediante el modelo de van Genuchten. El dominio de flujo está delimitado por una superficie de recarga, base de goteo y contornos laterales de flujo nulo. Las variables inciertas son el contenido de agua residual, el de saturación, los parámetros del modelo de van Genuchten (α y n) y la conductividad hidráulica saturada. El objetivo era evaluar la contribución de cada variable incierta al flujo probabilístico. Para arenas, las variables inciertas más importantes, en condiciones de humedad, son el contenido de agua residual y el de saturación en ausencia de humedad, lo son ambos parámetros del modelo de van Genuchten. Para margas arcillosas, las variables más significativas en condiciones húmedas son el parámetro n y la conductividad hidráulica saturada; en condiciones secas, el contenido de agua en saturación.

LES of a Jet Excited by the Localized Arc Filament Plasma Actuators

NASA Technical Reports Server (NTRS)

Brown, Clifford A.

2011-01-01

The fluid dynamics of a high-speed jet are governed by the instability waves that form in the free-shear boundary layer of the jet. Jet excitation manipulates the growth and saturation of particular instability waves to control the unsteady flow structures that characterize the energy cascade in the jet.The results may include jet noise mitigation or a reduction in the infrared signature of the jet. The Localized Arc Filament Plasma Actuators (LAFPA) have demonstrated the ability to excite a high-speed jets in laboratory experiments. Extending and optimizing this excitation technology, however, is a complex process that will require many tests and trials. Computational simulations can play an important role in understanding and optimizing this actuator technology for real-world applications. Previous research has focused on developing a suitable actuator model and coupling it with the appropriate computational fluid dynamics (CFD) methods using two-dimensional spatial flow approximations. This work is now extended to three-dimensions (3-D) in space. The actuator model is adapted to a series of discrete actuators and a 3-D LES simulation of an excited jet is run. The results are used to study the fluid dynamics near the actuator and in the jet plume.

Simulation d'ecoulements internes compressibles laminaires et turbulents par une methode d'elements finis

NASA Astrophysics Data System (ADS)

Rebaine, Ali

1997-08-01

Ce travail consiste en la simulation numerique des ecoulements internes compressibles bidimensionnels laminaires et turbulents. On s'interesse, particulierement, aux ecoulements dans les ejecteurs supersoniques. Les equations de Navier-Stokes sont formulees sous forme conservative et utilisent, comme variables independantes, les variables dites enthalpiques a savoir: la pression statique, la quantite de mouvement et l'enthalpie totale specifique. Une formulation variationnelle stable des equations de Navier-Stokes est utilisee. Elle est base sur la methode SUPG (Streamline Upwinding Petrov Galerkin) et utilise un operateur de capture des forts gradients. Un modele de turbulence, pour la simulation des ecoulements dans les ejecteurs, est mis au point. Il consiste a separer deux regions distinctes: une region proche de la paroi solide, ou le modele de Baldwin et Lomax est utilise et l'autre, loin de la paroi, ou une formulation nouvelle, basee sur le modele de Schlichting pour les jets, est proposee. Une technique de calcul de la viscosite turbulente, sur un maillage non structure, est implementee. La discretisation dans l'espace de la forme variationnelle est faite a l'aide de la methode des elements finis en utilisant une approximation mixte: quadratique pour les composantes de la quantite de mouvement et de la vitesse et lineaire pour le reste des variables. La discretisation temporelle est effectuee par une methode de differences finies en utilisant le schema d'Euler implicite. Le systeme matriciel, resultant de la discretisation spatio-temporelle, est resolu a l'aide de l'algorithme GMRES en utilisant un preconditionneur diagonal. Les validations numeriques ont ete menees sur plusieurs types de tuyeres et ejecteurs. La principale validation consiste en la simulation de l'ecoulement dans l'ejecteur teste au centre de recherche NASA Lewis. Les resultats obtenus sont tres comparables avec ceux des travaux anterieurs et sont nettement superieurs concernant les ecoulements turbulents dans les ejecteurs.

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

NASA Technical Reports Server (NTRS)

Raman, Venkatramanan; Pitsch, Heinz; Fox, Rodney O.

2003-01-01

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Exploring the large-scale structure of Taylor–Couette turbulence through Large-Eddy Simulations

NASA Astrophysics Data System (ADS)

Ostilla-Mónico, Rodolfo; Zhu, Xiaojue; Verzicco, Roberto

2018-04-01

Large eddy simulations (LES) of Taylor-Couette (TC) flow, the flow between two co-axial and independently rotating cylinders are performed in an attempt to explore the large-scale axially-pinned structures seen in experiments and simulations. Both static and dynamic LES models are used. The Reynolds number is kept fixed at Re = 3.4 · 104, and the radius ratio η = ri /ro is set to η = 0.909, limiting the effects of curvature and resulting in frictional Reynolds numbers of around Re τ ≈ 500. Four rotation ratios from Rot = ‑0.0909 to Rot = 0.3 are simulated. First, the LES of TC is benchmarked for different rotation ratios. Both the Smagorinsky model with a constant of cs = 0.1 and the dynamic model are found to produce reasonable results for no mean rotation and cyclonic rotation, but deviations increase for increasing rotation. This is attributed to the increasing anisotropic character of the fluctuations. Second, “over-damped” LES, i.e. LES with a large Smagorinsky constant is performed and is shown to reproduce some features of the large-scale structures, even when the near-wall region is not adequately modeled. This shows the potential for using over-damped LES for fast explorations of the parameter space where large-scale structures are found.

Improving Human Effectiveness Through Embedded Virtual Simulation (Amelioration de l’efficacite humaine grace a la simulation virtuelle integree)

DTIC Science & Technology

2014-01-01

mesurer leur expérience. Le groupe a rencontré des experts des technologies associées : les environnements virtuels, la réalité augmentée, les agents...virtuels, l’entraînement, l’ergonomie et la performance humaine. Les réflexions et conclusions issues de ces réunions et discussions sont résumées dans...en matière de EVS/ET (incluant les caractéristiques de l’utilisateur, la mission et l’environnement) et la gestion de l’entraînement est un

Calibration and Forward Uncertainty Propagation for Large-eddy Simulations of Engineering Flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Templeton, Jeremy Alan; Blaylock, Myra L.; Domino, Stefan P.

2015-09-01

The objective of this work is to investigate the efficacy of using calibration strategies from Uncertainty Quantification (UQ) to determine model coefficients for LES. As the target methods are for engineering LES, uncertainty from numerical aspects of the model must also be quantified. 15 The ultimate goal of this research thread is to generate a cost versus accuracy curve for LES such that the cost could be minimized given an accuracy prescribed by an engineering need. Realization of this goal would enable LES to serve as a predictive simulation tool within the engineering design process.

Description of the LASSO Alpha 1 Release

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gustafson, William I.; Vogelmann, Andrew M.; Cheng, Xiaoping

The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility began a pilot project in May 2015 to design a routine, high-resolution modeling capability to complement ARM’s extensive suite of measurements. This modeling capability has been named the Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) project. The availability of LES simulations with concurrent observations will serve many purposes. LES helps bridge the scale gap between DOE ARM observations and models, and the use of routine LES adds value to observations. It provides a self-consistent representation of the atmosphere and a dynamical context for the observations. Further,more » it elucidates unobservable processes and properties. LASSO will generate a simulation library for researchers that enables statistical approaches beyond a single-case mentality. It will also provide tools necessary for modelers to reproduce the LES and conduct their own sensitivity experiments. Many different uses are envisioned for the combined LASSO LES and observational library. For an observationalist, LASSO can help inform instrument remote-sensing retrievals, conduct Observation System Simulation Experiments (OSSEs), and test implications of radar scan strategies or flight paths. For a theoretician, LASSO will help calculate estimates of fluxes and co-variability of values, and test relationships without having to run the model yourself. For a modeler, LASSO will help one know ahead of time which days have good forcing, have co-registered observations at high-resolution scales, and have simulation inputs and corresponding outputs to test parameterizations. Further details on the overall LASSO project are available at http://www.arm. gov/science/themes/lasso.« less

Sensitivity of Shallow Convection in Large-Eddy Simulations to Forcing Datasets Across a Range of Days: Examining Results from the DOE LASSO Projec

NASA Astrophysics Data System (ADS)

Gustafson, W. I., Jr.; Vogelmann, A. M.; Li, Z.; Cheng, X.; Endo, S.; Krishna, B.; Toto, T.; Xiao, H.

2017-12-01

Large-eddy simulation (LES) is a powerful tool for understanding atmospheric turbulence and cloud development. However, the results are sensitive to the choice of forcing data sets used to drive the LES model, and the most realistic forcing data is difficult to identify a priori. Knowing the sensitivity of boundary layer and cloud processes to forcing data selection is critical when using LES to understand atmospheric processes and when developing associated parameterizations. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility has been developing the capability to routinely generate ensembles of LES based on a selection of plausible input forcing data sets. The LES ARM Symbiotic Simulation and Observation (LASSO) project is initially generating simulations for shallow convection days at the ARM Southern Great Plains site in Oklahoma. This talk will examine 13 days with shallow convection selected from the period May-August 2016, with multiple forcing sources and spatial scales used to generate an LES ensemble for each of the days, resulting in hundreds of LES runs with coincident observations from ARM's extensive suite of in situ and retrieval-based products. This talk will focus particularly on the sensitivity of the cloud development and its relation to forcing data. Variability of the PBL characteristics, lifting condensation level, cloud base height, cloud fraction, and liquid water path will be examined. More information about the LASSO project can be found at https://www.arm.gov/capabilities/modeling/lasso.

Large Eddy Simulation of Gravitational Effects on Transitional and Turbulent Gas-Jet Diffusion Flames

NASA Technical Reports Server (NTRS)

Givi, Peyman; Jaberi, Farhad A.

2001-01-01

The basic objective of this work is to assess the influence of gravity on "the compositional and the spatial structures" of transitional and turbulent diffusion flames via large eddy simulation (LES), and direct numerical simulation (DNS). The DNS is conducted for appraisal of the various closures employed in LES, and to study the effect of buoyancy on the small scale flow features. The LES is based on our "filtered mass density function"' (FMDF) model. The novelty of the methodology is that it allows for reliable simulations with inclusion of "realistic physics." It also allows for detailed analysis of the unsteady large scale flow evolution and compositional flame structure which is not usually possible via Reynolds averaged simulations.

The cause and resolution of log-layer mismatch in wall-modeled LES: a new perspective and its implication in complex flows

NASA Astrophysics Data System (ADS)

Park, George; Yang, Xiang; Moin, Parviz

2017-11-01

Log-layer mismatch (LLM) refers to the erroneous shifts of the mean velocity profile in the log-law region when wall models are coupled to the LES solution at the first off-wall grid points. It is often believed that the discretization error and subgrid-scale modeling error in the highly under resolved near-wall region contaminates the first off-wall LES solution, thereby providing inaccurate input to wall models resulting in inaccurate wall shear stress. Placing the LES/wall-model interface a couple of cells away from the wall has been recommended to avoid LLM. However, its non-local nature render this method impractical for flows involving complex geometry, by incurring significant overhead in LES mesh preparation and wall-model implementation. We propose an alternative remedy for LLM which warrants the removal of LLM while utilizing the first off-wall LES data. The method is based on filtering the wall-model input either in space or in time. It is simple, easy to implement, and would be particularly well suited for unstructured-grid LES involving complex geometries. We also demonstrate that LLM is caused by excessive correlation between the wall-model input and its wall shear stress output. This research is sponsored by NASA (NNX15AU93A) and ONR (FA9550-16-1-0319).

Film cooling from inclined cylindrical holes using large eddy simulations

NASA Astrophysics Data System (ADS)

Peet, Yulia V.

2006-12-01

The goal of the present study is to investigate numerically the physics of the flow, which occurs during the film cooling from inclined cylindrical holes, Film cooling is a technique used in gas turbine industry to reduce heat fluxes to the turbine blade surface. Large Eddy Simulation (LES) is performed modeling a realistic film cooling configuration, which consists of a large stagnation-type reservoir, feeding an array of discrete cooling holes (film holes) flowing into a flat plate turbulent boundary layer. Special computational methodology is developed for this problem, involving coupled simulations using multiple computational codes. A fully compressible LES code is used in the area above the flat plate, while a low Mach number LES code is employed in the plenum and film holes. The motivation for using different codes comes from the essential difference in the nature of the flow in these different regions. Flowfield is analyzed inside the plenum, film hole and a crossflow region. Flow inside the plenum is stagnating, except for the region close to the exit, where it accelerates rapidly to turn into the hole. The sharp radius of turning at the trailing edge of the plenum pipe connection causes the flow to separate from the downstream wall of the film hole. After coolant injection occurs, a complex flowfield is formed consisting of coherent vortical structures responsible for bringing hot crossflow fluid in contact with the walls of either the film hole or the blade, thus reducing cooling protection. Mean velocity and turbulent statistics are compared to experimental measurements, yielding good agreement for the mean flowfield and satisfactory agreement for the turbulence quantities. LES results are used to assess the applicability of basic assumptions of conventional eddy viscosity turbulence models used with Reynolds-averaged (RANS) approach, namely the isotropy of an eddy viscosity and thermal diffusivity. It is shown here that these assumptions do not hold for the film cooling flows. Comparison of film cooling effectiveness with experiments shows fair agreement for the centerline and laterally-averaged effectiveness. Lateral growth of the jet as judged from the lateral distribution of effectiveness is predicted correctly.

LES versus DNS: A comparative study

NASA Technical Reports Server (NTRS)

Shtilman, L.; Chasnov, J. R.

1992-01-01

We have performed Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of forced isotropic turbulence at moderate Reynolds numbers. The subgrid scale model used in the LES is based on an eddy viscosity which adjusts instantaneously the energy spectrum of the LES to that of the DNS. The statistics of the large scales of the DNS (filtered DNS field or fDNS) are compared to that of the LES. We present results for the transfer spectra, the skewness and flatness factors of the velocity components, the PDF's of the angle between the vorticity and the eigenvectors of the rate of strain, and that between the vorticity and the vorticity stretching tensor. The above LES statistics are found to be in good agreement with those measured in the fDNS field. We further observe that in all the numerical measurements, the trend was for the LES field to be more gaussian than the fDNS field. Future research on this point is planned.

Large eddy simulation for predicting turbulent heat transfer in gas turbines

PubMed Central

Tafti, Danesh K.; He, Long; Nagendra, K.

2014-01-01

Blade cooling technology will play a critical role in the next generation of propulsion and power generation gas turbines. Accurate prediction of blade metal temperature can avoid the use of excessive compressed bypass air and allow higher turbine inlet temperature, increasing fuel efficiency and decreasing emissions. Large eddy simulation (LES) has been established to predict heat transfer coefficients with good accuracy under various non-canonical flows, but is still limited to relatively simple geometries and low Reynolds numbers. It is envisioned that the projected increase in computational power combined with a drop in price-to-performance ratio will make system-level simulations using LES in complex blade geometries at engine conditions accessible to the design process in the coming one to two decades. In making this possible, two key challenges are addressed in this paper: working with complex intricate blade geometries and simulating high-Reynolds-number (Re) flows. It is proposed to use the immersed boundary method (IBM) combined with LES wall functions. A ribbed duct at Re=20 000 is simulated using the IBM, and a two-pass ribbed duct is simulated at Re=100 000 with and without rotation (rotation number Ro=0.2) using LES with wall functions. The results validate that the IBM is a viable alternative to body-conforming grids and that LES with wall functions reproduces experimental results at a much lower computational cost. PMID:25024418

Unstructured LES of Reacting Multiphase Flows in Realistic Gas Turbine Combustors

NASA Technical Reports Server (NTRS)

Ham, Frank; Apte, Sourabh; Iaccarino, Gianluca; Wu, Xiao-Hua; Herrmann, Marcus; Constantinescu, George; Mahesh, Krishnan; Moin, Parviz

2003-01-01

As part of the Accelerated Strategic Computing Initiative (ASCI) program, an accurate and robust simulation tool is being developed to perform high-fidelity LES studies of multiphase, multiscale turbulent reacting flows in aircraft gas turbine combustor configurations using hybrid unstructured grids. In the combustor, pressurized gas from the upstream compressor is reacted with atomized liquid fuel to produce the combustion products that drive the downstream turbine. The Large Eddy Simulation (LES) approach is used to simulate the combustor because of its demonstrated superiority over RANS in predicting turbulent mixing, which is central to combustion. This paper summarizes the accomplishments of the combustor group over the past year, concentrating mainly on the two major milestones achieved this year: 1) Large scale simulation: A major rewrite and redesign of the flagship unstructured LES code has allowed the group to perform large eddy simulations of the complete combustor geometry (all 18 injectors) with over 100 million control volumes; 2) Multi-physics simulation in complex geometry: The first multi-physics simulations including fuel spray breakup, coalescence, evaporation, and combustion are now being performed in a single periodic sector (1/18th) of an actual Pratt & Whitney combustor geometry.

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

NASA Technical Reports Server (NTRS)

Bui, Trong T.

1999-01-01

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Large Eddy Simulation of Engineering Flows: A Bill Reynolds Legacy.

NASA Astrophysics Data System (ADS)

Moin, Parviz

2004-11-01

The term, Large eddy simulation, LES, was coined by Bill Reynolds, thirty years ago when he and his colleagues pioneered the introduction of LES in the engineering community. Bill's legacy in LES features his insistence on having a proper mathematical definition of the large scale field independent of the numerical method used, and his vision for using numerical simulation output as data for research in turbulence physics and modeling, just as one would think of using experimental data. However, as an engineer, Bill was pre-dominantly interested in the predictive capability of computational fluid dynamics and in particular LES. In this talk I will present the state of the art in large eddy simulation of complex engineering flows. Most of this technology has been developed in the Department of Energy's ASCI Program at Stanford which was led by Bill in the last years of his distinguished career. At the core of this technology is a fully implicit non-dissipative LES code which uses unstructured grids with arbitrary elements. A hybrid Eulerian/ Largangian approach is used for multi-phase flows, and chemical reactions are introduced through dynamic equations for mixture fraction and reaction progress variable in conjunction with flamelet tables. The predictive capability of LES is demonstrated in several validation studies in flows with complex physics and complex geometry including flow in the combustor of a modern aircraft engine. LES in such a complex application is only possible through efficient utilization of modern parallel super-computers which was recognized and emphasized by Bill from the beginning. The presentation will include a brief mention of computer science efforts for efficient implementation of LES.

The Aliso Canyon Natural Gas Leak : Large Eddy Simulations for Modeling Atmospheric Dynamics and Interpretation of Observations.

NASA Astrophysics Data System (ADS)

Prasad, K.; Thorpe, A. K.; Duren, R. M.; Thompson, D. R.; Whetstone, J. R.

2016-12-01

The National Institute of Standards and Technology (NIST) has supported the development and demonstration of a measurement capability to accurately locate greenhouse gas sources and measure their flux to the atmosphere over urban domains. However, uncertainties in transport models which form the basis of all top-down approaches can significantly affect our capability to attribute sources and predict their flux to the atmosphere. Reducing uncertainties between bottom-up and top-down models will require high resolution transport models as well as validation and verification of dispersion models over an urban domain. Tracer experiments involving the release of Perfluorocarbon Tracers (PFTs) at known flow rates offer the best approach for validating dispersion / transport models. However, tracer experiments are limited by cost, ability to make continuous measurements, and environmental concerns. Natural tracer experiments, such as the leak from the Aliso Canyon underground storage facility offers a unique opportunity to improve and validate high resolution transport models, test leak hypothesis, and to estimate the amount of methane released.High spatial resolution (10 m) Large Eddy Simulations (LES) coupled with WRF atmospheric transport models were performed to simulate the dynamics of the Aliso Canyon methane plume and to quantify the source. High resolution forward simulation results were combined with aircraft and tower based in-situ measurements as well as data from NASA airborne imaging spectrometers. Comparison of simulation results with measurement data demonstrate the capability of the LES models to accurately model transport and dispersion of methane plumes over urban domains.

Numerical simulation of turbulent combustion: Scientific challenges

NASA Astrophysics Data System (ADS)

Ren, ZhuYin; Lu, Zhen; Hou, LingYun; Lu, LiuYan

2014-08-01

Predictive simulation of engine combustion is key to understanding the underlying complicated physicochemical processes, improving engine performance, and reducing pollutant emissions. Critical issues as turbulence modeling, turbulence-chemistry interaction, and accommodation of detailed chemical kinetics in complex flows remain challenging and essential for high-fidelity combustion simulation. This paper reviews the current status of the state-of-the-art large eddy simulation (LES)/prob-ability density function (PDF)/detailed chemistry approach that can address the three challenging modelling issues. PDF as a subgrid model for LES is formulated and the hybrid mesh-particle method for LES/PDF simulations is described. Then the development need in micro-mixing models for the PDF simulations of turbulent premixed combustion is identified. Finally the different acceleration methods for detailed chemistry are reviewed and a combined strategy is proposed for further development.

Comparison of Measured and WRF-LES Turbulence Statistics in a Real Convective Boundary Layer over Complex Terrain

NASA Astrophysics Data System (ADS)

Rai, R. K.; Berg, L. K.; Kosovic, B.; Mirocha, J. D.; Pekour, M. S.; Shaw, W. J.

2015-12-01

Resolving the finest turbulent scales present in the lower atmosphere using numerical simulations helps to study the processes that occur in the atmospheric boundary layer, such as the turbulent inflow condition to the wind plant and the generation of the wake behind wind turbines. This work employs several nested domains in the WRF-LES framework to simulate conditions in a convectively driven cloud free boundary layer at an instrumented field site in complex terrain. The innermost LES domain (30 m spatial resolution) receives the boundary forcing from two other coarser resolution LES outer domains, which in turn receive boundary conditions from two WRF-mesoscale domains. Wind and temperature records from sonic anemometers mounted at two vertical levels (30 m and 60 m) are compared with the LES results in term of first and second statistical moments as well as power spectra and distributions of wind velocity. For the two mostly used boundary layer parameterizations (MYNN and YSU) tested in the WRF mesoscale domains, the MYNN scheme shows slightly better agreement with the observations for some quantities, such as time averaged velocity and Turbulent Kinetic Energy (TKE). However, LES driven by WRF-mesoscale simulations using either parameterization have similar velocity spectra and distributions of velocity. For each component of the wind velocity, WRF-LES power spectra are found to be comparable to the spectra derived from the measured data (for the frequencies that are accurately represented by WRF-LES). Furthermore, the analysis of LES results shows a noticeable variability of the mean and variance even over small horizontal distances that would be considered sub-grid scale in mesoscale simulations. This observed statistical variability in space and time can be utilized to further analyze the turbulence quantities over a heterogeneous surface and to improve the turbulence parameterization in the mesoscale model.

Assessment of stretched vortex subgrid-scale models for LES of incompressible inhomogeneous turbulent flow

PubMed Central

Shetty, Dinesh A.; Frankel, Steven H.

2013-01-01

Summary The physical space version of the stretched vortex subgrid scale model [Phys. Fluids 12, 1810 (2000)] is tested in large eddy simulations (LES) of the turbulent lid driven cubic cavity flow. LES is carried out using a higher order finite-difference method [J. Comput. Phys. 229, 8802 (2010)]. The effects of different vortex orientation models and subgrid turbulence spectrums are assessed through comparisons of the LES predictions against direct numerical simulations (DNS) [Phys. Fluids 12, 1363 (2000)]. Three Reynolds numbers 12000, 18000, and 22000 are studied. Good agreement with the DNS data for the mean and fluctuating quantities is observed. PMID:24187423

Hybrid RANS-LES using high order numerical methods

NASA Astrophysics Data System (ADS)

Henry de Frahan, Marc; Yellapantula, Shashank; Vijayakumar, Ganesh; Knaus, Robert; Sprague, Michael

2017-11-01

Understanding the impact of wind turbine wake dynamics on downstream turbines is particularly important for the design of efficient wind farms. Due to their tractable computational cost, hybrid RANS/LES models are an attractive framework for simulating separation flows such as the wake dynamics behind a wind turbine. High-order numerical methods can be computationally efficient and provide increased accuracy in simulating complex flows. In the context of LES, high-order numerical methods have shown some success in predictions of turbulent flows. However, the specifics of hybrid RANS-LES models, including the transition region between both modeling frameworks, pose unique challenges for high-order numerical methods. In this work, we study the effect of increasing the order of accuracy of the numerical scheme in simulations of canonical turbulent flows using RANS, LES, and hybrid RANS-LES models. We describe the interactions between filtering, model transition, and order of accuracy and their effect on turbulence quantities such as kinetic energy spectra, boundary layer evolution, and dissipation rate. This work was funded by the U.S. Department of Energy, Exascale Computing Project, under Contract No. DE-AC36-08-GO28308 with the National Renewable Energy Laboratory.

Regularization method for large eddy simulations of shock-turbulence interactions

NASA Astrophysics Data System (ADS)

Braun, N. O.; Pullin, D. I.; Meiron, D. I.

2018-05-01

The rapid change in scales over a shock has the potential to introduce unique difficulties in Large Eddy Simulations (LES) of compressible shock-turbulence flows if the governing model does not sufficiently capture the spectral distribution of energy in the upstream turbulence. A method for the regularization of LES of shock-turbulence interactions is presented which is constructed to enforce that the energy content in the highest resolved wavenumbers decays as k - 5 / 3, and is computed locally in physical-space at low computational cost. The application of the regularization to an existing subgrid scale model is shown to remove high wavenumber errors while maintaining agreement with Direct Numerical Simulations (DNS) of forced and decaying isotropic turbulence. Linear interaction analysis is implemented to model the interaction of a shock with isotropic turbulence from LES. Comparisons to analytical models suggest that the regularization significantly improves the ability of the LES to predict amplifications in subgrid terms over the modeled shockwave. LES and DNS of decaying, modeled post shock turbulence are also considered, and inclusion of the regularization in shock-turbulence LES is shown to improve agreement with lower Reynolds number DNS.

Large eddy simulation for aerodynamics: status and perspectives.

PubMed

Sagaut, Pierre; Deck, Sébastien

2009-07-28

The present paper provides an up-to-date survey of the use of large eddy simulation (LES) and sequels for engineering applications related to aerodynamics. Most recent landmark achievements are presented. Two categories of problem may be distinguished whether the location of separation is triggered by the geometry or not. In the first case, LES can be considered as a mature technique and recent hybrid Reynolds-averaged Navier-Stokes (RANS)-LES methods do not allow for a significant increase in terms of geometrical complexity and/or Reynolds number with respect to classical LES. When attached boundary layers have a significant impact on the global flow dynamics, the use of hybrid RANS-LES remains the principal strategy to reduce computational cost compared to LES. Another striking observation is that the level of validation is most of the time restricted to time-averaged global quantities, a detailed analysis of the flow unsteadiness being missing. Therefore, a clear need for detailed validation in the near future is identified. To this end, new issues, such as uncertainty and error quantification and modelling, will be of major importance. First results dealing with uncertainty modelling in unsteady turbulent flow simulation are presented.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gustafson, William I.; Vogelmann, Andrew M.; Cheng, Xiaoping

The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility began a pilot project in May 2015 to design a routine, high-resolution modeling capability to complement ARM’s extensive suite of measurements. This modeling capability has been named the Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) project. The initial focus of LASSO is on shallow convection at the ARM Southern Great Plains (SGP) Climate Research Facility. The availability of LES simulations with concurrent observations will serve many purposes. LES helps bridge the scale gap between DOE ARM observations and models, and the use of routine LES addsmore » value to observations. It provides a self-consistent representation of the atmosphere and a dynamical context for the observations. Further, it elucidates unobservable processes and properties. LASSO will generate a simulation library for researchers that enables statistical approaches beyond a single-case mentality. It will also provide tools necessary for modelers to reproduce the LES and conduct their own sensitivity experiments. Many different uses are envisioned for the combined LASSO LES and observational library. For an observationalist, LASSO can help inform instrument remote sensing retrievals, conduct Observation System Simulation Experiments (OSSEs), and test implications of radar scan strategies or flight paths. For a theoretician, LASSO will help calculate estimates of fluxes and co-variability of values, and test relationships without having to run the model yourself. For a modeler, LASSO will help one know ahead of time which days have good forcing, have co-registered observations at high-resolution scales, and have simulation inputs and corresponding outputs to test parameterizations. Further details on the overall LASSO project are available at https://www.arm.gov/capabilities/modeling/lasso.« less

Dynamic large eddy simulation: Stability via realizability

NASA Astrophysics Data System (ADS)

Mokhtarpoor, Reza; Heinz, Stefan

2017-10-01

The concept of dynamic large eddy simulation (LES) is highly attractive: such methods can dynamically adjust to changing flow conditions, which is known to be highly beneficial. For example, this avoids the use of empirical, case dependent approximations (like damping functions). Ideally, dynamic LES should be local in physical space (without involving artificial clipping parameters), and it should be stable for a wide range of simulation time steps, Reynolds numbers, and numerical schemes. These properties are not trivial, but dynamic LES suffers from such problems over decades. We address these questions by performing dynamic LES of periodic hill flow including separation at a high Reynolds number Re = 37 000. For the case considered, the main result of our studies is that it is possible to design LES that has the desired properties. It requires physical consistency: a PDF-realizable and stress-realizable LES model, which requires the inclusion of the turbulent kinetic energy in the LES calculation. LES models that do not honor such physical consistency can become unstable. We do not find support for the previous assumption that long-term correlations of negative dynamic model parameters are responsible for instability. Instead, we concluded that instability is caused by the stable spatial organization of significant unphysical states, which are represented by wall-type gradient streaks of the standard deviation of the dynamic model parameter. The applicability of our realizability stabilization to other dynamic models (including the dynamic Smagorinsky model) is discussed.

Large-eddy simulation using the finite element method

DOE Office of Scientific and Technical Information (OSTI.GOV)

McCallen, R.C.; Gresho, P.M.; Leone, J.M. Jr.

1993-10-01

In a large-eddy simulation (LES) of turbulent flows, the large-scale motion is calculated explicitly (i.e., approximated with semi-empirical relations). Typically, finite difference or spectral numerical schemes are used to generate an LES; the use of finite element methods (FEM) has been far less prominent. In this study, we demonstrate that FEM in combination with LES provides a viable tool for the study of turbulent, separating channel flows, specifically the flow over a two-dimensional backward-facing step. The combination of these methodologies brings together the advantages of each: LES provides a high degree of accuracy with a minimum of empiricism for turbulencemore » modeling and FEM provides a robust way to simulate flow in very complex domains of practical interest. Such a combination should prove very valuable to the engineering community.« less

Flexible Inhibitor Fluid-Structure Interaction Simulation in RSRM.

NASA Astrophysics Data System (ADS)

Wasistho, Bono

2005-11-01

We employ our tightly coupled fluid/structure/combustion simulation code 'Rocstar-3' for solid propellant rocket motors to study 3D flows past rigid and flexible inhibitors in the Reusable Solid Rocket Motor (RSRM). We perform high resolution simulations of a section of the rocket near the center joint slot at 100 seconds after ignition, using inflow conditions based on less detailed 3D simulations of the full RSRM. Our simulations include both inviscid and turbulent flows (using LES dynamic subgrid-scale model), and explore the interaction between the inhibitor and the resulting fluid flow. The response of the solid components is computed by an implicit finite element solver. The internal mesh motion scheme in our block-structured fluid solver enables our code to handle significant changes in geometry. We compute turbulent statistics and determine the compound instabilities originated from the natural hydrodynamic instabilities and the inhibitor motion. The ultimate goal is to studdy the effect of inhibitor flexing on the turbulent field.

Towards Direct Numerical Simulation of mass and energy fluxes at the soil-atmospheric interface with advanced Lattice Boltzmann methods

NASA Astrophysics Data System (ADS)

Wang, Ying; Krafczyk, Manfred; Geier, Martin; Schönherr, Martin

2014-05-01

The quantification of soil evaporation and of soil water content dynamics near the soil surface are critical in the physics of land-surface processes on many scales and are dominated by multi-component and multi-phase mass and energy fluxes between the ground and the atmosphere. Although it is widely recognized that both liquid and gaseous water movement are fundamental factors in the quantification of soil heat flux and surface evaporation, their computation has only started to be taken into account using simplified macroscopic models. As the flow field over the soil can be safely considered as turbulent, it would be natural to study the detailed transient flow dynamics by means of Large Eddy Simulation (LES [1]) where the three-dimensional flow field is resolved down to the laminar sub-layer. Yet this requires very fine resolved meshes allowing a grid resolution of at least one order of magnitude below the typical grain diameter of the soil under consideration. In order to gain reliable turbulence statistics, up to several hundred eddy turnover times have to be simulated which adds up to several seconds of real time. Yet, the time scale of the receding saturated water front dynamics in the soil is on the order of hours. Thus we are faced with the task of solving a transient turbulent flow problem including the advection-diffusion of water vapour over the soil-atmospheric interface represented by a realistic tomographic reconstruction of a real porous medium taken from laboratory probes. Our flow solver is based on the Lattice Boltzmann method (LBM) [2] which has been extended by a Cumulant approach similar to the one described in [3,4] to minimize the spurious coupling between the degrees of freedom in previous LBM approaches and can be used as an implicit LES turbulence model due to its low numerical dissipation and increased stability at high Reynolds numbers. The kernel has been integrated into the research code Virtualfluids [5] and delivers up to 30% of the peak performance of modern General Purpose Graphics Processing Units (GPGPU, [6]) allowing the simulation of several minutes real-time for an LES LBM model. In our contribution we will present detailed profiles of the velocity distribution for different surface roughnesses, describe our multi-scale approach for the advection diffusion and estimate water vapour fluxes from transient simulations of the coupled problem. REFERENCES [1] J. Fröhlich and D. von Terzi. Hybrid LES/RANS methods for the simulation of turbulent flows. Progress in Aerospace Sciences, 44(5):349 - 377, 2008. [2] S. Chen and G. D. Doolen, Annual Review, of Fluid Mechanics 30, 329, 1998, [3] S. Seeger and K. H. Hoffmann, The cumulant method for computational kinetic theory, Continuum Mech. Thermodyn., 12:403-421, 2000. [4] S. Seeger and K. H. Hoffmann, The cumulant method applied to a mixture of Maxwell gases, Continuum Mech. Thermodyn., 14:321-335, 2002. [5] S. Freudiger, J. Hegewald and M. Krafczyk. A parallelisation concept for a mult-physics Lattice Boltzmann prototype based on hierarchical grids. Progress in Computational Fluid Dynamics, 8(1):168-178, 2008. [6] M. Schönherr, K. Kucher, M. Geier, M. Stiebler, S. Freudiger and M. Krafczyk, Multi- thread implementations of the Lattice Boltzmann method on non-uniform grids for CPUs and GPUs. Computers & Mathematics with Applications, 61(12):3730-3743, 2011.

Estimating the impact of internal climate variability on ice sheet model simulations

NASA Astrophysics Data System (ADS)

Tsai, C. Y.; Forest, C. E.; Pollard, D.

2016-12-01

Rising sea level threatens human societies and coastal habitats and melting ice sheets are a major contributor to sea level rise (SLR). Thus, understanding uncertainty of both forcing and variability within the climate system is essential for assessing long-term risk of SLR given their impact on ice sheet evolution. The predictability of polar climate is limited by uncertainties from the given forcing, the climate model response to this forcing, and the internal variability from feedbacks within the fully coupled climate system. Among those sources of uncertainty, the impact of internal climate variability on ice sheet changes has not yet been robustly assessed. Here we investigate how internal variability affects ice sheet projections using climate fields from two Community Earth System Model (CESM) large-ensemble (LE) experiments to force a three-dimensional ice sheet model. Each ensemble member in an LE experiment undergoes the same external forcings but with unique initial conditions. We find that for both LEs, 2m air temperature variability over Greenland ice sheet (GrIS) can lead to significantly different ice sheet responses. Our results show that the internal variability from two fully coupled CESM LEs can cause about 25 35 mm differences of GrIS's contribution to SLR in 2100 compared to present day (about 20% of the total change), and 100m differences of SLR in 2300. Moreover, only using ensemble-mean climate fields as the forcing in ice sheet model can significantly underestimate the melt of GrIS. As the Arctic region becomes warmer, the role of internal variability is critical given the complex nonlinear interactions between surface temperature and ice sheet. Our results demonstrate that internal variability from coupled atmosphere-ocean general circulation model can affect ice sheet simulations and the resulting sea-level projections. This study highlights an urgent need to reassess associated uncertainties of projecting ice sheet loss over the next few centuries to obtain robust estimates of the contribution of ice sheet melt to SLR.

Large-eddy simulation of flow in a plane, asymmetric diffuser

NASA Technical Reports Server (NTRS)

Kaltenbach, Hans-Jakob

1993-01-01

Recent improvements in subgrid-scale modeling as well as increases in computer power make it feasible to investigate flows using large-eddy simulation (LES) which have been traditionally studied with techniques based on Reynolds averaging. However, LES has not yet been applied to many flows of immediate technical interest. Preliminary results from LES of a plane diffuser flow are described. The long term goal of this work is to investigate flow separation as well as separation control in ducts and ramp-like geometries.

Eulerian Time-Domain Filtering for Spatial LES

NASA Technical Reports Server (NTRS)

Pruett, C. David

1997-01-01

Eulerian time-domain filtering seems to be appropriate for LES (large eddy simulation) of flows whose large coherent structures convect approximately at a common characteristic velocity; e.g., mixing layers, jets, and wakes. For these flows, we develop an approach to LES based on an explicit second-order digital Butterworth filter, which is applied in,the time domain in an Eulerian context. The approach is validated through a priori and a posteriori analyses of the simulated flow of a heated, subsonic, axisymmetric jet.

Marine boundary layer cloud regimes and POC formation in an LES coupled to a bulk aerosol scheme

NASA Astrophysics Data System (ADS)

Berner, A. H.; Bretherton, C. S.; Wood, R.; Muhlbauer, A.

2013-07-01

A large-eddy simulation (LES) coupled to a new bulk aerosol scheme is used to study long-lived regimes of aerosol-boundary layer cloud-precipitation interaction and the development of pockets of open cells (POCs) in subtropical stratocumulus cloud layers. The aerosol scheme prognoses mass and number concentration of a single log-normal accumulation mode with surface and entrainment sources, evolving subject to processing of activated aerosol and scavenging of dry aerosol by cloud and rain. The LES with the aerosol scheme is applied to a range of steadily-forced simulations idealized from a well-observed POC case. The long-term system evolution is explored with extended two-dimensional simulations of up to 20 days, mostly with diurnally-averaged insolation. One three-dimensional two-day simulation confirms the initial development of the corresponding two-dimensional case. With weak mean subsidence, an initially aerosol-rich mixed layer deepens, the capping stratocumulus cloud slowly thickens and increasingly depletes aerosol via precipitation accretion, then the boundary layer transitions within a few hours into an open-cell regime with scattered precipitating cumuli, in which entrainment is much weaker. The inversion slowly collapses for several days until the cumulus clouds are too shallow to efficiently precipitate. Inversion cloud then reforms and radiatively drives renewed entrainment, allowing the boundary layer to deepen and become more aerosol-rich, until the stratocumulus layer thickens enough to undergo another cycle of open-cell formation. If mean subsidence is stronger, the stratocumulus never thickens enough to initiate drizzle and settles into a steady state. With lower initial aerosol concentrations, this system quickly transitions into open cells, collapses, and redevelops into a different steady state with a shallow, optically thin cloud layer. In these steady states, interstitial scavenging by cloud droplets is the main sink of aerosol number. The system is described in a reduced two-dimensional phase plane with inversion height and boundary-layer average aerosol concentrations as the state variables. Simulations with a full diurnal cycle show similar evolutions, except that open-cell formation is phase-locked into the early morning hours. The same steadily-forced modeling framework is applied to the development and evolution of a POC and the surrounding overcast boundary layer. An initial aerosol perturbation applied to a portion of the model domain leads that portion to transition into open-cell convection, forming a POC. Reduced entrainment in the POC induces a negative feedback between areal fraction covered by the POC and boundary layer depth changes. This stabilizes the system by controlling liquid water path and precipitation sinks of aerosol number in the overcast region, while also preventing boundary-layer collapse within the POC, allowing the POC and overcast to coexist indefinitely in a quasi-steady equilibrium.

Impact of two chemistry mechanisms fully coupled with mesoscale model on the atmospheric pollutants distribution

NASA Astrophysics Data System (ADS)

Arteta, J.; Cautenet, S.; Taghavi, M.; Audiffren, N.

Air quality models (AQM) consist of many modules (meteorology, emission, chemistry, deposition), and in some conditions such as: vicinity of clouds or aerosols plumes, complex local circulations (mountains, sea breezes), fully coupled models (online method) are necessary. In order to study the impact of lumped chemical mechanisms in AQM simulations, we examine the ability of both different chemical mechanisms: (i) simplified: Condensed Version of the MOdèle de Chimie Atmosphérique 2.2 (CV-MOCA2.2), and (ii) reference: Regional Atmospheric Chemistry Model (RACM), which are coupled online with the Regional Atmospheric Modeling Systems (RAMS) model, on the distribution of pollutants. During the ESCOMPTE experiment (Expérience sur Site pour COntraindre les Modèles de Pollution et de Transport d'Emissions) conducted over Southern France (including urban and industrial zones), Intensive observation periods (IOP) characterized by various meteorological and mixed chemical conditions are simulated. For both configurations of modeling, numerical results are compared with surface measurements (75 stations) for primary (NO x) and secondary (O 3) species. We point out the impact of the two different chemical mechanisms on the production of species involved in the oxidizing capacity such as ozone and radicals within urban and industrial areas. We highlight that both chemical mechanisms produce very similar results for the main pollutants (NO x and O 3) in three-dimensional (3D) distribution, despite large discrepancies in 0D modeling. For ozone concentration, we found sometimes small differences (5-10 ppb) between the mechanisms under study according to the cases (polluted or not). The relative difference between the two mechanisms over the whole domain is only -7% for ozone from CV-MOCA 2.2 versus RACM. When the order of magnitude is needed rather than an accurate estimate, a reduced mechanism is satisfactory. It has the advantage of running faster (four times less than CPU time on SGI 3800 with 30 processors). Simplified mechanisms are really important to study cases for which an online coupling is necessary between meso-scale and chemistry models (clouds or aerosols plumes impacts, highly variable meteorology).

WRF nested large-eddy simulations of deep convection during SEAC4RS

NASA Astrophysics Data System (ADS)

Heath, Nicholas Kyle

Deep convection is an important component of atmospheric circulations that affects many aspects of weather and climate. Therefore, improved understanding and realistic simulations of deep convection are critical to both operational and climate forecasts. Large-eddy simulations (LESs) often are used with observations to enhance understanding of convective processes. This study develops and evaluates a nested-LES method using the Weather Research and Forecasting (WRF) model. Our goal is to evaluate the extent to which the WRF nested-LES approach is useful for studying deep convection during a real-world case. The method was applied on 2 September 2013, a day of continental convection having a robust set of ground and airborne data available for evaluation. A three domain mesoscale WRF simulation is run first. Then, the finest mesoscale output (1.35 km grid length) is used to separately drive nested-LES domains with grid lengths of 450 and 150 m. Results reveal that the nested-LES approach reasonably simulates a broad spectrum of observations, from reflectivity distributions to vertical velocity profiles, during the study period. However, reducing the grid spacing does not necessarily improve results for our case, with the 450 m simulation outperforming the 150 m version. We find that simulated updrafts in the 150 m simulation are too narrow to overcome the negative effects of entrainment, thereby generating convection that is weaker than observed. Increasing the sub-grid mixing length in the 150 m simulation leads to deeper, more realistic convection, but comes at the expense of delaying the onset of the convection. Overall, results show that both the 450 m and 150 m simulations are influenced considerably by the choice of sub-grid mixing length used in the LES turbulence closure. Finally, the simulations and observations are used to study the processes forcing strong midlevel cloud-edge downdrafts that were observed on 2 September. Results suggest that these downdrafts are forced by evaporative cooling due to mixing near cloud edge and by vertical perturbation pressure gradient forces acting to restore mass continuity around neighboring updrafts. We conclude that the WRF nested-LES approach provides an effective method for studying deep convection for our real-world case. The method can be used to provide insight into physical processes that are important to understanding observations. The WRF nested-LES approach could be adapted for other case studies in which high-resolution observations are available for validation.

Large eddy simulation for predicting turbulent heat transfer in gas turbines.

PubMed

Tafti, Danesh K; He, Long; Nagendra, K

2014-08-13

Blade cooling technology will play a critical role in the next generation of propulsion and power generation gas turbines. Accurate prediction of blade metal temperature can avoid the use of excessive compressed bypass air and allow higher turbine inlet temperature, increasing fuel efficiency and decreasing emissions. Large eddy simulation (LES) has been established to predict heat transfer coefficients with good accuracy under various non-canonical flows, but is still limited to relatively simple geometries and low Reynolds numbers. It is envisioned that the projected increase in computational power combined with a drop in price-to-performance ratio will make system-level simulations using LES in complex blade geometries at engine conditions accessible to the design process in the coming one to two decades. In making this possible, two key challenges are addressed in this paper: working with complex intricate blade geometries and simulating high-Reynolds-number (Re) flows. It is proposed to use the immersed boundary method (IBM) combined with LES wall functions. A ribbed duct at Re=20 000 is simulated using the IBM, and a two-pass ribbed duct is simulated at Re=100 000 with and without rotation (rotation number Ro=0.2) using LES with wall functions. The results validate that the IBM is a viable alternative to body-conforming grids and that LES with wall functions reproduces experimental results at a much lower computational cost. © 2014 The Author(s) Published by the Royal Society. All rights reserved.

Characterization and Modeling of Atmospheric Flow Within and Above Plant Canopies

NASA Astrophysics Data System (ADS)

Souza Freire Grion, Livia

The turbulent flow within and above plant canopies is responsible for the exchange of momentum, heat, gases and particles between vegetation and the atmosphere. Turbulence is also responsible for the mixing of air inside the canopy, playing an important role in chemical and biophysical processes occurring in the plants' environment. In the last fifty years, research has significantly advanced the understanding of and ability to model the flow field within and above the canopy, but important issues remain unsolved. In this work, we focus on (i) the estimation of turbulent mixing timescales within the canopy from field data; and (ii) the development of new computationally efficient modeling approaches for the coupled canopy-atmosphere flow field. The turbulent mixing timescale represents how quickly turbulence creates a well-mixed environment within the canopy. When the mixing timescale is much smaller than the timescale of other relevant processes (e.g. chemical reactions, deposition), the system can be assumed to be well-mixed and detailed modeling of turbulence is not critical to predict the system evolution. Conversely, if the mixing timescale is comparable or larger than the other timescales, turbulence becomes a controlling factor for the concentration of the variables involved; hence, turbulence needs to be taken into account when studying and modeling such processes. In this work, we used a combination of ozone concentration and high-frequency velocity data measured within and above the canopy in the Amazon rainforest to characterize turbulent mixing. The eddy diffusivity parameter (used as a proxy for mixing efficiency) was applied in a simple theoretical model of one-dimensional diffusion, providing an estimate of turbulent mixing timescales as a function of height within the canopy and time-of-day. Results showed that, during the day, the Amazon rainforest is characterized by well-mixed conditions with mixing timescales smaller than thirty minutes in the upper-half of the canopy, and partially mixed conditions in the lower half of the canopy. During the night, most of the canopy (except for the upper 20%) is either partially or poorly mixed, resulting in mixing timescales of up to several hours. For the specific case of ozone, the mixing timescales observed during the day are much lower than the chemical and deposition timescales, whereas chemical processes and turbulence have comparable timescales during the night. In addition, the high day-to-day variability in mixing conditions and the fast increase in mixing during the morning transition period indicate that turbulence within the canopy needs to be properly investigated and modeled in many studies involving plant-atmosphere interactions. Motivated by the findings described above, this work proposes and tests a new approach for modeling canopy flows. Typically, vertical profiles of flow statistics are needed to represent canopy-atmosphere exchanges in chemical and biophysical processes happening within the canopy. Current single-column models provide only steady-state (equilibrium) profiles, and rely on closure assumptions that do not represent the dominant non-local turbulent fluxes present in canopy flows. We overcome these issues by adapting the one-dimensional turbulent (ODT) model to represent atmospheric flows from the ground up to the top of the atmospheric boundary layer (ABL). The ODT model numerically resolves the one-dimensional diffusion equation along a vertical line (representing a horizontally homogeneous ABL column), and the presence of three-dimensional turbulence is added through the effect of stochastic eddies. Simulations of ABL without canopy were performed for different atmospheric stabilities and a diurnal cycle, to test the capabilities of this modeling approach in representing unsteady flows with strong non-local transport. In addition, four different types of canopies were simulated, one of them including the transport of scalar with a point source located inside the canopy. The comparison of all simulations with theory and field data provided satisfactory results. The main advantages of using ODT compared to typical 1D canopy-flow models are the ability to represent the coupled canopy-ABL flow with one single modeling approach, the presence of non-local turbulent fluxes, the ability to simulate transient conditions, the straightforward representation of multiple scalar fields, and the presence of only one adjustable parameter (as opposed to the several adjustable constants and boundary conditions needed for other modeling approaches). The results obtained with ODT as a stand-alone model motivated its use as a surface parameterization for Large-Eddy Simulation (LES). In this two-way coupling between LES and ODT, the former is used to simulate the ABL in a case where a canopy is present but cannot be resolved by the LES (i.e., the LES first vertical grid point is above the canopy). ODT is used to represent the flow field between the ground and the first LES grid point, including the region within and just above the canopy. In this work, we tested the ODT-LES model for three different types of canopies and obtained promising results. Although more work is needed in order to improve first and second-order statistics within the canopy (i.e. in the ODT domain), the results obtained for the flow statistics in the LES domain and for the third order statistics in the ODT domain demonstrate that the ODT-LES model is capable of capturing some important features of the canopy-atmosphere interaction. This new surface superparameterization approach using ODT provides a new alternative for simulations that require complex interactions between the flow field and near-surface processes (e.g. sand and snow drift, waves over water surfaces) and can potentially be extended to other large-scale models, such as mesoscale and global circulation models.

A new approach for turbulent simulations in complex geometries

NASA Astrophysics Data System (ADS)

Israel, Daniel M.

Historically turbulence modeling has been sharply divided into Reynolds averaged Navier-Stokes (RANS), in which all the turbulent scales of motion are modeled, and large-eddy simulation (LES), in which only a portion of the turbulent spectrum is modeled. In recent years there have been numerous attempts to couple these two approaches either by patching RANS and LES calculations together (zonal methods) or by blending the two sets of equations. In order to create a proper bridging model, that is, a single set of equations which captures both RANS and LES like behavior, it is necessary to place both RANS and LES in a more general framework. The goal of the current work is threefold: to provide such a framework, to demonstrate how the Flow Simulation Methodology (FSM) fits into this framework, and to evaluate the strengths and weaknesses of the current version of the FSM. To do this, first a set of filtered Navier-Stokes (FNS) equations are introduced in terms of an arbitrary generalized filter. Additional exact equations are given for the second order moments and the generalized subfilter dissipation rate tensor. This is followed by a discussion of the role of implicit and explicit filters in turbulence modeling. The FSM is then described with particular attention to its role as a bridging model. In order to evaluate the method a specific implementation of the FSM approach is proposed. Simulations are presented using this model for the case of a separating flow over a "hump" with and without flow control. Careful attention is paid to error estimation, and, in particular, how using flow statistics and time series affects the error analysis. Both mean flow and Reynolds stress profiles are presented, as well as the phase averaged turbulent structures and wall pressure spectra. Using the phase averaged data it is possible to examine how the FSM partitions the energy between the coherent resolved scale motions, the random resolved scale fluctuations, and the subfilter quantities. The method proves to be qualitatively successful at reproducing large turbulent structures. However, like other hybrid methods, it has difficulty in the region where the model behavior transitions from RANS to LES. Consequently the phase averaged structures reproduce the experiments quite well, and the forcing does significantly reduce the length of the separated region. Nevertheless, the recirculation length is significantly too large for all the cases. Overall the current results demonstrate the promise of bridging models in general and the FSM in particular. However, current bridging techniques are still in their infancy. There is still important progress to be made and it is hoped that this work points out the more important avenues for exploration.

Criteria for Modeling in LES of Multicomponent Fuel Flow

NASA Technical Reports Server (NTRS)

Bellan, Josette; Selle, Laurent

2009-01-01

A report presents a study addressing the question of which large-eddy simulation (LES) equations are appropriate for modeling the flow of evaporating drops of a multicomponent liquid in a gas (e.g., a spray of kerosene or diesel fuel in air). The LES equations are obtained from the direct numerical simulation (DNS) equations in which the solution is computed at all flow length scales, by applying a spatial low-pass filter. Thus, in LES the small scales are removed and replaced by terms that cannot be computed from the LES solution and instead must be modeled to retain the effect of the small scales into the equations. The mathematical form of these models is a subject of contemporary research. For a single-component liquid, there is only one LES formulation, but this study revealed that for a multicomponent liquid, there are two non-equivalent LES formulations for the conservation equations describing the composition of the vapor. Criteria were proposed for selecting the multicomponent LES formulation that gives the best accuracy and increased computational efficiency. These criteria were applied in examination of filtered DNS databases to compute the terms in the LES equations. The DNS databases are from mixing layers of diesel and kerosene fuels. The comparisons resulted in the selection of one of the multicomponent LES formulations as the most promising with respect to all criteria.

A Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

NASA Technical Reports Server (NTRS)

Yang, H. Q.; West, Jeff

2014-01-01

A capability to couple NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This paper summarizes the efforts in applying the installed coupling software to demonstrate/investigate fluid-structure interaction (FSI) between pressure wave and flexible inhibitor inside reusable solid rocket motor (RSRM). First a unified governing equation for both fluid and structure is presented, then an Eulerian-Lagrangian framework is described to satisfy the interfacial continuity requirements. The features of fluid solver, Loci/CHEM and structural solver, CoBi, are discussed before the coupling methodology of the solvers is described. The simulation uses production level CFD LES turbulence model with a grid resolution of 80 million cells. The flexible inhibitor is modeled with full 3D shell elements. Verifications against analytical solutions of structural model under steady uniform pressure condition and under dynamic condition of modal analysis show excellent agreements in terms of displacement distribution and eigen modal frequencies. The preliminary coupled result shows that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor.

Smoothed particle hydrodynamics method from a large eddy simulation perspective

NASA Astrophysics Data System (ADS)

Di Mascio, A.; Antuono, M.; Colagrossi, A.; Marrone, S.

2017-03-01

The Smoothed Particle Hydrodynamics (SPH) method, often used for the modelling of the Navier-Stokes equations by a meshless Lagrangian approach, is revisited from the point of view of Large Eddy Simulation (LES). To this aim, the LES filtering procedure is recast in a Lagrangian framework by defining a filter that moves with the positions of the fluid particles at the filtered velocity. It is shown that the SPH smoothing procedure can be reinterpreted as a sort of LES Lagrangian filtering, and that, besides the terms coming from the LES convolution, additional contributions (never accounted for in the SPH literature) appear in the equations when formulated in a filtered fashion. Appropriate closure formulas are derived for the additional terms and a preliminary numerical test is provided to show the main features of the proposed LES-SPH model.

Large-eddy simulation of turbulent flow with a surface-mounted two-dimensional obstacle

NASA Technical Reports Server (NTRS)

Yang, Kyung-Soo; Ferziger, Joel H.

1993-01-01

In this paper, we perform a large eddy simulation (LES) of turbulent flow in a channel containing a two-dimensional obstacle on one wall using a dynamic subgrid-scale model (DSGSM) at Re = 3210, based on bulk velocity above the obstacle and obstacle height; the wall layers are fully resolved. The low Re enables us to perform a DNS (Case 1) against which to validate the LES results. The LES with the DSGSM is designated Case 2. In addition, an LES with the conventional fixed model constant (Case 3) is conducted to allow identification of improvements due to the DSGSM. We also include LES at Re = 82,000 (Case 4) using conventional Smagorinsky subgrid-scale model and a wall-layer model. The results will be compared with the experiment of Dimaczek et al.

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

NASA Technical Reports Server (NTRS)

Fulton, Jesse A.; Edwards, Jack R.; Hassan, Hassan A.; Rockwell, Robert; Goyne, Christopher; McDaniel, James; Smith, Chad; Cutler, Andrew; Johansen, Craig; Danehy, Paul M.;

A survey of modelling methods for high-fidelity wind farm simulations using large eddy simulation.

PubMed

Breton, S-P; Sumner, J; Sørensen, J N; Hansen, K S; Sarmast, S; Ivanell, S

2017-04-13

Large eddy simulations (LES) of wind farms have the capability to provide valuable and detailed information about the dynamics of wind turbine wakes. For this reason, their use within the wind energy research community is on the rise, spurring the development of new models and methods. This review surveys the most common schemes available to model the rotor, atmospheric conditions and terrain effects within current state-of-the-art LES codes, of which an overview is provided. A summary of the experimental research data available for validation of LES codes within the context of single and multiple wake situations is also supplied. Some typical results for wind turbine and wind farm flows are presented to illustrate best practices for carrying out high-fidelity LES of wind farms under various atmospheric and terrain conditions.This article is part of the themed issue 'Wind energy in complex terrains'. © 2017 The Author(s).

A survey of modelling methods for high-fidelity wind farm simulations using large eddy simulation

PubMed Central

Sumner, J.; Sørensen, J. N.; Hansen, K. S.; Sarmast, S.; Ivanell, S.

2017-01-01

Large eddy simulations (LES) of wind farms have the capability to provide valuable and detailed information about the dynamics of wind turbine wakes. For this reason, their use within the wind energy research community is on the rise, spurring the development of new models and methods. This review surveys the most common schemes available to model the rotor, atmospheric conditions and terrain effects within current state-of-the-art LES codes, of which an overview is provided. A summary of the experimental research data available for validation of LES codes within the context of single and multiple wake situations is also supplied. Some typical results for wind turbine and wind farm flows are presented to illustrate best practices for carrying out high-fidelity LES of wind farms under various atmospheric and terrain conditions. This article is part of the themed issue ‘Wind energy in complex terrains’. PMID:28265021

Large eddy simulation for atmospheric boundary layer flow over flat and complex terrains

NASA Astrophysics Data System (ADS)

Han, Yi; Stoellinger, Michael; Naughton, Jonathan

2016-09-01

In this work, we present Large Eddy Simulation (LES) results of atmospheric boundary layer (ABL) flow over complex terrain with neutral stratification using the OpenFOAM-based simulator for on/offshore wind farm applications (SOWFA). The complete work flow to investigate the LES for the ABL over real complex terrain is described including meteorological-tower data analysis, mesh generation and case set-up. New boundary conditions for the lateral and top boundaries are developed and validated to allow inflow and outflow as required in complex terrain simulations. The turbulent inflow data for the terrain simulation is generated using a precursor simulation of a flat and neutral ABL. Conditionally averaged met-tower data is used to specify the conditions for the flat precursor simulation and is also used for comparison with the simulation results of the terrain LES. A qualitative analysis of the simulation results reveals boundary layer separation and recirculation downstream of a prominent ridge that runs across the simulation domain. Comparisons of mean wind speed, standard deviation and direction between the computed results and the conditionally averaged tower data show a reasonable agreement.

Large-eddy simulation of wind turbine wake interactions on locally refined Cartesian grids

NASA Astrophysics Data System (ADS)

Angelidis, Dionysios; Sotiropoulos, Fotis

2014-11-01

Performing high-fidelity numerical simulations of turbulent flow in wind farms remains a challenging issue mainly because of the large computational resources required to accurately simulate the turbine wakes and turbine/turbine interactions. The discretization of the governing equations on structured grids for mesoscale calculations may not be the most efficient approach for resolving the large disparity of spatial scales. A 3D Cartesian grid refinement method enabling the efficient coupling of the Actuator Line Model (ALM) with locally refined unstructured Cartesian grids adapted to accurately resolve tip vortices and multi-turbine interactions, is presented. Second order schemes are employed for the discretization of the incompressible Navier-Stokes equations in a hybrid staggered/non-staggered formulation coupled with a fractional step method that ensures the satisfaction of local mass conservation to machine zero. The current approach enables multi-resolution LES of turbulent flow in multi-turbine wind farms. The numerical simulations are in good agreement with experimental measurements and are able to resolve the rich dynamics of turbine wakes on grids containing only a small fraction of the grid nodes that would be required in simulations without local mesh refinement. This material is based upon work supported by the Department of Energy under Award Number DE-EE0005482 and the National Science Foundation under Award number NSF PFI:BIC 1318201.

Evaluation of Subgrid-Scale Models for Large Eddy Simulation of Compressible Flows

NASA Technical Reports Server (NTRS)

Blaisdell, Gregory A.

1996-01-01

The objective of this project was to evaluate and develop subgrid-scale (SGS) turbulence models for large eddy simulations (LES) of compressible flows. During the first phase of the project results from LES using the dynamic SGS model were compared to those of direct numerical simulations (DNS) of compressible homogeneous turbulence. The second phase of the project involved implementing the dynamic SGS model in a NASA code for simulating supersonic flow over a flat-plate. The model has been successfully coded and a series of simulations has been completed. One of the major findings of the work is that numerical errors associated with the finite differencing scheme used in the code can overwhelm the SGS model and adversely affect the LES results. Attached to this overview are three submitted papers: 'Evaluation of the Dynamic Model for Simulations of Compressible Decaying Isotropic Turbulence'; 'The effect of the formulation of nonlinear terms on aliasing errors in spectral methods'; and 'Large-Eddy Simulation of a Spatially Evolving Compressible Boundary Layer Flow'.

Developpement D'un Modele Climatique Regional: Fizr Simulation des Conditions de Janvier de la Cote Ouest Nord Americaine

NASA Astrophysics Data System (ADS)

Goyette, Stephane

1995-11-01

Le sujet de cette these concerne la modelisation numerique du climat regional. L'objectif principal de l'exercice est de developper un modele climatique regional ayant les capacites de simuler des phenomenes de meso-echelle spatiale. Notre domaine d'etude se situe sur la Cote Ouest nord americaine. Ce dernier a retenu notre attention a cause de la complexite du relief et de son controle sur le climat. Les raisons qui motivent cette etude sont multiples: d'une part, nous ne pouvons pas augmenter, en pratique, la faible resolution spatiale des modeles de la circulation generale de l'atmosphere (MCG) sans augmenter a outrance les couts d'integration et, d'autre part, la gestion de l'environnement exige de plus en plus de donnees climatiques regionales determinees avec une meilleure resolution spatiale. Jusqu'alors, les MCG constituaient les modeles les plus estimes pour leurs aptitudes a simuler le climat ainsi que les changements climatiques mondiaux. Toutefois, les phenomenes climatiques de fine echelle echappent encore aux MCG a cause de leur faible resolution spatiale. De plus, les repercussions socio-economiques des modifications possibles des climats sont etroitement liees a des phenomenes imperceptibles par les MCG actuels. Afin de circonvenir certains problemes inherents a la resolution, une approche pratique vise a prendre un domaine spatial limite d'un MCG et a y imbriquer un autre modele numerique possedant, lui, un maillage de haute resolution spatiale. Ce processus d'imbrication implique alors une nouvelle simulation numerique. Cette "retro-simulation" est guidee dans le domaine restreint a partir de pieces d'informations fournies par le MCG et forcee par des mecanismes pris en charge uniquement par le modele imbrique. Ainsi, afin de raffiner la precision spatiale des previsions climatiques de grande echelle, nous developpons ici un modele numerique appele FIZR, permettant d'obtenir de l'information climatique regionale valide a la fine echelle spatiale. Cette nouvelle gamme de modeles-interpolateurs imbriques qualifies d'"intelligents" fait partie de la famille des modeles dits "pilotes". L'hypothese directrice de notre etude est fondee sur la supposition que le climat de fine echelle est souvent gouverne par des forcages provenant de la surface plutot que par des transports atmospheriques de grande echelle spatiale. La technique que nous proposons vise donc a guider FIZR par la Dynamique echantillonnee d'un MCG et de la forcer par la Physique du MCG ainsi que par un forcage orographique de meso-echelle, en chacun des noeuds de la grille fine de calculs. Afin de valider la robustesse et la justesse de notre modele climatique regional, nous avons choisi la region de la Cote Ouest du continent nord americain. Elle est notamment caracterisee par une distribution geographique des precipitations et des temperatures fortement influencee par le relief sous-jacent. Les resultats d'une simulation d'un mois de janvier avec FIZR demontrent que nous pouvons simuler des champs de precipitations et de temperatures au niveau de l'abri beaucoup plus pres des observations climatiques comparativement a ceux simules a partir d'un MCG. Ces performances sont manifestement attribuees au forcage orographique de meso-echelle de meme qu'aux caracteristiques de surface determinees a fine echelle. Un modele similaire a FIZR peut, en principe, etre implante sur l'importe quel MCG, donc, tout organisme de recherche implique en modelisation numerique mondiale de grande echelle pourra se doter d'un el outil de regionalisation.

Skill of ship-following large-eddy simulations in reproducing MAGIC observations across the northeast Pacific stratocumulus to cumulus transition region

DOE PAGES

McGibbon, J.; Bretherton, C. S.

2017-03-17

During the Marine ARM GPCI Investigation of Clouds (MAGIC) in October 2011 to September 2012, a container ship making periodic cruises between Los Angeles, CA, and Honolulu, HI, was instrumented with surface meteorological, aerosol and radiation instruments, a cloud radar and ceilometer, and radiosondes. Here large-eddy simulation (LES) is performed in a ship-following frame of reference for 13 four day transects from the MAGIC field campaign. The goal is to assess if LES can skillfully simulate the broad range of observed cloud characteristics and boundary layer structure across the subtropical stratocumulus to cumulus transition region sampled during different seasons andmore » meteorological conditions. Results from Leg 15A, which sampled a particularly well-defined stratocumulus to cumulus transition, demonstrate the approach. The LES reproduces the observed timing of decoupling and transition from stratocumulus to cumulus and matches the observed evolution of boundary layer structure, cloud fraction, liquid water path, and precipitation statistics remarkably well. Considering the simulations of all 13 cruises, the LES skillfully simulates the mean diurnal variation of key measured quantities, including liquid water path (LWP), cloud fraction, measures of decoupling, and cloud radar-derived precipitation. The daily mean quantities are well represented, and daily mean LWP and cloud fraction show the expected correlation with estimated inversion strength. There is a –0.6 K low bias in LES near-surface air temperature that results in a high bias of 5.6 W m –2 in sensible heat flux (SHF). Altogether, these results build confidence in the ability of LES to represent the northeast Pacific stratocumulus to trade cumulus transition region.« less

Large Eddy Simulation of High Reynolds Number Complex Flows

NASA Astrophysics Data System (ADS)

Verma, Aman

Marine configurations are subject to a variety of complex hydrodynamic phenomena affecting the overall performance of the vessel. The turbulent flow affects the hydrodynamic drag, propulsor performance and structural integrity, control-surface effectiveness, and acoustic signature of the marine vessel. Due to advances in massively parallel computers and numerical techniques, an unsteady numerical simulation methodology such as Large Eddy Simulation (LES) is well suited to study such complex turbulent flows whose Reynolds numbers (Re) are typically on the order of 10. 6. LES also promises increasedaccuracy over RANS based methods in predicting unsteady phenomena such as cavitation and noise production. This dissertation develops the capability to enable LES of high Re flows in complex geometries (e.g. a marine vessel) on unstructured grids and provide physical insight into the turbulent flow. LES is performed to investigate the geometry induced separated flow past a marine propeller attached to a hull, in an off-design condition called crashback. LES shows good quantitative agreement with experiments and provides a physical mechanism to explain the increase in side-force on the propeller blades below an advance ratio of J=-0.7. Fundamental developments in the dynamic subgrid-scale model for LES are pursued to improve the LES predictions, especially for complex flows on unstructured grids. A dynamic procedure is proposed to estimate a Lagrangian time scale based on a surrogate correlation without any adjustable parameter. The proposed model is applied to turbulent channel, cylinder and marine propeller flows and predicts improved results over other model variants due to a physically consistent Lagrangian time scale. A wall model is proposed for application to LES of high Reynolds number wall-bounded flows. The wall model is formulated as the minimization of a generalized constraint in the dynamic model for LES and applied to LES of turbulent channel flow at various Reynolds numbers up to Reτ=10000 and coarse grid resolutions to obtain significant improvement.

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

NASA Technical Reports Server (NTRS)

Givi, P.; Frankel, S. H.; Adumitroaie, V.; Sabini, G.; Madnia, C. K.

1993-01-01

The primary objective of this research is to extend current capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first two years of this research have been concentrated on a priori investigations of single-point Probability Density Function (PDF) methods for providing subgrid closures in reacting turbulent flows. In the efforts initiated in the third year, our primary focus has been on performing actual LES by means of PDF methods. The approach is based on assumed PDF methods and we have performed extensive analysis of turbulent reacting flows by means of LES. This includes simulations of both three-dimensional (3D) isotropic compressible flows and two-dimensional reacting planar mixing layers. In addition to these LES analyses, some work is in progress to assess the extent of validity of our assumed PDF methods. This assessment is done by making detailed companions with recent laboratory data in predicting the rate of reactant conversion in parallel reacting shear flows. This report provides a summary of our achievements for the first six months of the third year of this program.

Heart rate and pulmonary function while wearing the launch-entry crew escape suit (LES) during + Gx acceleration and simulated Shuttle launch

NASA Technical Reports Server (NTRS)

Krutz, Robert W., Jr.; Bagian, James P.; Burton, Russell R.; Meeker, Larry J.

1990-01-01

Space shuttle crewmembers have been equipped with a launch-entry crew escape system (LES) since the Challenger accident in 1986. Some crewmembers, wearing the new pressure suit, have reported breathing difficulties and increased effort to achieve the desired range of motion. This study was conducted to quantify the reported increased physical workloads and breathing difficulty associated with wearing the LES. Both veteran astronauts and centrifuge panel members were exposed to various + Gx profiles (including simulated shuttle launch) + Gx on the USAF School of Aerospace Medicine (USAFSAM) human-use centrifuge. Maximum heart rate data showed no increased workload associated with arm and head movement in the LES when compared to the flight suit/helmet ensemble (LEH). However, the LES did impose a significant increase in breathing difficulty beginning at +2.5 Gx which was demonstrated by a decrease in forced vital capacity and subjected questionnaries.

Large eddy simulation of flows in industrial compressors: a path from 2015 to 2035

PubMed Central

Gourdain, N.; Sicot, F.; Duchaine, F.; Gicquel, L.

2014-01-01

A better understanding of turbulent unsteady flows is a necessary step towards a breakthrough in the design of modern compressors. Owing to high Reynolds numbers and very complex geometry, the flow that develops in such industrial machines is extremely hard to predict. At this time, the most popular method to simulate these flows is still based on a Reynolds-averaged Navier–Stokes approach. However, there is some evidence that this formalism is not accurate for these components, especially when a description of time-dependent turbulent flows is desired. With the increase in computing power, large eddy simulation (LES) emerges as a promising technique to improve both knowledge of complex physics and reliability of flow solver predictions. The objective of the paper is thus to give an overview of the current status of LES for industrial compressor flows as well as to propose future research axes regarding the use of LES for compressor design. While the use of wall-resolved LES for industrial multistage compressors at realistic Reynolds number should not be ready before 2035, some possibilities exist to reduce the cost of LES, such as wall modelling and the adaptation of the phase-lag condition. This paper also points out the necessity to combine LES to techniques able to tackle complex geometries. Indeed LES alone, i.e. without prior knowledge of such flows for grid construction or the prohibitive yet ideal use of fully homogeneous meshes to predict compressor flows, is quite limited today. PMID:25024422

Numerics and subgrid-scale modeling in large eddy simulations of stratocumulus clouds.

PubMed

Pressel, Kyle G; Mishra, Siddhartha; Schneider, Tapio; Kaul, Colleen M; Tan, Zhihong

2017-06-01

Stratocumulus clouds are the most common type of boundary layer cloud; their radiative effects strongly modulate climate. Large eddy simulations (LES) of stratocumulus clouds often struggle to maintain fidelity to observations because of the sharp gradients occurring at the entrainment interfacial layer at the cloud top. The challenge posed to LES by stratocumulus clouds is evident in the wide range of solutions found in the LES intercomparison based on the DYCOMS-II field campaign, where simulated liquid water paths for identical initial and boundary conditions varied by a factor of nearly 12. Here we revisit the DYCOMS-II RF01 case and show that the wide range of previous LES results can be realized in a single LES code by varying only the numerical treatment of the equations of motion and the nature of subgrid-scale (SGS) closures. The simulations that maintain the greatest fidelity to DYCOMS-II observations are identified. The results show that using weighted essentially non-oscillatory (WENO) numerics for all resolved advective terms and no explicit SGS closure consistently produces the highest-fidelity simulations. This suggests that the numerical dissipation inherent in WENO schemes functions as a high-quality, implicit SGS closure for this stratocumulus case. Conversely, using oscillatory centered difference numerical schemes for momentum advection, WENO numerics for scalars, and explicitly modeled SGS fluxes consistently produces the lowest-fidelity simulations. We attribute this to the production of anomalously large SGS fluxes near the cloud tops through the interaction of numerical error in the momentum field with the scalar SGS model.

Numerics and subgrid‐scale modeling in large eddy simulations of stratocumulus clouds

PubMed Central

Mishra, Siddhartha; Schneider, Tapio; Kaul, Colleen M.; Tan, Zhihong

2017-01-01

Abstract Stratocumulus clouds are the most common type of boundary layer cloud; their radiative effects strongly modulate climate. Large eddy simulations (LES) of stratocumulus clouds often struggle to maintain fidelity to observations because of the sharp gradients occurring at the entrainment interfacial layer at the cloud top. The challenge posed to LES by stratocumulus clouds is evident in the wide range of solutions found in the LES intercomparison based on the DYCOMS‐II field campaign, where simulated liquid water paths for identical initial and boundary conditions varied by a factor of nearly 12. Here we revisit the DYCOMS‐II RF01 case and show that the wide range of previous LES results can be realized in a single LES code by varying only the numerical treatment of the equations of motion and the nature of subgrid‐scale (SGS) closures. The simulations that maintain the greatest fidelity to DYCOMS‐II observations are identified. The results show that using weighted essentially non‐oscillatory (WENO) numerics for all resolved advective terms and no explicit SGS closure consistently produces the highest‐fidelity simulations. This suggests that the numerical dissipation inherent in WENO schemes functions as a high‐quality, implicit SGS closure for this stratocumulus case. Conversely, using oscillatory centered difference numerical schemes for momentum advection, WENO numerics for scalars, and explicitly modeled SGS fluxes consistently produces the lowest‐fidelity simulations. We attribute this to the production of anomalously large SGS fluxes near the cloud tops through the interaction of numerical error in the momentum field with the scalar SGS model. PMID:28943997

Simulation of turbulent separated flows using a novel, evolution-based, eddy-viscosity formulation

NASA Astrophysics Data System (ADS)

Castellucci, Paul

Currently, there exists a lack of confidence in the computational simulation of turbulent separated flows at large Reynolds numbers. The most accurate methods available are too computationally costly to use in engineering applications. Thus, inexpensive models, developed using the Reynolds-averaged Navier-Stokes (RANS) equations, are often extended beyond their applicability. Although these methods will often reproduce integrated quantities within engineering tolerances, such metrics are often insensitive to details within a separated wake, and therefore, poor indicators of simulation fidelity. Using concepts borrowed from large-eddy simulation (LES), a two-equation RANS model is modified to simulate the turbulent wake behind a circular cylinder. This modification involves the computation of one additional scalar field, adding very little to the overall computational cost. When properly inserted into the baseline RANS model, this modification mimics LES in the separated wake, yet reverts to the unmodified form at the cylinder surface. In this manner, superior predictive capability may be achieved without the additional cost of fine spatial resolution associated with LES near solid boundaries. Simulations using modified and baseline RANS models are benchmarked against both LES and experimental data for a circular cylinder wake at Reynolds number 3900. In addition, the computational tool used in this investigation is subject to verification via the Method of Manufactured Solutions. Post-processing of the resultant flow fields includes both mean value and triple-decomposition analysis. These results reveal substantial improvements using the modified system and appear to drive the baseline wake solution toward that of LES, as intended.

Structured Overlapping Grid Simulations of Contra-rotating Open Rotor Noise

NASA Technical Reports Server (NTRS)

Housman, Jeffrey A.; Kiris, Cetin C.

2015-01-01

Computational simulations using structured overlapping grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting tonal noise generated by a contra-rotating open rotor (CROR) propulsion system. A coupled Computational Fluid Dynamics (CFD) and Computational AeroAcoustics (CAA) numerical approach is applied. Three-dimensional time-accurate hybrid Reynolds Averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) CFD simulations are performed in the inertial frame, including dynamic moving grids, using a higher-order accurate finite difference discretization on structured overlapping grids. A higher-order accurate free-stream preserving metric discretization with discrete enforcement of the Geometric Conservation Law (GCL) on moving curvilinear grids is used to create an accurate, efficient, and stable numerical scheme. The aeroacoustic analysis is based on a permeable surface Ffowcs Williams-Hawkings (FW-H) approach, evaluated in the frequency domain. A time-step sensitivity study was performed using only the forward row of blades to determine an adequate time-step. The numerical approach is validated against existing wind tunnel measurements.

Large-Eddy Simulation (LES) of a Compressible Mixing Layer and the Significance of Inflow Turbulence

NASA Technical Reports Server (NTRS)

Mankbadi, Mina Reda; Georgiadis, Nicholas J.; Debonis, James R.

2017-01-01

In the context of Large Eddy Simulations (LES), the effects of inflow turbulence are investigated through the Synthetic Eddy Method (SEM). The growth rate of a turbulent compressible mixing layer corresponding to operating conditions of GeobelDutton Case 2 is investigated herein. The effects of spanwise width on the growth rate of the mixing layer is investigated such that spanwise width independence is reached. The error in neglecting inflow turbulence effects is quantified by comparing two methodologies: (1) Hybrid-RANS-LES methodology and (2) SEM-LES methodology. Best practices learned from Case 2 are developed herein and then applied to a higher convective mach number corresponding to Case 4 experiments of GeobelDutton.

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

NASA Technical Reports Server (NTRS)

Bisset, D. K.; Rogers, M. M.; Kega, Dennis (Technical Monitor)

1999-01-01

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Perceptions des jeunes victimes de violence sexuelle au sein de leurs relations amoureuses sur leur pire expérience

PubMed Central

Van Camp, Tinneke; Hébert, Martine; Fernet, Mylène; Blais, Martin; Lavoie, Francine

2016-01-01

Cette étude explore les pires expériences vécues dans les relations amoureuses de jeunes qui ont rapporté avoir vécu de la violence sexuelle dans une relation de couple récente. Quelles sont les situations perçues comme étant les plus difficiles par les jeunes et est-ce que celles-ci se limitent à des incidents violents? Le questionnaire sur les parcours amoureux des jeunes (PAJ) a été complété par des jeunes Québécois âgés de 14 à 18 ans. Au total, plus de 600 participants ont rapporté au moins un épisode de violence sexuelle (souvent en combinaison avec d’autres formes de violence). Nous présentons les résultats de l’analyse qualitative inductive fondée sur une question ouverte concernant la pire expérience vécue. Les observations suggèrent que, en plus des expériences de violence, les difficultés relationnelles, les ruptures amoureuses et les sentiments amoureux non réciproques sont des situations particulièrement difficiles selon les propos des jeunes. Ces différents enjeux vécus par les jeunes devraient être pris en considération dans l’offre de services d’intervention à leur intention. PMID:28191266

Detached-Eddy Simulations of Attached and Detached Boundary Layers

NASA Astrophysics Data System (ADS)

Caruelle, B.; Ducros, F.

2003-12-01

This article presents Detached-Eddy Simulations (DESs) of attached and detached turbulent boundary layers. This hybrid Reynolds Averaged Navier-Stokes (RANS) / Large Eddy Simulation (LES) model goes continuously from RANS to LES according to the mesh definition. We propose a parametric study of the model over two "academic" configurations, in order to get information on the influence of the mesh to correctly treat complex flow with attached and detached boundary layers.

Studying marine stratus with large eddy simulation

NASA Technical Reports Server (NTRS)

Moeng, Chin-Hoh

1990-01-01

Data sets from field experiments over the stratocumulus regime may include complications from larger scale variations, decoupled cloud layers, diurnal cycle, or entrainment instability, etc. On top of the already complicated turbulence-radiation-condensation processes within the cloud-topped boundary layer (CTBL), these complexities may sometimes make interpretation of the data sets difficult. To study these processes, a better understanding is needed of the basic processes involved in the prototype CTBL. For example, is cloud top radiative cooling the primary source of the turbulent kinetic energy (TKE) within the CTBL. Historically, laboratory measurements have played an important role in addressing the turbulence problems. The CTBL is a turbulent field which is probably impossible to generate in laboratories. Large eddy simulation (LES) is an alternative way of 'measuring' the turbulent structure under controlled environments, which allows the systematic examination of the basic physical processes involved. However, there are problems with the LES approach for the CTBL. The LES data need to be consistent with the observed data. The LES approach is discussed, and results are given which provide some insights into the simulated turbulent flow field. Problems with this approach for the CTBL and information from the FIRE experiment needed to justify the LES results are discussed.

Design optimization using adjoint of Long-time LES for the trailing edge of a transonic turbine vane

NASA Astrophysics Data System (ADS)

Talnikar, Chaitanya; Wang, Qiqi

2017-11-01

Adjoint-based design optimization methods have been applied to low-fidelity simulation methods like Reynolds Averaged Navier-Stokes (RANS) and are useful for designing fluid machinery components. But to reliably capture the complex flow phenomena involved in turbomachinery, high fidelity simulations like large eddy simulation (LES) are required. Unfortunately due to the chaotic dynamics of turbulence, the unsteady adjoint method for LES diverges and produces incorrect gradients. Using a viscosity stabilized unsteady adjoint method developed for LES, the gradient can be obtained with reasonable accuracy. In this paper, design of the trailing edge of a gas turbine inlet guide vane is performed with the objective to reduce stagnation pressure loss and heat transfer over the surface of the vane. Slight changes in the shape of trailing edge can significantly impact these quantities by altering the boundary layer development process and separation points. The trailing edge is parameterized using a linear combination of 5 convex designs. Bayesian optimization is used as a global optimizer with the objective function evaluated from the LES and gradients obtained using the viscosity adjoint method. Results from the optimization, performed on the supercomputer Mira, are presented.

Multi-scale properties of large eddy simulations: correlations between resolved-scale velocity-field increments and subgrid-scale quantities

NASA Astrophysics Data System (ADS)

Linkmann, Moritz; Buzzicotti, Michele; Biferale, Luca

2018-06-01

We provide analytical and numerical results concerning multi-scale correlations between the resolved velocity field and the subgrid-scale (SGS) stress-tensor in large eddy simulations (LES). Following previous studies for Navier-Stokes equations, we derive the exact hierarchy of LES equations governing the spatio-temporal evolution of velocity structure functions of any order. The aim is to assess the influence of the subgrid model on the inertial range intermittency. We provide a series of predictions, within the multifractal theory, for the scaling of correlation involving the SGS stress and we compare them against numerical results from high-resolution Smagorinsky LES and from a-priori filtered data generated from direct numerical simulations (DNS). We find that LES data generally agree very well with filtered DNS results and with the multifractal prediction for all leading terms in the balance equations. Discrepancies are measured for some of the sub-leading terms involving cross-correlation between resolved velocity increments and the SGS tensor or the SGS energy transfer, suggesting that there must be room to improve the SGS modelisation to further extend the inertial range properties for any fixed LES resolution.

Adaptability in Coalition Teamwork (Faculte d’adaptation au travail d’ equipe en coalition)

DTIC Science & Technology

2008-04-01

et des outils sont nécessaires au développement rapide d’équipes multiculturelles efficaces pour assurer le succès des missions, celles-ci étant...Les principaux résultats des 30 communications théoriques et de recherche ont été les suivants : • Les outils de formation (jeux, simulations...parmi les militaires ; • Le retour d’information sur le moral et les performances des équipes en opérations est un instrument qui est particulièrement

A priori and a posteriori analysis of the flow around a rectangular cylinder

NASA Astrophysics Data System (ADS)

Cimarelli, A.; Leonforte, A.; Franciolini, M.; De Angelis, E.; Angeli, D.; Crivellini, A.

2017-11-01

The definition of a correct mesh resolution and modelling approach for the Large Eddy Simulation (LES) of the flow around a rectangular cylinder is recognized to be a rather elusive problem as shown by the large scatter of LES results present in the literature. In the present work, we aim at assessing this issue by performing an a priori analysis of Direct Numerical Simulation (DNS) data of the flow. This approach allows us to measure the ability of the LES field on reproducing the main flow features as a function of the resolution employed. Based on these results, we define a mesh resolution which maximize the opposite needs of reducing the computational costs and of adequately resolving the flow dynamics. The effectiveness of the resolution method proposed is then verified by means of an a posteriori analysis of actual LES data obtained by means of the implicit LES approach given by the numerical properties of the Discontinuous Galerkin spatial discretization technique. The present work represents a first step towards a best practice for LES of separating and reattaching flows.

A unified RANS–LES model: Computational development, accuracy and cost

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gopalan, Harish, E-mail: hgopalan@uwyo.edu; Heinz, Stefan, E-mail: heinz@uwyo.edu; Stöllinger, Michael K., E-mail: MStoell@uwyo.edu

2013-09-15

Large eddy simulation (LES) is computationally extremely expensive for the investigation of wall-bounded turbulent flows at high Reynolds numbers. A way to reduce the computational cost of LES by orders of magnitude is to combine LES equations with Reynolds-averaged Navier–Stokes (RANS) equations used in the near-wall region. A large variety of such hybrid RANS–LES methods are currently in use such that there is the question of which hybrid RANS-LES method represents the optimal approach. The properties of an optimal hybrid RANS–LES model are formulated here by taking reference to fundamental properties of fluid flow equations. It is shown that unifiedmore » RANS–LES models derived from an underlying stochastic turbulence model have the properties of optimal hybrid RANS–LES models. The rest of the paper is organized in two parts. First, a priori and a posteriori analyses of channel flow data are used to find the optimal computational formulation of the theoretically derived unified RANS–LES model and to show that this computational model, which is referred to as linear unified model (LUM), does also have all the properties of an optimal hybrid RANS–LES model. Second, a posteriori analyses of channel flow data are used to study the accuracy and cost features of the LUM. The following conclusions are obtained. (i) Compared to RANS, which require evidence for their predictions, the LUM has the significant advantage that the quality of predictions is relatively independent of the RANS model applied. (ii) Compared to LES, the significant advantage of the LUM is a cost reduction of high-Reynolds number simulations by a factor of 0.07Re{sup 0.46}. For coarse grids, the LUM has a significant accuracy advantage over corresponding LES. (iii) Compared to other usually applied hybrid RANS–LES models, it is shown that the LUM provides significantly improved predictions.« less

DNS and LES/FMDF of turbulent jet ignition and combustion

NASA Astrophysics Data System (ADS)

Validi, Abdoulahad; Jaberi, Farhad

2014-11-01

The ignition and combustion of lean fuel-air mixtures by a turbulent jet flow of hot combustion products injected into various geometries are studied by high fidelity numerical models. Turbulent jet ignition (TJI) is an efficient method for starting and controlling the combustion in complex propulsion systems and engines. The TJI and combustion of hydrogen and propane in various flow configurations are simulated with the direct numerical simulation (DNS) and the hybrid large eddy simulation/filtered mass density function (LES/FMDF) models. In the LES/FMDF model, the filtered form of the compressible Navier-Stokes equations are solved with a high-order finite difference scheme for the turbulent velocity and the FMDF transport equation is solved with a Lagrangian stochastic method to obtain the scalar field. The DNS and LES/FMDF data are used to study the physics of TJI and combustion for different turbulent jet igniter and gas mixture conditions. The results show the very complex and different behavior of the turbulence and the flame structure at different jet equivalence ratios.

LES of Temporally Evolving Mixing Layers by an Eighth-Order Filter Scheme

NASA Technical Reports Server (NTRS)

Hadjadj, A; Yee, H. C.; Sjogreen, B.

2011-01-01

An eighth-order filter method for a wide range of compressible flow speeds (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) are employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) and Reynolds numbers. The high order filter method is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The value of Mc considered is for the TML range from the quasi-incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self similarity property, compressibility effects and the presence of large-scale structure with shocklets for high Mc) are considered for the LES study. The LES results using the same scheme parameters for all studied cases agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002).

LES, DNS and RANS for the analysis of high-speed turbulent reacting flows

NASA Technical Reports Server (NTRS)

Adumitroaie, V.; Colucci, P. J.; Taulbee, D. B.; Givi, P.

1995-01-01

The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES), direct numerical simulation (DNS), and Reynolds averaged Navier Stokes (RANS) methods for the computational analysis of high-speed reacting turbulent flows. In the second phase of this work, covering the period 1 Aug. 1994 - 31 Jul. 1995, we have focused our efforts on two programs: (1) developments of explicit algebraic moment closures for statistical descriptions of compressible reacting flows and (2) development of Monte Carlo numerical methods for LES of chemically reacting flows.

Uncertainty quantification in LES of channel flow

DOE PAGES

Safta, Cosmin; Blaylock, Myra; Templeton, Jeremy; ...

2016-07-12

Here, in this paper, we present a Bayesian framework for estimating joint densities for large eddy simulation (LES) sub-grid scale model parameters based on canonical forced isotropic turbulence direct numerical simulation (DNS) data. The framework accounts for noise in the independent variables, and we present alternative formulations for accounting for discrepancies between model and data. To generate probability densities for flow characteristics, posterior densities for sub-grid scale model parameters are propagated forward through LES of channel flow and compared with DNS data. Synthesis of the calibration and prediction results demonstrates that model parameters have an explicit filter width dependence andmore » are highly correlated. Discrepancies between DNS and calibrated LES results point to additional model form inadequacies that need to be accounted for.« less

Étude comparative des techniques d'analyse par fluorescence X à dispersion d'énergie (ED-XRF) et à dispersion de longueur d'onde (WD-XRF), et par spectrométrie d'émission atomique à source plasma couplée par induction (ICP-AES)

NASA Astrophysics Data System (ADS)

Rahmani, A.; Benyaïch, F.; Bounakhla, M.; Bilal, E.; Moutte, J.; Gruffat, J. J.; Zahry, F.

2004-11-01

Dans ce travail, nous présentons une étude comparative des techniques d'analyse par fluorescence X à dispersion d'énergie (ED-XRF) et à dispersion de longueur d'onde (WD-XRF), et par spectrométrie d'émission atomique à source plasma couplée par induction (ICP-AES). Les résultats de la calibration des spectromètres à dispersion d'énergie, à excitation par sources radioactives (55Fe, 109Cd et 241Am) et à excitation secondaire (cible secondaire Mo et Cu) du Centre National pour l'Energie, les Sciences et les Techniques Nucléaires (CNESTEN, Rabat, Maroc) sur des échantillons étalons de références de l'Agence International de l'Energie Atomique (AIEA) et du Community Bureau of Référence (BCR) ont été comparés aux résultats d'analyse des mêmes échantillons étalons par la spectrométrie X à dispersion de longueur d'onde (WD-XRF) et par spectrométrie d'émission atomique à source plasma couplé par induction (ICP-AES) au département GENERIC du centre SPIN à l'Ecole des Mines de Saint-Etienne (France). Les trois techniques d'analyse utilisées donnent des résultats comparables pour le dosage des éléments majeurs, alors que pour les traces on note des déviations importantes à cause des effets de matrice qui sont difficiles à corriger dans le cas de la fluorescence X.

Boundary Layer Simulation and Control in Wind Tunnels

DTIC Science & Technology

1988-04-01

de l’art dans le domaine de la simulation de la couche limite, ou le nombre de Reynolds n’est pas ou ne peut pas etre simule, examine les effets ...pour la definition de certains essais en soufflerie oil les effets visqueux sont d’une importance particuliere. CONTENTS Page PREFACE iii 1...transition associated with cylindrical bodies at high incidence in subsonic flow. Other relevant references are given therein. Figures 11 a-b, from Ref

GPU Accelerated DG-FDF Large Eddy Simulator

NASA Astrophysics Data System (ADS)

Inkarbekov, Medet; Aitzhan, Aidyn; Sammak, Shervin; Givi, Peyman; Kaltayev, Aidarkhan

2017-11-01

A GPU accelerated simulator is developed and implemented for large eddy simulation (LES) of turbulent flows. The filtered density function (FDF) is utilized for modeling of the subgrid scale quantities. The filtered transport equations are solved via a discontinuous Galerkin (DG) and the FDF is simulated via particle based Lagrangian Monte-Carlo (MC) method. It is demonstrated that the GPUs simulations are of the order of 100 times faster than the CPU-based calculations. This brings LES of turbulent flows to a new level, facilitating efficient simulation of more complex problems. The work at Al-Faraby Kazakh National University is sponsored by MoES of RK under Grant 3298/GF-4.

Consistent Large-Eddy Simulation of a Temporal Mixing Layer Laden with Evaporating Drops. Part 2; A Posteriori Modelling

NASA Technical Reports Server (NTRS)

Leboissertier, Anthony; Okong'O, Nora; Bellan, Josette

2005-01-01

Large-eddy simulation (LES) is conducted of a three-dimensional temporal mixing layer whose lower stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase equations are written in an Eulerian frame for two perfect gas species (carrier gas and vapour emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame. The effect of drop evaporation on the gas phase is considered through mass, species, momentum and energy source terms. The drop evolution is modelled using physical drops, or using computational drops to represent the physical drops. Simulations are performed using various LES models previously assessed on a database obtained from direct numerical simulations (DNS). These LES models are for: (i) the subgrid-scale (SGS) fluxes and (ii) the filtered source terms (FSTs) based on computational drops. The LES, which are compared to filtered-and-coarsened (FC) DNS results at the coarser LES grid, are conducted with 64 times fewer grid points than the DNS, and up to 64 times fewer computational than physical drops. It is found that both constant-coefficient and dynamic Smagorinsky SGS-flux models, though numerically stable, are overly dissipative and damp generated small-resolved-scale (SRS) turbulent structures. Although the global growth and mixing predictions of LES using Smagorinsky models are in good agreement with the FC-DNS, the spatial distributions of the drops differ significantly. In contrast, the constant-coefficient scale-similarity model and the dynamic gradient model perform well in predicting most flow features, with the latter model having the advantage of not requiring a priori calibration of the model coefficient. The ability of the dynamic models to determine the model coefficient during LES is found to be essential since the constant-coefficient gradient model, although more accurate than the Smagorinsky model, is not consistently numerically stable despite using DNS-calibrated coefficients. With accurate SGS-flux models, namely scale-similarity and dynamic gradient, the FST model allows up to a 32-fold reduction in computational drops compared to the number of physical drops, without degradation of accuracy; a 64-fold reduction leads to a slight decrease in accuracy.

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

NASA Technical Reports Server (NTRS)

Anderson, William E.; Lucht, Robert P.; Mongia, Hukam

2015-01-01

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

Numerical simulation of flow around the NREL S826 airfoil at moderate Reynolds number using delayed detached Eddy simulation (DDES)

NASA Astrophysics Data System (ADS)

Prytz, Erik R.; Huuse, Øyvind; Müller, Bernhard; Bartl, Jan; Sætran, Lars Roar

2017-07-01

Turbulent flow at Reynolds numbers 5 . 104 to 106 around the NREL S826 airfoil used for wind turbine blades is simulated using delayed detached eddy simulation (DDES). The 3D domain is built as a replica of the low speed wind tunnel at the Norwegian University of Science and Technology (NTNU) with the wind tunnel walls considered as slip walls. The subgrid turbulent kinetic energy is used to model the sub-grid scale in the large eddy simulation (LES) part of DDES. Different Reynoldsaveraged Navier-Stokes (RANS) models are tested in ANSYS Fluent. The realizable k - ∈ model as the RANS model in DDES is found to yield the best agreement of simulated pressure distributions with the experimental data both from NTNU and the Technical University of Denmark (DTU), the latter for a shorter spanwise domain. The present DDES results are in excellent agreement with LES results from DTU. Since DDES requires much fewer cells in the RANS region near the wing surface than LES, DDES is computationally much more efficient than LES. Whereas DDES is able to predict lift and drag in close agreement with experiment up to stall, pure 2D RANS simulations fail near stall. After testing different numerical settings, time step sizes and grids for DDES, a Reynolds number study is conducted. Near stall, separated flow structures, so-called stall cells, are observed in the DDES results.

Fuels and Combustion Technology for Advanced Aircraft Engines (Les Propergols et les Systemes de Combustion pour les Moteurs d’Aeronefs)

DTIC Science & Technology

1993-09-01

geometrical ceL.•er of the expressed as follows: control volume coupled with the use of linear interpolation for internodal variation usually leads to non ...2827) fuel using the such as copper, sulfur, and nitrogen. Note that F. is a experimental data in Figs. 1 and 2. It is assumed that non -depletins...combustor concepts One case, however, exhibited a very non -uniform is that proper control of fuel-air mixing is essential for distribution of fuel liquid

A method for obtaining a statistically stationary turbulent free shear flow

NASA Technical Reports Server (NTRS)

Timson, Stephen F.; Lele, S. K.; Moser, R. D.

1994-01-01

The long-term goal of the current research is the study of Large-Eddy Simulation (LES) as a tool for aeroacoustics. New algorithms and developments in computer hardware are making possible a new generation of tools for aeroacoustic predictions, which rely on the physics of the flow rather than empirical knowledge. LES, in conjunction with an acoustic analogy, holds the promise of predicting the statistics of noise radiated to the far-field of a turbulent flow. LES's predictive ability will be tested through extensive comparison of acoustic predictions based on a Direct Numerical Simulation (DNS) and LES of the same flow, as well as a priori testing of DNS results. The method presented here is aimed at allowing simulation of a turbulent flow field that is both simple and amenable to acoustic predictions. A free shear flow is homogeneous in both the streamwise and spanwise directions and which is statistically stationary will be simulated using equations based on the Navier-Stokes equations with a small number of added terms. Studying a free shear flow eliminates the need to consider flow-surface interactions as an acoustic source. The homogeneous directions and the flow's statistically stationary nature greatly simplify the application of an acoustic analogy.

Hybrid LES/RANS simulation of a turbulent boundary layer over a rectangular cavity

NASA Astrophysics Data System (ADS)

Zhang, Qi; Haering, Sigfried; Oliver, Todd; Moser, Robert

2016-11-01

We report numerical investigations of a turbulent boundary layer over a rectangular cavity using a new hybrid RANS/LES model and the traditional Detached Eddy Simulation (DES). Our new hybrid method aims to address many of the shortcomings from the traditional DES. In the new method, RANS/LES blending controlled by a parameter that measures the ratio of the modeled subgrid kinetic energy to an estimate of the subgrid energy based on the resolved scales. The result is a hybrid method automatically resolves as much turbulence as can be supported by the grid and transitions appropriately from RANS to LES without the need for ad hoc delaying functions that are often required for DES. Further, the new model is designed to improve upon DES by accounting for the effects of grid anisotropy and inhomogeneity in the LES region. We present comparisons of the flow features inside the cavity and the pressure time history and spectra as computed using the new hybrid model and DES.

Electromagnetic modelling of conductive or superconductive microstrip lines using spice as electromagnetic solver

NASA Astrophysics Data System (ADS)

Dinh, Thanh Vu; Cabon, Béatrice; Daoud, Nahla; Chilo, Jean

1992-11-01

This paper presents a simple and efficient method for calculating the propagating line parameters (actually, a microstrip one) and its magnetic fields, by simulating an original equivalent circuit with an electrical nodal simulator (SPICE). The losses in the normal conducting line (due to DC losses and to skin effect losses) and also in the superconducting one can be investigated. This allows us to integrate the electromagnetic solutions to the CAD softwares. Dans ce papier, une méthode simple et efficace pour calculer les paramètres de propagation d'une ligne microruban et les champs magnétiques qu'elle engendre est présentée; pour cela, nous simulons un circuit original équivalent à l'aide du simulateur nodal SPICE. Les pertes dans une ligne conductrice (pertes continues et par effet de peau) ainsi que dans une ligne supraconductrice peuvent être considérées. Les solutions électromagnétiques peuvent être intégrées dans les simulateurs de CAO.

Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder

DOE PAGES

Morgan, B. E.; Greenough, J. A.

2015-04-08

Two-dimensional numerical simulations of the Richtmyer–Meshkov unstable “shock-jet” problem are conducted using both large-eddy simulation (LES) and unsteady Reynolds-averaged Navier–Stokes (URANS) approaches in an arbitrary Lagrangian–Eulerian hydrodynamics code. Turbulence statistics are extracted from LES by running an ensemble of simulations with multimode perturbations to the initial conditions. Detailed grid convergence studies are conducted, and LES results are found to agree well with both experiment and high-order simulations conducted by Shankar et al. (Phys Fluids 23, 024102, 2011). URANS results using a k–L approach are found to be highly sensitive to initialization of the turbulence lengthscale L and to the timemore » at which L becomes resolved on the computational mesh. As a result, it is observed that a gradient diffusion closure for turbulent species flux is a poor approximation at early times, and a new closure based on the mass-flux velocity is proposed for low-Reynolds-number mixing.« less

A stochastic particle method for the investigation of turbulence/chemistry interactions in large-eddy simulations of turbulent reacting flows

NASA Astrophysics Data System (ADS)

Ferrero, Pietro

The main objective of this work is to investigate the effects of the coupling between the turbulent fluctuations and the highly non-linear chemical source terms in the context of large-eddy simulations of turbulent reacting flows. To this aim we implement the filtered mass density function (FMDF) methodology on an existing finite volume (FV) fluid dynamics solver. The FMDF provides additional statistical sub-grid scale (SGS) information about the thermochemical state of the flow - species mass fractions and enthalpy - which would not be available otherwise. The core of the methodology involves solving a transport equation for the FMDF by means of a stochastic, grid-free, Lagrangian particle procedure. Any moments of the distribution can be obtained by taking ensemble averages of the particles. The main advantage of this strategy is that the chemical source terms appear in closed form so that the effects of turbulent fluctuations on these terms are already accounted for and do not need to be modeled. We first validate and demonstrate the consistency of our implementation by comparing the results of the hybrid FV/FMDF procedure against model-free LES for temporally developing, non-reacting mixing layers. Consistency requires that, for non-reacting cases, the two solvers should yield identical solutions. We investigate the sensitivity of the FMDF solution on the most relevant numerical parameters, such as the number of particles per cell and the size of the ensemble domain. Next, we apply the FMDF modeling strategy to the simulation of chemically reacting, two- and three-dimensional temporally developing mixing layers and compare the results against both DNS and model-free LES. We clearly show that, when the turbulence/chemistry interaction is accounted for with the FMDF methodology, the results are in much better agreement to the DNS data. Finally, we perform two- and three-dimensional simulations of high Reynolds number, spatially developing, chemically reacting mixing layers, with the intent of reproducing a set of experimental results obtained at the California Institute of Technology. The mean temperature rise calculated by the hybrid FV/FMDF solver, which is associated with the amount of product formed, lies very close to the experimental profile. Conversely, when the effects of turbulence/chemistry coupling are ignored, the simulations clearly over predict the amount of product that is formed.

Techniques d'inspection par ondes guidees ultrasonores d'assemblages brases dans des reacteurs aeronautiques =

NASA Astrophysics Data System (ADS)

Comot, Pierre

L'industrie aeronautique, cherche a etudier la possibilite d'utiliser de maniere structurelle des joints brases, dans une optique de reduction de poids et de cout. Le developpement d'une methode d'evaluation rapide, fiable et peu couteuse pour evaluer l'integrite structurelle des joints apparait donc indispensable. La resistance mecanique d'un joint brase dependant principalement de la quantite de phase fragile dans sa microstructure. Les ondes guidees ultrasonores permettent de detecter ce type de phase lorsqu'elles sont couplees a une mesure spatio-temporelle. De plus la nature de ce type d'ondes permet l'inspection de joints ayant des formes complexes. Ce memoire se concentre donc sur le developpement d'une technique basee sur l'utilisation d'ondes guidees ultrasonores pour l'inspection de joints brases a recouvrement d'Inconel 625 avec comme metal d'apport du BNi-2. Dans un premiers temps un modele elements finis du joint a ete utilise pour simuler la propagation des ultrasons et optimiser les parametres d'inspection, la simulation a permis egalement de demontrer la faisabilite de la technique pour la detection de la quantite de phase fragile dans ce type de joints. Les parametres optimises sont la forme de signal d'excitation, sa frequence centrale et la direction d'excitation. Les simulations ont montre que l'energie de l'onde ultrasonore transmise a travers le joint aussi bien que celle reflechie, toutes deux extraites des courbes de dispersion, etaient proportionnelles a la quantite de phase fragile presente dans le joint et donc cette methode permet d'identifier la presence ou non d'une phase fragile dans ce type de joint. Ensuite des experimentations ont ete menees sur trois echantillons typiques presentant differentes quantites de phase fragile dans le joint, pour obtenir ce type d'echantillons differents temps de brasage ont ete utilises (1, 60 et 180 min). Pour cela un banc d'essai automatise a ete developpe permettant d'effectuer une analyse similaire a celle utilisee en simulation. Les parametres experimentaux ayant ete choisis en accord avec l'optimisation effectuee lors des simulations et apres une premiere optimisation du procede experimental. Finalement les resultats experimentaux confirment les resultats obtenus en simulation, et demontrent le potentiel de la methode developpee.

Microscale anthropogenic pollution modelling in a small tropical island during weak trade winds: Lagrangian particle dispersion simulations using real nested LES meteorological fields

NASA Astrophysics Data System (ADS)

Cécé, Raphaël; Bernard, Didier; Brioude, Jérome; Zahibo, Narcisse

2016-08-01

Tropical islands are characterized by thermal and orographical forcings which may generate microscale air mass circulations. The Lesser Antilles Arc includes small tropical islands (width lower than 50 km) where a total of one-and-a-half million people live. Air quality over this region is affected by anthropogenic and volcanic emissions, or saharan dust. To reduce risks for the population health, the atmospheric dispersion of emitted pollutants must be predicted. In this study, the dispersion of anthropogenic nitrogen oxides (NOx) is numerically modelled over the densely populated area of the Guadeloupe archipelago under weak trade winds, during a typical case of severe pollution. The main goal is to analyze how microscale resolutions affect air pollution in a small tropical island. Three resolutions of domain grid are selected: 1 km, 333 m and 111 m. The Weather Research and Forecasting model (WRF) is used to produce real nested microscale meteorological fields. Then the weather outputs initialize the Lagrangian Particle Dispersion Model (FLEXPART). The forward simulations of a power plant plume showed good ability to reproduce nocturnal peaks recorded by an urban air quality station. The increase in resolution resulted in an improvement of model sensitivity. The nesting to subkilometer grids helped to reduce an overestimation bias mainly because the LES domains better simulate the turbulent motions governing nocturnal flows. For peaks observed at two air quality stations, the backward sensitivity outputs identified realistic sources of NOx in the area. The increase in resolution produced a sharper inverse plume with a more accurate source area. This study showed the first application of the FLEXPART-WRF model to microscale resolutions. Overall, the coupling model WRF-LES-FLEXPART is useful to simulate the pollutant dispersion during a real case of calm wind regime over a complex terrain area. The forward and backward simulation results showed clearly that the subkilometer resolution of 333 m is necessary to reproduce realistic air pollution patterns in this case of short-range transport over a complex terrain area. Globally, this work contributes to enrich the sparsely documented domain of real nested microscale air pollution modelling. This study dealing with the determination of the proper resolution grid and proper turbulence scheme, is of significant interest to the near-source and complex terrain air quality research community.

Simulation of heat and mass transfer in turbulent channel flow using the spectral-element method: effect of spatial resolution

NASA Astrophysics Data System (ADS)

Ryzhenkov, V.; Ivashchenko, V.; Vinuesa, R.; Mullyadzhanov, R.

2016-10-01

We use the open-source code nek5000 to assess the accuracy of high-order spectral element large-eddy simulations (LES) of a turbulent channel flow depending on the spatial resolution compared to the direct numerical simulation (DNS). The Reynolds number Re = 6800 is considered based on the bulk velocity and half-width of the channel. The filtered governing equations are closed with the dynamic Smagorinsky model for subgrid stresses and heat flux. The results show very good agreement between LES and DNS for time-averaged velocity and temperature profiles and their fluctuations. Even the coarse LES grid which contains around 30 times less points than the DNS one provided predictions of the friction velocity within 2.0% accuracy interval.

Large Eddy Simulation of a Film Cooling Flow Injected from an Inclined Discrete Cylindrical Hole into a Crossflow with Zero-Pressure Gradient Turbulent Boundary Layer

NASA Technical Reports Server (NTRS)

Johnson, Perry L.; Shyam, Vikram

2012-01-01

A Large Eddy Simulation (LES) is performed of a high blowing ratio (M = 1.7) film cooling flow with density ratio of unity. Mean results are compared with experimental data to show the degree of fidelity achieved in the simulation. While the trends in the LES prediction are a noticeable improvement over Reynolds-Averaged Navier-Stokes (RANS) predictions, there is still a lack a spreading on the underside of the lifted jet. This is likely due to the inability of the LES to capture the full range of influential eddies on the underside of the jet due to their smaller structure. The unsteady structures in the turbulent coolant jet are also explored and related to turbulent mixing characteristics

Large eddy simulations in 2030 and beyond

PubMed Central

Piomelli, U

2014-01-01

Since its introduction, in the early 1970s, large eddy simulations (LES) have advanced considerably, and their application is transitioning from the academic environment to industry. Several landmark developments can be identified over the past 40 years, such as the wall-resolved simulations of wall-bounded flows, the development of advanced models for the unresolved scales that adapt to the local flow conditions and the hybridization of LES with the solution of the Reynolds-averaged Navier–Stokes equations. Thanks to these advancements, LES is now in widespread use in the academic community and is an option available in most commercial flow-solvers. This paper will try to predict what algorithmic and modelling advancements are needed to make it even more robust and inexpensive, and which areas show the most promise. PMID:25024415

Exploring the Effects of Cloud Vertical Structure on Cloud Microphysical Retrievals based on Polarized Reflectances

NASA Astrophysics Data System (ADS)

Miller, D. J.; Zhang, Z.; Platnick, S. E.; Ackerman, A. S.; Cornet, C.; Baum, B. A.

2013-12-01

A polarized cloud reflectance simulator was developed by coupling an LES cloud model with a polarized radiative transfer model to assess the capabilities of polarimetric cloud retrievals. With future remote sensing campaigns like NASA's Aerosols/Clouds/Ecosystems (ACE) planning to feature advanced polarimetric instruments it is important for the cloud remote sensing community to understand the retrievable information available and the related systematic/methodical limitations. The cloud retrieval simulator we have developed allows us to probe these important questions in a realistically relevant test bed. Our simulator utilizes a polarized adding-doubling radiative transfer model and an LES cloud field from a DHARMA simulation (Ackerman et al. 2004) with cloud properties based on the stratocumulus clouds observed during the DYCOMS-II field campaign. In this study we will focus on how the vertical structure of cloud microphysics can influence polarized cloud effective radius retrievals. Numerous previous studies have explored how retrievals based on total reflectance are affected by cloud vertical structure (Platnick 2000, Chang and Li 2002) but no such studies about the effects of vertical structure on polarized retrievals exist. Unlike the total cloud reflectance, which is predominantly multiply scattered light, the polarized reflectance is primarily the result of singly scattered photons. Thus the polarized reflectance is sensitive to only the uppermost region of the cloud (tau~<1) where photons can scatter once and still escape before being scattered again. This means that retrievals based on polarized reflectance have the potential to reveal behaviors specific to the cloud top. For example cloud top entrainment of dry air, a major influencer on the microphysical development of cloud droplets, can be potentially studied with polarimetric retrievals.

Analysis of a Simulation Experiment on Optimized Crewing for Damage Control

DTIC Science & Technology

2012-06-01

base donnaient un rendement supérieur à l’automatisation moyenne pour l’intervention en cas d’inondation. À partir de ces analyses, les auteurs du...et l’analyse ultérieures de données aux fins d’expériences de simulation semblables. Enfin, les auteurs du rapport ont établi des pistes...DRDC Toronto. [6] Floyd, J., Hunt, S., Williams, F., & Tatem, P. (2004). Fire + Smoke Simulator (FSSIM), Version 1 - Theory manual (NRL/MR/6180-04

Toward Development of a Stochastic Wake Model: Validation Using LES and Turbine Loads

DOE PAGES

Moon, Jae; Manuel, Lance; Churchfield, Matthew; ...

2017-12-28

Wind turbines within an array do not experience free-stream undisturbed flow fields. Rather, the flow fields on internal turbines are influenced by wakes generated by upwind unit and exhibit different dynamic characteristics relative to the free stream. The International Electrotechnical Commission (IEC) standard 61400-1 for the design of wind turbines only considers a deterministic wake model for the design of a wind plant. This study is focused on the development of a stochastic model for waked wind fields. First, high-fidelity physics-based waked wind velocity fields are generated using Large-Eddy Simulation (LES). Stochastic characteristics of these LES waked wind velocity field,more » including mean and turbulence components, are analyzed. Wake-related mean and turbulence field-related parameters are then estimated for use with a stochastic model, using Multivariate Multiple Linear Regression (MMLR) with the LES data. To validate the simulated wind fields based on the stochastic model, wind turbine tower and blade loads are generated using aeroelastic simulation for utility-scale wind turbine models and compared with those based directly on the LES inflow. The study's overall objective is to offer efficient and validated stochastic approaches that are computationally tractable for assessing the performance and loads of turbines operating in wakes.« less

Toward Development of a Stochastic Wake Model: Validation Using LES and Turbine Loads

DOE Office of Scientific and Technical Information (OSTI.GOV)

Moon, Jae; Manuel, Lance; Churchfield, Matthew

Wind turbines within an array do not experience free-stream undisturbed flow fields. Rather, the flow fields on internal turbines are influenced by wakes generated by upwind unit and exhibit different dynamic characteristics relative to the free stream. The International Electrotechnical Commission (IEC) standard 61400-1 for the design of wind turbines only considers a deterministic wake model for the design of a wind plant. This study is focused on the development of a stochastic model for waked wind fields. First, high-fidelity physics-based waked wind velocity fields are generated using Large-Eddy Simulation (LES). Stochastic characteristics of these LES waked wind velocity field,more » including mean and turbulence components, are analyzed. Wake-related mean and turbulence field-related parameters are then estimated for use with a stochastic model, using Multivariate Multiple Linear Regression (MMLR) with the LES data. To validate the simulated wind fields based on the stochastic model, wind turbine tower and blade loads are generated using aeroelastic simulation for utility-scale wind turbine models and compared with those based directly on the LES inflow. The study's overall objective is to offer efficient and validated stochastic approaches that are computationally tractable for assessing the performance and loads of turbines operating in wakes.« less

Investigation of extinction and re-ignition in piloted turbulent non-premixed methane-air flames using LES and high-speed OH-LIF

NASA Astrophysics Data System (ADS)

Prasad, Vinayaka N.; Juddoo, Mrinal; Masri, Assaad R.; Jones, William P.; Luo, Kai H.

2013-06-01

Extinction and re-ignition processes observed experimentally in thin reaction zones of piloted turbulent non-premixed methane flames approaching blow-off are analysed using Large Eddy Simulation (LES) along with the Eulerian stochastic field method representing the unresolved sub-grid turbulence-chemistry interactions. Eight stochastic fields in conjunction with a reduced chemical mechanism involving 19 species are employed to perform simulations of the Sydney flames L, B and M, which exhibit increasing levels of extinction. The agreement of the flame statistics of the velocities, mixture fraction and selected reactive species were found to be encouraging and highlight the ability of the method to capture quantitatively the effects of increasing jet velocity in this series. In a subsequent analysis of the flame structure using the LES simulation data, the strong three-dimensionality of the flame was emphasised. Quantitative comparisons with recent measurements using high-speed Planar Laser-Induced Fluorescence of OH (OH-PLIF) were found to be in reasonably good agreement with LES simulations and confirm the previous observations that the rates of flame breakages are greater than those of flame closures. This study, which also represents the first successful numerical attempt to describe the entire flame series, highlights the potential and complementary capabilities of a hybrid LES and high-speed imaging approach to resolve issues such as the role of out-of-plane motion in the investigation of transient processes such as flame breakages and re-ignition.

Operation and Equivalent Loads of Wind Turbines in Large Wind Farms

NASA Astrophysics Data System (ADS)

Andersen, Soren Juhl; Sorensen, Jens Norkaer; Mikkelsen, Robert Flemming

2017-11-01

Wind farms continue to grow in size and as the technology matures, the design of wind farms move towards including dynamic effects besides merely annual power production estimates. The unsteady operation of wind turbines in large wind farms has been modelled with EllipSys3D(Michelsen, 1992, and Sørensen, 1995) for a number of different scenarios using a fully coupled large eddy simulations(LES) and aero-elastic framework. The turbines are represented in the flow fields using the actuator line method(Sørensen and Shen, 2002), where the aerodynamic forces and deflections are derived from an aero-elastic code, Flex5(Øye, 1996). The simulations constitute a database of full turbine operation in terms of both production and loads for various wind speeds, turbulence intensities, and turbine spacings. The operating conditions are examined in terms of averaged power production and thrust force, as well as 10min equivalent flapwise bending, yaw, and tilt moment loads. The analyses focus on how the performance and loads change throughout a given farm as well as comparing how various input parameters affect the operation and loads of the wind turbines during different scenarios. COMWIND(Grant 2104-09- 067216/DSF), Nordic Consortium on Optimization and Control of Wind Farms, Eurotech Greentech Wind project, Winds2Loads, and CCA LES. Ressources Granted on SNIC and JESS. The Vestas NM80 turbine has been used.

A novel approach to predict the stability limits of combustion chambers with large eddy simulation

NASA Astrophysics Data System (ADS)

Pritz, B.; Magagnato, F.; Gabi, M.

2010-06-01

Lean premixed combustion, which allows for reducing the production of thermal NOx, is prone to combustion instabilities. There is an extensive research to develop a reduced physical model, which allows — without time-consuming measurements — to calculate the resonance characteristics of a combustion system consisting of Helmholtz resonator type components (burner plenum, combustion chamber). For the formulation of this model numerical investigations by means of compressible Large Eddy Simulation (LES) were carried out. In these investigations the flow in the combustion chamber is isotherm, non-reacting and excited with a sinusoidal mass flow rate. Firstly a combustion chamber as a single resonator subsequently a coupled system of a burner plenum and a combustion chamber were investigated. In this paper the results of additional investigations of the single resonator are presented. The flow in the combustion chamber was investigated without excitation at the inlet. It was detected, that the mass flow rate at the outlet cross section is pulsating once the flow in the chamber is turbulent. The fast Fourier transform of the signal showed that the dominant mode is at the resonance frequency of the combustion chamber. This result sheds light on a very important source of self-excited combustion instabilities. Furthermore the LES can provide not only the damping ratio for the analytical model but the eigenfrequency of the resonator also.

The Interaction of the Space Shuttle Launch and Entry Suits and Sustained Weightless on Astronaut Egress Locomotion

NASA Technical Reports Server (NTRS)

Greenisen, M. C.; Bishop, P. A.; Sothmann, M.

2008-01-01

The purpose of this study was to determine the consequences of extended periods of weightlessness during space missions on astronauts f ability to perform a simulated contingency egress while wearing either of the Launch and Entry suits immediately after space flight. In our previous lab-based study of simulated contingency egress, we found only 4 of 12 non-astronauts wearing the Launch and Entry Suit (LES) successfully completed the simulated egress. However, 4 of 4 of the previous failures (when tested wearing the LES), were then successful in completing the test wearing the Advanced Crew Escape Suit (ACES). Therefore, this study tested 21 Astronaut Volunteers wearing either the LES or ACES while performing a simulated egress on a treadmill (TM) onboard the Crew Transportation Vehicle immediately after space flight at either the Kennedy Space Center or Edwards AFB. Astronauts walked for 400 meters at 1.6m/sec with g-suit inflation level set to preflight testing levels, visor down, breathing from the suit emergency O2 supply. Metabolic, heartrate, and perceived exertion data were collected during these post-flight tests. Exactly the same preflight simulated egress tests on a TM were performed in the lab at NASA/JSC by each crewmember at L-60. Preflight testing found 2 of the 21 crewmembers were unable to complete the simulated contingency egress. Postflight, 9 crew (8 ACES, 1 LES) completed the simulated contingency egress of 400 meters at 1.6m/sec. and 12 failed to meet that standard (7 ACES, 5 LES). Preflight physiological response tests failed to identify crew capable of performing the egress vs. those who failed. However, 18 of the 21 crew did make at least 2.67 minutes into the postflight egress testing. At that point in time, heartrate was higher (P <=.20) for the failures compared to the finishers. These findings indicate that NASA fs switch to the ACES for space flight crews should be expedited.

Effet des Plastifiants sur la Temperature de Transition Vitreuse, Tg, du Butadiene a Terminaison Hydroxy et du Polyazoture de Glycidyle - Simulations et Experiences (Effect of the Plasticizers on the Vitreous Transition Temperature, Tg, of the Butadiene to Termination Hydroxy and the Polyazoture of Glycidyle - Simulations and Experiences)

DTIC Science & Technology

2008-09-01

diverses temperatures 26 a) HTPB pur b) HTPB-DOA (polymere et plastifiant) c) GAP pur d) Gpl pur e)Gap-Gpl Liste des tableaux Tableau 1...Composition des mailles amorphes construites 11 Tableau 2. Proprietes des polymeres et plastifiants utilises 11 Tableau 3. Comparaisons entre les Tt...obtenues experimentalement, les T% publiees dans les ecrits scientifiques et celles predites a partir des 7"gdes composes purs 19 Tableau 4. Comparaison

Dynamic Load Predictions for Launchers Using Extra-Large Eddy Simulations X-Les

NASA Astrophysics Data System (ADS)

Maseland, J. E. J.; Soemarwoto, B. I.; Kok, J. C.

2005-02-01

Flow-induced unsteady loads can have a strong impact on performance and flight characteristics of aerospace vehicles and therefore play a crucial role in their design and operation. Complementary to costly flight tests and delicate wind-tunnel experiments, unsteady loads can be calculated using time-accurate Computational Fluid Dynamics. A capability to accurately predict the dynamic loads on aerospace structures at flight Reynolds numbers can be of great value for the design and analysis of aerospace vehicles. Advanced space launchers are subject to dynamic loads in the base region during the ascent to space. In particular the engine and nozzle experience aerodynamic pressure fluctuations resulting from massive flow separations. Understanding these phenomena is essential for performance enhancements for future launchers which operate a larger nozzle. A new hybrid RANS-LES turbulence modelling approach termed eXtra-Large Eddy Simulations (X-LES) holds the promise to capture the flow structures associated with massive separations and enables the prediction of the broad-band spectrum of dynamic loads. This type of method has become a focal point, reducing the cost of full LES, driven by the demand for their applicability in an industrial environment. The industrial feasibility of X-LES simulations is demonstrated by computing the unsteady aerodynamic loads on the main-engine nozzle of a generic space launcher configuration. The potential to calculate the dynamic loads is qualitatively assessed for transonic flow conditions in a comparison to wind-tunnel experiments. In terms of turn-around-times, X-LES computations are already feasible within the time-frames of the development process to support the structural design. Key words: massive separated flows; buffet loads; nozzle vibrations; space launchers; time-accurate CFD; composite RANS-LES formulation.

A Hybrid Numerical Method for Turbulent Mixing Layers. Degree awarded by Case Western Reserve Univ.

NASA Technical Reports Server (NTRS)

Georgiadis, Nicholas J.

2001-01-01

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modern day aircraft and also those of hypersonic vehicles currently under development. The method configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. Closure for the RANS equations was obtained using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The wall-function approach enabled a continuous computational grid from the RANS regions to the LES region. The LES equations were closed using the Smagorinsky subgrid scale model. The hybrid RANS-LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Vortex shedding from the base region of a splitter plate separating the upstream flows was observed to eventually transition to turbulence. The location of the transition, however, was much further downstream than indicated by experiments. Actual LES calculations, performed in three spatial directions, also indicated vortex shedding, but the transition to turbulence was found to occur much closer to the beginning of the mixing section. which is in agreement with experimental observations. These calculations demonstrated that LES simulations must be performed in three dimensions. Comparisons of time-averaged axial velocities and turbulence intensities indicated reasonable agreement with experimental data.

Annual Research Briefs, 1998

NASA Technical Reports Server (NTRS)

Spinks, Debra (Compiler)

1998-01-01

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Calcul numérique accéléré de mesures de performance sur les modèles Markoviens appliqués aux systèmes informatiques et de télécommunications

NASA Astrophysics Data System (ADS)

Ghaddar, A.; Sinno, N.

2005-05-01

La complexité du phénomène de files d'attente dans les systèmes informatiques et télécommunications nécessite leur simulation par des modèles Markoviens pour les mesures de performance, mesure des délais d'attente au niveau des routeurs pour le modèle informatique et l'étude de la gestion des appels téléphoniques pour le modèle des circuits téléphoniques. L'optimisation des méthodes numériques de résolution des équations relatives à ces deux modèles va permettre d' ídentifier les critères de convergence rapide vers les états stationnaires correspondant à ces mesures.

Large Eddy Simulation of a Supercritical Turbulent Mixing Layer

NASA Astrophysics Data System (ADS)

Sheikhi, Reza; Hadi, Fatemeh; Safari, Mehdi

2017-11-01

Supercritical turbulent flows are relevant to a wide range of applications such as supercritical power cycles, gas turbine combustors, rocket propulsion and internal combustion engines. Large eddy simulation (LES) analysis of such flows involves solving mass, momentum, energy and scalar transport equations with inclusion of generalized diffusion fluxes. These equations are combined with a real gas equation of state and the corresponding thermodynamic mixture variables. Subgrid scale models are needed for not only the conventional convective terms but also the additional high pressure effects arising due to the nonlinearity associated with generalized diffusion fluxes and real gas equation of state. In this study, LES is carried out to study the high pressure turbulent mixing of methane with carbon dioxide in a temporally developing mixing layer under supercritical condition. LES results are assessed by comparing with data obtained from direct numerical simulation (DNS) of the same layer. LES predictions agree favorably with DNS data and represent several key supercritical turbulent flow features such as high density gradient regions. Supported by DOE Grant SC0017097; computational support is provided by DOE National Energy Research Scientific Computing Center.

Velocity-Resolved LES (VR-LES) technique for simulating turbulent transport of high Schmidt number passive scalars

NASA Astrophysics Data System (ADS)

Verma, Siddhartha; Blanquart, Guillaume; P. K. Yeung Collaboration

2011-11-01

Accurate simulation of high Schmidt number scalar transport in turbulent flows is essential to studying pollutant dispersion, weather, and several oceanic phenomena. Batchelor's theory governs scalar transport in such flows, but requires further validation at high Schmidt and high Reynolds numbers. To this end, we use a new approach with the velocity field fully resolved, but the scalar field only partially resolved. The grid used is fine enough to resolve scales up to the viscous-convective subrange where the decaying slope of the scalar spectrum becomes constant. This places the cutoff wavenumber between the Kolmogorov scale and the Batchelor scale. The subgrid scale terms, which affect transport at the supergrid scales, are modeled under the assumption that velocity fluctuations are negligible beyond this cutoff wavenumber. To ascertain the validity of this technique, we performed a-priori testing on existing DNS data. This Velocity-Resolved LES (VR-LES) technique significantly reduces the computational cost of turbulent simulations of high Schmidt number scalars, and yet provides valuable information of the scalar spectrum in the viscous-convective subrange.

Large-eddy simulations of mechanical and thermal processes within boundary layer of the Graciosa Island

NASA Astrophysics Data System (ADS)

Sever, G.; Collis, S. M.; Ghate, V. P.

2017-12-01

Three-dimensional numerical experiments are performed to explore the mechanical and thermal impacts of Graciosa Island on the sampling of oceanic airflow and cloud evolution. Ideal and real configurations of flow and terrain are planned using high-resolution, large-eddy resolving (e.g., Δ < 100 meter) simulations. Ideal configurations include model initializations with ideal dry and moist temperature and wind profiles to capture flow features over an island-like topography. Real configurations will use observations from different climatological background states over the Eastern Northern Atlantic, Atmospheric Radiation Measurement (ENA-ARM) site on Graciosa Island. Initial small-domain large-eddy simulations (LES) of dry airflow produce cold-pool formation upstream of an ideal two-kilometer island, with von Kármán like vortices propagation downstream. Although the peak height of Graciosa is less than half kilometer, the Azores island chain has a mountain over 2 km, which may be leading to more complex flow patterns when simulations are extended to a larger domain. Preliminary idealized low-resolution moist simulations indicate that the cloud field is impacted due to the presence of the island. Longer simulations that are performed to capture diurnal evolution of island boundary layer show distinct land/sea breeze formations under quiescent flow conditions. Further numerical experiments are planned to extend moist simulations to include realistic atmospheric profiles and observations of surface fluxes coupled with radiative effects. This work is intended to produce a useful simulation framework coupled with instruments to guide airborne and ground sampling strategies during the ACE-ENA field campaign which is aimed to better characterize marine boundary layer clouds.

A resolvable subfilter-scale model specific to large-eddy simulation of under-resolved turbulence

NASA Astrophysics Data System (ADS)

Zhou, Yong; Brasseur, James G.; Juneja, Anurag

2001-09-01

Large-eddy simulation (LES) of boundary-layer flows has serious deficiencies near the surface when a viscous sublayer either does not exist (rough walls) or is not practical to resolve (high Reynolds numbers). In previous work, we have shown that the near-surface errors arise from the poor performance of algebraic subfilter-scale (SFS) models at the first several grid levels, where integral scales are necessarily under-resolved and the turbulence is highly anisotropic. In under-resolved turbulence, eddy viscosity and similarity SFS models create a spurious feedback loop between predicted resolved-scale (RS) velocity and modeled SFS acceleration, and are unable to simultaneously capture SFS acceleration and RS-SFS energy flux. To break the spurious coupling in a dynamically meaningful manner, we introduce a new modeling strategy in which the grid-resolved subfilter velocity is estimated from a separate dynamical equation containing the essential inertial interactions between SFS and RS velocity. This resolved SFS (RSFS) velocity is then used as a surrogate for the complete SFS velocity in the SFS stress tensor. We test the RSFS model by comparing LES of highly under-resolved anisotropic buoyancy-generated homogeneous turbulence with a corresponding direct numerical simulation (DNS). The new model successfully suppresses the spurious feedback loop between RS velocity and SFS acceleration, and greatly improves model predictions of the anisotropic structure of SFS acceleration and resolved velocity fields. Unlike algebraic models, the RSFS model accurately captures SFS acceleration intensity and RS-SFS energy flux, even during the nonequilibrium transient, and properly partitions SFS acceleration between SFS stress divergence and SFS pressure force.

Quantifying phytoplankton productivity and photoinhibition in the Ross Sea Polynya with large eddy simulation of Langmuir circulation

NASA Astrophysics Data System (ADS)

Smyth, Robyn L.; Akan, Cigdem; Tejada-Martínez, Andrés.; Neale, Patrick J.

2017-07-01

Southern Ocean phytoplankton assemblages acclimated to low-light environments that result from deep mixing are often sensitive to ultraviolet and high photosynthetically available radiation. In such assemblages, exposures to inhibitory irradiance near the surface result in loss of photosynthetic capacity that is not rapidly recovered and can depress photosynthesis after transport below depths penetrated by inhibitory irradiance. We used a coupled biophysical modeling approach to quantify the reduction in primary productivity due to photoinhibition based upon experiments and observations made during the spring bloom in Ross Sea Polynya (RSP). Large eddy simulation (LES) was used to generate depth trajectories representative of observed Langmuir circulation that were passed through an underwater light field to yield time series of spectral irradiance representative of what phytoplankton would have experienced in situ. These were used to drive an assemblage-specific photosynthesis-irradiance model with inhibition determined from a biological weighting function and repair rate estimated from shipboard experiments on the local assemblage. We estimate the daily depth-integrated productivity was 230 mmol C m-2. This estimate includes a 6-7% reduction in daily depth-integrated productivity over potential productivity (i.e., effects of photoinhibition excluded). When trajectory depths were fixed (no vertical transport), the reduction in productivity was nearly double. Relative to LES estimates, there was slightly less depth-integrated photoinhibition with random walk trajectories and nearly twice as much with circular rotations. This suggests it is important to account for turbulence when simulating the effects of vertical mixing on photoinhibition due to the kinetics of photodamage and repair.

Simulating effects of a wind-turbine array using LES and RANS: Simulating turbines using LES and RANS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vanderwende, Brian J.; Kosović, Branko; Lundquist, Julie K.

2016-08-27

Growth in wind power production has motivated investigation of wind-farm impacts on in situ flow fields and downstream interactions with agriculture and other wind farms. These impacts can be simulated with both large-eddy simulations (LES) and mesoscale wind-farm parameterizations (WFP). The Weather Research and Forecasting (WRF) model offers both approaches. We used the validated generalized actuator disk (GAD) parameterization in WRF-LES to assess WFP performance. A 12-turbine array was simulated using the GAD model and the WFP in WRF. We examined the performance of each scheme in both convective and stable conditions. The GAD model and WFP produced qualitatively similarmore » wind speed deficits and turbulent kinetic energy (TKE) production across the array in both stability regimes, though the magnitudes of velocity deficits and TKE production levels were underestimated and overestimated, respectively. While wake growth slowed in the latter half of the WFP array as expected, wakes did not approach steady state by the end of the array as simulated by the GAD model. A sensitivity test involving the deactivation of explicit TKE production by the WFP resulted in turbulence levels within the array well that were below those produced by the GAD in both stable and unstable conditions. Finally, the WFP overestimated downwind power production deficits in stable conditions because of the lack of wake stabilization in the latter half of the array.« less

Stochasticity of convection in Giga-LES data

NASA Astrophysics Data System (ADS)

De La Chevrotière, Michèle; Khouider, Boualem; Majda, Andrew J.

2016-09-01

The poor representation of tropical convection in general circulation models (GCMs) is believed to be responsible for much of the uncertainty in the predictions of weather and climate in the tropics. The stochastic multicloud model (SMCM) was recently developed by Khouider et al. (Commun Math Sci 8(1):187-216, 2010) to represent the missing variability in GCMs due to unresolved features of organized tropical convection. The SMCM is based on three cloud types (congestus, deep and stratiform), and transitions between these cloud types are formalized in terms of probability rules that are functions of the large-scale environment convective state and a set of seven arbitrary cloud timescale parameters. Here, a statistical inference method based on the Bayesian paradigm is applied to estimate these key cloud timescales from the Giga-LES dataset, a 24-h large-eddy simulation (LES) of deep tropical convection (Khairoutdinov et al. in J Adv Model Earth Syst 1(12), 2009) over a domain comparable to a GCM gridbox. A sequential learning strategy is used where the Giga-LES domain is partitioned into a few subdomains, and atmospheric time series obtained on each subdomain are used to train the Bayesian procedure incrementally. Convergence of the marginal posterior densities for all seven parameters is demonstrated for two different grid partitions, and sensitivity tests to other model parameters are also presented. A single column model simulation using the SMCM parameterization with the Giga-LES inferred parameters reproduces many important statistical features of the Giga-LES run, without any further tuning. In particular it exhibits intermittent dynamical behavior in both the stochastic cloud fractions and the large scale dynamics, with periods of dry phases followed by a coherent sequence of congestus, deep, and stratiform convection, varying on timescales of a few hours consistent with the Giga-LES time series. The chaotic variations of the cloud area fractions were captured fairly well both qualitatively and quantitatively demonstrating the stochastic nature of convection in the Giga-LES simulation.

Nested high-resolution large-eddy simulations in WRF to support wind power

NASA Astrophysics Data System (ADS)

Mirocha, J.; Kirkil, G.; Kosovic, B.; Lundquist, J. K.

2009-12-01

The WRF model’s grid nesting capability provides a potentially powerful framework for simulating flow over a wide range of scales. One such application is computation of realistic inflow boundary conditions for large eddy simulations (LES) by nesting LES domains within mesoscale domains. While nesting has been widely and successfully applied at GCM to mesoscale resolutions, the WRF model’s nesting behavior at the high-resolution (Δx < 1000m) end of the spectrum is less well understood. Nesting LES within msoscale domains can significantly improve turbulent flow prediction at the scale of a wind park, providing a basis for superior site characterization, or for improved simulation of turbulent inflows encountered by turbines. We investigate WRF’s grid nesting capability at high mesh resolutions using nested mesoscale and large-eddy simulations. We examine the spatial scales required for flow structures to equilibrate to the finer mesh as flow enters a nest, and how the process depends on several parameters, including grid resolution, turbulence subfilter stress models, relaxation zones at nest interfaces, flow velocities, surface roughnesses, terrain complexity and atmospheric stability. Guidance on appropriate domain sizes and turbulence models for LES in light of these results is provided This work is performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344 LLNL-ABS-416482

Effects of Submesoscale Turbulence on Reactive Tracers in the Upper Ocean

NASA Astrophysics Data System (ADS)

Smith, Katherine Margaret

In this dissertation, Large Eddy Simulations (LES) are used to model the coupled turbulence-reactive tracer dynamics within the upper mixed layer of the ocean. Prior work has shown that LES works well over the spatial and time scales relevant to both turbulence and reactive biogeochemistry. Additionally, the code intended for use is able to carry an arbitrary number of tracer equations, allowing for easy expansion of the species reactions. Research in this dissertation includes a study of 15 idealized non-reactive tracers within an evolving large-scale temperature front in order determine and understand the fundamental dynamics underlying turbulence-tracer interaction in the absence of reactions. The focus of this study, in particular, was on understanding the evolution of biogeochemically-relevant, non-reactive tracers in the presence of both large ( 5 km) submesoscale eddies and smallscale ( 100 m) wave-driven Langmuir turbulence. The 15 tracers studied have different initial, boundary, and source conditions and significant differences are seen in their distributions depending on these conditions. Differences are also seen between regions where submesoscale eddies and small-scale Langmuir turbulence are both present, and in regions with only Langmuir turbulence. A second study focuses on the examination of Langmuir turbulence effects on upper ocean carbonate chemistry. Langmuir mixing time scales are similar to those of chemical reactions, resulting in potentially strong tracer-flow coupling effects. The strength of the Langmuir turbulence is varied, from no wave-driven turbulence (i.e., only shear-driven turbulence), to Langmuir turbulence that is much stronger than that found in typical upper ocean conditions. Three different carbonate chemistry models are also used in this study: time-dependent chemistry, equilibrium chemistry, and no-chemistry (i.e., non-reactive tracers). The third and final study described in this dissertation details the development of a reduced-order biogeochemical model with 17 state equations that can accurately reproduce the Bermuda Atlantic Time-series Study (BATS) ecosystem behavior, but that can also be integrated within high-resolution LES.

Wildfire simulation using LES with synthetic-velocity SGS models

NASA Astrophysics Data System (ADS)

McDonough, J. M.; Tang, Tingting

2016-11-01

Wildland fires are becoming more prevalent and intense worldwide as climate change leads to warmer, drier conditions; and large-eddy simulation (LES) is receiving increasing attention for fire spread predictions as computing power continues to improve (see, e.g.,). We report results from wildfire simulations over general terrain employing implicit LES for solution of the incompressible Navier-Stokes (N.-S.) and thermal energy equations with Boussinesq approximation, altered with Darcy, Forchheimer and Brinkman extensions, to represent forested regions as porous media with varying (in both space and time) porosity and permeability. We focus on subgrid-scale (SGS) behaviors computed with a synthetic-velocity model, a discrete dynamical system, based on the poor man's N.-S. equations and investigate the ability of this model to produce fire whirls (tornadoes of fire) at the (unresolved) SGS level. Professor, Mechanical Engineering and Mathematics.

Parametric modeling studies of turbulent non-premixed jet flames with thin reaction zones

NASA Astrophysics Data System (ADS)

Wang, Haifeng

2013-11-01

The Sydney piloted jet flame series (Flames L, B, and M) feature thinner reaction zones and hence impose greater challenges to modeling than the Sanida Piloted jet flames (Flames D, E, and F). Recently, the Sydney flames received renewed interest due to these challenges. Several new modeling efforts have emerged. However, no systematic parametric modeling studies have been reported for the Sydney flames. A large set of modeling computations of the Sydney flames is presented here by using the coupled large eddy simulation (LES)/probability density function (PDF) method. Parametric studies are performed to gain insight into the model performance, its sensitivity and the effect of numerics.

The coupling of high-speed high resolution experimental data and LES through data assimilation techniques

NASA Astrophysics Data System (ADS)

Harris, S.; Labahn, J. W.; Frank, J. H.; Ihme, M.

2017-11-01

Data assimilation techniques can be integrated with time-resolved numerical simulations to improve predictions of transient phenomena. In this study, optimal interpolation and nudging are employed for assimilating high-speed high-resolution measurements obtained for an inert jet into high-fidelity large-eddy simulations. This experimental data set was chosen as it provides both high spacial and temporal resolution for the three-component velocity field in the shear layer of the jet. Our first objective is to investigate the impact that data assimilation has on the resulting flow field for this inert jet. This is accomplished by determining the region influenced by the data assimilation and corresponding effect on the instantaneous flow structures. The second objective is to determine optimal weightings for two data assimilation techniques. The third objective is to investigate how the frequency at which the data is assimilated affects the overall predictions. Graduate Research Assistant, Department of Mechanical Engineering.

Filter and Grid Resolution in DG-LES

NASA Astrophysics Data System (ADS)

Miao, Ling; Sammak, Shervin; Madnia, Cyrus K.; Givi, Peyman

2017-11-01

The discontinuous Galerkin (DG) methodology has proven very effective for large eddy simulation (LES) of turbulent flows. Two important parameters in DG-LES are the grid resolution (h) and the filter size (Δ). In most previous work, the filter size is usually set to be proportional to the grid spacing. In this work, the DG method is combined with a subgrid scale (SGS) closure which is equivalent to that of the filtered density function (FDF). The resulting hybrid scheme is particularly attractive because a larger portion of the resolved energy is captured as the order of spectral approximation increases. Different cases for LES of a three-dimensional temporally developing mixing layer are appraised and a systematic parametric study is conducted to investigate the effects of grid resolution, the filter width size, and the order of spectral discretization. Comparative assessments are also made via the use of high resolution direct numerical simulation (DNS) data.

Large-Eddy Simulation of Chemically Reactive Pollutant Transport from a Point Source in Urban Area

NASA Astrophysics Data System (ADS)

Du, Tangzheng; Liu, Chun-Ho

2013-04-01

Most air pollutants are chemically reactive so using inert scalar as the tracer in pollutant dispersion modelling would often overlook their impact on urban inhabitants. In this study, large-eddy simulation (LES) is used to examine the plume dispersion of chemically reactive pollutants in a hypothetical atmospheric boundary layer (ABL) in neutral stratification. The irreversible chemistry mechanism of ozone (O3) titration is integrated into the LES model. Nitric oxide (NO) is emitted from an elevated point source in a rectangular spatial domain doped with O3. The LES results are compared well with the wind tunnel results available in literature. Afterwards, the LES model is applied to idealized two-dimensional (2D) street canyons of unity aspect ratio to study the behaviours of chemically reactive plume over idealized urban roughness. The relation among various time scales of reaction/turbulence and dimensionless number are analysed.

Investigations of Flow Over a Hemisphere Using Numerical Simulations (Postprint)

DTIC Science & Technology

2015-06-22

ranging from missile defense, remote sensing , and imaging . An important aspect of these applications is determining the effective beam-on-target...Stokes (URANS), detached eddy simulation (DES), and hybrid RANS/LES. The numerical results were compared with the experiment conducted at Auburn...turret. Using the DES and hybrid RANS/LES turbulence models, Loci-Chem was able to capture the unsteady flow structures, such as the shear layer

Skill of ship-following large-eddy simulations in reproducing MAGIC observations across the northeast Pacific stratocumulus to cumulus transition region

NASA Astrophysics Data System (ADS)

McGibbon, J.; Bretherton, C. S.

2017-06-01

During the Marine ARM GPCI Investigation of Clouds (MAGIC) in October 2011 to September 2012, a container ship making periodic cruises between Los Angeles, CA, and Honolulu, HI, was instrumented with surface meteorological, aerosol and radiation instruments, a cloud radar and ceilometer, and radiosondes. Here large-eddy simulation (LES) is performed in a ship-following frame of reference for 13 four day transects from the MAGIC field campaign. The goal is to assess if LES can skillfully simulate the broad range of observed cloud characteristics and boundary layer structure across the subtropical stratocumulus to cumulus transition region sampled during different seasons and meteorological conditions. Results from Leg 15A, which sampled a particularly well-defined stratocumulus to cumulus transition, demonstrate the approach. The LES reproduces the observed timing of decoupling and transition from stratocumulus to cumulus and matches the observed evolution of boundary layer structure, cloud fraction, liquid water path, and precipitation statistics remarkably well. Considering the simulations of all 13 cruises, the LES skillfully simulates the mean diurnal variation of key measured quantities, including liquid water path (LWP), cloud fraction, measures of decoupling, and cloud radar-derived precipitation. The daily mean quantities are well represented, and daily mean LWP and cloud fraction show the expected correlation with estimated inversion strength. There is a -0.6 K low bias in LES near-surface air temperature that results in a high bias of 5.6 W m-2 in sensible heat flux (SHF). Overall, these results build confidence in the ability of LES to represent the northeast Pacific stratocumulus to trade cumulus transition region.