Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface stress”

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 37 records · Page 2

Geothermal Hydraulic Stimulation: Overview of Methods and Best Practices

This paper outlines important information related to the hydraulic stimulation of geothermal wells and considerations for regulators responsible for evaluating applications to conduct these operations. Hydraulic stimulation applied to geothermal wells is an evolving technology that is different from the hydraulic fracturing process widely used to complete unconventional oil and gas wells. Important differences relate to operating pressures, variations in rock and fracture processes, and characteristics and disposal of hydraulic fracture fluids. Hydraulic stimulation is not widely used in the geothermal industry; however, future development of enhanced geothermal systems (EGS) will require hydraulic stimulation in order to enhance and create fracture permeability to allow wells to deliver sufficient heat and fluid to power plants. Gigawatts of EGS potential have been identified (U.S. Department of Energy 2019), so it is expected that the use of hydraulic stimulation of geothermal wells will be more common in the future. At present, there are no formal regulations guiding drilling programs or sundry notices that propose hydraulic stimulation of geothermal wells. On federal lands, important constraints and oversight are embedded in the process for obtaining Geothermal Drilling Permits outlined in the Code of Federal Regulations (43 CFR Section 3262.11); in Geothermal Resource Order 2 (GRO 2) guidance for assessment and mitigation of impacts of geothermal operations as well as construction and testing of geothermal wells; and in a Bureau of Land Management (BLM) Induced Seismicity Instruction Memorandum (BLM 2018). Additionally, local knowledge from BLM field offices and expertise of BLM engineers, or that of state regulators for operations on private and state lands, will be applied to fully assess operator applications to conduct geothermal hydraulic stimulation activities. Significant research effort has been directed toward EGS, so geothermal hydraulic stimulation techniques, impacts, and results are evolving and becoming better understood, including adaptation of decades of oil and gas industry experience with hydraulic fracturing of unconventional reservoirs. Recent EGS activities in the United States provide details about hydraulic stimulation with respect to oversight, stimulation design, execution, and results. Based on experience to date, the important issues to address when hydraulically stimulating a geothermal well include the following: 1. Wellbore construction and integrity must be appropriate to protect groundwater and manage stimulation pressures. 2. Understanding of lithology, faults, fractures, and subsurface stress state is necessary to design stimulation plans and predict results. 3. Seismic monitoring allows for observation and mitigation of induced seismicity.

geothermal↗

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗

A Non-Invasive Approach for Elucidating the Spatial Distribution of In-Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Final Report)

The primary objectives of this project were to: develop a new, non-invasive method to define the orientation and magnitude of in-situ stresses in deep geologic formations (depths that exceed 1500 meters) considered for carbon dioxide (CO 2 ) sequestration; to demonstrate the method at a real field site; and, to extend current technical capabilities and reduce costs of acquiring information about deep-stress fields. This investigation was segmented into five tasks, including Project Management (Task 1) and four technical tasks (Tasks 2 through 5). The results of the technical tasks are described in this document.

02 PETROLEUM↗

Stress engineering for crack and dendrite prevention in solid electrolytes via ion implantation

Solid-state batteries represent a promising technology that offers safer and more densely packed energy storage. A primary cause of failure in solid-state cells is the penetration of metal dendrites through the solid electrolyte. Here, we report that fluorine-ion implantation can enhance the mechanical resistance of solid electrolyte Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 to dendrite propagation by inducing residual compressive stress in the subsurface of the electrolyte. Ion implantation modifies subsurface residual stress and also alters the electrolyte’s surface, resulting in multifunctional enhancement. The combined chemical and mechanical effects of ion implantation enable reversible lithium metal stripping and plating at elevated current densities while enhancing air stability and mitigating the formation of a harmful carbonate layer at the surface. This study provides new insights into a scalable dendrite-suppression strategy for designing solid-state batteries suitable for cycling at room temperature and low stack pressures.

25 ENERGY STORAGE↗

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↗

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗

Geothermal Deep Direct-Use for Low-Carbon Heating: A Case Study at Cornell University

