Search NASA⌕ Search

SEARCH · Search NASA

Results for “Batch processing”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Testing Destruction of Clay Minerals By Silica-Poor Brines in Gale Crater, Mars Using Laboratory and Field Experiments

The Mars Science Laboratory (MSL) Curiosity rover landed in Gale crater to investigate the cause for mineralogical changes seen from orbit, including the clay-sulfate mineral transition observed on its central mound, Aeolis Mons. The lower strata show spectral signatures of clay minerals, hematite, and hydrated Mg sulfates, and the upper strata are dominated by hydrated Mg sulfate signatures. This significant transition has been observed in other locations on Mars, and therefore processes related to the clay mineral-sulfate transition are likely important broadly across Mars. The bulk mineralogy of ~3.7 Ga sedimentary rocks in Gale crater is a product of variations in sediment sources and depositional processes, with overprinting by episodes of early and late diagenetic alteration. Sedimentary analyses of the smectite clay minerals observed in the Glen Torridon (GT) region of Mount Sharp suggest formation in a lacustrine aqueous environment. Examination of the Vera Rubin ridge (VRR) strata indicate that it was formed at the same time as GT and later altered by extensive diagenesis. Recent work suggests that the observed mineral transition from smectites and talc-serpentinites in GT to iron oxides and oxyhydroxides and VRR may have resulted from extensive diagenetic alteration by density driven, silica-poor brines sourced from the overlying sulfur-bearing unit on Mount Sharp. This process of alteration by silica-poor brines is observed in limited locations on Earth. We are therefore working to demonstrate the impacts of silica poor brines on clay mineral alteration in laboratory column and terrestrial analog burial experiments. Column experiments allow the investigation of dissolution and precipitation under conditions that represent a closer analog to natural environments, while burial experiments allow examination of weathering reactions of a known starting material after a specific period of time in a natural field environment. Both methods are valuable for interpreting the water-rock interactions between clay minerals and brines. In preparation for these longer-term column and burial experiments, we are performing preliminary batch dissolution experiments on materials relevant to Glen Torridon: nontronite and serpentinite. We are using magnesium sulfate brine solutions relevant to those proposed on Mars. The results of this study can provide constraints on the characteristics of the diagenetic fluids that formed the altered minerals in Gale crater. They can also provide insights into the role of brine density-driven diagenesis in past aqueous environments on Mars, which have been proposed to be widespread.

J E Valdueza↗
DOE OSTI2026

Carbon Tetrachloride Degradation Results for 200-ZP-1 Operable Unit

Carbon tetrachloride (CT) contamination in the 200-ZP-1 Operable Unit (OU) at the Hanford Site originated from large-volume discharges to the subsurface during plutonium production operations between 1955 and 1973. Contamination migrated through more than 70 meters of unsaturated sediment to reach the underlying unconfined aquifer, where it persists as a large and complex groundwater plume. The 200-ZP-1 OU Record of Decision (ROD) requires that groundwater CT concentrations be reduced to 3.4 µg/L within 125 years. Current groundwater modeling projections estimate that the existing pump-and-treat, even when combined with monitored natural attenuation (specifically hydrolysis), will not achieve this target within the designated timeframe. A fundamental contributor to this shortfall is the extremely slow rate of CT hydrolysis under Hanford aquifer conditions, which has been estimated to have a half-life of 630 years. If faster-acting biotic and abiotic degradation processes are operating within the aquifer, their contribution to CT mass reduction could have a meaningful impact. However, site-specific measurements of these processes and their rates have not previously been performed. This report documents the results of a two-phase laboratory investigation designed to characterize and quantify the capacity of site-specific 200-ZP-1 OU sediments and groundwater to support natural attenuation of CT through biotic and abiotic pathways. In this context, degradation capacity is defined as the intrinsic potential of the subsurface matrix to transform CT under optimized, controlled conditions. System capacity is evaluated in two ways: (1) as rate-limited capacity, which establishes the maximum kinetic velocity of CT transformation and is measured using half-lives and first order rate constants; and (2) as mass limited capacity, which defines the total contaminant mass the batch experimental system can degrade before reactants are exhausted, representing the maximum amount of contaminant the microbial community and reactive mineral phases can transform under the experimental conditions.

