Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BUBBLES”

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

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)↗
DOE OSTI2012

THE FERMI BUBBLES: SUPERSONIC ACTIVE GALACTIC NUCLEUS JETS WITH ANISOTROPIC COSMIC-RAY DIFFUSION

The Fermi Gamma-Ray Space Telescope reveals two large bubbles in the Galaxy, which extend nearly symmetrically {approx}50 Degree-Sign above and below the Galactic center. Using three-dimensional (3D) magnetohydrodynamic simulations that self-consistently include the dynamical interaction between cosmic rays (CRs) and thermal gas and anisotropic CR diffusion along the magnetic field lines, we show that the key characteristics of the observed gamma-ray bubbles and the spatially correlated X-ray features in the ROSAT 1.5 keV map can be successfully reproduced by recent jet activity from the central active galactic nucleus. We find that after taking into account the projection of the 3D bubbles onto the sky the physical heights of the bubbles can be much smaller than previously thought, greatly reducing the formation time of the bubbles to about a Myr. This relatively small bubble age is needed to reconcile the simulations with the upper limit of bubble ages estimated from the cooling time of high-energy electrons. No additional physical mechanisms are required to suppress large-scale hydrodynamic instabilities because the evolution time is too short for them to develop. The simulated CR bubbles are edge-brightened, which is consistent with the observed projected flat surface brightness distribution. Furthermore, we demonstrate that the sharp edges of the observed bubbles can be due to anisotropic CR diffusion along magnetic field lines that drape around the bubbles during their supersonic expansion, with suppressed perpendicular diffusion across the bubble surface. Possible causes of the slight bends of the Fermi bubbles to the west are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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↗
DOE OSTI2025

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 ↗
DOE OSTI2019

Surfactant stabilized bubbles flowing in a Newtonian fluid

Bubbles suspended in a fluid cause the suspension to have different rheological properties than the base fluid. In general, the viscosity of the suspension increases as the volume fraction of the bubbles is increased. A current application, and motivation for this study, is in wellbore cements used for hydrocarbon extraction and carbon sequestration. In these settings, the gas bubbles are dispersed into the cement to reduce the density as well as improve the properties for specific conditions or wellbore issues. In this paper, we use Stokesian dynamics to numerically simulate the behavior of a large number of bubbles suspended in a Newtonian fluid. Going beyond prior work on simulating particles in suspension, we account for the nature of bubbles by allowing for slip on the bubble surface, the deflection on the bubble surface, and a bubble–bubble pairwise interaction that represents the surfactant physics; we do not account for bubble compressibility. We incorporate these interactions and simulate bubble suspensions of monodisperse size at several volume fractions. We find that the bubbles remain better dispersed compared with hard spherical particles that show a greater tendency to structure or cluster.

Rosenbaum, Eilis↗
DOE OSTI2020

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↗
DOE OSTI2024

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↗
DOE OSTI2016

MOLECULAR LINES OF 13 GALACTIC INFRARED BUBBLE REGIONS

We investigated the physical properties of molecular clouds and star formation (SF) processes around infrared bubbles, which are essentially expanding H ii regions. We performed observations of 13 galactic infrared bubble fields containing 18 bubbles. We observed five molecular lines—{sup 12}CO (J=1→0), {sup 13}CO (J=1→0), C{sup 18}O (J=1→0), HCN (J=1→0), and HCO{sup +} (J=1→0)—and several publicly available surveys were used for comparison: Galactic Legacy Infrared Mid-Plane Survey Extraordinaire, Multiband Imaging Photometer for Spitzer Galactic Plane Survey, APEX Telescope Large Area Survey of the Galaxy, Bolocam Galactic Plane Survey, Very Large Array (VLA) Galactic Plane Survey, Multi-Array Galactic Plane Imaging Survey, and NRAO VLA Sky Survey. We find that these bubbles are generally connected with molecular clouds, most of which are giant. Several bubble regions display velocity gradients and broad-shifted profiles, which could be due to the expansion of bubbles. The masses of molecular clouds within bubbles range from 100 to 19,000 M {sub ⊙}, and their dynamic ages are about 0.3–3.7 Myr, which takes into account the internal turbulence pressure of surrounding molecular clouds. Clumps are found in the vicinity of all 18 bubbles, and molecular clouds near four of these bubbles with larger angular sizes show shell-like morphologies, indicating that either collect-and-collapse or radiation-driven implosion processes may have occurred. Due to the contamination of adjacent molecular clouds, only six bubble regions are appropriate to search for outflows, and we find that four have outflow activities. Three bubbles display ultra-compact H ii regions at their borders, and one is probably responsible for its outflow. In total, only six bubbles show SF activities in the vicinity, and we suggest that SF processes might have been triggered.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

