Search NASA⌕ Search

SEARCH · Search NASA

Results for “Processing maps”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Deformation Twinning of a Silver Nanocrystal under High Pressure

Within a high-pressure environment, crystal deformation is controlled by complex processes such as dislocation motion, twinning, and phase transitions, which change materials' microscopic morphology and alter their properties. Understanding a crystal's response to external stress provides a unique opportunity for rational tailoring of its functionalities. It is very challenging to track the strain evolution and physical deformation from a single nanoscale crystal under high-pressure stress. Here, we report an in situ three-dimensional mapping of morphology and strain evolutions in a single-crystal silver nanocube within a high-pressure environment using the Bragg Coherent Diffractive Imaging (CDI) method. We observed a continuous lattice distortion, followed by a deformation twining process at a constant pressure. The ability to visualize stress-introduced deformation of nanocrystals with high spatial resolution and prominent strain sensitivity provides an important route for interpreting and engineering novel properties of nanomaterials.

coherent diffractive imaging↗
DOE OSTI2018

Mechanism of Na{sub 2}SO{sub 4} Promoting Nickel Extraction from Sulfide Concentrates by Sulfation Roasting–Water Leaching

As a more environmentally friendly and energy-efficient route, the sulfation roasting–water leaching technique has been developed for highly effective extraction of non-ferrous metals from nickel sulfide concentrate in the presence of a Na{sub 2}SO{sub 4} additive. The effects of several important roasting parameters—the roasting temperature, the addition of Na{sub 2}SO{sub 4}, the holding time, and the heating rate in particular—have been investigated. The results suggest that about 90 pct Ni, 92 pct Co, 95 pct Cu, and < 1 pct Fe can be leached from the calcine roasted under the optimum conditions. Furthermore, the behavior and mechanism of the Na{sub 2}SO{sub 4} additive in the roasting process have been well addressed by detailed characterization of the roasted product and leaching residue using quantitative phase analysis (QPA) and energy dispersive spectroscopy (EDS) mapping. The Na{sub 2}SO{sub 4} additive was observed to play a noticeable role in promoting the sulfation degree of valuable metals by forming liquid phases [Na{sub 2}Me(SO{sub 4}){sub 2}] at the outermost layer, which can create a suitable dynamic environment for sulfation. Thus, addition of Na{sub 2}SO{sub 4} might be conducive to an alternative metallurgical process involving complex sulfide ores.

36 MATERIALS SCIENCE↗
DOE OSTI2016

A butterfly diagram and carrington maps for century-long Ca II K spectroheliograms from the Kodaikanal observatory

The century-long (1907–2007) Ca ii K spectroheliograms from the Kodaikanal Solar Observatory (KSO) are calibrated, processed, and analyzed to follow the evolution of the bright on-disc structures called plages, possible representatives of magnetic activity on the Sun. This is the longest data set studied in Ca ii K to date, covering about 9.5 cycles of 11 yr periods. Plages are segmented with area ⩾1 arcmin{sup 2} using global thresholds for individual full disc images and subsequent application of a morphological closing operation. The plage index is calculated and is seen to have a close positive correlation with the fractional disc area covered by plages. The newly generated plage area cycle (from KSO) was compared with the same from the Mount Wilson Observatory (correlation 95.6%) for the overlapping years, i.e., 1915–2000. This study illustrates the time–latitude distribution of plage centroids by rendering a butterfly diagram (as observed for sunspots). The 3D visualization of the diagram shows one-to-one mapping between plage location, time, and area. This work further delineates the positional correlation between magnetic patches and plage regions through the comparison of synoptic maps derived from both KSO Ca ii K images and space-based full disc line-of-sight magnetograms. Regular synoptic magnetograms from ground-based observatories are available only after 1970s. Thus the long term Ca ii K data from KSO can be used as a proxy for estimating magnetic activity locations and their strengths at earlier times.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Particle-image Velocimetry Sensitivity through Automatic Differentiation

