Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial accelerator”

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 685 records · Page 38

Towards performance portability in the Spark astrophysical magnetohydrodynamics solver in the Flash-X simulation framework

Simulations of core-collapse supernovae, and other astrophysical phenomena, are quintessential extreme-scale computing challenges. For core-collapse supernova simulations to be carried out by the ExaStar project under the Exascale Computing Project umbrella, a robust, efficient, and state-of-the-art magnetohydrodynamics solver is a critical requirement. In Flash-X, the primary software instrument for ExaStar, a new magnetohydrodynamics solver has been designed and implemented from the ground up to achieve accuracy and efficiency for simulations of complex astrophysical flows. This new solver, dubbed Spark, uses high-order spatial reconstruction, Runge-Kutta time integration, and an efficient cell-centered approach to satisfying the divergence-free condition for the magnetic fields. Spark was written to be optimized for data locality in cache hierarchy of CPUs. Since data locality optimizations for cache hierarchy are not directly compatible with those of accelerators, we have taken the approach of using program synthesis to avoid massive amounts of code replication that would be necessary if we were to maintain two different versions of the solver. Our program synthesis relies on a simple key-dictionary approach, implemented in python, that enables us to assemble the version of the solver suitable for the target hardware from code fragments identified by specific keys. In this work, we describe the data locality optimizations of the solver for CPUs and accelerators and the program synthesis tools that enable this portability. We also detail the parallel performance of Spark for both CPUs and accelerators.

97 MATHEMATICS AND COMPUTING↗

Final Seismic Shake Table Test Plan

The Spent Fuel Waste Disposition (SFWD) program is planning to conduct a full-scale seismic shake table test on the dry storage systems of spent nuclear fuel (SNF) to close the gap related to seismic loads on fuel assemblies in dry storage systems. This test will allow for quantifying the strains and accelerations on surrogate fuel assembly hardware and cladding during earthquakes of different magnitudes and frequency content. Full-scale testing is needed because a dry storage system is a complex and highly nonlinear system making it hard to predict (model) the responses to seismic excitations. The non-linearity arises from the multiple spatial gaps in the system – between fuel rods and the basket, between the basket and dry storage canister, between the dry storage canister and the storage cask (overpack), and ventilation gaps. The non-linearities pose significant limitations on the value of tests with scaled systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A three dimensional multigrid multiblock multistage time stepping scheme for the Navier-Stokes equations

A general multiblock method for the solution of the three-dimensional, unsteady, compressible, thin-layer Navier-Stokes equations has been developed. The convective and pressure terms are spatially discretized using Roe's flux differencing technique while the viscous terms are centrally differenced. An explicit Runge-Kutta method is used to advance the solution in time. Local time stepping, adaptive implicit residual smoothing, and the Full Approximation Storage (FAS) multigrid scheme are added to the explicit time stepping scheme to accelerate convergence to steady state. Results for three-dimensional test cases are presented and discussed.

Elmiligui, Alaa↗

Plasma and Field Observations at the Day-Side, Equatorial Magnetopause, Boundary Layers and Magnetosphere

The Polar spacecraft's orbit has precessed in latitude to an orientation that places it at the dayside magnetopause every 18 hours. In this configuration the various regions near the magnetopause(LLBL, turbulent boundary layer, magnetosphere, and magnetosheath) are sampled with high temporal and spatial resolution. These observational periods-ranging from several minutes to more than two hours-provide an unprecedented look at plasma conditions in these regions. Initial analysis of the low-energy ion data from TIDE reveal plasmaspheric-like ions within the turbulent boundary layer. Within this layer, circularly polarized waves accelerate these ions to 30-40 kilometers per second perpendicular to the local magnetic field. These relatively high velocities allow the H(+) to be observed above the -2V spacecraft potential. They also put the low-density O(+) in the higher-energy, higher sensitivity channels such that densities of order 10e-2 can be observed.

Chandler, M. O.↗

Polar Spacecraft Based Comparisons of Intense Electric Fields and Poynting Flux Near and Within the Plasma Sheet-Tail Lobe Boundary to UVI Images: An Energy Source for the Aurora

