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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

The interpretation of crustal dynamics data in terms of plate interactions and active tectonics of the Anatolian Plate and surrounding regions in the Middle East

The long-term objective of this project is to interpret NASA's Crustal Dynamics measurements (SLR) in the Eastern Mediterranean region in terms of relative plate movements and intraplate deformation. The approach is to combine realistic modeling studies with analysis of available geophysical and geological observations to provide a framework for interpreting NASA's measurements. This semi-annual report concentrates on recent results regarding the tectonics of Anatolia and surrounding regions from ground based observations. Also reported on briefly is progress in the use of the Global Positioning System to densify SLR observations in the Eastern Mediterranean. Reference is made to the previous annual report for a discussion of modeling results.

Toksoz, M. Nafi↗
DOE OSTI2023

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗
DOE OSTI2023

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗
DOE OSTI2024

Enhancing Regional Seismic Velocity Models With Higher‐Resolution Local Results Using Sparse Dictionary Learning

Abstract We use sparse dictionary learning to develop transformations between seismic velocity models of different resolution and spatial extent. Starting with data in the common region of both models, the method can enhance a regional lower‐resolution model to match the style and resolution of local higher‐resolution results while preserving its regional coverage. The method is demonstrated by applying it to two‐dimensional V S and three‐dimensional V P and V S regional and local velocity models in southern California. The enhanced reconstructed regional results exhibit clear visual improvements, especially in the reconstructed V P / V S ratios, and better correlations with geological features. Moreover, the reconstructed regional V P , V S models outperform the original ones in comparison of simulated earthquake waveforms to observations. The improved fitting to observed waveforms extends beyond the domain of the overlapping region. The developed dictionary learning approach provides physically interpretable results and offers a powerful tool for additional applications of data enhancement in earth sciences.

Geochemistry & Geophysics↗
NASA NTRS1997

Hydrothermal hydration of Martian crust: illustration via geochemical model calculations

If hydrothermal Systems existed on Mars, hydration of crustal rocks may have had the potential to affect the water budget of the planet. We have conducted geochemical model calculations to investigate the relative roles of host rock composition, temperature, water-to-rock ratio, and initial fluid oxygen fugacity on the mineralogy of hydrothermal alteration assemblages, as well as the effectiveness of alteration to store water in the crust as hydrous minerals. In order to place calculations for Mars in perspective, models of hydrothermal alteration of three genetically related Icelandic volcanics (a basalt, andesite, and rhyolite) are presented, together with results for compositions based on SNC meteorite samples (Shergotty and Chassigny). Temperatures from 150 degrees C to 250 degrees C, water-to-rock ratios from 0.1 to 1000, and two initial fluid oxygen fugacities are considered in the models. Model results for water-to-rock ratios less than 10 are emphasized because they are likely to be more applicable to Mars. In accord with studies of low-grade alteration of terrestrial rocks, we find that the major controls on hydrous mineral production are host rock composition and temperature. Over the range of conditions considered, the alteration of Shergotty shows the greatest potential for storing water as hydrous minerals, and the alteration of Icelandic rhyolite has the lowest potential.

NASA Discipline Exobiology↗
DOE OSTI2020

Carbon Limitation Leads to Thermodynamic Regulation of Aerobic Metabolism

Organic matter (OM) metabolism in freshwater ecosystems is a critical source of uncertainty in global biogeochemical cycles, yet the processes regulating aerobic respiration in aquatic environments remain poorly understood. Aerobic respiration is typically predicted through kinetic controls such as organic carbon and oxygen concentrations, while assuming thermodynamic regulation has no influence. The role of OM thermodynamics has been mainly explored under anaerobic conditions; however, recent observations have suggested that thermodynamic properties of OM may influence aerobic respiration rates. Here, we explicitly test OM thermodynamics as a key regulator of aerobic respiration. We pair controlled microcosm experiments with ultrahigh-resolution OM characterization to show that aerobic respiration increases with OM thermodynamic favorability under carbon limitation. We also demonstrate a shift in the regulation of aerobic respiration from OM thermodynamics to nitrogen content when carbon is in excess, highlighting a central role for OM thermodynamics particularly in carbon-limited ecosystems. This illuminates a structural gap in aquatic biogeochemical models by demonstrating that both kinetic and thermodynamic constraints can modulate aerobic respiration. Our work therefore proposes a new conceptual model in which thermodynamic and nutrient limitations dually control aerobic respiration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

