Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atmospheric simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗

An Ocean-Atmosphere Simulation for Studying Air-Sea Interactions

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interactivity between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. To remove this limitation, we aim to couple the km-scale GEOS atmosphere simulation to the km-scale MIT ocean simulation.

Strobach, Ehud↗

Air-Sea Interactions in a High-Resolution Ocean-Atmosphere Simulation

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interaction between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. We present here results from a coupled GEOS-MIT "nature run" simulation, wherein we have coupled a cubed-sphere-720 (~ 1/8) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere simulation to equivalent atmosphere-only and ocean-only simulations. A particular focus of the comparisons is the coupled versus uncoupled differences in interactions between Sea Surface Temperature (SST) and ocean surface wind. We discuss, in particular, a several-day mode of temporal variability in the SST-wind cycle and how it is represented in the different model simulations and in observationally-based products. A mechanism for the cycle, which is driven by SST-wind feedback, is proposed.

Strobach, Ehud↗

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗

Mapping Distant Continental Influences in the Remote Pacific Atmosphere; Simulations of CO Relevant to the Photochemistry of Oxidants

An animated sequence of maps of simulated carbon monoxide concentrations graphically portrays the extent of residual continental influence upon the tropical Pacific Ocean as studied by NASA aircraft during the PEM-Tropics B intensive sampling campaign. We used the MM5 at a 90 km resolution in a globally wrapped grid to simulate the meteorology of transport, and our GRACES model to follow the basic chemistry. The CO we simulate derives from different sources, and so we distinguish anthropogenic, natural terpenoid oxidation, biomass burning, and pervasive CH4-oxidation influences. "Influence" is always judged with an implicit timescale, and these maps describe influence on the 15-45 day timescale appropriate for CO oxidation. In consequence, the maps are useful in assessing the origins of slowly reacting compounds like acetone, methanol, and the lightest hydrocarbons. At 8 km altitude, The Eastern South Pacific to ca. 130 W (eastern Polynesia) was frequently affected by continental influences but NASA's DC-8's flight path did not happen to take it into these regions very often. Near the surface, continentally influenced air crossed into t he Western South Pacific, in the region northwest of the Southern Pacific Convergence Zone but south of the Intertropical Convergence Zone. This air originated from the NE Pacific, and partly from North America. Comparisons are made to CO and other compounds measured aboard the DC-8 and the P-3 aircraft. We will also use tracers to describe the influence of marine convection in the upper troposphere. As time allows, we will discuss the "age" of ozone within the very cleanest region sampled in portions of the near-equatorial Western South Pacific, using a simple chemical mechanism for ozone levels. These simulations describe the chemistry of an atmosphere with very low ozone.

Chatfield Robert B.↗

gdess: A framework for evaluating simulated atmospheric CO 2 in Earth System Models

Atmospheric carbon dioxide (CO 2 ) plays a key role in the global carbon cycle and global warming. Climate-carbon feedbacks are often studied and estimated using Earth System Models (ESMs), which couple together multiple model components—including the atmosphere, ocean, terrestrial biosphere, and cryosphere—to jointly simulate mass and energy exchanges within and between these components. Despite tremendous advances, model intercomparisons and benchmarking are aspects of ESMs that warrant further improvement (Fer et al., 2021; Smith et al., 2014). Such benchmarking is critical because comparing the value of state variables in these simulations against observed values provides evidence for appropriately refining model components; moreover, researchers can learn much about Earth system dynamics in the process (Randall et al., 2019). We introduce `gdess` (a.k.a., Greenhouse gas Diagnostics for Earth System Simulations), which parses observational datasets and ESM simulation output, combines them to be in a consistent structure, computes statistical metrics, and generates diagnostic visualizations. In its current incarnation, `gdess` facilitates evaluating a model's ability to reproduce observed temporal and spatial variations of atmospheric CO 2 . The diagnostics implemented modularly in `gdess` support more rapid assessment and improvement of model-simulated global CO 2 sources and sinks associated with land and ocean ecosystem processes. We intend for this set of automated diagnostics to form an extensible, open source framework for future comparisons of simulated and observed concentrations of various greenhouse gases across Earth system models.

