Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface data models”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Groundwater Drought: Environmental Controls and Monitoring

Arising from a long-term deficit of precipitation, groundwater drought often lags meteorological drought by months to years, depending on climate and subsurface hydrogeological conditions. Due to the paucity of in situ groundwater data in many parts of the world, global monitoring of groundwater storage variations and drought may be best achieved using satellite observations and/or groundwater time series simulated by hydrological models. The Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO) missions have greatly benefitted the modeling and monitoring of groundwater storage changes and drought at the global scale. In this chapter, we first review environmental controls on the temporal variability of groundwater using in situ data. We then describe an approach that infuses GRACE and GRACE-FO observations into a land surface model for assessing groundwater storage changes and drought globally. We also discuss characteristics of simulated groundwater drought and the limitations of current groundwater drought monitoring approaches.

Bailing Li↗
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2024

On the transferability of residence time distributions in two 10-km long river sections with similar hydromorphic units

Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface and their residence time distributions (RTDs) in the subsurface are important for managing the water quality and ecosystem health in dynamic river corridors. However, direct simulating high-spatial resolution HEFs and RTDs can be time-consuming, especially for watershed-scale modeling. Efficient surrogate models linking RTDs to hydromorphic units (HUs) can be alternatives for simulating RTDs in large-scale models. A common concern of these surrogate models, though, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this issue, this work evaluates the HEFs and resulting RTD-HU relationships for two 10-km long river corridors along the Columbia River leveraging a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework we previously developed. Applying such a framework at the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. Finally, our comparison shows that the similarity and transferability of the RTD-HU relationship is very low for the two investigated river sections, which suggests that devising a general algorithm to estimate RTDs based solely on surface water hydrodynamics and short-distance river channel topography data, as well as HU classification, might be nearly impossible.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

Mars Exploration Rovers as Virtual Instruments for Determination of Terrain Roughness and Physical Properties

The two 2003 Mars Exploration Rovers (MERs), in combination with the Athena Payload, will be used as virtual instrument systems to infer terrain properties during traverses, in addition to using the rover wheels to excavate trenches, exposing subsurface materials for remote and in-situ observations. The MERs are being modeled using finite element-based rover system transfer functions that utilize the distribution of masses associated with the vehicle, together with suspension and wheel dynamics, to infer surface roughness and mechanical properties from traverse time series data containing vehicle yaw, pitch, roll, encoder counts, and motor currents. These analyses will be supplemented with imaging and other Athena Payload measurements. The approach is being validated using Sojourner data, the FIDO rover, and experiments with MER testbed vehicles. In addition to conducting traverse science and associated analyses, trenches will be excavated by the MERs to depths of approximately 10-20 cm by locking all but one of the front wheels and rotating that wheel backwards so that the excavated material is piled up on the side of the trench away from the vehicle. Soil cohesion and angle of internal friction will be determined from the trench telemetry data. Emission spectroscopy and in-situ observations will be made using the Athena payload before and after imaging. Trenching and observational protocols have been developed using Sojourner results; data from the FIDO rover, including trenches dug into sand, mud cracks, and weakly indurated bedrock; and experiments with MER testbed rovers. Particular attention will be focused on Mini-TES measurements designed to determine the abundance and state of subsurface water (e.g. hydrated, in zeolites, residual pore ice?) predicted to be present from Odyssey GRS/NS/HEND data.

R E Arvidson↗
NASA NTRSDate not supplied

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain-Corrected Satellite Sea Surface Salinity Using the GMOA S2s Forecast System

