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 73 records · Page 4

Nonhydrostatic effects in numerical modeling of mesoscale convective systems and baroclinic waves

The present investigation is concerned with the role of convection upon mesoscale modeling results, particularly when the grid resolution becomes small enough that there is not a clear scale separation between the explicitly resolved circulations and the parameterized clouds. In those situations, the vertical accelerations in explicitly resolved circulations become strong enough that the hydrostatic assumption may no longer be valid. These concerns arise from interests in improving mesoscale modeling per se and in improving the subgrid-scale parameterizations in global models. The hydrostatic and the nonhydrostatic options of the Colorado State University Regional Atmospheric Modeling System were used to simulate dry gravity currents in two dimensions, using several different horizontal grid sizes. With horizontal grid intervals of 10 km or less, nonhydrostatic simulations produce wider and colder heads and weaker but wider forced updrafts than do the hydrostatic simulations. Comparing the hydrostatic and nonhydrostatic models show that the difference between the vertical mass fluxes is much less than the difference between the vertical velocities. When the grid is fine enough to resolve the head of the gravity current, horizontal convergence at the gust front extends upwards almost to the head of the cold air. Vertical mass flux in the forced updraft at the front varies with horizontal grid size mainly as a function of the height of the simulated head. For coarser grids, which do not resolve the head, vertical mass flux at all heights decreases with increasing horizontal grid size. A comparison on nonhydrostatic simulations with horizontal grid intervals of 1 km and 2 km illustrates how decreasing the grid size does not necessarily increase the intensity of the resolved circulation. The smaller grid enables the simulated gravity current to entrain a bubble of warm air behind the head, which results in a weaker circulation with a shorter head and weaker updraft.

Cohen, Charles↗

Search for subgrid scale parameterization by projection pursuit regression

The dependence of subgrid-scale stresses on variables of the resolved field is studied using direct numerical simulations of isotropic turbulence, homogeneous shear flow, and channel flow. The projection pursuit algorithm, a promising new regression tool for high-dimensional data, is used to systematically search through a large collection of resolved variables, such as components of the strain rate, vorticity, velocity gradients at neighboring grid points, etc. For the case of isotropic turbulence, the search algorithm recovers the linear dependence on the rate of strain (which is necessary to transfer energy to subgrid scales) but is unable to determine any other more complex relationship. For shear flows, however, new systematic relations beyond eddy viscosity are found. For the homogeneous shear flow, the results suggest that products of the mean rotation rate tensor with both the fluctuating strain rate and fluctuating rotation rate tensors are important quantities in parameterizing the subgrid-scale stresses. A model incorporating these terms is proposed. When evaluated with direct numerical simulation data, this model significantly increases the correlation between the modeled and exact stresses, as compared with the Smagorinsky model. In the case of channel flow, the stresses are found to correlate with products of the fluctuating strain and rotation rate tensors. The mean rates of rotation or strain do not appear to be important in this case, and the model determined for homogeneous shear flow does not perform well when tested with channel flow data. Many questions remain about the physical mechanisms underlying these findings, about possible Reynolds number dependence, and, given the low level of correlations, about their impact on modeling. Nevertheless, demonstration of the existence of causal relations between sgs stresses and large-scale characteristics of turbulent shear flows, in addition to those necessary for energy transfer, provides important insight into the relation between scales in turbulent flows.

Meneveau, C.↗

Chapman Conference on the Hydrologic Aspects of Global Climate Change, Lake Chelan, WA, June 12-14, 1990, Selected Papers

The present conference on the hydrological aspects of global climate change discusses land-surface schemes for future climate models, modeling of the land-surface boundary in climate models as a composite of independent vegetation, a land-surface hydrology parameterizaton with subgrid variability for general circulation models, and conceptual aspects of a statistical-dynamical approach to represent landscape subgrid-scale heterogeneities in atmospheric models. Attention is given to the impact of global warming on river runoff, the influence of atmospheric moisture transport on the fresh water balance of the Atlantic drainage basin, a comparison of observations and model simulations of tropospheric water vapor, and the use of weather types to disaggregate the prediction of general circulation models. Topics addressed include the potential response of an Arctic watershed during a period of global warming and the sensitivity of groundwater recharge estimates to climate variability and change.