97 MATHEMATICS AND COMPUTING↗

An Ocean-Atmosphere Simulation for Studying Air-Sea Interactions

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interaction between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. To remove this limitation, we aim to perform a coupled simulation using the km-scale GEOS atmosphere and the km-scale MIT ocean models. The initial attempt at the km-scale coupled simulation resulted in computational issues which will be presented here. As a preliminary step towards the km-scale objective, we present results from a high resolution but not yet km-scale simulation, wherein we have coupled a cubed-sphere-720 (~ 1/8) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. A particular focus of the comparisons is the differences in interactions between Sea Surface Temperature (SST) and ocean surface wind for the coupled and uncoupled simulations. We discuss observed and modeled high temporal variability (~days) SST-wind cycle and how it is represented in the different systems. A mechanism for the cycle, which is driven by SST-wind feedback, is proposed.

Strobach, Ehud↗

Generation of realistic input parameters for simulating atmospheric point-spread functions at astronomical observatories

High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear – correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all other sources of correlated image distortions must be modeled or removed. One such source is the correlated blurring of images due to optical turbulence in the atmosphere, which dominates the point-spread function (PSF) for ground-based instruments. In this work, we leverage data from weather forecasting models to produce wind speeds and directions, and turbulence parameters, that are realistically correlated with altitude. To study the resulting correlations in the size and shape of the PSF, we generate simulated images of the PSF across a ~10 square-degree field of view – the size of the camera focal plane for the Vera C. Rubin Observatory in Chile – using weather data and historical seeing for a geographic location near the Observatory. We make quantitative predictions for two-point correlation functions (2PCF) that are used in analyses of cosmic shear. We observe a strong anisotropy in the two-dimensional 2PCF, which is expected based on observations in real images, and study the dependence of the orientation of the anisotropy on dominant wind directions near the ground and at higher altitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of interactive and prescribed agricultural ammonia emissions for simulating atmospheric composition in CAM-chem

Abstract. Ammonia (NH3) plays a central role in the chemistry of inorganic secondary aerosols in the atmosphere. The largest emission sector for NH3 is agriculture, where NH3 is volatilized from livestock wastes and fertilized soils. Although the NH3 volatilization from soils is driven by the soil temperature and moisture, many atmospheric chemistry models prescribe the emission using yearly emission inventories and climatological seasonal variations. Here we evaluate an alternative approach where the NH3 emissions from agriculture are simulated interactively using the process model FANv2 (Flow of Agricultural Nitrogen, version 2) coupled to the Community Atmospheric Model with Chemistry (CAM-chem). We run a set of 6-year global simulations using the NH3 emission from FANv2 and three global emission inventories (EDGAR, CEDS and HTAP) and evaluate the model performance using a global set of multi-component (atmospheric NH3 and NH4+, and NH4+ wet deposition) in situ observations. Over East Asia, Europe and North America, the simulations with different emissions perform similarly when compared with the observed geographical patterns. The seasonal distributions of NH3 emissions differ between the inventories, and the comparison to observations suggests that both FANv2 and the inventories would benefit from more realistic timing of fertilizer applications. The largest differences between the simulations occur over data-scarce regions. In Africa, the emissions simulated by FANv2 are 200 %–300 % higher than in the inventories, and the available in situ observations from western and central Africa, as well as NH3 retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) instrument, are consistent with the higher NH3 emissions as simulated by FANv2. Overall, in simulating ammonia and ammonium concentrations over regions with detailed regional emission inventories, the inventories based on these details (HTAP, CEDS) capture the atmospheric concentrations and their seasonal variability the best. However these inventories cannot capture the impact of meteorological variability on the emissions, nor can these inventories couple the emissions to the biogeochemical cycles and their changes with climate drivers. Finally, we show with sensitivity experiments that the simulated time-averaged nitrate concentration in air is sensitive to the temporal resolution of the NH3 emissions. Over the CASTNET monitoring network covering the US, resolving the NH3 emissions hourly instead monthly reduced the positive model bias from approximately 80 % to 60 % of the observed yearly mean nitrate concentration. This suggests that some of the commonly reported overestimation of aerosol nitrate over the US may be related to unresolved temporal variability in the NH3 emissions.

