Results for “Local sensitivity analysis” — Page 5

DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Temporal model of an optically pumped co-doped solid state laser

Currently, research is being conducted on the optical properties of materials associated with the development of solid state lasers in the two micron region. In support of this effort, a mathematical model describing the energy transfer in a holmium laser sensitized with thulium is developed. In this paper, we establish some qualitative properties of the solution of the model, such as non-negativity, boundedness, and integrability. A local stability analysis is then performed from which conditions for asymptotic stability are attained. Finally, we report on our numerical analysis of the system and how it compares with experimental results.

Wangler, T. G.↗

NASA NTRS1993

On the global dynamics of adaptive systems - A study of an elementary example

The inherent nonlinear character of adaptive systems poses serious theoretical problems for the analysis of their dynamics. On the other hand, the importance of their dynamic behavior is directly related to the practical interest in predicting such undesirable phenomena as nonlinear oscillations, abrupt transients, intermittence or a high sensitivity with respect to initial conditions. A geometrical/qualitative description of the phase portrait of a discrete-time adaptive system with unmodeled disturbances is given. For this, the motions in the phase space are referred to normally hyperbolic (structurally stable) locally invariant sets. The study is complemented with a local stability analysis of the equilibrium point and periodic solutions. The critical character of adaptive systems under rather usual working conditions is discussed. Special emphasis is put on the causes leading to intermittence. A geometric interpretation of the effects of some commonly used palliatives to this problem is given. The 'dead-zone' approach is studied in more detail. The predicted dynamics are compared with simulation results.

Espana, Martin D.↗

NASA NTRS2013

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

NASA NTRS2014

Characterization of MgB2 Superconducting Hot Electron Bolometers

Hot-Electron Bolometer (HEB) mixers have proven to be the best tool for high-resolution spectroscopy at the Terahertz frequencies. However, the current state of the art NbN mixers suffer from a small intermediate frequency (IF) bandwidth as well as a low operating temperature. MgB2 is a promising material for HEB mixer technology in view of its high critical temperature and fast thermal relaxation allowing for a large IF bandwidth. In this work, we have fabricated and characterized thin-film (approximately 15 nanometers) MgB2-based spiral antenna-coupled HEB mixers on SiC substrate. We achieved the IF bandwidth greater than 8 gigahertz at 25 degrees Kelvin and the device noise temperature less than 4000 degrees Kelvin at 9 degrees Kelvin using a 600 gigahertz source. Using temperature dependencies of the radiation power dissipated in the device we have identified the optical loss in the integrated microantenna responsible as a cause of the limited sensitivity of the current mixer devices. From the analysis of the current-voltage (IV) characteristics, we have derived the effective thermal conductance of the mixer device and estimated the required local oscillator power in an optimized device to be approximately 1 microwatts.

Superconducting Devices↗

NASA NTRS2019

Radio Follow-up of a Candidate γ-Ray Transient in the Sky Localization Area of GW170608

After the identification of a candidate γ-ray transient in the error region of the binary black hole (BBH) merger GW150914 by the Fermi satellite, the question of whether BBH mergers can be associated with electromagnetic counterparts remains highly debated. Here, we present radio follow-up observations of GW170608, a BBH merger that occurred during the second observing run (O2) of the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). Our radio follow-up focused on a specific field contained in the GW170608 sky localization area, where a candidate high-energy transient was detected by the Fermi Large Area Telescope (LAT). We make use of data collected at 1.4 GHz with the Karl G. Jansky Very Large Array (VLA), as well as with the VLA Lowband Ionosphere and Transient Experiment (VLITE). Our analysis is sensitive to potential radio afterglows with luminosity densities L(1.4GHz) ≳ 6 x 10^(28) erg/sHz. In the most optimistic theoretical models, ≈20% of BBH events occurring in massive hosts could be associated with outflows as radio luminous as this. Although we find no evidence for the presence of a radio counterpart associated with the LAT candidate in the GW170608 error region, our analysis demonstrates the feasibility of future radio follow-up observations of well-localized BBHs. Comparing our radio upper limits with theoretical expectations for the radio afterglows potentially associated with jets launched in BBH mergers, we find that for jets of energy ≈10^(49) erg seen on-axis, only jet angles θ(jet)≳40° are compatible with the observations.

Kyle Artkop↗