Lettenmaier, Dennis P.↗

Coupled land surface/hydrologic/atmospheric models

The topics covered include the following: prototype land cover characteristics data base for the conterminous United States; surface evapotranspiration effects on cumulus convection and implications for mesoscale models; the use of complex treatment of surface hydrology and thermodynamics within a mesoscale model and some related issues; initialization of soil-water content for regional-scale atmospheric prediction models; impact of surface properties on dryline and MCS evolution; a numerical simulation of heavy precipitation over the complex topography of California; representing mesoscale fluxes induced by landscape discontinuities in global climate models; emphasizing the role of subgrid-scale heterogeneity in surface-air interaction; and problems with modeling and measuring biosphere-atmosphere exchanges of energy, water, and carbon on large scales.

Pielke, Roger↗

Evapotranspiration from nonuniform surfaces - A first approach for short-term numerical weather prediction

Observations of surface heterogeneity of soil moisture from scales of meters to hundreds of kilometers are discussed, and a relationship between grid element size and soil moisture variability is presented. An evapotranspiration model is presented which accounts for the variability of soil moisture, standing surface water, and vegetation internal and stomatal resistance to moisture flow from the soil. The mean values and standard deviations of these parameters are required as input to the model. Tests of this model against field observations are reported, and extensive sensitivity tests are presented which explore the importance of including subgrid-scale variability in an evapotranspiration model.

Wetzel, Peter J.↗

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

The length scale for sub-grid-scale parameterization with anisotropic resolution

Use of the Smagorinsky eddy-viscosity formulation and related schemes for subgrid-scale parameterization of large eddy simulation models requires specification of a single length scale, earlier related by Lilly to the scale of filtering and/or numerical resolution. An anisotropic integration of the Kolmogoroff enstrophy spectrum allows generalization of that relationship to anisotropic resolution. It is found that the Deardorff assumption is reasonably accurate for small anisotropies and can be simply improved for larger values.

Lilly, Douglas K.↗

A description of rainfall interception over large areas

A new scheme is developed for describing interception at spatial scales comparable to the typical resolution of climate models. The scheme is based on the Rutter model of interception and statistical description of the subgrid-scale spatial variability of canopy storage and rainfall. The interception loss simulated by the new scheme is significantly smaller than those simulated by other schemes that do not include considerations for spatial variability. The explanation of this result is partly in the enhancement of spatially averaged canopy drainage due to the large local drainage from the few buckets of large canopy storage. The relative reduction in interception loss simulated by the new scheme may explain the overestimation of interception loss by climate models that do not include the effects of spatial variability on interception processes.

Eltahir, E. A. B.↗

On the Two-Dimensionalization of Homogeneous Rotating Turbulence

Large-eddy simulation of the incompressible Navier-Stokes equations has been used to examine the long-time development of initially isotropic turbulence subjected to solid-body rotation. The simulations were carried out using a pseudo-spectral method with 128 x 128 x 512 collocation points in a computational domain that is four times larger along the rotation axis than in the other directions; subgrid-scale motions were parameterized using a spectral eddy viscosity model modified for system rotation. Simulation results show that the correlation length along the rotation am's of velocities orthogonal to the rotation vector exhibits rapid growth while the integral length-scale of velocities aligned with the rotation axis is relatively unaffected by rotation. Examination of the energy spectrum of two-dimensional, two-component motions indicates the presence of an inverse cascade of energy. System rotation also causes an alignment of vorticity along the rotation axis with relatively stronger cyclonic vorticity than anticyclonic. The onset of anisotropic effects are well characterized by Rossby numbers defined in terms of both macroscopic and microscopic quantities.

Squires, K. D.↗