Cornell University is developing a geothermal deep direct-use system to provide reliable, low carbon, and cost-competitive heating for its main campus in Ithaca, New York. This geothermal system, also referred to as Earth Source Heat, would provide the majority of the campus heating load and represents a key pillar for the university to become carbon neutral by 2035. Over the last 10 years, numerous desktop, seismic, and geophysical studies have been undertaken at Cornell to characterize the local subsurface, predict long-term reservoir performance, design the integration into the existing district heating network, and estimate capital and levelized costs. Recently, the university received funding for drilling an observatory well to obtain in-situ subsurface data (e.g., reservoir temperatures, pressures, fluid chemistry, stress field), significantly lowering subsurface uncertainty and de-risking the project. Drilling of the test well is planned for fall 2021. This paper provides an overview and status of the project, including a summary of past studies and a description of the proposed well and planned testing activities.

Cornell University↗

Modeling Validation of Mechanical and Thermal Borehole Breakouts From Polyaxial Laboratory Tests to Determine Maximum Horizontal in Situ Stress

To accurately assess rock behavior in the deep subsurface, it is necessary to measure the in situ stress directions and magnitudes. The current methods for measuring the in situ stress state in the deep subsurface primarily include hydraulic fracturing tests (i.e., the minimum horizontal stress) or occasionally observing existing compressive borehole breakouts (i.e., the maximum horizontal stress) that have occurred naturally from drilling. If there are no existing compressive breakouts in a borehole, the maximum horizontal in situ stress cannot be estimated with much confidence. In response to this data gap, a new thermal breakout technology is being developed that will provide a method for thermally inducing borehole breakouts and obtaining consistent measurements of the maximum horizontal stress magnitude. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which is directly correlated to the maximum horizontal stress magnitude. In support of developing the thermal breakout technology, polyaxial laboratory tests have been performed on small-scale boreholes within rock blocks where mechanically- and thermally-induced borehole breakouts have been created. Numerical models along with the principal of superposition were created and used to analyze the polyaxial laboratory tests to predict the maximum horizontal stress, given the same data that would be obtained in an actual subsurface borehole thermal breakout test. Multiple failure criteria were used to evaluate the best prediction of the breakout onset and maximum horizontal stress. The maximum horizontal stress predictions were compared to the actual maximum horizontal stress applied at breakout using acoustic events recorded from emission sensors in the polyaxial tests. The results showed consistent results that could be used to refine the modeling approach and failure criterion that are used to make the maximum horizontal stress predictions. This study provided insight and validation for the thermal breakout stress measurement concept.

Jones, Matt↗

Subsurface Cooling and Sea Surface Temperature Pattern Formation Over the Equatorial Pacific

Abstract The equatorial cold tongue region has not warmed up in response to historical radiative forcing in the real world, contrary to the strong warming often simulated by climate models. Here we demonstrate that climate models fail to represent one or both of the key processes driving observed sea surface temperature (SST) pattern formation: a realistic surface wind stress pattern shaping subsurface cooling through wind‐driven circulation changes, and effective connectivity between subsurface and surface temperatures via upwelling and mixing. Consequently, none of the models approximate the observed lack of cold tongue SST warming and strengthening of zonal SST gradient across the equatorial Pacific. Furthermore, those that come closest achieve this due to interhemispheric warming differences rather than equatorial dynamics as observed. Addressing different origins of subsurface cooling in observations and simulations, and how they connect to SST, will lead to improved understanding of tropical Pacific SST changes to date and how they will evolve in the future.

