Search NASA⌕ Search

SEARCH · Search NASA

Results for “pressure dependence”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spiers Memorial Lecture: Theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our ability to model and predict the kinetics of such reactions has grown in leaps and bounds through the synergy of ever more sophisticated experimental studies with increasingly high level theoretical analyses. Herein, this review begins with a brief historical overview of the progress from the Lindemann mechanism to the master equation, which now provides the cornerstone for most theoretical analyses. The current status of ab initio transition state theory based implementations of the master equation is then reviewed, beginning with a discussion of the energy resolved chemical conversion rates, followed by a review of the pressure dependence of the thermal kinetics. The latter discussion focusses on the collisional energy and angular momentum transfer rates as well as the complexities of non-thermal effects and of multiple-well multiple-channel reactions. The synergy with recent state-of-the-art experiments is used to motivate these discussions. Attempts to automate the calculations are also briefly reviewed. Throughout the discussion, we provide our perspective on current understanding and continuing challenges and opportunities. One key challenge relates to the coupling of sequences of reactions, which frequently leads to deviations from the classic presumption of thermal reactants.

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

The role of collisional energy transfer in the thermal and prompt dissociation of 1-methyl allyl

Collisional energy transfer plays a pivotal role in a-priori calculations of the pressure-dependent rate constants obtained from a master equation. However, accurate determinations of collisional energy transfer parameters are rare and so most kinetic studies rely on best-fits of these parameters to experimental measurements (if available) or estimates guided from literature studies. In this work, we have quantified the effect of the uncertainty in energy transfer parameters on the thermal and prompt dissociation kinetics of a resonance stabilized radical, 1-methyl allyl (1MA), of relevance to the combustion of 1- and 2-butene isomers. Simulations using literature kinetics models were performed to assess the impact of these uncertainties on flame propagation and speciation data in laminar flames of 1- and 2-butene. Analyses of the uncertainty propagated by the energy transfer parameters for 1MA dissociation to the flame simulations, in particular the laminar flame speed, indicate an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Ab-initio based theoretical calculations were also performed to obtain pressure-dependent kinetics for the reaction of 1MA with H-atoms. Lastly, theoretically calculated energy transfer parameters were used to best characterize the kinetics and branching between the chain propagating 1MA dissociation and the chain-terminating reaction 1MA + H.

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

Oxidation and pyrolysis of methyl propyl ether

The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.

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

Probing the high-pressure densification of amorphous silica nanomaterials using SBA-15: An investigation into the paradoxical nature of the first sharp diffraction peak

The densification and X-ray scattering of mesoporous silica (SBA-15) were measured simultaneously under gigapascal (GPa) pressures. The results are compared to previous work on amorphous silica (aSiO 2 ) and demonstrate the feasibility of measuring the densification of aSiO 2 nanomaterials with small angle X-ray scattering (SAXS) in-situ in a diamond anvil cell. Compared to fused silica, the position of the SBA-15 first sharp diffraction peak (FSDP) is 7 times more sensitive to pressure and has a transition in its pressure dependance at a lower pressure (~2 GPa vs.~13 GPa). SBA-15 has two densification regimes, low-density amorphous and highdensity amorphous, which have equations of state comparable to low-density amorphous and high-density amorphous fused silica. Further, the transition between these two regimes occurs at a lower pressure than for fused silica (~1.5 GPa vs.~13 GPa). The results suggest that there is no direct relationship between the FSDP position and the aSiO 2 density during compression.

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

X-Ray Diffraction and Electron Microscopy Studies of the Size Effects on Pressure-Induced Phase Transitions in CdS Nanocrystals