abiotic degradation↗
DOE OSTI2026

Affinity and Treatment of LDR Organics in Low Activity Waste by Cementitious Materials – 26109

The physiochemical interactions between organic contaminants, minerals, and polymer additives in grout/cementitious materials are being identified. This research supports a proposed Resource Conservation and Recovery Act (RCRA) treatment variance relevant to the Land Disposal Restriction (LDR) organics present in Hanford tank waste. It provides the technical basis to address the aspects documented by the Environmental Protection Agency in a draft memo necessary to establish that organics can be immobilized/stabilized in a cementitious matrix. The organic species may interact with cementitious minerals, including slag, fly ash, cement, and other components/dopants like activated carbon. A list of 132 reasonably expected LDR organics associated with Hanford tank waste (RPP-RPT-63493, Rev 1a) was screened by functional groups, octanol-water partitioning coefficients, and detection frequency in Hanford tank waste samples. A subset of 15 compounds spanning the identified properties was tested using sorption and leachate tests. These compounds were subjected to traditional batch sorption tests and leaching tests on/from Cast Stone cementitious material with and without activated carbon (a potential organic adsorption additive). The interactions between organic compounds were observed in traditional batch sorption and modified TCLP leachate tests for surrogate Cast Stone material with and without activated carbon addition. Specifically, phthalic acid (negatively charged) and 4-chloroaniline (polar) sorbed strongly to unmodified Cast Stone, while a group of three phenolic compounds (2,4,6-Trichlorophenol, Pentachlorophenol, o-Cresol) sorbed to and were retained by Cast Stone with a 1% activated carbon amendment. Mixtures of these compounds demonstrated visible retention of the organics within the matrix with a quantifiable affinity of the organics to the cement. This novel testing supports the notion that organics within Hanford tank waste would be retained onto the cement phase with the addition of simple materials easily added to the formulation as part of the treatment process. To further provide additional support to this claim, models employing Flory-Huggins theory in conjunction with macromolecular and organic molecule solubility parameters and species activity coefficients are being utilized to better understand these interactions. Within this context, Van Krevelen (VK) and Hansen organic and polymer solubility parameters were used to inform the Flory-Huggins interaction parameter, allowing for a quantifiable method of evaluating polymer-solvent compatibility and, by extension, interactions with proposed binding additives (e.g., activated carbon). These models and measurement should test the hypothesis for our systems understanding and applicability of Flory-Huggins theory, and whether the prediction of organics in the polymer phase versus the aqueous phase can be confidently made using this strategy. The estimation of the organic molecular activity coefficient in these high salt-saturated aqueous systems is another factor that will be assessed. This project aims to predict and understand complicated organic interactions in these heterogenous systems, guiding optimal treatment and stabilization methodologies.