NASA NTRSDate not supplied

CHP-PRA Proof-of-Concept [Simulation] Sensitivity Assessment

An effort is underway to establish a Crew Health and Performance system (CHP) tradespace tool using a Probabilistic Risk Assessment (PRA) modeling and simulation system. The goal of the CHP-PRA effort is to provide a means of quantifying the integrated influence of CHP functions and capabilities on risk outcome metrics associated with health, performance, and long-term health. These metrics can then be used to establish potential risk-based trades on CHP system designed functionality and capabilities. Previously, our team demonstrated a proof-of-concept PRA approach that estimated the integrated influence of exercise countermeasures on 8 human system risks (Figure 1a) with outcomes associated with health and medical risk metrics. We reported that the change in the integrated relative health risk was small (Figure 1b) and that the small change in overall risk resulted from compounding and competing contribution levels of the individual risks. This interesting observation illustrates the emergent complexity of even straightforward representations of the human health and performance risk space and the ability of PRA models to capture this balance of global risk concerns. A key question that is not addressed in the initial analysis is “even though the global risk is relatively nominal, do any of the local risks become unacceptable?” In essence, we seek to determine what relative change in the human system risks are contributing to the relatively muted sensitivity of the proof-of-concept model combined risk assessments. Evaluations at the component risk level should elucidate if any individual risk reaches a high level that is subsequentially balanced by reductions in other areas. To further understand the relative changes in the component risks in the proof-of-concept model, and to elucidate how future refinements can be targeted, a means of establishing the contributions of the robustness of the proof-of-concept approach will be demonstrated.

astronaut health↗

NASA NTRS1995

The Thermal Ion Dynamics Experiment and Plasma Source Instrument

The Thermal Ion Dynamics Experiment (TIDE) and the Plasma Source Instrument (PSI) have been developed in response to the requirements of the ISTP Program for three-dimensional (3D) plasma composition measurements capable of tracking the circulation of low-energy (0-500 eV) plasma through the polar magnetosphere. This plasma is composed of penetrating magnetosheath and escaping ionospheric components. It is in part lost to the downstream solar wind and in part recirculated within the magnetosphere, participating in the formation of the diamagnetic hot plasma sheet and ring current plasma populations. Significant obstacles which have previously made this task impossible include the low density and energy of the outflowing ionospheric plasma plume and the positive spacecraft floating potentials which exclude the lowest-energy plasma from detection on ordinary spacecraft. Based on a unique combination of focusing electrostatic ion optics and time of flight detection and mass analysis, TIDE provides the sensitivity (seven apertures of about 1 cm squared effective area each) and angular resolution (6 x 18 degrees) required for this purpose. PSI produces a low energy plasma locally at the POLAR spacecraft that provides the ion current required to balance the photoelectron current, along with a low temperature electron population, regulating the spacecraft potential slightly positive relative to the space plasma. TIDE/PSI will: (a) measure the density and flow fields of the solar and terrestrial plasmas within the high polar cap and magnetospheric lobes; (b) quantify the extent to which ionospheric and solar ions are recirculated within the distant magnetotail neutral sheet or lost to the distant tail and solar wind; (c) investigate the mass-dependent degree energization of these plasmas by measuring their thermodynamic properties; (d) investigate the relative roles of ionosphere and solar wind as sources of plasma to the plasma sheet and ring current.

Moore, T. E.↗

NASA NTRS2018

Test and Analysis of Full-Scale 27.5-Foot-Diameter Stiffened Metallic Launch Vehicle Cylinders

The Shell Buckling Knockdown Factor Project (SBKF) was established with the goal of developing improved (i.e., less-conservative, more robust) shell buckling knockdown factors (KDFs) for modern launch-vehicle structures. To this end, SBKF has engaged in several activities to support the development, validation, and implementation of the new design factors, including subscale and full-scale structural testing. Tests on eight different subscale, 8-foot-diameter, integrally stiffened aluminum-lithium 2195 (Al-Li 2195) cylinders were conducted in order to obtain the majority of the required validation data. In addition, two full-scale, 27.5-foot-diameter, Al-Li 2195, integrally stiffened cylinders were tested to provide additional validation data and to determine structural scaling trends. Presented herein are the details of a recent analysis model development and test and analysis correlation effort on the full-scale test articles. The effects of selected modeling assumptions and approaches are discussed, and results from a modeling sensitivity study are presented. It was found that simplified finite element models, that assume nominal test article geometry and material properties, can predict the overall response characteristics well. However, several discrepancies in the test and analysis results were observed. A sensitivity study was performed to determine the effects of several modeling assumptions and address the observed discrepancies. The results from the study indicated that the evolution of local skin pocket buckling and the presence of residual stresses due to the manufacturing process can have a significant influence on the predicted buckling response of the cylinders considered.

