Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Volume”

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

Analysis of Degradation and Electrochemical Behavior of Lithium–Oxygen Batteries under Lean Electrolyte Conditions

Lithium-oxygen batteries have garnered much attention in the past decades due to their high energy density and active material abundancy. Decreasing the electrolyte volume is critical to achieving a high energy density for practical applications. Recent works have demonstrated a serious performance limitation at lower electrolyte fillings. Here, in this work, we study specifically the difference between flooded and lean electrolyte lithium-oxygen batteries by systematically tracing the source of performance degradation. We find that the key contributor to the decrease in performance under lean electrolyte conditions stems from the lithium metal anode rather than the cathode. Solid electrolyte interphase (SEI) compositionally appears to be the same from X-ray photoelectron spectroscopy but the coverage is less uniform, and impedance is much higher under lean electrolyte conditions. Such an effect likely produced the observed poorer cycle performance at 10 μL cm -2 (similar to 28% cathode's pore volume filled) vs 100 μL cm -2 (278% of cathode's pore volume) at a capacity of 1.0 mAh cm -2 (0.1 mAh cm -2 ) using an electrolyte composition of 1 M LiTFSI in TEGDME.

Córdoba, Daniel [Argonne National Laboratory (ANL)↗
DOE OSTI2024

Colloidal Stability of PFSA-Ionomer Dispersions. Part I. Single-Ion Electrostatic Interaction Potential Energies

Charged colloidal particles neutralized by a single counterion are increasingly important for many emerging technologies. Attention here is paid specifically to hydrogen fuel cells and water electrolyzers whose catalyst layers are manufactured from a perfluorinated sulfonic acid polymer (PFSA) suspended in aqueous/alcohol solutions. Partially dissolved PFSA aggregates, known collectively as ionomers, are stabilized by the electrostatic repulsion of overlapping diffuse double layers consisting of only protons dissociated from the suspended polymer. We denote such double layers containing no added electrolyte as "single ion". Size-distribution predictions build upon interparticle interaction potential energies from the Derjaguin-Landau-Verwey-Overbeek (DLVO) formalism. However, when only a single counterion is present in solution, classical DLVO electrostatic potential energies no longer apply. Accordingly, here a new formulation is proposed to describe how single-counterion diffuse double layers interact in colloidal suspensions. Part II (Srivastav, H.; Weber, A. Z.; Radke, C. J. Langmuir 2024 DOI: 10.1021/acs.langmuir.3c03904) of this contribution uses the new single-ion interaction energies to predict aggregated size distributions and the resulting solution pH of PFSA in mixtures of n-propanol and water. A single-counterion diffuse layer cannot reach an electrically neutral concentration far from a charged particle. Consequently, nowhere in the dispersion is the solvent neutral, and the diffuse layer emanating from one particle always experiences the presence of other particles (or walls). Thus, in addition to an intervening interparticle repulsive force, a backside osmotic force is always present. With this new construction, we establish that single-ion repulsive pair interaction energies are much larger than those of classical DLVO electrostatic potentials. The proposed single-ion electrostatic pair potential governs dramatic new dispersion behavior, including dispersions that are stable at a low volume fraction but unstable at a high volume fraction and finite volume-fraction dispersions that are unstable with fine particles but stable with coarse particles. Finally, the proposed single-counterion electrostatic pair potential provides a general expression for predicting colloidal behavior for any charged particle dispersion in ionizing solvents with no added electrolyte.

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

Electrolyte-phobic Surface for the Next-Generation Nanostructured Battery Electrodes

Nanostructured electrodes are among the most important candidates rationally designed to enable high capacity battery chemistry. Nanostructures can solve issues such as the volume change and mechanical fragmentations. However, the high surface area they usually possess would decrease the Coulombic efficiencies, since the side chemical reactions scale with surface area. Moreover, electrodes comprised of nanomaterials have significant intakes of liquid electrolytes, which reduces the overall energy density and increases the cost of the battery. Here we present a new strategy of limiting effective surface area by introducing an “electrolyte-phobic surface”. In this study, a porous Si anode functions as a model material to demonstrate this concept. Silicon boasts high theoretical capacity, but experiences large volume change during its lithiation and delithiation processes. Porous Si can address this volume change problem with the buffer effect of its inner pores. However, porous silicon shows low initial Coulombic efficiencies and high irreversible lithium loss, owing to its intrinsic high surface area. In this report, a covalently linked perfluorinated surface coating layer on porous Si particles serves as an electrolyte-phobic protection layer, minimizing the accessible surface area for the electrolytes, decreasing the side reactions between the electrolyte and Si surface, and thus significantly enhancing the initial Coulombic efficiencies, up to ~88% compared to ~60% for the pristine porous silicon. Meanwhile, the electrolyte-phobic protection layer of Si particles keeps the silicon surface compatible with the conventional polyvinylidene fluoride (PVDF) binder, which helps to stabilize the Si electrode for long-term battery cycling.

Qian, Chenxi↗
DOE OSTI2021

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

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

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion↗
DOE OSTI2025

Unique Conductivity Behavior in Water-In-Salt Electrolytes Driven by Ion Clusters

Understanding and predicting ion transport in aqueous electrolytes are crucial for advanced energy storage and biophysics, and many emergent technologies yet remain elusive. Herein, we introduce a unified framework to quantitatively describe and predict electrolyte conductivity that shifts from conventional molar concentration-based metrics to a volume fraction-based approach. Through analyzing a variety of electrolyte solutions via this perspective, we observe a universal conductivity peak at a 37% volume fraction. Small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations reveal that nanometer-scale ion clusters drive this general behavior. Moreover, key geometric features of the ion transport pathwayssuch as pore size, tortuosity, and connectivityfollow a consistent dependence with respect to the volume fraction, reinforcing the argument for the universal conductivity trend. This paradigm shift opens new avenues for designing high-performance electrolytes and provides transformative insights for advancing studies in many fields, wherein molecular aggregates dictate transport properties.

Nguyen, Huong T. D.↗
DOE OSTI2020

Collar Properties and Measurement Time Confer Minimal Bias Overall on Annual Soil Respiration Estimates in a Global Database

Abstract Measuring the soil‐to‐atmosphere carbon dioxide (CO 2 ) flux (soil respiration, R S ) is important to understanding terrestrial carbon balance and to forecasting climate change. Such measurements are frequently made using measurement collars permanently inserted into the soil surface. However, differences in measurement duration and frequency, as well as collar properties, may lead to biases in the estimation of annual R S . Using a newly updated global R S database (SRDB‐V5), we investigated the annual R S bias associated with five methodological factors: collar height, collar coverage area, collar insertion depth, measurement duration, and measurement frequency. We found that annual R S was negatively correlated with collar insertion depth, consistent with the idea that collar insertion cuts roots and thus reduces R S . Annual R S was also negatively related with collar height and collar coverage area, perhaps because uniform head‐space mixing is difficult to achieve in larger volume chambers; however, these effects were quantitatively small (bias of ~2% to 10% of mean R S ). We found no correlation of measurement duration or measurement frequency with annual R S . These findings suggest that variation in R S methodology generally introduces minimal bias overall. Therefore, compilations of minimally adjusted annual R S measurements provide a reliable resource for synthesis studies, global annual R S modeling, and investigation of how soil carbon responds to climate change.

Jian, Jinshi↗
DOE OSTI2023

Scaling of Floods With Geomorphologic Characteristics and Precipitation Variability Across the Conterminous United States

Abstract Accurate flood risk assessment requires a comprehensive understanding of flood sensitivity to regional drivers and climate factors. This paper presents the scaling of floods (duration, peak, volume) with geomorphologic characteristics of the basin (i.e., drainage area, slope, elevation) and precipitation patterns (rainfall accumulation, variability). Long‐term daily streamflow observations over the 20th and early 21st centuries from Hydro‐Climatic Data Network streamgages across the conterminous United States are used to create a flood event database based on their flood stage information. Antecedent daily rainfall accumulation and variability corresponding to these floods are computed using Global Historical Climatology Network daily data set. Two Bayesian scaling models are developed, and the spatial organization of scaling exponents is investigated. The baseline model quantifies the scaling of floods to geomorphologic characteristics. The dynamic model quantifies the scaling of floods to antecedent precipitation distribution which is further conditioned on geomorphologic characteristics. Results show that small and low‐elevation basins have a stronger response to antecedent rainfall distribution in amplifying flood peaks, while high‐elevation steeper basins have a lower response for flood duration and volume. The dynamic models demonstrate that there are significant variations in the flood scaling rates, with the largest rates up to 40% and 4.5% for flood duration, 64% and 44% for peak, and 98% and 40% for volume found across the Northeast, Coastal Southeast, and Northwest with intensifying rainfall accumulation and variability, respectively. This study advances flood predictions by better informing the flood attributes in the context of dynamical land‐atmosphere perturbations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Verifying Ray Tracing Amplitude Methods for Global Magnetospheric Modeling

Abstract Ray tracing is a commonly used method for modeling the propagation of electromagnetic waves in Earth's magnetosphere. To apply ray tracing results to global models of wave‐particle interaction such as energetic electron scattering, it is useful to map the discrete rays to a volume filling mesh. However, some methods have inherent losses of energy from the wave source, or do not account for the full range of wave properties within a sample volume. We have developed and tested a 3D magnetospheric ray tracing code “MESHRAY” which resolves these issues. MESHRAY uses the conservation of Poynting flux through ray triplets with finite volume to determine the local field amplitudes. Electromagnetic wave energy density from all ray data points is mapped to a mesh and verified against the wave source power for energy conservation varying time step length, number of rays, and total time steps. We find that the method is self‐consistent and numerically robust. We further investigate whether the neglect of phase information and superposition has a significant impact on the accuracy of mapping wave intensity to a mesh. We find excellent agreement between the analytic solution for waves emitted by a line source in a plane‐stratified medium and an equivalent ray tracing solution. When phase information is excluded, ray tracing reproduces an average amplitude spread over regions of coherent constructive and destructive interference. This may be an important consideration for interpolating ray tracing results of longer wavelength waves such as magnetosonic, electromagnetic ion cyclotron, or ULF waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Brain ageing in schizophrenia: evidence from 26 international cohorts via the ENIGMA Schizophrenia consortium

Schizophrenia (SZ) is associated with an increased risk of life-long cognitive impairments, age-related chronic disease, and premature mortality. We investigated evidence for advanced brain ageing in adult SZ patients, and whether this was associated with clinical characteristics in a prospective meta-analytic study conducted by the ENIGMA Schizophrenia Working Group. The study included data from 26 cohorts worldwide, with a total of 2803 SZ patients (mean age 34.2 years; range 18–72 years; 67% male) and 2598 healthy controls (mean age 33.8 years, range 18–73 years, 55% male). Brain-predicted age was individually estimated using a model trained on independent data based on 68 measures of cortical thickness and surface area, 7 subcortical volumes, lateral ventricular volumes and total intracranial volume, all derived from T1-weighted brain magnetic resonance imaging (MRI) scans. Deviations from a healthy brain ageing trajectory were assessed by the difference between brain-predicted age and chronological age (brain-predicted age difference [brain-PAD]). On average, SZ patients showed a higher brain-PAD of +3.55 years (95% CI: 2.91, 4.19; I 2 = 57.53%) compared to controls, after adjusting for age, sex and site (Cohen’s d = 0.48). Among SZ patients, brain-PAD was not associated with specific clinical characteristics (age of onset, duration of illness, symptom severity, or antipsychotic use and dose). This large-scale collaborative study suggests advanced structural brain ageing in SZ. Longitudinal studies of SZ and a range of mental and somatic health outcomes will help to further evaluate the clinical implications of increased brain-PAD and its ability to be influenced by interventions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Super-resolving microscopy images of Li-ion electrodes for fine-feature quantification using generative adversarial networks

Abstract For a deeper understanding of the functional behavior of energy materials, it is necessary to investigate their microstructure, e.g., via imaging techniques like scanning electron microscopy (SEM). However, active materials are often heterogeneous, necessitating quantification of features over large volumes to achieve representativity which often requires reduced resolution for large fields of view. Cracks within Li-ion electrode particles are an example of fine features, representative quantification of which requires large volumes of tens of particles. To overcome the trade-off between the imaged volume of the material and the resolution achieved, we deploy generative adversarial networks (GAN), namely SRGANs, to super-resolve SEM images of cracked cathode materials. A quantitative analysis indicates that SRGANs outperform various other networks for crack detection within aged cathode particles. This makes GANs viable for performing super-resolution on microscopy images for mitigating the trade-off between resolution and field of view, thus enabling representative quantification of fine features.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Topology-enhanced mechanical stability of swelling nanoporous electrodes

Abstract Materials like silicon and germanium offer a 10-fold improvement in charge capacity over conventional graphite anodes in lithium-ion batteries but experience a roughly threefold volume increase during lithiation, which challenges ensuring battery integrity. Nanoporous silicon, created by liquid-metal-dealloying, is a potentially attractive anode design to mitigate this challenge, exhibiting both higher capacity and extended cycle lifetimes. However, how nanoporous structures accommodate the large volume change is unknown. Here, we address this question by using phase-field modeling to produce nanoporous particles and to investigate their elastoplastic swelling behavior and fracture. Our simulations show that enhanced mechanical stability results from the network topology consisting of ligaments connected by bulbous, sphere-like nodes. The ligaments forcefully resist elongation while the nodes, behaving like isolated spherical particles, experience large stresses driving fracture. However, being smaller compared to a sphere of the same volume as the entire nanoporous particle, the nodes are more protected against fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pit growth kinetics in aluminum: effects of salt loading and relative humidity

Abstract The growth kinetics of localized corrosion, e.g. pits, in corrosive environments often controls the service life of metallic components. Yet, our understanding of these kinetics is largely based on coupon-level, e.g. mass-loss, studies which provide limited insights into the evolution of individual damage events. It is critical to relate observed cumulative loss trends, such as links between changing humidity and mass loss rates, to the growth kinetics of individual pits. Towards this goal, we leverage in-situ X-ray computed tomography to measure the growth rates of over sixty pits in aluminum in four different humid, chloride environments over ≈3 days of exposure. Pit growth rates and final volumes increased with increasing droplet volume, which was observed to increase with increasing humidity and salt loading. Two factors, droplet spreading and oxide jacking, dramatically increased pit growth rates and final volumes.

Noell, Philip J. (ORCID:0000000346998999)↗
DOE OSTI2021

In situ synchrotron diffraction and modeling of non-equilibrium solidification of a MnFeCoNiCu alloy

The solidification mechanism and segregation behavior of laser-melted Mn 35 Fe 5 Co 20 Ni 20 Cu 20 was firstly investigated via in situ synchrotron x-ray diffraction at millisecond temporal resolution. The transient composition evolution of the random solid solution during sequential solidification of dendritic and interdendritic regions complicates the analysis of synchrotron diffraction data via any single conventional tool, such as Rietveld refinement. Therefore, a novel approach combining a hard-sphere approximation model, thermodynamic simulation, thermal expansion measurement and microstructural characterization was developed to assist in a fundamental understanding of the evolution of local composition, lattice parameter, and dendrite volume fraction corresponding to the diffraction data. This methodology yields self-consistent results across different methods. Via this approach, four distinct stages were identified, including: (I) FCC dendrite solidification, (II) solidification of FCC interdendritic region, (III) solid-state interdiffusion and (IV) final cooling with marginal diffusion. It was found out that in Stage I, Cu and Mn were rejected into liquid as Mn 35 Fe 5 Co 20 Ni 20 Cu 20 solidified dendritically. During Stage II, the lattice parameter disparity between dendrite and interdendritic region escalated as Cu and Mn continued segregating into the interdendritic region. After complete solidification, during Stage III, the lattice parameter disparity gradually decreases, demonstrating a degree of composition homogenization. The volume fraction of dendrites slightly grew from 58.3 to 65.5%, based on the evolving composition profile across a dendrite/interdendritic interface in diffusion calculations. Postmortem metallography further confirmed that dendrites have a volume fraction of 64.7 ± 5.3% in the final microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

3D isotope density measurements by energy-resolved neutron imaging

Abstract Tools for three-dimensional elemental characterization are available on length scales ranging from individual atoms, using electrons as a probe, to micrometers with X-rays. However, for larger volumes up to millimeters or centimeters, quantitative measurements of elemental or isotope densities were hitherto only possible on the surface. Here, a novel quantitative elemental characterization method based on energy-resolved neutron imaging, utilizing the known neutron absorption cross sections with their ‘finger-print’ absorption resonance signatures, is demonstrated. Enabled by a pixilated time-of-flight neutron transmission detector installed at an intense short-pulsed spallation neutron source, for this demonstration 3.25 million state-of-the-art nuclear physics neutron transmission analyses were conducted to derive isotopic densities for five isotopes in 3D in a volume of 0.25 cm 3 . The tomographic reconstruction of the isotope densities provides elemental maps similar to X-ray microprobe maps for any cross section in the probed volume. The bulk isotopic density of a U-20Pu-10Zr-3Np-2Am nuclear transmutation fuel sample was measured, agrees well with mass-spectrometry and is evidence of the accuracy of the method.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multicolor multifocal 3D microscopy using in-situ optimization of a spatial light modulator

Abstract Multifocal microscopy enables high-speed three-dimensional (3D) volume imaging by using a multifocal grating in the emission path. This grating is typically designed to afford a uniform illumination of multifocal subimages for a single emission wavelength. Using the same grating for multicolor imaging results in non-uniform subimage intensities in emission wavelengths for which the grating is not designed. This has restricted multifocal microscopy applications for samples having multicolored fluorophores. In this paper, we present a multicolor multifocal microscope implementation which uses a Spatial Light Modulator (SLM) as a single multifocal grating to realize near-uniform multifocal subimage intensities across multiple wavelength emission bands. Using real-time control of an in-situ-optimized SLM implemented as a multifocal grating, we demonstrate multicolor multifocal 3D imaging over three emission bands by imaging multicolored particles as well as Escherichia coli ( E. coli ) interacting with human liver cancer cells, at $$\sim 2.5$$ ∼ 2.5 multicolor 3D volumes per second acquisition speed. Our multicolor multifocal method is adaptable across SLM hardware, emission wavelength band locations and number of emission bands, making it particularly suited for researchers investigating fast processes occurring across a volume where multiple species are involved.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Dynamic compression of Ce and Pr with millisecond time-resolved X-ray diffraction

Abstract Both cerium (Ce) and praseodymium (Pr) undergo a volume collapse transition under compression that originate from similar electronic mechanisms. Yet the outcome could not be more different. In the case of Ce with one affected 4f electron the volume collapse leaves the crystal symmetry intact, whereas for Pr with two 4f electrons the crystal symmetry changes from a distorted face centered cubic structure to a lower symmetry orthorhombic structure. In this paper, we present a study of the effect of strain/compression rate spanning nearly 4 orders of magnitude on the volume collapse phase transitions in Ce and Pr. These dynamic compression experiments in a diamond anvil cell also reveal kinetic differences between the phase transformations observed in these two materials. The transition cannot be overdriven in pressure in Ce, which indicates a fast kinetic process, whereas fast compression rates in Pr lead to a shift of the phase boundary to higher pressures, pointing to slower kinetics possibly due to the realization of a new crystal structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations

Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$N_\text {V}$$ N V ), He atoms ( $$N_\text {He}$$ N He ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$N_\text {V}$$ N V vacancies and $$N_\text {He}$$ N He helium atoms, the $$N_\text {He}/N_\text {V}$$ N He / N V ratio that causes the event, the resulting increase in $$N_\text {V}$$ N V , and the associated shift of the bubble depth are formulated as a function of $$N_\text {V}$$ N V and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$N_\text {He}/N_\text {V}$$ N He / N V ratio in order to burst. The burst depth and $$N_\text {He}/N_\text {V}$$ N He / N V are also modeled as a function of $$N_\text {V}$$ N V and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.

36 MATERIALS SCIENCE↗