Search NASA⌕ Search

SEARCH · Search NASA

Results for “helium bubble”

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 91 records · Page 5

Incorporating bubble growth volume feedback to improve simulation of the response of a structure containing liquid and gas to sudden energy input

The SNS target module is a stainless-steel structure that contains and directs mercury. The mercury is struck with short, intense proton pulses to create neutrons. The pulses deposit energy and cause loads on the target structure and mercury cavitation. Helium bubbles are introduced into the mercury to reduce loads and cavitation. This gas complicates the prediction of the target module’s physical response. The current method used to simulate the structural response to the proton pulse incorporates a simple but effective approach to approximate mercury cavitation effects. This method cannot account for the changes from injected non-condensable gas. The ability to model the target vessel response is a prerequisite for estimating its life. Known approaches for simulating a bubbly fluid mixture’s response are computationally expensive and impractical for detailed engineering models. A constitutive model is proposed for target mercury with injected gas bubbles, with the intent of providing more accurate mercury–vessel response behavior without the computing costs required to track individual bubbles. The model’s assumptions and a computational model embodiment are described. The bubbly mercury computational model is subjected to test problems to assess its potential utility. Here, the model was able to better predict the effect of gas on wave propagation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupled Thermal-Hydraulic Analysis and Species Mass Transport in a Versatile Experimental Salt Irradiation Loop (VESIL) Using CTF

With the resurgence of interest in molten salt reactors, there is a need for new experiments and modeling capabilities to characterize the unique phenomena present in this fluid fuel system. A Versatile Experimental Salt Irradiation Loop (VESIL) is currently under investigation at Idaho National Laboratory to be placed in the Advanced Test Reactor (ATR). One of the key phenomena this proposed experiment plans to elucidate is fission product speciation in the fuel-salt and the subsequent effects this has on the fuel-salt properties, source term generation, and corrosion control. Specifically, noble gases (Xe & Kr) will bubble out to a plenum or off-gas system, and noble metals (Mo, Tc, Te, etc.) will precipitate and deposit in specific zones in the loop. This work extends the mass transfer and species interaction models in CTF (Coolant-Boiling in Rod Arrays—Two Fluids) and applies these models to give a preliminary estimation of fission product behavior in the proposed VESIL design. A noble metal–helium bubble mass transfer model is coupled with the thermal-hydraulic results from CTF to determine the effectiveness of this insoluble fission product (IFP) extraction method for VESIL. Amounts of IFP species extracted to the off-gas system and species distributions in VESIL after a 60-day ATR cycle are reported.

Walker, Samuel A.↗

Effects of helium on irradiation response of reduced-activation ferritic-martensitic steels: Using nickel isotopes to simulate fusion neutron response

Understanding the effects of helium on microstructures and mechanical properties of reduced-activation ferritic-martensitic steels is important to use of these steels in fusion reactor structures. The 9Cr-2WVTa steels were doped with 58Ni and 60Ni isotopes at 2 weight percent to control the rate of transmutation helium generation. The samples were irradiated in the High Flux Isotope Reactor to ~24 displacements per atom at nominal temperatures of 300, 400, and 500°C, producing 228 and 7 atomic parts-per-million helium in the 58 Ni- and 60 Ni-doped samples, respectively. Transmission electron microscopy revealed a variety of precipitates and the radiation-induced dislocation loops and cavities (voids or helium bubbles). Tensile tests of the irradiated samples at the irradiation temperatures showed radiation-induced hardening at 300°C and radiation-induced softening at 400°C. Analysis indicates that the hardening primarily originated from the loops and cavities. The 58 Ni-doped samples had greater strengthening contributions from loops and cavities, leading to higher hardening with lower ductility than the 60 Ni-doped samples. The greater helium production of 58 Ni did not show pronounced reductions in ductility of the samples.

36 MATERIALS SCIENCE↗

Non-dilute helium-related defect interactions in the near-surface region of plasma-exposed tungsten

