Search NASA⌕ Search

Engineering topics

Bell, J.

Publications and source records attributed to Bell, J..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Porting WarpX to GPU-accelerated platforms

WarpX is a general purpose electromagnetic particle-in-cell code that was originally designed to run on many-core CPU architectures. We describe the strategy, based on the AMReX library, followed to allow WarpX to use the GPU-accelerated nodes on OLCF’s Summit supercomputer, a strategy we believe will extend to the upcoming machines Frontier and Aurora. We summarize the challenges encountered, lessons learned, and give current performance results on a series of relevant benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling of a chain of three plasma accelerator stages with the WarpX electromagnetic PIC code on GPUs

The fully electromagnetic particle-in-cell code WarpX is being developed by a team of the U.S. DOE Exascale Computing Project (with additional non-U.S. collaborators on part of the code) to enable the modeling of chains of tens to hundreds of plasma accelerator stages on exascale supercomputers, for future collider designs. The code is combining the latest algorithmic advances (e.g., Lorentz boosted frame and pseudo-spectral Maxwell solvers) with mesh refinement and runs on the latest computer processing unit and graphical processing unit (GPU) architectures. In this paper, we summarize the strategy that was adopted to port WarpX to GPUs, report on the weak parallel scaling of the pseudo-spectral electromagnetic solver, and then present solutions for decreasing the time spent in data exchanges from guard regions between subdomains. In Sec. IV, we demonstrate the simulations of a chain of three consecutive multi-GeV laser-driven plasma accelerator stages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Warp-X: A new exascale computing platform for beam–plasma simulations

Turning the current experimental plasma accelerator state-of-the-art from a promising technology into mainstream scientific tools depends critically on high-performance, high-fidelity modeling of complex processes that develop over a wide range of space and time scales. As part of the U.S. Department of Energy's Exascale Computing Project, a team from Lawrence Berkeley National Laboratory, in collaboration with teams from SLAC National Accelerator Laboratory and Lawrence Livermore National Laboratory, is developing a new plasma accelerator simulation tool that will harness the power of future exascale supercomputers for high-performance modeling of plasma accelerators. We present the various components of the codes such as the new Particle-In-Cell Scalable Application Resource (PICSAR) and the redesigned adaptive mesh refinement library AMReX, which are combined with redesigned elements of the Warp code, in the new WarpX software. Lastly, the code structure, status, early examples of applications and plans are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

MAESTRO, CASTRO, and SEDONA -- Petascale Codes for Astrophysical Applications

Performing high-resolution, high-fidelity, three-dimensional simulations of Type Ia supernovae (SNe Ia) requires not only algorithms that accurately represent the correct physics, but also codes that effectively harness the resources of the most powerful supercomputers. We are developing a suite of codes that provide the capability to perform end-to-end simulations of SNe Ia, from the early convective phase leading up to ignition to the explosion phase in which deflagration/detonation waves explode the star to the computation of the light curves resulting from the explosion. In this paper we discuss these codes with an emphasis on the techniques needed to scale them to petascale architectures. We also demonstrate our ability to map data from a low Mach number formulation to a compressible solver.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Theoretical UV absorption spectra of hydrodynamically escaping O{sub 2}/CO{sub 2}-rich exoplanetary atmospheres

Characterizing Earth- and Venus-like exoplanets' atmospheres to determine if they are habitable and how they are evolving (e.g., equilibrium or strong erosion) is a challenge. For that endeavor, a key element is the retrieval of the exospheric temperature, which is a marker of some of the processes occurring in the lower layers and controls a large part of the atmospheric escape. We describe a method to determine the exospheric temperature of an O{sub 2}- and/or CO{sub 2}-rich transiting exoplanet, and we simulate the respective spectra of such a planet in hydrostatic equilibrium and hydrodynamic escape. The observation of hydrodynamically escaping atmospheres in young planets may help constrain and improve our understanding of the evolution of the solar system's terrestrial planets' atmospheres. We use the dependency of the absorption spectra of the O{sub 2} and CO{sub 2} molecules on the temperature to estimate the temperature independently of the total absorption of the planet. Combining two observables (two parts of the UV spectra that have a different temperature dependency) with the model, we are able to determine the thermospheric density profile and temperature. If the slope of the density profile is inconsistent with the temperature, then we infer the hydrodynamic escape. We address the question of the possible biases in the application of the method to future observations, and we show that the flare activity should be cautiously monitored to avoid large biases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

CARBON DEFLAGRATION IN TYPE Ia SUPERNOVA. I. CENTRALLY IGNITED MODELS