A study of the accuracy of neutrally buoyant bubbles used as flow tracers in air

Research has been performed to determine the accuracy of neutrally buoyant and near neutrally buoyant bubbles used as flow tracers in air. Theoretical, computational, and experimental results are presented to evaluate the dynamics of bubble trajectories and factors affecting their ability to trace flow-field streamlines. The equation of motion for a single bubble was obtained and evaluated using a computational scheme to determine the factors which affect a bubble's trajectory. A two-dimensional experiment was also conducted to experimentally determine bubble trajectories in the stagnation region of NACA 0012 airfoil at 0 deg angle of attack using a commercially available helium bubble generation system. Physical properties of the experimental bubble trajectories were estimated using the computational scheme. These properties included the density ratio and diameter of the individual bubbles. the helium bubble system was then used to visualize and document the flow field about a 30 deg swept semispan wing with simulated glaze ice. Results were compared to Navier-Stokes calculations and surface oil flow visualization. The theoretical and computational analysis have shown that neutrally buoyant bubbles will trace even the most complex flow patterns. Experimental analysis revealed that the use of bubbles to trace flow patterns should be limited to qualitative measurements unless care is taken to ensure neutral buoyancy. This is due to the difficulty in the production of neutrally buoyant bubbles.

Kerho, Michael F.↗
NASA NTRS1999

Buoyancy Driven Shear Flows of Bubble Suspensions

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Koch, D. L.↗
NASA NTRS2003

Prediction of Bubble Diameter at Detachment from a Wall Orifice in Liquid Cross Flow Under Reduced and Normal Gravity Conditions

Bubble formation and detachment is an integral part of the two-phase flow science. The objective of the present work is to theoretically investigate the effects of liquid cross-flow velocity, gas flow rate embodied in the momentum flux force, and orifice diameter on bubble formation in a wall-bubble injection configuration. A two-dimensional one-stage theoretical model based on a global force balance on the bubble evolving from a wall orifice in a cross liquid flow is presented in this work. In this model, relevant forces acting on the evolving bubble are expressed in terms of the bubble center of mass coordinates and solved simultaneously. Relevant forces in low gravity included the momentum flux, shear-lift, surface tension, drag and inertia forces. Under normal gravity conditions, the buoyancy force, which is dominant under such conditions, can be added to the force balance. Two detachment criteria were applicable depending on the gas to liquid momentum force ratio. For low ratios, the time when the bubble acceleration in the direction of the detachment angle is greater or equal to zero is calculated from the bubble x and y coordinates. This time is taken as the time at which all the detaching forces that are acting on the bubble are greater or equal to the attaching forces. For high gas to liquid momentum force ratios, the time at which the y coordinate less the bubble radius equals zero is calculated. The bubble diameter is evaluated at this time as the diameter at detachment from the fact that the bubble volume is simply given by the product of the gas flow rate and time elapsed. Comparison of the model s predictions was also made with predictions from a two-dimensional normal gravity model based on Kumar-Kuloor formulation and such a comparison is presented in this work.

Nahra, Henry K.↗
NASA NTRS2002

Mathematical model of diffusion-limited gas bubble dynamics in unstirred tissue with finite volume