We report a systematic energetic analysis of helium-related defect interactions that mediate helium (He) segregation on surfaces of plasma-exposed tungsten at different levels of He ion implantation. We focus on high He fluences that increase the He content in the plasma-exposed material well beyond the dilute limit of He concentration and employ atomic configurations generated by large-scale molecular dynamics simulations of He-implanted tungsten. We perform systematic molecular statics computations of cluster– defect interaction energetics in the highly defect-rich near-surface region of plasma-exposed tungsten for small mobile helium clusters as a function of the cluster distance from the surface. In this region, mobile clusters are also subjected to the stress fields generated by defects such as helium bubbles and other clusters, which govern cluster–defect interactions in addition to the cluster–surface interaction. Based on systematic investigation of individual cluster–defect interactions, we develop a mathematical framework to describe the interaction energy landscapes consisting of elastic interaction potential perturbations to finite-width square-well potentials, where the potential well accounts for cluster trapping by the defect at close range and subsequent coalescence and the perturbation potential is parameterized according to elastic inclusion theory. Superposition of all the relevant interaction potentials provides a comprehensive description of the interaction energy landscape that would be experienced by a small mobile cluster along its migration path toward the plasma-exposed surface at high He fluence. Such descriptions are particularly important for developing atomistically-informed, hierarchical multi-scale models of helium cluster dynamics in plasma-facing materials.

36 MATERIALS SCIENCE↗

Microstructural characterization of cold-worked 316 stainless steel flux thimble tubes irradiated up to 100 dpa in a commercial Pressurized Water Reactor

Two flux thimble tubes (FTT) made of 15% cold-worked 316 stainless steel (SS) were harvested from Ringhals Pressurized Water Reactor (PWR) Unit 2, with peak damages of 76 and 100 displacements per atom (dpa) after 29 and 34 years’ service, respectively. Specimens sectioned from parent tubes were comprehensively characterized with nominal damage levels of ~0, ~41, ~74, 76, and 100 dpa at a nominal temperature range of 285-323 °C. Both FTTs contained helium and hydrogen gases as transmutation products. The helium follows a production rate of ~9.8 appm/dpa, while environmental factors complicate hydrogen production obscuring an exact H/dpa ratio. Irradiation-induced dislocation loops, nano-cavities, solute clusters, and microsegregation were all observed. The dislocation loops and nano-cavities indicated saturation at 41 dpa. The solute clusters continued to evolve with Ni-Si clusters formed at 41 dpa, and Ni-Si-Mn-P clusters formed at 74 and 100 dpa, but neither clusters exhibited distinct diffraction patterns at any damage levels. Solute clusters were observed to frequently be co-located with dislocation loops, but fully decorated loops were rarely detected. Significant radiation-induced segregation (RIS) was observed around grain boundaries at all damage levels. The modified inverse Kirkendall (MIK) model captured the RIS behavior of major elements. Large cavities within or around an Mn-S rich region were observed for the first time. Through all the damage levels, void swelling is always below 0.05%, making significant dimensional change unlikely in core internals when used at similar conditions. Meanwhile, the role of overwhelming nanocavities, presumably helium bubbles, should be considered in other potential degradation mechanisms, including irradiation-assisted stress corrosion cracking, embrittlement, and loss of fracture toughness, which remain the concerns for extended operation of nuclear power plants.

Post Irradiation Examination, Isotope Dilution Mas↗

Segregation behavior and phase instability of Eurofer97 after neutron irradiation to 72 dpa

Here, a reduced activation ferritic/martensitic steel, Eurofer97, was neutron irradiated in the vicinity of 300 °C in the High Flux Isotope Reactor (HFIR) up to 72 dpa. Advanced analytical scanning transmission electron microscopy and conventional transmission electron microscopy were applied to investigate the radiation-induced segregation and phase instability behavior after neutron irradiation. Amorphization was observed in M 23 C 6 carbides. Cr-rich clusters were seen within the matrix, near the lath boundaries and close to the M 23 C 6 carbides. Cr enrichment and Fe depletion were detected at both prior austenite grain boundaries and lath boundaries, despite different segregation magnitude. In addition, the enrichment of Ni, the depletion of V, and tiny cavities (presumably helium bubbles) are also found at lath boundaries. This work interrogates the evolution of microstructures after neutron irradiation, which provides detailed understanding on the microstructural aspects controlling the mechanical integrity of Eurofer97 under high-dose neutron damage.

36 MATERIALS SCIENCE↗

No ball milling needed: Alternative ODS steel manufacturing with gas atomization reaction synthesis (GARS) and friction-based processing

Oxide dispersion strengthened (ODS) steels are promising structural materials for future fusion reactors. The high-density (~10 23 /m 3 ) of highly stable Y-(Ti)-O nano-oxides provide high sink strength for radiation resistance and high-temperature (> 650 °C) creep strength. Concomitantly, helium management is enabled by trapping high density (~10 23 /m 3 ) of small (< 3 nm) helium bubbles in the vicinity of nano-oxides. However, conventional route of making ODS steels involves prolonged ball milling, canning, degassing, and laborious thermo-mechanical processing (TMP). Such route, especially the batch-by-batch ball milling step, faces persistent challenge with scalability and high costs. Gas atomization reaction synthesis (GARS) method has demonstrated the potential of making precursor ODS steel powders without ball milling, but the nano-oxide density was around 10 21 /m 3 in the final consolidated form by conventional TMP. Taking advantage of GARS precursor powder, we use friction-based processing, including friction consolidation and extrusion, to manufacture ODS steel with further improved nano-oxide characteristics. Here, preliminary results showed that Y/Ti/O species were intimately mixed and rapidly reacted to form nano-oxides with a number density of ~10 22 /m 3 .

