Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subgrid-scale processes”

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

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

Parameterization of Mixed Layer and Deep-Ocean Mesoscales Including Nonlinearity

In 2011, Chelton et al. carried out a comprehensive census of mesoscales using altimetry data and reached the following conclusions: "essentially all of the observed mesoscale features are nonlinear" and "mesoscales do not move with the mean velocity but with their own drift velocity," which is "the most germane of all the nonlinear metrics." Accounting for these results in a mesoscale parameterization presents conceptual and practical challenges since linear analysis is no longer usable and one needs a model of nonlinearity. A mesoscale parameterization is presented that has the following features: 1) it is based on the solutions of the nonlinear mesoscale dynamical equations, 2) it describes arbitrary tracers, 3) it includes adiabatic (A) and diabatic (D) regimes, 4) the eddy-induced velocity is the sum of a Gent and McWilliams (GM) term plus a new term representing the difference between drift and mean velocities, 5) the new term lowers the transfer of mean potential energy to mesoscales, 6) the isopycnal slopes are not as flat as in the GM case, 7) deep-ocean stratification is enhanced compared to previous parameterizations where being more weakly stratified allowed a large heat uptake that is not observed, 8) the strength of the Deacon cell is reduced. The numerical results are from a stand-alone ocean code with Coordinated Ocean-Ice Reference Experiment I (CORE-I) normal-year forcing.

Baroclinic flows; Mesoscale processes; Mixing; Sub↗

Simulation of Boundary-Layer Cumulus and Stratocumulus Clouds using a Cloud-Resolving Model With Low- and Third-Order Turbulence Closures

The effects of subgrid-scale condensation and transport become more important as the grid spacings increase from those typically used in large-eddy simulation (LES) to those typically used in cloud-resolving models (CRMs). Incorporation of these effects can be achieved by a joint probability density function approach that utilizes higher-order moments of thermodynamic and dynamic variables. This study examines how well shallow cumulus and stratocumulus clouds are simulated by two versions of a CRM that is implemented with low-order and third-order turbulence closures (LOC and TOC) when a typical CRM horizontal resolution is used and what roles the subgrid-scale and resolved-scale processes play as the horizontal grid spacing of the CRM becomes finer. Cumulus clouds were mostly produced through subgrid-scale transport processes while stratocumulus clouds were produced through both subgrid-scale and resolved-scale processes in the TOC version of the CRM when a typical CRM grid spacing is used. The LOC version of the CRM relied upon resolved-scale circulations to produce both cumulus and stratocumulus clouds, due to small subgrid-scale transports. The mean profiles of thermodynamic variables, cloud fraction and liquid water content exhibit significant differences between the two versions of the CRM, with the TOC results agreeing better with the LES than the LOC results. The characteristics, temporal evolution and mean profiles of shallow cumulus and stratocumulus clouds are weakly dependent upon the horizontal grid spacing used in the TOC CRM. However, the ratio of the subgrid-scale to resolved-scale fluxes becomes smaller as the horizontal grid spacing decreases. The subcloud-layer fluxes are mostly due to the resolved scales when a grid spacing less than or equal to 1 km is used. The overall results of the TOC simulations suggest that a 1-km grid spacing is a good choice for CRM simulation of shallow cumulus and stratocumulus.

Xu, Kuan-Man↗

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

El Nino and the Global Ocean Observing System

