Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STEP requirements”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
DOE OSTI2025

Neutrino oscillations

Neutrino oscillation measurements provide an important window on what lies beyond the Standard Model of particle physics. These measurements may unlock the mechanism by which the universe came to become matter-dominated, and may offer hints of another mechanism to generate particles with mass. Measurements of oscillations from muon neutrinos to electron neutrinos (and their antineutrino counterpart) as a function of time provide critical inputs to understanding both mechanisms. These are challenging measurements and a variety of techniques and strategies are required to get the complete picture. This article describes the current status of our understanding of neutrino masses and how neutrinos oscillate between one flavour and another as they propagate through space and time, and what remains to be understood. Taking the next steps in this field requires a variety of approaches and a better understanding of how neutrinos interact with nuclei. This article describes two of those next steps, highlighting where Canadian groups are active, and concludes with a discussion of the broad range of additional physics that becomes accessible by having two very different large sensitive neutrino detectors making these measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Alternative approach to quantum imaginary time evolution

There is increasing interest in quantum algorithms (QAs) that are based on the imaginary time evolution (ITE), a successful classical numerical approach to obtain ground states. However, most of the proposals so far require heavy postprocessing computational steps on a classical computer, such as solving linear equations. Here we provide an alternative approach to implement ITE. A key feature in our approach is the use of an orthogonal basis set: the propagated state is efficiently expressed in terms of orthogonal basis states at every step of the evolution. We argue that the number of basis states needed at those steps to achieve an accurate solution can be kept on the order of n , the number of qubits, by controlling the precision (number of significant digits) and the imaginary time increment. The number of quantum gates per imaginary time step is estimated to be polynomial in n . Additionally, while in many QAs the locality of the Hamiltonian is a key assumption, in our algorithm this restriction is not required. This characteristic of our algorithm renders it useful for studying highly nonlocal systems, such as the occupation-representation nuclear shell model. Here, we illustrate our algorithm through numerical implementation on an IBM quantum simulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Demonstration and analysis of volumetric additive manufacturing via sub-orbital spaceflight testing

Computed Axial Lithography (CAL) represents a significant advancement in the emerging field of Volumetric Additive Manufacturing (VAM). CAL addresses key limitations of traditional photopolymer additive manufacturing technologies, by eliminating the need for layering and support structures. Unlike conventional methods, CAL prints components by illuminating all points within a desired geometry simultaneously, using tomographic reconstruction to form the object in a single step. This unique approach eliminates the relative motion between the object and the precursor material, enabling faster printing speeds and reducing the waste associated with support structures. However, CAL parts require post-processing steps before they can be utilized. CAL's core attributes make it particularly suited for In-Space Manufacturing (ISM), due to its fast fabrication times, wide breadth of materials it can use, and minimized footprint. CAL has been successfully demonstrated in microgravity during parabolic flight experiments. However to fully validate and understand CAL's behaviour in microgravity, all manufacturing and post-processing steps must be integrated. In June 2024, we conducted SpaceCAL Mission 3, testing this entire workflow on a suborbital flight aboard Virgin Galactic's SpaceShipTwo. During ~140 s of microgravity, the system autonomously manufactured and post-processed four parts using PEGDA700 resin. Post-flight analysis showed that 2/4 parts were recognisable, while others were distorted due to bubble formation from residual water droplets, off-axis optical aberrations, and non-uniform solvent rinsing. Despite these limitations, this study represents the first integrated CAL workflow in space, providing an initial experimental demonstration and analysis for closed-loop in-space manufacturing.