36 MATERIALS SCIENCE↗

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗

Deciphering the influence of high-flux helium plasma exposure on tungsten nanostructure growth using real-time, in situ spectroscopic ellipsometry

In situ analysis of surfaces during high-flux plasma exposure represents a long-standing challenge in the study of plasma–material interactions. While post-mortem microscopy can provide a detailed picture of structural and compositional changes, in situ techniques can capture the dynamic evolution of the surface. In this study, we demonstrate how spectroscopic ellipsometry can be applied to the real-time characterization of W nanostructure (also known as “fuzz”) growth during exposure to low temperature, high-flux He plasmas. Strikingly, over a wide range of sample temperatures and helium fluences, the measured ellipsometric parameters (Ψ,Δ) collapse onto a single curve that can be directly correlated with surface morphologies characterized by ex situ helium ion microscopy. The initial variation in the (Ψ,Δ) parameters appears to be governed by small changes in surface roughness (<50 nm) produced by helium bubble nucleation and growth, followed by the emergence of 50 nm diameter W tendrils. This basic behavior appears to be reproducible over a wide parameter space, indicating that the spectroscopic ellipsometry may be of general practical use as a diagnostic to study surface morphologies produced by high-flux He implantation in refractory metals. An advantage of the methods outlined here is that they are applicable at low incident ion energies, even below the sputtering threshold. As an example of this application, we apply in situ ellipsometry to examine how W fuzz growth is affected both by varying ion energy and the temperature of the surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Compact Gas Liquid Separator for the Spallation Neutron Source Mercury Process Loop

Upgrades at the spallation neutron source (SNS) accelerator at Oak Ridge National Laboratory are underway to double its proton beam power from 1.4 to 2.8 MW. About 2 MW will go to the current first station while the rest will go to the future Second Target Station. The increase of beam power to the first target station is especially challenging for its mercury target. When the short proton beam hits the target, strong pressure waves are generated, causing cavitation erosion and challenging stresses for the target's weld regions. SNS has successfully operated reliably at 1.4 MW by mitigating the pressure wave with the injection of small Helium bubbles into the mercury. To operate reliably at 2 MW, more gas will be injected into mercury to mitigate the pressure wave further. However, the mercury process loop was not originally designed for gas injection, and the accumulation of gas in the pipes is a concern. Due to space constraints, a custom gas liquid separator (GLS) was designed to fit a 90-deg horizontal elbow space in the SNS mercury loop. Simulations and experiments were performed, and a successful design was developed that has the desired efficiency while keeping the pressure losses acceptable.

42 ENGINEERING↗

Kinetic Monte Carlo simulations of aging in δ -Pu

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

Actinides↗

Comprehensive Characterization of Helium-Induced Degradation of the Friction Stir Weld on Neutron-Irradiated 304L Stainless Steel

The report describes new experimental results on the mechanical performance of the friction stir welds made on neutron-irradiated 304L stainless steel with helium. The report focuses on helium-related issues (e.g., the helium-induced degradation in the welded joint), aiming at the need to repair irradiated components of nuclear power plants. The friction stir welds analyzed here were previously produced at the US Department of Energy’s Oak Ridge National Laboratory, and initial characterization work was performed, mostly addressing the microstructure and if macroscopic cracks and helium bubbles present. The present work attempts to perform a more comprehensive study to assess mechanical performance (i.e., microhardness distribution, tensile properties, tensile deformation behaviors, and fractography analysis).

36 MATERIALS SCIENCE↗

2024 Second Half Semi Annual Report: Modeling plasticity-mediated flow in metals with pressurized cavities

The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.

36 MATERIALS SCIENCE↗

Transport of Highly Volatile Gases Related to Noble Gases and Tritium in the MSRE