Until a decade ago, an often-quoted expression in oceanography is that very few observations are recorded throughout the ocean. Now, the sentiment is no longer valid in the uppermost 10% of the tropical Pacific Ocean nor at the surface of the global ocean. One of the remarkable legacies of the 1985-1994 Tropical Oceans Global Atmosphere (TOGA) Program is an in situ marine meteorological and upper oceanographic measurement array throughout the equatorial Pacific to monitor the development and maintenance of El Nino episodes. The TOGA Observing System, which initially consisted of moored- and drifting-buoy arrays, a network of commercial ships, and coastal and island stations, now includes a constellation of satellites and data-assimilating models to simulate subsurface oceanographic conditions. The El Nino and La Nina tropical Pacific Ocean observing system represents the initial phase of an integrated global ocean observing system. Remarkable improvements have been made in ocean model simulation of subsurface currents, but some problems persist. For example, the simulation of the South Equatorial Current (SEC) remains an important challenge in the 2S-2N Pacific equatorial wave guide. During El Nino the SEC at the equator is reduced and sometimes the direction is reversed, becoming eastward. Both conditions allow warm water stored in the western Pacific to invade the eastern region, creating an El Nino episode. Assimilation of data is a tenet of faith to correct simulation errors caused by deficiencies in surface fluxes (especially wind stress) and parameterizations of subgrid-scale physical processes. In the first of two numerical experiments, the Pacific SEC was simulated with and without assimilation of subsurface temperature data. Along the equator, a very weak SEC occurred throughout the eastern Pacific, independent of assimilation of data. However, as displayed in the diagram, in the western Pacific there was no satisfactory agreement between the two simulations. To help determine reliability of the simulated SEC in the western Pacific, current measurements recorded during the 9-19 October 1994 voyage of the French research vessel L'Atalante are also shown in the diagram. With data assimilation, the simulated SEC was in much better agreement with L'Atalante observations. The simulated SEC with data assimilation was far from perfect, in part because of the sparsity of subsurface temperature observations. In the next experiment, TOPEX/POSEIDON sea surface height data in combination with subsurface temperatures will be assimilated to assess further improvement of the simulation of the SEC.

Halpern, David↗

How Difficult is it to Reduce Low-Level Cloud Biases With the Higher-Order Turbulence Closure Approach in Climate Models?

Low-level clouds cover nearly half of the Earth and play a critical role in regulating the energy and hydrological cycle. Despite the fact that a great effort has been put to advance the modeling and observational capability in recent years, low-level clouds remains one of the largest uncertainties in the projection of future climate change. Low-level cloud feedbacks dominate the uncertainty in the total cloud feedback in climate sensitivity and projection studies. These clouds are notoriously difficult to simulate in climate models due to its complicated interactions with aerosols, cloud microphysics, boundary-layer turbulence and cloud dynamics. The biases in both low cloud coverage/water content and cloud radiative effects (CREs) remain large. A simultaneous reduction in both cloud and CRE biases remains elusive. This presentation first reviews the effort of implementing the higher-order turbulence closure (HOC) approach to representing subgrid-scale turbulence and low-level cloud processes in climate models. There are two HOCs that have been implemented in climate models. They differ in how many three-order moments are used. The CLUBB are implemented in both CAM5 and GDFL models, which are compared with IPHOC that is implemented in CAM5 by our group. IPHOC uses three third-order moments while CLUBB only uses one third-order moment while both use a joint double-Gaussian distribution to represent the subgrid-scale variability. Despite that HOC is more physically consistent and produces more realistic low-cloud geographic distributions and transitions between cumulus and stratocumulus regimes, GCMs with traditional cloud parameterizations outperform in CREs because tuning of this type of models is more extensively performed than those with HOCs. We perform several tuning experiments with CAM5 implemented with IPHOC in an attempt to produce the nearly balanced global radiative budgets without deteriorating the low-cloud simulation. One of the issues in CAM5-IPHOC is that cloud water content is much higher than in CAM5, which is combined with higher low-cloud coverage to produce larger shortwave CREs in some low-cloud prevailing regions. Thus, the cloud-radiative feedbacks are exaggerated there. The turning exercise is focused on microphysical parameters, which are also commonly used for tuning in climate models. The results will be discussed in this presentation.

Xu, Kuan-Man↗

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

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

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↗

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↗

Insights on TTL Dehydration Mechanisms from Microphysical Modelling of Aircraft Observations

