Search NASA⌕ Search

SEARCH · Search NASA

Results for “explainability”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analyzing the impact of design factors on solar module thermomechanical durability using interpretable machine learning techniques

Solar modules in utility-scale systems are expected to maintain decades of lifetime to rival conventional energy sources. However, cyclic thermomechanical loading often degrades their long-term performance, highlighting the importance of effective design to mitigate thermal expansion mismatches between module materials. Given the complex composition of solar modules, isolating the impact of individual components on overall durability remains a challenging task. In this work, we analyze a comprehensive data set that comprises bill-of-materials (BOM) and thermal cycling power loss from 251 distinct module designs to identify the predominant design factors and their impacts on the thermomechanical durability of modules. The methodology of our analysis combines machine learning modeling (random forest) and Shapley additive explanation (SHAP) to correlate design factors with power loss and interpret the model’s decision-making. The interpretation reveals that silicon type (monocrystalline or polycrystalline), encapsulant thickness, busbar numbers, and wafer thickness predominantly influence the degradation. With lower power loss of around 0.6% on average in the SHAP analysis, monocrystalline cells present better durability than polycrystalline cells. This finding is further substantiated by statistical testing on our raw data set. The SHAP analysis also demonstrates that while thicker encapsulants lead to reduced power loss, further increasing their thickness over around 0.6 to 0.7 mm does not yield additional benefits, particularly for the front side one. In addition, other important BOM features such as the number of busbars are analyzed. This study provides a blueprint for utilizing explainable machine learning techniques in a complex material system and can potentially guide future research on optimizing the design of solar modules.

14 SOLAR ENERGY↗
DOE OSTI2025

Caught in headlights: Captive white-tailed deer responses to variations in vehicle lighting during imminent collision scenarios

Vehicle collisions with deer (Odocoileus spp.) cause billions of dollars in damages and injure thousands of drivers every year in the United States, and few mitigation methods have proven effective. However, recent research suggests that vehicle lighting might influence white-tailed deer (Odocoileus virginianus; hereafter, deer) responses to oncoming vehicles. Most new vehicles are manufactured with light emitting diode (LED) headlights which differ in total radiance and wavelength of light emitted compared to the previous industry standard of tungstenhalogen (halogen) headlights. Also, frontal vehicle illumination through rear-facing lighting has shown promise in enhancing deer responses to vehicles, but its effectiveness has not been tested under various headlight conditions (headlight type or intensity). As such, it remains unclear how these aspects of vehicle lighting affect deer responses to an approaching vehicle. We used 23 captive, wild-type deer to investigate how variations in vehicle lighting affect deer responses to an approaching vehicle at night, when most collisions occur. We released deer into a 95 m long, 3 m wide chute and approached them from the opposite end with an electric golf cart outfitted with two versions of stock 2017–2020 Ford Fusion headlights (LED and halogen) and a 51 cm rear-facing lightbar to test how vehicle lighting affected deer avoidance behaviors in an imminent, head-on collision scenario. Each deer received eight lighting treatments consisting of unique combinations of headlight type (LED vs. halogen), light intensity (low vs. high beam), and rear-facing lighting (lightbar off vs. on). We measured deer alert and flight behavior using infrared videography. We found that the halogen, high beam, lightbar off treatment had the greatest probability of evoking an alert response. Furthermore, when the lightbar was off, high beams appeared to increase alert probability for halogen headlights. Also, we found evidence that high beam, halogen headlights tend to increase alert probability over high beam, LED headlights, when the lighbar was off. We found no effect of our lighting treatments on deer alert distance, flight probability, or flight initiation distance. Across all behavioral responses, the random effect deer ID explained 0.86–9.19 × more variation than our lighting treatments, reflecting large differences in responses among deer. Overall, we found that vehicle lighting can impact deer behavior during an imminent, head-on collision scenario, although lighting was ineffective at increasing favorable flight behaviors. Future research should investigate how vehicle lighting treatments affect free-ranging, wild deer in a variety of real-world scenarios and at longer approach distances.

White-tailed deer (Odocoileus virginianus) Deer be↗
DOE OSTI2025

Burn parameters affect PAH emissions at conditions relevant for prescribed fires

