Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microcracking”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An Oxygen‐Scavenger Sulfide Coating Enabling Long‐Term Stable Nickel‐Rich Cathodes

Oxygen release is a major issue associated with layer-structured metal oxide cathodes in lithium batteries, which can further cause a series of problems, such as irreversible phase transition, microcracking, and electrolyte decomposition. Eventually, these issues jointly result in cell performance degradation and safety hazards. Thus, it is very significant to tackle oxygen release for achieving long-term stable cyclability, but very challenging. Although intensive efforts have been invested to date, there still lacks a feasible solution. In this study, nanoscale ZrS 2 coatings are applied on prefabricated LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) cathodes directly via atomic layer deposition (ALD). Very encouragingly, we reveal that this ALD-deposited conformal ZrS 2 nanocoating can serve as an exceptional oxygen scavenger and then convert into a stable sulfate (Zr(SO 4 ) 2 ) coating. Such an in situ conversion is very beneficial and effective for protecting the electrolyte from decomposition. In addition, the resultant Zr(SO 4 ) 2 coating further inhibits undesirable reactions, stabilizes the interface between NMC811 and the electrolyte, suppresses microcracking, mitigates transition metal dissolution, and maintains the structural stability of the NMC811 cathode. Consequently, the ZrS 2 -coated NMC811 cathode has demonstrated extraordinary performance. Thus, this study advances the understanding of interface engineering while paving a new technical pathway for commercializing NMC811 cathodes.

Nickel-rich cathodes↗
DOE OSTI2022

Multiscale investigation of the microstructural mechanisms driving ratchet growth in PBX 9502

The high explosive PBX 9502 undergoes irreversible expansion during thermal cycling (“ratchet growth”). Recent innovations in thermomechanical modeling via homogenization strategies are beginning to incorporate mesoscale information such as grain size, total porosity, and spatial distribution of voids and cracks. To generate a complete experimental data set to challenge and inform these models, PBX 9502 pellets were thermally cycled, cross-sectioned using ion polishing, and imaged in high resolution with scanning electron microscopy. Ratchet growth was found to drive expansion through microcracking. Microcracks were affected by agglomeration of crystals within the PBX. Virgin material showed greater ratchet growth than recycled material.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗
DOE OSTI2023

Microscale drivers and mechanisms of fracture in post-processed additively manufactured Ti–6Al– 4 V

Herein, we focus on understanding the microstructure-fracture correlations in a Ti–6Al– 4 V alloy additively manufactured via electron beam melting (EBM) and subjected to various post-process heat-treatments. Specifically, the as fabricated material is subjected to a sub-transus heat-treatment followed by air-cooling and a super-transus heat-treatment followed by either air- or furnace-cooling. Next, a series of in-situ single edge notch tension (SENT) tests are carried out under a high-resolution digital optical microscope. The panoramic high-resolution images captured during the in-situ tests are then used to characterize the planar deformation on the specimen surface using microstructure-based digital image correlation (DIC). Additionally, the results of the in-situ SENT tests together with DIC and post-mortem fractographic analyses provided us with a better understanding of the microstructure-fracture correlations in these materials. Our results show that the fracture mechanism of the as fabricated and sub-transus heat-treated materials is essentially the same, while the changes in the microstructure following the super-transus heat-treatments significantly affects the fracture mechanism. In this case, several microcracks of hundreds of microns in length first nucleate away from the deformed notch following extreme plastic deformation at discrete locations. Furthermore, the location of these microcracks in the super-transus heat-treated materials is extremely sensitive to the details of the underlying microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shrinkage-induced deformations and creep of structural concrete: 1-year measurements and numerical prediction

Highlights: • Extensive experimental study on drying shrinkage, creep and microcracking of concrete • All specimens prepared from a single batch of ordinary-strength structural concrete • 1st year of measurements of (not only) non-uniformly drying beams with span up to 3-m • The database is downloadable from free-to-use research data repository • Modified MPS model for concrete creep used in blind prediction of all experiments The material models for creep and shrinkage operating on the material point level in FEM are usually intended for challenging complex applications and structures, where the average cross-sectional approach does not suffice. The identification of the growing number of material parameters induced by increasing model capabilities relies on very specific and narrow-oriented yet interconnected experiments which are scarce. The presented comprehensive experiments aim to provide a clearer image of the complicated interaction among the basic phenomena: drying, shrinkage, creep, and microcracking. The cornerstone of this ongoing research is a unique set of 30 partially-sealed unreinforced concrete beams with span 1.75–3.0 m subjected to drying. To minimize material variation, all specimens in this study were cast from a single batch of ordinary strength structural concrete with slag-blended binder. The resulting experimental database will be suitable both for validation and development of the constitutive models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deterioration of concrete due to ASR: Experiments and multiscale modeling