This study aims to comprehend mass transfer in the closed-loop circulation of highly volatile gases, including noble gases and tritium. We explore the impact of steady-state xenon-135 and tritium on the MSRE and reveal their isotopic distributions using online noble gas stripping of fuel salts. The MSRE was engineered to extract fission product gases from fuel salts and efficiently eliminate inert gases with the help of helium bubbles within a circulating fuel pump. These reactors introduce significant theoretical challenges in estimating interfacial area and mass transfer coefficients, crucial for modeling mass transfer processes. An essential component of our analysis is the mass transfer coefficient. These coefficients are important for understanding how radionuclides move during various phase transitions within a nuclear reactor. Xenon-135 and tritium are found in both liquid and gas phases within the reactor system. In the liquid phase, they dissolve in molten salts, while in the gas phase, they manifest as bubbles. These elements have significant adverse effects on reactor operation due to their strong neutron absorption properties, influencing both safety and performance. The Mole code, which predicts the behavior of chemical species under steady-state conditions, facilitates multiphysics coupling with Griffin to update species distributions and address inherent MSR safety.

Lee, Kyoung↗

Corrigendum and Addendum: Helium flux effects on bubble growth and surface morphology in plasma-facing tungsten from large-scale molecular dynamics simulations (2019 Nucl. Fusion 59 066035)

Two of the simulations discussed in a prior article (Hammond et al 2019 Nucl. Fusion 59 066035) were affected by a simulation glitch. We repeated the affected calculations and discuss them here. Here, the overall conclusions are essentially unchanged, though the details are different. In particular, observations that we referred to as ‘concerted bursting’ were caused primarily by non-physical heating and cooling applied by the thermostat after most atoms’ velocities were deleted (for reasons that are not known for certain). The phenomenon of one bubble bursting and causing another nearby bubble to burst does exist, though its effects are much less spectacular in the absence of non-physical driving forces. The observation of an interconnected network of sub-surface cavities formed by burst bubbles is real, and the observation of holes on the surface 1–2 nm in diameter is also confirmed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Microstructural Effects of High Dose Helium Implantation in ErD 2

Metal hydrides can store hydrogen isotopes with high volumetric density. In metal tritides, tritium beta decay can result in accumulation of helium within the solid, in some cases exceeding 10 at.% helium after only 4 years of aging. Helium is insoluble in most materials, but often does not readily escape, and instead coalesces to form nanoscale bubbles when helium concentrations are near 1 at.%. Blistering or spallation often occurs at higher concentrations. Radioactive particles shed during this process present a potential safety hazard. This study investigates the effects of high helium concentrations on erbium deuteride (ErD 2 ), a non-radioactive surrogate material for erbium tritide (ErT 2 ). To simulate tritium decay in the surrogate, high doses of 120 keV helium ions were implanted into ErD 2 films at room temperature. Scanning and transmission electron microscopy indicated spherical helium bubble formation at a critical concentration of 1.5 at.% and bubble linkage leading to nanoscale crack formation at a concentration of 7.5 at.%. Additionally, crack propagation occurred through the nanocrack region, resulting in spallation extending from the implantation peak to the surface. Electron energy loss spectroscopy was utilized to confirm the presence of high-pressure helium in the nanocracks, suggesting that helium gas plays a predominant role in deformation. This work improves the overall understanding of helium behavior in ErD 2 by using modern characterization techniques to determine: the critical helium concentration required for bubble formation, the material failure mechanism at high concentration, and the nanoscale mechanisms responsible for material failure in helium implanted ErD 2 .

36 MATERIALS SCIENCE↗

Experimental observations of shock-wave-induced bubble collapse and hot-spot formation in nitromethane liquid explosive

We studied the collapse of individual helium gas bubbles in the homogeneous explosive nitromethane (NM) to investigate effects of hot-spot formation on the detonation process. A bubble was injected into a NM sample, and a shock wave from an explosive detonator compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and the pre- and post-bubble collapse processes. An infrared camera image showed the residual radiance temperature distribution after the bubble collapse, and an optical fiber pyrometer measured the time-resolved thermal radiance. We measured the optical spectra of light emitted from detonating NM without a bubble and from a collapsing bubble in shocked, undetonated NM. We estimated temperatures of the detonation fronts and of the hot spots formed by bubble collapse. To study the incipient detonation process, we performed all bubble collapse experiments at pressures below the threshold for creating a sustained detonation. Where the bubble collapsed, we observed an opaque, thermally emissive region believed to be chemical reaction products. Chemical reactions in NM can be produced with lower shock pressures (~1 GPa) when a helium bubble is present than without a bubble (~10 GPa). We used hydrodynamic modeling to predict shock wave propagation, extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. Simulations using a temperature-dependent Arrhenius burn model gave much better results than reactive burn models that depend only on pressure and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