Search NASA⌕ Search

SEARCH · Search NASA

Results for “He 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 73 records · Page 4

The Investigation of the Effects of Gravity on Single Bubble Sonoluminescence

In single bubble following it's rapid collapse each cycle of oscillation of an ultrasonic field. Since widely varying length and time scales affect the bubble dynamics and optical emission processes, it is difficult to anticipate the importance of the effects of gravity present for observations on earth. Our bubble is driven in an acoustically resonating cavity at it's first harmonic mode. The acoustical radiation pressure (Bjerknes force) will then keep it suspended in the center near the pressure antinode. When driven in a region where the diffusive processes balance the bubble it acts in a nonlinear but regular way, emitting a short (approx. 200ps) burst of light each acoustic cycle. Balancing the Bjerknes force with buoyancy, as in, we can see that the bubble should be displaced from the velocity node approximately 20m at normal gravity. Therefore, water flows past the bubble at the time of collapse. Gravitation also changes the ambient pressure at the bubble's location, as Delta.P = rho.g.h this gives a change of approximately -0.5% in our experiment when going from 1.8g to 0g. Studies of ambient pressure changes were also done in order to assess these effects. Inside a pressure sealed chamber a spherical glass cell is filled with distilled water which has been degassed to 120mmHg. A bubble is then trapped in the center and driven by a piezoelectric transducer at 32.2kHz attached to the side of the cell. An optical system is then set up to take strobbed video images along and light emission data simultaneously. Temperature, pressure, drive voltage, and listener voltage are also monitored. PMT output in Volts The radii of the bubbles for both experiment s are fit using the Rayleigh-Plesset equation and the acoustic drive amplitude and the ambient bubble radius are found. There is little change in the acoustic drive amplitude as we expect, since we are not varying the drive voltage. However. the ambient bubble radius goes up considerably. These changes (increased light output, increased maximum bubble radius, and increased ambient bubble radius) are also observed when the ambient pressure is varied in the laboratory by an amount similar to that due to gravitation. The changes in the ambient bubble radius and light output with a change in ambient pressure are predicted by the "dissociation hypothesis" and have been observed by other groups in the laboratory. It seems clear that buoyancy's effect on light output and bubble radius, are at best on the same order as the effects of ambient pressure.

Dzikowicz, Ben↗

The Isolated Bubble Regime in Pool Nucleate Boiling

We consider an isolated bubble boiling regime in which vapour bubbles are intermittently produced at a prearranged set of nucleation site on an upward facing overheated wall plane. In this boiling regime, the bubbles depart from the wall and move as separate entities. Except in the matter of rise velocity, the bubbles do not interfere and are independent of one another. However, the rise velocity is dependent on bubble volume concentration in the bulk. Heat transfer properties specific to this regime cannot be described without bubble detachment size, and we apply our previously developed dynamic theory of vapour bubble growth and detachment to determine this size. Bubble growth is presumed to be thermally controlled. Two limiting cases of bubble evolution are considered: the one in which buoyancy prevails in promoting bubble detachment and the one in which surface tension prevails. We prove termination of the isolated regime of pool nucleate boiling to result from one of the four possible causes, depending on relevant parameters values. The first cause consists in the fact that the upward flow of rising bubbles hampers the downward liquid flow, and under certain conditions, prevents the liquid from coming to the wall in an amount that would be sufficient to compensate for vapour removal from the wall. The second cause is due to the lateral coalescence of growing bubbles that are attached to their corresponding nucleation sites, with ensuing generation of larger bubbles and extended vapour patches near the wall. The other two causes involve longitudinal coalescence either 1) immediately in the wall vicinity, accompanied by the establishment of the multiple bubble boiling regime, or 2) in the bulk, with the formation of vapour columns. The longitudinal coalescence in the bulk is shown to be the most important cause. The critical wall temperature and the heat flux density associated with isolated bubble regime termination are found to be functions of the physical and operating parameters and are discussed in detail.

Buyevich, Y. A.↗

