Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subgrid-Scale 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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Improved subgrid-scale models for large-eddy simulation

The paper analyzes models for subgrid-scale turbulence. The analysis indicates that there is sufficient information in the resolved scales to determine some characteristics of the complete flow field. The kinetic energy of the small-scale motions can be decomposed into two parts: one results from the large scales and is correlated with them, and the other part is uncorrelated which leads to a two-component eddy-viscosity model. The 'production equals dissipation' argument does not apply to the small scales in the decay of turbulence because it does not account for the uncorrelated component. The exchange between the large and small scales takes place mainly between the smallest scales of the former and the largest scales of the latter; this argument is the basis of a new model shown to be superior to the Smagorinsky model (1963).

Bardina, J.↗
NASA NTRS1987

On the subgrid-scale modeling of compressible turbulence

A subgrid-scale model recently derived for use in the large-eddy simulation of compressible turbulent flows is examined from a fundamental theoretical and computational standpoint. It is demonstrated that this model, which is applicable only to compressible turbulent flows in the limit of small density fluctuations, correlates somewhat poorly with the results of direct numerical simulations of compressible isotropic turbulence at low Mach numbers. An alternative model, based on Favre-filtered fields, is suggested which appears to reduce these limitations.

Speziale, C. G.↗
NASA NTRS1979

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.↗
NASA NTRS1980

Tests of subgrid-scale models in strained turbulence

Strained and sheared turbulence is computed by direct simulation and it is shown that the results are in good qualitative agreement with experiments. It is found that after large amounts of strain have been applied to turbulence, the energy flow to the small scales is reduced and, in some cases, reversed. Eddy viscosity models are shown to be very poor in strained turbulence and, when they are used, the mean strain should not be included in them. Finally, new models proposed by Bardina et al. have been tested and found to offer considerable promise for the future.

Mcmillan, O. J.↗
NASA NTRS1987

Toward the large-eddy simulations of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model (in terms of Favre-filtered fields) is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The three dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 to the third power grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. Future applications of these compressible subgrid-scale models to the large-eddy simulation of supersonic aerodynamic flows are discussed briefly.

Erlebacher, G.↗
NASA NTRS1977

Evaluation of subgrid-scale turbulence models using a fully simulated turbulent flow

An exact turbulent flow field was calculated on a three-dimensional grid with 64 points on a side. The flow simulates grid-generated turbulence from wind tunnel experiments. In this simulation, the grid spacing is small enough to include essentially all of the viscous energy dissipation, and the box is large enough to contain the largest eddy in the flow. The method is limited to low-turbulence Reynolds numbers, in our case R sub lambda = 36.6. To complete the calculation using a reasonable amount of computer time with reasonable accuracy, a third-order time-integration scheme was developed which runs at about the same speed as a simple first-order scheme. It obtains this accuracy by saving the velocity field and its first-time derivative at each time step. Fourth-order accurate space-differencing is used.

Clark, R. A.↗
DOE OSTI2026

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗
NASA NTRS1987

Numerical simulation of homogeneous, isotropic, compressible turbulence

A new numerical algorithm is developed and applied to the direct simulation of compressible, homogeneous turbulent flows at low Mach numbers. A split method in time first solves a subset of the equations explicitly, followed by an implicit treatment of the pressure terms. The definition of an average sound speed over the entire field allows the implicit equations to be solved analytically, while relaxing the severe time limit imposed by the large sound speeds. Results from direct simulations on 96(3) grids provide a data base against which a new subgrid-scale model for compressible homogeneous turbulence is tested. This model reduces to the linear combination model in the absence of compressibility.

Erlebacher, G.↗
NASA NTRS1981

Trends and pacing items in computational aerodynamics

A perspective is presented of trends in computational aerodynamics, and of important technology development items that pace future advanced applications. From a survey of AIAA Journal papers published during the past two decades, the growth trends and the progressively increasing emphasis on code development for viscous, compressible, turbulent flow are illustrated. These trends are reflected in the chronology of introduction by the aerospace industry of new computational methods in aircraft design. Key pacing items outlined are: automatic grid generation for nonlinear inviscid computations; advanced computers, improved efficiency of numerical methods, and improved turbulence models for Reynolds-averaged Navier-Stokes computations; advanced computers, time-dependent three-dimensional law-of-the-wall, code development, improved efficiency of numerical methods, and improved subgrid-scale turbulence modeling for large eddy simulations.

Chapman, D. R.↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Simulation of mesoscale convective response

Results are reported from a transfer of thermodynamic and dynamic data downscale from a two-dimensional mesoscale model to a two-dimensional cloud model and a subsequent examination of the differential convective response of the cloud model to a mesoscale structure. The mesoscale model included a high resolution PBL formulation, with convective forces expressed in a profile of exchange coefficients over the height of the PBL and the magnitude of the surface heat flux. Account was also taken of forcing by long- and short-wave radiation, surface forcing by the surface energy budget, the terrain shape, and possible wave reflection at the top boundary. The slab-symmetric cloud model possessed subgrid-scale features, five types of moisture terms, and a subroutine for accretion processes. The mesoscale environment was observed to have a significant impact on convective response, i.e., the vertical velocity and the cloud water. Various factors which were not included in the study and which must be considered are discussed.

Mcnider, R. T.↗
NASA NTRS1987

Mesoscale lake-effect snowstorms in the vicinity of Lake Michigan - Linear theory and numerical simulations