In recent years, investigations of the phase transition behavior of semiconducting nanoparticles under high pressure has attracted increasing attention due to their potential applications in sensors, electronics, and optics. However, current understanding of how the size of nanoparticles influences this pressure-dependent property is somewhat lacking. In particular, phase behaviors of semiconducting CdS nanoparticles under high pressure have not been extensively reported. Therefore, in this work, CdS nanoparticles of different sizes are used as a model system to investigate particle size effects on high-pressure-induced phase transition behaviors. In particular, 7.5, 10.6, and 39.7 nm spherical CdS nanoparticles are synthesized and subjected to controlled high pressures up to 15 GPa in a diamond anvil cell. Analysis of all three nanoparticles using in-situ synchrotron wide-angle X-ray scattering (WAXS) data shows that phase transitions from wurtzite to rocksalt occur at higher pressures than for bulk material. Bulk modulus calculations not only show that the wurtzite CdS nanomaterial is more compressible than rocksalt, but also that the compressibility of CdS nanoparticles depends on their particle size. Furthermore, sintering of spherical nanoparticles into nanorods was observed for the 7.5 nm CdS nanoparticles. Our results provide new insights into the fundamental properties of nanoparticles under high pressure that will inform designs of new nanomaterial structures for emerging applications.

Meng, Lingyao↗
DOE OSTI2022

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗
DOE OSTI2026

A study of silver acetate under extreme conditions

With the aim of exploring chemical systems that may undergo metallization when irradiated with hard X-rays, we selected silver acetate (AgC 2 H 3 O 2 ) for the subject of this study. X-ray-induced decomposition of silver acetate under ambient and high-pressure conditions was observed in a diamond anvil cell (DAC), leading to the formation of metallic nanograins of silver at ambient pressure and 1.65 GPa. At 4 GPa, no decomposition was observed. The Avrami kinetics equation also provides information about novel structural formation at ambient pressure and 1.65 GPa. By modeling of the XRD data, it was found that the size of the silver nanocrystallites formed at 1.65 GPa pressure steadily increased to ∼5 nm after 600 min of X-ray irradiation as determined by applying the Scherrer equation to the diffraction peak widths. Time-resolved X-ray diffraction (XRD) revealed pressure-dependent kinetics, demonstrating that coupling pressure with irradiation enables controlled photochemical pathways in this model system. Concurrent with previous studies, the application of high pressure (HP) can be considered as a means of controlling X-ray synthetic photochemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-material hydrodynamics with algebraic sharp interface capturing

A finite volume method for Eulerian multi-material hydrodynamics with sharp interface capturing is presented here. The pressure-temperature non-equilibrium multi-material system with finite-rate pressure relaxation in mixed-cells is considered here. This pressure closure facilitates material-property-dependent pressure relaxation, rather than instantaneous pressure equilibration, which in turn allows the use of unsplit high-order time-integrators. A modified tangent of hyperbola for interface capturing (THINC) method is used to reconstruct multi-material (> 2) interfaces, on three-dimensional unstructured meshes. A simple modification which extends the THINC reconstruction to interfaces between more than two materials is proposed. It is demonstrated that the modified THINC can capture multi-material interfaces within 2–4 tetrahedral cells. Since no geometric reconstructions are required by the THINC method, the presented multi-material method is algorithmically simple, and computationally efficient. Consistent reconstructions of conserved quantities at material interfaces ensure that conservation and closure laws are satisfied at the discrete level. Through a suite of test problems solved on unstructured meshes, it is demonstrated that the presented method is a promising candidate for accurate and efficient multi-material hydrodynamics computations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Pressure-induced suppression of charge density phases across the entire rare-earth tritellurides by optical spectroscopy

The rare-earth tritellurides (RTe 3 ) are a distinct class of 2D layered materials that recently gained significant attention due to hosting such quantum collective phenomena as superconductivity or charge density waves (CDWs). Many members of this van der Waals (vdW) family crystals exhibit CDW behavior at room temperature, i.e., RTe 3 compound where R = La, Ce, Pr, Nd, Sm, Gd, and Tb. Here, our systematic studies establish the CDW properties of RTe 3 when the vdW spacing/interaction strength between adjacent RTe 3 layers is engineered under extreme hydrostatic pressures. Using a non-destructive spectroscopy technique, pressure-dependent Raman studies first establish the pressure coefficients of phonon and CDW amplitude modes for a variety of RTe 3 materials, including LaTe 3 , CeTe 3 , PrTe 3 , NdTe 3 , SmTe 3 , GdTe 3 , and TbTe 3 . Results further show that the CDW phase is eventually suppressed at high pressures when the interlayer spacing is reduced and interaction strength is increased. Comparison between different RTe 3 materials shows that LaTe 3 with the largest thermodynamic equilibrium interlayer spacing (smallest chemical pressure) exhibits the most stable CDW phases at high pressures. In contrast, CDW phases in late RTe 3 systems with the largest internal chemical pressures are suppressed easily with applied pressure. Overall results provide comprehensive insights into the CDW response of the entire RTe 3 series under extreme pressures, offering an understanding of CDW formation/engineering in a unique class of vdW RTe 3 material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Effect of Pressure on NO x Entitlement and Reaction Timescales in a Premixed Axial Jet-In-Crossflow

