Search NASA⌕ Search

SEARCH · Search NASA

Results for “engineered matrix composites”

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 811 records · Page 45

High-Pressure Lightweight Thrusters

Returning samples of Martian soil and rock to Earth is of great interest to scientists. There were numerous studies to evaluate Mars Sample Return (MSR) mission architectures, technology needs, development plans, and requirements. The largest propulsion risk element of the MSR mission is the Mars Ascent Vehicle (MAV). Along with the baseline solid-propellant vehicle, liquid propellants have been considered. Similar requirements apply to other lander ascent engines and reaction control systems. The performance of current state-ofthe- art liquid propellant engines can be significantly improved by increasing both combustion temperature and pressure. Pump-fed propulsion is suggested for a single-stage bipropellant MAV. Achieving a 90-percent stage propellant fraction is thought to be possible on a 100-kg scale, including sufficient thrust for lifting off Mars. To increase the performance of storable bipropellant rocket engines, a high-pressure, lightweight combustion chamber was designed. Iridium liner electrodeposition was investigated on complex-shaped thrust chamber mandrels. Dense, uniform iridium liners were produced on chamber and cylindrical mandrels. Carbon/carbon composite (C/C) structures were braided over iridium-lined mandrels and densified by chemical vapor infiltration. Niobium deposition was evaluated for forming a metallic attachment flange on the carbon/ carbon structure. The new thrust chamber was designed to exceed state-of-the-art performance, and was manufactured with an 83-percent weight savings. High-performance C/Cs possess a unique set of properties that make them desirable materials for high-temperature structures used in rocket propulsion components, hypersonic vehicles, and aircraft brakes. In particular, more attention is focused on 3D braided C/Cs due to their mesh-work structure. Research on the properties of C/Cs has shown that the strength of composites is strongly affected by the fiber-matrix interfacial bonding, and that weakening interface realizes pseudo-plastic behavior with significant increase in the tensile strength. The investigation of high-temperature strength of C/Cs under high-rate heating (critical for thrust chambers) shows that tensile and compression strength increases from 70 MPa at room temperature to 110 MPa at 1,773 K, and up to 125 MPa at 2,473 K. Despite these unique properties, the use of C/Cs is limited by its high oxidation rate at elevated temperatures. Lining carbon/carbon chambers with a thin layer of iridium or iridium and rhenium is an innovative way to use proven refractory metals and provide the oxidation barrier necessary to enable the use of carbon/ carbon composites. Due to the lower density of C/Cs as compared to SiC/SiC composites, an iridium liner can be added to the C/C structure and still be below the overall thruster weight. Weight calculations show that C/C, C/C with 50 microns of Ir, and C/C with 100 microns of Ir are of less weight than alternative materials for the same construction.

Holmes, Richard↗

Evaluation of Equiatomic CrMnFeCoNiCu System and Subsequent Derivation of a Non-Equiatomic MnFeCoNiCu Alloy

Investigation into non-equiatomic high-entropy alloys has grown in recent years due to questions about the role of entropy stabilization in forming single-phase solid solutions. Non-equiatomic alloys have been shown to retain the outstanding mechanical properties exhibited by their equiatomic counterparts and even improve electrical, thermal, and magnetic properties, albeit with relaxed composition bounds. However, much remains to understand the processing–structure–property relationships in all classes of so-called high-entropy alloys (HEAs). Here, we are motivated by the natural phenomena of crystal growth and equilibrium conditions to introduce a method of HEA development where controlled processing conditions determine the most probable and stable composition. This is demonstrated by cooling an equiatomic CrMnFeCoNiCu alloy from the melt steadily over 3 days (cooling rate ~4 °C/h). The result is an alloy containing large Cr-rich precipitates and an almost Cr-free matrix exhibiting compositions within the MnFeCoNiCu system (with trace amounts of Cr). From this juncture, it is argued that the most stable composition is within the CrMnFeCoNiCu system rather than the CrMnFeCoNi system. With further optimization and evaluation, a unique non-equiatomic alloy, Mn 17 Fe 21 Co 24 Ni 24 Cu 14 , is derived. The alloy solidifies and recrystallizes into a single-phase face-centered cubic (FCC) polycrystal. In addition to possible applications where Invar is currently utilized, this alloy can be used in fundamental studies that contrast its behavior with its equiatomic counterpart and shed light on the development of HEAs.