The Tropical Tropopause Layer (TTL), a transition layer between the upper troposphere and lower stratosphere in the tropics, serves as the entryway of various trace gases into the stratosphere. Of particular interest is the transport of water vapor through the TTL, as WV is an important greenhouse gas and also plays a significant role in stratospheric chemistry by affecting polar stratospheric cloud formation and the ozone budget. While the dominant control of stratospheric water vapor by tropical cold point temperatures via the "freeze-drying" process is generally well understood, the details of the TTL dehydration mechanisms, including the relative roles of deep convection, atmospheric waves and cloud microphysical processes, remain an active area of research. The dynamical and microphysical processes that influence TTL water vapor concentrations are investigated in simulations of cloud formation and dehydration along air parcel trajectories. We first confirm the validity of our Lagrangian models in a case study involving measurements from the Airborne Tropical TRopopause EXperiment (ATTREX) flights over the central and eastern tropical Pacific in Oct-Nov 2011 and Jan-Feb 2013. ERA-Interim winds and seasonal mean heating rates from Yang et al. (2010) are used to advance parcels back in time from the flight tracks, and time-varying vertical profiles of water vapor along the diabatic trajectories are calculated in a one-dimensional cloud model as in Jensen and Pfister (2004) but with more reliable temperature field, wave and convection schemes. The simulated water vapor profiles demonstrate a significant improvement over estimates based on the Lagrangian Dry Point, agreeing well with aircraft observations when the effects of cloud microphysics, subgrid-scale gravity waves and convection are included. Following this approach, we examine the dynamical and microphysical control of TTL water vapor in the 30ºS-30ºN latitudinal belt and elucidate the dominant processes in the winter and summer seasons. Implications of the TTL dehydration processes for the regulation of global stratospheric humidity will be discussed.

stratosphere↗

Multi-Layer Arctic Mixed-Phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity Experiments

A cloud-resolving model (CRM) is used to simulate the multiple-layer mixed-phase stratiform (MPS) clouds that occurred during a three-and-a-half day subperiod of the Department of Energy-Atmospheric Radiation Measurement Program s Mixed-Phase Arctic Cloud Experiment (M-PACE). The CRM is implemented with an advanced two-moment microphysics scheme, a state-of-the-art radiative transfer scheme, and a complicated third-order turbulence closure. Concurrent meteorological, aerosol, and ice nucleus measurements are used to initialize the CRM. The CRM is prescribed by time-varying large-scale advective tendencies of temperature and moisture and surface turbulent fluxes of sensible and latent heat. The CRM reproduces the occurrences of the single- and double-layer MPS clouds as revealed by the M-PACE observations. However, the simulated first cloud layer is lower and the second cloud layer thicker compared to observations. The magnitude of the simulated liquid water path agrees with that observed, but its temporal variation is more pronounced than that observed. As in an earlier study of single-layer cloud, the CRM also captures the major characteristics in the vertical distributions and temporal variations of liquid water content (LWC), total ice water content (IWC), droplet number concentration and ice crystal number concentration (nis) as suggested by the aircraft observations. However, the simulated mean values differ significantly from the observed. The magnitude of nis is especially underestimated by one order of magnitude. Sensitivity experiments suggest that the lower cloud layer is closely related to the surface fluxes of sensible and latent heat; the upper cloud layer is probably initialized by the large-scale advective cooling/moistening and maintained through the strong longwave (LW) radiative cooling near the cloud top which enhances the dynamical circulation; artificially turning off all ice-phase microphysical processes results in an increase in LWP by a factor of 3 due to interactions between the excessive LW radiative cooling and extra cloud water; heating caused by phase change of hydrometeors could affect the LWC and cloud top height by partially canceling out the LW radiative cooling. It is further shown that the resolved dynamical circulation appears to contribute more greatly to the evolution of the MPS cloud layers than the parameterized subgrid-scale circulation.

Luo, Yali↗

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

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

Dynamic subgrid-scale modeling for high-speed transitional boundary layers

The subgrid scales are modeled dynamically in a large-eddy simulation of transitional boundary-layer flow along a hollow cylinder at a Mach number of 4.5. The behavior of the dynamic-model coefficients, which is determined from local information in the resolved field, is investigated through both an a priori test with direct numerical simulation data for the same case and a complete large-eddy simulation. Both contractions proposed by Germano et al. (1991) and Lilly (1992) are used for the unique determination of the coefficients of the dynamic model, and their results are compared and assessed. The behavior, as well as the energy cascade of the subgridscale field structure, is investigated at various stages of the transition process.

El-Hady, Nabil M.↗

Evaluating and Understanding Parameterized Convective Processes and their Role in the Development of Mesoscale Precipitation Systems

Research efforts during the second year have centered on improving the manner in which convective stabilization is achieved in the Penn State/NCAR mesoscale model MM5. Ways of improving this stabilization have been investigated by (1) refining the partitioning between the Kain-Fritsch convective parameterization scheme and the grid scale by introducing a form of moist convective adjustment; (2) using radar data to define locations of subgrid-scale convection during a dynamic initialization period; and (3) parameterizing deep-convective feedbacks as subgrid-scale sources and sinks of mass. These investigations were conducted by simulating a long-lived convectively-generated mesoscale vortex that occurred during 14-18 Jul. 1982 and the 10-11 Jun. 1985 squall line that occurred over the Kansas-Oklahoma region during the PRE-STORM experiment. The long-lived vortex tracked across the central Plains states and was responsible for multiple convective outbreaks during its lifetime.