Jolin, Will [Savannah River National Laboratory (S↗
NASA NTRS2004

The Lewice Console

Lewice (LEWis ICE accretion program) is software used by literally hundreds of users in the aeronautics community for predicting ice shapes, collections efficiencies, and anti-icing heat requirements for aircraft. LEWICE performs its analysis in minutes on a desktop PC, allowing the user to run several parameter studies for design purposes. The ice shape predictions have been used to assess performance degradation both as an input to a CFD program or experimentally in flight or in a wind tunnel. This information is important to ensure an airplane s safe passage through an icing cloud. Currently, Lewice runs as a DOS program that accepts many different inputs such as cloud conditions, wing shapes, and thermal deicing inputs. Usually, such experimental data is stored in spreadsheets. However, Lewice inputs are text files; therefore, they must be generated by the user. Lewice s outputs (collection efficiency, ice shapes and thicknesses) are also text files; to plot the data, users must generate a spreadsheet with this output. Because all Lewice J/O is in the form of text files, using Lewice can be tricky and time-consuming. Our goal was to improve Lewice s usability by creating a user interface that would automatically generate Lewice input from a spreadsheet and automatically put Lewice output into spreadsheets with charts. Additionally, this user interface would automatically convert units (as Lewice only accepts input in certain units) and offer several output options. I call this program the Lewice Console. The Lewice Console is an easy to use interface for Lewice written in Visual Basic. It allows users to run Lewice given a spreadsheet listing experimental conditions. It automatically generates the input to Lewice, does necessary unit conversions, runs Lewice, and produces a spreadsheet with charts plotting the data. It allows users to import data from previously generated Lewice inputs into a spreadsheet. It also allows users to batch run Lewice on several different inputs to automatically generate multiple output spreadsheets. You can also generate plots of actual data vs. experimental data. These capabilities are just the beginning for the Lewice Console. Lewice is capable of running a full deicing experiment given a geometry and heating apparatus information. However, users find it difficult to run such experiments due to the number of inputs and the difficult input file format. The Lewice Console would simplify experiment generation by allowing the user to interactively draw a geometry, place heating apparatus, and specify information about each part. The input to Lewice would be automatically generated from the experiment the user draws on the screen. The Lewice Console would simplify the experiment building process. Currently, Lewice runs as a DOS program that accepts many different inputs such as

Armstrong, Aaron E.↗
NASA NTRS1979

A Study of Linear Approximation Techniques for SAR Azimuth Processing

The application of the step transform subarray processing techniques to synthetic aperture radar (SAR) was studied. The subarray technique permits the application of efficient digital transform computational techniques such as the fast Fourier transform to be applied while offering an effective tool for range migration compensation. Range migration compensation is applied at the subarray level, and with the subarray size based on worst case range migration conditions, a minimum control system is achieved. A baseline processor was designed for a four-look SAR system covering approximately 4096 by 4096 SAR sample field every 2.5 seconds. Implementation of the baseline system was projected using advanced low power technologies. A 20 swath is implemented with approximately 1000 circuits having a power dissipation of from 70 to 195 watts. The baseline batch step transform processor is compared to a continuous strip processor, and variations of the baseline are developed for a wide range of SAR parameters.

Martinson, L. W.↗
DOE OSTI2025

Pyrogenic Organic Matter Laboratory Experiment: Aerobic Respiration and Geochemistry from Variably Inundated Stream Sediments (v3)

This dataset supports a broader study examining the effects of variable inundation and pyrogenic organic matter on ecosystem respiration. The dataset provides data generated from a laboratory batch experiment investigating the interaction between variable inundation conditions (wet and dry sediment) and pyrogenic organic matter (burned and unburned treatments). The contents include time series dissolved oxygen, sediment geochemistry data, and field metadata (including qualitative information on instream and river corridor characteristics). This data package was originally published in November 2025. It was updated in April 2026 (v2; new and modified files) and May 2026 (v3; modified files). See the change history section in the readme for more details For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) international generic sample number (IGSN) mapping file; (5) readme; (6) field protocol; (7) sample name metadata; (8) an environmental context picture for the dry and inundated sampling locations; and (9) a subfolder with sample data from the sediment incubation experiment. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC); (2) total nitrogen (TN); (3) gravimetric moisture; (4) partial pressure and production rates of carbon dioxide, methane, and nitrous oxide; (5) field wet sediment mass, dry sediment mass, water mass, and field wet sediment volume in incubation and sediment NPOC/TN vials; (6) methods codes; (7) respiration rates, pH, and temperature from after the incubation, raw time series dissolved oxygen and temperature, and a subfolder containing associated plots and scripts; (8) ions; (9) FTICR-MS methods; and (10) a subfolder of 12 Tesla (12T) FTICR-MS data. This folder contains the CoreMS processed data and three subfolders, one containing the .xml files, one containing the CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .xml, .html, .Rmd, .py, .cal, .json, or .jpg.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

