Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2025

Characterizing leaf-scale fluorescence with spectral invariants

Sun-induced chlorophyll fluorescence (SIF) is increasingly recognized as a non-destructive probe for tracking terrestrial photosynthesis. Emerging developments in spectral invariants theory provide an innovative and efficient approach for representing SIF radiative transfer processes at the canopy scale. However, modeling leaf-scale fluorescence based on the spectral invariants properties (SIP) remains underexplored. In this study, the spectral invariants theory is employed for the first time to model the leaf-scale total, backward and forward fluorescence (leaf-SIP SIF). The leaf-SIP SIF model separates the leaf-scale radiative transfer process into two distinct components: the wavelength-dependent one associated with leaf biochemical properties, and the wavelength-independent component linked to leaf structural characteristics. The leaf structure-related effects are characterized by two spectrally invariant parameters: the photon recollision probability (p) and the scattering asymmetry parameter (q), which are parameterized using the directly measurable leaf dry matter. Evaluation against field measurements shows that the proposed leaf-SIP SIF model has a good performance, with coefficient of determination (R 2 ) of 0.89, 0.89, 0.90 and root mean squared errors (RMSE) of 1.28, 0.69, 0.74 Wm -2 µm -1 sr -1 , respectively for the total, backward, and forward fluorescence (660–800 nm). The leaf-SIP SIF model with a more concise formulation demonstrates comparable performance with the widely used Fluspect model. Further, the leaf-SIP SIF model provides a simple and efficient approach for simulating leaf-scale fluorescence, with the potential to be integrated into a unified SIP-based model framework for simulating the radiative transfer processes across the soil-leaf-canopy-atmosphere continuum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Contrasting Effects of Catecholate and Hydroxamate Siderophores on Molybdenite Dissolution

Molybdenum (Mo) is essential for many enzymes but is often sequestered within minerals, rendering it not readily bioavailable. Metallophores, metabolites secreted by microorganisms and plants, promote mineral dissolution to increase metal bioavailability. However, interactions between metallophores and Mo-bearing minerals remain unclear. Here, in this study, catecholate protochelin and hydroxamate desferrioxamine B (DFOB) were utilized to examine their effects on dissolution of the common Mo-bearing mineral, molybdenite (MoS 2 ), under both oxic and anoxic conditions. Protochelin promoted molybdenite dissolution under oxic conditions with the formation of MoO 3 on the surface and Mo-siderophore complexes in solution. This was attributed to air-oxidation of both molybdenite and protochelin, as evidenced by lack of dissolution under anoxic conditions but enhanced dissolution by either pre-oxidized protochelin or pre-oxidized molybdenite. Liquid chromatography-mass spectroscopy (LC-MS), X-ray photoelectron spectroscopy (XPS), and Time-of-flight secondary ion mass spectrometry (ToF-SIMS) analyses revealed degradation of protochelin and adsorptions of its byproducts on molybdenite surface to promote dissolution. Conversely, DFOB inhibited molybdenite dissolution under both oxic and anoxic conditions, likely attributed to surface adsorption of DFOB and its weak complexation with Mo(VI) at the circumneutral pH. This work highlights the need to consider the balance between promoting and inhibitory effects of different metallophores on Mo-mineral dissolution.

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

Using SF 6 and Xe to Monitor Gas Migration Through Explosion-Generated Fracture Networks

Here we describe a field study where tracer gas was injected into a subsurface cavity created by a small chemical explosion beneath the water table. The main objective of the study is to compare the migration of sulfur hexafluoride (SF 6 ) and xenon (Xe) through an explosion-generated fracture network and to study the influence of ground water on gas transport. A mixture of tracer gases (50% of SF 6 and 50% of Xe) was injected on 31 October 2018 and gas sampling continued until 8 November 2018. We observe similar trends in SF 6 and Xe concentrations at four ground surface sampling sites. The changes in the SF 6 /Xe ratios with time show that more SF 6 than Xe is observed during the barometric pressure lows when the absolute measured concentrations are highest. Conversely, the ratio SF 6 /Xe is less than 1 during the high-pressure intervals when absolute measured concentrations are low. The results of the experiment suggest that during barometric pressure lows the tracer is migrating to the surface primarily by advective gas phase transport, whereas during barometric pressure highs, advection is suppressed and near-surface evaporation of interstitial pore fluid with tracer dissolved in it becomes more important. Thus, the results of the experiment show that the gas concentrations at the surface are controlled by the combined effects of the gas dissolution into pore water and the barometric pressure fluctuations.

