Search NASA⌕ Search

SEARCH · Search NASA

Results for “Working Fluid demonstrates”

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 505 records · Page 28

JSC Materials Laboratory Reproduction and Failure Analysis of Cracked Orbiter Reaction Control System Niobium Thruster Injectors

In April 2004 a Space Shuttle Orbiter Reaction Control System (RCS) thruster was found to be cracked while undergoing a nozzle (niobium/C103 alloy) retrofit. As a failure resulting from an in-flight RCS thruster burn-through (initiated from a crack) could be catastrophic, an official Space Shuttle Program flight constraint was issued until flight safety could be adequately demonstrated. This paper describes the laboratory test program which was undertaken to reproduce the cracking in order to fully understand and bound the driving environments. The associated rationale developed to justify continued safe flight of the Orbiter RCS system is also described. The laboratory testing successfully reproduced the niobium cracking, and established specific bounding conditions necessary to cause cracking in the C103 thruster injectors. Each of the following conditions is necessary in combination together: 1) a mechanically disturbed / cold-worked free surface, 2) an externally applied sustained tensile stress near yield strength, 3) presence of fluorine-containing fluids on exposed tensile / cold-worked free surfaces, and 4) sustained exposure to temperatures greater than 400 F. As a result of this work, it was concluded that fluorine-containing materials (e.g. HF acid, Krytox , Brayco etc.) should be carefully controlled or altogether eliminated during processing of niobium and its alloys.

Castner, Willard L.↗

Spin-on deposition of amorphous zeolitic imidazolate framework films for lithography applications

Amorphous zeolitic imidazolate framework (aZIF) films have been recently introduced as resists for electron beam and extreme ultraviolet lithography. aZIFs are also being considered for separation applications, including thin film membranes. However, the reported methods for aZIF deposition are currently based on highly empirical trial-and-error approaches that hinder control of film composition, thickness and uniformity as well as scale-up and transferability to different coating geometries. This work presents a method for depositing aZIF films with controllable thickness using dilute precursors mixed immediately before encountering the substrate. Importantly, the method is amenable to quantitative analysis by computational fluid dynamics to extract intrinsic deposition rates and limiting reactant transport diffusivities, enabling predictive physics-based modeling of the deposition process. This allows the deposition method to be adapted for spin coating on silicon wafers to prepare high-quality aZIF films with consistently controlled thickness. Using this approach, high-resolution resist performance and wafer-scale use for beyond extreme-ultraviolet lithography of aZIF films is demonstrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ultralight Fabric Reflux Tube (UFRT) Thermal/Vacuum Test

Spacecraft thermal control systems are essential to provide the necessary environment for the crew and equipment to function adequately on space missions. The Ultralight Fabric Reflux Tube (UFRT) was developed by Pacific Northwest Laboratory (PNL) as a lightweight radiator concept to be used on planetary-type missions (e.g., Moon, Mars). The UFRT consists of a thin-walled tube (acting as the fluid boundary), overwrapped with a low-mass ceramic fabric (acting as the primary pressure boundary). The tubes are placed in an array in the vertical position with the evaporators at the lower end. Heat is added to the evaporators, which vaporizes the working fluid. The vapor travels to the condenser end above and cools as heat is radiated to the environment. The fluid condensed on the tube wall is then returned to the evaporator by gravity. The primary objectives for the fiscal year 1994 program included the design and fabrication of prototype UFRTs and thermal/vacuum chamber testing of these test articles. Six UFRTS, with improved titanium liners, were successfully manufactured and provided to the Johnson Space Center in July 1994. Five were tested in a thermal/vacuum chamber in September 1994. Data obtained to characterize the performance of the UFRTs under simulated lunar conditions demonstrated the design concept successfully. In addition, a trade study showed that an optimized/improved UFRT could achieve as much as a 25% mass savings in the heat rejection subsystem of future planetary-type thermal control systems.

Hurlbert, K. M.↗

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗

2D kinetic-ion simulations of inverted corona fusion targets

Laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. The scheme consists of a hollow or gas-filled spherical shell made of deuterated plastic. The shell has one or more laser entrance holes (LEH), resembling a spherical hohlraum. The laser passes through the LEH’s and illuminates the interior surface of the shell, ablating a plasma that travels inward towards the target center. Long ion mean free paths in the converging plasma can lead to significant interpenetration, atomic mix, and other kinetic effects. Here, in this work we report on numerical simulations of inverted corona targets using the kinetic-ion, fluid–electron hybrid particle-in-cell (PIC) approach in 2D RZ geometry. 2D simulations suggest that shape effects do not have a significant impact on plasma evolution and observed yield trends are primarily the result of 1D kinetic mix mechanisms. Simulations are also compared against available experimental data recorded at the OMEGA laser facility. In particular, synthetic x-ray emission images show good qualitative agreement with experimental results, albeit with an apparent timing discrepancy for the two-sided vacuum target. More generally, we demonstrate the potential of hybrid-PIC simulations for full-system modeling and experimental design, including collisional absorption of laser energy, plasma evolution, mix, and fusion burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model fusion with physics-guided machine learning: Projection-based reduced-order modeling

The unprecedented amount of data generated from experiments, field observations, and large-scale numerical simulations at a wide range of spatiotemporal scales has enabled the rapid advancement of data-driven and especially deep learning models in the field of fluid mechanics. Although these methods are proven successful for many applications, there is a grand challenge of improving their generalizability. This is particularly essential when data-driven models are employed within outer-loop applications like optimization. In this work, we put forth a physics-guided machine learning (PGML) framework that leverages the interpretable physics-based model with a deep learning model. Leveraging a concatenated neural network design from multi-modal data sources, the PGML framework is capable of enhancing the generalizability of data-driven models and effectively protects against or inform about the inaccurate predictions resulting from extrapolation. We apply the PGML framework as a novel model fusion approach combining the physics-based Galerkin projection model and long- to short-term memory (LSTM) network for parametric model order reduction of fluid flows. We demonstrate the improved generalizability of the PGML framework against a purely data-driven approach through the injection of physics features into intermediate LSTM layers. Our quantitative analysis shows that the overall model uncertainty can be reduced through the PGML approach, especially for test data coming from a distribution different than the training data. Moreover, we demonstrate that our approach can be used as an inverse diagnostic tool providing a confidence score associated with models and observations. The proposed framework also allows for multi-fidelity computing by making use of low-fidelity models in the online deployment of quantified data-driven models.

42 ENGINEERING↗

Tokamak divertor plasma emulation with machine learning

Abstract Future tokamak devices that aim to create conditions relevant to power plant operations must consider strategies for mitigating damage to plasma facing components in the divertor. One of the goals of MAST-U tokamak operations is to inform these considerations by researching advanced divertor configurations that aid stable plasma detachment. Machine design, scenario planning and detachment control would all greatly benefit from tools that enable rapid calculation of scenario-relevant quantities given some input parameters. This paper presents a method for generating large, simulated scrape-off layer data sets, which was applied to generate a data set of steady-state Hermes-3 simulations of the MAST-U tokamak. A machine learning model was constructed using a Bayesian approach to hyperparameter optimisation to predict diagnosable output quantities given control-relevant input features. The resulting best-performing model, which is based on a feedforward neural network, achieves high accuracy when predicting electron temperature at the divertor target and carbon impurity radiation front position and runs in around 1 ms in inference mode. Techniques for interpreting the predictions made by the model were applied, and a high-resolution parameter scan of upstream conditions was performed to demonstrate the utility of rapidly generating accurate predictions using the emulator. This work represents a step forward in the design of machine learning-driven emulators of tokamak exhaust simulation codes in operational modes relevant to divertor detachment control and plasma scenario design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced System Thermal Fluids Solver Development for SAM