GDSA Framework Development and Process Model Integration FY2023

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is disposal system modeling (Sassani et al. 2021). The SFWST Geologic Disposal Safety Assessment (GDSA) work package is charged with developing a disposal system modeling and analysis capability for evaluating generic disposal system performance for nuclear waste in geologic media.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1994

The chemical reactivity of the Martian soil and implications for future missions

Possible interpretations of the results of the Viking Biology Experiments suggest that greater than 1 ppm of a thermally labile oxidant, perhaps H2O2, and about 10 ppm of a thermally stable oxidant are present in the martian soil. We reexamine these results and discuss implications for future missions, the search for organics on Mars, and the possible health and engineering effects for human exploration. We conclude that further characterization of the reactivity of the martian regolith materials is warrented-although if our present understanding is correct the oxidant does not pose a hazard to humans. There are difficulties in explaining the reactivity of the Martian soil by oxidants. Most bulk phase compounds that are capable of oxidizing H2O to O2 per the Gas Exchange Experiment (GEx) are thermally labile or unstable against reduction by atmospheric CO2. Models invoking trapped O2 or peroxynitrates (NOO2(-)) require an unlikely geologic history for the Viking Lander 2 site. Most suggested oxidants, including H2O2, are expected to decompose rapidly under martian UV. Nonetheless, we conclude that the best model for the martian soil contains oxidants produced by heterogeneous chemical reactions with a photochemically produced atmospheric oxidant. The GEx results may be due to catalytic decomposition of an unstable oxidizing material by H2O. We show that interfacial reaction sites covering less than 1% of the available soil surfaces could explain the Viking Biology Experiments results.

Zent, Aaron P.↗
NASA NTRS2001

A Model for Formation of Dust, Soil and Rock Coatings on Mars: Physical and Chemical Processes on the Martian Surface

This model is one of many possible scenarios to explain the generation of the current surface material on Mars using chemical, magnetic and spectroscopic data From Mars and geologic analogs from terrestrial sites. One basic premise of this model is that the dust/soil units are not derived exclusively from local rocks, but are rather a product of global, and possibly remote, weathering processes. Another assumption in this model is that there are physical and chemical interactions of the atmospheric dust particles and that these two processes create distinctly different results on the surface. Physical processes distribute dust particles on rocks and drift units, forming physically-aggregated layers; these are reversible processes. Chemical reactions of the dust/soil particles create alteration rinds on rock surfaces and cohesive, crusted surface units between rocks, both of which are relatively permanent materials. According to this model the dominant components of the dust/soil particles are derived from alteration of volcanic ash and tephra, and contain primarily nanophase and poorly crystalline ferric oxides/oxyhydroxide phases as well as silicates. These phases are the alteration products that formed in a low moisture environment. These dust/soil particles also contain a smaller amount of material that was exposed to more water and contains crystalline ferric oxides/oxyhydroxides, sulfates and clay silicates. These components could have formed through hydrothermal alteration at steam vents or fumeroles, thermal fluids, or through evaporite deposits. Wet/dry cycling experiments are presented here on mixtures containing poorly crystalline and crystalline ferric oxides/oxyhydroxides, sulfates and silicates that range in size from nanophase to 1-2 pm diameter particles. Cemented products of these soil mixtures are formed in these experiments and variation in the surface texture was observed for samples containing smectites, non-hydrated silicates or sulfates. Reflectance spectra were measured of the initial particulate mixtures, the cemented products and ground versions of the cemented material. The spectral contrast in the visible/near-infrared and mid-infrared regions is significantly reduced for the cemented material compared to the initial soil, and somewhat reduced for the ground, cemented soil compared to the initial soil. The results of this study suggest that diurnal and seasonal cycling on Mars will have a profound effect on the texture and spectral properties of the dust/soil particles on the surface. The model developed in this study provides an explanation for the generation of cemented or crusted soil units and rock coatings on Mars and may explain albedo variations on the surface observed near large rocks or crater rims.

Bishop, Janice↗
NASA NTRS2000

The faint young Sun paradox: an observational test of an alternative solar model

We report the results of deep observations at radio (3.6 cm) wavelengths of the nearby solar-type star pi 01 Ursa Majoris with the Very Large Array (VLA) intended to test an alternative theory of solar luminosity evolution. The standard model predicts a solar luminosity only 75% of the present value and surface temperatures below freezing on Earth and Mars at 4 Ga, seemingly in conflict with geologic evidence for liquid water on these planets. An alternative model invokes a compensatory mass loss through a declining solar wind that results in a more consistent early luminosity. The free-free emission from an enhanced wind around nearby young Sun-like stars should be detectable at microwave frequencies. Our observations of pi 01 UMa, a 300 million year-old solar-mass star, place an upper limit on the mass loss rate of 4-5 x 10(-11) M(solar) yr-1. Total mass loss from such a star over 4 Gyr would be less than 6%. If this star is indeed an analog of the early Sun, it casts doubt on the alternative model as a solution to the faint young Sun paradox, particularly for Mars.