58 GEOSCIENCES↗
DOE OSTI2020

Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones

Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one-way coupling between the surface and subsurface domains. Such one-way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam, was developed using the open-source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod-type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three-dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.

58 GEOSCIENCES↗
DOE OSTI2020

Semiquantitative Estimates of Rainfall Variability During the 8.2 kyr Event in California Using Speleothem Calcium Isotope Ratios

A multiproxy record from a fast-growing stalagmite reveals variable hydroclimate on the California coast across the 8.2 kyr event and a precursor event likely caused by initial drainage of proglacial Lake Agassiz. Using speleothem δ 44 Ca, we develop the first semiquantitative estimates of paleorainfall variability for California through calibration with measurements of the modern climate and cave environment. We find that the magnitude of rainfall variability during the 8.2 kyr event approached the multiyear variability observable in the recent past (1950–2019) and the magnitude of variability during the precursor event likely exceeded this range. Additionally, we observe other instances of multidecadal variability comparable in magnitude to the precursor event during the record. Finally, our work suggests that speleothem calcium isotope ratios are a powerful semiquantitative means to reconstruct paleorainfall, although numerous factors must be assessed in each cave system before applying this approach.

58 GEOSCIENCES↗
DOE OSTI2025

Using Radiogenic Noble Gas Nuclides to Identify and Characterize Rock Fracturing

Abstract Fracture‐released radiogenic noble gas nuclides are used to identify locations and constrain the volume of new fracture creation during subsurface detonations. Real‐time, in situ noble gases and reactive gases were monitored using a field‐deployed mass spectrometer and automated sampling system in a multilevel borehole array. Released gases were measured after two different detonations having distinct energy, pressure, and gas volume characteristics. Explosive‐derived gases (N 2 O, CO 2 ) and excess radiogenic 4 He and 40 Ar above atmospheric background are used to identify locations of gas transport and new fracture creation after each detonation. Fracture‐released radiogenic 4 He is used to constrain the volume of newly created fractures with a model of helium release from fracturing. Explosive by‐product gas was observed in multiple locations both near and distal to the shot locations for both detonations. Radiogenic 4 He and 40 Ar release from rock damage was observed in locations near the detonation after the second, more powerful detonation. Observed 4 He response is consistent with a model of diffusive release from newly created fractures. Volume of new fractures estimated from the 4 He release ranges from 1 to 5 m 2 with apertures ranging from 0.1 to 1 m. Our results provide evidence that radiogenic noble gases released during fracture creation can be identified at the field scale in real time and used to identify timing and location of fracture creation during deformation events. This technique could be useful in subsurface science and engineering problems where the location and amount of newly created rock fracturing is of interest including fault rupture, mine safety, subsurface detonation monitoring and reservoir stimulation.

58 GEOSCIENCES↗
DOE OSTI2025

A Unified Analytical Model for Pressure Solution With Fully Coupled Diffusion and Reaction

Abstract Geophysical models for pressure solution are typically developed for diffusion‐controlled or reaction‐controlled scenarios. We present a unified analytical model that considers fully coupled diffusion and reaction during pressure solution. The model recovers the diffusion‐controlled and reaction‐controlled models in the literature as specific limiting cases. When diffusion and reaction exhibit comparable influences, we validate the proposed model against independent numerical simulations. The proposed model is then employed in interpreting experimental measurements, demonstrating a better agreement compared to previous interpretations.