54 ENVIRONMENTAL SCIENCES↗

Atmospheric simulator and calibration system for remote sensing radiometers

A system for calibrating the MAPS (measurement of air pollution from satellites) instruments was developed. The design of the system provides a capability for simulating a broad range of radiant energy source temperatures and a broad range of atmospheric pressures, temperatures, and pollutant concentrations for a single slab atmosphere. The system design and the system operation are described.

Holland, J. A.↗

Simulating Atmospheric Turbulence

Dimensionless time series for gusts and gradients calculated with model. Single-sided one-dimensional spectra obtained from Fourier analysis of simulated vertical gust. Each of four plots represents spectrum in one of four altitude bands ranging from 0 to 10 km. Model differs from most models; is nondimensional for greater generality and uses spectra with finite wave number limits. Developed to predict effects of atmospheric turbulence in horizontal and near-horizontal flight also applicable to commercial aviation.

Tatom, F. B.↗

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗

Impact of surface coupling grids on tropical cyclone extremes in high-resolution atmospheric simulations

This article discusses the sensitivity of tropical cyclone climatology to surface coupling strategy in high-resolution configurations of the Community Earth System Model. Using two supported model setups, we demonstrate that the choice of grid on which the lowest model level wind stress and surface fluxes are computed may lead to differences in cyclone strength in multi-decadal climate simulations, particularly for the most intense cyclones. Using a deterministic framework, we show that when these surface quantities are calculated on an ocean grid that is coarser than the atmosphere, the computed frictional stress is misaligned with wind vectors in individual atmospheric grid cells. This reduces the effective surface drag, and results in more intense cyclones when compared to a model configuration where the ocean and atmosphere are of equivalent resolution. Our results demonstrate that the choice of computation grid for atmosphere–ocean interactions is non-negligible when considering climate extremes at high horizontal resolution, especially when model components are on highly disparate grids.

54 ENVIRONMENTAL SCIENCES↗

Long-Term Simulated Atmospheric Nitrogen Deposition Alters Leaf and Fine Root Decomposition

Atmospheric nitrogen deposition increases forest carbon sequestration across broad parts of the Northern Hemisphere. Slower organic matter decomposition and greater soil carbon accumulation could contribute to this increase in carbon sequestration. We investigated the effects of chronic simulated nitrogen deposition on leaf litter and fine root decomposition at four sugar maple (Acer saccharum)-dominated northern hardwood forests. At these sites, we previously observed that nitrogen additions increased soil organic carbon and altered litter chemistry. We conducted a 3-year decomposition study with litter bags. Litter production of leaves and fine roots were combined with decomposition dynamics to estimate how fine roots and leaf litter contribute to soil organic carbon. We found that nitrogen additions marginally stimulated early-stage decomposition of leaf litter, an effect associated with previously documented changes in litter chemistry. In contrast, nitrogen additions inhibited the later stages of fine root decomposition, which is consistent with observed decreases in lignin-degrading enzyme activities with nitrogen additions at these sites. At the ecosystem scale, slower fine root decomposition led to additional root mass retention (g m -2 ), and this greater retention of root residues was estimated to explain 5–51% of previously documented carbon accumulation in the surface soil due to nitrogen additions. Our results demonstrated that simulated nitrogen deposition created contrasting effects on the decomposition of leaf litter and fine roots. Although previous nitrogen deposition studies have focused on leaf litter, this work suggests that slower fine root decomposition is a major driver of soil organic carbon accumulation under elevated nitrogen deposition.

54 ENVIRONMENTAL SCIENCES↗

Verification of interim simulated atmospheric turbulence

During horizontal and near-horizontal flight of the space shuttle orbiter, the effects of atmospheric turbulence on the vehicle are crucial in establishing the space shuttle design, control, and human pilot effects. The effects of both turbulent gusts and gust gradients (shears) must be taken into account. By means of a nonrecursive turbulence simulation model, twelve simulated turbulence time series for gusts and gust gradients were generated as part of a previous investigation. The current investigation was concerned with the spectral analysis of the time series based on the use of Fourier transform theory. A comparison between the computed spectra and the corresponding von Karman theoretical spectra is also included. A basic difference is explained between the computed and theoretical spectra which was identified as a result of the comparison. Conclusions and recommendations are presented.

