Search NASASearch

SEARCH · Search NASA

Results for “subgrid”

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.

At least 37 records · Page 2

Large-Eddy Simulations with Explicit Equations for Subgrid-Scale Quantities

Various alternative formulations of the LES equations have been explored in which additional evolution equations for variables such as the acceleration, the subgrid-scale stress tensor, or the subgrid-scale force are explicitly carried. Statistics of the velocity field obtained from the equation for the acceleration are shown to depend strongly on the initial conditions. This feature, which is independent of LES modeling issues, seems to prove that the velocity-acceleration formulation of the Navier-Stokes is not useful for numerical simulation. Equations for the subgrid-scale quantities appear to be much more stable. However, models required by this formulation of the LES problem still require additional study.

Carati, Daniele

Subgrid Combustion Modeling for the Next Generation National Combustion Code

In the first year of this research, a subgrid turbulent mixing and combustion methodology developed earlier at Georgia Tech has been provided to researchers at NASA/GRC for incorporation into the next generation National Combustion Code (called NCCLES hereafter). A key feature of this approach is that scalar mixing and combustion processes are simulated within the LES grid using a stochastic 1D model. The subgrid simulation approach recovers locally molecular diffusion and reaction kinetics exactly without requiring closure and thus, provides an attractive feature to simulate complex, highly turbulent reacting flows of interest. Data acquisition algorithms and statistical analysis strategies and routines to analyze NCCLES results have also been provided to NASA/GRC. The overall goal of this research is to systematically develop and implement LES capability into the current NCC. For this purpose, issues regarding initialization and running LES are also addressed in the collaborative effort. In parallel to this technology transfer effort (that is continuously on going), research has also been underway at Georgia Tech to enhance the LES capability to tackle more complex flows. In particular, subgrid scalar mixing and combustion method has been evaluated in three distinctly different flow field in order to demonstrate its generality: (a) Flame-Turbulence Interactions using premixed combustion, (b) Spatially evolving supersonic mixing layers, and (c) Temporal single and two-phase mixing layers. The configurations chosen are such that they can be implemented in NCCLES and used to evaluate the ability of the new code. Future development and validation will be in spray combustion in gas turbine engine and supersonic scalar mixing.

Menon, Suresh

Subgrid Scale Modeling in Solar Convection Simulations using the ASH Code

The turbulent solar convection zone has remained one of the most challenging and important subjects in physics. Understanding the complex dynamics in the solar con- vection zone is crucial for gaining insight into the solar dynamo problem. Many solar observatories have generated revealing data with great details of large scale motions in the solar convection zone. For example, a strong di erential rotation is observed: the angular rotation is observed to be faster at the equator than near the poles not only near the solar surface, but also deep in the convection zone. On the other hand, due to the wide range of dynamical scales of turbulence in the solar convection zone, both theory and simulation have limited success. Thus, cutting edge solar models and numerical simulations of the solar convection zone have focused more narrowly on a few key features of the solar convection zone, such as the time-averaged di erential rotation. For example, Brun & Toomre (2002) report computational finding of differential rotation in an anelastic model for solar convection. A critical shortcoming in this model is that the viscous dissipation is based on application of mixing length theory to stellar dynamics with some ad hoc parameter tuning. The goal of our work is to implement the subgrid scale model developed at CTR into the solar simulation code and examine how the differential rotation will be a affected as a result. Specifically, we implement a Smagorinsky-Lilly subgrid scale model into the ASH (anelastic spherical harmonic) code developed over the years by various authors. This paper is organized as follows. In x2 we briefly formulate the anelastic system that describes the solar convection. In x3 we formulate the Smagorinsky-Lilly subgrid scale model for unstably stratifed convection. We then present some preliminary results in x4, where we also provide some conclusions and future directions.

Young, Y.-N.

On the TFNS Subgrid Models for Liquid-Fueled Turbulent Combustion

This paper describes the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models used to emulate the major processes occurring in the turbulence-chemistry interaction. These two subgrid models are termed as LEM-like model and EUPDF-like model (Eulerian probability density function), respectively. Two-phase turbulent combustion in a single-element lean-direct-injection (LDI) combustor is calculated by employing the TFNS/LEM-like approach as well as the TFNS/EUPDF-like approach. Results obtained from the TFNS approach employing these two different subgrid models are compared with each other, along with the experimental data, followed by more detailed comparison between the results of an updated calculation using the TFNS/LEM-like model and the experimental data.

Liu, Nan-Suey

A Dynamic Nonlinear Subgrid-Scale Model for Large-Eddy Simulation of Complex Turbulent Flows