A commonly encountered difficulty in modern experimental fluid dynamics experiments is the need for full-flow visualization and measurement; traditional techniques such as hot-wire probes and laser Doppler velocimetry cannot provide such information without difficult scanning procedures. Particle-image velocimetry (PIV) is a measurement technique which can overcome this spatial limitation to provide real-time three-dimensional fluid velocity vectors for an entire plane within the flow, as opposed to a single point that is provided by other measurement techniques. To accomplish this, a PIV system seeds the fluid flow with tracer particles; these must have nearly the same density as that of the fluid and a high drag coefficient to ensure that they follow the fluid flow without deviation. A laser is sent through optics to expand its beam into a sheet which illuminates the plane of interest in the fluid. The particles reflect more light than the fluid, and their image is recorded using a high speed camera. Due to limitations in bandwidth, video cannot be taken fast enough to be useful, and instead pairs of images are captured with a known separation in time. The illuminated particles will move between the two images taken by the camera, and their relative motion can be used to compute the fluid velocity at that point. When a single camera is used, only two-dimensional velocity vectors can be computed, but with the addition of another camera this can be improved to full three-dimensional vectors. The flow's time evolution can be measured by taking series of image pairs and processing each pair to get the velocity field at the time between each image pair. Only the two-dimensional version will be considered in this work. To compute the velocity vectors of the flow, the images are broken up into interrogation regions, usually in a uniform grid. In the most basic PIV algorithm, corresponding regions in the image pairs are cross-correlated to compute the most likely displacement from one image to the next. This is used since the cross correlation will achieve a maximum value when the two interrogation regions line up the best. The location of the peak in the correlation map and scaling information for the optical setup is used to compute the particle displacement from one image to the next. This process is similar to holding two semi-transparent images and moving one until they match as well as possible. Sub-pixel displacements between images can be obtained by fitting a Gaussian function to the correlation peak. Knowing the displacement vector for each set of interrogation regions in the image pair, it is a simple matter to divide by the time gap between the two images to compute the fluid velocity at the point in time and space. Since there are many interrogation regions in each image pair, the velocity field for a sheet is sampled according to the arrangement of the interrogation regions. The sensitivity of the final computed vector to the input images is difficult to compute, and has led most considerations of uncertainty in the PIV field to apply Monte Carlo methods and correlations to quantify the systematic and random errors which result from application of various PIV algorithms. In this work, automatic differentiation (AD) has been used to trace the dependence of velocity vectors computed using the correlation-peak PIV algorithm back to the input values given to an artificial image generator. These initial results form the foundation for continued work in the application of AD to algorithms in PIV. They demonstrate that AD methods can be successfully applied to PIV algorithms without undue complexity. The derivatives of each individual vector with regard to the coefficients of the velocity field used to create the artificial images can be computed and the statistics of them considered collectively, or correlated between each other to quantify the effectiveness of a PIV algorithm. With the initial implementation complete, adding extra algorithms is relatively simple and considerable data can be generated for analysis. It is hoped that application of these methods can lead to improved PIV algorithms which demonstrate better sensitivity to small variations in particle displacement. While these methods can be used as an analysis tool for refinement of new algorithms, computation of actual resultant vector sensitivities to input images' individual pixels is the next logical step of this work. Such information will be analyzed with experimental data to produce a better understanding of how output vectors' uncertainty relates to the actual image pairs, and lead to automatic bad vector removal. (authors)

42 ENGINEERING↗
DOE OSTI2013

High-Throughput Mapping of the Electrochemical Properties of (Ni-Fe-Co-Ce)O x Oxygen-Evolution Catalysts