Jiang, Feng [Lamont‐Doherty Earth Observatory Colu↗

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↗

Refining Principal Stress Measurements in Reservoir Underburden in Regions of Induced Seismicity through Seismological Tools, Laboratory Experiments - Final Technical Report

This project developed methodologies to measure the in-situ principal stress in the deep subsurface through use of multiple independent, but complementary, seismic methods, laboratory verification, and development of theoretical frameworks. By leveraging existing regional and local datasets we developed, tested, and refined a set of diagnostic tools for determining the in-situ stress state with reduced uncertainty at and below reservoir depths (1.5-6 km). A set of novel tools was produced that are scale independent, such that their utility is equivalent on regional, field scale, and near borehole monitoring of principal stresses in reservoir underburden for carbon storage projects. During a 4-year Department of Energy (DOE) and Southern Company funded project, carried out by the Electric Power Research Institute (EPRI), Lawrence Livermore National Laboratory (LLNL), the University of Oklahoma (OU), and the U.S. Geological Survey (USGS), the project team developed methodologies to measure the far-field in-situ principal stress in the deep subsurface, leveraging induced seismicity data from waste-water disposal projects. These methodologies consisted in the use of well-established and technically advanced seismic processing methods, such as virtual seismometer method-moment tensor (VSM-MT) and shear wave splitting (SWS), that are adept at recovering the stress orientation and certain components of the stress tensor. These methods were applied to robust seismicity catalogs created with matched filter techniques near sites of active fluid disposal—a proxy for carbon storage sites where such datasets are more limited. Estimates of the stress orientation made with seismic processing tools were considered along with laboratory acoustic emission experiments conducted on rock samples from the region of interest. Stress orientations in the studied region do not vary significantly across distances of ~100 km, nor are they found to rotate through time as a consequence of local wastewater disposal, as previously speculated. Finally, the project team investigated the trade-offs among the different seismic methods and evaluated the range of uncertainty that is generated with these methodologies, which led to a practical use and refinement of the VSM-MT technique when it is applied to field datasets. Understanding the trade-offs between these different methods highlighted the potential benefits of improved quantification of uncertainties on stress field estimations.

58 GEOSCIENCES↗

Historical Subsurface Cooling in the Tropical Pacific and Its Dynamics

Abstract Understanding how the tropical Pacific responds to rising greenhouse gases in recent decades is of paramount importance given its central role in global climate systems. Extensive research has explored the long-term trends of tropical Pacific sea surface temperatures (SSTs) and the overlying atmosphere, yet the historical change in the upper ocean has received far less attention. Here, we present compelling evidence of a prominent subsurface cooling pattern along the thermocline in the central-to-eastern tropical Pacific since 1958. This subsurface cooling has been argued to be contributing to the observed cooling or lack of warming of the equatorial cold tongue SST. We further demonstrate that different mechanisms are responsible for different parts of the subsurface cooling. In the central-to-eastern equatorial Pacific and the southeastern off-equatorial Pacific, where zonal wind stress strengthens, a pronounced subsurface cooling trend emerges just above the thermocline that is closely tied to increased Ekman pumping. In the eastern equatorial Pacific where zonal wind stress weakens, the westward surface current and eastward Equatorial Undercurrent weaken as well, resulting in reduced vertical current shear and increased ocean stability, which suppresses vertical mixing and leads to local cooling. We conclude that the historical subsurface cooling is primarily linked to dynamical adjustments of ocean currents to tropical surface wind stress changes.

Jiang, Feng↗

Initial Field Testing in the Deep Subsurface for the Thermal Breakout Project for Measuring In Situ Stress

The current state-of-the-art technology for in situ stress measurements involves an integrated approach that combines borehole breakout observations, drilling-induced tensile fractures, and hydraulic fracturing tests (i.e., mini-fracs); however, this methodology has several limitations that often prevent successful in situ stress measurements. One major limitation is that breakouts do not appear in every borehole and are generally a natural occurrence that cannot easily be controlled. Because borehole breakouts are used to measure the maximum horizontal in situ stress magnitude, the absence of borehole breakouts presents a major data gap for in situ stress measurements. As a response to this data gap a new thermal breakout technology is being developed that will provide a method for thermally inducing borehole breakouts and obtaining consistent measurements of the maximum horizontal stress magnitude. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which is directly correlated to the maximum horizontal stress magnitude. An initial component of the thermal breakout project involves developing a prototype downhole tool and performing field testing in the deep subsurface. The primary objective of the initial tool development and field testing is to provide a physical proof of concept for the thermal breakout technology. Three thermal breakout field tests were performed at approximately 4850 feet deep in the Sanford Underground Research Facility. Each test was performed in a borehole that had been previously used for overcoring stress measurements. The third and final test successfully created two diametrically opposed breakouts after deploying the thermal breakout tool for approximately 3 hours of continuous heating. Post-analysis of the thermally induced breakouts confirmed that the thermal breakout orientation directly corresponds to the known in situ stress orientation and magnitude.

Voegeli, Sam↗

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↗