Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sequential ranging”

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.

197 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Extreme-scale stochastic optimization and simulation via learning-enhanced decomposition and parallelization (Final Technical Report)

Stochastic optimization and simulation models ubiquitously arise in designing and operating complex service/engineering systems. They can be extreme in scale due to high-dimensional data and decisions, and can also involve decisions made sequentially in response to newly revealed data, both causing significant computational challenge. The objective of this research is to explore a unified framework that integrates machine learning with discrete optimization and risk-averse modeling, to improve the efficiency of decomposition paradigms for stochastic optimization and simulations at extreme scale. The models we consider represent a broad class of complex decision-making problems, where 0-1 or continuous decisions are made before and/or after knowing multiple sources of uncertainties that could be correlated. We will employ machine learning methods to dynamically decide and prioritize computational procedures, including cut generation, branching, and bounding of the optimal objective. Furthermore, the research will shed new lights on the traditional decomposition algorithms for extreme-scale computing. Deliverables of the research include new modeling and computational methods for advancing the state-of-the-art research in optimization and simulation, bringing many relevant risk-averse, data-driven optimization problems in practice within the range of tractability. Examples include distributed computing server scheduling and sensor deployment for monitoring critical infrastructures. Success in this effort will enable progress in solving multiple extreme-scale problems in the complex system design and operations arising from DoE missions in energy, environment, and national security.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Autonomous Hydrogen Fueling Station

This project “Autonomous Hydrogen Fueling Station” covered the autonomous refueling with both gaseous hydrogen and liquid hydrogen. The part on gaseous hydrogen focused on the development of an autonomous robotic fueling arm that would couple to a fuel cell engine for hydrogen refueling without guidance from the forklift operator and budget period. Research was also covered for the robotic fueling with a commercial vehicle. The second phase of the project created the baseline for an autonomous liquid hydrogen transfer system that would minimize boil off losses by operating at thermodynamically efficient state points. For the development of the robotic fueling arm, testing was conducted to establish a baseline measurement of the accuracy and repeatability of a human operator positioning a lift truck in front of a dispenser. The goal was to establish the range of motion required for an autonomous fueling mechanism to mate a hydrogen nozzle with a receptacle on a fuel cell system installed in a forklift. The final design comprised a selective compliance articulated robot arm (SCARA)-type mechanism with two arms for horizontal motion and a ball screw for vertical movement and color and LIDAR cameras were used for marker identification and proximity awareness to guide the robotic arm to its target receptacle. Initial tests resulted in 199 out of 200 successful attempts at autonomous coupling of the dispensing coupler and a fuel cell engine, without hydrogen. The dispenser prototype was modified to include tubing for both hydrogen fuel and air purge lines, but subsequent tests were confounded by the shoulder motor over current errors which limited the robot from getting to the fully inserted position to achieve a positive latch. Robotic hydrogen refueling was successfully demonstrated over 1.5 hours of testing, Plug completed 29 successful latches with an average number of 4 sequential latches before failure. However, a robot capable of placing the nozzle with more force is required for higher reliability. For budget period two, a small scale (10 kg / transfer) automated control system was designed that would operate valves to control pressure and flow of liquid nitrogen between a source and receiving tank with an aim to minimize boil off losses by operating at the most thermodynamically efficient state points. Control system logic flow and a P&ID were developed prior to system safety characterization via HAZOP. A control narrative and system state points were defined. Delays in approval for a change of project objective and procurement issues precluded the construction and test of the final prototype system.

08 HYDROGEN↗
DOE OSTI2026

Assessing Structural, Thermal, and Functional Characteristics of Marigold Flower Protein as a Sustainable Food Ingredient