The ENSO phenomenon has a significant global socio-economic impact and has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help improve the initialization of the thermocline, while satellite SST aids in constraining surface heat-fluxes, leading to improved coupled system sub-seasonal to seasonal forecasts. However, few studies have focused on improving the near-surface density and mixing through satellite sea surface salinity (SSS) assimilation. The few ocean models that assimilate satellite SSS, bias correct to normalize towards the near-surface Argo data for expediency. This assumption is likely inadequate in rainy regions, where buoyant water forms a fresh surface lens. In previous work, we showed that adjusting SSS to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 improves the near-surface density and mixed layer depth, leading to deeper thermocline and improved NINO3.4 SST forecasts. We now utilize the Soil Moisture and Ocean Salinity rain-corrected (SMOS_RC) SSS, available in SMOS-CATDS products, to represent Sb more accurately at the first model layer (e.g., 5 m). Rather than a diffusivity model as RIM, SMOS_RC uses a statistical correction dependent on Integrated Multi-satellitE Retrievals for GPM (IMERG) rain rates, established on observed SMOS SSS decreases related to Sb in the presence of rain (Supply et al., 2020). For all experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is, and a separate reanalysis assimilates SMOS_RC. We assess the impact on near-surface and subsurface dynamics by validating against observations and explore how SSS assimilation (SMOS vs SMOS_RC) impacts ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (S2S-v3, Molod et al. 2020). We show that improved estimates of density and near-surface mixing led to more accurate coupled air/sea interaction and better ENSO forecasts. The increased SSS, resulting from the removal of the instantaneous rain effect, modifies the ocean state by enhancing mixing and deepening the thermocline.

Veronica Ruiz Xomchuk↗
DOE OSTI2020

Microbial Community Composition in Deep‐Subsurface Reservoir Fluids Reveals Natural Interwell Connectivity

Abstract The identification of natural fractures and the wells they connect is crucial for the development of geological reservoirs because it may have an important impact on reservoir model construction and hydraulic fracture propagation. In this study we investigated the use of a novel data source, the microbial community composition in the reservoir formation fluids, for identification of interwell connectivity caused by natural fractures. We verified this concept at a newly developed mesoscale enhanced geothermal system (EGS) field testbed located 4,850 ft (1,478 m) beneath ground surface at the Sanford Underground Research Facility in Lead, SD. Fluids produced at or near the EGS testbed were sampled and subjected to high‐throughput 16S rRNA gene amplicon sequencing to analyze the microbial community profile therein. Despite the typically substantial heterogeneity across the community profiles of samples spatially distributed (10 m to 1.9 km apart) throughout the site, samples from two wells at the EGS testbed showed highly similar microbial community composition, suggesting the two wells intersected the same natural fracture. This evidence of natural connectivity between the two wells at the EGS testbed was later corroborated by core log analysis and sewer camera surveys into the boreholes. Besides the field case described in this study, microbial community composition as a reservoir diagnostic tool would be applicable in a much broader context such as unconventional hydrocarbon exploitation, groundwater reservoir characterization, and environmental remediation, adding valuable “hard” data capable of pinpointing the origins of fluids unambiguously.

Zhang, Yuran↗
DOE OSTI2022

Structural Characterization of a Toppling Rock Slab From Array-Based Ambient Vibration Measurements and Numerical Modal Analysis

Accurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array-based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross-correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500-m-long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Further, statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

58 GEOSCIENCES↗
DOE OSTI2023

ParFlow Sand Tank: A tool for groundwater exploration