This paper investigates the pressure dependency of a lean premixed jet injected into a lean vitiated crossflow with an experimentally verified detailed chemistry computational fluid dynamics (CFD) model and 53 species considered. Experimental data were taken in an axially staged combustor with an optically accessible test section, allowing the use of particle image velocimetry (PIV) and CH* chemiluminescence techniques as well as point measurement of species concentration, temperature, and pressure. The experimental data cases at one, three, and five atmospheres were selected to describe the flame stabilization dependency on pressure and gain the required knowledge for an extrapolation to engine condition. Simulated exit nitrogen oxide levels were validated with experimental emission data, and a global emission trend for the NO reduction at elevated pressure and constant turbine inlet temperature level was defined. The nitrogen oxide benefit at elevated operating pressure was justified with the significantly smaller flame surface area: the analysis of the simulated spanwise and top-view profiles showed a relatively short receded core flame with nitrogen oxide production in the center at high pressure relative to a longer and larger shear layer flame at atmospheric condition that produced NO toward the inner and outer side of the flame. Furthermore, decomposition of the Damköhler number revealed the strong influence of the reaction timescales with higher reaction rates at elevated pressure, along with a moderate influence of the turbulent timescales, showing higher turbulence intensity in the lee-side recirculation zone at lower pressure.

03 NATURAL GAS↗
DOE OSTI2021

Computational Fluid Dynamics Simulation of Compressible Non-Newtonian Biomass in a Compression-Screw Feeder

Compression-screw feeders play a critical role in biorefineries to transport lignocellulosic biomass feedstocks from the feeding hoppers to biomass-conversion reactors in order to pretreat and convert biomass to hydrocarbon liquid biofuels and other power and energy resources. One of the main challenges in the operation of screw feeder is plugging and jamming of compressed biomass with high concentration of insoluble solids. The focus of this paper is to numerically investigate the screw feeder at these challenging operating conditions and help with the optimization of the screw feeder design to avoid operation failure. In this work a customized CFD model based on open-source OpenFOAM package [1] was developed to simulate the concentrated biomass as a highly viscous non-Newtonian fluid in the screw feeder. The biomass is modeled as a single-phase compressible Bingham fluid with a plastic viscosity as well as a density-dependent yield stress. The compressibility formulation (pressure-dependent density) and the density-dependent yield stress formulation in the governing equations follow the suggestions from a recent study by Duncan et al. [2]. A pilot-scale hopper/screw feeding system at NREL [3] is used to compare the experimental observations with our simulation results. The auger is 280 mm long and tapered with outer diameter changing from 80 mm to 35 mm. The auger rotates from 10 to 60 rpm in a conical throat which contains anti-rotational bars. The simulations predicted the required torque for the screw feeder and the pressure increase at the exit for biomass feedstocks with various fluid viscosity properties and auger rotating speeds. The analysis of the stress forces helped to identify the critical conditions were the screw feeder excessive wear or jamming could occur.

biofuels↗
DOE OSTI2026

Hall effect in LaCrGe 3 at pressures up to 24 kbar and a 𝐻−𝑇 phase diagram of the pressure-induced magnetic phase