The demand for sustainable and alternative protein sources has been on the rise, driving interest in the valorization of underutilized plants. This study evaluated Calendula officinalis (marigold), a common floral waste, as a sustainable alternative protein source for the food industry. The primary objective of this study was to investigate the physicochemical properties of protein fractions from Calendula officinalis flower to evaluate their potential as a novel protein ingredient. Extraction of the Calendula officinalis flower yielded 92.17% of the crude protein. A sequential extraction of albumin, globulin, glutelin, and prolamin from marigold flower revealed albumin as the dominant fraction (65.47%) and exhibited the highest protein functionality, including water-holding capacity (2.37 g/g), oil-holding capacity (2.49 g/g), and emulsifying capacity (65.22 mL/g). Compared with other protein fractions, glutelin showed a relatively high emulsifying and foaming capacity (EC: 59.13 mL/g; FC: 16.23%). Differential scanning calorimetry revealed high thermal stability for albumin (T p = 105.28 °C) and glutelin (T p = 97.6 °C). Sodium Dodecyl Sulfate–Polyacrylamide Gel Electrophoresis (SDS-PAGE) and Liquid Chromatography–Mass Spectrometry (LC-MS) confirmed the presence of abundant low-molecular-weight polypeptides (<37 kDa), which enhanced emulsification, while scanning electron microscopy revealed porous structures aligned with hydration properties. Antioxidant activity was higher in albumin and glutelin, linked to surface hydrophobicity. LC-MS/MS identified 33 short-chain proteins, including oxidoreductase proteins and lipid-transfer proteins. Findings highlight marigold flower proteins as a sustainable, functional ingredient for a diverse range of food applications.

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

Development of NDE/NDT Tools for High-Volume & High-Speed Inspection of CFRP Structures in Automotive Manufacturing

Main advantages of the air-coupled ultrasound testing (ACUT) and electromagnetic testing (EMT) techniques for NDE of CFRP composites were non-contact sensing, scalability for high-speed inspection, cost-effectiveness, and non-hazardous operation. Despite these advantages, no systems that would satisfy the project requirements were commercially available. Hence, one of the major efforts of the Michigan State University (MSU) team at the initial stage of the project was to close this technological gap by developing, optimizing, and validating array sensors that would provide sufficient sensitivity, spatial coverage, and resolution for robust defect detection. Optimization of the ACUT and EMT sensor designs was performed using experimentally validated finite element models. Initial experiments using array probes were conducted on relatively flat CFRP samples. In parallel, the MSU team designed and assembled a portable platform with two robotic arms. The robots were equipped with newly designed sensors that enabled high-speed NDE of curved CFRP parts. Presently, the developed robotic platform can be used as a demo/template NDE system, which is easily adaptable to manufacturing environments and in-line NDE. The ACUT NDE system developed by the MSU team used a high-power 4-channel pulser receiver for parallel data acquisition. The array probes were designed by stacking commercially available ACUT transducers, which operated in the frequency range between 100 kHz and 500 kHz. MSU optimized the excitation procedure and developed wave focusing cones so as to reduce the crosstalk between the transducers and to provide higher pulse repletion frequency (PRF). The through-transmission (TT) and single-side access (SSA) inspection modes were successfully implemented. In the TT-ACUT, structural defects in CFRP were detected by passing ultrasonic waves through the test part. Hence, the ACUT transmitters and receivers needed to be placed on the opposite sides of the test part. In the SSA-ACUT, guided waves (GW) were excited in the test part using the transmitters and were sensed by the receivers from the same side. Multi-channel TT-ACUT and SSA-ACUT provided high-speed NDE, and were successfully validated on CFRP test samples with interlaminar delaminations and other embedded defects The EM techniques developed by the MSU team included: 1) eddy current testing (ECT), 2) capacitive imaging (CI) and hybrid dual-mode imaging. In ECT, structural damage was detected in CFRP using coils sensor arrays. In ECT, the excitation magnetic field is generated by passing an alternating current through a coil, which is placed above the test sample. The excitation field penetrates the conductive sample and induces the eddy currents in its transect. In turn, the eddy currents generate the reaction field, which affects the total field sensed by a coil. Hence, the presence of structural flaws will alter the eddy current flow and the picked-up signal. ECT is mostly sensitive to local changes of the electric conductivity of the test sample, and CFRPs are mostly conductive in the direction of carbon fibers. Hence, ECT was well suited for the detection of fiber damage/fiber irregularities. The MSU team developed printed circuit boards (PCB) with coil sensor arrays optimized for NDE of CFRP. Unlike most commercial probes designed for ECT of metallic structures, the MSU array probes were designed for operation in [1-10] MHz frequency range, which was optimal for low-conductive CFRP. Multiple sensing topologies (coil groups excitation/sensing arrangements) were implemented and successfully validated. Capacitive Imaging (CI) technique developed by MSU was complementary to ECT. In contrast to ECT, which was sensitive to local changes of the electrical conductivity, the CI was sensitive to local changes of the dielectric constant. Therefore, CI could provide information about matrix damage/matrix irregularities in CFRP. The MSU CI sensor arrays were made of multiple circular or rectangular open-plate capacitors printed on PCB. Sensors of this type are not commercially available. In addition to ECT and CI, the MSU team developed a hybrid (dual-mode) inductive/capacitive measurement technique that synergistically combined the benefits of inductive and capacitive sensing for rapid NDE of fiber reinforced polymer (FRP) composite structures. Fiber damage and fiber irregularities in FRPs were detected by configuring hybrid sensors as coil sensors. Similarly, matrix damage, matrix irregularities and interlaminar delaminations were detected by configuring hybrid sensors as capacitive sensors. ECT and CI were performed sequentially by means of electronic switching. Hence, eliminating the need for mounting two separate sensor arrays on the probe. Portable robotic platform was developed by MSU for multi-technique high-speed NDE of CFRP test parts. The platform had two 6-axis robots, which enabled inspection of curved parts in approximately a 6×6×6 ft 3 active scan area. On the software side, the MSU team integrated scripts for NDE hardware control with scripts for robot motion control. MSU also implemented automated path planning for the robots, reconstruction of part’s surfaces via stereovision, 3D rendering of inspection data, and image processing algorithms for enhanced defect detection. Automotive composite parts manufactured by Plasan Composites from Phase I were used to validate the ACUT and EMT techniques on representative testbeds. Among those parts were three X-braces for a Dodge Viper, one composite calibration plaque with known defects at known locations, and four other test sections, including sections from a front splitter, a corner section from a composite hood, and a high-pressure RTM panel made using non crimp fabric. Other test samples included CFRP and GFRP calibration plates with fiber/matrix defects fabricated at MSU/CVRC.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The role of information structures in game-theoretic multi-agent learning