We present a new dynamic nonlinear subgrid-scale (SGS) model for large-eddy simulations (LES) and apply it to compute a flow involving pressure gradients, surface curvature and separation, for which data from a direct numerical simulation are available for comparison. The model, inspired by the triple model idea of Bardina et al. (“Improved Turbulence Models Based on Large Eddy Simulation of Homogeneous, Incompressible, Turbulent Flows,” Report No. TF-19, Thermosciences Division, Department of Mechanical Engineering, Stanford University, 1983), includes a Galilean-invariant term called the modified Leonard stress tensor, and two nonlinear terms comprised of the products of the strain-rate and rotation-rate tensors for an improved representation of the subgrid-scale dissipation, backscatter and anisotropy effects. The model does not employ any ad hoc averaging or clipping procedures, and does not require the specification of a characteristic length scale; hence, it naturally avoids the ambiguities associated with defining a proper length scale for anisotropic grids. Results from the wall-resolved LES of flow past a Gaussian bump using the new model demonstrate improved prediction of skin-friction, flow separation, mean flow profiles and turbulent quantities when compared to implicit LES as well as explicit LES using the Vreman SGS model on the same grid.

Turbulence

Evaluation of subgrid-scale models using an accurately simulated turbulent flow

A calculation of periodic homogeneous isotropic turbulence is used to simulate the experimental decay of grid turbulence. The calculation is found to match the experiment in a number of important aspects and the computed flow field is then treated as a realization of a physical turbulent flow. From this flow, a calculation is conducted of the large eddy field and the various averages of the subgrid-scale turbulence that occur in the large eddy simulation equations. These quantities are compared with the predictions of the models that are usually applied in large eddy simulation. The results show that the terms which involve the large-scale field are accurately modeled but the subgrid-scale Reynolds stresses are only moderately well modeled. It is also possible to use the method to predict the constants of the models without reference to experiment. Attempts to find improved models have not met with success.

Clark, R. A.

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej

Land surface hydrology parameterization for atmospheric general circulation models including subgrid scale spatial variability

Parameterizations are developed for the representation of subgrid hydrologic processes in atmospheric general circulation models. Reasonable a priori probability density functions of the spatial variability of soil moisture and of precipitation are introduced. These are used in conjunction with the deterministic equations describing basic soil moisture physics to derive expressions for the hydrologic processes that include subgrid scale variation in parameters. The major model sensitivities to soil type and to climatic forcing are explored.

Entekhabi, D.

Stratified turbulence in the atmosphere and oceans - A new subgrid model

Turbulence in stably and unstably stratified media is studied, and the results are used to construct a subgrid simulation (SGS) model for use in large eddy simulation (LES). It was found that, although the assumption of inertiality of the subgrid scales is reasonable for the case of unstable stratification, it ceases to be so in the case of stable stratification, where gravity removes kinetic energy from the eddies: the generation of gravity waves becomes the dominant physical process, with dissipation relegated to higher wavenumbers. Preliminary results show that the total kinetic energy dissipation length scale increases with stability, in accordance with LES results but in disagreement with Deardorff's model that suggests a decrease of all dissipation scales in presence of stratification.

Canuto, V. M.

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Horiuti, K.

Large-eddy simulation of laminar-turbulent breakdown at high speeds with dynamic subgrid-scale modeling

The laminar-turbulent breakdown of a boundary-layer flow along a hollow cylinder at Mach 4.5 is investigated with large-eddy simulation. The subgrid scales are modeled dynamically, where the model coefficients are determined from the local resolved field. The behavior of the dynamic-model coefficients is investigated through both an a priori test with direct numerical simulation data for the same case and a complete large-eddy simulation. Both formulations proposed by Germano et al. and Lilly are used for the determination of unique coefficients for the dynamic model and their results are compared and assessed. The behavior and the energy cascade of the subgrid-scale field structure are investigated at various stages of the transition process. The investigations are able to duplicate a high-speed transition phenomenon observed in experiments and explained only recently by the direct numerical simulations of Pruett and Zang, which is the appearance of 'rope-like' waves. The nonlinear evolution and breakdown of the laminar boundary layer and the structure of the flow field during the transition process were also investigated.

El-Hady, Nabil M.

The importance of mesoscale circulations generated by subgrid-scale landscape heterogeneities in general circulation models

A mesoscale atmospheric model was used to evaluate the impact of subgrid-scale landscape discontinuities on the vertical profiles of resolved temperature, moisture, and moist static energy in the planetary boundary layer (PBL) of general circulation models (GCMs). These profiles were produced with a 3D version of the model (using a horizontal grid resolution of 7.5 km and 13 vertical layers in the PBL) by averaging horizontally the various atmospheric variables over a 180 x 180 sq km domain-about the size of the horizontal domain represented by a single grid element in a GCM. They were compared to corresponding vertical profiles produced with a 1D version of the model, which simulates the PBL, as in a GCM, over a single horizontal grid element. Differences obtained between the horizontally averaged atmospheric variables produced with the 3D simulations and the 1D simulations emphasize the impact of subgrid-scale landscape discontinuities on GCM-resolved variables. Various types of landscape discontinuities, characterized by horizontal contrasts of surface wetness and size of land patches, were simulated under various background-wind conditions. Differences of temperature, specific humidity, and moist static energy as large as 4 K, 6 g/kg, and 10 kJ/kg were obtained in some cases. These differences were not affected significantly by moderate winds but were sensitive to the spatial distribution of surface wetness. These results emphasize the need to parameterize mesoscale processes induced by landscape discontinuities in GCMss.