Wildfire smoke is a health hazard as it contains carcinogenic volatile compounds and fine particulate matter. In particular, exposure to polycyclic aromatic hydrocarbons (PAHs) is a major concern, since these compounds have been recognized as important contributors to the overall carcinogenic risk. In this work, gas and particle-phase PAH emissions from combustion of Eastern White Pine (Pinus strobus) were quantified using time-of-flight mass spectrometry over a range of burn conditions representative of wildfires and prescribed fires, including fuel moisture, heat flux, and oxygen concentration. We found that changing the burn environment lead to a variability of up to 77% in phenanthrene/anthracene emissions. This could explain a large part of the variability in PAH emission factors from biomass combustion reported in the literature. Here, we found that optimal conditions for fuel moisture content of 20–30%, sample heat load of 60 - 70 kW m -2 , and oxygen concentrations of 5–15% can significantly reduce the emissions of heavy molar weight PAHs. Our analysis showed that the relative carcinogenic risk from PAH exposure can be reduced by more than 50% under optimal conditions. In light of the increasing use of prescribed fire for forest management, the relationship between emissions and burn conditions that we have established provides a guidance for assessing the expected health impact from prescription burns, and can inform strategies to reduce PAH emissions from prescribed fire activities.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Manganese concentration influences nitrogen cycling in agricultural soil

Manganese (Mn) can modulate nitrogen (N) transformations in soil, yet its role in agroecosystems remains understudied. We conducted a 51-day microcosm incubation with agricultural soils differing in long-term N history (N 0 , no added N; N 1 , added 225 kg N ha −1 for a duration of 27 years) and amended with soluble Mn at 0 (M 0 ), 50 (M 1 ), or 250 (M 2 ) mg kg −1 and Glyceria striata (Lam.) residue. In N 1 soils, Mn additions (both M 1 and M 2 ) lowered total mineral N by 25 % relative to N 1 M 0 at day 51 and reduced 51-day cumulative N₂O by 32 % (N 1 M 1 ) and 42 % (N 1 M 2 ) vs. N 1 M 0 , whereas effects in N 0 were negligible. Mn also depressed ammonia-oxidizing bacterial amoA gene transcripts at day 15 in N 1 M 2 vs. N 1 M 0 (2.5 fold change). This reduction was likely due to increased N loss via complete denitrification to N₂ through microbial pathways such as nitrate/nitrite-dependent manganese oxidation (NDMO), where bacteria directly used the added Mn 2+ to reduce nitrate (NO 3 − ) and nitrite (NO 2 − ) to N 2 or as Mn dependent-ammonia oxidation (Mnammox) where bacteria oxidized ammonium (NH 4 + ) to N 2 , using Mn oxides as electron acceptors that formed from the oxidation of the added Mn. Other contributing mechanisms may include Mn-induced N immobilization, toxicity, and changes in the microbial community. These mechanistic results indicate that background Mn availability and redox dynamics can shape nitrification–denitrification pathways under N-rich conditions. In conclusion, we highlight how native Mn pools and redox state may help explain observed variability in N losses and greenhouse gas production across agricultural soils.

Decomposition↗
DOE OSTI2026

Response of soil nutrient pools and microbiomes to recurrent wildfire disturbance and varying burn severities in a mixed conifer forest

Wildfire is a pervasive disturbance in mixed-conifer forests, yet the relative influence of fire recurrence versus burn severity on soil biogeochemistry and microbial communities remains poorly quantified. We examined a natural gradient of fire history (0–3 prior fires) and burn severity (low–high) spanning 50 yr in a mixed-conifer ecosystem to assess how repeated fire shapes soil carbon (C) and nitrogen (N) pools, their isotopic signatures, mineral and particulate fractions, microbial community composition, carbon-use, CO₂ fluxes, and vegetation cover. Successive fires produced progressively higher bare-ground percentages and lower tree cover, which were tightly linked to declines in microbial diversity and reductions bulk %C, and %N. δ 13 C increased with fire frequency, indicating preferential loss of labile C through combustion or enhanced microbial oxidation, thereby explaining the observed net soil-C decline. Conversely, δ 15 N decreased and pH increased as tree density declined, reflecting altered N cycling and reduced acidification in post-fire soils. Fire recurrence, more than severity, corresponded with a marked shift in the bacterial community: for example, Xanthobacteraceae—key N-fixers and C-cyclers—diminished, while N-fixing Bacillaceae increased, underscoring the tightly coupled nature of soil nutrient dynamics and microbiome composition after repeated burns. Our results demonstrate that fire recurrence appears to be a stronger driver of post-fire soil ecosystem responses in this mixed-conifer forest, influencing both abiotic nutrient pools and the functional potential of the soil microbiome. These findings provide a more enhanced assessment and understanding to date of the biogeochemical consequences of repeated wildfire disturbance that can be used to inform management strategies aimed at preserving soil health in fire-prone landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Strain-modulated intercalated phases of Pb monolayer with dual periodicity in SiC(0001)-graphene interface