36 MATERIALS SCIENCE↗

Liquid Interfacial Electron Microscopy Identifies Nanogalvanic Corrosion in Pearlitic Steel

The nanoscale mechanisms of localized corrosion in low carbon steels have remained elusive due to the complexity of studying the degradative material behavior at nanoscale solid-liquid interfaces. We identified various steps in the nanogalvanic corrosion processes using in-situ liquid-cell scanning transmission electron microscopy (STEM) using a microfluidic holder by Hummingbird Scientific. Initial work, performed at low magnification, identified the initiation point on a 1018 low-carbon steel surface. This initiation point was determined to be a triple junction of two ferrite grains bridging a cementite grain in contact with a baseline electrolyte of 6 uM CO2 dissolved in a buffered (2.78 uM Na2SO4) aqueous solution, pH 6.1. The pre-etched low-carbon steel surface was prepared using focused ion beam lift-out procedures to extract a cross-section of the low-carbon steel surface, which then was thinned to about 150 nm and transferred to a SiN membrane microfluidic window. The transfer was made using a lift-out needle to attach the low-carbon steel lamella to the corner of the SiN window, and then Pt/C deposition held the lamella in contact with the window while it was released from the lift out needle. To identify the triple point on the low carbon steel lamella, prior to attachment on the SiN window, the sample was characterized for compositional variations with energy dispersive x-ray spectroscopy mapping, grain orientation and phase mapping with precession electron diffraction, and thickness mapping with energy filtered transmission electron microscopy. This pre-characterization prior to the in-situ experiment provided a map of the multiphase and multigrain structure, where the in-situ liquid cell imaging provided a clear understanding of the initiation point on the sample. These data were cross-correlated to paint a holistic picture of the triple junction site, enabling low electron-fluence in-situ snapshot imaging to avoid dominating the native corrosion reactions with effects from the incident electron beam. This initial result identified that localized, nanogalvanic corrosion at the phase interface was the dominant corrosion process in the low-carbon steel, so we next targeted the observation of an array of these nanogalvanic features phase boundaries in a pearlite grain. Near-surface ferrite/cementite phase interfaces that typify pearlitic low-carbon steel were extracted, pre-characterized, and imaged for the in-situ corrosion processes. The sample was a cross-section from a pearlite grain, with alternating ferrite and cementite grains that extended microns down from the pre-etched low-carbon steel pipe surface. After contact with a buffered aqueous solution, the phase boundaries between the ferrite and cementite began to dissolve, with observable material loss and thickness changes in the dark-field and bright-field STEM images. Within minutes, the corrosion front proceeded deeper into the material, claiming a thin layer of ferrite around all exposed phase boundaries before progressing laterally into the ferrite matrix, converting the ferrite to corrosion product normal to each buried cementite grain. Formation of the corrosion product causes a volumetric expansion, creating a lateral wedging force that mechanically ejects the cementite grains from their grooves and leaves behind percolation channels into the steel substructure. Rapid and deleterious, this nanogalvanic corrosion pathway represents an important target for understanding and preventing run-away degradation in this common building material. Observation of this corrosion mechanism was enabled by the combination of pre-characterization using standard structural, grain, and compositional analysis in the TEM, which provides maps for understanding the reaction propagation captured in low-dose, in-situ, liquid-cell STEM.

corrosion↗

Scalability of Cohesive Fatigue Analyses Using Explicit Solvers

A cohesive fatigue law has been integrated into a constitutive material model compatible with an explicit finite element solver. The cohesive fatigue model response is based on engineering approximations of the endurance limit and the Goodman diagram. This approach can predict stress-life diagrams for crack initiation, the Paris law regime, and transient effects of crack initiation and stable tearing. Simplified cyclic loading is utilized so that the applied load(or displacement) corresponds to the peak load of a fatigue cycle. Loads are held constant during fatigue while damage develops with increasing solution increments. An automatically-calculated ratio of fatigue cycles per solution increment controls the rate of damage growth, ensuring that damage growth is modeled with a sufficient minimum number of increments and damage growth advances to a minimum desired extent within the explicit analysis step time. The compatibility with an explicit finite element solver enables the analysis of structures that are computationally intractable for implicit finite element solvers. Scalability studies are conducted for geometrically nonlinear problems involving fiber-reinforced composite structures that exhibit fatigue damage growth of interacting matrix cracks and delaminations