In this paper, we present measurements from two passes of the Polar spacecraft of intense electric and magnetic field structures associated with Alfven waves at and within the outer boundary of the plasma sheet at geocentric distances of 4-6 R(sub E), near local midnight. The electric field variations have maximum values exceeding 100 mV/m and are typically polarized approximately normal to the plasma sheet boundary. The electric field structures investigated vary over timescales (in the spacecraft frame.) ranging front 1 to 30 s. They are associated with strong magnetic field fluctuations with amplitudes of 10-40 nT which lie predominantly ill the plane of the plasma sheet and are perpendicular to the local magnetic field. The Poynting flux associated with the perturbation fields measured at these altitudes is about 1-2 ergs per square centimeters per second and is directed along the average magnetic field direction toward the ionosphere. If the measured Poynting flux is mapped to ionospheric altitudes along converging magnetic field lines. the resulting energy flux ranges up to 100 ergs per centimeter squared per second. These strongly enhanced Poynting fluxes appear to occur in layers which are observed when the spacecraft is magnetically conjugate (to within a 1 degree mapping accuracy) to intense auroral structures as detected by the Polar UV Imager (UVI). The electron energy flux (averaged over a spatial resolution of 0.5 degrees) deposited in the ionosphere due to auroral electron beams as estimated from the intensity in the UVI Lyman-Birge-Hopfield-long filters is 15-30 ergs per centimeter squared per second. Thus there is evidence that these electric field structures provide sufficient Poynting flux to power the acceleration of auroral electrons (as well as the energization of upflowing ions and Joule heating of the ionosphere). During some events the phasing and ratio of the transverse electric and magnetic field variations are consistent with earthward propagation of Alfven surface waves with phase velocities of 4000-10000 kilometers per second. During other events the phase shifts between electric and magnetic fields suggest interference between upward and downward propagating Alfven waves. The E/B ratios are about an order of magnitude larger than typical values of C/SIGMA(sub p), where SIGMA(sub p), is the height integrated Pedersen conductivity. The contribution to the total energy flux at these altitudes from Poynting flux associated with Alfven waves is comparable to or larger than the contribution from the particle energy flux and 1-2 orders of magnitude larger than that estimated from the large-scale steady state convection electric field and field-aligned current system.

Wygant, J. R.↗

The Focusing Optics X-Ray Solar Imager: FOXSI

The Focusing Optics x-ray Solar Imager (FOXSI) is a sounding rocket payload funded under the NASA Low Cost Access to Space program to test hard x-ray (HXR) focusing optics and position-sensitive solid state detectors for solar observations. Today's leading solar HXR instrument, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) provides excellent spatial (2 arcseconds) and spectral (1 keV) resolution. Yet, due to its use of an indirect imaging system, the derived images have a low dynamic range (typically <10) and sensitivity. These limitations make it difficult to study faint x-ray sources in the solar corona which are crucial for understanding the particle acceleration processes which occur there. Grazing-incidence x-ray focusing optics combined with position-sensitive solid state detectors can overcome both of these limitations enabling the next breakthrough in understanding impulsive energy release on the Sun. The FOXSI project is led by the Space Sciences Laboratory at the University of California, Berkeley. The NASA Marshall Space Flight Center is responsible for the grazing-incidence optics, while the Astro-H team at JAXA/ISAS has provided double-sided silicon strip detectors. FOXSI is a pathfinder for the next generation of solar hard x-ray spectroscopic imagers. Such observatories will be able to image the non-thermal electrons within the solar flare acceleration region, trace their paths through the corona, and provide essential quantitative measurements such as energy spectra, density, and energy content in accelerated electrons.

Krucker, Saem↗

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗

Robust and optimal alignment of high-dimensional data using maximum likelihood estimation through a random sample consensus framework

Abstract Correcting spatial orientations of groups of high-dimensional data sets such that they are all in a consistent coordinate system is often a time-consuming and error-prone process. Automation of this process can be accomplished by using Generalized Procrustes Analysis to estimate the relative orientations among a population of high-dimensional data sets. A least squares Procrustes solution is applied through a maximum likelihood estimation and random sample consensus framework for robustness. The likelihood model is comprised of a mixture distribution where inliers are modeled using t -distribution and outliers from a uniform distribution. Applications will focus on a synthetic data set that emulates triaxial acceleration data and also real shock data from a population of triaxial accelerometers. Outliers represent either non-rigid body responses, environmental noise, and/or sensor and data acquisition issues. The intended application for the methodology is to robustly automate the rotation of populations of experimentally collected triaxial accelerometer data sets to a single global coordinate system.

