Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1994

Hotspot evolution and Venusian tectonic style

Because hotspots represent an important manifestation of heat loss on Venus, their geological evolution is of fundamental importance for any attempt to understand Venusian tectonics. Eistla Regio is an approximately 7500-km-long, moderately elevated region inferred to overlie one or more large mantle upwellings or hotspots. It also contains many shield volcanoes and coronae believed due to the rise of thermal plumes in the mantle. Central Eistla Regio includes two large volcanoes, Sappho and Anala, and several coronae in close proximity. Detailed mapping in this region results in two conclusions of tectonic significance: (1) Sappho and Anala occur near the intersection of two major extensional deformation zones, and (2) the coronae are older than the large volcanoes. Several of the coronae occur as a chain along Guor Linea, one of the major extensional deformation zones. Stratigraphic relationships indicate that the coronae began forming very soon after the emplacement of the widespread regional plains materials. Thus Central Eistla Regio was the site of a swarm of plumes that first formed coronae and then later formed shield volcanoes. The expected result of such a swarm would be thermal thinning of the elastic lithosphere with time. However, model results, geological observations, and gravity data suggest that the change from coronae to shield volcanoes was accompanied by a thickening of the lithosphere with time. This thickening is interpreted to be the result of global cooling of the lithosphere following the most recent episode of near-global resurfacing. The global cooling must have occurred faster than local heating of the lithosphere due to the impingement of thermal plumes.

Mcgill, George E.↗
DOE OSTI2025

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

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

Geothermal Direct-Use Applications for the District Energy System in Bucharest, Romania

The city of Bucharest, Romania, hosts the second-largest district energy system (DES) in the world. Geothermal resources can be considered as a supplementary heat source to support the demand for domestic hot water and space heating in the winter and shoulder seasons. The National Laboratory of the Rockies (NLR) has conducted a study that considers geothermal energy to serve a fraction of the existing district heating network operated by Electrocentrale Bucure?ti (ELCEN), the utility operating the DES. Lower Cretaceous and Jurassic limestones make up the main geothermal aquifer underlying Bucharest, which hosts temperatures suitable for district heating (up to 90 degrees C to the north of the city). Anomalous geothermal gradients have been observed to the north of the city, where a pumped well has produced 82 degrees C brine at the wellhead to feed the Therme Bucharest Spa. An anomalous gradient has also been reported to the southeast of the city (35 degrees C/km). NLR modeled the building heating loads of a small portion of the DES (a block of nine prototypical buildings) in its Urban Renewable Building and Neighborhood Optimization (URBANopt ) platform. To simulate meeting a baseload benchmark of 20 MWth deliverable to a small portion of the DES, the NLR team used GEOPHIRES to model production scenarios for (1) hydrothermal systems coupled with heat pumps targeting the main geothermal aquifer in the north, (2) enhanced geothermal systems targeting hot dry rock in the southeast, and (3) huff-and-puff systems targeting a gradient of 25 degrees C/km. Finally, NLR conducted a high-level sensitivity study around the techno-economics of these systems. The outcomes of this work indicate that hot dry rock geothermal resources that can deliver at least 90 degrees C hot water to the Geothermal District Energy System (GeoDES) offer a possible solution for supplemental geothermal heat delivered to the existing DES.

15 GEOTHERMAL ENERGY↗
NASA NTRS2001

Impact Cratering Calculations

We examined the von Mises and Mohr-Coulomb strength models with and without damage effects and developed a model for dilatancy. The models and results are given in O'Keefe et al. We found that by incorporating damage into the models that we could in a single integrated impact calculation, starting with the bolide in the atmosphere produce final crater profiles having the major features found in the field measurements. These features included a central uplift, an inner ring, circular terracing and faulting. This was accomplished with undamaged surface strengths of approximately 0.1 GPa and at depth strengths of approximately 1.0 GPa. We modeled the damage in geologic materials using a phenomenological approach, which coupled the Johnson-Cook damage model with the CTH code geologic strength model. The objective here was not to determine the distribution of fragment sizes, but rather to determine the effect of brecciated and comminuted material on the crater evolution, fault production, ejecta distribution, and final crater morphology.

Ahrens, Thomas J.↗
DOE OSTI2022

A Bayesian Approach for In-Situ Stress Prediction and Uncertainty Quantification for Subsurface Engineering