Frank A Leone↗

Humans to Mars, but How Many? Using Training Requirements Modeling to Inform Crew Size

Missions to Mars will differ from previous human spaceflight missions in that the onboard crew of astronauts will be required to operate in an Earth-independent manner due to the long communication delays. Without a systematic, repeatable process to determine the number and composition of crew necessary to successfully accomplish these missions, NASA increases the risk that crew sizes may be too small to meet primary mission objectives under nominal conditions and, more consequentially, that crewmembers may not have the expertise needed to successfully respond to unforeseen failures without the real-time expertise of the Mission Control Central (MCC) team NASA currently relies upon. The NASA Engineering and Safety Center (NESC) is developing a methodology for assessing the trade space of factors that affect the number of crew for future missions. This methodology includes the consideration of results from three human performance models developed using the Improved Performance Research and Integration Tool (IMPRINT) modeling platform as well as a custom-built model on expertise trained within the crew. The IMPRINT results will be presented in the modeling and simulation sub-tag. Here we present results of a model based on NASA’s crew qualification and responsibility matrix (CQRM), a tool used to identify the crew qualifications for each area of responsibility (operation, system, and payload) for a mission. The model outputs an optimized allocation of training assignments along with a flight-assigned CQRM that can be used to consider the expertise that can be trained within a crew of a given size. We discuss the CQRM model result implications on the trade space for Mars mission crew size.

Mars↗

Humans to Mars, but How Many? Using Training Requirements Modeling to Inform Crew Size

Missions to Mars will differ from previous human spaceflight missions in that the onboard crew of astronauts will be required to operate in an Earth-independent manner due to the long communication delays. Without a systematic, repeatable process to determine the number and composition of crew necessary to successfully accomplish these missions, NASA increases the risk that crew sizes may be too small to meet primary mission objectives under nominal conditions and, more consequentially, that crewmembers may not have the expertise needed to successfully respond to unforeseen failures without the real-time expertise of the Mission Control Central (MCC) team NASA currently relies upon. The NASA Engineering and Safety Center (NESC) is developing a methodology for assessing the trade space of factors that affect the number of crew for future missions. This methodology includes the consideration of results from three human performance models developed using the Improved Performance Research and Integration Tool (IMPRINT) modeling platform as well as a custom-built model on expertise trained within the crew. The IMPRINT results will be presented in the modeling and simulation sub-tag. Here we present results of a model based on NASA’s crew qualification and responsibility matrix (CQRM), a tool used to identify the crew qualifications for each area of responsibility (operation, system, and payload) for a mission. The model outputs an optimized allocation of training assignments along with a flight-assigned CQRM that can be used to consider the expertise that can be trained within a crew of a given size. We discuss the CQRM model result implications on the trade space for Mars mission crew size.

Mars↗

Tunable PA6 Polymer System for Thermoplastic Reinforced Body Panels

Project Goal: Overall, a redesigned polyamide system will combine the ductility and processing of a PA6 matrix, the strength and modulus of CF, and the facile recyclability of esters. Approach: PEA can be synthesized under conventional, industrially relevant pathways. Material properties and processing will be analyzed and compared to commercial PA6 material. Results: Polycondensation, traditional for polyesters, was successful in synthesized a PEA. PEA exhibited fast crystallization and bimodal melting behavior suggesting beneficial crystallization kinetics. PEA exhibit reduced melting temperatures owed to the incorporation of ester units in the polymer structure.

42 ENGINEERING↗

Advanced Multi-Tube Mixer Combustion for 65% Efficiency (Final Report)