Discovering improved electrocatalysts is critical for many technologically important processes and for the development of new clean-energy technologies. In this paper, high-throughput methods for measuring fundamental electrochemical properties are demonstrated through the investigation of oxygen-evolution catalysis by using 665 oxide compositions containing nickel, iron, cobalt, and cerium. The behavior of each composition is characterized in 1.0 M NaOH(aq) by using a scanning drop three-electrode cell to perform chronopotentiometry (CP) and cyclic voltammetry experiments. CP measurements at different current densities identify different composition–performance trends, owing to underlying variations in fundamental electrochemical behavior. We report systematic, coincident, composition-dependent trends in the Tafel slopes and the reversible redox potentials of the catalysts. Finally, applying high-throughput electrochemical methods provides insight into composition–property–performance relationships and motivates new directions for the study of catalyst mechanisms by using informatics and theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mapping the Temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, IR-IR photodissociation mass spectrometry to the D3O+?(HDO)(D2O)19 isotopologue of the “magic number” protonated water cluster, H+?(H2O)21. The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure. The OH frequency is observed to evolve over tens of milliseconds in the temperature range (90-120 K). Starting at 100 K, large “jumps” are observed between two OH frequencies separated by ~300 cm-1 indicating migration of the OH group from the bound OH site at 3350 cm-1 to the free position at 3686 cm-1. Increasing the temperature to 110 K leads to partial interconversion among many sites. All sites are observed to interconvert at 120 K such that the distribution of the unique OH group among them adopts the form one would expect for a canonical ensemble. The spectral dynamics displayed by the clusters thus offer an unprecedented view into the molecular-level processes that drive spectral diffusion in an extended network of water molecules.

Yang, Nan↗
DOE OSTI2020

Mapping the Concentration Distribution of Contaminant Plumes to the Computational Grid of the Plateau to River Model (Ver. 8.3)

This environmental calculation file (ECF) outlines the process conducted for defining the initial condition groundwater plumes for select contaminants of interest (COIs) to the Plateau River model version 8.3 (P2Rv8.3) domain (CP-57037, Model Package Report: Plateau to River Groundwater Model Version 8.3). Various data sources were combined into a single set of files defining plume conditions at time 1/1/2018. The resulting files are intended for use as inputs for transport modeling to support forecasts of the Central Plateau contaminant plume fate and transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Dereplicating and spatial mapping of secondary metabolites from fungal cultures in situ

Ambient ionization techniques coupled to mass spectrometry have recently become prevalent in natural product research due to their ability to examine secondary metabolites in situ. Identifying, mapping, and monitoring secondary metabolites directly on an organism provides invaluable spatial and temporal details that are lost through traditional extraction processes. Most ambient ionization techniques do not collect mutually supportive data, such as chromatographic retention times and/or UV/VIS spectra, and this can limit the ability to identify certain metabolites, such as differentiating isomers. To overcome this, the droplet liquid microjunction surface sampling probe (droplet LMJ SSP) was coupled with UPLC PDA HRMS MS/MS, thus providing separation, retention times, and UV/VIS data used in traditional dereplication protocols. By capturing these mutually supportive data, the identity of secondary metabolites could be confidently and rapidly assigned in situ. Using the droplet LMJ SSP, a protocol was constructed to analyze the secondary metabolite profile of fungal cultures directly without any sample preparation. The results demonstrate that fungal cultures can be dereplicated from the Petri dish, thus identifying secondary metabolites, including isomers, and confirming them against reference standards. As a result, heat maps, similar to mass spectrometry imaging, can be used to ascertain the location and relative concentration of secondary metabolites directly on the surface and/or surroundings of a fungal culture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Acquisition, preprocessing, and reconstruction of ultralow dose volumetric CT scout for organ-based CT scan planning

