SEARCH · Search NASA
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.
The Allan variance for an averaging time tau can be estimated either from all available phase samples or from a subgrid of samples with spacing tau. This paper gives a set of computational recipes that yield the variance of both estimators, with less than 2 percent error, for the five power-law components of the classical continuous-time clock noise model.
In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.
One of the most exciting recent developments in the field of large eddy simulation (LES) is the dynamic subgrid-scale model. The dynamic model concept is a general procedure for evaluating model constants by sampling a band of the smallest scales actually resolved in the simulation. To date, the procedure has been used primarily in conjunction with the Smagorinsky model. The dynamic procedure has the advantage that the value of the model constant need not be specified a priori, but rather is calculated as a function of space and time as the simulation progresses. This feature makes the dynamic model especially attractive for flows in complex geometries where it is difficult or impossible to calibrate model constants. The dynamic model was highly successful in benchmark tests involving homogeneous and channel flows. Having demonstrated the potential of the dynamic model in these simple flows, the overall direction of the LES effort at CTR shifted toward an evaluation of the model in more complex situations. The current test cases are basic engineering-type flows for which Reynolds averaged approaches were unable to model the turbulence to within engineering accuracy. Flows currently under investigation include a backward-facing step, wake behind a circular cylinder, airfoil at high angles of attack, separated flow in a diffuser, and boundary layer over a concave surface. Preliminary results from the backward-facing step and cylinder wake simulations are encouraging. Progress on the LES of a boundary layer on a concave surface is discussed. Although the geometry of a concave wall is not very complex, the boundary layer that develops on its surface is difficult to model due to the presence of streamwise Taylor-Gortler vortices. These vortices arise as a result of a centrifugal instability associated with the convex curvature.
Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. To better understand the magnitude and possible causes of this problem, we analyze the disagreements between the Arctic System Reanalysis version 2 (ASR) and data taken during the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign conducted over the Beaufort Sea. Previous analysis of grid scale-averaged meteorological variables revealed a lack of cloud water produced in ASR, related to warm/dry biases, a failure to produce sufficient cloud water in conditions that are observed to be favorable for large cloud water values, and sampling issues with the ARISE flight paths. However, ARISE collected meteorological observations at one second intervals, offering spatial data resolution comparable to the subgrid scale of ASR. Subgrid scale properties and processes must be parameterized in conventional atmospheric models and reanalyses such as ASR, and often rely on assumed distributions of meteorological variables about the grid box mean state. To enable an examination of the realism of ASR’s subgrid scale properties, we use high resolution data from ARISE determine the “subgrid scale” distributions of cloud water and thermodynamic variables in the real atmosphere. Preliminary results show that the subgrid scale distribution of thermodynamic variables about the grid scale mean resembles the triangular distribution assumed in many bulk microphysical parameterization schemes, but also with tails resembling Gaussian distributions. The main goal is to understand the partitioning of the total available atmospheric water into water vapor and cloud water. Models typically covert some percentage of total water in a grid box to cloud water based on the estimated amount of supersaturation from the assumed vapor distribution. The ARISE measurements can reveal the relationship between subgrid scale supersaturation and the partition between cloud water and water vapor.
Predicting realistic cloud subgrid variability remains a challenge for Global Climate Models (GCMs) even though it can play an important role for proper representation of processes pertaining to cloud microphysics, precipitation, and radiation, but also for comparisons with satellite observations which are of much higher resolution than model grids. The diagnostic approach of subgrid cloud variability is commonly handled by subcolumn cloud generators. For a specific GCM, one ideally wants the subgrid variability used for comparisons with satellite observations to be created by the same generator and with the same rules as the one used for model integration. With this in mind, we have embarked in an effort to test and improve cloud subcolumn generators appropriate for GCMs. For this purpose, we use cloud (hydrometeor) products from active observations by the CloudSat radar (CPR) and the CALIPSO lidar (CALIOP). Cloud products from active sensors while suffering significant sampling and coverage drawbacks have the advantage of resolving both horizontal and vertical variability. The main question is: given a profile of cloud condensate mean and variance, can we create a subgrid cloud field that is statistically similar to the observed subgrid cloud field? By “statistically”, we suggest that we do not aspire to reproduce “well” each individual subgrid cloud field of a GCM-scale region, but that our generator performs well for a large ensemble of cases. We simulate radar, passive imager, and radiation flux fields from the observed 2D cloud fields and create one-point statistics; we use the profiles of cloud condensate mean and variance as input to the generator to create subgrid cloud fields; we compare the statistics; we adjust the rules of the generator to create the best possible agreement between the radar, imager and radiation field statistics. In this process, the active observations have actually a dual role: they provide the actual subgrid cloud field which can be reduced to a profile of mean and variance used by the subcolumn generator, but they also provide the rules needed by the generator, such as measures (e.g, decorrelation length) of the vertical overlap of cloud fraction and of the condensate horizontal variability. This presentation will show our progress using Cloudsat and CALIPSO products in this dual fashion.
Predicting realistic cloud subgrid variability remains a challenge for Global Climate Models (GCMs) even though it can play an important role for proper representation of processes pertaining to cloud microphysics, precipitation, and radiation, but also for comparisons with satellite observations which are of much higher resolution than model grids. The diagnostic approach of subgrid cloud variability is commonly handled by subcolumn cloud generators. For a specific GCM, one ideally wants the subgrid variability used for comparisons with satellite observations to be created by the same generator and with the same rules as the one used for model integration. With this in mind, we have embarked in an effort to test and improve cloud subcolumn generators appropriate for GCMs. For this purpose, we use cloud (hydrometeor) products from active observations by the CloudSat radar (CPR) and the CALIPSO lidar (CALIOP). Cloud products from active sensors while suffering signicant sampling and coverage drawbacks have the advantage of resolving both horizontal and vertical variability. The main question is: given a profile of cloud condensate mean and variance, can we create a subgrid cloud field that is statistically similar to the observed subgrid cloud field? By “statistically”, we suggest that we do not aspire to reproduce “well” each individual subgrid cloud field of a GCM-scale region, but that our generator performs well for a large ensemble of cases. We simulate radar, passive imager, and radiation flux fields from the observed 2D cloud fields and create one-point statistics; we use the profiles of cloud condensate mean and variance as input to the generator to create subgrid cloud fields; we compare the statistics; we adjust the rules of the generator to create the best possible agreement between the radar, imager and radiation field statistics. In this process, the active observations have actually a dual role: they provide the actual subgrid cloud field which can be reduced to a profile of mean and variance used by the subcolumn generator, but they also provide the rules needed by the generator, such as measures (e.g, decorrelation length) of the vertical overlap of cloud fraction and of the condensate horizontal variability. This presentation will show our progress using Cloudsat and CALIPSO products in this dual fashion.
Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales
Abstract Machine‐learning (ML) parameterizations of subgrid processes (here of turbulence, convection, and radiation) may one day replace conventional parameterizations by emulating high‐resolution physics without the cost of explicit simulation. However, uncertainty about the relationship between offline and online performance (i.e., when integrated with a large‐scale general circulation model) hinders their development. Much of this uncertainty stems from limited sampling of the noisy, emergent effects of upstream ML design decisions on downstream online hybrid simulation. Our work rectifies the sampling issue via the construction of a semi‐automated, end‐to‐end pipeline for size ensembles of hybrid simulations, revealing important nuances in how systematic reductions in offline error manifest in changes to online error and online stability. For example, removing dropout and switching from a Mean Squared Error to a Mean Absolute Error loss both reduce offline error, but they have opposite effects on online error and online stability. Other design decisions, like incorporating memory, converting moisture input from specific humidity to relative humidity, using batch normalization, and training on multiple climates do not come with any such compromises. Finally, we show that ensemble sizes of may be necessary to reliably detect causally relevant differences online. By enabling rapid online experimentation at scale, we can empirically settle debates regarding subgrid ML parameterization design that would have otherwise remained unresolved in the noise.
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.
One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.
One fundamental question about atmospheric moist convection processes that remains debated is whether or under what conditions a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using two different reference frames: static cloudy-updraft grid cells versus tracked cumulus thermals. In both frameworks we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy-updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy-updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy-updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.
The NASA Goddard Earth Observing System (GEOS-5) Chemistry and Climate Model (CCM) consists of a global atmospheric general circulation model and the combined stratospheric and tropospheric chemistry package from the NASA Global Modeling Initiative (GMI) chemical transport model. The subgrid process of convective tracer transport is represented through the Relaxed Arakawa-Schubert parameterization in the GEOS-5 CCM. However, substantial uncertainty for tracer transport is associated with this parameterization, as is the case with all global and regional models. We have designed a project to comprehensively evaluate this parameterization from the point of view of tracer transport, and determine the most appropriate improvements that can be made to the GEOS-5 convection algorithm, allowing improvement in our understanding of the role of convective processes in determining atmospheric composition. We first simulate tracer transport in individual observed convective events with a cloud-resolving model (WRF). Initial condition tracer profiles (CO, CO2, O3) are constructed from aircraft data collected in undisturbed air, and the simulations are evaluated using aircraft data taken in the convective anvils. A single-column (SCM) version of the GEOS-5 GCM with online tracers is then run for the same convective events. SCM output is evaluated based on averaged tracer fields from the cloud-resolving model. Sensitivity simulations with adjusted parameters will be run in the SCM to determine improvements in the representation of convective transport. The focus of the work to date is on tropical continental convective events from the African Monsoon Multidisciplinary Analyses (AMMA) field mission in August 2006 that were extensively sampled by multiple research aircraft.