Intercalation of metal atoms at the SiC(0001)-graphene (Gr) interface can provide confined 2D metal layers with interesting electronic properties. The intercalated Pb monolayer (ML) has shown the coexistence of the Gr(10 x 10)-moiré and a stripe phase, which still lacks understanding. Using density functional theory calculation and thermal annealing with ab initio molecular dynamics as motivated by experiment, we have studied the formation energy of Gr/Pb/SiC(0001) for different Pb coverages. Near the coverage of a Pb(111)-like ML mimicking the (10 x 10)-moiré, we find a slightly more stable stripe structure, where one half of the structure has compressive strain with Pb occupying the Si-top sites and the other half has tensile strain with Pb off the Si-top sites. This stripe structure along the Gr zigzag direction has a periodicity of 2.3 nm across the [1$\overline{2}$10] direction agreeing with the previous observations using scanning tunneling microscopy. Analysis with electron density difference and density of states show the tensile region has a more metallic character than the compressive region, while both are dominated by charge transfer from Pb ML to SiC(0001). As a result, the small energy difference between the stripe and Pb(111)-like structures means the two phases are almost degenerate and can coexist, which explains the experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Source attribution of near-surface ozone pollution in Jiangsu Province of China over 2013–2019

Near-surface ozone (O3) is one of the most severe air pollutants in China, particularly over densely populated Jiangsu Province in the Yangtze River Delta. In this study, an O3 source tagging technique is utilized in a chemistry-climate model to quantify the source contributions of various emission sectors and regions for nitrogen oxides (NOx) and volatile organic compounds (VOCs) to O3 concentrations in Jiangsu Province during 2013–2019. The results show that the near-surface O3 in Jiangsu Province is mainly contributed by surrounding and remote anthropogenic NOx emissions through long-range transport. Local anthropogenic NOx emissions account for only 13 % and 18 % of the annual and summertime mean near-surface O3 in Jiangsu Province, respectively. Anthropogenic NOx emissions from the surface transportation, industry, and energy sectors account for 21 %, 22 % and 20 % of the annual mean near-surface O3 concentration in Jiangsu Province, respectively. Biogenic and anthropogenic VOCs emissions each explains one-third of the annual mean near-surface O3 concentration in Jiangsu, while methane and stratospheric chemical production contribute 21 % and 6 %, respectively. The sources from stratospheric production, aircraft, lightning, and foreign emissions are the primary contributors to O3 in the mid- and high troposphere. During high pollution days in Jiangsu Province, the near-surface O3 concentrations increase with the maximum exceeding 20 ppb, which is attributed to both the enhanced photochemical production and regional transport in favorable meteorological conditions.

ozone pollution↗
DOE OSTI2026

Ammonia in northeast Colorado is increasing, rising most quickly in regions close to confined animal feeding operations