Multi-agent learning (MAL) studies how agents learn to behave optimally and adaptively from their experience when interacting with other agents in dynamic environments. The outcome of a MAL process is jointly determined by all agents’ decision-making. Hence, each agent needs to think strategically about others’ sequential moves, when planning future actions. The strategic interactions among agents makes MAL go beyond the direct extension of single-agent learning to multiple agents. With the strategic thinking, each agent aims to build a subjective model of others decision-making using its observations. Such modeling is directly influenced by agents’ perception during the learning process, which is called the information structure of the agent’s learning. As it determines the input to MAL processes, information structures play a significant role in the learning mechanisms of the agents. This review creates a taxonomy of MAL and establishes a unified and systematic way to understand MAL from the perspective of information structures. Here we define three fundamental components of MAL: the information structure (i.e., what the agent can observe), the belief generation (i.e., how the agent forms a belief about others based on the observations), as well as the policy generation (i.e., how the agent generates its policy based on its belief). In addition, this taxonomy enables the classification of a wide range of state-of-the-art algorithms into four categories based on the belief-generation mechanisms of the opponents, including stationary, conjectured, calibrated, and sophisticated opponents. We introduce Value of Information (VoI) as a metric to quantify the impact of different information structures on MAL. Finally, we discuss the strengths and limitations of algorithms from different categories and point to promising avenues of future research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Process Intensification for Direct Conversion of Biomass-Based Syngas to High Octane Gasoline