Purpose: Traditionally, 2D radiographic preparatory scan images (scout scans) are used to plan diagnostic CT scans. However, a 3D CT volume with a full 3D organ segmentation map could provide superior information for customized scan planning and other purposes. A practical challenge is to design the volumetric scout acquisition and processing steps to provide good image quality (at least good enough to enable 3D organ segmentation) while delivering a radiation dose similar to that of the conventional 2D scout. Methods: The authors explored various acquisition methods, scan parameters, postprocessing methods, and reconstruction methods through simulation and cadaver data studies to achieve an ultralow dose 3D scout while simultaneously reducing the noise and maintaining the edge strength around the target organ. Results: In a simulation study, the 3D scout with the proposed acquisition, preprocessing, and reconstruction strategy provided a similar level of organ segmentation capability as a traditional 240 mAs diagnostic scan, based on noise and normalized edge strength metrics. At the same time, the proposed approach delivers only 1.25% of the dose of a traditional scan. In a cadaver study, the authors’ pictorial-structures based organ localization algorithm successfully located the major abdominal-thoracic organs from the ultralow dose 3D scout obtained with the proposed strategy. Conclusions: The authors demonstrated that images with a similar degree of segmentation capability (interpretability) as conventional dose CT scans can be achieved with an ultralow dose 3D scout acquisition and suitable postprocessing. Furthermore, the authors applied these techniques to real cadaver CT scans with a CTDI dose level of less than 0.1 mGy and successfully generated a 3D organ localization map.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Autonomous 3D and Radiation Mapping in Tunnel Environments Using Aerial Robots - 18156

In this paper we present the algorithms and systems developed to enable autonomous 3D and radiation mapping, as well as radiological source localization in tunnel environments using aerial robots. The developed technologies aim to address the need of autonomous robotic inspection of nuclear sites such as facilities in the process of being decommissioned. As a motivational example, the PUREX tunnels at the Hanford Site are considered and the robot is designed to be able to conduct multi-modal 3D/radiation characterization of the tunnel environment without any prior knowledge or map provided beforehand. In order to enable the necessary GPS-denied simultaneous localization and mapping in dark, visually-degraded tunnel environments, a multi-modal mapping unit was developed which fuses the data from stereo vision cameras synchronized with flashing LEDs, inertial sensors, and three time-of-flight 3D depth sensors. The robot's pose and surrounding map are estimated at 20 Hz through a depth-enhanced visual- inertial localization and mapping framework. Provided this ability to estimate the robot's pose and develop a dense map in real-time, a receding horizon sampling-based autonomous exploration path planning algorithm was utilized and equips the robot with the ability to explore its environment given no prior map. The aforementioned key robotic functionalities are further extended through the integration of miniaturized fully solid-state scintillation detectors, utilizing advanced scintillators and Silicon Photomultipliers. The derived sensor data is employed for radiation mapping, which when fused with the 3D maps and the spectroscopic source localization methods, allow the aerial robotic platform to localize the radiation source and identify the specific source isotope. The source location is estimated through a geometric algorithm which can localize a source given only three measurements and short dwell times. The developed systems were tested using two small aerial robot platforms of the F450 and F550 frame classes with weights of 2.6 kg and 2.9 kg respectively. All localization and planning as well as high-level control algorithms are executed onboard. The scintillation detector is integrated on the robot and its data is collected through a dedicated microcontroller. A set of field experiments were conducted inside the dark and GPS-denied environment of a 4.25 m wide, 7.5 m tall, and 137 m long railroad tunnel. A Cs-137 source was placed within the tunnel, the associated 3D/radiation maps were derived, and the source location was estimated. Through these efforts, the goal is to develop the robotic technology that would enable the autonomous systematic inspection or emergency response in nuclear waste decommissioning facilities. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Behavioral simulations and experimental evaluations of stress induced spatial nonuniformity of dynamic bulk modulus in additive friction stir deposited AA 6061