Mesoscale lake-effect snowstorms in the vicinity of Lake Michigan are studied by a linear steady-state analytic model and a nonlinear time-dependent numerical model with parameterized subgrid-scale physics. The solutions of the linear model show that the orientation of the mean wind field to the surface heating pattern is crucial to the shapes of the disturbances. The results indicate that the relative warmth of the lake surface can induce three updraft centers under a westerly wind, two updraft centers/bands under a northwesterly wind, and a convergence band under a northerly wind. Such convergences are caused by the interaction between the mean wind and the local circulations forced by the curved thermal contrasts. The numerical results from the nonlinear model not only produce convergence patterns very similar to those found in the linear theory and in other numerical studies, but also capture the transient property in some of the lake-induced disturbances. All of these results are qualitatively confirmed by satellite images.

Hsu, Hsiao-Ming↗
NASA NTRS1984

A Vertically Resolved Planetary Boundary Layer

Increase of the vertical resolution of the GLAS Fourth Order General Circulation Model (GCM) near the Earth's surface and installation of a new package of parameterization schemes for subgrid-scale physical processes were sought so that the GLAS Model GCM will predict the resolved vertical structure of the planetary boundary layer (PBL) for all grid points.

Helfand, H. M.↗
DOE OSTI2025

Ship‐Based Lidar Evaluation of Southern Ocean Low Clouds in the Storm‐Resolving General Circulation Model ICON and the ERA5 and MERRA‐2 Reanalyses

Global storm resolving models (GSRMs) represent the next generation of global climate models. One of them is a 5-km Icosahedral Nonhydrostatic Weather and Climate Model (ICON). Its high resolution means that parameterizations of convection and clouds, including subgrid-scale clouds, are omitted, relying on explicit simulation but necessarily utilizing microphysics and turbulence parameterizations. Standard-resolution (10–100 km) models, which use convection and cloud parameterizations, have substantial cloud biases over the Southern Ocean (SO), adversely affecting radiation and sea surface temperature. The SO is dominated by low clouds, which cannot be observed accurately from space due to overlapping clouds, attenuation, and ground clutter. We evaluated SO clouds in ICON and the ERA5 and MERRA-2 reanalyzes using approximately 2400 days of lidar observations and 2300 radiosonde profiles from 31 voyages and a Macquarie Island station during 2010–2021, compared to the model and reanalyzes using a ground-based lidar simulator. We found that ICON and the reanalyzes underestimate the total cloud fraction by about 10% and 20%, respectively. ICON and ERA5 overestimate the cloud occurrence peak at about 500 m, associated with underestimated lower tropospheric stability and overestimated lifting condensation level. The reanalyzes strongly underestimate fog and very low-level clouds, and MERRA-2 underestimates cloud occurrence at almost all heights. Outgoing shortwave radiation is overestimated in MERRA-2, implying a “too few, too bright” cloud problem. SO cloud and fog biases are a substantial issue in the analyzed model and reanalyzes and result in shortwave and longwave radiation biases.

Kuma, Peter [Stockholm Univ. (Sweden); Univ. of Ca↗
DOE OSTI2025

Understanding drivers and uncertainty in projected African precipitation

We investigate the drivers of projected summer precipitation changes and their uncertainties across Africa in the second half of the 21st century under the SSP2-4.5 and SSP5-8.5 scenarios using CMIP6 models. Our results reveal distinct regional precipitation changes, particularly under SSP5-8.5, with robust increases of 75%, 24%, and 17% over the Sahara, South Eastern Africa, and South Central Africa, respectively, and a decline of up to 5% over West Southern Africa (WSAF). In most regions, precipitation increases are driven by enhanced vertical thermodynamic processes associated with temperature-induced moisture increases and enhanced moisture convergence. In contrast, the WSAF decrease is associated with vertical dynamic processes driven by a weakening of the Hadley circulation’s ascending branch. Model uncertainty accounts for over 85% of total projection uncertainty across all regions and is largely due to subgrid-scale parameterizations. Overall, this study enhances our understanding of climate change impacts on African precipitation.

Taguela, Thierry N↗
NASA NTRS1976

Spectral and spatial energetics of the GISS model atmosphere

A calculation of the spectral and spatial energetics of the Goddard Institute for Space Studies (GISS) 9-level model of the global atmosphere is given. Results are presented from a model simulation of January 1973 and compared with previous model studies and an observational analysis for the Northern Hemisphere. The spectral analysis permits detailed study of those physical processes of scale larger than the grid spacing and allows indirect conclusions concerning unresolved subgrid-scale processes. Major conclusions for larger-than-grid-scale processes are that (1) the GISS model with a relatively coarse resolution of 4 deg latitude by 5 deg longitude produces significant nonlinear energy transfers; and (2) baroclinic conversion is large near planetary wave-number 6, but both the GISS model and another orographic model show a less peaked structure than nonorographic models and observations. A linear fit to the kinetic energy spectrum yields a power law of -2.6 in the wavenumber region 8 to 15. For subgrid-scale processes, the turbulence theory result of very small energy flux cascade is obtained.

Tenenbaum, J.↗
DOE OSTI2025

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗
NASA NTRS1977

Numerical simulation of the Gulf Stream and Mid-Ocean eddies

Oceanographic field measurements have revealed intense, transient mesoscale motions in many parts of the world ocean. The circulation of the western North Atlantic, considered in the present study as a rectangular basin, is simulated with a primitive equation model that has five levels and a horizontal grid size of 37 km. The model ocean is driven by a 2.5 gyre pattern of steady zonal wind stress and by a Newtonian-type surface heating. Two cases are considered: the first uses a Laplacian formulation for the subgrid-scale lateral diffusions of heat and momentum; the second uses a highly scale-selective biharmonic formulation for these diffusions. An analysis of the heat transport, in the biharmonic experiment, shows that the horizontal transport of heat by eddies is much larger than the subgrid-scale horizontal heat diffusion. In the Gulf Stream region, the eddy heat transport is comparable to the effect of a lateral diffusion coefficient of 10 to the 7th sq cm/s.

Semtner, A. J., Jr.↗