Models of gas bubble dynamics for studying decompression sickness have been developed by considering the bubble to be immersed in an extravascular tissue with diffusion-limited gas exchange between the bubble and the surrounding unstirred tissue. In previous versions of this two-region model, the tissue volume must be theoretically infinite, which renders the model inapplicable to analysis of bubble growth in a finite-sized tissue. We herein present a new two-region model that is applicable to problems involving finite tissue volumes. By introducing radial deviations to gas tension in the diffusion region surrounding the bubble, the concentration gradient can be zero at a finite distance from the bubble, thus limiting the tissue volume that participates in bubble–tissue gas exchange. It is shown that these deviations account for the effects of heterogeneous perfusion on gas bubble dynamics, and are required for the tissue volume to be finite. The bubble growth results from a difference between the bubble gas pressure and an average gas tension in the surrounding diffusion region that explicitly depends on gas uptake and release by the bubble. For any given decompression, the diffusion region volume must stay above a certain minimum in order to sustain bubble growth.

NASA Discipline Environmental Health↗
NASA NTRS2010

Bubble Departure from Metal-Graphite Composite Surfaces and Its Effects on Pool Boiling Heat Transfer

The formation and growth processes of a bubble in the vicinity of graphite micro-fiber tips on metal-graphite composite boiling surfaces and their effects on boiling behavior are investigated. It is discovered that a large number of micro bubbles are formed first at the micro scratches and cavities on the metal matrix in pool boiling. By virtue of the non-wetting property of graphite, once the growing micro bubbles touch the graphite tips, the micro bubbles are sucked by the tips and merged into larger micro bubbles sitting on the end of the tips. The micro bubbles grow rapidly and coalesce to form macro bubbles, each spanning several tips. The necking process of a detaching macro bubble is analyzed. It is revealed that a liquid jet is produced by sudden break-off of the bubble throat. The composite surfaces not only have higher temperatures in micro- and macrolayers but also make higher frequency of the bubble departure, which increase the average heat fluxes in both the bubble growth stage and in the bubble departure period. Based on these analyses, the enhancement mechanism of pool boiling heat transfer on composite surfaces is clearly revealed.

Chao, David F.↗
NASA NTRS2004

Using an Ultrasonic Instrument to Size Extravascular Bubbles

In an ongoing development project, microscopic bubbles in extravascular tissue in a human body will be detected by use of an enhanced version of the apparatus described in Ultrasonic Bubble- Sizing Instrument (MSC-22980), NASA Tech Briefs, Vol. 24, No. 10 (October 2000), page 62. To recapitulate: The physical basis of the instrument is the use of ultrasound to excite and measure the resonant behavior (oscillatory expansion and contraction) of bubbles. The resonant behavior is a function of the bubble diameter; the instrument exploits the diameter dependence of the resonance frequency and the general nonlinearity of the ultrasonic response of bubbles to detect bubbles and potentially measure their diameters. In the cited prior article, the application given most prominent mention was the measurement of gaseous emboli (essentially, gas bubbles in blood vessels) that cause decompression sickness and complications associated with cardiopulmonary surgery. According to the present proposal, the instrument capabilities would be extended to measure extravascular bubbles with diameters in the approximate range of 1 to 30 m. The proposed use of the instrument could contribute further to the understanding and prevention of decompression sickness: There is evidence that suggests that prebreathing oxygen greatly reduces the risk of decompression sickness by reducing the number of microscopic extravascular bubbles. By using the ultrasonic bubble-sizing instrument to detect and/or measure the sizes of such bubbles, it might be possible to predict the risk of decompression sickness. The instrument also has potential as a tool to guide the oxygen-prebreathing schedules of astronauts; high-altitude aviators; individuals who undertake high-altitude, low-opening (HALO) parachute jumps; and others at risk of decompression sickness. For example, an individual at serious risk of decompression sickness because of high concentrations of extravascular microscopic bubbles could be given a warning to continue to prebreathe oxygen until it was safe to decompress.