Tatom, F. B.↗

Simulating Atmospheric Processes in Earth System Models and Quantifying Uncertainties With Deep Learning Multi‐Member and Stochastic Parameterizations

Abstract Deep learning is a powerful tool to represent subgrid processes in climate models, but many application cases have so far used idealized settings and deterministic approaches. Here, we develop stochastic parameterizations with calibrated uncertainty quantification to learn subgrid convective and turbulent processes and surface radiative fluxes of a superparameterization embedded in an Earth System Model (ESM). We explore three methods to construct stochastic parameterizations: (a) a single Deep Neural Network (DNN) with Monte Carlo Dropout; (b) a multi‐member parameterization; and (c) a Variational Encoder Decoder with latent space perturbation. We show that the multi‐member parameterization improves the representation of convective processes, especially in the planetary boundary layer, compared to individual DNNs. The respective uncertainty quantification illustrates that methods (b) and (c) are advantageous compared to a dropout‐based DNN parameterization regarding the spread of convective processes. Hybrid simulations with our best‐performing multi‐member parameterizations remained challenging and crash within the first days. Therefore, we develop a pragmatic partial coupling strategy relying on the superparameterization for condensate emulation. Partial coupling reduces the computational efficiency of hybrid Earth‐like simulations but enables model stability over 5 months with our multi‐member parameterizations. However, our hybrid simulations exhibit biases in thermodynamic fields and differences in precipitation patterns. Despite this, the multi‐member parameterizations enable improvements in reproducing tropical extreme precipitation compared to a traditional convection parameterization. Despite these challenges, our results indicate the potential of a new generation of multi‐member machine learning parameterizations leveraging uncertainty quantification to improve the representation of stochasticity of subgrid effects.

Behrens, Gunnar [Deutsches Zentrum für Luft‐ und R↗

Expanding the computational frontier of multi-scale atmospheric simulation to advance understanding of low cloud / climate feedbacks: Final Technical Report

Marine boundary layer clouds are a challenge for climate models. They cover much of the oceans, but they are driven by small-scale turbulent eddies only a few hundred meters across, which in turn respond to the cloud formation. Global climate models have a grid spacing that is far too large to simulate such fluid motions, so they parameterize these cloud formation processes. Different climate modeling groups have designed different parameterizations in which the clouds turn out to be differently sensitive to a warming climate. This drives uncertainties in our best guess at the sensitivity of global warming to greenhouse gas increases. If these cloud-forming eddies could be directly simulated using the well-known equations of fluid motion, they would no longer need to be parameterized, removing a major source of climate modeling uncertainty. In this project, we overcame software engineering challenges to successfully implemented ‘ultraparameterization’ (UP), the first global model that does this, and we tested how well it works. We simulated five-year periods with present-day temperatures and with a warmer climate, and we investigated how the UP-simulated clouds responded to climate – the ‘cloud feedback’ problem. We found little response of clouds at all latitudes to the imposed climate change, which is within the range of predictions of conventional global climate models. UP is a variation on superparameterization, in which small cloud-resolving models (CRMs) are embedded in each column of the global model. In UP, the CRM grid is fine enough (250 m horizontal × 20 m vertical) to explicitly capture boundary-layer turbulent eddies and associated clouds. Because only one small columnar patch is simulated within each climate model grid cell, this is a million-fold more efficient than simulating the entire globe on this same CRM grid., but achieves much of the same effect for the cloud properties. It doesn’t work perfectly. For instance, like conventional climate models, UP simulates too little subtropical stratocumulus cloud, a bias that we are continuing to work to reduce. However, because it directly simulates the turbulent cloud-forming processes, UP is inherently more plausible for simulating how clouds will change in a perturbed climate. In future, we hope to apply UP to another key climate modeling issue: cloud-aerosol interaction and the effect of human-produced aerosols on the climate change we have already experienced and that which is likely to come.

54 ENVIRONMENTAL SCIENCES↗