The conversion of lignocellulosic-derived carbon sources via the creation of a CO2 rich syngas stream provides means to produce liquid hydrocarbon fuels, with advancements intended to create a cost-competitive method of sequential conversion via methanol and dimethyl ether intermediates. Previously developed methods with Cu/Beta zeolites enabled the conversion of dimethyl ether to linear and branched hydrocarbons, with high selectivities towards products with high-octane fuel properties. Through process intensification, use of a single reactor to convert syngas to methanol and DME, then further to hydrocarbons, allows for reduced capitol and operation cost for the same chemistries that typically use 3 reactors instead of one. Utilization of commercially available methanol synthesis catalyst (Megamax 800, Clariant, CZA) and methanol dehydration catalyst with the Cu/Beta catalyst allowed for conversion into hydrocarbon products in this single reactor set up. Reactor configuration with CZA and alumina catalysts mixed with or stacked in a bed physically above the Cu/Beta catalyst provided a system for the direct conversion of the syngas to hydrocarbons with improved yields and selectivities with higher net C1 conversions for stacked bed configurations. Through control of the process conditions, greater conversion of C1 oxygenate intermediates (i.e., methanol and dimethyl ether) was achieved with concomitant increase in selectivities towards gasoline and jet fuel range olefin and paraffin products. Lower hydrocarbon number products and less naphthenes were observed for the conversion of syngas compared to DME conversion on only the Cu/Beta system with prominent selectivities observed for C4, C5, and C7 hydrocarbons. The DME and methanol formation rates, with understanding of the equilibria for both processes, were determined as important factors to allow for improved performance of the Cu/Beta. Determination of the factors which allow for tuning of selectivities and yields created an intensified process which allows for a "market-responsive" biorefinery design, which can produce high octane gasoline or jet fuel range hydrocarbons to meet demands for a more sustainable route to liquid fuels.

BIOMASS FUELS↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2023

Rates and reversibility of CO 2 hydrogenation on Cu-based catalysts

Kinetics of reaction pathways involved in the conversion of CO 2 to methanol and CO on Cu/ZnO/Al 2 O 3 are resolved using in situ chemical titration, steady-state kinetic measurements, and mathematical formalisms for reversibility to probe salient species governing methanol selectivity and yield during CO 2 hydrogenation. Across a range of H 2 :CO 2 = 1:1 to H 2 :CO 2 = 80.5:1, active site density determined from in situ chlorine uptake remained invariant; hence, observed trends in rates can be interpreted as only arising from reaction kinetics and not from changing active site density. Kinetic and thermodynamic contributions to rates are decoupled to evaluate forward and reverse rates of methanol synthesis and reverse water-gas shift (RWGS) reactions. These kinetic analyses show that the forward rates of methanol synthesis exhibit persistent first order dependence on hydrogen pressure and are inhibited by water more significantly than the forward rates of RWGS. In contrast, the reverse rates of methanol synthesis and RWGS are both inhibited by H 2 . Consequently, without any modifications to the Cu/ZnO/Al 2 O 3 catalyst formulation, methanol selectivity can be increased to > 80% by increasing inlet H 2 partial pressure and methanol yield can be enhanced by ~20% by adding water adsorbents even under conditions far from equilibrium. Here, the kinetic treatments presented herein demonstrate a dearth of H* species during catalysis, provide thermodynamic constraints precluding sequential RWGS and CO hydrogenation as the pathway for methanol synthesis, reveal P H2 and P H2O as salient in determining methanol selectivity and yield by impacting both the forward and reverse rates of CO 2 hydrogenation on Cu/ZnO/Al 2 O 3 , and explicate the fundamentals of novel sorption-enhanced methanol synthesis, which not only alleviates equilibrium constraints but also alters the intrinsic rate at which the system approaches equilibrium.

42 ENGINEERING↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2025

Effects of dose rate on the microstructure and deuterium retention in γ-LiAlO 2 pellets