The role of irradiation-enhanced interstitial diffusion in over-pressurizing fission gas bubbles in UO 2

Fission gas bubbles in UO 2 nuclear fuel have been observed to exhibit pressures in excess of the equilibrium bubble pressure; however, the cause of bubble over-pressurization has not yet been demonstrated. The mechanical interaction between a bubble and the surrounding matrix or grain boundary depends on the internal pressure of the bubble and local stress state, such that over-pressurized bubbles are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp. Here, in this work, we investigate the role of U interstitials, produced through irradiation, in over-pressurizing bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the capture of interstitials and vacancies by bubbles have been determined from MD as a function of the ratio of gas atoms to vacancies that make up the bubble. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to predict bubble over-pressurization as a function of temperature for typical fission rates. It was found that there is a transition from low pressure bubbles (at high temperatures) to high pressure bubbles (at lower temperatures). The cause of this behavior was shown to be the creation of irradiation-induced interstitials that are highly mobile relative to vacancies at low temperature; whereas, vacancies are sufficiently mobile at high temperatures to limit bubble pressures. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and that this behavior is highly sensitive to the local pellet temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Guest investigator program study: Physics of equatorial plasma bubbles

Plasma bubbles are large-scale (10 to 100 km) depletions in plasma density found in the night-time equatorial ionosphere. Their formation has been found to entail the upward transport of plasma over hundreds of kilometers in altitude, suggesting that bubbles play significant roles in the physics of many of the diverse and unique features found in the low-latitude ionosphere. In the simplest scenario, plasma bubbles appear first as perturbations in the bottomside F layer, which is linearly unstable to the gravitationally driven Rayleigh-Taylor instability. Once initiated, bubbles develop upward through the peak of the F layer into its topside (sometimes to altitudes in excess of 1000 km), a behavior predicted by the nonlinear form of the same instability. While good general agreement has been found between theory and observations, little is known about the detailed physics associated with plasma bubbles. Our research activity centered around two topics: the shape of plasma bubbles and associated electric fields, and the day-to-day variability in the occurrence of plasma bubbles. The first topic was pursued because of a divergence in view regarding the nonlinear physics associated with plasma bubble development. While the development of perturbations in isodensity contours in the bottomside F layer into plasma bubbles is well accepted, some believed bubbles to be cylinder-like closed regions of depleted plasma density that floated upward leaving a turbulent wake behind them (e.g., Woodman and LaHoz, 1976; Ott, 1978; Kelley and Ott, 1978). Our results, summarized in a paper submitted to the Journal of Geophysical Research, consisted of incoherent scatter radar measurements that showed unambiguously that the depleted region is wedgelike and not cylinderlike, and a case study and modeling of SM-D electric field instrument (EFI) measurements that showed that the absence of electric-field perturbations outside the plasma-depleted region is a distinct signature of wedge-shaped plasma bubbles. The second topic was pursued because the inability to predict the day-to-day occurrence of plasma bubbles indicated inadequate knowledge of the physics of plasma bubbles. An understanding of bubble formation requires an understanding of the roles of the various terms in the linearized growth rate of the collisional Rayleigh-Taylor instability. In our study, we examined electric-field perturbations found in SM-D EFI data and found that the seeding is more likely to be produced in the E region rather than the F region. The results of this investigation are presented in the Appendix of this report and will be submitted for publication in the Journal of Geophysical Research.

Tsunoda, Roland T.↗

Generation of Bubbly Suspensions in Low Gravity