This project targeted advanced low NOx combustion for advanced gas turbines capable of 65%, or greater, efficiency in combined cycle application. This technology advancement has further potential to benefit gas turbines used in coal based IGCC applications with pre-combustion carbon capture and hydrogen as the resulting fuel. The program developed and synthesized the most advanced combustion system capable of achieving low NOX emissions up to turbine inlet temperatures of 3100F while also supporting the load-following needs of a modern grid. The combustion system contributes to the overall gas turbine efficiency goal by setting the maximum cycle temperature achievable for a given NOX level and by minimizing the through-combustor air flow pressure drop. Focus areas for this project targeted maximizing the turbine inlet temperature entitlement, as constrained by emissions considerations. The design also minimized parasitic air flow pressure drop by using advanced cooling techniques and performance materials selections and by minimizing hot surface area. These two technology objectives (maximum, emissions-compliant cycle temperature and minimum air flow pressure drop) were integrated into a prototype design. The primarily analytical project sought to identify the most promising technologies to meet these objectives. Additional critical “jugular” data were obtained from multi-tube mixer tests to realize the potential of leveraging “micro flames” for minimizing overall hot surface area. This data was used, in conjunction with an understanding of advanced material and cooling design technologies, to analytically develop multiple design concepts. Phase I focused on in-depth engineering analysis and design, with minimal supporting laboratory testing to enable a selection of the top three combustion architectures for achieving these overall objectives. Phase II of the program developed the selected design through a combination of sub-scale testing and analytical efforts. Early tests included a cold-flow cascade to establish aerodynamic performance characteristics and a sub-scale fired test at GE Global Research in Niskayuna, NY, to establish cooling and heat transfer characteristics in conjunction with combustion performance. The data from these tests validated the analytical models to ultimately design a full-scale, test article to evaluate at prototypical pressure and temperature conditions at GE Gas Power’s Gas Turbine Technology Laboratory in Greenville, SC. GE Gas Power also developed, tested, and recommended a suitable seal design to be applied to the unique features of the combustor. To assess the technology challenges from prospective future production of the combustor from a ceramic matrix composite material, screening tests of Environmental Barrier Coatings were completed.

20 FOSSIL-FUELED POWER PLANTS↗

Coal-derived graphene foam and micron-sized silicon composite anodes for lithium-ion batteries

Silicon-based materials demonstrate significant potential as lithium-ion batteries (LIBs) anode, but their expansion and degradation present engineering design challenges for commercial application. In this study, a porous three-dimensional (3D) graphene and micron-sized silicon composite anode (Si@G foam) was synthesized using humic acid (HA) derived from coal as a graphene precursor. In-situ formation of graphene structure through reducing HA was confirmed by Raman spectra. The reduced HA (rHA) shows similar electrical conductivity compared to the commercial conductive carbon. SEM images depict a 3D skeleton of coal-derived graphene, with silicon particles distributed on the 3D graphene foam's internal surface. The Si@G composite-anode displays a good reversible capacity of ~656 mAh/g at a current density of 50 mA/g, as well as a high-rate capability of ~433 mAh/g at a current density of 800 mA/g, and outstanding cycling stability –89.8% capacity retention after 300 cycles, which is significantly higher than that of other foam structures. During lithiation and de-lithiation, the graphene foam serves as a matrix of electrical conductors and a volume expansion support for silicon. This 3D graphene network will be beneficial for developing advanced silicon-based anodes for high-performance LIBs.

25 ENERGY STORAGE↗

Particle‐Size‐Dependent Lithium‐Ion Transport in PEO/LLZO Composite Electrolytes

Lithium-metal batteries with solid electrolytes can deliver higher energy density and improved safety than conventional Li-ion batteries. Among solid electrolyte candidates, polymer/ceramic composite electrolytes are attractive because they combine polymer flexibility with the high ionic conductivity of ceramics. However, whether ceramic fillers synergistically reduce polarization losses in the polymer matrix remains unclear. A central unknown is the critical polymer/ceramic interfacial resistance (Rint,crit), below which adding ceramics lowers electrolyte overpotential. Here, we present the first macroscale model framework to quantify R int,crit for composite electrolytes based on polyethylene oxide (PEO) and Ta-doped Li 7 La 3 Zr 2 O 12 (LLZO). A 1D model for DC-polarization of tri-layer cells (PEO-LiTFSI/LLZO/PEO-LiTFSI) shows that LLZO surface functionalization reduces the PEO/LLZO interfacial resistance, consistent with electrochemical impedance measurements. Extending to a 2D composite model, we show notably that Rint,crit scales linearly with LLZO particle diameter and shifts toward experimentally accessible values (e.g., 28.8 Ωcm 2 ) as particle size increases. At fixed ceramic volume fraction, larger LLZO particles reduce the number of interfacial crossings, driving more current through the ceramic phase and lowering concentration polarization. In contrast, R int,crit is largely independent of ceramic volume fraction. These results demonstrate that ceramic filler-size engineering can enable synergistic, energy-efficient transport in polymer/ceramic composite electrolytes.