The process of ASR (Alkali-Silica Reaction) induced expansion and damage in pavement concrete specimens is investigated using laboratory experiments and computational modeling. In the experimental program, the concrete specimens are subject to CS-CPT (climate simulation concrete prism test) to obtain ASR induced expansion with and without external supply of alkali. The dissolution rates of the granodiorite used in the concrete mix and the gel formation rates are determined under concrete-like conditions (pH 13.8, with/without Ca(OH){sub 2} and NaCl) at different temperatures. A micromechanics based computational model with aggregate-scale diffusion and reaction kinetics coupled to an Eigenstrain based micromechanics damage model is developed for the simulation of ASR induced expansion and damage. Data from the experimental program are used to calibrate and validate the computational model. Model predictions show that for the given concrete mixture, ASR induced expansion in the specimen exposed to water is predominantly governed by microcracking in the aggregate, while the expansion in the specimen subjected to external alkali supply is governed by microcracking in both the aggregates and the cement paste.

36 MATERIALS SCIENCE↗
DOE OSTI2021

DEM analysis of the effect of interface transition zone on dynamic splitting tensile behavior of high-strength concrete based on multi-phase model

Highlights: • Crushable aggregates with the real geometry shape are generated by using clump-cluster method. • Flat-Joint Model is adopted as the constitutive law. • The tensile strength is almost not affected by ITZ strength when its strength exceeds 0.8 times that of mortar. • A modified formula that can reflect the functional relationship between the ITZ strength and the TDIF is proposed. The tensile strength is one of the main parameters in the design of concrete structures and has an important influence on the cracking performance and durability of concrete materials. As the weak region in concrete, the strength of the interface transition zone (ITZ) directly affects the dynamic mechanical response of concrete. Meanwhile, considering that the tensile performance of concrete is far weaker than its compressive performance, the study of the strength of ITZ is of great significance to the dynamic tensile fracture characteristics and mechanical properties of concrete. In this paper, a method combining 3D scanning technology and Clump-Cluster was proposed to generate a real geometry shape model of crushable aggregate. On this basis, the 3D discrete element model (DEM) of high-strength concrete was constructed using the Flat-Joint Model (FJM). The effects of the ITZ strength on the failure mode, number of microcracks, tensile strength and strain rate effect of high-strength concrete under strain rates from 0.01/s to 100/s were studied by splitting tensile test. The simulation results indicate that the number of microcracks generated decreases with the increase of ITZ strength, and the proportion of ITZ cracks decreases obviously. When the ITZ strength exceeds the critical value of 0.8 times that of mortar, the tensile strength is almost not affected by ITZ strength. In addition, based on the numerical simulation results, a modified formula that can reflect the functional relationship between the ITZ strength and the dynamic increasing factor of the concrete tensile strength (TDIF) is proposed, and the rationality of the TDIF modified formula is verified.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accelerating high-strain continuum-scale brittle fracture simulations with machine learning

Failure in brittle materials under dynamic loading conditions is a result of the propagation and coalescence of microcracks. Simulating this discrete crack evolution at the continuum level is computationally expensive or, in some cases, intractable, resulting in the need to make broad assumptions or neglect key physics. In this work, we have developed an approach using machine learning that overcomes the current inability to represent meso-scale physics at the macro-scale. Our approach leverages damage and stress data from a computationally expensive high-fidelity model that explicitly resolves microcrack behavior to build an inexpensive machine learning emulator. Once trained, the machine learning emulator is used to predict the evolution of crack length statistics, which then informs a continuum-scale constitutive model. This results in a significant speed-up of the workflow by four orders of magnitude. Both the machine learning emulator and the continuum-scale model are validated against the high-fidelity model and experimental data, respectively, showing excellent agreement. There are two key findings. The first is that we can reduce the dimensionality of the problem, establishing that the machine learning emulator only needs the length of the longest crack and one of the maximum stress components to capture the necessary physics. Another compelling finding is that the emulator can be trained in one experimental setting and transferred successfully to predict behavior in a different setting.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of gas storage and transport in Eagle Ford shale using pressure pulse decay measurements with He, Kr and CO 2

Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.

58 GEOSCIENCES↗
DOE OSTI2025

Exceptional radiation resistance of hardened amorphous SiC under high-temperature hydrogen ion implantation

This study provides a compelling comparison of the structural and mechanical responses of single-crystal silicon carbide (sc-SiC), nanocrystalline silicon carbide (nc-SiC), and amorphous silicon carbide (am-SiC) to hydrogen ion implantation at 650 °C across three distinct fluences: low, medium, and high (LF, MF, HF). While both sc-SiC and nc-SiC exhibit blistering, microcracking, and exfoliation, am-SiC remains free of blisters, demonstrating superior resilience. Notably, nc-SiC, with its high density of stacking faults (SFs), requires a higher fluence to initiate blistering compared to sc-SiC. At high fluences, hydrogen accumulation at grain boundaries (GBs) in nc-SiC elevates internal gas pressure, intensifying exfoliation. In sc-SiC, blistering leads to increased hardness, whereas in nc-SiC, the degradation of the SF structure results in a reduction in hardness. In contrast, am-SiC undergoes structural relaxation during irradiation, resulting in a significant increase in hardness while maintaining its structural integrity, with only the formation of nano-sized spherical bubbles observed. Furthermore, these findings highlight the exceptional suitability of am-SiC for nuclear applications, where resistance to radiation-induced microcracking is critical.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of sample bias on adhesion of magnetron sputtered Cr coatings on SiC

Swelling of SiC at 300 C due to in-service neutron irradiation causes tensile residual stresses in coatings which are expected to adversely affect the performance of coated SiC composite fuel cladding for light water reactors. Matching the coating swelling with the substrate, a solution common for thermal expansion, is not practical in the case of neutron irradiation. Biasing samples during magnetron sputtering deposition induces compressive residual stress which may counteract this. In this study, chromium coatings were deposited on SiC by DC magnetron sputtering with no external heating at bias voltages of –50V, –75V, and –100V. The effects of the bias voltage on morphology, residual stress, microstrain, texture, and adhesion are shown. The low deposition temperature resulted in the coating microstructure evolution following an energetic particle bombardment dominated trend. Additionally, at the two lower bias voltages knock-on implantation dominated increasing the residual stress and microstrain while at the highest bias voltage, thermal spike migration allowed for defect relaxation. When the knock-on induced compressive residual stress exceeded 0.8 GPa microcrack formation in the SiC substrate decreased coating adhesion. While no microcracks formed at the lowest bias voltage, insufficient atomic mobility during coating growth lead to voids forming in the coating. A balance is needed to form void-free coatings that have high compressive residual stress.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of weave architecture on mechanical response of SiC f -SiC m tubular composites

Due to their high resistance to radiation damage and elevated temperature, silicon carbide fiber-reinforced, silicon carbide matrix composites (SiC f -SiC m ) are identified as potential cladding structures for use in nuclear reactors. In this study, four composite architectures with varying ply numbers along the thickness direction and different biaxial or triaxial plain weave orientations at either 45° or 60°, were systematically evaluated under various stress states to assess the influence of weave architecture on mechanical performance. Experiments were conducted on SiC f -SiC m composite tubes under tensile hoop, axial compression, and rotating flexural loading to evaluate the mechanical response and investigate the failure modes using high-speed imaging and digital image correlation (DIC) techniques. It was found that for tensile hoop burst and flexural loading, the braiding angle had the most significant influence on the strength of the composite, whereas the effect of fiber angle was more limited for compression testing. Under axial compression a unique failure mode where a microcrack nucleates and grows only to a length equal to the thickness of a single yarn was identified. This crack growth behavior is reflected as periodic oscillations in the load-displacement response. For both axial and hoop loading, regardless of weave angle and number of plies, failure always initiated parallel to the tube axis in a single yarn and the cumulative interaction of these microcracks lead to either axial burst or fracture at an angle to the tube axis along a yarn. Furthermore, these results point to the importance of customizing the design of tube architecture for enhanced performance in specified nuclear applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The comparative strength and fracture toughness properties of commercial 95W-3.5Ni1.5Fe and 95W-3.5Ni1.5Cu tungsten heavy alloys

Due to the sharp surface cracks resulting from high and cyclic thermal loads, fracture toughness is a crucial property for fusion divertor materials. Tungsten heavy alloys (WHA) have emerged as a leading candidate for plasma facing component applications. Both the WNiFe and WNiCu WHAs have been considered. However, there is no fracture toughness data in the literature on the WNiCu WHA. Therefore, maximum load elastic–plastic fracture toughness (K Jm ) of a commercial 95W-3.5Ni1.5Cu (NiCu) WHA plate is reported here for the first time. The NiCu WHA K Jm is compared to a corresponding K Jm for commercial 95W-3.5Ni1.5Fe (NiFe) WHA, along with comparisons of their respective microstructures, microhardness, and tensile properties. Although their room temperature yield strengths are similar (≈ 600 MPa), the NiFe WHA has higher ultimate stress (818 MPa vs. 642 MPa) and total elongation (8% vs. 1.2%) compared to the NiCu WHA. Most notably, the room temperature elastic–plastic K Jm for the NiFe WHA of ≈ 80 ± 8 MPa$\sqrt{m}$, is nearly twice that of the NiCu WHA at ≈ 47 ± 4 MPa$\sqrt{m}$. Further, while the NiFe WHA experienced stable ductile tearing, the crack propagated unstably (approximately elastically) after reaching the maximum load in the NiCu WHA. The higher K Jm of the NiFe WHA is attributed to its higher flow stress and ductility, along with the energy dissipation and microcrack dilatational stress shielding, provided by a wide distribution of arrested process-zone microcracks, which are largely absent in the NiCu WHA. Since, the 95W-NiFe WHA has distinctly better mechanical properties than 95W-NiCu WHA, it is a more suitable candidate for the fusion reactor application.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tribological behavior of spark plasma sintered Ti 3 SiC 2 MAX phase composites