Additive friction stir deposition (AFSD) is a rapidly developing solid-state additive manufacturing (AM) technology that enables the fabrication of customized products on a relatively large scale. The unique forge-liked process offers high quality products with refined equiaxed microstructure and exceptional mechanical performance. However, in AFSD, the spatially distributed mechanical properties within a deposited layer have not been established so far, due to the lack of a suitable large-scale mapping method. In the present work, an ultrasonic dynamic-elasticity mapping method was employed to map the distribution of physical properties in a AFSD 6061 Al. Strong spatial asymmetries in elastographies were mainly due to the residual stresses. The origin of the asymmetries was investigated through thermo-mechanical simulations of the AFSD process. It was observed that majority of the asymmetries were associated with the mechanical residual stresses instead of the symmetric thermal residual stresses during linear deposition of the AFSD process. However, the thermal residual stresses can be manipulated to become asymmetric with a 90° turn in the linear deposition path, due to the difference between the levels of cooling along the edges of the inner and outer corners. Hence, the mechanical and thermo-mechanical residual stresses can be additive or subtractive depending on the direction of the turns. Here, such numerically predicted behaviors were also experimentally observed on the AFSD fabricated samples by ultrasonic elastography. Furthermore, additional preliminary studies on multi-layer AFSD fabrications show the asymmetric mechanical stresses generated during the deposition process on different layers generate constructive or destructive interference which results to spatially variable dynamic bulk modulus depending on the deposition orientations of layers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Autonomous scanning probe microscopy investigations over WS2 and Au{111}

Abstract Individual atomic defects in 2D materials impact their macroscopic functionality. Correlating the interplay is challenging, however, intelligent hyperspectral scanning tunneling spectroscopy (STS) mapping provides a feasible solution to this technically difficult and time consuming problem. Here, dense spectroscopic volume is collected autonomously via Gaussian process regression, where convolutional neural networks are used in tandem for spectral identification. Acquired data enable defect segmentation, and a workflow is provided for machine-driven decision making during experimentation with capability for user customization. We provide a means towards autonomous experimentation for the benefit of both enhanced reproducibility and user-accessibility. Hyperspectral investigations on WS 2 sulfur vacancy sites are explored, which is combined with local density of states confirmation on the Au{111} herringbone reconstruction. Chalcogen vacancies, pristine WS 2 , Au face-centered cubic, and Au hexagonal close-packed regions are examined and detected by machine learning methods to demonstrate the potential of artificial intelligence for hyperspectral STS mapping.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Model Data Archive for Manuscript Titled "Evaluation of a Coupled Surface–Subsurface Hydrologic Model Using Dense Water‑Level Sensors in a Mixed Urban–Rural Watershed"

This archive provides scripts, input files, and datasets used for the implementation and evaluation of a fully coupled surface–subsurface hydrologic model in the Neches River Basin, southeast Texas. The study uses the Advanced Terrestrial Simulator (ATS) to simulate coupled surface–subsurface hydrologic processes over a mixed urban–rural watershed and evaluates model performance using a dense network of 136 in situ water-level sensors, nine U.S. Geological Survey (USGS) stream gauges, and SSEBop-derived evapotranspiration estimates during the period October 2014–June 2024. The workflow is implemented primarily in Python 3 using the Watershed Workflow package. The Jupyter notebooks can be executed using open-source software such as Anaconda JupyterLab or Visual Studio Code. Other data files include TXT, CSV, XML, SHP, TIF, NetCDF, HDF5, and ExodusII files, which can be processed using the provided Python scripts. ATS input files are provided in XML format and can be edited using any commonly used text editor. This archive contains: *Scripts and input files used to generate the ATS model setup, including watershed discretization, mesh generation, parameter mapping, and model configuration. *Jupyter notebooks used for preprocessing observational data, evaluating streamflow, water levels, and evapotranspiration, computing performance metrics, and generating the figures presented in the manuscript. *ATS simulation outputs and processed observational datasets, including OneRain and DD6 water-level sensors, USGS streamflow observations, GIS data, and supporting spatial datasets used throughout the study.

Dense water-level sensor network↗
DOE OSTI2013

Recent Improvement of Measurement Instrumentation to Supervise Nuclear Operations and to Contribute Input Data to 3D Simulation Code - 13289

