Engineering topics
Nielsen, Eric
Publications and source records attributed to Nielsen, Eric.
An Optimized Multicolor Point-Implicit Solver for Unstructured Grid Applications on Graphics Processing Units
In the field of computational fluid dynamics, the Navier-Stokes equations are often solved using an unstructuredgrid approach to accommodate geometric complexity. Implicit solution methodologies for such spatial discretizations generally require frequent solution of large tightly-coupled systems of block-sparse linear equations. The multicolor point-implicit solver used in the current work typically requires a significant fraction of the overall application run time. In this work, an efficient implementation of the solver for graphics processing units is proposed. Several factors present unique challenges to achieving an efficient implementation in this environment. These include the variable amount of parallelism available in different kernel calls, indirect memory access patterns, low arithmetic intensity, and the requirement to support variable block sizes. In this work, the solver is reformulated to use standard sparse and dense Basic Linear Algebra Subprograms (BLAS) functions. However, numerical experiments show that the performance of the BLAS functions available in existing CUDA libraries is suboptimal for matrices representative of those encountered in actual simulations. Instead, optimized versions of these functions are developed. Depending on block size, the new implementations show performance gains of up to 7x over the existing CUDA library functions.
GMI Capabilities
We describe the capabilities of the Global Modeling Initiative (GMI) chemical transport model (CTM) with a special focus on capabilities related to the Atmospheric Tomography Mission (ATom). Several science results based on GMI hindcast simulations and preliminary results from the ATom simulations are highlighted. We also discuss the relationship between GMI and GEOS-5.
FUN3D v12.7 Training: Session 9: Adjoint-Based Design for Steady Flows
No abstract available
FUN3D v12.7 Training: Session 4: Gridding, Solution and Visualization Basics
No abstract available
Atmospheric Constituents in GEOS-5: Components for an Earth System Model
The GEOS-S model is being developed for weather and climate processes, including the implementation of "Earth System" components. While the stratospheric chemistry capabilities are mature, we are presently extending this to include predictions of the tropospheric composition and chemistry - this includes CO2, CH4, CO, nitrogen species, etc. (Aerosols are also implemented, but are beyond the scope of this paper.) This work will give an overview of our chemistry modules, the approaches taken to represent surface emissions and uptake of chemical species, and some studies of the sensitivity of the atmospheric circulation to changes in atmospheric composition. Results are obtained through focused experiments and multi-decadal simulations.
Transport of Cs-137 from Boreal Biomass Burning in Summer of 2010
The summer of 2010 was a severe fire season in western Russia. Wildfires were detected in the Bryansk region, raising concerns that radionuclide contamination from the Chernobyl accident could be resuspended in the atmosphere. We simulate the transport of passive and particulate tracers of biomass burning from this region using the GEOS5 GOCART model driven by assimilated meteorology. Biomass burning emissions are based on MODIS fire detections. We validate the model against aerosol optical depth from MODIS. Using a range of estimates for Cs-137 emissions during wildfires, we estimate the downwind concentration and deposition of Cs-137 based on the emission ratios of Cs-137 to the simulated tracers. We discuss the sensitivity of our results to the location of the fires and the fraction of Cs-137 resuspended.
The Influence of the 2006 Indonesian Biomass Burning Aerosols on Tropical Dynamics Studied with the GEOS-5 AGCM
The direct and semi-direct effects of aerosols produced by Indonesian biomass burning (BB) during August November 2006 on tropical dynamics have been examined using NASA's Goddard Earth Observing System, Version 5 (GEOS-5) atmospheric general circulation model (AGCM). The AGCM includes CO, which is transported by resolved and sub-grid processes and subject to a linearized chemical loss rate. Simulations were driven by two sets of aerosol forcing fields calculated offline, one that included Indonesian BB aerosol emissions and one that did not. In order to separate the influence of the aerosols from internal model variability, the means of two ten-member ensembles were compared. Diabatic heating from BB aerosols increased temperatures over Indonesia between 150 and 400 hPa. The higher temperatures resulted in strong increases in upward grid-scale vertical motion, which increased water vapor and CO over Indonesia. In October, the largest increases in water vapor were found in the mid-troposphere (~25%) while the largest increases in CO occurred just below the tropopause (80 ppbv or ~50%). Diabatic heating from the Indonesian BB aerosols caused CO to increase by 9% throughout the tropical tropopause layer in November and 5% in the lower stratosphere in December. The results demonstrate that aerosol heating plays an important role in the transport of BB pollution and troposphere-to-stratosphere transport. Changes in vertical motion and cloudiness induced by aerosol heating can also alter the transport and phase of water vapor in the upper troposphere/lower stratosphere.
The GEOS-Chemistry-Climate Model (CCM), Versions 1,2 and 3
This paper examines evolution of key ozone and climate aspects in our range of Goddard Earth Observing System Models, Versions 1-3. Version 1 (Pawson et al., 2008; )GR) used GSFC stratospheric ozone chemistry in GEOS-4 GCM. A chronic high bias in polar ozone at low chlorine loading persisted into Version 2, which used the GEQS-5 GCM and the same chemistry. This is much improved in Version 3, which uses the GMI-COMBO stratosphere-troposphere chemistry. A dynamical problem with Version 1, the overactive nature of the Antarctic polar vortex, is corrected in Versions 2 and 3, when GEQS-5 replaces GEO5-4. Other changes between the various model versions will be documented.
Aerothermodynamic Analyses of Towed Ballutes
A ballute (balloon-parachute) is an inflatable, aerodynamic drag device for application to planetary entry vehicles. Two challenging aspects of aerothermal simulation of towed ballutes are considered. The first challenge, simulation of a complete system including inflatable tethers and a trailing toroidal ballute, is addressed using the unstructured-grid, Navier-Stokes solver FUN3D. Auxiliary simulations of a semi-infinite cylinder using the rarefied flow, Direct Simulation Monte Carlo solver, DSV2, provide additional insight into limiting behavior of the aerothermal environment around tethers directly exposed to the free stream. Simulations reveal pressures higher than stagnation and corresponding large heating rates on the tether as it emerges from the spacecraft base flow and passes through the spacecraft bow shock. The footprint of the tether shock on the toroidal ballute is also subject to heating amplification. Design options to accommodate or reduce these environments are discussed. The second challenge addresses time-accurate simulation to detect the onset of unsteady flow interactions as a function of geometry and Reynolds number. Video of unsteady interactions measured in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel and CFD simulations using the structured grid, Navier-Stokes solver LAURA are compared for flow over a rigid spacecraft-sting-toroid system. The experimental data provides qualitative information on the amplitude and onset of unsteady motion which is captured in the numerical simulations. The presence of severe unsteady fluid - structure interactions is undesirable and numerical simulation must be able to predict the onset of such motion.
BOREAS TE-22 Allometric Forest Survey Data
The Boreal Ecosystem-Atmospheric Study (BOREAS) TE-22 (Terrestrial Ecology) team collected data sets in support of its efforts to characterize and interpret information on the forest structure of boreal vegetation in the Southern and Northern Study Areas (SSA and NSA) during the 1994 growing season. The data are stored in tabular ASCII files. The data files are available on a CD-ROM (see document number 20010000884), or from the Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC).