SEARCH · Search NASA

Results for “Observation geometry”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Tradeoffs and Variability of Mechanical Properties in Cold Worked 21-9-6 Stainless Steel Bars

Sub-standard yield stress in a batch of 21-9-6 stainless steel (Nitronic 40) bars was observed while other batches passed requirements. This material requires enhanced property uniformity due to its application in critical components and therefore has tighter restrictions on alloying elements. Using tensile testing, this study investigates the distribution of yield and ductility across bar geometry for three lots of material and is corroborated alongside Vickers microhardness profiles.

Estrada, Micah Garnet Watson↗

DOE OSTI2025

Design Principles in Engineering of Multigrain Nanocatalysts via Multiscale Electronic Structure Characterization

Engineering grain boundary (GB) strain provides a promising pathway to tune the catalytic properties of nanocrystals. However, structural heterogeneity from random grain orientation and geometry has limited clear structure–property correlations. Here, we utilize a multigrain Co3O4/Mn3O4 core/shell nanocrystal platform as a model system to systematically investigate how geometric misfit strain at GBs serves as catalytically active sites for the oxygen reduction reaction. Through precise subnanometer-level control over grain morphology and by integrating multiscale electronic structure characterization, we identify the electronic structural signature of GB defects and establish a direct correlation between localized strain fields and modified electronic states. Strain modulation at GBs alters the eg orbital energy levels, with elongation along the z-axis combined with shear strain stabilizing the eg states, in contrast to the destabilization observed under pure shear strain. This stabilization mechanism enhances the electrocatalytic activity and selectivity of strained GBs compared with strain-relaxed grain surfaces. Furthermore, we reveal that GBs exhibit a radial strain gradient, producing a spatial energy shift that further modulates local electronic structures, as resolved through the classification of electron energy loss spectroscopy data. Together, these findings demonstrate that geometric misfit strain enables precise tuning of grain geometry and the resulting electronic structures, offering a robust strategy for engineering next-generation nanocatalysts.

Cho, Min Gee↗

DOE OSTI2026

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗

DOE OSTI2024

Hierarchical Bayesian modeling for Inverse Uncertainty Quantification of system thermal-hydraulics code using critical flow experimental data

The best estimate plus uncertainty methodology in nuclear system thermal-hydraulic studies necessitates a comprehensive understanding of uncertainties in system code predictions. The forward uncertainty quantification (UQ) process involves the propagation of input uncertainties through the computational models to obtain uncertainties in the outputs. To this end, achieving an accurate estimation of input uncertainties is important, which is the focus of inverse UQ (IUQ). Traditionally, research in Bayesian IUQ within the nuclear engineering domain has largely relied on single-level Bayesian inference. While being effective for relatively small datasets, this approach encounters limitations for cases with large datasets. The use of a single-level model may prove inefficient, as the resultant posterior distributions can significantly differ when distinct subsets of data are employed. To address this issue, we employ an hierarchical Bayesian model for IUQ. Furthermore, this approach involves organizing observations into different groups based on the test conditions, thereby accommodating varying calibration parameters across these distinct groups. In this study, we developed and implemented a hierarchical Bayesian IUQ method to consider the grouping effect of critical flow measurement data from various geometries. Comparing the outcomes of IUQ under different selections of test data using hierarchical Bayesian IUQ against those obtained from single-level Bayesian IUQ, the forward propagation of hierarchical Bayesian IUQ results demonstrates a notably improved agreement with the experimental data.

42 ENGINEERING↗

DOE OSTI2025

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗

DOE OSTI2025

Macro-level mechanical interlocking: A rapid joining approach for additively manufactured compression molded composite panels