Generation of a uniform monodisperse bubbly suspension in low gravity is a rather difficult task because bubbles do not detach as easily as on Earth. Under microgravity, the buoyancy force is not present to detach the bubbles as they are formed from the nozzles. One way to detach the bubbles is to establish a detaching force that helps their detachment from the orifice. The drag force, established by flowing a liquid in a cross or co-flow configuration with respect to the nozzle direction, provides this additional force and helps detach the bubbles as they are being formed. This paper is concerned with studying the generation of a bubbly suspension in low gravity in support of a flight definition experiment titled "Behavior of Rapidly Sheared Bubbly Suspension." Generation of a bubbly suspension, composed of 2 and 3 mm diameter bubbles with a standard deviation <10% of the bubble diameter, was identified as one of the most important engineering/science issues associated with the flight definition experiment. This paper summarizes the low gravity experiments that were conducted to explore various ways of making the suspension. Two approaches were investigated. The first was to generate the suspension via a chemical reaction between the continuous and dispersed phases using effervescent material, whereas the second considered the direct injection of air into the continuous phase. The results showed that the reaction method did not produce the desired bubble size distribution compared to the direct injection of bubbles. However, direct injection of air into the continuous phase (aqueous salt solution) resulted in uniform bubble-diameter distribution with acceptable bubble-diameter standard deviation.

Nahra, Henry K.↗

Modeling of Vapor Bubble Growth Under Nucleate Boiling Conditions in Reduced Gravity

A dynamic model is developed to describe the evolution of a vapor bubble growing at a nucleation site on a superheated surface under arbitrary gravity. The bubble is separated from the surface by a thin microlayer and grows due to the evaporation from the microlayer interface. The average thickness of the microlayer increases as the bubble expands along the surface if the evaporation rate is lower than some critical value. The corresponding threshold value of the surface temperature has to be associated with the burn-out crisis. Two main reasons make for bubble separation, which are the buoyancy force and a force caused by the vapor momentum that comes to the bubble with vapor molecules. The latter force is somewhat diminished if condensation takes place at the upper bubble surface in subcooled liquids. The action of the said forces is opposed by inertia of the additional mass of liquid as the bubble center rises above the surface and by inertia of liquid being expelled by the growing bubble in radial directions. An extra pressure force arises due to the liquid inflow into the microlayer with a finite velocity. The last force helps in holding the bubble close to the surface during an initial stage of bubble evolution. Two limiting regimes with distinctly different properties can be singled out, depending on which of the forces that favor bubble detachment dominates. Under conditions of moderately reduced gravity, the situation is much the same as in normal gravity, although the bubble detachment volume increases as gravity diminishes. In microgravity, the buoyancy force is negligible. Then the bubble is capable of staying near the surface for a long time, with intensive evaporation from the microlayer. It suggests a drastic change in the physical mechanism of heat removal as gravity falls below a certain sufficiently low level. Inferences of the model and conclusions pertaining to effects caused on heat transfer processes by changes in bubble hydrodynamics induced by gravity are discussed in connection with experimental evidence, both available in current and in as yet unpublished literature.

Buyevich, Yu A.↗

Pressure fields produced by single-bubble collapse near a corner

Damage produced by repeated bubble collapse to neighboring rigid objects in hydraulic systems is an important consequence of cavitation. Although bubble collapse near a single wall has received significant attention in the past, few studies exist on the dynamics of bubbles collapsing near a corner, i.e., two flat rigid surfaces intersecting at a right angle. In this work we quantify the pressure fields produced by a single bubble collapsing near two perpendicular rigid walls. Using a high-order accurate shock- and interface-capturing method to solve the three-dimensional compressible Navier-Stokes equations for gas and liquid flows, we simulate the dynamics of a single bubble collapsing at different initial stand-off distances from the two walls. In contrast to a bubble collapsing near a single wall, the collapse of bubbles within a critical stand-off distance is not symmetric about the bisecting plane due to the interaction between the bubble and the second wall. The second wall affects the pressure produced during the collapse in the following ways: (i) For bubbles initially located sufficiently close to both walls, the reentrant jet produced during collapse no longer points in the direction normal to the closest wall but at an angle toward the corner, (ii) the part of the emitted shock with the highest amplitude propagates in line with the jet, and (iii) the bubble migrates in that same direction during its collapse with a dependence on the stand-off distance, consistent with predictions made using Kelvin impulse. The location of maximum pressure along the walls is measured for the different initial stand-off distances. Using acoustic arguments, we find a semiempirical relationship to predict the initial stand-off distances for which the maximum pressure occurs in the corner. Here we find that when the bubble is sufficiently close to equidistant from each boundary, the maximum pressure is observed in the corner due to the water-hammer and implosion shocks reflecting off the boundaries and intersecting in the corner. We also show that when bubbles are initially attached to either wall the wall pressure produced can be significantly increased compared to bubbles detached from either wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A phase field study of the thermal migration of gas bubbles in UO 2 nuclear fuel under temperature gradient