25 ENERGY STORAGE↗

Crafting Multifunctional Materials with Tailored Mechanical and Magnetic Properties by Solid-State Non-equilibrium Processing

Aluminum matrix composites hold promise for creating lightweight structural components with multifunctional properties. Nonetheless, achieving desired magnetic or electrical characteristics without compromising their fundamental mechanical properties presents a significant challenge. Here, in this study, we explore mechanical stir-based extrusion as a solid-state metalworking technique to develop fine-grained microstructures and highly refined metallic composites. Our study focuses on creating a SmCo 5 -reinforced aluminum matrix composite using solid stir extrusion (SSE) and assessing its microstructure, mechanical, and magnetic properties through a comprehensive suite of characterization tools, including x-ray diffraction, electron microscopy, and magnetometry. Our findings reveal that incorporating a small amount (3.76 wt.%) of SmCo 5 phase in a non-heat-treatable, non-magnetic aluminum alloy via SSE yields a composite that exhibits good hard magnetic characteristics with a large coercivity (H ci = 13.5 kOe) and improved mechanical properties. This study underscores non-equilibrium processing via solid-state high-speed stirring as a compelling method for crafting multifunctional materials with tailored mechanical and magnetic properties.

36 MATERIALS SCIENCE↗

Growth of Cd 0.9 Zn 0.1 Te 1-y Se y Single Crystals for Room Temperature Gamma-Ray Detection

Quaternary Cd 0.9 Zn 0.1 Te 1-y Sey (CZTS) single crystals, a novel room-temperature nuclear radiation detector semiconductor material, have been grown using a modified vertical Bridgman method (VBM) and a travelling heater method (THM). The percentage concentration of selenium in the VBM-grown crystal was 3% and that in the THM-grown crystal was 2%. While the THM Frisch collar detector (4.4 × 4.4 × 10.7 mm 3 ) produced a highly resolved pulse height spectra (PHS) with a resolution of ~1.06% for 662-keV gamma rays without any correction, the VBM-grown detector (10 × 10 × 2 mm 3 ) offered a high energy resolution of ~2% after application of a digital biparametric correction. The high-resolution performance of these detectors has been attributed to the addition of Se in the Cd 0.9 Zn 0.1 Te (CZT) matrix. Ab-initio calculations based on density functional theory (DFT) also confirmed that the addition of Se in the CZT matrix helps to reduce the formation of TeCd and the TeZn anti-sites. The VBM-grown crystals were characterized using powder x-ray diffraction (XRD) and energy dispersive x-ray spectroscopy (EDS). While XRD results revealed sharp diffraction peaks confirming the crystalline nature of the grown crystal, the EDS results confirmed the targeted stoichiometry of the elemental composition. The bulk resistivity of the grown crystal was calculated to be ~3 × 10 10 Ω-cm from current-voltage characteristics recorded at room temperature in a planar configuration, ensuring that the grown CZTS crystals have low dark current as required for detector-grade crystals.

42 ENGINEERING↗

The ATIC Long Duration Balloon Project

Long Duration Balloon (LDB) scientific experiments, launched to circumnavigate the south pole over Antarctica, have particular advantages compared to Shuttle or other Low Earth Orbit (LEO) missions in terms of cost, weight, scientific 'duty factor' and work force development. The Advanced Thin Ionization Calorimeter (ATIC) cosmic ray astrophysics experiment is a good example of a university-based project that takes full advantage of current LDB capability. The ATIC experiment is currently being prepared for its first LDB science flight that will investigate the charge composition and energy spectra of primary cosmic rays over the energy range from about 10(exp 10) to 10(exp 14) eV. The instrument is built around a fully active, Bismuth Germanate (BGO) ionization calorimeter to measure the energy deposited by the cascades formed by particles interacting in a thick carbon target. A highly segmented silicon matrix, located above the target, provides good incident charge resolution plus rejection of the 'backscattered' particles from the interaction. Trajectory reconstruction is based on the cascade profile in the BGO calorimeter, plus information from the three pairs of scintillator hodoscope layers in the target section above it. A full evaluation of the experiment was performed during a test flight occurring between 28 December 2000 and 13 January 2001 where ATIC was carried to an altitude of approx. 37 km above Antarctica by an approx. 850,000 cu m helium filled balloon for one circumnavigation of the continent. All systems behaved well, the detectors performed as expected, more than 43 gigabytes of engineering and cosmic ray event data was returned and these data are now undergoing preliminary data analysis. During the coming 2002-2003 Antarctica summer season, we are preparing for a ATIC science flight with approx. 15 to 30 days of continuous data collection in the near-space environment of LDB float altitudes.