A leading model for Type Ia supernovae (SNe Ia) begins with a white dwarf near the Chandrasekhar mass that ignites a degenerate thermonuclear runaway close to its center and explodes. In a series of papers, we shall explore the consequences of ignition at several locations within such dwarfs. Here we assume central ignition, which has been explored before, but is worth revisiting, if only to validate those previous studies and to further elucidate the relevant physics for future work. A perturbed sphere of hot iron ash with a radius of {approx}100 km is initialized at the middle of the star. The subsequent explosion is followed in several simulations using a thickened flame model in which the flame speed is either fixed-within the range expected from turbulent combustion-or based on the local turbulent intensity. Global results, including the explosion energy and bulk nucleosynthesis (e.g., {sup 56}Ni of 0.48-0.56 M{sub Sun }) turn out to be insensitive to this speed. In all completed runs, the energy released by the nuclear burning is adequate to unbind the star, but not enough to give the energy and brightness of typical SNe Ia. As found previously, the chemical stratification observed in typical events is not reproduced. These models produce a large amount of unburned carbon and oxygen in central low velocity regions, which is inconsistent with spectroscopic observations, and the intermediate mass elements and iron group elements are strongly mixed during the explosion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

CASTRO: A NEW COMPRESSIBLE ASTROPHYSICAL SOLVER. III. MULTIGROUP RADIATION HYDRODYNAMICS

We present a formulation for multigroup radiation hydrodynamics that is correct to order O(v/c) using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses a Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts: one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A Multifaceted Mathematical Approach for Complex Systems

Applied mathematics has an important role to play in developing the tools needed for the analysis, simulation, and optimization of complex problems. These efforts require the development of the mathematical foundations for scientific discovery, engineering design, and risk analysis based on a sound integrated approach for the understanding of complex systems. However, maximizing the impact of applied mathematics on these challenges requires a novel perspective on approaching the mathematical enterprise. Previous reports that have surveyed the DOE's research needs in applied mathematics have played a key role in defining research directions with the community. Although these reports have had significant impact, accurately assessing current research needs requires an evaluation of today's challenges against the backdrop of recent advances in applied mathematics and computing. To address these needs, the DOE Applied Mathematics Program sponsored a Workshop for Mathematics for the Analysis, Simulation and Optimization of Complex Systems on September 13-14, 2011. The workshop had approximately 50 participants from both the national labs and academia. The goal of the workshop was to identify new research areas in applied mathematics that will complement and enhance the existing DOE ASCR Applied Mathematics Program efforts that are needed to address problems associated with complex systems. This report describes recommendations from the workshop and subsequent analysis of the workshop findings by the organizing committee.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

ROCK AND SOIL PHYSICAL PROPERTIES AT THE MER GUSEV CRATER AND MERIDIANI PLANUM LANDING SITES

Following the successful landings of both Mars Exploration Rover (MER) vehicles at Gusev Crater and Meridiani Planum, respectively, their Athena suite of instruments is being used to study the geologic history of these two very different landing sites on Mars that had been selected on the basis of showing different types of evidence for aqueous processes in the planet s past. Utilizing the on-board instruments as well as the rovers mobility system, a wide range of physical properties investigations is carried out as well - the subject of this abstract - that provide additional information on the geology and processes at the sites. Results of the mission in general as well as of the physical properties studies thus far greatly exceed expectations in that observations and measurements by both vehicles show a rich variety in materials and processes: the Gusev site in the vicinity of the lander is remarkably flat and generally devoid of large rocks along traverses up to the time of this writing (approx.Sol 50) and suggestive of a deflated surface with generally only thin veneers of bright dust while exhibiting evidence of a widespread occurrence of a crust from cemented fines that has been observed to fail in the form of blocky clods when disturbed by vehicle rolling action; numerous small and shallow depressions - presumably created by impacts - are observed at the site which are infilled with bright, fine-grained material that likewise appears indurated and which was studied by a trenching experiment; small ripple bedforms are scattered across the site and were characterized in terms of particle size distributions. At the Meridiani site, studies so far - up to approx.Sol 33 - have focussed on soils and the rock outcrop encountered within the approx.20 m diameter crater that the spacecraft came to rest in: from a physical properties point of view, a mantle of dark, well-sorted, apparently basaltic sand with small to moderate cohesion has been of interest - and has been studied by a trenching experiment - as well as a fine-grained unit underlying the mantle at least locally within the crater. Rock grindings were accomplished successfully at both sites at the time of this writing, suggesting different strengths of the two targets (the basaltic rock nicknamed Adirondack at Gusev and the Meridiani rock outcrop) in addition to enabling compositional measurements below the original rock surfaces.

Ming, Douglas W.↗