Ti₃SiC₂ MAX phases are considered promising candidates for tribological applications; but, their low hardness (5-6 GPa) can lead to increased abrasive wear and higher wear rates. This study investigates the effect of incorporating hard SiC particles on microstructure, mechanical, and tribological behavior of Ti₃SiC₂-SiC-based MAX phase composites prepared using spark plasma sintering at 1400 °C, 40 MPa, 15-minute holding time. Phase and microstructural characterization of composites confirmed the formation of Ti₃SiC₂ MAX phase (90%) along with TiC as a minor phase (10%). For Ti₃SiC₂-SiC, homogeneous distribution of SiC grains within Ti₃SiC₂ matrix resulted in increased hardness from ~11.6 GPa to ~14.8 GPa; however, the flexural strength decreased from ~615 MPa to ~597 MPa due to coefficient of thermal expansion mismatch. Tribological behavior of Ti₃SiC₂-SiC MAX phase composites was assessed using unidirectional sliding ball-on-flat tests with a 52100 steel ball at different loads and a constant speed of 0.05 m/s. SiC reinforced composites exhibited a decreased friction coefficient from 0.39 to 0.29 and wear rates comparable to Ti₃SiC₂ (7-10×10⁻³ mm³/N·m). SEM and EDS analysis of the worn surfaces indicated material transfer from the steel ball counter body and its oxidation due to frictional heating, at lower loads (2N and 5N), whereas at higher load (10N), fracture of the transfer layer dominated, with the presence of microcracks, delamination, grain pull-outs, and wear debris. Improved mechanical properties and good adhesion of SiC with Ti₃SiC₂ matrix resulted in reduced microcracks and grain pull-outs, making Ti₃SiC₂-SiC/Steel tribo-pair more suitable for tribological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Material Characterization-Based Wear Mechanism Investigation for Biomass Hammer Mills

Biomass, as harvested, is composed of inorganic compounds both intrinsically and extrinsically and can be abrasive. The present study investigates the wear modes and mechanisms of two types of blades of hammer mills used in biomass size reduction (impacting the particle size and distribution) and densification (impacting the size, shape, and density). The dominant wear modes for the stage 1 steel blades are determined to be erosive and polishing wear. For the stage 2 blades with a carbide weld overlay, the main wear mechanisms are erosion and fracture. Partial replacement of Co by Fe in the outer layer of the carbide grits, likely induced by diffusion during high-temperature welding, has been correlated to the observed microcracking. Finally, the microcracking is believed to weaken the grit strength and fracture toughness to make the overlay prone to fracture and erosion due to repetitive contact with the inorganic contents in chopping biomass.

42 ENGINEERING↗
DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗
DOE OSTI2021

Damage Modeling of Solid Oxide Fuel Cells Accounting for Redox Effects

This work applies a multiscale mechanistic damage model developed for brittle ceramics and implemented in commercial finite element (FE) packages via user subroutines to study progressive damage in solid oxide fuel cells (SOFC) subjected to thermomechanical loading under normal operating and shutdown conditions including redox effects. The damage model captures the micromechanics of stiffness reduction due to material porosity change and microcracking and integrates the as-obtained stiffness reduction law into a continuum damage mechanics (CDM) formulation for the evolution of microcracks up to fracture. The volumetric “swelling” that occurs during redox is treated in constitutive modeling similarly to thermal expansion, but swelling strains are irreversible. Furthermore, this damage model was first validated through predictions of strength and stress-strain response for the SOFC electrode materials. Next, it has been applied to predict the potential for degradation in a generic planar SOFC stack with large active area cells. Multicell stack models were simulated in both co-flow and counter-flow configurations. In addition, a constant temperature redox cycle was also simulated to capture overall cell electrode damage due to volumetric swelling of the nickel (Ni)-based anode in the anode-supported cells.

25 ENERGY STORAGE↗