Guzik, T. G.↗

A Microchannel Heat Exchanger Produced From a Metal Matrix Composite by Hybrid Laser Powder Bed Fusion and Inkjet Printing

Abstract This paper explores the production of an oxide dispersion strengthened (ODS) 304L stainless steel microchannel heat exchanger (HX) using a hybrid additive manufacturing process of laser powder bed fusion and inkjet printing. The study investigates the capabilities and economics of the hybrid inkjet-laser powder bed fusion (LPBF) process and evaluates the dimensional accuracy, functionality, and mechanical properties of the resulting ODS alloy. The effectiveness and pressure drop of the ODS heat exchangers produced by the hybrid LPBF tool are also determined. Results show that the inkjet-doped samples have a lower mean channel height with higher standard deviation than samples produced by LPBF alone. This is attributed to greater absorption of laser energy for the powder coated with the oxide precursor. The economic analysis shows that the hybrid process has a potential for reducing the unit cost of the heat exchanger based on cost modeling assumptions.

Engineering↗

Additive manufacturing of carbon fiber-reinforced thermoset composites via in-situ thermal curing

Fiber-reinforced polymer composites are lightweight structural materials widely used in the transportation and energy industries. Current approaches for the manufacture of composites require expensive tooling and long, energy-intensive processing, resulting in a high cost of manufacturing, limited design complexity, and low fabrication rates. Here, we report rapid, scalable, and energy-efficient additive manufacturing of fiber-reinforced thermoset composites, while eliminating the need for tooling or molds. Use of a thermoresponsive thermoset resin as the matrix of composites and localized, remote heating of carbon fiber reinforcements via photothermal conversion enables rapid, in-situ curing of composites without further post-processing. Rapid curing and phase transformation of the matrix thermoset, from a liquid or viscous resin to a rigid polymer, immediately upon deposition by a robotic platform, allows for the high-fidelity, freeform manufacturing of discontinuous and continuous fiber-reinforced composites without using sacrificial support materials. This method is applicable to a variety of industries and will enable rapid and scalable manufacture of composite parts and tooling as well as on-demand repair of composite structures.

36 MATERIALS SCIENCE↗

Probabilistic Micromechanics and Macromechanics for Ceramic Matrix Composites

The properties of ceramic matrix composites (CMC's) are known to display a considerable amount of scatter due to variations in fiber/matrix properties, interphase properties, interphase bonding, amount of matrix voids, and many geometry- or fabrication-related parameters, such as ply thickness and ply orientation. This paper summarizes preliminary studies in which formal probabilistic descriptions of the material-behavior- and fabrication-related parameters were incorporated into micromechanics and macromechanics for CMC'S. In this process two existing methodologies, namely CMC micromechanics and macromechanics analysis and a fast probability integration (FPI) technique are synergistically coupled to obtain the probabilistic composite behavior or response. Preliminary results in the form of cumulative probability distributions and information on the probability sensitivities of the response to primitive variables for a unidirectional silicon carbide/reaction-bonded silicon nitride (SiC/RBSN) CMC are presented. The cumulative distribution functions are computed for composite moduli, thermal expansion coefficients, thermal conductivities, and longitudinal tensile strength at room temperature. The variations in the constituent properties that directly affect these composite properties are accounted for via assumed probabilistic distributions. Collectively, the results show that the present technique provides valuable information about the composite properties and sensitivity factors, which is useful to design or test engineers. Furthermore, the present methodology is computationally more efficient than a standard Monte-Carlo simulation technique; and the agreement between the two solutions is excellent, as shown via select examples.

Murthy, Pappu L. N.↗

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↗