This work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key is the implementation of a high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. Leveraging the existing capabilities of the SAM code, significant code coverages were established in the finite volume method code. This in turn allows for a suite of test problems with different problem sizes and levels of complexity to be used to quantify the performance improvement of the finite volume method code. As evidently shown in this study, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE for the wide range of selected problems. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. In addition, for a complex reactor model, transient simulation was performed using the finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development. In this work, short-term priority development and testing items were identified, and long-term code adoption and integration plans were made for the eventual deployment of the finite volume method in the SAM code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Algorithms for the automatic generation of 2-D structured multi-block grids

Two different approaches to the fully automatic generation of structured multi-block grids in two dimensions are presented. The work aims to simplify the user interactivity necessary for the definition of a multiple block grid topology. The first approach is based on an advancing front method commonly used for the generation of unstructured grids. The original algorithm has been modified toward the generation of large quadrilateral elements. The second method is based on the divide-and-conquer paradigm with the global domain recursively partitioned into sub-domains. For either method each of the resulting blocks is then meshed using transfinite interpolation and elliptic smoothing. The applicability of these methods to practical problems is demonstrated for typical geometries of fluid dynamics.

Schoenfeld, Thilo↗

Calcium Isotope Constraints on Recycled Carbonates in Subduction-related Magmas

Calcium isotope ratios are readily mass-fractionated by low-temperature depositional processes that can be recorded in sedimentary rocks and, therefore, have the potential to track distinct geochemical signatures of recycled materials in mantle-derived igneous rocks. In this study, I report calcium isotopic compositions for well-characterized lavas from the Central American volcanic arc that exhibit a range of trace element and radiogenic isotope signatures that have been used to indicate variable amounts of subducted marine carbonate and hemipelagic sedimentary rocks, fresh carbonatite lavas from Oldoinyo Lengai Volcano, Tanzania, and an intrusive carbonatite clast erupted 12.8 ka from the Laacher See Volcano, Germany. I observed no calcium isotope evidence for recycled sedimentary rocks in the Central American arc magmas or in the Oldoinyo Lengai carbonatite volcanic rocks. They all exhibit isotopic compositions similar to rocks dominated by a primitive, mantle-like, bulk silicate Earth (BSE) composition. The exception in this work is the calcium isotope composition measured in the intrusive Laacher See carbonatite that is resolvable from BSE (δ^44 Ca = 0.0). Although this calcium isotopic signature might be related to ancient carbonate recycling, the magnitude and relatively light calcium isotopic composition (δ^44 Ca = –0.4) of this carbonatite could also reflect an origin that involved partial melting of altered lithospheric mantle, from which isotopically heavy Ca-bearing fluids have been lost. The decoupled signatures between trace elements and their radiogenic isotopes and those recorded by calcium isotope data are important because they likely reflect different sources and processes, and demonstrate that sediment subduction is not a bulk mixing process.

Ca isotopes↗

Semi-Annual Report for Horizontal Compact High Temperature Gas Reactor (HC-HTGR) Development during Performance Period October 2022-March 2023

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. Additionally, Argonne is providing analysis of the primary heat transport system to ensure temperatures in the reactor systems, structures and components with significant safety margins during normal operation and design basis accident scenarios. This third semi-annual report summarized the progress made at Argonne on the two tasks during the first half of FY23. As a part of the RCCS design task, a design process for the water panel and the system configuration was performed to improve the thermal performance of the RCCS. The updated water panel design had achieved enhanced thermal performance with 0.96 MWth capability with major design updates made in structural interfaces with the RPV and initial configuration of the water tanks. The loop configuration of the RCCS has been proposed to have two independently working loops for system redundancy. A water panel material study was performed focusing on the use of carbon steel in water systems. Additional modeling strategies of primary system thermal fluids analyses were developed to meet modeling needs that are not well suited for the 1D-3D assembly level approach. The first of these is a reduced order assembly model, in which fuel centered unit cells are used to represent a fuel assembly. The 2D approach used in this model is much more computationally efficient, making this model useful for full core transient scenarios where fuel to coolant heat transfer is still the dominant flow path. The next model is the 3D core conduction model to be used to analyze decay heat removal in loss of primary system flow scenarios. Because these scenarios require a large domain to be modeled, a homogenized core model is being pursued to reduce the required computational costs. To accurately model decay heat scenarios it is necessary to couple a RCCS model to the 3D core conduction model. A simplified case is presented to demonstrate how the coupling methodology will be applied to the full core 3D conduction model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High-Performance Java Codes for Computational Fluid Dynamics