The Colorado Front Range urban corridor and nearby agricultural operations are important source regions of atmospheric ammonia (NH 3 ). Upslope flows periodically transport these emissions into Rocky Mountain National Park (RMNP), located 50 km west of the urban corridor, where wet and dry deposition of excess reactive nitrogen (N) impacts ecosystems. Here, we use a combination of in situ passive NH 3 measurements and NH 3 vertical column density retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) to assess variability and changes in NH 3 across three land use categories in the northeast Colorado source region (agricultural, urban, and remote) during the period 2013-2023. A strong seasonal cycle is present across the region with increased NH 3 during summer months. Elevated NH 3 is spatially correlated with the number of permitted animal units in confined animal feeding operations (CAFOs) within 12 km. Ground-level NH 3 concentrations are strongly positively correlated with monthly gridded IASI satellite column densities. Satellite retrievals reveal an increasing trend in NH 3 column amounts of ∼3% per year in agricultural and ∼2% per year in urban sub-regions. The magnitude of the trend observed in NH 3 columns averaged over the agricultural sub-region is > 3 times larger than observed near and over Denver. The largest increases in NH 3 are closely aligned with the distribution of CAFOs. Reductions in particle sulfate associated with declining sulfur dioxide (SO 2 ) emissions could account for only ∼0.1% per year increase in gaseous NH 3 . Wildfire smoke across the region has increased but appears unlikely to explain the majority of the observed NH 3 increase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Can we mitigate tropical cyclone formation using aerosols? A review of cyclogenesis and aerosol effects as a theoretical basis

Tropical cyclogenesis is a complex phenomenon, involving processes across multiple scales, with various theories proposed to explain its mechanisms. This paper provides an up-to-date summary of the established physical mechanisms and controlling factors of cyclogenesis. Understanding the nature of tropical cyclone formation is necessary to improve tropical cyclone forecasts and inform tropical cyclone mitigation strategies targeting early-stage intervention. Despite the disparities among theories, the critical role of persistent, intensified convection in tropical cyclone initiation is universally acknowledged. Therefore, our hypothesis is that interrupting convective development using aerosol injection is a potential strategy for mitigating tropical cyclone formation. Targetted injection of coarse and fine aerosols at the formation stage are two proposed potential approaches for mitigating tropical cyclones. These methods have completely different microphysical effects but each may ultimately interfere with the formation of the nascent cyclone. In conclusion, these potential interventions open new areas of further research, which are important and necessary to establish a solid scientific foundation for cyclogenesis mitigation.

Tropical cyclone formation↗
DOE OSTI2025

N-terminal domain swapping: A new paradigm for spermidine/spermine N -acetyltransferase (SSAT) protein structures?

Enterococcus faecalis is a multi-drug-resistant human pathogen that is found in a variety of environments and is challenging to treat. Under stress conditions, some bacteria regulate intracellular polyamine concentrations via polyamine acetyltransferases to reduce their toxicity. The E. faecalis genome encodes two polyamine acetyltransferases: PmvE and BltD. Both of these proteins belong to the Gcn5-related N-acetyltransferase (GNAT) superfamily. It is unclear why there are two enzymes with similar substrate specificities in this organism. To better understand the structure/function relationship of the E. faecalis BltD enzyme, we determined its crystal structure and performed additional assays to explore its oligomeric state and enzymatic activity. The goal was to determine whether there were structural or catalytic differences between this enzyme and other polyamine acetyltransferases that could explain this redundancy and be exploited for future development of targeted inhibitors for this important human pathogen. We found the BltD enzyme was structurally unique due to its N-terminal domain swapped dimer. However, this enzyme adopts a catalytically active monomer rather than dimer in solution. This indicates the crystal structure we obtained may represent a state that forms at high protein and salt concentrations and at low pH used during crystallization. The BltD dimer found in the crystal may represent a unique view of how an inhibitory peptide or molecule could be designed to occupy its active site. Additionally, this structure shows the extensive flexibility of the N-terminal portion of the E. faecalis BltD enzyme.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Nanoparticle-mediated antagonism of sustained endosomal signaling of the calcitonin receptor-like receptor provides enhanced and persistent relief of oral cancer pain