Phase field models are developed to study the gas bubble migration in uranium dioxide nuclear fuel in which a large temperature gradient exists during the operation. In this work, thermal diffusion mechanism for nanosized gas bubbles and vapor transport process for micron-sized gas bubbles are considered, respectively. In both cases, gas bubbles migrate to the high-temperature area. Due to the velocity difference between leading and trailing edges of the gas bubbles, nanosized gas bubbles are elongated along the temperature gradient direction when thermal diffusion is dominated. Micron-sized gas bubbles are either compressed along temperature gradient direction to form lenticular shape bubbles or elongated along temperature gradient direction, depending on the location of the gas bubbles within the fuel pellet. Initial gas bubble radius has no significant effect on the gas bubble migration velocity for both thermal diffusion and vapor transport mechanisms. We notice that the shape change of the gas bubble due to vapor transport mechanism has no significant effect on the migration velocity. Furthermore, the center cavity formation is also captured by our model which is due to the migration and accumulation of lenticular gas bubbles at the center of the fuel pellet. The modeling results compare well with experimental observations and theoretical analysis in the literature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗

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↗

Measurement of the Shear Lift Force on a Bubble in a Channel Flow

Two-phase flow systems play vital roles in the design of some current and anticipated space applications of two-phase systems which include: thermal management systems, transfer line flow in cryogenic storage, space nuclear power facilities, design and operation of thermal bus, life support systems, propulsion systems, In Situ Resource Utilization (ISRU), and space processes for pharmaceutical applications. The design of two-phase flow systems for space applications requires a clear knowledge of the behaviors of the dispersed phase (bubble), its interaction with the continuous phase (liquid) and its effect on heat and mass transfer processes, The need to understand the bubble generation process arises from the fact that for all space applications, the size and distribution of bubbles are extremely crucial for heat and mass transfer control. One important force in two-phase flow systems is the lift force on a bubble or particle in a liquid shear flow. The shear lift is usually overwhelmed by buoyancy in normal gravity, but it becomes an important force in reduced gravity. Since the liquid flow is usually sheared because of the confining wall, the trajectories of bubbles and particles injected into the liquid flow are affected by the shear lift in reduced gravity. A series of experiments are performed to investigate the lift force on a bubble in a liquid shear flow and its effect on the detachment of a bubble from a wall under low gravity conditions. Experiments are executed in a Poiseuille flow in a channel. An air-water system is used in these experiments that are performed in the 2.2 second drop tower. A bubble is injected into the shear flow from a small injector and the shear lift is measured while the bubble is held stationary relative to the fluid. The trajectory of the bubble prior, during and after its detachment from the injector is investigated. The measured shear lift force is calculated from the trajectory of the bubble at the detachment point. These values for the shear lift are then compared with the theoretical predictions from various published works on shear lift in the open literature, which include asymptotic solutions at low bubble Reynolds number, potential flow predictions and numerical studies that deal with intermediate bubble Reynolds numbers.

Nahra, Henry K.↗

The effect of stress on the migration of He gas bubbles under a thermal gradient in Fe by phase-field modeling