Wang, Ziyan [Energy Geosciences Division Lawrence ↗
DOE OSTI2026

A Mountain Glacier Perspective on the Bipolar Seesaw

A global record of mountain glacier terminations during the last deglaciation (∼19–11 ka) dated by a large, uncurated data set of cosmogenic-nuclide exposure ages highlights a statistically significant asynchrony in termination ages between the Northern and Southern Hemispheres. This interhemispheric offset in the timing of glacier terminations is consistent with previously correlated ice core records that show a systematic interhemispheric lag in the timing of abrupt climate events, with the Southern Hemisphere leading the Northern Hemisphere by ∼300–3,500 years. Our analysis (a) aggregates cosmogenic-nuclide exposure ages from a global data set of moraines to discern climatically driven peaks in moraine emplacement events, and (b) utilizes a Monte Carlo simulation based on a null hypothesis that moraine emplacement is interhemispherically synchronous to estimate the statistical significance of the observed offset. The observed lag of Northern Hemisphere emplacement events compared to the Southern Hemisphere is statistically significant and is consistent with the “bipolar seesaw” pattern observed in ice core records.

58 GEOSCIENCES↗
DOE OSTI2026

Process‐Oriented Calibration of a Turbulence Scheme in the DOE's Global Storm‐Resolving Model Using Machine Learning

A process‐oriented calibration framework is developed for the Simplified Higher‐Order Closure (SHOC) turbulence scheme in DOE's Simple Cloud Resolving E3SM Atmospheric Model (SCREAM). This framework leverages machine learning surrogates and observational constraints to efficiently calibrate SHOC adjustable parameters across two convective regimes: clear‐sky dry convective boundary layer and fair‐weather shallow cumulus clouds from ARM observations. We use perturbed‐parameter ensembles of a doubly periodic version of SCREAM to train surrogates and apply Markov Chain Monte Carlo sampling guided by cost functions based on benchmarking large‐eddy simulations and observations to identify optimized parameter sets that perform well in both regimes. The calibrated SHOC parameters substantially improve boundary‐layer turbulence and cloud boundaries, and modeled cloud fraction and radiative effects align better with observations than the default. These results demonstrate that combining multiple process‐specific convective regimes with machine‐learning surrogates can reduce parametric uncertainties and yield a model more faithful to cloud–turbulence interactions.

58 GEOSCIENCES↗
DOE OSTI2025

The Energy Exascale Earth System Model Version 3: 1. Overview of the Atmospheric Component

This paper describes the atmospheric component of the US Department of Energy's Energy Exascale Earth System Model (E3SM) version 3. Significant updates have been made to the atmospheric physics compared to earlier versions. Specifically, interactive gas chemistry has been implemented, along with improved representations of aerosols and dust emissions. A new stratiform cloud microphysics scheme more physically treats ice processes and aerosol‐cloud interactions. The deep convection parameterization has been largely improved with sophisticated microphysics for convective clouds, making model convection sensitive to large‐scale dynamics, and incorporating the dynamical and physical effects of organized mesoscale convection. Improvements in aerosol wet removal processes and parameter re‐tuning of key aerosol and cloud processes have improved model aerosol radiative forcing. The model's vertical resolution has increased from 72 to 80 layers with the extra eight layers added in the lower stratosphere to better simulate the Quasi‐Biennial Oscillation. These improvements have enhanced E3SM's capability to couple aerosol, chemistry, and biogeochemistry and reduced some long‐standing biases in simulating tropical variability. Compared to its predecessors, the model shows a much stronger signal for the Madden‐Julian Oscillation, Kelvin waves, mixed Rossby‐gravity waves, and eastward inertia‐gravity waves. Aerosol radiative forcing has been considerably reduced and is now better aligned with community best estimates, leading to significantly improved skill in simulating historical temperature records. Its simulated mean‐state climate is largely comparable to E3SMv2, but with some notable degradation in shortwave cloud radiative effect, precipitable water, and surface wind stress, which will be addressed in future updates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗
DOE OSTI2021

Standing shock prevents propagation of sparks in supersonic explosive flows

Abstract Volcanic jet flows in explosive eruptions emit radio frequency signatures, indicative of their fluid dynamic and electrostatic conditions. The emissions originate from sparks supported by an electric field built up by the ejected charged volcanic particles. When shock-defined, low-pressure regions confine the sparks, the signatures may be limited to high-frequency content corresponding to the early components of the avalanche-streamer-leader hierarchy. Here, we image sparks and a standing shock together in a transient supersonic jet of micro-diamonds entrained in argon. Fluid dynamic and kinetic simulations of the experiment demonstrate that the observed sparks originate upstream of the standing shock. The sparks are initiated in the rarefaction region, and cut off at the shock, which would limit their radio frequency emissions to a tell-tale high-frequency regime. We show that sparks transmit an impression of the explosive flow, and open the way for novel instrumentation to diagnose currently inaccessible explosive phenomena.

58 GEOSCIENCES↗
DOE OSTI2023

Assessing reservoir performance for geologic carbon sequestration in offshore saline reservoirs

Carbon capture and storage (CCS) is a technological strategy to reduce CO 2 emissions from hydrocarbon-based energy production and non-energy industrial emitters, like fertilizer and cement production. To expand the geographic opportunities for CCS, offshore geologic settings are increasingly considered for saline CO 2 storage, particularly where legacy oil and gas infrastructure may be repurposed for CO 2 transport and injection. In this context, the Gulf of Mexico (GoM), United States (U.S.), may offer tremendous opportunity for CCS because the surrounding states are among the largest CO 2 producers in the U.S. and there has been extensive oil and gas development over the last half century. Nevertheless, there remains significant uncertainty in how GoM reservoir/seal systems will respond to industrial-scale CO 2 injections, particularly in the context of meso-scale (1–10s of meters) permeability variations that are difficult to identify prior to injection. Such permeability variations can have significant impacts on CO 2 migration and fluid pressure propagation, both of which govern the overall effectiveness of the CCS project. This research uses ensemble simulation methods to quantify the influence of spatially variable and a priori unknown permeability fields on CO 2 plume development and fluid pressure propagation when CO 2 is injected at 1 MMt per year for ten years. Results show a number of characteristic patterns that reveal the substantial influence of both near- and far-field permeability on reservoir performance. By focusing on basin characteristics that are common among offshore basins worldwide, results from this study suggest that offshore CCS may be a feasible carbon management strategy in regions where offshore basins are proximal to industrial CO 2 sources.

42 ENGINEERING↗
DOE OSTI2021

Network analysis of earthquake ground motion spatial correlation: a case study with the San Jacinto seismic nodal array

SUMMARY The spatial correlation of earthquake ground motion intensity can be measured from strong motion data; however, the data used in past studies is sparsely sampled in space, and only the interstation distance was considered as a correlation variable. These limitations mean that we have only weak constraints on the true correlation structure of ground motion and that potentially important aspects of spatial correlation are unconstrained. In this study, we combine a large-N seismic array and graph analytics to explore this issue at a local scale using small local and regional earthquakes. Our result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements, of seismic wavefields. Aftershock sequences of major earthquakes would provide particularly data-rich targets of opportunity.

58 GEOSCIENCES↗
DOE OSTI2023

Transforming microseismic clouds into near real-time visualization of the growing hydraulic fracture

SUMMARY Microseismic observations during unconventional reservoir stimulation are typically seen as a proxy for clusters of hydraulic fractures and the extent of the stimulated reservoir. Such straightforward interpretation is often misleading and fails to provide a physically reasonable image of the fracturing process. This paper demonstrates the application of a physics-based machine learning algorithm which enables a rapid and accurate fracture mapping from the microseismic data. Our training and validation data set relies on a history-matched geomechanical modelling workflow implemented in GEOS software for the Hydraulic Fracturing Test Site 1 (HFTS-1) project. For this study we augmented the simulated fracture growth through geostatistical modelling of induced seismicity, so that the synthetic microseismic catalogue matches the main statistical properties of the field observations. We formulated the problem of mapping the actual fracture in the clutter of events to parallel common video segmentation workflows: several past video frames (microseismic density snapshots) are passed through a deep convolutional network to classify whether a given voxel is associated with a fracture or intact rock. We found that for accurate fracture mapping, the network’s input and architecture must be augmented to incorporate the fluid injection parameters (pressure, rate, concentration of proppant, and location of the perforation within the cluster). The error rate for the network reached as little as 10 per cent of the fracture area, while a conventional microseismic interpretation approach yielded ∼300 per cent. Our approach also yields must faster predictions than conventional methods (minutes instead of weeks), and could enable engineers to make rapid decisions regarding engineering parameters (pumping rate, viscosity) in real time during stimulation.

58 GEOSCIENCES↗
DOE OSTI2025

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