Comparison of EKF and UKF for Robotic Missions to Mars

The batched Kalman lter works well for a wide variety of orbit determination scenarios but requires linearized approximations to the underlying dynamical system. Furthermore, implementation depends on the computation of partial derivatives for all measurements and parameters used in the lter, a non-trivial task even when analytical derivations are possible. To address these limitations, many extensions and modi - cations to the linear Kalman lter have been proposed over the years since Kalman's initial publication. Our investigation examines two popular variants of the Kalman lter, namely the Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF), within the context of spacecraft missions to Mars. The EKF relies on the same mathematical basis as the classic version of the lter, however it sequentially updates the linearization as the ltering process is conducted. In contrast, the UKF is built upon the unscented transform which samples the nearby solution space and numerically propagates perturbed trajectories as part of the update / estimation process. We assess the performance of the ltering algorithms under a variety of mission scenarios with the goal of establishing decision criteria for which implementation to use under di ering circumstances.

Ely, Todd↗
NASA NTRS1997

Density and composition of microorganisms during long-term (418 day) growth of potato using biologically reclaimed nutrients from inedible plant biomass

A study evaluating alternative methods for long term operation of biomass production systems was recently completed at the Kennedy Space Center (KSC). The 418-day study evaluated repeated batch versus mixed-aged production of potato grown on either standard 1/2-strength Hoagland's nutrient solution or solutions including nutrients recycled from inedible plant material. The long term effects of closure and recycling on microbial dynamics were evaluated by monitoring the microbial communities associated with various habitats within the plant growth system (i.e., plant roots, nutrient solution, biofilms within the hydroponic systems, atmosphere, and atmospheric condensate). Plate count methods were used to enumerate and characterize microorganisms. Microscopic staining methods were used to estunate total cell densities. The primary finding was that the density and composition of microbial communities associated with controlled environmental plant growth systems are stable during long term operation. Continuous production resulted in slightly greater stability. Nutrient recycling, despite the addition of soluble organic material from the waste processing system, did not significantly increase microbial density in any of the habitats.

NASA Center KSC↗
NASA NTRS1990

Attitude determination error analysis - General model and specific application

This paper presents a comprehensive approach to filter and dynamics modeling for attitude determination error analysis. The discussion includes models for both batch least-squares and sequential estimators, a specific dynamic model for attitude determination error analysis of a three-axis stabilized spacecraft equipped with strapdown gyros, and the incorporation of general attitude sensor observations. An analyst using this approach to perform an error analysis chooses a subset of the spacecraft parameters to be 'solve-for' parameters, which are to be estimated, and another subset to be 'consider' parameters, which are assumed to have errors but not to be estimated. The result of the error analysis is an indication of overall uncertainties in the 'solve-for' parameters, as well as the contributions of the various error sources to these uncertainties, including those of errors in the a priori 'solve-for' estimates, of measurement noise, of dynamic noise (also known as process noise or plant noise), and of 'consider' parameter uncertainties. The analysis of attitude, star tracker alignment, and gyro bias uncertainties for the Gamma Ray Observatory spacecraft provide a specific example of the use of a general-purpose software package incorporating these models.

Markley, F. Landis↗
NASA NTRS1993

Comparison of TOPEX/Poseidon orbit determination solutions obtained by the Goddard Space Flight Center Flight Dynamics Division and Precision Orbit Determination Teams

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using the Goddard Trajectory Determination System (GTDS) and a real-time extended Kalman filter estimation system to process Tracking Data and Relay Satellite (TDRS) System (TDRSS) measurements in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. GTDS is the operational orbit determination system used by the FDD, and the extended Kalman fliter was implemented in an analysis prototype system, the Real-Time Orbit Determination System/Enhanced (RTOD/E). The Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generates an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the Geodynamics (GEODYN) orbit determination system with laser ranging tracking data. The TOPEX/Poseidon trajectories were estimated for the October 22 - November 1, 1992, timeframe, for which the latest preliminary POD results were available. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch cases were assessed using overlap comparisons, while the sequential cases were assessed with covariances and the first measurement residuals. The batch least-squares and forward-filtered RTOD/E orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 10 meters (m) for the batch least squares and less than 18 m for the sequential estimation solutions. The differences among the POD, GTDS, and RTOD/E solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C.↗
NASA NTRS2012