The CEA has developed many compact characterization tools to follow sensitive operations in a nuclear environment. Usually, these devices are made to carry out radiological inventories, to prepare nuclear interventions or to supervise some special operations. These in situ measurement techniques mainly take place at different stages of clean-up operations and decommissioning projects, but they are also in use to supervise sensitive operations when the nuclear plant is still operating. In addition to this, such tools are often associated with robots to access very highly radioactive areas, and thus can be used in accident situations. Last but not least, the radiological data collected can be entered in 3D calculation codes used to simulate the doses absorbed by workers in real time during operations in a nuclear environment. Faced with these ever-greater needs, nuclear measurement instrumentation always has to involve on-going improvement processes. Firstly, this paper will describe the latest developments and results obtained in both gamma and alpha imaging techniques. The gamma camera has been used by the CEA since the 1990's and several changes have made this device more sensitive, more compact and more competitive for nuclear plant operations. It is used to quickly identify hot spots, locating irradiating sources from 50 keV to 1500 keV. Several examples from a wide field of applications will be presented, together with the very latest developments. The alpha camera is a new camera used to see invisible alpha contamination on several kinds of surfaces. The latest results obtained allow real time supervision of a glove box cleaning operation (for {sup 241}Am contamination). The detection principle as well as the main trials and results obtained will be presented. Secondly, this paper will focus on in situ gamma spectrometry methods developed by the CEA with compact gamma spectrometry probes (CdZnTe, LaBr{sub 3}, NaI, etc.). The radiological data collected is used to quantify the activity of hot spots and can also then be entered in 3D models of nuclear plants to simulate intervention scenarios. Recent developments and results will be presented regarding this. Finally, thanks to a large amount of feedback, the interest of using complementary measurements will be discussed. In fact, the recent use of 3D simulation codes requires very accurate knowledge of nuclear plant radiological data. The use of coupled devices such as imaging devices, (gamma and alpha cameras), gamma spectrometry, dose rate mapping, collimated / un-collimated measurements and many other physical values gives an approach to the radiological knowledge of a process or plant with the lowest possible uncertainty. In line with this, the paper will conclude with the future developments and trials that could be assessed in that field of application. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

The G2 erosion model: An algorithm for month-time step assessments

Highlights: • G2 is a soil erosion model for developing monthly erosion maps at regional scale. • G2 enhances the spatio-temporal variability of cover management and erosivity factors. • The V-factor models the impact of vegetation cover, land use and imperviousness layer. • G2 uses data from ESDAC, LUCAS, COPERNICUS, CORINE, SENTINEL and EU-DEM databases. • G2 has been applied in Strymonas Catchment, Crete, Cyprus and Albania catchments. A detailed description of the G2 erosion model is presented, in order to support potential users. G2 is a complete, quantitative algorithm for mapping soil loss and sediment yield rates on month-time intervals. G2 has been designed to run in a GIS environment, taking input from geodatabases available by European or other international institutions. G2 adopts fundamental equations from the Revised Universal Soil Loss Equation (RUSLE) and the Erosion Potential Method (EPM), especially for rainfall erosivity, soil erodibility, and sediment delivery ratio. However, it has developed its own equations and matrices for the vegetation cover and management factor and the effect of landscape alterations on erosion. Provision of month-time step assessments is expected to improve understanding of erosion processes, especially in relation to land uses and climate change. In parallel, G2 has full potential to decision-making support with standardised maps on a regular basis. Geospatial layers of rainfall erosivity, soil erodibility, and terrain influence, recently developed by the Joint Research Centre (JRC) on a European or global scale, will further facilitate applications of G2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

Computational Particle Dynamic Simulations on Multicore Processors (CPDMu) Final Report Phase I