LaCrGe 3 is a metallic, itinerant ferromagnet with 𝑇 𝐶 ≈ 86K and an easy magnetic direction along its hexagonal 𝑐 axis. Here we present studies of the 𝐻 ∥ 𝑐 Hall effect in LaCrGe 3 at ambient pressure and under pressures up to 24 kbar. At ambient pressure the anomalous Hall conductivity in the ferromagnetic state reaches a large value of ∼520 Ω − 1 cm −1 . The overall data are consistent with the behavior expected for a ferromagnet with magnetic field applied along its easy axis. Measurements under pressure yield an 𝐻 = 0, 𝑃−𝑇 phase diagram consistent with earlier publications with 𝑇 𝐶 being suppressed with pressure and nonferromagnetic phases appearing for 𝑃 >15 kbar. The data suggest that the observed 𝑇 𝐶 suppression is due to the pressure-dependent exchange interactions rather than suppression of the size of the magnetic moment. Field-dependent Hall resistivity in the higher-pressure, nonferromagnetic state is consistent with the net moment being zero for 𝐻=0. In a finite field there is a sharp anomaly comparable to the one associated with a spin-flip metamagnetic transition. The 𝐻−𝑇 phase diagrams corresponding to a pressure-induced magnetic phase are constructed for several pressures. In conclusion, the Hall data suggest the possibility of the pressure-induced nonferromagnetic phase being a combination of multiple small regions of antiferromagnetic phases with different wave vectors.

Bud'ko, Sergey L. [Ames Laboratory, and Iowa State↗
DOE OSTI2021

Dependency of active pressure and equation of state on stiffness of wall

Autonomous motion and motility are hallmarks of active matter. Active agents, such as biological cells and synthetic colloidal particles, consume internal energy or extract energy from the environment to generate self-propulsion and locomotion. These systems are persistently out of equilibrium due to continuous energy consumption. It is known that pressure is not always a state function for generic active matter. Torque interaction between active constituents and confinement renders the pressure of the system a boundary-dependent property. The mechanical pressure of anisotropic active particles depends on their microscopic interactions with a solid wall. Using self-propelled dumbbells confined by solid walls as a model system, we perform numerical simulations to explore how variations in the wall stiffness influence the mechanical pressure of dry active matter. In contrast to previous findings, we find that mechanical pressure can be independent of the interaction of anisotropic active particles with walls, even in the presence of intrinsic torque interaction. Particularly, the dependency of pressure on the wall stiffness vanishes when the stiffness is above a critical level. In such a limit, the dynamics of dumbbells near the walls are randomized due to the large torque experienced by the dumbbells, leading to the recovery of pressure as a state variable of density.

42 ENGINEERING↗
DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗
DOE OSTI2022

Pressure-induced changes of valence fluctuation in β - YbAlB 4 probed by x-ray absorption spectroscopy

Here we report detailed valence fluctuation phenomena in the unique quantum-critical compound β-YbAlB 4 under pressures and at low temperatures. We directly observed the drastic change in the pressure dependence of the mean Yb valence at about 3 GPa with no pressure-induced structural transition. Below 3 GPa, the anomaly in the temperature dependences of the mean Yb valence was observed at about 55 K under ambient pressure, corresponding to a previously observed reduced Kondo coherence temperature, and this temperature rapidly decreases to 22 K at 2.2 GPa, in contrast to the monotonous decrease of the high Kondo temperature (~250K at ambient pressure). The present results reveal that this reduced Kondo coherence temperature is responsible for the unconventional quantum critical behaviors in β-YbAlB 4 . Above 3 GPa, it was found that the temperature dependence of the mean Yb valence exhibits a discontinuous change at low temperature (10.5 and 9 K under 4.8 and 5.5 GPa, respectively).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Pressure-Induced Phase Transitions in Bilayer La 3 Ni 2 O 7

La 3 Ni 2 O 7 differs from expectations in part because at ambient pressure, it exists in two polymorphs: an unconventional crystal structure with alternating layers of single- and triple-layered nickel−oxygen octahedra (i.e., La 2 NiO 4 plus La 4 Ni 3 O 10 , with the so-called 1313 crystal structure) and a classical double-layer Ruddlesden−Popper phase (i.e., La 3 Ni 2 O 7 , with the so-called 2222 crystal structure, which dominates at high pressure). In this study, we report the pressure-dependent structural and electrical resistive properties of single crystals of the classical double-layered La 3 Ni 2 O 7 , grown at slightly elevated pressure using the floating-zone method. Structural characterization under pressures up to 15.4 GPa reveals a gradual transition from orthorhombic to tetragonal symmetry that is completed at 12−14 GPa. Additionally, we present pressure and field-dependent electrical resistance measurements up to 27.4 GPa, from which we construct a phase diagram. Although a resistive transition at 80 K is observed at high pressures, the resistive hallmarks of superconductivity are not.

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

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