SE-FIT

The mathematical theory of capillary surfaces has developed steadily over the centuries, but it was not until the last few decades that new technologies have put a more urgent demand on a substantially more qualitative and quantitative understanding of phenomena relating to capillarity in general. So far, the new theory development successfully predicts the behavior of capillary surfaces for special cases. However, an efficient quantitative mathematical prediction of capillary phenomena related to the shape and stability of geometrically complex equilibrium capillary surfaces remains a significant challenge. As one of many numerical tools, the open-source Surface Evolver (SE) algorithm has played an important role over the last two decades. The current effort was undertaken to provide a front-end to enhance the accessibility of SE for the purposes of design and analysis. Like SE, the new code is open-source and will remain under development for the foreseeable future. The ultimate goal of the current Surface Evolver Fluid Interface Tool (SEFIT) development is to build a fully integrated front-end with a set of graphical user interface (GUI) elements. Such a front-end enables the access to functionalities that are developed along with the GUIs to deal with pre-processing, convergence computation operation, and post-processing. In other words, SE-FIT is not just a GUI front-end, but an integrated environment that can perform sophisticated computational tasks, e.g. importing industry standard file formats and employing parameter sweep functions, which are both lacking in SE, and require minimal interaction by the user. These functions are created using a mixture of Visual Basic and the SE script language. These form the foundation for a high-performance front-end that substantially simplifies use without sacrificing the proven capabilities of SE. The real power of SE-FIT lies in its automated pre-processing, pre-defined geometries, convergence computation operation, computational diagnostic tools, and crash-handling capabilities to sustain extensive computations. SE-FIT performance is enabled by its so-called file-layer mechanism. During the early stages of SE-FIT development, it became necessary to modify the original SE code to enable capabilities required for an enhanced and synchronized communication. To this end, a file-layer was created that serves as a command buffer to ensure a continuous and sequential execution of commands sent from the front-end to SE. It also establishes a proper means for handling crashes. The file layer logs input commands and SE output; it also supports user interruption requests, back and forward operation (i.e. undo and redo), and others. It especially enables the batch mode computation of a series of equilibrium surfaces and the searching of critical parameter values in studying the stability of capillary surfaces. In this way, the modified SE significantly extends the capabilities of the original SE.

Chen, Yongkang↗
DOE OSTI2024

Characterization of Carbon Foam Products Made from Coal at Atmospheric Pressure

Carbon foam products are currently manufactured on a commercial scale utilizing a high-pressure autoclave in batch mode. In an effort to increase manufacturing speed and reduce cost, carbon foam products (both panels and aggregates) were produced at atmospheric pressure in a continuous belt kiln. Here, the microstructure, properties, and performance of these materials were characterized and assessed in various applications. Lightweight aggregate was explored for use as a proppant, as well as it use in concrete and asphalt. The thermal conductivity, water permeability, and fire resistance of panels were evaluated, and its performance with respect to screw pullout was compared against conventional building materials. Some interesting performance aspects and potential for future work was revealed. This work was supported by a cooperative agreement with the Department of Energy – National Energy Technology Laboratory (DE-FOA-0002185, Area of Interest 4: Coal-Derived Carbon Foam Produced via a Continuous Process, award number DE-FE-0031992).

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Performance Comparison of Machine Learning Models for Ultrasonic Nondestructive Evaluation of Alkali-Silica Reaction in Concrete