The implementation and validation of improved land-surface hydrology in an atmospheric general circulation model

New land-surface hydrologic parameterizations are implemented into the NASA Goddard Institute for Space Studies (GISS) General Circulation Model (GCM). These parameterizations are: 1) runoff and evapotranspiration functions that include the effects of subgrid-scale spatial variability and use physically based equations of hydrologic flux at the soil surface and 2) a realistic soil moisture diffusion scheme for the movement of water and root sink in the soil column. A one-dimensional climate model with a complete hydrologic cycle is used to screen the basic sensitivities of the hydrological parameterizations before implementation into the full three-dimensional GCM. Results of the final simulation with the GISS GCM and the new land-surface hydrology indicate that the runoff rate, especially in the tropics, is significantly improved. As a result, the remaining components of the heat and moisture balance show similar improvements when compared to observations. The validation of model results is carried from the large global (ocean and land-surface) scale to the zonal, continental, and finally the regional river basin scales.

Johnson, Kevin D.↗

Initialization of soil-water content in regional-scale atmospheric prediction models

The purpose of this study is to demonstrate the feasibility of determining the soil-water content fields required as initial conditions for land surface components within atmospheric prediction models. This is done using a model of the hydrologic balance and conventional meteorological observations, land cover, and soils information. A discussion is presented of the subgrid-scale effects, the integration time, and the choice of vegetation type on the soil-water content patterns. Finally, comparisons are made between two The Pennsylvania State University/National Center for Atmospheric Research mesoscale model simulations, one using climatological fields and the other one using the soil-moisture fields produced by this new method.

Smith, Christopher B.↗

Inertial transfers in the helical decomposition

The directions of the triadic energy transfers assumed in a previous paper agree with the spectral closures, in a similarity range, if the exponent of the power-law energy spectrum is less than unity. The helical interactions showing a strong local energy transfer when the triad is nonlocal, sum up to a reverse cascade unless the spectrum falls off faster than a -7/3 power of the wave number. The energy cascades from each type of helical interaction are calculated for a -5/3 inertial range using the eddy damped quasinormal Markovian model. One type of interaction is responsible for 86 percent of the cascade. The contributions of the two classes of helical interactions to the subgrid-scale eddy viscosity are presented, together with the contributions from the forward and reverse cascading interactions. An application of the assumption on the triadic energy transfers to turbulence under rapid rotation gives a simple argument for the tendency toward nonlinear two-dimensionalization of the flow.

Waleffe, Fabian↗

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

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

Preliminary large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids

Large-eddy simulation (LES) has matured to the point where application to complex flows is desirable. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential. In the past two years, important building blocks have been put into place making possible a careful study of LES on unstructured grids. The first building block was an efficient mesh generator which allowed the placement of points according to smooth variation of physical length scales. This variation of length scales is in all three directions independently, which allows a large reduction in points when compared to structured-grid methods, which can only vary length scales in one direction at a time. The second building block was the development of a dynamic model appropriate for unstructured grids. The principle obstacle was the development of an unstructured-grid filtering operator. In the past year, some of the new filters developed by Jansen have been implemented into a highly parallelized finite element code based on the Galerkin/least-squares finite element method. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation for a variety of reasons. First, it is a problem of significant interest since it would be the first LES of an aircraft component. Second, this flow has been the subject of three experimental studies. The third reason for considering this flow is the variety of flow features which provide an important test of the dynamic model. Only the dynamic model can be expected to perform satisfactorily in this variety of situations: from the laminar regions where it must not modify the flow at all to the turbulent boundary layers and wake where it must represent a wide variety of subgrid-scale structures. The flow configuration we have chosen is that of Wadcock (1987) at Reynolds number based on chord Re(sub c) = u(sub infinity)c/v = 1.64 x 10(exp 6), Mach number M = 0.2, and 12 deg angle of attack.

Jansen, Kenneth↗

The GISS model of the global atmosphere

