Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics and Computing, Geosciences”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multiscale Dynamics of Reactive Fronts in the Subsurface

Understanding and predicting flow and reactive transport in rocks (i.e. geologic porous media) is critical to many technologies at the heart of the energy transition, including CO 2 sequestration and H2 storage. However, accurate modeling and prediction of these systems is very complex because physico-chemical processes that occur at very small spatial scales, i.e. in the pores of the rocks, can dramatically control the system performance at the field scale (km). For example, precipitation reactions at the pore-scale can lead to large permeability changes at the field scale and dramatically alter the migration of stored gases. Properly accounting for multi-scale coupling effects is critical to achieve predictivity and confidence in model outputs, which then can guide design and optimization at the system-scale. This can be achieved through the development of rigorous mathematical models that can appropriately account for fine-scale effects at the large scale. The final report of the Early Career award DE-SC0019075 “Multiscale dynamics of reactive fronts in the subsurface” summarizes the mathematical, numerical and experimental advancements to study and predict reactive transport in geologic porous media across scales.

58 GEOSCIENCES↗
DOE OSTI2024

FECM/NETL CO 2 Saline Storage Cost Model (2024): User’s Manual

The U.S. Department of Energy's (DOE) Office of Fossil Energy and Carbon Management (FECM), in collaboration with the National Energy Technology Laboratory (NETL), has developed the FECM/NETL CO 2 Saline Storage Cost Model (CO2_S_COM). This Excel-based tool provides a comprehensive framework for estimating the costs and breakeven prices associated with storing carbon dioxide (CO 2 ) in deep saline formations. Designed from the perspective of a CO 2 storage site owner, the CO2_S_COM incorporates four integrated modules—project management, financial analysis, activity cost estimation, and geological evaluation—to deliver fast, robust, and actionable insights for evaluating project finances. This is the user's manual for CO2_S_COM. The model may be accessed at this link: FECM/NETL CO2 Saline Storage Cost Model CO2_S_COM 2024 (v4) - Submissions - EDX

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Model Card for WaveDenoiser

This study used STEAD to train and evaluate this model because STEAD is among the best benchmark datasets available for local to regional data. STEAD is a global dataset with over 1 million 60 s long seismic waveforms that originated from approximately 450,000 earthquakes and background noise captured by more than 2,500 seismic stations. Each waveform in STEAD was attached with metadata such as earthquake locations, station locations, and signal arrival times when available

58 GEOSCIENCES↗
DOE OSTI2026

Dynamic Networks Experiment 2: Measuring Associator Sensitivity to Signal Detection Errors

Using the Dynamic Networks Experiment 2 (DNE2) human-analyst event bulletin picks as a baseline signal detection dataset, we generate 47 additional datasets by gradually reducing their accuracy and completeness by randomly removing DNE2 picks, changing the initial phase labels from P to S and vice-versa, and injecting noise detections to simulate real-world signal detection algorithms.

58 GEOSCIENCES↗
DOE OSTI2026

Strategic Petroleum Reserve Cavern Leaching Monitoring CY24

This report provides an analysis of the effects of raw water leaching at the U.S. SPR using the Sandia Solution Mining Code (SANSMIC). A new version of the code has been established in the past year (implemented in Python and C++ code) and was used for this annual report. Additionally, the setup, running, and post-processing of leaching modeling has now been streamlined into a single Python notebook environment for each cavern run.

58 GEOSCIENCES↗
DOE OSTI2026

Superstructure Optimization for Brine Valorization from Brackish Water Desalination

This poster presents preliminary results from a superstructure optimization framework developed to identify cost-optimal brine valorization configurations for brackish water desalination plants across diverse U.S. regional feed chemistries. The study uses brackish groundwater compositions from Arizona, California, Florida, New Mexico, and Texas. Using Pyomo Generalized Disjunctive Programming (GDP) within the WaterTAP modeling environment, the optimization framework simultaneously evaluates thousands of candidate treatment configurations, spanning nanofiltration, reverse osmosis, and chemical precipitation, to minimize the levelized cost of water (LCOW) while meeting water recovery targets and product recovery constraints. Results across eight representative feed clusters demonstrate water recovery rates of 57–87% and net LCOW values ranging from -$0.032/m³ (net revenue-positive) to $0.80/m. Notably, no single process configuration was optimal across all feed types, underscoring the necessity of feed-specific optimization. Products targeted include calcium carbonate (CaCO₃) at $0.01/kg and sodium chloride (NaCl) at $0.10/kg, both at 95% purity, with product revenues offsetting treatment costs in several scenarios. The work advances NAWI's process systems engineering capabilities for multi-configuration screening.

58 GEOSCIENCES↗
DOE OSTI2025

Designing Yield Grids for Airborne Nuclear Explosions to Maintain Decrements in Rayleigh Waveform Amplitudes

This document details a method for selecting grid points that sample the yields of airborne, low yield nuclear explosions that constrains the relative size of seismic waveforms sourced by explosions that are parameterized by this grid. In detail, we impose this constraint so that the vertical component of ground displacement that is sourced by explosions with a fixed height of burst (HoB) at grid point k change by a fixed factor of a when compared to the waveform amplitudes sourced by explosions at grid point k + 1, as observed at a fixed standoff distance. This gridding thereby preserves amplitude ratios along the grid so that they are agnostic to k. We assert that this method is useful in tests that apply algorithms against synthetic or augmented datasets for which the performance of the algorithm must be measured against source size, and waveform amplitudes must change uniformly (e.g., Berg et al, 2025).

58 GEOSCIENCES↗
DOE OSTI2026

Subsurface mafic and ultramafic rock mapping and analysis for carbon mineralization in the United States (SubMAP-CO 2 )