The computational science community is reluctant to write large-scale computationally -intensive applications in Java due to concerns over Java's poor performance, despite the claimed software engineering advantages of its object-oriented features. Naive Java implementations of numerical algorithms can perform poorly compared to corresponding Fortran or C implementations. To achieve high performance, Java applications must be designed with good performance as a primary goal. This paper presents the object-oriented design and implementation of two real-world applications from the field of Computational Fluid Dynamics (CFD): a finite-volume fluid flow solver (LAURA, from NASA Langley Research Center), and an unstructured mesh adaptation algorithm (2D_TAG, from NASA Ames Research Center). This work builds on our previous experience with the design of high-performance numerical libraries in Java. We examine the performance of the applications using the currently available Java infrastructure and show that the Java version of the flow solver LAURA performs almost within a factor of 2 of the original procedural version. Our Java version of the mesh adaptation algorithm 2D_TAG performs within a factor of 1.5 of its original procedural version on certain platforms. Our results demonstrate that object-oriented software design principles are not necessarily inimical to high performance.

Riley, Christopher↗

Online Alpha Monitoring of High Cs-137 Hanford and SRS High Level Waste with Tensioned Metastable Fluid Detectors

The Department of Energy’s Hanford and Savannah River Sites maintain millions of gallons of caustic supernate and salt high activity waste in their high-level waste (HLW) tank farm inventories. The Savannah River Site is currently treating this waste with a calixarene-based solvent extraction of Cs-137 to reduce these inventories. Hanford is employing an at-tank crystalline silico titanate (CST) based solid phase extraction methodology to reduce their liquid HLW inventories. Due to the high solubility of Cs-137 and the relative insolubility of the actinides in these caustic waste forms, the beta to alpha radioactivity ratio can often exceed six orders of magnitude in the feed solutions to these treatment processes. This unique characteristic leads to significant technical challenges in making rapid gross alpha measurements in the presence of the overwhelming beta, gamma, as well as dissolved sodium salt in these HLW matrices. Conventional radioanalytical techniques, such as liquid scintillation analysis or gas flow proportional counting require significant radiochemistry preparation prior to the radiometric measurements for gross alpha activity. These required pretreatments render these technologies untenable for rapid quantification of gross alpha activity that could be required to support on or at-line measurements ensuring a waste stream will meet regulatory requirements. The radiation measurement properties of Tensioned Metastable Fluid Detectors (TMFDs) have been studied by Purdue University’s Taleyarkhan research group for well over a decade. Fluids tensioned to the appropriate degree will rupture when struck by radiation, resulting in a measurable cavitation event. The negative pressure generating this tension can be adjusted by centrifugal rotation or by acoustic means in such a way that these cavitation events can be generated from alpha radiation but will not be generated by beta or gamma radiation. Purdue University and the Savannah River National Laboratory are currently collaborating to develop a gamma/beta blind, spectroscopic alpha measurement system based on the Tensioned Metastable Fluid Detector technology to provide a potential solution for performing rapid gross alpha measurements on these high gamma/beta sample matrices. Measurements using the Indirect Drive Acoustically Tensioned Metastable Fluid Detectors developed as part of this collaboration were performed with an alpha emitting radionuclide. Successful determination of gross alpha activity was observed, indicating a potential pathway for rapid gross alpha measurements in remote-handled shielded cells or in process situations requiring online alpha monitoring. Measurements using this system have been conducted on high beta activity solutions, demonstrating the beta blind capability of this system. Measurements are currently underway to test the system’s capability to measure gross alpha activity on Savannah River Site high level waste high Cs-137 samples that have been measured by the SRNL radiochemistry team. This work was supported by the DOE EM Technology Development program.