Alkali-silica reaction (ASR) causes concrete degradation, leading to cracking, rebar corrosion, and reduced structural integrity, which raises safety concerns. Ultrasonic nondestructive evaluation (NDE) effectively assesses concrete properties and monitors ASR progression. However, its deployment and analysis require specialized expertise and subjective interpretation. As computational power increases, artificial intelligence (AI) and machine learning (ML) algorithms are increasingly being used to automate NDE data analysis across various industries for AI-assisted automation. Regulatory agencies are adapting to this technological shift, prompting a need to evaluate current ML technologies’ capabilities and limitations in assessing concrete material properties and damage. This report presents a comparative analysis of four ML regression models for predicting concrete material damage induced by ASR expansion using long-term ultrasonic data monitoring. The models investigated include linear regression (LR), support vector regression (SVR), shallow neural networks (NN), and deep neural networks (DNN). LR, SVR, and shallow NN models use features extracted from ultrasonic signals, whereas the DNN model processes time-domain ultrasonic signals and frequency spectra directly. The study systematically compared the models’ performance from various perspectives, including model input, prediction performance, and generalization ability. The findings indicate significant variability in model performance, with some ML algorithms achieving very high or very low prediction accuracy depending on the preprocessing and feature engineering (extraction and selection) applied. Key insights include the observation that shallow ML models (LR, SVR, and shallow NNs) require meticulous preprocessing and feature extraction to achieve high accuracy. In contrast, the DNN model, although it bypasses the need for feature engineering, necessitates extensive preprocessing to mitigate noise and computational demands. The SVR model emerged as the top performer among the shallow models, and the DNN model exhibited superior performance on specific datasets but struggled with generalization across specimens from different batches. Additionally, the SVR model is sensitive to temperature variations, whereas the DNN model is robust in this regard. Using recurrent neural networks is recommended for future ASR expansion prediction studies. Recurrent neural networks’ inherent ability to capture temporal dependencies and long-term patterns makes them well suited for analyzing sequential ultrasonic monitoring data. Overall, the results and conclusions of this study could provide insights into the capabilities and effectiveness of ML when applied to ultrasonic NDE data and help identify best practices for using ML for ultrasonic NDE of concrete material properties.

36 MATERIALS SCIENCE↗
NASA NTRS1994

High-efficiency, deep-junction, epitaxial InP solar cells on (100) and (111)B InP substrates

We report on the development and performance of deep-junction (approximately 0.25 micron), graded-emitter-doped, n(sup +)-p InP solar cells grown by metallorganic chemical vapor deposition (MOCVD). A novel, diffusion-transport process for obtaining lightly-doped p-type base regions of the solar cell is described. The I-V data and external quantum-efficiency response of these cells are presented. The best active-area AMO efficiency for these deep-junction cells on (100)-oriented InP substrates is 16.8 percent, with a J(sub SC) of 31.8 mA/sq cm, a V(sub OC) of 0.843 V, and a fill-factor of 0.85. By comparison, the best cell efficiency on the (111)B-oriented InP substrates was 15.0 percent. These efficiency values for deep-junction cells are encouraging and compare favorably with performance of thin-emitter (0.03 micron) epitaxial cells as well as that of deep-emitter diffused cells. The cell performance and breakdown voltage characteristics of a batch of 20 cells on each of the orientations are presented, indicating the superior breakdown voltage properties and other characteristics of InP cells on the (111)B orientation. Spectral response, dark I-V data, and photoluminescence (PL) measurements on the InP cells are presented with an analysis on the variation in J(sub SC) and V(sub OC) of the cells. It is observed, under open-circuit conditions, that lower-V(sub OC) cells exhibit higher band-edge PL intensity for both the (100) and (111)B orientations. This anomalous behavior suggests that radiative recombination in the heavily-doped n(sup +)-InP emitter may be detrimental to achieving higher V(sub OC) in n(sup +)-p InP solar cells.

Venkatasubramanian, R.↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Carbon and oxygen isotope study of carbonates from highly shocked clasts of the polymict breccia of the Haughton Crater (Canada)