Waddell, Taylor [Department of Mechanical Engineer↗
DOE OSTI2020

Testing of Crystalline Silicotitanate to Support Tank-Side Cesium Removal System Operations - 20458

The direct feed of waste from the Hanford tank farms to the Low-Activity Waste Facility at the Hanford Waste Treatment and Immobilization Plant requires an intermediate treatment step that prepares the waste for vitrification. Washington River Protection Solutions (WRPS) selected a near-tank treatment system known as the Tank-Side Cesium Removal (TSCR) system to perform the required intermediate treatment functions. The approach in the TSCR system is to filter waste supernatant with a sintered metal dead-end filter and then use a series of ion exchange columns to remove cesium; this system has several similarities to the Tank Closure Cesium Removal system that has been deployed on the Savannah River Site. The ion exchange media proposed for use is crystalline silicotitanate (CST), which is a non-elutable media with a high affinity for cesium. Because the media is non-elutable, loaded columns will be blown down with compressed air, moved into interim storage, and replaced with new columns when cesium capacity is reached. The assumed extent and rate of drying that could be achieved in the TSCR column geometry lacked confirmation by experimental data. In addition, once the columns are loaded with cesium, they generate flammable gases (primarily hydrogen) via radiolysis acting on any resident liquids, i.e., moisture, remaining in the media bed. An assessment of available information on CST concluded that gas generation data that bounds expected operational conditions were needed to support the TSCR system safety basis and planned operations. Working with WRPS, Pacific Northwest National Laboratory (PNNL) designed and conducted testing with CST media to address the need for (a) representative in-column drying data, and (b) bounding flammable gas generation data. Drying testing was performed using a full-height column with a diameter of approximately 2 inches with injected air at temperatures of 18 deg. C and 30 deg. C. Gas generation testing was conducted in an engineered bunker that allowed simultaneous irradiation of up to eight samples at a time. These two experimental approaches used at PNNL for the CST testing are summarized and the major outcomes are presented. (authors)

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

DOE JGI Metagenome Workflow

The DOE Joint Genome Institute (JGI) Metagenome Workflow performs metagenome data processing, including assembly; structural, functional, and taxonomic annotation; and binning of metagenomic data sets that are subsequently included into the Integrated Microbial Genomes and Microbiomes (IMG/M) (I.-M. A. Chen, K. Chu, K. Palaniappan, A. Ratner, et al., Nucleic Acids Res, 49:D751–D763, 2021, https://doi.org/10.1093/nar/gkaa939) comparative analysis system and provided for download via the JGI data portal (https://genome.jgi.doe.gov/portal/). This workflow scales to run on thousands of metagenome samples per year, which can vary by the complexity of microbial communities and sequencing depth. Here, we describe the different tools, databases, and parameters used at different steps of the workflow to help with the interpretation of metagenome data available in IMG and to enable researchers to apply this workflow to their own data. We use 20 publicly available sediment metagenomes to illustrate the computing requirements for the different steps and highlight the typical results of data processing. The workflow modules for read filtering and metagenome assembly are available as a workflow description language (WDL) file (https://code.jgi.doe.gov/BFoster/jgi_meta_wdl). The workflow modules for annotation and binning are provided as a service to the user community at https://img.jgi.doe.gov/submit and require filling out the project and associated metadata descriptions in the Genomes OnLine Database (GOLD) (S. Mukherjee, D. Stamatis, J. Bertsch, G. Ovchinnikova, et al., Nucleic Acids Res, 49:D723–D733, 2021, https://doi.org/10.1093/nar/gkaa983).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Structure-function studies reveal ComEA contains an oligomerization domain essential for transformation in gram-positive bacteria

An essential step in bacterial transformation is the uptake of DNA into the periplasm, across the thick peptidoglycan cell wall of Gram-positive bacteria, or the outer membrane and thin peptidoglycan layer of Gram-negative bacteria. ComEA, a DNA-binding protein widely conserved in transformable bacteria, is required for this uptake step. Here we determine X-ray crystal structures of ComEA from two Gram-positive species, Bacillus subtilis and Geobacillus stearothermophilus, identifying a domain that is absent in Gram-negative bacteria. X-ray crystallographic, genetic, and analytical ultracentrifugation (AUC) analyses reveal that this domain drives ComEA oligomerization, which we show is required for transformation. We use multi-wavelength AUC (MW-AUC) to characterize the interaction between DNA and the ComEA DNA-binding domain. Finally, we present a model for the interaction of the ComEA DNA-binding domain with DNA, suggesting that ComEA oligomerization may provide a pulling force that drives DNA uptake across the thick cell walls of Gram-positive bacteria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Gridded Sub-daily Climate Forcings for North America Based on Daymet and ERA5 (Daymet-ERA5)

To support high spatial and temporal resolution land surface modeling, this dataset provides hourly time step historic weather forcing at 1-km spatial resolution for the entire North America. The latest Daymet V4 data provides gridded historic daily weather observations at 1-km spatial resolution from 1980 to 2024. Using sub-daily temporal information from the ECMWF ERA5 reanalysis, Daymet was further temporally downscaled to hourly time steps and provided in the format required for land surface model simulations. The process of temporal downscaling preserves the relative magnitude in each hourly time step from ERA5 while maintaining the total and average values from Daymet for each day. This results in a blended 1980-2024 Daymet-ERA5 dataset. Available variables include surface air temperature, precipitation, specific humidity, shortwave and longwave radiation, wind speed, and pressure. These data can be used as a high-resolution meteorological forcing dataset to support high-resolution land surface modeling where accurate meteorological forcing datasets built from historic observations and/or reanalysis datasets are desirable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Pressure stability in explicitly coupled simulations of poromechanics with application to CO 2 sequestration

We study in detail the pressure stabilizing effects of the non-iterated fixed-stress splitting in poromechanical problems which are nearly undrained and incompressible. When applied in conjunction with a spatial discretization which does not satisfy the discrete inf–sup condition, namely a mixed piecewise linear–piecewise constant spatial discretization, the explicit fixed-stress scheme can have a pressure stabilizing effect in transient problems. This effect disappears, however, upon time step refinement or the attainment of steady state. The interpretation of the scheme as an Augmented Lagrangian method similar to Uzawa iteration for incompressible flow helps explain these results. Moreover, due to the slowly evolving solution within undrained seal regions, we show that the explicit fixed-stress scheme requires very large time steps to reveal its pressure stabilizing effect in examples of geologic CO 2 sequestration. We note that large time steps can result in large errors in drained regions, such as the aquifer or reservoir regions of these examples, and can prevent convergence of nonlinear solvers in the case of multiphase flows, which can make the explicit scheme an unreliable source of pressure stabilization. We conclude by demonstrating that pressure jump stabilization is as effective in the explicit fixed-stress setting as in the fully implicit setting for undrained problems, while maintaining the stability and convergence of the fixed-stress split for drained problems.

58 GEOSCIENCES↗
DOE OSTI2020

Coal direct chemical looping process: 250 kW pilot-scale testing for power generation and carbon capture

Chemical looping combustion (CLC) is an energy conversion technology that can produce concentrated CO 2 stream without the need for a gas separation step, and thus, has the potential to drastically reduce the energy consumption and cost associated with CO 2 capture in power generation. The coal-direct chemical looping (CDCL) process is a CLC technology that uses a moving bed reducer configuration that can directly consume coal as a feedstock without requiring an upstream gasification step. An integrated 250 kW th CDCL pilot unit using iron-based oxygen carriers was constructed and demonstrated for over 1000 h of testing. The principles for the CDCL pilot unit design and operation are summarized in this article. During the 288-hour continuous operation testing, the CDCL pilot unit achieved >96% coal conversion with a CO 2 purity of >97%. Low carbon carryover into the combustor, i.e. <2%, was also confirmed during the test, which shows the capability of the moving bed reactor to retain and convert coal using the oxygen available on the iron-based oxygen carrier. The results from the pilot unit testing confirms the CDCL concept as a promising coal combustion technology for heat and power generation with CO 2 capture.

10 SYNTHETIC FUELS↗
DOE OSTI2023

Polymer-based microfluidic device for on-chip counter-diffusive crystallization and in situ X-ray crystallography at room temperature

Proteins are long chains of amino acid residues that perform a myriad of functions in living organisms, including enzymatic reactions, signalling, and maintaining structural integrity. Protein function is determined directly by the protein structure. X-ray crystallography is the primary technique for determining the 3D structure of proteins, and facilitates understanding the effects of protein structure on function. The first step towards structure determination is crystallizing the protein of interest. We have developed a centrifugally-actuated microfluidic device that incorporates the fluid handling and metering necessary for protein crystallization. Liquid handling takes advantage of surface forces to control fluid flow and enable metering, without the need for any fluidic or pump connections. Our approach requires only the simple steps of pipetting the crystallization reagents into the device followed by either spinning or shaking to set up counter-diffusive protein crystallization trials. The use of thin, UV-curable polymers with a high level of X-ray transparency allows for in situ X-ray crystallography, eliminating the manual handling of fragile protein crystals and streamlining the process of protein structure analysis. We demonstrate the utility of our device using hen egg white lysozyme as a model system, followed by the crystallization and in situ, room temperature structural analysis of the hub domain of calcium–calmodulin dependent kinase II (CaMKIIβ).

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

Low synchronization Gram–Schmidt and generalized minimal residual algorithms

The Gram–Schmidt process uses orthogonal projection to construct the A = QR factorization of a matrix. When Q has linearly independent columns, the operator P = I - Q(QTQ)-1QT defines an orthogonal projection onto Q⊥. In finite precision, Q loses orthogonality as the factorization progresses. A family of approximate projections is derived with the form P = I - QTQT, with correction matrix T. When T = (QTQ)-1, and T is triangular, it is postulated that the best achievable orthogonality is $\mathcal{O}(ε)\mathcal{K}(A)$. We present new variants of modified (MGS) and classical Gram–Schmidt algorithms that require one global reduction step. An interesting form of the projector leads to a compact WY representation for MGS. In particular, the inverse compact WY MGS algorithm is equivalent to a lower triangular solve. Our main contribution is to introduce a backward normalization lag into the compact WY representation, resulting in a $\mathcal{O}(ε)\mathcal{K}[r_0, AV_m])$ stable Generalized Minimal Residual Method (GMRES) algorithm that requires only one global reduce per iteration. Finally, further improvements in performance are achieved by accelerating GMRES on GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Conceptual Basis and Techno-Economic Modeling for Integrated Algal Biorefinery Conversion of Microalgae to Fuels and Products (2019 NREL TEA Update: Highlighting Paths to Future Cost Goals via a New Pathway for Combined Algal Processing)

The report documents the conceptual basis for a new potential Combined Algal Processing design strategy which may allow more flexibility in accommodating different algal biomass feedstock compositions, by enabling upgrading of both protein and carbohydrates in a single step, without a strict requirement for either component to be in soluble or monomeric form, while maintaining effective wet lipid extraction techniques to enable high lipid recoveries. In light of previously-established constraints around algal biomass costs (which are significantly higher than lignocellulosic terrestrial biomass), the present CAP processing strategy reflects an integrated biorefinery concept producing both fuels and value-added chemical coproducts as a means to improve profitability and generate coproduct revenues to help drive down the minimum fuel selling price (MFSP) towards economically viable levels. Namely, this report highlights an integrated CAP biorefinery process and associated technical targets that would be required to achieve U.S. Department of Energy target MFSP goals of $2.5/gallon gasoline equivalent by 2030. This is accomplished by a process involving low-cost seasonal storage of algal biomass during high-growth seasons, rapid flash hydrolysis pretreatment of the biomass, solvent extraction of pretreated biomass, cleanup and fractionation of lipids into triglyceride and free fatty acid fractions, and a series of thermochemical conversion steps to upgrade carbohydrates and protein to hydrocarbon fuels. These steps include mild oxidative treatment (MOT), a process originally investigated at NREL for upgrading lignin, followed by catalytic ketonization and hydrotreating of MOT products to fuels. Isolated triglycerides are sent to a coproduct train, with the base case focused on upgrading to polyurethane foams as a high-value, high-market-volume coproduct.

09 BIOMASS FUELS↗
DOE OSTI2025

Automation of Vulnerability and Patch Management: Information Extraction, Association, and Optimization

Vulnerability and patch management is an integral part of a robust cybersecurity program, yet it grows increasingly complex due to the sheer amount of data that must be analyzed. Particularly in Operational Technology (OT) environments, analysis must be done manually because of the lack of automated solutions. Additionally, there are many steps in this process, from the initial discovery of the vulnerability to the implementation of its remediation, and each step in the process requires different data in order to be performed effectively. In this work, we provide approaches and strategies to assist operators in industrial or OT environments throughout the vulnerability management cycle. Security advisories provide key information about mitigation strategies, or actions that can be taken when a patch is unavailable or cannot be installed. Details of these strategies are not shared in public vulnerability databases and must be found manually. We approach this problem by designing a solution to automatically identify that information within vendor security advisories and retrieve it for operator use. We start with an approach that requires domain-specific knowledge of certain frequently-seen reference websites. Next, an approach that can work on an arbitrary website but relies on certain keywords. Finally, an approach that uses Natural Language Processing (NLP) methods and does not require specific knowledge or keywords. Each of these approaches is more general than its predecessor; we demonstrate high accuracy for all approaches Advisories also often contain details of affected products in non-standard or natural language formats. While this information can be easily understood when read by an operator, the non-standard format acts as a barrier to effective automation. We provide an approach for the first step in this process: identifying vendors in security advisories and mapping them to a standard framework for representing digital assets and software products. We evaluate five established string similarity algorithms, plus one of our own design that combines string similarity and information theory, on the task of mapping vendors to their corresponding entries in the Common Platform Enumeration (CPE) repository. Our results show that our proposed metric outperforms all others. Due to the constraints on time, finances, and personnel for organizations, Large Language Models (LLMs) may seem like attractive opportunities for security operators to speed up information gathering; however, it is still not clear whether LLMs can handle vulnerability management tasks well. To answer this question, we perform an empirical study of LLMs’ ability to provide consistent, accurate information about vulnerabilities in order to guide organizations in their adoption of LLMs. We observe poor performance for all models tested, suggesting that these models are not well-suited to the consistent retrieval of accurate vulnerability information. Finally, once vulnerabilities have been identified and any additional information has been obtained, operators must decide which remediation actions to implement based on their available resources. This already-complex problem becomes even more so when we consider that a vulnerability may have multiple avenues for remediation. We formulate this scenario as two knapsack problems and provide solutions, which we then compare against several existing strategies for vulnerability prioritization seen in real operational environments.

McClanahan, Kylie↗
DOE OSTI2022

Oak Ridge National Laboratory FY 2023 Site Sustainability Plan With FY 2022 Performance Data

At the close of each fiscal year, the US Department of Energy (DOE) Sustainability Performance Division (SPD) issues guidance documents and technical resource aids/tools necessary for DOE sites and national laboratories to complete sustainability reporting requirements. SPD is part of the DOE Office of Asset Management. As required by DOE Order 436.1, Departmental Sustainability, “each site will develop and commit to an annual Site Sustainability Plan (SSP) that identifies its respective contribution toward meeting the DOE’s sustainability goals.” SPD collects and compiles information reported by each site to develop an agency-wide Sustainability Report and Implementation Plan, which is used to report DOE sustainability progress to the federal government as required by all major federal agencies. DOE launched a formal Sustainability Office and annual SSP process in 2011. Each year, Oak Ridge National Laboratory (ORNL), in concert with the Office of Science (SC), provides the resources essential to fulfill its commitment to deliver a complete and accurate SSP report and quality performance data for entry into the DOE Sustainability Dashboard as managed by SPD. The performance data entered by each DOE site are then combined to disclose the progress of each DOE Program Office and are further combined to show comprehensive progress for the agency. The Office of Asset Management provides assistance to program offices in sustaining their missions, freeing up resources by reducing waste, avoiding excess expenditure on utilities, maximizing productivity, and improving the efficiency of facilities and processes. By focusing on mission needs, programs and associated DOE sites can help the agency meet its sustainability goals, as outlined in federal statutory and regulatory requirements. In FY 2022, the SSP guidance was updated to capture requirements from Executive Order (EO) 14008, Tackling the Climate Crisis at Home and Abroad, the Energy Act of 2020 (EAct 20), actions outlined in DOE’s Climate Adaptation & Resilience Plan and Sustainability Plan, and EO 14057, Catalyzing Clean Energy Industries and Jobs Through Federal Sustainability. Updates in SSP guidance help to minimize and streamline reporting while simultaneously addressing updated federal requirements. Per DOE, each SSP report should provide an overview of the site’s planned actions, as well as an overview of efforts and accomplishments during the reporting period. SPD collects and compiles information reported by each site to develop DOE’s Annual Sustainability Report, Climate Adaptation & Resilience Plan, and Annual Energy Management Report to Congress. The agency goal has been to lower the reporting burden for sites and increase and improve the consistency of information available to decision makers, allowing them to better identify projects and potential for increased efficiency, as well as to reduce waste, lower emissions, and enhance operational resilience. Sites may elect to produce a more polished publication for their leadership and stakeholders, but this step is no longer required. The ORNL SSP narrative report (this document) and the reporting of DOE SPD Sustainability Dashboard performance data is a collaborative effort of approximately 30 subject matter experts (SMEs) from ORNL facility management and research divisions. Annually, these associates come together to provide a report that can be used by DOE to demonstrate continued agency progress in energy efficiency and sustainable federal operations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm

Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales $\textit{O}$ (10 0–6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions $\textit{O}$ (10 0–2 years) precisely within longer $\textit{O}$ (10 3–6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution $\textit{O}$ (10 0–2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution $\textit{O}$ (10 3–6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ~ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Toward designing effective exascale scientific computing workflows: experiences and best practices

Many fields within scientific computing have embraced advances in big-data analysis and machine learning, which often requires the deployment of large, distributed and complicated workflows that may combine training neural networks, performing simulations, running inference, and performing database queries and data analysis in asynchronous, parallel and pipelined execution frameworks. Such a shift has brought into focus the need for scalable, efficient workflow management solutions with reproducibility, error and provenance handling, traceability, and checkpoint-restart capabilities, among other needs. Here, we discuss challenges and best-practices for deploying exascale-generation computational science workflows on resources at the Oak Ridge Leadership Computing Facility (OLCF). We present our experiences with large-scale deployment of distributed workflows on the Summit supercomputer, including for bioinformatics and computational biophysics, materials science, and deep learning model optimization. We also present problems and solutions created by working within a Python-centric software base on traditional HPC systems, and discuss steps that will be required before the convergence of HPC, AI, and data science can be fully realized. Our results point to a wealth of exciting new possibilities for harnessing this convergence to tackle new scientific challenges.

Coletti, Mark↗
DOE OSTI2023

Evaluating Retention Index Score Assumptions to Refine GC–MS Metabolite Identification

As metabolomics grows into a high-throughput and high demand research field, current metrics for the identification of small molecules in gas chromatography–mass spectrometry (GC–MS) still require manual verification. Though steps have been taken to improve scoring metrics by combining spectral similarity (SS) and retention index (RI), the problem persists. A large body of literature has analyzed and refined SS scores, but few studies have explicitly studied improvements to RI scores. Here, in this study, we examined whether uninvestigated assumptions of the RI score are valid and propose ways to improve them. Query RIs were matched to library RI with a generous window of ±35 to avoid unintentional removal of valid compound identifications. Each match was manually verified as a true positive (TP), true negative, or unknown. Metabolites with at least 30 TP identifications were included in downstream analyses, resulting in a total of 87 metabolites from samples of varying complexity and type (e.g., amino acid mixtures, human urine, fungal species, and so on.). Our results showed that the RI score assumptions of normality, consistent variance across metabolites, and a mean error centered at 0 are often violated. We demonstrated through a cross-validation analysis that modifying these underlying assumptions according to empirical metabolite-specific distributions improved the TP and negative rankings. Further, we statistically determined the minimum number of samples required to estimate distributional parameters for scoring metrics. Overall, this work proposes a robust statistical pipeline to reduce the time bottleneck of metabolite identification by improving RI scores and thus minimize the effort to complete manual verification.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