Many subsurface engineering applications require accurate knowledge of the in-situ state of stress for their safe design and operation. Existing methods to meet this need primarily include field measurements for estimating one or more of the principal stresses from a borehole, or optimization methods for constructing a 3D geomechanical model in terms of geophysical measurements. These methods, however, often contain considerable uncertainty in estimating the state of stress. Here, in this paper, we build on a Bayesian approach to quantify uncertainty in stress estimations for subsurface engineering applications. This approach can provide an estimate of the 3D distribution of stress throughout the volume of interest and provide an estimate of the uncertainty arising from the stress measurement, the rheology parameters, and a paucity of measurements. The value of this approach is demonstrated using stress measurements from the In Salah carbon storage site, which was one of the world’s first industrial carbon capture and storage projects. This demonstration shows the application of this Bayesian approach for estimating the initial state of stress for In Salah and quantifying the uncertainty in the estimated stress. Also, an assessment of a maximum injection pressure to prevent geomechanical risks from CO 2 injection pressures is provided in terms of the probability distribution of the minimum principal stress quantified by the approach. With the In Salah case study, this paper demonstrates that using the Bayesian approach can provide additional insights for site explorations and/or project operations to make informed-site decisions for subsurface engineering applications.

58 GEOSCIENCES↗
DOE OSTI2022

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗
DOE OSTI2021

Can artificial intelligence and data-driven machine learning models match or even replace process-driven hydrologic models for streamflow simulation?: A case study of four watersheds with different hydro-climatic regions across the CONUS

With recent developments in computational techniques, Data-driven Machine Learning Models (DMLs) have shown great potential in simulating streamflow and capturing the rainfall-runoff relationship in given watersheds, which are traditionally fulfilled by Process-based Hydrologic Models (PHMs). There are debates on whether the DMLs can outperform and possibly replace the classical PHMs for streamflow simulation and river forecasting, but no clear conclusions have been made. This study aims to investigate whether the newer DMLs have any potential in further improving the simulation accuracy of classical PHMs, and vice versa. To do this, we compared a few popular PHMs and DMLs over four watersheds across the Continental US (CONUS) that are associated with different input, climate, and regional conditions. A total of five hydrologic models were chosen, including (1) two classical lumped models, i.e., the Sacramento Soil Moisture Accounting (SAC-SMA) and Xinanjiang (XAJ); (2) one modern distributed model, termed Coupled Routing and Excess Storage (CREST); (3) and two DMLs including an Artificial Neural Networks (ANN) and a deep learning model, termed Long Short Term Memory (LSTM). Our results demonstrated that the DMLs still significantly biased when using the baseline input scenario with the PHMs. However, the DMLs fed with delayed input scenarios had great potential and can reach high simulation accuracy. The DMLs, especially the ANN, outperformed other employed models under the rainfall-runoff relationship in which rainfall dominantly drives. Furthermore, the DMLs also showed better performance in the high-flow regime, while the PHMs had a better performance for the low-flow regime, implying both PHMs and DMLs have their own merits and are worthy of joint development. In general, our study indicated a great potential of using DMLs to simulate streamflow, but further studies are still needed to verify the transferability and scalability of DMLs in large-scale experiments, such as the Distributed Model Intercomparison Projects 1&2 conducted by National Weather Services but to compare modern DMLs and PHMs.

58 GEOSCIENCES↗
DOE OSTI2025

Understanding the Hydronium Cation in the Solid State: A Study in Synthetic Hydronium Uranyl Phosphate and Arsenate Mineral Systems and Their Irradiation Stability

Here, we present neutron diffraction, electronic structure calculations, and optical vibrational spectroscopic characterization of hydronium uranyl phosphate (HUP), the synthetic analog of chernikovite ((H 3 O)(UO 2 )(PO 4 )·3H 2 O) to gain insight into structural features of the hydronium cation in the solid state. HUP crystallizes in P2 1 /c instead of previously reported P4/ncc 1 at room temperature and Pccn below 302 K. Lower symmetry was required due to hydronium and interstitial water positions within the structure. Crystallographic positions of hydronium cations determined from powder neutron diffraction data are bolstered by Raman and attenuated total reflectance infrared spectroscopic measurements and are further informed by density functional theory with phonon eigenvector analysis for spectral assignments. Finally, HUP and its arsenate analog (trogerite, (H 3 O)(UO 2 )(AsO 4 )·3H 2 O), (HUAs) were studied using He 2+ irradiation as an analog for α radiolysis to investigate the irradiation stability of these phases and the stability of hydronium cations in the solid state. Dose studies were employed wherein each sample was irradiated to 5, 10, 15, 25, and 50 MGy. Structural insight regarding irradiated materials is gained using optical vibrational spectroscopy and powder X-ray diffraction. Surprising irradiation stability of HUAs was found up to 50 MGy of dose, which may have important implications for understanding and modeling the geologic stability of legacy U.