LOSAC↗

AI-NERD: Elucidation of relaxation dynamics beyond equilibrium through AI-informed X-ray photon correlation spectroscopy

Abstract Understanding and interpreting dynamics of functional materials in situ is a grand challenge in physics and materials science due to the difficulty of experimentally probing materials at varied length and time scales. X-ray photon correlation spectroscopy (XPCS) is uniquely well-suited for characterizing materials dynamics over wide-ranging time scales. However, spatial and temporal heterogeneity in material behavior can make interpretation of experimental XPCS data difficult. In this work, we have developed an unsupervised deep learning (DL) framework for automated classification of relaxation dynamics from experimental data without requiring any prior physical knowledge of the system. We demonstrate how this method can be used to accelerate exploration of large datasets to identify samples of interest, and we apply this approach to directly correlate microscopic dynamics with macroscopic properties of a model system. Importantly, this DL framework is material and process agnostic, marking a concrete step towards autonomous materials discovery.

36 MATERIALS SCIENCE↗

Laser-Produced Coherent X-ray Sources. Final report

This is the final report for the DOE-AMO DE-FG02-05ER15663 grant. Using a high-peak-power laser system at the Extreme Light Laboratory of the University of Nebraska-Lincoln, we develop novel femtosecond radiation sources and use them to investigate ultrafast photo-induced processes. Pulses of high energy electrons and x-rays are generated when an optical pulse is focused to ultra-high intensity. The radiation is used to probe the evolution of matter under highly non-equilibrium conditions with atomic-scale temporal and spatial resolution. Moreover, the facility is small enough to fit in a university laboratory and is operated by graduate students and postdocs.

74 ATOMIC AND MOLECULAR PHYSICS↗

Observational data needs for plasma phenomena

Bright comets display a rich variety of interesting plasma phenomena which occur over an enormous range of spatial scales, and which require different observational techniques to be studied effectively. Wide-angle photography of high time resolution is probably the best method of studying the phenomenon of largest known scale: the plasma tail disconnection event (DE), which has been attributed to magnetic reconnection at interplanetary sector boundary crossings. These structures usually accelerate as they recede from the head region and observed velocities are typically in the range 50 V km/s. They are often visible for several days following the time of disconnection, and are sometimes seen out past 0.2 AU from the cometary head. The following areas pertaining to plasma phenomena in the ionoshere are addressed: the existence, size, and heliocentric distance variations of the contact surface, and the observational signatures of magnetic reconnection at sector boundary crossings.

Niedner, M. B., Jr.↗

Features of steady magnetospheric convection

The large-scale patterns of ionospheric convection and particle precipitation are described during two intervals of steady magnetospheric convection (SMC) on November 24, 1981. The unique data set used in the analysis includes recordings from the worldwide network of magnetometers and all-sky cameras, global auroral images from the Dynamics Explorer (DE) 1 spacecraft, and particle precipitation data from low-altitude National Oceanic and Atmospheric Administration (NOAA) 6 and NOAA 7 spacecraft. The data show that intense magnetospheric convection continued during more than 10 hours under the steady southward interplanetary magnetic field without any distinct substorm signatures. All data sets available confirmed the stable character of the large-scale magnetospheric configuration during this period. In particular, the magnetic flux threading the polar cap was stable (within 10%) during 3.5 hours of continued DE 1 observations. The dayside cusp was located at an unusually low latitude (70 deg CGL). The nightside auroral pattern consisted of two distinct regions. The diffuse aurora in the equatorward half of the expanded (10 deg wide) auroral oval was well-separated from the bright, active auroral forms found in the vicinity of the poleward boundary of the oval. The twin-vortex convection pattern had no signature of the Harang discontinuity; its nightside 'convection throat' was spatially coincident with the poleward active auroras. This region of the auroral oval was identified as the primary site of the short-lived transient activations during the SMC intervals. The energetic particle observations show that the auroral precipitation up to its high-latitude limit is on closed field lines and that particle acceleration up to greater than 30-keV energy starts close to this limit. The isotropic boundaries of the greater than 30-keV protons and electrons were found close to each other, separating regions of discrete and diffuse precipitation. This suggests that these precipitation types originate on the very taillike and very dipolelike field lines, respectively.