The ParFlow Sand Tank model is an open source application designed to allow users to interactively simulate and visualize groundwater movement through the subsurface. The app is designed for both research and education; teaching hydrogeology concepts and making it easy explore and run sophisticated groundwater simulations. Our goal is to support increased accessibility and usability of research grade hydrology tools for research and teaching. The Sand Tank application simulates groundwater and surface water fluxes as well as contaminant transport in real time using the integrated physical hydrology model ParFlow (Kollet & Maxwell, 2006; Maxwell & Miller, 2005; Osei-Kuffuor et al., 2014) and the particle tracking code EcoSlim (Maxwell et al., 2019). ParFlow is a numerical hydrology model that simulates spatially distributed groundwater and surface water flow. It is a well established research tool with more than 90 publications documenting its development use to advance our understanding of groundwater dynamics and groundwater surface water interactions from the hillslope to the continental scale e.g. (Condon et al., 2020; Condon & Maxwell, 2019; Maxwell & Condon, 2016). It is designed for efficient parallel computation and has been run on many platforms spanning from laptops to supercomputers. However, one of the challenges of ParFlow is that it requires significant training and hydrologic expertise to develop simulations. The Sand Tank application makes this model accessible to anyone for education and exploration. Our application uses ParFlow for its simulation backend and ParaView for the data loading and processing. The communication infrastructure relies on the ParaViewWeb framework. We use model templates deployed in Docker images to setup the Sand Tank framework. Users can build the application locally or interact with it through our web deployment. When interacting with a template users can interactively change model parameters like subsurface processes or pump/inject water into the subsurface and watch the system respond to their changes in real time as the simulation runs. Additionally, our template setup will allow more advanced users to build custom templates of increasing complexity for both research and educational purposes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Revisiting Wagner’s Criteria to Predict Establishment and Retention of Alumina Scales on Ternary NiCrAl Alloys

Despite the advent of numerical modeling approaches and high-performance computing infrastructure, the design and development of corrosion-resistant high temperature alloys (> 500 °C) continue to be largely empirical and typically involve extensive experimentation. This is mainly due to the lack of a single unified physics-based model that can address the impact of multiple competing factors such as time, environment, alloy composition, microstructure, and geometry. The classical Wagner’s criteria have been foundational to estimate the minimum concentrations required of an oxide-forming element to establish and sustain a protective oxide scale. However, the formulation is primarily limited to lower-order alloy systems (binary alloys) and ignores the time dependence of subsurface compositional changes in the alloy. The lack of key data on the temperature and composition dependence of the solubility and transport of oxidants in multicomponent-multiphase alloys further exacerbates the problem. In the present work, a few of these limitations were addressed using a flux-based approach (FLAP) which tracks the spatiotemporal evolution of the fundamental flux balance between oxygen and the oxide-forming elements to enable the prediction of the formation of an external alumina scale in ternary NiCrAl alloys. In conclusion, the modeling results were validated with the literature findings and additional experimental work conducted in the present work.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Model Development for Thermal-Hydrology Simulations of a Full-Scale Heater Experiment in Opalinus Clay

Disposal of commercial spent nuclear fuel in a geologic repository is studied. In situ heater experiments in underground research laboratories provide a realistic representation of subsurface behavior under disposal conditions. Here, this study describes process model development and modeling analysis for a full-scale heater experiment in opalinus clay host rock. The results of thermal-hydrology simulation, solving coupled nonisothermal multiphase flow, and comparison with experimental data are presented. The modeling results closely match the experimental data.

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

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Radar observations at 3.5 and 12.6 cm wavelength of asteroid 433 Eros

A study of the asteroid 433 Eros using 3.5 and 12.6 cm radar waves indicates that the surface is very much rougher than any planetary or lunar surface observed by this method. A surface completely covered with sharp edges, pits, subsurface holes, or embedded chunks with scale sizes on the order of our wavelengths seems to be indicated. A model based on a rough rotating triaxial ellipsoid having radii in the rotation equator of 18.6 and 7.9 km agrees well with our data, although a strong wobble in the apparent center frequency of the spectra as rotation progresses indicates that one side may be more reflective than the other, or more likely, that the projected axis of rotation does not equally divide the projected area.

Jurgens, R. F.↗
NASA NTRS2003

Lunar and Planetary Science XXXV: Mars Impact Cratering