Composite joining typically involves multiple steps, such as drilling and surface treatment, as part of the manufacturing process, which leads to low throughput and long cycle times. In the present study, we demonstrated a macro-level mechanical interlocking (MI) based, rapid joining technique to assemble additively manufactured compression molded (AMCM) panels, enabling the production of parts larger than the mold dimensions. Composite panels made of 20 wt% short carbon fiber reinforced acrylonitrile butadiene styrene (CF/ABS) were joined using MI features of various geometries, namely tree (TR), dovetail (Dov), rectangle 2 (Rect2), and rectangle 1 (Rect1), and their in-plane strength was evaluated. The resultant strength of the tested MI joints reached up to 74 % of the baseline tensile strength (i.e., the ‘no joint’ case). Observations from optical and scanning electron microscopy revealed inadequate polymer diffusion between the adherends, indicating that the joint strength was primarily derived from mechanical interlocking. Additionally, the fracture surfaces exhibited stress-whitening marks, which were characterized using differential scanning calorimetry (DSC). The increase in melting enthalpy suggested local stretching of polymer chains due to MI. Finite element analysis (FEA) indicated that the Rect1 MI feature, which generated the lowest stress concentration, outperformed the others in terms of joint strength, achieving 42 MPa. As a demonstration of the MI joining method, a battery box tray measuring 108 cm × 34 cm using a mold with an effective dimension of 36 cm × 34 cm successfully manufactured, resulting in a part with an area three times larger than the mold. In conclusion, this study presents a promising approach to improving composite joining techniques while minimizing production complexities.

In-plane joining↗

DOE OSTI2024

Probing degradation at solid-state battery interfaces using machine-learning interatomic potential

Solid-state batteries featuring fast ion-conducting solid electrolytes are promising next-generation energy storage technologies, yet challenges remain for practical deployment due to electro-chemo-mechanical instabilities at solid-solid interfaces. These interfaces, which include homogeneous/internal interfaces such as grain boundaries (GBs) and heterogeneous/external interfaces between solid-electrolyte and electrode materials, can impede Li-ion transport, deteriorate performance, and eventually lead to cell failure. Here, in this study, we leverage large-scale molecular simulations, enabled by validated machine-learning interatomic potentials, to directly probe the onset of interfacial degradation at the garnet Li 7 La 3 Zr 2 O 12 (LLZO) solid-electrolyte/LiCoO 2 (LCO) cathode interface. By surveying different interfacial geometries and compositions, it is found that Li-deficient interfaces can lead to severe interfacial disordering with cation mixing and Co interdiffusion from LCO into LLZO. By contrast, Li-sufficient interfaces are less disordered, although elemental segregation with local ordering is observed. As a consequence of Co interdiffusion, Co-rich regions are formed at the GBs of LLZO due to cation segregation and trapping effects. This behavior is independent of the GB tilting axis, degree of disorder at the GBs, and Co concentration, which implies Co clustering at GBs is a general phenomenon in polycrystalline LLZO and can dictate its overall transport and mechanical properties. Our findings elucidate the underlying fundamental mechanisms that give rise to experimentally observed physicochemical properties and provide guidelines for interface design that can mitigate interfacial degradation and improve cycling performance.

25 ENERGY STORAGE↗

DOE OSTI2026

Substrate birefringence as a source of artifacts in spintronic THz emission spectra

Terahertz (THz) emission from spintronic THz emitters (STEs) has been extensively studied, both for its potential application in THz technologies and as a tool for probing spin dynamics, conductivity, and crystal anisotropy via time-domain THz spectroscopy. However, substrate birefringence can modify the emitted THz polarization, distorting the measured spectra and leading to potential misinterpretation of the results. While some aspects of substrate birefringence on THz emission have been studied in the time domain, its effects on THz spectra have not been discussed in detail. In the THz emission spectrum from STEs, we observe distinct spectral signatures that depend on sample orientation. To understand these signatures, we use STEs grown on c-cut sapphire and systematically vary the orientation of a (100) rutile TiO 2 window in both pump-through and THz-through geometries to identify the birefringence-induced effects. These findings demonstrate that substrate birefringence plays a critical role in the emitted THz spectrum and must be carefully considered in the analysis of time-domain THz spectroscopy data. We also show that artifacts from smaller birefringence can have a more misleading effect on the THz spectra.