This report presents experimental findings obtained from November 2024 to September 2025. Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by two competing processes: defect production, driven by dose rate, and defect recovery, controlled by defect diffusion, interaction, and annihilation. At a given dose, the resulting microstructural evolution depends on both the dose rate and irradiation temperature. As a continuation of our FY24 tritium science project, which investigated temperature effects at a fixed dose rate, this study focuses on dose-rate effects at a fixed temperature to provide deeper insights into the irradiated microstructure and compositional changes in γ-LiAlO 2 pellets. The study aims to elucidate the impact of dose rate on microstructure, precipitate morphology, deuterium retention, and lithium volatility in γ-LiAlO 2 pellets under sequential 120 keV He + and 80 keV D 2 + ion irradiation. Three dose rates of 7.3×10 4 , 2.9×10 4 and 6.8×10 5 dpa/s were applied to achieve the same total ion fluence of 2×10 17 (He + +D + )/cm 2 at 500 °C, corresponding to a maximum combined dose of 7.55 dpa at ~255 nm. The irradiated pellets were subsequently characterized using scanning transmission electron microscopy (STEM) and time-of-flight secondary ion mass spectrometry (ToF-SIMS). The microstructural response of γ-LiAlO 2 to ion irradiation was found to be strongly dose-rate dependent. At medium and high dose rates, irradiation produced a surface amorphized layer and an underlying crystalline layer containing LiAl 5 O 8 precipitates, with implanted gases accumulating and forming blisters at the crystalline-amorphous interface. It remains to be investigated whether the amorphized layer was produced by He + ion irradiation prior to D 2 + ion irradiation. In contrast, at low dose rates, the material remained crystalline, with cavities, likely gas-filled, distributed around precipitates, within the γ-LiAlO 2 matrix, and along grain boundaries. While precipitate morphology exhibited anisotropy, their size showed little sensitivity to dose rate in the applied range of this study. This result, however, does not rule out the possibility that further reductions in dose rate could influence precipitate size. High dose-rate irradiation enhanced protonium–deuterium isotopic exchange and lithium depletion in the amorphized region. Collectively, the results show that dose rate governs amorphization, gas redistribution, isotopic exchange, and lithium depletion, providing important insights into the mechanisms underlying structural evolution in γ-LiAlO 2 pellets under reactor-relevant irradiation conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Multicriteria-Based Selection of Microalgae Biorefineries: Biomass Composition Defines the Most Suitable Product Portfolios

The supply of microalgae-derived biofuels and bioproducts will be vital in a global-scale bioeconomy. As the diversity of microalgae strains and their compositional plasticity may yield a wide array of products of market interest, the design of effective biorefineries is a central aspect in the path to making microalgae-based products available in the market. In this way, the conversion of microalgae biomass should be planned to employ mature technologies, while maximizing economic and environmental benefits obtained from a diverse product portfolio. This study applied sequential, hybrid Multicriteria Decision Analyses (MCDA) to aid the decision-making process of outlining the most suitable biorefining pathways for specific compositional profiles. For this, the methodology simultaneously considered multiple technical, economic, and environmental criteria, such as market aspects for the main algae-derived products, potential reduction in greenhouse gas (GHG) emissions provided by the biobased alternatives, and technology readiness level of conversion routes, among others. The framework was tested using productivity and compositional data for 13 high-productivity summer strains cultivated under varying nutrient availability. The analysis pointed to compositional profile being a key driver in defining the core biorefining strategy of algae biomass, with lipid-to-carbohydrate ratios higher than roughly 1 warranting the preferential processing of lipids into hydrocarbon fuels with the remainder of the algae biomass compounds being sent to higher-value applications, such as carboxylic acids and renewable thermoplastic substitutes. An in-depth process simulation, techno-economic assessment (TEA), and life-cycle analysis (LCA) effort was carried out as a closing step to corroborate the results stemming from the proposed framework. This study validates the use of MCDAs as a screening method prior to implementing more time-consuming analysis techniques and makes a compelling case for this approach as a product selection tool on a wide range of algae species, thus helping establish species-agnostic (but composition-driven) biorefineries.

biofuels↗
DOE OSTI2026

Mechanism of Tyrosine-Driven Deprotonation in Photosystem II Revealed by Multiscale Simulations