Lovejoy, Andrew E.↗

DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗

NASA NTRS1994

Shortwave feedbacks and El Nino-Southern Oscillation: Forced ocean coupled ocean-atmosphere experiments

Changes in tropical sea surface temperature (SST) can produce changes in cloudiness that modify incoming solar shortwave (SW) radiation, which in turn affects SST. The effects of this negative feedback on Pacific interannual variability are examined in forced ocean model and hybrid coupled ocean-atmosphere model simulations. Two empirical schemes are used to model the large-scale, low-frequency response of surface SW to SST anomalies. The first scheme attempts to account for the nonlocal nature of the atmospheric response to SST based patterns of covariability analyzed through singular value decomposition. In the observations the primary coupled mode of variability is composed of a SW anomaly in the central Pacific that covaries with anomalous SST in the eastern Pacific. This is applied in the model as a nonlocal feedback. The second scheme examines the effects of a purely local feedback with a spatially varying coefficient of magnitude chosen similar to the first scheme. In almost all cases the second scheme behaved similarly to the first, presumably because the correlation scale of SST is large enough for El Nino-Southern Oscillation (ENSO) dynamics that there is little sensitivity to the local approximation in the SW feedback. In simulations forced by time series of observed wind stress the SW feedback induced very minor SST damping. Results for a simpified heat budget analysis showed that while the SW feedback increased the local heat flux damping on SST, it also induced a mean shallowing of the mixed layer. The resulting changes in both the local mean vertical temperature gradient and the zonal velocity response to the wind stress acted to oppose the local heat flux damping effects. When the observed SW anomalies were applied to forced simulations, the simulated SST anomalies were modified as expected, and agreement with observed SST improved. In coupled simulations the SW feedbacks had greater impact than in the case of specified stress. The main effects were to decrease the magnitude of the warm and cold SST anomalies in the central Pacific, while leaving the pattern and evolution of ENSO anomalies essentially unchanged elsewhere. The SW feedbacks thus produce a modest improvement of the model ENSO SST pattern compared with observations, although they tended to shorten the period of the model ENSO cycle. Overall the results suggest that large-scale SW feedbacks are of quantitative importance to simulating some aspects of the ENSO cycle but are not critical to the overall occurrence of the phenomenon.

Waliser, Duane E.↗

NASA NTRS2001

Palmitoylation of caveolin-1 in endothelial cells is post-translational but irreversible

Caveolin-1 is a palmitoylated protein involved in assembly of signaling molecules in plasma membrane subdomains termed caveolae and in intracellular cholesterol transport. Three cysteine residues in the C terminus of caveolin-1 are subject to palmitoylation, which is not necessary for caveolar targeting of caveolin-1. Protein palmitoylation is a post-translational and reversible modification that may be regulated and that in turn may regulate conformation, membrane association, protein-protein interactions, and intracellular localization of the target protein. We have undertaken a detailed analysis of [(3)H]palmitate incorporation into caveolin-1 in aortic endothelial cells. The linkage of palmitate to caveolin-1 was hydroxylamine-sensitive and thus presumably a thioester bond. However, contrary to expectations, palmitate incorporation was blocked completely by the protein synthesis inhibitors cycloheximide and puromycin. In parallel experiments to show specificity, palmitoylation of aortic endothelial cell-specific nitric-oxide synthase was unaffected by these reagents. Inhibitors of protein trafficking, brefeldin A and monensin, blocked caveolin-1 palmitoylation, indicating that the modification was not cotranslational but rather required caveolin-1 transport from the endoplasmic reticulum and Golgi to the plasma membrane. In addition, immunophilin chaperones that form complexes with caveolin-1, i.e. FK506-binding protein 52, cyclophilin A, and cyclophilin 40, were not necessary for caveolin-1 palmitoylation because agents that bind immunophilins did not inhibit palmitoylation. Pulse-chase experiments showed that caveolin-1 palmitoylation is essentially irreversible because the release of [(3)H]palmitate was not significant even after 24 h. These results show that [(3)H]palmitate incorporation is limited to newly synthesized caveolin-1, not because incorporation only occurs during synthesis but because the continuous presence of palmitate on caveolin-1 prevents subsequent repalmitoylation.