By improving the delivery and tumor retention of chemotherapeutics, nanomedicines hold potential for cancer treatment. The usefulness of nanoparticle (NP)-encapsulated analgesics for the cancer pain treatment is comparatively unexplored. We investigated whether NPs encapsulating olcegepant (OCP), an antagonist of the calcitonin receptor-like receptor (CLR) for the calcitonin gene-related peptide (CGRP), effectively relieved oral cancer pain in mice. Because persistent endosomal CLR signaling in Schwann cells mediates craniofacial pain, we reasoned that the predisposition of NPs to accumulate in endosomes could be leveraged to effectively relieve oral cancer pain. By expressing biosensors for activated CLR, Gα proteins and β-arrestins in HEK293T and Schwann cells, we found that CGRP activates CLR signaling first at the plasma membrane and then in early, late and recycling endosomes and the cis- and trans-Golgi apparatus. We synthesized biocompatible NPs encapsulating OCP and fluorophores by integrating hydrophobic ion pairing nanoformulation with Flash NanoPrecipitation. NPs slowly released OCP and accumulated in early endosomes, leading to sustained inhibition of endosomal CLR signaling in HEK293T and Schwann cells. Oral cancers were established in mice, which led to heightened pain-like responses. After intra-tumoral injection, NPs were retained in tumors for at least one week. OCP-loaded NPs almost completely reversed allodynia and hyperalgesia for a prolonged period, whereas unencapsulated OCP had small and transient effects. The NP accumulation in endosomal sites of pain signaling, the sustained release of antagonist, and the retention of NPs in tumors explain their beneficial actions. Thus, NP-encapsulation holds promise for the relief of painful cancers that are inadequately treated by opioids.

Calcitonin gene-related peptide↗
DOE OSTI2024

Mechanistic study of direct coupling of CO 2 and C 2 H 4 over atomically dispersed metal at graphene edges

Direct coupling of CO 2 and ethylene (hereinafter DCCE) to acrylic acid is valuable for valorizing CO 2 to manufacture acrylate-derived products. However, previous studies in DCCE have been limited on molecular catalysts with challenges in improving catalytic performance. In this work, we employed density functional theory calculations and ab initio molecular dynamics simulations to investigate the heterogeneous catalysis of DCCE over atomically dispersed metal centers at nitrogen-doped zigzag edge of graphene. Based on competitive adsorption and structural stability, Mo, Cr, V, Ru, and Ni active sites are chosen to explore the reaction kinetics. Here, we find that the activation barriers are determined by the charge redistribution at transition states, which explains the trend of activity for the C-C coupling and the hydrogen transfer, two key steps in DCCE. Furthermore, we show that the intramolecular hydrogen transfer (rate-limiting step) is hindered due to the lack of local coordinate at the active sites. We thus propose to use co-adsorbed water as a “proton-exchanger” following a water-assisted route, and show that the activation barriers are reduced over all metal centers. Particularly, water promotes the hydrogen transfer over metals with strong CO 2 -ethylene co-activation and facile C-C coupling kinetics, which could be considered promising for DCCE. In both mechanisms, the stability of metallactone intermediate can be used to predict the catalytic activity. It is anticipated that the insights from this work can provide guidelines for mimicking well-defined multifunctional active sites in molecular catalysts to design heterogeneous catalysts for such C-C coupling, which advances catalytic utilization of CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Dual aging pathways of Cu-SSZ-13 SCR catalysts: Hydrothermal vs. sulfur-induced deactivation

Hydrothermal aging (HTA) and chemical poisoning are two primary factors contributing to the real-world degradation of Cu-SSZ-13 SCR catalysts. Investigating field-returned samples offers valuable insights into performance degradation caused by these mechanisms. However, the simultaneous presence of both deactivation pathways complicates the isolation of their individual effects in post-mortem analyses. In this study, we separately prepared model Cu-SSZ-13 SCR catalysts subjected to hydrothermal-aging and sulfur-induced chemical poisoning. Using various characterization techniques, we elucidated the specific role of each aging process in catalyst deactivation and compared the results to real-world field-aged catalysts. Our findings show that hydrothermal aging at 650 °C for 100 h caused dealumination of the zeolite framework but no significant CuO x cluster formation. In contrast, sulfur aging (via sulfur exposure, calcination at 550 °C, and desulfation up to 750 °C) led to CuO x formation without any observable dealumination. On model catalysts, sulfur poisoning was found to reduce Cu mobility and the amount of active Cu sites, thus degrading catalyst activity. Although some activity was recovered upon desulfation, a portion of the initial catalyst activity remained irreversibly lost due to CuO x formation. We demonstrate that this occurs because sulfated species impede the ability of multi-nuclear Cu species (e.g., Cu dimers) to split back into their isolated form, leading to CuSO 4 -clusters that oxidatively desulfate to CuO x species. This degradation pathway explains the significant reduction in activity of field-aged samples, where substantial CuSO 4 -cluster accumulation leads to reduced active Cu and subsequent conversion to CuO x . Furthermore, the conclusions from model catalysts were extended directly to field-aged commercial samples, elucidating the decline in activity and chemical properties during field deployment.