Photosystem II (PSII) drives light-induced water oxidation via stepwise redox transitions of its oxygen-evolving complex (OEC), a Mn 4 Ca cluster advancing through five intermediate S-states (S 0 –S 4 ). The S 2 → S 3 transition involves a redox event in which a Mn ion donates an electron to the redox-active tyrosine YZ, coupled to deprotonation of an OEC-bound water ligand─yet the underlying coupling mechanism remains unresolved. Time-resolved serial femtosecond crystallography (TR-SFX) has revealed transient electron density shifts near the redox-active tyrosine Y Z , interpreted as sequential oxidation and reduction, with reduction initiating ∼1 μs after excitation and substantially progressed by 30 μs. However, this interpretation conflicts with kinetics from photothermal beam deflection (PBD), time-resolved X-ray absorption spectroscopy (TR-XAS), and electron paramagnetic resonance (EPR), which place electron transfer at 190–400 μs and proton transfer around 30 μs. Here, we reconcile these discrepancies using quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations. We show that oxidation of P680 and Y Z breaks the symmetry of the nearby hydrogen bonds involving water molecule W4, displacing Y Z and replicating the TR-SFX features of Y Z and Q165 observed at 1 μs. This local perturbation propagates through a hydrogen-bond network, transmitting the electrostatic signal from Y Z to the E65-E312 dyad and triggering redox-coupled deprotonation via the Cl1 channel. By 30 μs, the hydrogen-bond symmetry is restored through deprotonation of W2 (or alternatively W1), reproducing the disappearance of TR-SFX density differences around Y Z and Q165 without requiring Y Z reduction. Our proposed mechanism also gives molecular insights into the O6* density, assigning it to water reorganization rather than a discrete Ca-bound hydroxide species. Here, our results reveal a detailed atomistic mechanism linking Y Z oxidation to long-range proton release and suggest a functional role for the nearby Cl – ion in proton transfer. More broadly, this study underscores the importance of hydrogen-bond dynamics in mediating redox-driven proton transport and demonstrates how integrative simulations can resolve mechanistic ambiguities.

Deprotonation↗
DOE OSTI2023

A Deep Learning Filter for the Intraseasonal Variability of the Tropics

Abstract This paper presents a novel application of convolutional neural network (CNN) models for filtering the intraseasonal variability of the tropical atmosphere. In this deep learning filter, two convolutional layers are applied sequentially in a supervised machine learning framework to extract the intraseasonal signal from the total daily anomalies. The CNN-based filter can be tailored for each field similarly to fast Fourier transform filtering methods. When applied to two different fields (zonal wind stress and outgoing longwave radiation), the index of agreement between the filtered signal obtained using the CNN-based filter and a conventional weight-based filter is between 95% and 99%. The advantage of the CNN-based filter over the conventional filters is its applicability to time series with the length comparable to the period of the signal being extracted. Significance Statement This study proposes a new method for discovering hidden connections in data representative of tropical atmosphere variability. The method makes use of an artificial intelligence (AI) algorithm that combines a mathematical operation known as convolution with a mathematical model built to reflect the behavior of the human brain known as artificial neural network. Our results show that the filtered data produced by the AI-based method are consistent with the results obtained using conventional mathematical algorithms. The advantage of the AI-based method is that it can be applied to cases for which the conventional methods have limitations, such as forecast (hindcast) data or real-time monitoring of tropical variability in the 20–100-day range.

Stan, Cristiana↗
DOE OSTI2026

Implementing a unified solver for nonlinearly constrained optimization

SQP and interior-point methods (also referred to as Lagrange-Newton methods) typically share key algorithmic components, such as strategies for computing descent directions and mechanisms that promote global convergence. Building on this insight, we introduce a unifying framework with eight building blocks that abstracts the workflows of Lagrange-Newton methods. We then present Uno, a modular C++ solver that implements our unifying framework and allows the automatic combination of a wide range of strategies with no programming effort from the user. Uno is meant to (1) organize mathematical optimization strategies into a coherent hierarchy; (2) offer a wide range of efficient and robust methods that can be compared for a given instance; (3) enable researchers to experiment with novel optimization strategies; and (4) reduce the cost of development and maintenance of multiple optimization solvers. Uno’s software design allows user to compose new customized solvers for emerging optimization areas such as robust optimization or optimization problems with complementarity constraints, while building on reliable nonlinear optimization techniques. We demonstrate that Uno is highly competitive against state-of-the-art solvers filterSQP, IPOPT, SNOPT, MINOS, LANCELOT, LOQO, and CONOPT on a subset of 429 small problems from the CUTE collection. Uno is available as open-source software under the MIT license at https://github.com/cvanaret/Uno and via its C, Julia, Python, Fortran, and AMPL interfaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