SEARCH · Search NASA
Results for “eddying model”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A priori study of models for large eddy simulations of drop-laden flows
Large Eddy Simulation (LES) models are evaluated on a Direct Numerical Simulation (DNS) database representing a three-dimensional temporal mixing layer with evaporating drops.
A priori study of models for large eddy simulations of drop-laden flows
Large Eddy simulation (LES) models are presented and evaluated on a database obtained from Direct Numerical Simulation (DNS) of a three-dimensional temporal mixing layer with evaporating drops.
Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion
Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.
CGILS: Results from the First Phase of an International Project to Understand the Physical Mechanisms of Low Cloud Feedbacks in Single Column Models
1] CGILS-the CFMIP-GASS Intercomparison of Large Eddy Models (LESs) and single column models (SCMs)-investigates the mechanisms of cloud feedback in SCMs and LESs under idealized climate change perturbation. This paper describes the CGILS results from 15 SCMs and 8 LES models. Three cloud regimes over the subtropical oceans are studied: shallow cumulus, cumulus under stratocumulus, and well-mixed coastal stratus/stratocumulus. In the stratocumulus and coastal stratus regimes, SCMs without activated shallow convection generally simulated negative cloud feedbacks, while models with active shallow convection generally simulated positive cloud feedbacks. In the shallow cumulus alone regime, this relationship is less clear, likely due to the changes in cloud depth, lateral mixing, and precipitation or a combination of them. The majority of LES models simulated negative cloud feedback in the well-mixed coastal stratus/stratocumulus regime, and positive feedback in the shallow cumulus and stratocumulus regime. A general framework is provided to interpret SCM results: in a warmer climate, the moistening rate of the cloudy layer associated with the surface-based turbulence parameterization is enhanced; together with weaker large-scale subsidence, it causes negative cloud feedback. In contrast, in the warmer climate, the drying rate associated with the shallow convection scheme is enhanced. This causes positive cloud feedback. These mechanisms are summarized as the "NESTS" negative cloud feedback and the "SCOPE" positive cloud feedback (Negative feedback from Surface Turbulence under weaker Subsidence-Shallow Convection PositivE feedback) with the net cloud feedback depending on how the two opposing effects counteract each other. The LES results are consistent with these interpretations
Aspects of turbulent combustion
A complete second-order closure program for the investigation of chemically reacting turbulent flowfields is being developed. Many of the turbulence models being used have been developed and tested during previous studies of a variety of incompressible and compressible shear flows. A number of additional scalar correlations have to be modeled for chemically reacting flows, and a 'typical eddy' model has been developed for the joint probability density function for all the scalars. A simpler model in which higher-order scalar correlations are set zero is also being used. The program has been recently extended to multi-species, multi-step chemical reactions. Computations have been performed using the second-order closure program to study the turbulence-chemistry interaction in DF chemical laser flows and in a hydrogen-air diffusion flame. The results demonstrate the importance of the inclusion of the concentration and temperature fluctuation correlations in turbulent reacting flows.
A second-order closure analysis of turbulent diffusion flames
A complete second-order closure computer program for the investigation of compressible, turbulent, reacting shear layers was developed. The equations for the means and the second order correlations were derived from the time-averaged Navier-Stokes equations and contain third order and higher order correlations, which have to be modeled in terms of the lower-order correlations to close the system of equations. In addition to fluid mechanical turbulence models and parameters used in previous studies of a variety of incompressible and compressible shear flows, a number of additional scalar correlations were modeled for chemically reacting flows, and a typical eddy model developed for the joint probability density function for all the scalars. The program which is capable of handling multi-species, multistep chemical reactions, was used to calculate nonreacting and reacting flows in a hydrogen-air diffusion flame.
Comparison of Forecast and Observed Energetics
An energetics analysis scheme was developed to compare the observed kinetic energy balance over North America with that derived from forecast cyclone case. It is found that: (1) the observed and predicted kinetic energy and eddy conversion are in good qualitative agreement, although the model eddy conversion tends to be 2 to 3 times stronger than the observed values. The eddy conversion which is stronger in the 12 h forecast than in observations and may be due to several factors is studied; (2) vertical profiles of kinetic energy generation and dissipation exhibit lower and upper tropospheric maxima in both the forecast and observations; and (3) a lag in the observational analysis with the maximum in the observed kinetic energy occurring at 0000 GMT 14 January over the same region as the maximum Eddy conversion 12 h earlier is noted.
Massive separation and dynamic stall on a cusped trailing-edge airfoil
The cross-over from a predominantly attached two-dimensional flow to the bluff body form of separation is modeled via the interacting boundary layer approximation. The initial breakdown of the predominantly attached flow on a cusped trailing edge airfoil is examined using the Hilbert integral form of the unsteady interacting boundary layer equations. In addition, an interacting boundary layer technique is developed for calculating bluff body separation. This new model eliminates the severe scaling problems associated with bluff body separation through the use of a realistic inviscid eddy model, based on the infinite eddy Kirchhoff free-streamline description of separation. Brief consideration is given to the cross-over from bluff body separation to a predominantly attached flow, the extension to finite eddies and cascade flows, and the possible coupling with full, or parabolized, Navier-Stokes calculations.
Developments and applications of dynamic models for large eddy simulation of complex flows
The dynamic modeling procedure for large eddy simulation of turbulent flows is reviewed and recent developments in the theoretical aspects and applications are described. Methods for inclusion of backscatter of energy from small to large scale motions are presented. New formulations of the dynamic procedure are proposed which are optimized based on the subgrid scale flux vector or the energy dissipation rate instead of the subgrid scale stress tensor. Recent results from application of the model to forced isotropic turbulence with an inertial subrange, flow over a backward facing step at Reynolds number of 28000, and flow over a concave curved surface are presented.
A multiple-scale model for compressible turbulent flows
A multiple-scale model for compressible turbulent flows is proposed. It is assumed that turbulent eddy shocklets are formed primarily by the 'collisions' of large energetic eddies. The extra straining of the large eddy, due to their interactions with shocklets, enhances the energy cascade to smaller eddies. Model transport equations are developed for the turbulent kinetic energies and the energy transfer rates of the different scale. The turbulent eddy viscosity is determined by the total turbulent kinetic energy and the rate of energy transfer from the large scale to the small scale, which is different from the energy dissipation rate. The model coefficients in the modeled turbulent transport equations depend on the ratio of the turbulent kinetic energy of the large scale to that of the small scale, which renders the model more adaptive to the characteristics of individual flow. The model is tested against compressible free shear layers. The results agree satisfactorily with measurements.
Numerical Study of Flow Past a Circular Cylinder Using RANS, Hybrid RANS/LES and PANS Formulations
Two multiscale type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(sub epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. Solutions from these models are compared to RANS results and experimental data for a stationary and rotating cylinder. The parameter f(sub k) varies between zero and one and has the characteristic to be equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation example are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The models are compared to RANS results and experimental data for turbulent separated flows (TS) and laminar separated flows (LS) around stationary and rotating cylinders. For a stationary cylinder, the TS case is accurately simulated using the general two-equation k(sub epsilon) turbulence model (eddy-viscosity model). PAB3D accurately predicts the drag coefficient (CD), lift coefficient (CL) and the Strouhal number (St). The LS case was a challenge for the RANS computation with an eddy-viscosity turbulence model. The RANS/LES and PANS performed well and showed marked improvements over the RANS solution. The modified PANS model was the most accurate. For the rotating cylinder, the spin ratio varied from zero to one, and the PANS results were in good agreement with published experimental data. RANS/LES and PANS capture both temporal and spatial fluctuations and produce large-scale structures that do not occur in the RANS simulation. The current results show promise for the capability of PANS in simulating unsteady and complex flow phenomena.
Assessment of Turbulence-Chemistry Interaction Models in the National Combustion Code (NCC) - Part I
This paper describes the implementations of the linear-eddy model (LEM) and an Eulerian FDF/PDF model in the National Combustion Code (NCC) for the simulation of turbulent combustion. The impacts of these two models, along with the so called laminar chemistry model, are then illustrated via the preliminary results from two combustion systems: a nine-element gas fueled combustor and a single-element liquid fueled combustor.
A Priori Analysis of Subgrid-Scale Models for Large Eddy Simulations of Supercritical Binary-Species Mixing Layers
Models for large eddy simulation (LES) are assessed on a database obtained from direct numerical simulations (DNS) of supercritical binary-species temporal mixing layers. The analysis is performed at the DNS transitional states for heptane/nitrogen, oxygen/hydrogen and oxygen/helium mixing layers. The incorporation of simplifying assumptions that are validated on the DNS database leads to a set of LES equations that requires only models for the subgrid scale (SGS) fluxes, which arise from filtering the convective terms in the DNS equations. Constant-coefficient versions of three different models for the SGS fluxes are assessed and calibrated. The Smagorinsky SGS-flux model shows poor correlations with the SGS fluxes, while the Gradient and Similarity models have high correlations, as well as good quantitative agreement with the SGS fluxes when the calibrated coefficients are used.
Baroclinic Control of Southern Ocean Eddy Upwelling Near Topography
In the Southern Ocean, mesoscale eddies contribute to the upwelling of deep waters along sloping isopycnals, helping to close the upper branch of the meridional overturning circulation. Eddy energy (EE) is not uniformly distributed along the Antarctic Circumpolar Current (ACC). Instead, “hotspots” of EE that are associated with enhanced eddy-induced upwelling exist downstream of topographic features. This study shows that, in idealized eddy-resolved simulations, a topographic feature in the ACC path can enhance and localize eddy-induced upwelling. However, the upwelling systematically occurs in regions where eddies grow through baroclinic instability, rather than in regions where EE is large. Across a range of parameters, along-stream eddy growth rate is a more reliable indicator of eddy upwelling than traditional parameterizations such as eddy kinetic energy, eddy potential energy, or isopycnal slope. Ocean eddy parameterizations should consider metrics specific to the growth of baroclinic instability to accurately model eddy upwelling near topography.
Wall-Modeled Large-Eddy Simulations for High-Lift Configurations using FUN3D
Wall-modeled large-eddy simulation (WMLES) capability has recently been implemented into FUN3D, an unstructured, node-centered, finite-volume solver developed at the NASA Langley Research Center. In this paper, WMLES is assessed for two configurations that are representative for high-lift applications. The first configuration is a nominal two-dimensional multielement airfoil that has been extensively studied in the literature. WMLES solutions are computed for four angles of attack and compared with previously reported solutions. Good agreement of integrated forces, surface pressures, and boundary-layer velocity profiles is shown with available experimental data especially at lower angles of attack. WMLES solutions are also computed for the NASA High-Lift Common Research Model over a large range of angles of attack. Forces, pitching moments, and pressure distributions are favorably compared with the experimental data up to the maximum lift, including the angle of attack where the maximum lift is obtained experimentally. Eddy visualization techniques of q-criterion and density-gradient magnitude illustrate the resolved content.
Two self-similar Reynolds-stress transport models with anisotropic eddy viscosity
Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.
Modelling Unsteady Wall Pressures Beneath Turbulent Boundary Layers
As a structural entity of turbulence, hairpin vortices are believed to play a major role in developing and sustaining the turbulence process in the near wall region of turbulent boundary layers and may be regarded as the simplest conceptual model that can account for the essential features of the wall pressure fluctuations. In this work we focus on fully developed typical hairpin vortices and estimate the associated surface pressure distributions and their corresponding spectra. On the basis of the attached eddy model, we develop a representation of the overall surface pressure spectra in terms of the eddy size distribution. Instantaneous wavenumber spectra and spatial correlations are readily derivable from this representation. The model is validated by comparison of predicted wavenumber spectra and cross-correlations with existing emperical models and experimental data.