Shrestha, Shreya [Univ. of Delaware, Newark, DE (U↗

DOE OSTI2026

Observation of flow vector fluctuations in p–Pb collisions at $\sqrt{{s}_{NN}}=5.02$ TeV

Measurements of transverse momentum (pT) and pseudorapidity (η) dependent flow vector fluctuations in p–Pb collisions at $$\sqrt{{s}_{NN}}=5.02$$ TeV at the CERN Large Hadron Collider are presented. By studying long-range two-particle correlations with a template fit method, potential biases from non-flow effects such as jets and resonance decays are effectively suppressed. Significant pT- and η-dependent fluctuations of the second-harmonic flow vector are observed with more than 5σ confidence in p–Pb collisions, similar to the observations in Pb–Pb collisions. The influence of residual non-flow effects has been evaluated and cannot account for the observed fluctuations, thereby confirming the observation of flow vector fluctuations in small collision systems at the LHC. Comparisons to model calculations from 3DGlauber+MUSIC+UrQMD and the parton transport model from AMPT are also presented. The measurements provide constraints on the theoretical modelling of the three-dimensional initial geometry and its event-by-event fluctuations, offering critical insights into the origin of collective flow in small collision systems at the LHC.

Abdallah, D A H↗

DOE OSTI2024

Assay-based background projection for the Majorana Demonstrator using Monte Carlo uncertainty propagation

The background index (BI) is an important quantity to project and calculate the half-life sensitivity of neutrinoless double-𝛽 decay (0⁢𝜈⁢𝛽⁢𝛽) experiments. An analysis framework is presented to calculate the BI using the specific activities, masses, and simulated efficiencies of an experiments components as distributions. This Bayesian framework includes a unified approach to combine specific activities from assay. Monte Carlo uncertainty propagation is used to build a BI distribution from the specific activity, mass, and efficiency distributions. This method is applied to the M AJORANA D EMONSTRATOR , which deployed arrays of high-purity Ge detectors enriched in 76 Ge to search for 0⁢𝜈⁢𝛽⁢𝛽. The original assay-based projection is requantified in the new framework, using the as-built geometry of the Demonstrator and additional assay information. While 47% higher than the original projection, the resulting BI of [8.95±0.36]×10 −4 cts/(keVkgyr) from the 232 Th and 238 U decay chains does not account for the higher-than-expected BI observed by the D EMONSTRATOR . Finally, this method enables us to demonstrate the statistical incompatibility between the D EMONSTRATOR 's observed background and the assay results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Why the US Needs an Enduring Plutonium Critical Assembly

Jezebel (operated from 1954-1977) has been the primary experiment for fast 239 Pu nuclear data validation for the last 70 years. Validation has included not only k eff , but also spectral indicies, Rossi-α, reactivity coefficients, and neutron leakage spectra. While this has been incredibly valuable to the community, there are three major issues. The first issue is that documentation on many of these experiments were lacking, leading to large uncertainties or (even worse) incorrect assumptions. The second issue is that while there have been many advances in research, there is no way to test those new advances today. The last issue is that since the assembly only operated for 23 years (and at time when 30 other critical assemblies were operating at the same facility and nuclear weapons testing was occurring), there were limited opportunities to observe how any system parameters changed as a function of time. Note that this work is not suggesting that a "Jezebel re placement" used for a limited experiment campaign would have great value. A new enduring (100 year target) plutonium (Pu) assembly with simple geometry and low uncertainties, however, would be extremely valuable. Such a capability would have a transformative impact on many research areas including nuclear data validation, analytical methods validation, dosimetry, reactor kinetics, and materials. The need for a new plutonium assembly is not new: it has been in the DOE Nuclear Criticality Safety Program (NCSP) Mission and Vision for over 10 years and was also the top priority established at the 2022 National Criticality Experiments Research Center (NCERC) Futures meeting. This work will discuss how such a new capability would help ensure that the US retains international leadership in plutonium research. Last, a brief overview of Lilith, a project aimed to design a new enduring plutonium assembly for operation at the NCERC will be given.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Isotope effects and Alfvén eigenmode stability in JET H, D, T, DT, and He plasmas

While much about Alfvén eigenmode (AE) stability has been explored in previous and current tokamaks, open questions remain for future burning plasma experiments, especially regarding exact stability threshold conditions and related isotope effects; the latter, of course, requiring good knowledge of the plasma ion composition. In the JET tokamak, eight in-vessel antennas actively excite stable AEs, from which their frequencies, toroidal mode numbers, and net damping rates are assessed. The effective ion mass can also be inferred using measurements of the plasma density and magnetic geometry. Thousands of AE stability measurements have been collected by the Alfvén Eigenmode Active Diagnostic in hundreds of JET plasmas during the recent Hydrogen, Deuterium, Tritium, DT, and Helium-4 campaigns. In this novel AE stability database, spanning all four main ion species, damping is observed to decrease with increasing Hydrogenic mass, but increase for Helium, a trend consistent with radiative damping as the dominant damping mechanism. These data are important for confident predictions of AE stability in both non-nuclear (H/He) and nuclear (D/T) operations in future devices. In particular, if radiative damping plays a significant role in overall stability, some AEs could be more easily destabilized in D/T plasmas than their H/He reference pulses, even before considering fast ion and alpha particle drive. Active MHD spectroscopy is also employed on select HD, HT, and DT plasmas to infer the effective ion mass, thereby closing the loop on isotope analysis and demonstrating a complementary method to typical diagnosis of the isotope ratio.

Alfvén eigenmodes↗

DOE OSTI2024

Advanced Instrumentation for Metal Additive Manufacturing

Laser powder bed fusion (LPBF) is the most widely used process for metal additive manufacturing (AM), particularly where complex geometries provide performance advantages unattainable with traditional manufacturing techniques. However, LPBF is highly sensitive to innate variability in both the powder spreading and fusion steps, often leading to defects such as pores that are difficult to detect yet significantly impair component mechanical properties and fatigue life. This thesis presents a range of novel instruments enabling both precise assessment of powder layer characteristics and in-situ thermal metrology of metal AM to advance the quality control of LPBF. First, leveraging a custom X-ray microscope and a radiation-transport model developed through this work, transmission X-ray imaging is used to study spreading of thin metal powder layers. Effective layer depth is directly mapped at a process-relevant size scale, surpassing optical techniques that can only estimate local deposition from layer surface topography. Layer packing density and quality are shown to be influenced by powder flowability and particle size relative to nominal powder layer thickness. Layer quality is additionally connected to the geometry of the spreading implement and its velocity. This technique and its presented findings enable pairing feedstocks with spreading strategies that create layers with consistent packing density and uniformity. Second, a twofold approach is employed to optically interrogate the laser fusion step of LPBF for observing signatures of defect formation. Aperture division multiplexing is conceptualized, providing for simultaneous laser delivery and high-fidelity infrared (IR) process monitoring through a common optic. In-situ microscopy at 50 μm spatial resolution and at mid-wave IR wavelengths is proven readily achievable with the first purpose-built optic of this type. Next, a bespoke imaging spectrometer, along with a temperature-emissivity separation technique, is used to retrieve accurate process temperatures over a 1000 K range. Data from these instruments are correlated to porosity as fine as 4.3 μm in two LPBF test artifacts, as verified using computed tomography (CT), establishing the viability of robust optically-based component qualification.

Penny, Ryan↗

DOE OSTI2025

High-Pressure Effects on an Octa-Hydrated Curium Complex: An Experimental and Theoretical Investigation

An octa-hydrated curium compound [Cm(H 2 O) 8 ](Hdtp)(dtp)·H 2 O (Cm1,H 2 dtp = 2,3-di(tetrazol-5-yl)pyrazine) along with its lanthanide analogues [Ln(H 2 O) 8 ](Hdtp)-(dtp)·H 2 O (Ln1, Ln 3+ = La 3+ –Nd 3+ , Sm 3+ –Lu 3+ ) were synthesized and characterized using single crystal X-ray diffraction and spectroscopic methods. Bond length analysis of VIII Cm(III)–OH 2 (where VIII refers to the coordination number) was compared to VIII Ln(III)–OH 2 (Ln 3+ = Nd 3+ and Sm 3+ ), indicating similar VIII M(III)–OH 2 bond lengths owing to the similar eight-coordinate ionic radii of these VIII M(III) cations. Owing to the reduced coordination number, the VIII Cm(III)–OH 2 bond lengths were shorter than previously reported IX Cm(III)–OH 2 bonds in [Cm(H 2 O) 9 ](CF 3 SO 3 ) 3 . The octa-aquo complexes were also characterized by solid-state UV–vis–NIR spectroscopy in addition to variable-temperature and variable-pressure photoluminescence. Variable-pressure absorption spectra of Cm1 were compared with Ln1 and show that the Cm(III) f → f transitions have a stronger dependence on pressure than that observed in Ln1 (Ln 3+ = Nd 3+ and Sm 3+ ). The experimental and computational analyses reveal that the monotonic decrease in the computed energy difference between the ground state and the first excited state corresponds to the observed red shift of the photoluminescence peak. This is accompanied by a gradual reduction in the average Cm(III)–OH 2 bond length and a delocalization of spin densities, alongside an intensified interaction involving the 5f orbitals under increasing pressure. Furthermore, these changes accommodate the new geometry and collectively modify the energy landscape, resulting in peak broadening and quenching.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Simulations of nozzle gas flow and gas-puff Z-pinch implosions on the Weizmann Z-pinch

We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.

Physics↗

DOE OSTI2026

NASA’s Pandora SmallSat Mission: Simulating the Impact of Stellar Photospheric Heterogeneity and Its Correction

Stellar photospheric heterogeneity is a dominant astrophysical systematic impacting exoplanet transmission spectroscopy. NASA’s Pandora SmallSat Mission is designed to address this challenge through contemporaneous visible-band photometry and near-infrared spectroscopy of exoplanet host stars. Here, we present an end-to-end simulation study quantifying Pandora’s ability to infer stellar photospheric properties and correct stellar contamination using out-of-transit observations. We construct eight representative stellar activity scenarios and generate 160 simulated Pandora datasets, incorporating time-dependent stellar spectra, instrument response, and noise. Given accurate models, Bayesian retrievals of joint visible photometry (0.4–0.7 μm) and near-infrared spectroscopy (0.9–1.6 μm, R ≈ 120) recover photospheric temperatures with typical uncertainties of ≈30 K, with no significant bias. Models with two spectral components (i.e., a quiescent photosphere and spots) are strongly favored in 95% of cases; one-component models are preferred when true spot filling factors fall below a detection threshold of ≈0.3%. We propagate the true and inferred stellar parameters to compute true, inferred, and residual contamination signals under physically motivated spot geometries. For simple spot distributions, contamination signals of 10 2 –10 3 ppm are reduced to ≲10 ppm—well below Pandora’s expected transmission spectroscopy precision (30–100 ppm). For more complex spot distributions, geometric degeneracies limit deterministic corrections, leaving residual contamination at the 10 3 ppm level that must be mitigated using additional constraints, such as spot-crossing events and joint stellar–planetary retrievals of transmission spectra. These results define regimes in which stellar contamination can be corrected from stellar observations alone and show how Pandora stellar observations can identify cases where additional information is required.

Astronomy and AstroPhysics↗

DOE OSTI2026

From MAST-U to STEP: Power exhaust control challenges and opportunities

This work explores the challenges and opportunities for power exhaust control in STEP, informed by dedicated MAST-U experiments. The STEP system primarily relies on hydrogenic and argon gas injection into the divertors, compensating for transients originating in the plasma core. Fast transients (>10 Hz) that evolve too quickly for actuators to compensate must be absorbed passively; the enhanced transient buffering of the foreseen long-legged divertor provides a clear advantage over conventional geometries, as demonstrated in the MAST-U Super-X divertor. Experiments further indicate that STEP’s tight divertor baffling enables near-independent control of the upper and lower divertors. The most challenging transients are expected from power-sharing fluctuations, as MAST-U experiments observed extremely fast dynamics. Fluctuations from core pellet fuelling remain relatively benign due to the small pellet size according to simulations. Exhaust control in a reactor like STEP requires an integrated approach befitting its machine-critical nature. A predictive control element, integrated with core plasma control, allows pre-emptive preparation of the divertor for incoming transients. An observer-based diagnostic approach is foreseen to monitor the divertor in the challenging reactor environment, supported by an extended diagnostic set in the non-nuclear phase to validate the required dynamic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite↗