Achieving durable, gigaton-scale carbon dioxide (CO₂) removal will require the development of storage strategies that are both permanent and scalable. Among the available options, carbon mineralization in mafic and ultramafic rocks offers a uniquely robust pathway for long-term CO₂ sequestration, converting CO₂ into thermodynamically stable carbonate minerals. While surface exposures of these lithologies have been extensively studied, particularly in ophiolitic and volcanic settings, the potential for carbon mineralization in the subsurface remains comparatively underexplored at national scale. This gap limits the ability to assess storage capacity, accessibility, and deployment pathways for mineralization-based carbon management in the United States.

08 HYDROGEN↗
DOE OSTI2026

MultiPEM Toolbox: User Manual

This document explains use of the Multi-Phenomenology Explosion Monitoring (Multi PEM) Toolbox, a collection of R scripts for estimating the unknown device parameters of a new event with uncertainty quantification. The methodology and application used for il lustration in this user manual are published in Geophysical Journal International and fully documented in a Los Alamos National Laboratory technical report hereafter designated “WPA” for reference. Additional details on the application are found in a recent journal article. Two assessment types are available: rapid and complete.

58 GEOSCIENCES↗
DOE OSTI2026

Yield Estimates of three Historical Atmospheric Nuclear Explosions at Lop Nor, from Reduced Order Models of Regionally Recorded Rayleigh Waveforms and One Seismic Station: A Brief Communication

We use a reduced order model (ROM) for Rayleigh waveforms sourced by large atmospheric explosions along with data collected from a single seismic station (TLG) to estimate the yields of three atmospheric nuclear tests with order-megaton (MT) yields. These historical tests were conducted by China between 1973 and 1980 (CHIC 15, CHIC 16, and CHIC 26). We use our ROM to estimate yields of 2.34 MT (CHIC 15), 0.51 MT (CHIC 16), and 1.15 (CHIC 26) in the 20.5 s period band that are consistent with publicly accepted values of 3 MT (CHIC 15), 0.50 MT (CHIC 16), and 0.730 MT (CHIC 26).

58 GEOSCIENCES↗
DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

c-HAND: near real-time coastal flood mapping

The Texas Gulf Coast region contains significant centers of population, infrastructure, and economy and is threatened by intensifying tropical storms. The flooding from these tropical storms often has multiple compounding drivers. This characteristic presents a complex numerical problem where a simulation must consider multiple hydrologic forcings. While several procedures exist for addressing this problem numerically, they tend to be resource-intensive and cannot be conducted in near real-time. We extend GeoFlood, a reduced physics approach for fluvial flood forecasting, to rapidly predict coastal and compound fluvial-coastal inundation. This method is validated against a numerical ocean circulation model (ADCIRC) simulation of Hurricane Ike, a major coastal flooding event that happened on the Texas Gulf Coast in 2008. We show that the inundation map generated by coastal HAND (c-HAND) has reasonable agreement with the ADCIRC simulation while taking about 1.7% of the time currently needed to run ADCIRC on a supercomputer. While our model correctly predicts 99% of ADCIRC-inundated DEM cells, it also overpredicts inundated area by a factor of approximately 27%. We combine c-HAND with the GeoFlood framework for fluvial flood forecasting to create a compound fluvial-coastal inundation mapping workflow that can be run in near real-time. c-HAND's fast wall-clock time and low CPU requirements can support decision making by first response personnel. The method provides timely and convenient access to crucial information, such as the locations of flooded roads and inundated coastal areas.

42 ENGINEERING↗
DOE OSTI2024

Unconventional Wells Interference: Supervised Machine Learning for Detecting Fracture Hits

The primary objective of the study was development of a machine learning (ML)-based workflow for fracture hit (“frac hit”) detection and monitoring using shale oil-field data such as drilling surveys, production history (oil and produced water), pressure, and fracking start time and duration records. The ML method takes advantage of long short-term memory (LSTM) and multilayer perceptron (MLP) neural networks to identify the frac hits due to hydraulic communication between the fracking child well(s) and the producing parent well(s) within the same pad (intra-pad interaction) and/or on different pads (inter-pad interaction). It utilizes time series of pressure and production data from within a pad and from adjacent pads. The workflow can capture time variable features of frac hits when the model architecture is deep and wide enough, with enough trainable parameters for deep learning and feature extraction, as demonstrated in this paper by using training and testing subsets of the field data from selected neighboring pads with over a couple of hundred wells. The study was focused on frac-hit interaction among paired wells and demonstrated that the ML model, once trained, can predict the frac-hit probability.

58 GEOSCIENCES↗
DOE OSTI2020

Power and Wind Shear Implications of Large Wind Turbine Scenarios in the US Central Plains

Continued growth of wind turbine physical dimensions is examined in terms of the implications for wind speed, power and shear across the rotor plane. High-resolution simulations with the Weather Research and Forecasting model are used to generate statistics of wind speed profiles for scenarios of current and future wind turbines. The nine-month simulations, focused on the eastern Central Plains, show that the power scales broadly as expected with the increase in rotor diameter (D) and wind speeds at hub-height (H). Increasing wind turbine dimensions from current values (approximately H = 100 m, D = 100 m) to those of the new International Energy Agency reference wind turbine (H = 150 m, D = 240 m), the power across the rotor plane increases 7.1 times. The mean domain-wide wind shear exponent (α) decreases from 0.21 (H = 100 m, D = 100 m) to 0.19 for the largest wind turbine scenario considered (H = 168 m, D = 248 m) and the frequency of extreme positive shear (α > 0.2) declines from 48% to 38% of 10-min periods. Thus, deployment of larger wind turbines potentially yields considerable net benefits for both the wind resource and reductions in fatigue loading related to vertical shear.

17 WIND ENERGY↗