Catalyst Deactivation↗
DOE OSTI2025

Persistent fatty acid catabolism during plant oil synthesis

Plant lipids are an essential energy source for diets and are a sustainable alternative to petroleum-based fuels and feedstocks. Fatty acid breakdown during seed germination is crucial for seedling establishment but unexpected during seed filling. Here, we demonstrate that the simultaneous biosynthesis and degradation of fatty acids begins early and continues across all phases of oil filling and throughout the photoperiod. Tests in camelina, rapeseed, and an engineered high-oil tobacco line confirmed that concomitant synthesis and breakdown in oil-producing tissues over development is the rule rather than the exception. Furthermore, we show that transgenics, designed to elevate fatty acid biosynthesis, failed to achieve anticipated increases in storage lipid levels due to increased degradation, potentially explaining the underperformance of engineered lines compared to expectations more generally.

coenzyme A↗
DOE OSTI2024

Mechanical degradation of duplex SiC-fiber reinforced SiC matrix composite tubes under a controlled high-temperature steam environment

This paper describes results from the mechanical evaluation of unirradiated SiC fiber–reinforced SiC matrix composite tubes under a controlled high-temperature steam environment. The experiments were aimed at identifying key material degradation behavior under environments relevant to loss-of-coolant accidents of light water reactors. Mechanical tests of the SiC composite tubes at 1000 °C under steam and inert environments were conducted using a unique test capability. Further, the material tested was a duplex tube with a thick monolithic SiC layer on the outer surface. The tubes were subjected to preloading at ~100 MPa in tension before exposure to high-temperature steam with up to 75% of the preload at a constant displacement. In the presence of matrix cracks, the steam exposure caused embrittlement of the SiC composite tubes and failure at a stress level below the pretest stress. The material degradation was explained by a fiber oxidation model, which can be applied to various SiC cladding concepts. The embrittlement could be a limiting factor for using SiC cladding subjected to loss-of-coolant accident conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pressure-stabilized fixed-stress iterative solutions of compositional poromechanics

We consider the numerical behavior of the fixed-stress splitting method for coupled poromechanics as undrained regimes are approached. We explain that pressure stability is related to the splitting error of the scheme, not the fact that the discrete saddle point matrix never appears in the fixed-stress approach. This observation reconciles previous results regarding the pressure stability of the splitting method. Using examples of compositional poromechanics with application to geological CO sequestration, we see that solutions obtained using the fixed-stress scheme with a low order finite element-finite volume discretization which is not inherently inf-sup stable can exhibit the same pressure oscillations obtained with the corresponding fully implicit scheme. Moreover, pressure jump stabilization can effectively remove these spurious oscillations in the fixed-stress setting, while also improving the efficiency of the scheme in terms of the number of iterations required at every time step to reach convergence.

42 ENGINEERING↗
DOE OSTI2024

Pressure stability in explicitly coupled simulations of poromechanics with application to CO 2 sequestration

We study in detail the pressure stabilizing effects of the non-iterated fixed-stress splitting in poromechanical problems which are nearly undrained and incompressible. When applied in conjunction with a spatial discretization which does not satisfy the discrete inf–sup condition, namely a mixed piecewise linear–piecewise constant spatial discretization, the explicit fixed-stress scheme can have a pressure stabilizing effect in transient problems. This effect disappears, however, upon time step refinement or the attainment of steady state. The interpretation of the scheme as an Augmented Lagrangian method similar to Uzawa iteration for incompressible flow helps explain these results. Moreover, due to the slowly evolving solution within undrained seal regions, we show that the explicit fixed-stress scheme requires very large time steps to reveal its pressure stabilizing effect in examples of geologic CO 2 sequestration. We note that large time steps can result in large errors in drained regions, such as the aquifer or reservoir regions of these examples, and can prevent convergence of nonlinear solvers in the case of multiphase flows, which can make the explicit scheme an unreliable source of pressure stabilization. We conclude by demonstrating that pressure jump stabilization is as effective in the explicit fixed-stress setting as in the fully implicit setting for undrained problems, while maintaining the stability and convergence of the fixed-stress split for drained problems.

58 GEOSCIENCES↗