Magari, Patrick J.↗
DOE OSTI2025

Driven probe particle dynamics in a bubble and pattern forming system

We numerically examine the dynamics of a probe particle driven at a constant force through an assembly of particles with competing long-range repulsion and short-range attraction that forms a bubble or stripe state. In the bubble regime, we identify several distinct types of motion, including an elastic or pinned regime where the probe particle remains inside a bubble and drags all other bubbles with it. There is also a plastic bubble phase where the bubble in which the probe particle is trapped is able to move past the adjacent bubbles. At larger drives, there is a breakthrough regime where the probe particle jumps from bubble to bubble and, in some cases, can induce correlated rotations or plastic rearrangements of the particles within the bubbles. At the highest drives, the probe particle moves sufficiently rapidly that the background particles undergo only small distortions. The distinctive dynamic flow states and the transitions between them are accompanied by signatures in the effective drag on the driven particle, jumps in the velocity–force curves, and changes in the time-dependent velocity fluctuations. We map the dynamic phase diagram for this system for varied interaction lengths, bubble sizes, and densities.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Measurement of interfacial structures in horizontal air-water bubbly flows

In order to predict multi-dimensional phenomena in nuclear reactor systems, methods relying on computational fluid dynamics (CFD) codes are essential. However, to be applicable in assessing thermal-hydraulic safety, these codes must be able to accurately predict the development of two-phase flows. Therefore, before practical application these codes must be assessed using experimental databases that capture multi-dimensional phenomena. While a large database exists that can be employed to assess predictions in vertical flows, the available database for horizontal flows is significantly lacking. Therefore, the current work seeks to develop an additional database in air-water horizontal bubbly flow through a 38.1 mm ID test section with a total development length of approximately 250 diameters. The experimental conditions are chosen to cover a wide range of the bubbly flow regime based upon flow visualization using a high-speed video camera. A database of local time-averaged void fraction, bubble velocity, interfacial area concentration, and bubble Sauter mean diameter are acquired throughout the pipe cross-section using a four-sensor conductivity probe. To investigate the evolution of the flow, measurements are made at axial locations of 44, 116, and 244 diameters downstream of the inlet. In the current work, only measurements obtained at L/D = 244 are presented. It is found that increasing the liquid superficial velocity tends to reduce both the bubble size and the degree of bubble packing near the upper wall. However, it is observed that the position of the maximum void fraction value remains nearly constant and is located approximately one bubble diameter away from the upper wall. It is also found that the bubble velocity exhibits a power law behavior resembling a single phase liquid turbulent velocity profile. Moreover, the local bubble velocity tends to decrease as the local void fraction increases. Conversely, increasing the gas superficial velocity is found to have little effect on the bubble size or velocity, but increases the bubble packing near the upper pipe wall. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Helium bubble evolution in a Zr–Sn–Nb–Fe–Cr alloy during post-annealing: An in-situ investigation

The formation of helium bubbles is considered to be detrimental to the mechanical performance of the nuclear materials. The growth behaviors of helium bubbles in a helium ion implanted Zr–Sn–Nb–Fe–Cr alloy with respect to the helium fluence and subsequently annealing procedure were investigated by in-situ transmission electron microscopy. In the as-implanted sample, the measured size distributions of the helium bubbles are consistent with the simulated helium concentrations. Moreover, the mean size of the helium bubbles increases with the increase of the irradiation temperatures and the helium fluence. The in-situ heating study performed in a transmission electron microscope indicates that the mean size of the helium bubbles increase slowly below 923 K and dramatically above 923 K. The coarsening mechanism of the helium bubbles in the alloy is suggested based on the study. - Highlights: • Helium bubble growth in zirconium with annealing was in-situ investigated in TEM. • The mean helium bubble size increase with helium fluence and annealing temperature. • Helium bubble size distribution is same as that of helium concentration by SRIM. • Mean bubble size increases slowly and quickly with temperature below and above 923 K. • The growth mechanism of the helium bubbles in Zr alloy has been discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