Yahnin, A.↗

On the Formation of Massive Stars

We calculate numerically the collapse of slowly rotating, nonmagnetic, massive molecular clumps of masses 30,60, and 120 Stellar Mass, which conceivably could lead to the formation of massive stars. Because radiative acceleration on dust grains plays a critical role in the clump's dynamical evolution, we have improved the module for continuum radiation transfer in an existing two-dimensional (axial symmetry assumed) radiation hydrodynamic code. In particular, rather than using "gray" dust opacities and "gray" radiation transfer, we calculate the dust's wavelength-dependent absorption and emission simultaneously with the radiation density at each wavelength and the equilibrium temperatures of three grain components: amorphous carbon particles. silicates, and " dirty ice " -coated silicates. Because our simulations cannot spatially resolve the innermost regions of the molecular clump, however, we cannot distinguish between the formation of a dense central cluster or a single massive object. Furthermore, we cannot exclude significant mass loss from the central object(s) that may interact with the inflow into the central grid cell. Thus, with our basic assumption that all material in the innermost grid cell accretes onto a single object. we are able to provide only an upper limit to the mass of stars that could possibly be formed. We introduce a semianalytical scheme for augmenting existing evolutionary tracks of pre-main-sequence protostars by including the effects of accretion. By considering an open outermost boundary, an arbitrary amount of material could, in principal, be accreted onto this central star. However, for the three cases considered (30, 60, and 120 Stellar Mass originally within the computation grid), radiation acceleration limited the final masses to 3 1.6, 33.6, and 42.9 Stellar Mass, respectively, for wavelength-dependent radiation transfer and to 19.1, 20.1, and 22.9 Stellar Mass. for the corresponding simulations with gray radiation transfer. Our calculations demonstrate that massive stars can in principle be formed via accretion through a disk. The accretion rate onto the central source increases rapidly after one initial free-fall time and decreases monotonically afterward. By enhancing the nonisotropic character of the radiation field, the accretion disk reduces the effects of radiative acceleration in the radial direction - a process we call the "flashlight effect." The flashlight effect is further amplified in our case by including the effects of frequency-dependent radiation transfer. We conclude with the warning that a careful treatment of radiation transfer is a mandatory requirement for realistic simulations of the formation of massive stars.

Yorke, Harold W.↗

Understanding the eco-geomorphologic feedback of coastal marsh under sea level rise: vegetation dynamic representations, processes interaction, and parametric sensitivity

A growing number of coastal eco-geomorphologic modeling studies have been conducted to understand coastal marsh evolution under sea-level rise (SLR). Although these models quantify marsh topographic change as a function of sedimentation and erosion, their representations of vegetation dynamics that control organic sedimentation differ. How vegetation dynamic schemes contribute to simulation outcomes is not well quantified. Additionally, the sensitivity of modeling outcomes to parameter selection in the available formulations has not been rigorously tested to date, especially under the influence of an accelerating SLR. In this paper, we used a coastal eco-geomorphologic model with different vegetation dynamic schemes to investigate the eco-geomorphologic feedbacks of coastal marshes and parametric sensitivity under SLR scenarios. We found that marsh platform relief increased with sea level rise rate. The simulations with different vegetation schemes exhibited different spatial-temporal variations in elevation and biomass. The nonlinear Spartina scheme presented the most resilient prediction with generally the highest marsh accretion and vegetation biomass, and the least elevation relief under SLR. But the linear Spartina scheme predicts the lowest unvegetated-vegetated ratio. We also found that vegetation-related parameters and sediment diffusivity, which were not well measured or discussed in previous studies, were identified as some of the most critical parameters. Additionally, the model sensitivity to vegetation-related parameters increased with SLR rates. The identified most sensitive parameters may inform how to appropriately choose modeling representations of key processes and parameters for different coastal marsh landscapes under SLR, and demonstrate the importance of future field measurements of these key parameters.