Statement of Problem - Department of Energy has many legacy codes for simulation of computational particle dynamics and computational fluid dynamics applications that are designed to run on sequential processors and are not easily parallelized. Emerging high-performance computing architectures employ massively parallel multicore architectures (e.g., graphics processing units) to increase throughput. Parallelization of legacy simulation codes is a high priority, to achieve compatibility, efficiency, accuracy, and extensibility. General Statement of Solution - A legacy simulation application designed for implementation on mainly-sequential processors has been represented as a graph G. Mathematical transformations, applied to G, produce a graph representation {und G} for a high-performance architecture. Key computational and data movement kernels of the application were analyzed/optimized for parallel execution using the mapping G {yields} {und G}, which can be performed semi-automatically. This approach is widely applicable to many types of high-performance computing systems, such as graphics processing units or clusters comprised of nodes that contain one or more such units. Phase I Accomplishments - Phase I research decomposed/profiled computational particle dynamics simulation code for rocket fuel combustion into low and high computational cost regions (respectively, mainly sequential and mainly parallel kernels), with analysis of space and time complexity. Using the research team's expertise in algorithm-to-architecture mappings, the high-cost kernels were transformed, parallelized, and implemented on Nvidia Fermi GPUs. Measured speedups (GPU with respect to single-core CPU) were approximately 20-32X for realistic model parameters, without final optimization. Error analysis showed no loss of computational accuracy. Commercial Applications and Other Benefits - The proposed research will constitute a breakthrough in solution of problems related to efficient parallel computation of particle and fluid dynamics simulations. These problems occur throughout DOE, military and commercial sectors: the potential payoff is high. We plan to license or sell the solution to contractors for military and domestic applications such as disaster simulation (aerodynamic and hydrodynamic), Government agencies (hydrological and environmental simulations), and medical applications (e.g., in tomographic image reconstruction). Keywords - High-performance Computing, Graphic Processing Unit, Fluid/Particle Simulation. Summary for Members of Congress - Department of Energy has many simulation codes that must compute faster, to be effective. The Phase I research parallelized particle/fluid simulations for rocket combustion, for high-performance computing systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Aerial Lidar and Geomorphic Change Detection Software- Tools for Monitoring Canyon Sediments at LANL - 18668

The Los Alamos National Laboratory Environmental Remediation program has conducted extensive geomorphic surveys of sediment deposits in the canyon bottoms within and surrounding the facility. The goal of the geomorphic surveys was to provide basic information needed to characterize chemical inventories of the sediments downstream of historical release sites. Sediment units were mapped based on depositional setting and age. Representative characterization samples were collected for each unit and submitted for chemical analysis. Chemical inventories were estimated by calculating the total volume of each depositional unit and multiplying by average chemical concentration. This method of characterization provides a spatial distribution of chemical data based on fluvial processes. In general, for the canyon bottom sediments, the higher concentrations of industrial chemicals were associated with fine-grained floodplain deposits closest to the source. Since initial characterization, many of the reaches have undergone changes in response to post-wildfire floods and severe rain events. These events impacted the channel and floodplain deposits. Much of the elevated concentrations are confined to the over-bank deposits distal from the primary channel. In order to ensure that the sediment units containing elevated chemical concentrations are stable, repeat aerial light detection and ranging (LiDAR) surveys are conducted and changes are assessed using the Geomorphic Change Detection (GCD) plug-in for ArcGIS. These tools provide high confidence uncertainty estimates and extensive areal coverage. This is a valuable method for long-term monitoring of surface sediments with low-level but widespread contaminant distributions. All of the monitored canyon bottom reaches pass risk assessment and present technically infeasible sediment removal options. Following the 2014 monsoon season, aerial LiDAR data has been collected annually and the resulting digital elevation models (DEMs) processed using the GCD method. Changes are identified using GCD and are verified by field walk downs. If erosion is identified in areas which are vulnerable to mobilizing contaminated sediment, enhanced controls can be installed to stabilize banks or other areas with high erosion potential. Annual sediment monitoring samples are collected where there is evidence for flow and deposition within the past year. Since 2014, no severe floods have been catalogued using the LiDAR surveys and GCD software. The recent rainfall/runoff events have not caused major changes to the canyon bottom channels and floodplains. Geomorphic processes identified in the monitoring results include channel aggradation and incision that over the course of a monsoon season results in minor changes to the system. Other active processes contributing to observed changes are minor slides, flows, slumps, and falls of unconsolidated sediment on steep bedrock or soil surfaces. A comprehensive baseline has been established by geomorphic mapping and characterization sampling. Periodic monitoring consists of LiDAR surveys and GCD analysis which offers an effective and defensible methodology for documenting stability of surface sediments. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