NASA Discipline Exobiology↗
DOE OSTI2022

Physically rigorous reduced-order flow models of fractured subsurface environments without explosive computational cost

Fractured media models comprise discontinuities of multiple lengths (e.g. fracture lengths and apertures, wellbore area) that fall into the relatively insignificant length scales spanning millimeter-scale fractures to centimeter-scale wellbores in comparison to the extensions of the field of interest, and challenge the conventional discretization methods imposing highly-fine meshing and formidably large numerical cost. By utilizing the recent developments in the finite element analysis of electromagnetics that allow to represent material properties on a hierarchical geometry, this project develops computational capabilities to model fluid flow, heat conduction, transport and induced polarization in large-scale geologic environments that possess geometrically-complex fractures and man-made infrastructures without explosive computational cost. The computational efficiency and robustness of this multi-physics modeling tool are demonstrated by considering various highly-realistic complex geologic environments that are common in many energy and national security related engineering problems.

42 ENGINEERING↗
DOE OSTI2022

Gas diffusion through variably-water-saturated zeolitic tuff: Implications for transport following a subsurface nuclear event

Noble gas transport through geologic media has important applications in the characterization of underground nuclear explosions (UNEs). Without accurate transport models, it is nearly impossible to distinguish between xenon signatures originating from civilian nuclear facilities and UNEs. Understanding xenon transport time through the earth is a key parameter for interpreting measured xenon isotopic ratios. One of the most challenging aspects of modeling gas transport time is accounting for the effect of variable water saturation of geological media. In this study, we utilize bench-scale laboratory experiments to characterize the diffusion of krypton, xenon, and sulfur hexafluoride (SF6) through intact zeolitic tuff under different saturations. Here, we demonstrate that the water in rock cores with low partial saturation dramatically affects xenon transport time compared to that of krypton and SF6 by blocking sites in zeolitic tuff that preferentially adsorb xenon. This leads to breakthrough trends that are strongly influenced by the degree of the rock saturation. Xenon is especially susceptible to this phenomenon, a finding that is crucial to incorporate in subsurface gas transport models used for nuclear event identification. We also find that the breakthrough of SF6 diverges significantly from that of noble gases within our system. When developing field scale models, it is important to understand how the behavior of xenon deviates from chemical tracers used in the field, such as SF6 (Carrigan et al., 1996). These new insights demonstrate the critical need to consider the interplay between rock saturation and fission product sorption during transport modeling, and the importance of evaluating specific interactions between geomedia and gases of interest, which may differ from geomedia interactions with chemical tracers.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
NASA NTRS1994

Sediment oxygen profiles in a super-oxygenated antarctic lake

Perennially ice-covered lakes are found in the McMurdo Dry Valleys of southern Victoria Land, Antarctica. In contrast to temperate lakes that have diurnal photic periods, antarctic (and arctic) lakes have a yearly photic period. An unusual feature of the antarctic lakes is the occurrence of O2 at supersaturated levels in certain portions of the water column. Here we report the first sediment O2 profiles obtained using a microelectrode from a perennially ice-covered antarctic lake. Sediment cores collected in January and October 1987 from Lake Hoare in Taylor Valley show oxygenation down to 15, and in some cases, 25 cm. The oxygenation of sediments several centimeters below the sediment-water interface is atypical for lake sediments and may be characteristic of perennially ice-covered lakes. There is a significant difference between the observed January and October sediment O2 profiles. Several explanations may account for the difference, including seasonality. A time-dependent model is presented which tests the feasibility of a seasonal cycle resulting from the long photoperiod and benthic primary production in sediments overlain by a highly oxygenated water column.

NASA Center ARC↗
DOE OSTI2021

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗
DOE OSTI2022

Task 5 Field Scale Stress Modeling Scripts and Data

Project data and Matlab scripts for Task 5 of the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage”. Project includes the static earth models, grids and other data used for simulation and calibration of the two case studies addressed in the project, referred to herein as the “FutureGen2” model (Illinois) and the “Perch” model (Otsego county, Michigan). Large Dataset, please contact EDXSupport@netl.doe.gov

FutureGen2↗
DOE OSTI2022

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