This document discusses the following topics: High Resolution Digital Elevation Models of Pristine Explosion Craters; Crater Degradation in the Martian Highlands: Morphometric Analysis of the Sinus Sabaeus Region and Simulation Modeling Suggest Fluvial Processes; Analysis of Impact Crater Preservation on Mars Using THEMIS Data; Atmospheric Entry Studies and the Smallest Impact Craters on Mars; Updating the Crater Count Chronology System for Mars; Control of Impact Crater-related Fracture Systems on the Subsurface Hydrology and Ground Collapse; Quantitative Analyses of Terrestrial Crater Deposits: Constraining Formation and Sediment Transport Processes on Mars; and Predicted Effects of Surface Processes on Martian Impact Crater Depth/Diameter Relationships

Stephen J Mackwell↗
NASA NTRS2008

Transient Thermal Model and Analysis of the Lunar Surface and Regolith for Cryogenic Fluid Storage

A transient thermal model of the lunar surface and regolith was developed along with analytical techniques which will be used to evaluate the storage of cryogenic fluids at equatorial and polar landing sites. The model can provide lunar surface and subsurface temperatures as a function of latitude and time throughout the lunar cycle and season. It also accounts for the presence of or lack of the undisturbed fluff layer on the lunar surface. The model was validated with Apollo 15 and Clementine data and shows good agreement with other analytical models.

Christie, Robert J.↗
DOE OSTI2020

Challenges in simulating ground interacting nuclear explosions

This paper summarizes recent above-ground nuclear explosion simulations as part of a broader effort to better characterize conditions within a fireball that may influence the chemical evolution of bomb materials and other materials entrained from the local explosion environment. A critical component of this work is validation against historic footage of atmospheric testing, requiring that we understand how the frequency-dependent sensitivity of the utilized film footage influences data captured in such images. We focus first on the early physics of a nuclear explosion in the atmosphere before discussing some of the technical challenges we seek to capture in late-time models that include more complex emplacement conditions and subsurface features. We discuss required physics packages (compressible hydrodynamics, radiation transport, as well as necessary ancillary tables such as equations of state (EOS) and opacities). Additionally, we note reasonable “shortcuts” one may make and their limitations, e.g., using ideal gas EOS, replacing spectrally resolved radiation with spectrally averaged radiation, and exchanging deterministic transport with diffusion. We then discuss an approach to achieving an equilibrated initial stress state for problems where buoyancy and subsurface lithostatic stress are important. Our methodology is presented in the context of LLNL’s ALE3D multiphysics code but may readily be implemented in other codes. In this paper, we start with a description of the challenges of NUDET simulations, followed by a presentation of the simulated intensity (flux) as it would appear on an analysis of the Dixie test. We then progressively introduce additional complexity in subsequent sections (near-surface burst and gravity initialization) before concluding.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
NASA NTRS1971

Thin highly conducting layer in the moon - Consistent interpretation of dayside and nightside electromagnetic responses.

The vacuum transient response of the moon to a time-varying spatially uniform magnetic field is determined for a lunar electrical conductivity model that was based on the harmonic analysis of Apollo 12 and Explorer 35 dayside magnetometer data. The transient response of the model is found to provide a plausible explanation of the behavior of the local vertical-surface magnetic field for an Apollo 12 magnetometer darkside transient event. A model containing a conducting core and a highly conducting thin subsurface layer is presented, and its transient behavior is discussed.

Schubert, G.↗
NASA NTRS1974

Reinterpretation of Mariner 9 IRIS data on the basis of a simulation of radiative-conductive convective transfer in the dust laden Martian atmosphere

Time dependent vertical distributions of atmospheric temperature and static stability were determined by a radiative-convective-conductive heat transfer model attuned to Mariner 9 IRIS radiance data. Of particular interest were conditions of both the dust-laden and dust-free atmosphere in the middle latitudes on Mars during the late S.H. summer season. The numerical model simulates at high spatial and temporal resolution (52 atmospheric and 30 subsurface levels; with a time-step of 7.5 min.) the heat transports in the ground-atmosphere system. The algorithm is based on the solution of the appropriate heating rate equation which includes radiative, molecular-conductive and convective heat transfer terms. Ground and atmosphere are coupled by an internal thermal boundary condition.

Pallman, A. J.↗