actinides↗
DOE OSTI2024

Ocean cavity regime shift reversed West Antarctic grounding line retreat in the late Holocene

Recent geologic and modeled evidence suggests that the grounding line of the Siple Coast of the West Antarctic Ice Sheet (WAIS) retreated hundreds of kilometers beyond its present position in the middle to late Holocene and readvanced within the past 1.7 ka. This grounding line reversal has been attributed to both changing rates of isostatic rebound and regional climate change. Here, we test these two hypotheses using a proxy-informed ensemble of ice sheet model simulations with varying ocean thermal forcing, global glacioisostatic adjustment (GIA) model simulations, and coupled ice sheet-GIA simulations that consider the interactions between these processes. Our results indicate that a warm to cold ocean cavity regime shift is the most likely cause of this grounding line reversal, but that GIA influences the rate of ice sheet response to oceanic changes. This implies that the grounding line here is sensitive to future changes in sub-ice shelf ocean circulation.

58 GEOSCIENCES↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
DOE OSTI2024

LYNM-PE1 Seismic Parameters from Borehole Log, Laboratory, and Tabletop Measurements

The goal of this work is to provide a database of quality-checked seismic parameters that can be integrated with the Geologic Framework Model (GFM) for the LYNM-PE1 (Low Yield Nuclear Monitoring – Physical Experiment 1) testbed. We integrated data from geophysical borehole logs, tabletop measurements on collected core, and laboratory measurements. We reviewed for internal consistency among each measurement type, documented the caveats of measurement conditions, and integrated lithologic logs to check the validity of outlier values. The resulting consolidated parameter tables can be used as inputs for modeling and analysis codes and are designed to interface with the GFM, which is being actively developed.

58 GEOSCIENCES↗
DOE OSTI2026

Improving Production in the Emerging Paradox Oil Play (Final Technical Report)

This report documents the results of the DOE-funded project Improving Production in the Emerging Paradox Oil Play, focused on integrating geological characterization, petrophysical analysis, reservoir simulation, and field development planning to improve hydrocarbon recovery from the Paradox Basin. Project activities included core and log analysis, assessment of reservoir heterogeneity and compartmentalization, evaluation of completion and stimulation strategies, and development of an uncertainty-aware field development framework. Findings highlight the importance of selective completions, geology-driven well placement, and targeted stimulation in vertically heterogeneous carbonate reservoirs.

02 PETROLEUM↗
NASA NTRS1982

An inventory of undiscovered Canadian mineral resources

Unit regional value (URV) and unit regional weight are area standardized measures of the expected value and quantity, respectively, of the mineral resources of a region. Estimation and manipulation of the URV statistic is the basis of an approach to mineral resource evaluation. Estimates of the kind and value of exploitable mineral resources yet to be discovered in the provinces of Canada are used as an illustration of the procedure. The URV statistic is set within a previously developed model wherein geology, as measured by point counting geologic maps, is related to the historical record of mineral resource production of well-developed regions of the world, such as the 50 states of the U.S.A.; these may be considered the training set. The Canadian provinces are related to this training set using geological information obtained in the same way from geologic maps of the provinces. The desired predictions of yet to be discovered mineral resources in the Canadian provinces arise as a consequence. The implicit assumption is that regions of similar geology, if equally well developed, will produce similar weights and values of mineral resources.

Labovitz, M. L.↗
NASA NTRS1984

Ice-covered volcanic water flows on Ganymede

The most geologically reasonable model for the origin of grooved terrain on Ganymede involves flooding of a planetary graben system by water or water-ice magmas that rose to the surface along normal faults in the rift zones. To investigate this process further, the thermal evolution of a water flow on Ganymede's surface was examined. The model is a variation of earlier work on ice-covered water flows. It is based on the simplifying assumptions that laminar flow and a concomitant solid ice cover are achieved relatively soon after eruption; and that the amount of energy lost to the substratum is negligible.

Allison, M. L.↗