Lynn, Barry H.

On why dynamic subgrid-scale models work

Dynamic subgrid models have proved to be remarkably successful in predicting the behavior of turbulent flows. Part of the reasons for their success are well understood. Since they are constructed to generate an effective viscosity which is proportional to some measure of the turbulent energy at the high wavenumber end of the spectrum, their eddy viscosity vanishes as the flow becomes laminar. This alone would justify their use over simpler models. But beyond this obvious advantage, which is confined to inhomogeneous and evolving flows, the reason why they also work better in simpler homogeneous cases, and how they do it without any obvious adjustable parameter, is not clear. This lack of understanding of the internal mechanisms of a useful tool is disturbing, not only as an intellectual challenge, but because it raises the doubt of whether it will work in all cases. This note is an attempt to clarify those mechanisms. We will see why dynamic models are robust and how they can get away with even comparatively gross errors in their formulations. This will suggest that they are only particular cases of a larger family of robust models, all of which would be relatively insensitive to large simplifications in the physics of the flow. We will also construct some such models, although mostly as research tools. It will turn out, however, that the standard dynamic formulation is not only robust to errors, but also behaves as if it were substantially well formulated. The details of why this is so will still not be clear at the end of this note, specially since it will be shown that the 'a priori' testing of the stresses gives, as is usual in most subgrid models, very poor results. But it will be argued that the basic reason is that the dynamic formulation mimics the condition that the total dissipation is approximately equal to the production measured at the test filter level.

Jimenez, J.

A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations

Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

Winckelmans, G. S.

A Priori Subgrid Analysis of Temporal Mixing Layers with Evaporating Droplets

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5. In the DNS, the gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. The Large Eddy Simulation (LES) equations corresponding to the DNS are first derived, and key assumptions in deriving them are first confirmed by computing the terms using the DNS database. Since LES of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be the sum of the filtered variables and a correction based on the filtered standard deviation; this correction is then computed from the Subgrid Scale (SGS) standard deviation. This model predicts the unfiltered variables at droplet locations considerably better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: the Smagorinsky approach, the Gradient model and the Scale-Similarity formulation. When the proportionality constant inherent in the SGS models is properly calculated, the Gradient and Scale-Similarity methods give results in excellent agreement with the DNS.

Okongo, Nora

Large-eddy Simulation of Stratocumulus-topped Atmospheric Boundary Layers with Dynamic Subgrid-scale Models

The objective of the present study is to evaluate the dynamic procedure in LES of stratocumulus topped atmospheric boundary layer and assess the relative importance of subgrid-scale modeling, cloud microphysics and radiation modeling on the predictions. The simulations will also be used to gain insight into the processes leading to cloud top entrainment instability and cloud breakup. In this report we document the governing equations, numerical schemes and physical models that are employed in the Goddard Cumulus Ensemble model (GCEM3D). We also present the subgrid-scale dynamic procedures that have been implemented in the GCEM3D code for the purpose of the present study.

Senocak, Inane

A Priori Subgrid Analysis of Temporal Mixing Layers with Evaporating Droplets

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5.

subgrid

Large eddy simulation of lean blow-off in swirl-stabilized flame with the subgrid dissipation concept

The lean blow-off mechanism of the premixed swirl flame is numerically investigated by large eddy simulation (LES) with the subgrid dissipation concept (SDC) combustion model. Three simulated cases cover stable, near blow-off, and transient conditions. Compared with the experiment, the LES-SDC approach captures the flow and combustion features for stable and near blow-off conditions. More importantly, the predictions of the blowoff procedure and duration agree satisfactorily with the experiment, indicating that the LES-SDC approach is a promising tool for predicting strong, unsteady turbulent combustion processes. Further, the numerical results are used to investigate the blow-off mechanism. Two stages in the blow-off procedure are specified. The first is the necking and extinction of the downstream flame surface, and the second is the shrinking of the upstream flame surface. The blow-off mechanism is well explained by the theory of stretched flame extinction. At the end of the recirculation zone, the large negative radial velocity pushes the flame to the central line. The combustion process here can be abstracted as the stretched counter-flame of the reactant-reactant configuration. The excessive flame stretch dominates the flame extinction and triggers the blow-off event. The upstream flame resists the intense stretch with the help of hot product recirculation, and the combustion here can be idealized as the counter-flame of the reactant-product configuration. The alignment of the temperature gradient and flow velocity, together with the excessive stretch, clearly indicates the tendency of flame local extinction. A Damköhler number-based blow-off criterion is raised from the mean flow strain rate and laminar flame bulk extinction strain rate.

42 ENGINEERING