A description and numerical results are presented for a global atmospheric circulation model developed at the Goddard Institute for Space Studies (GISS). The model version described is a 9-level primitive-equation model in sigma coordinates. It includes a realistic distribution of continents, oceans and topography. Detailed calculations of energy transfer by solar and terrestrial radiation make use of cloud and water vapor fields calculated by the model. The model hydrologic cycle includes two precipitation mechanisms: large-scale supersaturation and a parameterization of subgrid-scale cumulus convection. Results are presented both from a comparison of the 13th to the 43rd days (January) of one integration with climatological statistics, and from five short-range forecasting experiments. In the extended integration, the near-equilibrium January-mean model atmosphere exhibits an energy cycle in good agreement with observational estimates, together with generally realistic zonal mean fields of winds, temperature, humidity, transports, diabatic heating, evaporation, precipitation, and cloud cover.

Somerville, R. C. J.↗

Development and analysis of prognostic equations for mesoscale kinetic energy and mesoscale (subgrid scale) fluxes for large-scale atmospheric models

Generated by landscape discontinuities (e.g., sea breezes) mesoscale circulation processes are not represented in large-scale atmospheric models (e.g., general circulation models), which have an inappropiate grid-scale resolution. With the assumption that atmospheric variables can be separated into large scale, mesoscale, and turbulent scale, a set of prognostic equations applicable in large-scale atmospheric models for momentum, temperature, moisture, and any other gaseous or aerosol material, which includes both mesoscale and turbulent fluxes is developed. Prognostic equations are also developed for these mesoscale fluxes, which indicate a closure problem and, therefore, require a parameterization. For this purpose, the mean mesoscale kinetic energy (MKE) per unit of mass is used, defined as E-tilde = 0.5 (the mean value of u'(sub i exp 2), where u'(sub i) represents the three Cartesian components of a mesoscale circulation (the angle bracket symbol is the grid-scale, horizontal averaging operator in the large-scale model, and a tilde indicates a corresponding large-scale mean value). A prognostic equation is developed for E-tilde, and an analysis of the different terms of this equation indicates that the mesoscale vertical heat flux, the mesoscale pressure correlation, and the interaction between turbulence and mesoscale perturbations are the major terms that affect the time tendency of E-tilde. A-state-of-the-art mesoscale atmospheric model is used to investigate the relationship between MKE, landscape discontinuities (as characterized by the spatial distribution of heat fluxes at the earth's surface), and mesoscale sensible and latent heat fluxes in the atmosphere. MKE is compared with turbulence kinetic energy to illustrate the importance of mesoscale processes as compared to turbulent processes. This analysis emphasizes the potential use of MKE to bridge between landscape discontinuities and mesoscale fluxes and, therefore, to parameterize mesoscale fluxes generated by such subgrid-scale landscape discontinuities in large-scale atmospheric models.

Avissar, Roni↗

The impact of land-surface wetness heterogeneity on mesoscale heat fluxes

Vertical heat fluxes associated with mesoscale circulations generated by land-surface wetness discontinuities are often stronger than turbulent fluxes, especially in the upper part of the atmospheric planetary boundary layer. As a result, they contribute significantly to the subgrid-scale fluxes in large-scale atmospheric models. Yet they are not considered in these models. To provide some insights into the possible parameterization of these fluxes in large-scale models, a state-of-the-art mesoscale numerical model was used to investigate the relationships between mesoscale heat fluxes and atmospheric and land-surface characteristics that play a key role in the generation of mesoscale circulations. The distribution of land-surface wetness, the wavenumber and the wavelength of the land-surface discontinuities, and the large-scale wind speed have a significant impact on the mesoscale heat fluxes. Empirical functions were derived to characterize the relationships between mesoscale heat fluxes and the spatial distribution of land-surface wetness. The strongest mesoscale heat fluxes were obtained for a wavelength of forcing corresponding approximately to the local Rossby deformation radius. The mesoscale heat fluxes are weakened by large-scale background winds but remain significant even with moderate winds.

Chen, Fei↗