Here a phase-field model is parameterized to study the effect of elastic stresses on the migration of He gas bubbles in Fe under a temperature gradient. Stresses caused by the gas bubble pressure and residual stress in the Fe matrix are considered. The dependence of He bubble migration velocity on the magnitude of the residual stress, average temperature, temperature gradient, and bubble size is measured. In agreement with a theoretical model based on surface diffusion, simulation results demonstrate that He bubbles move towards the high temperature region with velocities in Fe that are orders of magnitude faster than previously reported in UO 2 . It is found that local stresses in the matrix caused by the He bubble have negligible effect on the bubble migration process; however, residual stresses in the Fe matrix, potentially caused by processing or irradiation, can modestly modify bubble kinetics through pressure dependence of the He diffusion coefficients. Compressive residual stress decreases diffusion coefficients for bulk and surface diffusion mechanisms, thus reducing the migration velocity of the gas bubble. In contrast, tensile residual stress increases the diffusion coefficients, resulting in an increase in the gas bubble migration velocity. This pressure dependence is also consistent with a theoretical model. This phase field model lays the foundation for analysis of bubble coalescence-induced fracture in He bubble-containing steels.

36 MATERIALS SCIENCE↗

Atomistic simulations of He bubbles in Beryllium

Formation of He bubbles can have a significant effect on the microstructural evolution and properties of irradiated materials. In this report, we use atomistic simulations based on machine learning potentials to investigate the fundamental behavior of He bubbles in Be, with a specific focus on the shape, stability, and diffusivity of bubbles. Stability of He bubbles is quantified in terms of formation energies, which are determined as a function of the ratio of He/V. We find that He bubbles become unstable with respect to plastic deformation through punch-out dislocations around the bubble when the He/V ratio is larger than ~1.25, and the punch-out process induces the change of the regular bubble shape. In general, the bubble shape of He in Be is found to be ellipsoid-like. It is also found that for a fixed He/V ratio, the bubble attracts vacancies to become larger in size. If the bubble size is constant, the bubble attracts additional He atoms until the punch-out reaction occurs. The dominant diffusion mechanism of He bubbles changes from surface diffusion to volume diffusion as the temperature is increased, with a crossover occurring at about 900 K.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Inertial collapse of a gas bubble in a shear flow near a rigid wall

Despite the extensive research on bubble collapse near rigid walls, the bubble collapse dynamics in the presence of shear flow near a rigid wall is poorly understood. We conduct direct simulations of the Navier–Stokes equations to explore the bubble dynamics and pressures during bubble collapse near a rigid, flat wall under linear shear flow conditions. We examine the dependence of the bubble collapse morphology and wall pressures on the initial bubble location and shear rate. We find that shear distorts the bubble, generating two re-entrant jets – one developing from the side opposite to the mean flow and the other from the far end toward the wall. Upon impact of the jet on the opposite side of the bubble, water-hammer shocks are produced, which propagate outward and interact with the convoluted bubble shape. The shock stretches the bubble towards the wall, resulting in a closer impact location for the jet originating from the far end compared with the case with no shear flow. The water-hammer pressure location can be approximated as the theoretical distance travelled by a particle initialised at the bubble centre with the corresponding constant shear flow velocity. The maximum wall pressures can thus be predicted by considering the distance between the far jet impingement location and the wall along the wall-normal direction. As the shear rate is increased, the maximum wall pressure increases, although only marginally. We determine the critical initial stand-off distance from the wall at which the bubble morphology is shear dominated, i.e. characterised by converging re-entrant jets.

Bhola, Sahil (ORCID:0000000236766302)↗

Sliding dynamics for bubble phases on periodic modulated substrates

We analyze a bubble-forming system composed of particles with competing long-range repulsive and short-range attractive interactions driven over a quasi-one-dimensional periodic substrate. We find various pinned and sliding phases as a function of substrate strength and drive amplitude. When the substrate is weak, a pinned bubble phase appears that depins elastically into a sliding bubble lattice. For stronger substrates, we find anisotropic bubbles, disordered bubbles, and stripe phases. Plastic depinning occurs via the hopping of individual particles from one bubble to the next in a pinned bubble lattice, and as the drive increases, there is a transition to a state where all of the bubbles are moving but are continuously shedding and absorbing individual particles. This is followed at high drives by a moving bubble lattice in which the particles can no longer escape their individual bubbles. The transition between the plastic and elastic sliding phases can be detected via signatures in the velocity-force curves, differential conductivity, and noise. When the bubbles shrink due to an increase in the attractive interaction term, they fit better inside the pinning troughs and become more strongly pinned, leading to a reentrant pinning phase. For weaker attractive terms, the size of the bubbles becomes greater than the width of the pinning troughs and the depinning becomes elastic with a reduced depinning threshold. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bubble Transport through a Porous Lattice with an Applied Inlet Flow