58 GEOSCIENCES↗

The acceleration of charged particles in interplanetary shock waves

Consideration of the theoretical and observational literature on energetic ion acceleration in interplanetary shock waves is the basis for the present discussion of the shock acceleration of the solar wind plasma and particle transport effects. It is suggested that ISEE data be used to construct data sets for shock events that extend continuously from solar wind to galactic cosmic ray energies, including data for electrons, protons, alphas and ions with Z values greater than 2.0, and that the temporal and spatial evolution of two- and three-dimensional particle distribution functions be studied by means of two or more spacecraft.

Pesses, M. E.↗

Thermal and nonthermal phenomena in solar flare loops at 20 cm wavelength and in X-rays

X-ray images from the P78-1 satellite are presented for a long-lasting burst at 20 cm wavelength mapped with the Very Large Array on May 19, 1979. The decimeter wave observations were originally interpreted in terms of two models, one invoking thermal electrons radiating at low harmonics of the gyrofrequency, and the other invoking mildly relativistic electrons emitting gyrosynchrotron radiation. If the 20 cm source is thermal, it should also be visible in soft X-rays, while if it is nonthermal, the soft X-ray emission should be weak or spatially or temporally distinct from the 20 cm burst. It is found that only one of the three 20 cm sources was approximately co-spatial with the soft X-ray source, and that it was only partially thermal. The 20 cm burst is therefore primarily decimeter type IV emission from mildly relativistic electrons of the post-flare phase. The over 2 hr lifetime and smooth temporal variation of the burst belie its nonthermal nature and suggest continuous acceleration as well as long-term storage of energetic electrons.

Schmahl, E. J.↗

The Advanced Energetic Pair Telescope for Gamma-Ray Polarimetry

The Advanced Energetic Pair Telescope (AdEPT), a future NASA/GSFC MIDEX mission, is being developed to perform high-sensitivity medium-energy (5-200 MeV) astronomy and revolutionary gamma-ray polarization measurements. The enabling technology for AdEPT is the Three- Dimensional Track Imager (3-DTI), a large volume gaseous time projection chamber with 2- dimentional micro-well detector (MWD) readout. The low density and high spatial resolution of the 3-DTI allows AdEPT to achieve high angular resolution (~0.5 at 67.5 MeV) and, for the first time, exceptional gamma-ray polarization sensitivity. These capabilities enable a wide range of scientific discovery potential for AdEPT. The key science goals of the AdEPT mission include: 1) Explore fundamental processes of particle acceleration in active astrophysical objects, 2) Reveal the magnetic field configuration of the most energetic accelerators in the Universe, 3) Explore the origins and acceleration of cosmic rays and the Galactic MeV diffuse emission, 4) Search for dark matter in the Galactic center, and 5) Test relativity with polarization measurements. We report on the latest developments of the MWDs for the 3-DTI.

Hunter, Stanley D.↗

Snow Albedo Feedbacks Enhance Snow Impurity-Induced Radiative Forcing in the Sierra Nevada

This study employs a fully coupled meteorology-chemistry-snow model to investigate the impacts of light-absorbing particles (LAPs) on snow darkening in the Sierra Nevada. After comprehensive evaluation with spatially and temporally complete satellite retrievals, the model shows that LAPs in snow reduce snow albedo by 0.013 (0-0.045) in the Sierra Nevada during the ablation season (April-July), producing a midday mean radiative forcing (RF) of 4.5 W m -2 which increases to 15-22 W m -2 in July. LAPs in snow accelerate snow aging processes and reduce snow cover fraction, which doubles the albedo change and radiative forcing caused by LAPs. The impurity-induced snow darkening effects decrease snow water equivalent and snow depth by 20 mm and 70 mm in June in Sierra Nevada bighorn sheep habitat. The earlier snowmelt reduces root zone soil moisture by 20%, deteriorating the forage productivity and playing a negative role in the survival of bighorn sheep.

54 ENVIRONMENTAL SCIENCES↗