DiPrete, David [Savannah River National Laboratory↗

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Aeroelastic Stability Assessment Methodology and Application to Slat Noise Treatments on the High-Lift Common Research Model

This work supports a larger effort at NASA to reduce airframe noise and, thus, environmental noise around airports. The leading-edge slat of conventional high-lift systems is a prominent source of airframe noise, and two technologies were previously identified as promising for noise reduction without significant aerodynamic or weight penalty: the slat cove filler (SCF) and the slat gap filler (SGF). NASA fabricated a 10%-scale semispan aircraft model for wind tunnel testing based upon the High-Lift Common Research Model (CRM-HL). NASA embarked on a test campaign with the 10% CRM-HL to determine the noise reduction effectiveness and other performance metrics of 3D SCF and SGF treatments in flow conditions representative of flight. The highly flexible SCF and SGF treatments posed a risk to the model and aeroelastic instability would compromise research integrity, so aeroelastic stability assessment was required. This work briefly describes the 10% CRM-HL, design of the SCF and SGF treatments to the model, and supporting fluid structure interaction (FSI) research. Relations for static and dynamic aeroelastic similitude are then developed for convenience of reference, and simplifications are introduced of particular utility to the cases considered in this study. The similitude relations are employed to assess the stability of the treatments under wind tunnel flow conditions based upon known behaviors from the supporting research. It is shown that the treatments are safe for test under the anticipated wind tunnel flow conditions with significant margins/factors of safety. It is furthermore demonstrated that the SCF is more susceptible to aeroelastic instability than the SGF, which has implications for flight feasibility.

Similitude↗

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗

Field Results from New Tensor Borehole Optical Fiber Strainmeter Installations in Oklahoma and Utah

The time evolving strain field contains a wealth of information that can be used to interpret subsurface behavior. For example, injecting or removing fluids from reservoirs or aquifers causes deformation that can be used as a diagnostic signal in some cases, while it can interfere with geodetic interpretations in other cases. We've previously demonstrated the feasibility of measuring the strain tensor at a depth of 30m caused by injection into a reservoir at 530m. The observed strain signals were interpreted using four independent analytic and numerical methods that resulted in estimates of the poroelastic properties and geometry of the reservoir that was consistent with data from well logs. However, studies like these are only possible if these deformations can be reliably measured. Years of lab and field work has culminated in the development of a novel borehole strainmeter capable of resolving multiple components of strain using embedded optical fibers configured as Michelson interferometers. It features four horizontal gauges separated by 45° to resolve the horizontal strain tensor as well as a vertical strain gauge and a sixth null component for state-of-health monitoring. The downhole sensing package also includes an open pipe through its center for grout circulation during single-trip deployments and a fully welded stainless steel exterior for robustness and longevity. These instruments have a resolution of 2x10-13 strain that can easily measure the solid earth tides. Preliminary data are available from four strainmeters in shale at our Oklahoma site and four in compacted sand and gravel in Utah. These are deployed from 40-60m, except one of the strainmeters in Oklahoma is deployed at 500m. The data include strains from the initial grout curing, comparisons to predicted earth tide models and in-situ calibration results, barometric pressure admittances and spectral analyses as well as signals from underground injections and surface waves from teleseismic events. Preliminary analyses indicate behavior consistent with other strainmeter deployments, and comparison to data from a Gladwin strainmeter at the Oklahoma site validate the performance of the new design. Analyses from a suite of six well tests at the Oklahoma site show for the first time how the strain tensor field varies with location during well testing.

DeWolf, Scott↗