Within gas-evolving electrochemical systems, bubbles negatively impact performance by covering electrode active sites for reactions, blocking electric field lines, and obstructing liquid electrolyte flow causing pressure buildup. Recent additive manufacturing advances have enabled tuned porous electrode microstructures to be created, but producing systems that maximize electrochemical throughput and minimize bubble impact remains challenging. Thus, improved physical understanding of and modeling capabilities for bubble behavior are critical to improve electrolyzer design. To address this need, this study examines rising stage bubbles within a lattice with an applied liquid flow—an underexplored regime that strongly influences an electrochemical bubble’s fate. Notably, theoretical predictions and resolved bubble simulations are complemented by experiments from a 3D-printed visualization cell that matches the simulation geometry. The minimum threshold flow rate to achieve bubble breakthrough is found to be larger for higher porosities and for smaller bubbles. Different-sized bubbles decrease expected electrochemical performance in different ways; smaller bubbles tend to stay stuck but cover less solid surface, while larger bubbles more readily break through but cover more surface while in the lattice. The bubble trajectory, deformation, and contact area provide insight into these different behaviors. These findings provide design guidelines toward creating more effective electrolyzers.

Guo, Jack [Lawrence Livermore National Laboratory ↗

Investigating bubble pressures in irradiated UO 2 for understanding fragmentation of high burnup structure

This report summarized work carried out to investigate mechanisms relevant to fragmentation of high burnup UO 2 fuel. High pressure bubbles that form in high burnup structure are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp during a transient. The over-pressurization of the bubbles during normal operation is a pre-requisite for this mechanism to occur. Therefore, in this work we investigate the role of interstitials, produced through irradiation, in over-pressurizing the bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the annihilation of interstitials and vacancies at bubbles have been determined from MD as a function of bubble pressure. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to determine the steady-state defect concentrations and bubble pressure under irradiation. It was found that there is a transition from low pressure bubbles, at high temperature, to high pressure bubbles, at low temperatures. This indicates the formation of over-pressurized bubbles in the low temperature periphery of the pellet, where high burnup structure forms. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and can then contribute to fragmentation during a transient. Future work will examine the influence of bubble pressure on Xe diffusion and study the impact of other sinks (e.g. grain boundaries and dislocations) that are relevant to high burnup structure. Ultimately, the goal is to provide predictions of the bubble pressure as a function of burnup that can be used to derive a more mechanistic fragmentation threshold for use in BISON.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Sliding Dynamics for Bubble Phases on Periodic Modulated Substrates

We analyze a bubble-forming system composed of particles with competing long-range repulsive and short range attractive interactions driven over a quasi-one-dimensional periodic substrate. We find various pinned and sliding phases as a function of substrate strength and drive amplitude. When the substrate is weak, a pinned bubble phase appears that depins elastically into a sliding bubble lattice. For stronger substrates, we find anisotropic bubbles, disordered bubbles, and stripe phases. Plastic depinning occurs via the hopping of individual particles from one bubble to the next in a pinned bubble lattice, and as the drive increases, there is a transition to a state where all of the bubbles are moving but are continuously shedding and absorbing individual particles. This is followed at high drives by a moving bubble lattice in which the particles can no longer escape their individual bubbles. The transition between the plastic and elastic sliding phases can be detected via signatures in the velocity-force curves, differential conductivity, and noise. When the bubbles shrink due to an increase in the attractive interaction term, they fit better inside the pinning troughs and become more strongly pinned, leading to a reentrant pinning phase. For weaker attractive terms, the size of the bubbles becomes greater than the width of the pinning troughs and the depinning becomes elastic with a reduced depinning threshold.

36 MATERIALS SCIENCE↗