Non-NASA Center↗

DOE OSTI2025

A Type II Hamiltonian Variational Principle and Adjoint Systems for Lie Groups

We present a novel Type II variational principle on the cotangent bundle of a Lie group which enforces Type II boundary conditions, i.e., fixed initial position and final momentum. In general, such Type II variational principles are only globally defined on vector spaces or locally defined on general manifolds; however, by left translation, we are able to define this variational principle globally on cotangent bundles of Lie groups. Type II boundary conditions are particularly important for adjoint sensitivity analysis, which is our motivating application. As such, we additionally discuss adjoint systems on Lie groups, their properties, and how they can be used to solve optimization problems subject to dynamics on Lie groups.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

NASA NTRS2018

Analysis of Observed Contamination Through SAGE III's First Year on Orbit

SAGE III is a payload on the International Space Station that conducts measurements of ozone and other atmospheric constituents through the use of a moderate resolution spectrometer with an operating wavelength range of 290 nm to 1550 nm. Because of the optically sensitive nature of the payload, a suite of eight Thermoelectric Quartz Crystal Microbalances (TQCMs) were included to monitor the operating environment. During the rst year of operation, the SAGE III TQCMs were instrumental in detecting several periods of higher contamination and localizing their sources. A clear window made from quartz crystal covers the instrument assembly's aperture. Under nominal operating conditions, this window is only open during science gathering activities. However, if the rates of contamination accumulation are detected to be above the background rate, the window will be kept closed during science gathering to protect the optically sensitive instrument mirror. An analysis of the signal transmissions through the window for the wavelengths of 290 nm to 1550 nm has been conducted to determine any possible degradation of the window and potential in uence on science data collected to date, and established a baseline for future analysis.

Thompson, Michael P.↗

NASA NTRS2021

Local Ultrasonic Resonance Spectroscopy of Lithium Metal Batteries for Aerospace Applications

As next-generation aircraft and vehicles continue to develop, so do their associated energy demands. Lithium metal batteries are a leading candidate to fulfill this energy requirement, but these batteries are prone to internal dendrite defects that can lead to catastrophic thermal runaway events. Current battery management systems are capable of mitigating such risks, but are unable to detect such defects until thermal runaway has already begun. Various nondestructive evaluation (NDE) techniques, particularly ultrasonic NDE, can directly monitor internal battery parameters giving them the potential to detect critical defects prior to catastrophic failure. However, most of the current battery NDE research has focused on improved battery state-of-charge (SOC) and stateof- health (SOH) monitoring with little emphasis on critical defect detection. Thus, a measurement technique sensitive to subtle battery defects is needed. In addition, the complex mechanics of ultrasound in porous, thin, multilayered batteries prompt the use of physics-based simulation to guide inspections. In this work, an ultrasonic NDE technique has been developed utilizing frequency domain analysis of local battery resonances to detect the presence of battery defects. This technique is a practical extension of local ultrasonic resonance spectroscopy (LURS) – which previously required non-contact laser ultrasonics – to measurements with piezoelectric contact and immersion scan transducers. To extend the technique to work with piezoelectric transducers, ultrasonic battery measurements were compared to a sans-battery calibration measurement. Then, a linear systems deconvolution was used to eliminate the transfer functions of extraneous factors such as the transducer and electronics, leaving only the frequency-dependent battery reflection coefficient. The LURS technique was first validated on stainless steel and aluminum plates, producing reflection coefficients in line with analytical and numerical finite element modeling (FEM) results. Functioning Li-metal pouch cells were then seeded with lithium chip defects prior to LURS measurements. The presence of these defects is shown to cause a measurable shift in the battery’s through-thickness local resonances. 2D, frequency-domain poroelastic models of ultrasonic propagation in a single-cell lithium metal pouch battery were created and corroborated these findings. Thus, this work has both extended and proven the feasibility of the LURS technique in the detection of local battery defects.

Ultrasound↗