Fritsch, J. Michael↗

SST Variation Due to Interactive Convective-Radiative Processes

The recent linking of Cloud-Resolving Models (CRMs) to Ocean-Mixed Layer (OML) models has provided a powerful new means of quantifying the role of cloud systems in ocean-atmosphere coupling. This is due to the fact that the CRM can better resolve clouds and cloud systems and allow for explicit cloud-radiation interaction. For example, Anderson (1997) applied an atmospheric forcing associated with a CRM simulated squall line to a 3-D OML model (one way or passive interaction). His results suggested that the spatial variability resulting from the squall forcing can last at least 24 hours when forced with otherwise spatially uniform fluxes. In addition, the sea surface salinity (SSS) variability continuously decreased following the forcing, while some of the SST variability remained when a diurnal mixed layer capped off the surface structure. The forcing used in the OML model, however, focused on shorter time (8 h) and smaller spatial scales (100-120 km). In this study, the 3-D Goddard Cumulus Ensemble Model (GCE; 512 x 512 x 23 cu km, 2-km horizontal resolution) is used to simulate convective active episodes occurring in the Western Pacific warm pool and Eastern Atlantic regions. The model is integrated for seven days, and the simulated results are coupled to an OML model to better understand the impact of precipitation and changes in the planetary boundary layer upon SST variation. We will specifically examine and compare the results of linking the OML model with various spatially-averaged outputs from GCE simulations (i.e., 2 km vs. 10-50 km horizontal resolutions), in order to help understand the SST sensitivity to multi-scale influences. This will allow us to assess the importance of explicitly simulated deep and shallow clouds, as well as the subgrid-scale effects (in coarse-model runs) upon SST variation. Results using both 1-D and 2-D OML models will be evaluated to assess the effects of horizontal advection.

Tao, W.-K.↗

Large-Eddy Simulation Code Developed for Propulsion Applications

A large-eddy simulation (LES) code was developed at the NASA Glenn Research Center to provide more accurate and detailed computational analyses of propulsion flow fields. The accuracy of current computational fluid dynamics (CFD) methods is limited primarily by their inability to properly account for the turbulent motion present in virtually all propulsion flows. Because the efficiency and performance of a propulsion system are highly dependent on the details of this turbulent motion, it is critical for CFD to accurately model it. The LES code promises to give new CFD simulations an advantage over older methods by directly computing the large turbulent eddies, to correctly predict their effect on a propulsion system. Turbulent motion is a random, unsteady process whose behavior is difficult to predict through computer simulations. Current methods are based on Reynolds-Averaged Navier- Stokes (RANS) analyses that rely on models to represent the effect of turbulence within a flow field. The quality of the results depends on the quality of the model and its applicability to the type of flow field being studied. LES promises to be more accurate because it drastically reduces the amount of modeling necessary. It is the logical step toward improving turbulent flow predictions. In LES, the large-scale dominant turbulent motion is computed directly, leaving only the less significant small turbulent scales to be modeled. As part of the prediction, the LES method generates detailed information on the turbulence itself, providing important information for other applications, such as aeroacoustics. The LES code developed at Glenn for propulsion flow fields is being used to both analyze propulsion system components and test improved LES algorithms (subgrid-scale models, filters, and numerical schemes). The code solves the compressible Favre-filtered Navier- Stokes equations using an explicit fourth-order accurate numerical scheme, it incorporates a compressible form of Smagorinsky s model for the subgrid-scale turbulence, and it uses generalized curvilinear coordinates to allow analysis of a wide range of geometries. The code runs in parallel on shared memory multiprocessor computers and is written in Fortran 90 with dynamic memory allocation. A sample result for a Mach-1.4 round jet is presented in the figure. Instantaneous Mach number contours in several cross-planes downstream of the nozzle exit are shown, illustrating how an LES captures the large unsteady three-dimensional turbulent structures present in the jet.

DeBonis, James R.↗