It is known that the release of volatiles on impact is an important controlling factor in cratering processes in carbonate terranes and in the mobility of chemical elements. In order to assess the nature and the role of carbon- and oxygen-bearing volatiles during impact-induced metamorphism of sedimentary rocks, the C-13/C-12 and O-18/O-16 ratios and carbonate contents were determined for 30 shocked clasts from the Haughton Crater polymict breccia as well as for some unshocked carbonates from the sedimentary cover adjacent to the crater. Shock-induced CO2 loss during decarbonation of calcite is known to be a function of peak pressure and ambient partial pressure of the volatile species. In our clast samples, shocked from 20 to 60 GPa, we expect about 20 to 100 percent CO2 loss and preferential depletion in C-13 and O-18 in the residual carbonate. Rayleigh model (progressive loss of CO2) and batch model (single-step loss of CO2) curves for this depletion are shown. The magnitudes of the C-13 and O-18 depletions increase with the increase of the CO2 loss. In addition, the isotopic depletions should be correlated with an enrichment in CaO and MgO in the residual solid.

Agrinier, P.↗
NASA NTRS2007

Tracking and Data Relay Satellite (TDRS) Orbit Estimation Using an Extended Kalman Filter

Alternatives to the Tracking and Data Relay Satellite (TDRS) orbit estimation procedure were studied to develop a technique that both produces more reliable results and is more amenable to automation than the prior procedure. The Earth Observing System (EOS) Terra mission has TDRS ephemeris prediction 3(sigma) requirements of 75 meters in position and 5.5 millimeters per second in velocity over a 1.5-day prediction span. Meeting these requirements sometimes required reruns of the prior orbit determination (OD) process, with manual editing of tracking data to get an acceptable solution. After a study of the available alternatives, the Flight Dynamics Facility (FDF) began using the Real-Time Orbit Determination (RTOD(Registered TradeMark)) Kalman filter program for operational support of TDRSs in February 2007. This extended Kalman filter (EKF) is used for daily support, including within hours after most thrusting, to estimate the spacecraft position, velocity, and solar radiation coefficient of reflectivity (C(sub R)). The tracking data used are from the Bilateration Ranging Transponder System (BRTS), selected TDRS System (TDRSS) User satellite tracking data, and Telemetry, Tracking, and Command (TT&C) data. Degraded filter results right after maneuvers and some momentum unloads provided incentive for a hybrid OD technique. The results of combining EKF strengths with the Goddard Trajectory Determination System (GTDS) Differential Correction (DC) program batch-least-squares solutions, as recommended in a 2005 paper on the chain-bias technique, are also presented.

Ward, Douglas T.↗
DOE OSTI2024

Core Design Optimization of the Westinghouse Lead Fast Reactor

Westinghouse is pursuing an advanced Nuclear Power Plant design based on Lead Fast Reactor (LFR) technology for global commercialization. To achieve an optimal combination of key attributes, such as safety, sustainability, and economic competitiveness, Westinghouse and ANL partnered in developing and applying a formalized core design optimization strategy. An LFR analysis workflow was developed to automate a suite of reactor physics, fuels performance, safety, and economics simulations on a selected LFR concept. The workflow streamlines analysis of a wide range of LFR designs with different dimensions and fuel types to assess their viability and economic performance, significantly reducing human processing time and risks of processing errors. The LFR optimization exercise was defined, resulting in selection of the design constraints (geometric, neutronics, thermo-mechanical, safety, thermal-hydraulics, and economics) and performance metrics researched (minimization of both the fuels LCOE and the first core inventory cost). A total of 14 varied design parameters were considered, including assembly dimensions, coolant temperature, and enrichment distribution throughout the core. The LFR analysis workflow was connected to DAKOTA for sensitivity and optimization analyses. Due to the extremely large size of the potential LFR optimization solution space relative to the computing time required to characterize one LFR solution, a multi-stage optimization approach was proposed to breakdown the problem into several stages with more reasonable sizes. This optimization approach enabled finding various viable core solutions with different cost tradeoffs that were considered by Westinghouse and justify selection of a smaller core with multi-batch 2-year cycle length.

Stauff, Nicolas E.↗