Engineering Papers⌕ Search

Engineering topics

Nielsen, Eric

Publications and source records attributed to Nielsen, Eric.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Understanding the Impacts of Stellar Companions on Planet Formation and Evolution: A Survey of Stellar and Planetary Companions within 25 pc

We explore the impact of outer stellar companions on the occurrence rate of giant planets detected with radial velocities. We searched for stellar and planetary companions to a volume-limited sample of solar-type stars within 25 pc. Using adaptive optics imaging observations from the Lick 3 m and Palomar 200″ Telescopes, we characterized the multiplicity of our sample stars, down to the bottom of the main sequence. With these data, we confirm field star multiplicity statistics from previous surveys. We additionally combined three decades of radial velocity (RV) data from the California Planet Search with newly collected RV data from Keck/HIRES and the Automated Planet Finder/Levy Spectrometer to search for planetary companions in these same systems. Using an updated catalog of both stellar and planetary companions, as well as detailed injection/recovery tests to determine our sensitivity and completeness, we measured the occurrence rate of planets among the single- and multiple-star systems. We found that planets with masses in the range of 0.1–10 M {sub J} and with semimajor axes of 0.1–10 au have an occurrence rate of 0.18{sub −0.03}{sup +0.04} planets per star when they orbit single stars and an occurrence rate of 0.12 ± 0.04 planets per star when they orbit a star in a binary system. Breaking the sample down by the binary separation, we found that only one planet-hosting binary system had a binary separation <100 au, and none had a separation <50 au. These numbers yielded planet occurrence rates of 0.20{sub −0.06}{sup +0.07} planets per star for binaries with separation a {sub B} > 100 au and 0.04{sub −0.02}{sup +0.04} planets per star for binaries with separation a {sub B} < 100 au. The similarity in the planet occurrence rate around single stars and wide primaries implies that wide binary systems should actually host more planets than single-star systems, since they have more potential host stars. We estimated a system-wide planet occurrence rate of 0.3 planets per wide binary system for binaries with separations a {sub B} > 100 au. Finally, we found evidence that giant planets in binary systems have a different semimajor-axis distribution than their counterparts in single-star systems. The planets in the single-star sample had a significantly higher occurrence rate outside of 1 au than inside 1 au by nearly 4σ, in line with expectations that giant planets are most common near the snow line. However, the planets in the wide binary systems did not follow this distribution, but rather had equivalent occurrence rates interior and exterior to 1 au. This may point to binary-mediated planet migration acting on our sample, even in binaries wider than 100 au.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Obliquity Constraints on an Extrasolar Planetary-mass Companion

We place the first constraints on the obliquity of a planetary-mass companion outside of the solar system. Our target is the directly imaged system 2MASS J01225093–2439505 (2M0122), which consists of a 120 Myr 0.4 M {sub ⊙} star hosting a 12–27 M {sub J} companion at 50 au. We constrain all three of the system’s angular-momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (v sin i) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from Hubble Space Telescope to obtain spin-axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b’s orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin–orbit resonances, and Kozai–Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk—the scenario favored for brown dwarf companions to stars—appears promising.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Individual, Model-independent Masses of the Closest Known Brown Dwarf Binary to the Sun

At a distance of ∼2 pc, our nearest brown dwarf neighbor, Luhman 16 AB, has been extensively studied since its discovery 3 years ago, yet its most fundamental parameter—the masses of the individual dwarfs—has not been constrained with precision. In this work, we present the full astrometric orbit and barycentric motion of Luhman 16 AB and the first precision measurements of the individual component masses. We draw upon archival observations spanning 31 years from the European Southern Observatory (ESO) Schmidt Telescope, the Deep Near-Infrared Survey of the Southern Sky (DENIS), public FORS2 data on the Very Large Telescope (VLT), and new astrometry from the Gemini South Multiconjugate Adaptive Optics System (GeMS). Finally, we include three radial velocity measurements of the two components from VLT/CRIRES, spanning one year. With this new data sampling a full period of the orbit, we use a Markov chain Monte Carlo algorithm to fit a 16-parameter model incorporating mutual orbit and barycentric motion parameters and constrain the individual masses to be 27.9{sub −1.0}{sup +1.1} M{sub J} for the T dwarf and 34.2{sub −1.1}{sup +1.3} M{sub J} for the L dwarf. Our measurements of Luhman 16 AB’s mass ratio and barycentric motion parameters are consistent with previous estimates in the literature utilizing recent astrometry only. The GeMS-derived measurements of the Luhman 16 AB separation in 2014–2015 agree closely with Hubble Space Telescope (HST) measurements made during the same epoch, and the derived mutual orbit agrees with those measurements to within the HST uncertainties of 0.3–0.4 mas.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

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.

Zubair, Mohammad↗
DOE OSTI2016

Bringing "The Moth" to light: A planet-sculpting scenario for the HD 61005 debris disk

Here, the HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2–2.3 μm that further constrains its outer morphology (projected separations of 27–135 au). We also present complementary Gemini Planet Imager 1.6 μm total intensity and polarized light detections that probe down to projected separations less than 10 au. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40–52 au and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 au to a Jupiter mass at 5 au.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

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.

modeling↗
NASA NTRS2011

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.

Pawson, Steven↗
NASA NTRS2010

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.

Strode, Sarah↗
NASA NTRS2009

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.

Ott, Lesley↗
NASA NTRS2009

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.

Pawson, Steven↗