Spatial Distribution Of Bubbles Research Articles

On the linear viscoelastic behavior of semidilute polydisperse bubble suspensions in Newtonian media

In this work, we investigated the linear viscoelasticity of semidilute polydisperse bubble suspensions via small amplitude oscillatory shear (SAOS) tests performed in a rheo-optical setup. For all tested suspensions, the measured viscoelastic moduli (G′, G″) aligned with the theoretical predictions of the Jeffreys model for average dynamic capillary numbers (⟨Cd⟩) greater than unity. But at lower ⟨Cd⟩ values, experimental G′ values exceeded theoretical predictions. To investigate this, we considered the effects of suspension polydispersity and various SAOS measurement artifacts, including bubble rise, coalescence, and changes in suspension microstructure over time. Polydispersity could not cause the observed deviation, because the viscoelastic trends deviate from those of classic single-mode relaxation only for bimodal bubble size distributions with equal volume fractions of very small and very large bubbles; in any other case, the polydisperse suspension behaves as monodisperse with a bubble radius equal to the volume-weighted mean radius. Furthermore, bubble rise proved to play a minor role, while SAOS rheo-optical experiments revealed that bubble size and organization varied negligibly during our measurements. The G′ deviation at low ⟨Cd⟩ values was linked to bubble fluid dynamic interactions induced by the bubble spatial distribution. Image analysis showed that at low bubble volume fractions, stronger and prolonged preshearing reduces these interactions by increasing the average interbubble distance. But this effect is negligible in denser suspensions, which show similar G′ trends for any applied preshearing. Finally, a multimode Jeffreys model fitted to the experimental data showed that bubble interactions complicate the relaxation process, introducing multiple relaxation modes.

Shock wave induced by the collapse of a bubble cluster with each bubble distributed randomly

The phenomenon of bubble collapse is ubiquitous in both natural and industrial settings, often manifesting in the form of bubble clusters. Despite its prevalence, there has been limited research on the dynamics of bubble clusters. In order to gain a better understanding of shock wave emission, liquid jet propagation, and the resulting pressure loads during bubble cluster collapse, we have developed a high-fidelity numerical approach aimed at accurately capturing shock waves and phase interfaces. The interface compression technique is introduced into the compressible two-phase model to improve the sharpness of the phase interface. We solve the physical model using the BVD (Boundary Variation Diminishing) principle. Both the 5th order WENO (Weighted Essentially Non-Oscillatory) scheme and the interface capturing function THINC (Tangent of Hyperbola for INterface Capturing) are utilized to reconstruct primary variables at cell boundaries, effectively minimizing numerical dissipation. For further resolution improvement, we implement these numerical methods on the AMR (Adaptive Mesh Refinement) mesh. Our numerical results demonstrate that the distance between the wall and the bubble cluster, as well as the nearest bubble to the wall, significantly influence the wall peak pressure. However, they have a minimal impact on the overall evolution of the bubble cluster and the peak pressure induced in the flow field. Remarkably, the spatial distribution of bubbles and the total number of bubbles play a crucial role in shaping the dynamic behavior of the bubble cluster and the pressure loads exerted both on the solid wall and in the flow field. The crucial findings offer profound insights into a range of applications, including cavitation erosion, underwater explosions, and drug delivery.

Study on the fluidity of foamed alkali-activated slag cementitious material (AASCM)

This study aims to investigate the evolution patterns of fluidity and rheological properties of AASCM under varying dosages of foaming agent and particle sizes of filling aggregate. The flow characteristics of AASCM are significantly affected by the filling aggregate's size and the foaming agent's dosage. Specifically, an increase in filling aggregate size (D(4,3) ϵ [26 μm, 69 μm]) enhances the fluidity of foamed AASCM, while an increase in foaming agent dosage reduces fluidity. These observed variations can be attributed to the presence of particle voids, the specific surface area of the aggregate, as well as the quantity and spatial distribution of bubbles within the slurry. A bubble-particle packing model is established, and by calibrating the simulation error coefficient to 1.1, the study investigates the evolution of water film thickness (WFT) in foamed AASCM with slurry expansion degree. It is observed that bubbles in the slurry affect the fluidity by altering the overall compactness and specific surface area of the foamed slurry, subsequently modifying the WFT.

Air-bubble entrainment by translating turbulent jets in stagnant water

This paper presents an experimental study on air-bubble entrainment in quiescent water by translating circular turbulent jets. The jet diameters, plunging heights, impact velocities, and translating velocities were varied during the experiments to investigate their relative effect on the bubble characteristics. The experimental observations reveal that the jet translation affects the air bubble entrainment mechanism and the distribution of bubble size. A rotating cavity forms around the plunging jet due to the translation of the jet. Depending on the translating velocity, the air bubble emanates from the cusp of the cavity and the downstream water surface meniscus with the jet. The bubble swarm produced by the translating jet exhibits vortex shedding with Strouhal numbers between 0.22 and 0.27, comparable to a circular cylinder in cross-flow. The peak of the bubble size distribution varies between 0.5 and 1.5 mm, while the Sauter mean diameter varies between 1.8 and 2.8 mm. The maximum penetration depth of bubbles is found to be a function of the jet impact to translating velocity ratio, and the Capillary number of the air–water interface. The spatial distribution of bubbles along the plume cross section exhibits Gaussian distributions. Finally, the terminal rising velocity of the bubbles shows no obvious effect of the jet translation.

Effects of pressure and particle size on bubble behaviors in a pseudo 2D pressured fluidized bed with Geldart A/B, B and D particles

Bubble behaviors of three types of particles from Geldart A/B, B and D groups were comprehensively investigated in a pseudo 2D pressurized fluidized bed via digital image analysis. Based on the analysis of a large quantity of experimental data, the effects of pressure and particle size on the spatial distribution of bubbles, bubble diameter, rise velocity, bubble fraction and flow regimes were obtained. Results show that bubbles are smaller and uniformly dispersed throughout the bed of Geldart A/B particles, while they become densely populated and grow into much larger sizes in the bed of Geldart B particles. A single peak of large bubbles in the center and emptying of bubble near the walls is formed in the bed of Geldart D particles due to the preferable rising path of bubbles towards the center and the significant coalescence of bubbles. Bubble diameter decreases with increasing pressure, and the effect of pressure is more pronounced on Geldart D particles. Bubble rise velocity increases with the increasing particle size and decreases as pressure rises. Large bubbles/slugs hardly appear in the bed of Geldart A/B particles, while they are clearly visible in the bed of Geldart B and D particles. The flow regime transition from slugging to fast fluidization shifts to the lower gas velocity for Geldart D particles compared with Geldart B particles.

An Experimental Model of Unconfined Bubbly Lava Flows: Importance of Localized Bubble Distribution

AbstractMost lava flows carry bubbles and crystals in suspension. From earlier works, it is known that spherical bubbles increase the effective viscosity while bubbles deformed by rapid flow decrease it. Changes in the spatial distribution of bubbles can lead to variable rheology and flow localization and thus modify the resulting lava flow structure and morphology. To understand the roles of bubble and solid phase crystal distributions, we conducted a series of analog experiments of high bubble fraction suspensions. We poured the analog lava on an inclined slope, observed its shape, calculated the velocity field, and monitored its local thickness. A region of localized rapid flow and low vesicularity, whose thickness is thinner than the surrounding area, develops at the center of the bubbly flows. These features suggest that the locally higher liquid fraction decreases the effective viscosity, increases the fluid density, and accelerates the flow. We also found that a halted particle‐bearing bubbly flow can resume flowing. We interpret this to result from the upward vertical separation of bubbles, which generates a liquid‐rich layer at the bottom of the flow. In our experiment, bubbles are basically spherical and decrease the flow velocity, while our estimate suggests that bubbles in natural lava flows could increase or decrease flow velocity. Downstream decreases in flow velocity stops the bubble deformation and can cause a sudden increase of effective viscosity. The vertical segregation of the liquid phase at the slowed flow front may be a way to generate a cavernous shelly paho’eho’e.

Ocean bubbles under high wind conditions – Part 1: Bubble distribution and development

Abstract. The bubbles generated by breaking waves are of considerable scientific interest due to their influence on air–sea gas transfer, aerosol production, and upper ocean optics and acoustics. However, a detailed understanding of the processes creating deeper bubble plumes (extending 2–10 m below the ocean surface) and their significance for air–sea gas exchange is still lacking. Here, we present bubble measurements from the HiWinGS expedition in the North Atlantic in 2013, collected during several storms with wind speeds of 10–27 m s−1. A suite of instruments was used to measure bubbles from a self-orienting free-floating spar buoy: a specialised bubble camera, acoustical resonators, and an upward-pointing sonar. The focus in this paper is on bubble void fractions and plume structure. The results are consistent with the presence of a heterogeneous shallow bubble layer occupying the top 1–2 m of the ocean, which is regularly replenished by breaking waves, and deeper plumes which are only formed from the shallow layer at the convergence zones of Langmuir circulation. These advection events are not directly connected to surface breaking. The void fraction distributions at 2 m depth show a sharp cut-off at a void fraction of 10−4.5 even in the highest winds, implying the existence of mechanisms limiting the void fractions close to the surface. Below wind speeds of 16 m s−1 or a wind-wave Reynolds number of RHw=2×106, the probability distribution of void fraction at 2 m depth is very similar in all conditions but increases significantly above either threshold. Void fractions are significantly different during periods of rising and falling winds, but there is no distinction with wave age. There is a complex near-surface flow structure due to Langmuir circulation, Stokes drift, and wind-induced current shear which influences the spatial distribution of bubbles within the top few metres. We do not see evidence for slow bubble dissolution as bubbles are carried downwards, implying that collapse is the more likely termination process. We conclude that the shallow and deeper bubble layers need to be studied simultaneously to link them to the 3D flow patterns in the top few metres of the ocean. Many open questions remain about the extent to which deep bubble plumes contribute to air–sea gas transfer. A companion paper (Czerski et al., 2022) addresses the observed bubble size distributions and the processes responsible for them.

Formation of Bubble-Loop Complexes During Helium Radiation in Fe-9Cr Steel

In the present study, the Fe-9Cr model alloy was irradiated with 240 keV He2+ at 550°C with a dose of 0.5 dpa at the peak damage region. The depth distribution of bubbles in Fe-9Cr alloy was investigated by transmission electron microscopy (TEM). The experimental results revealed that the spatial distribution of bubbles along the depth is different. In the region with higher helium concentration and irradiation dose, the bubbles inclined to be situated inside the plane of loops, forming a structure of “bubble-loop complex.” However, in regions where the helium concentration and irradiation dose are relatively low, the number of “bubble-loop complexes” significantly decreased. In addition, the Burgers vector of “bubble-loop complexes” was identified as <100> type. Radiation-induced enrichment of Cr atoms at the “bubble-loop complexes” was also quantitatively estimated by energy-dispersive X-ray spectroscopy (EDS) in the scanning TEM mode.

Effects of temperature on helium bubble behaviour in Fe–9Cr alloy

In the present study, Fe–9Cr model alloys were irradiated by 275 keV He+ ions at 500 to 800 °C to a peak damage level and implanted helium content of 0.5 dpa and 0.75 at% He, respectively. Effects of temperature on He bubble behaviour were studied by transmission electron microscopy (TEM) to reveal the temperature effects. The bubble spatial distribution changes from homogenous at 500 and 600 °C to heterogeneous at 700 and 800 °C, due to preferential bubble precipitation on dislocations and grain boundaries at high temperature. Faceted bubbles with low energy faces along the {001} planes were present at all four temperatures and are attributed to surface energy minimization. Bubble precipitation on the GBs became more pronounced with increasing temperature. Meanwhile, the cavity swelling increased with increasing temperature, from 0.59% at 500 °C to 8.21% at 800 °C with a sharp increase above 700 °C. By analysing the temperature dependence of bubble size and density, a transition temperature between low and high temperature regimes was obtained in terms of different activation energies which are related to different bubble nucleation and formation mechanisms.

Experimental investigation on the motion of bubbles in the moonpool with a sonar model inside

Bubbles are becoming a threat to the measurement of equipment and human activities in the moonpool, while little is known about the characteristics of bubbles in the moonpool. We aim to investigate the motion patterns and characteristic of bubbles in the moonpool. To support the investigation, experiments were conducted in a water tunnel. Images of bubbles in and under the moonpool were captured by two high-speed cameras. The software Tema and program PolyParticleTracker were used to track the bubbles in the moonpool. We obtain information of bubbles in the moonpool, including movement patterns of bubbles, generation of bubbles, relationship between amount of bubbles and the freesurface, temporal and spatial distribution of bubbles, the size and the velocities of bubbles.

Development of an SP3 neutron transport solver for the analysis of the Molten Salt Fast Reactor

The aim of this paper is the extension of a multiphysics OpenFOAM solver for the analysis of the Molten Salt Fast Reactor (MSFR), developed in previous works (Cervi et al., 2017, 2018). In particular, the neutronics sub-solver is improved by implementing a new module based on the SP3 approximation of the neutron transport equation. The new module is successfully tested against a Monte Carlo model of the MSFR, in order to assess its correct implementation. Then, a neutronics analysis of the MSFR is carried out on a simplified axial-symmetric model of the reactor. Particular focus is devoted to the analysis of the MSFR helium bubbling system and its effect on reactivity. The presence of bubbles inside the reactor is handled with a two-fluid thermal-hydraulics module, previously implemented into the solver. The void reactivity coefficient is evaluated on the basis of the bubble spatial distribution calculated by the multiphysics solver. Then, the results are compared to simulations carried out with uniform bubble distributions, highlighting significant differences between the two approaches. The outcomes of this work constitute a step forward in the multiphysics analysis of the Molten Salt Fast Reactor and represent a useful starting point for the optimization of the MSFR helium bubbling system, as well as for the development of appropriate control strategies.

Visualizing Probabilistic Multi‐Phase Fluid Simulation Data using a Sampling Approach

AbstractEulerian Method of Moment (MoM) solvers are gaining popularity for multi‐phase CFD simulation involving bubbles or droplets in process engineering. Because the actual positions of bubbles are uncertain, the spatial distribution of bubbles is described by scalar fields of moments, which can be interpreted as probability density functions. Visualizing these simulation results and comparing them to physical experiments is challenging, because neither the shape nor the distribution of bubbles described by the moments lend themselves to visual interpretation.In this work, we describe a visualization approach that provides explicit instances of the bubble distribution and produces bubble geometry based on local flow properties. To facilitate animation, the instancing of the bubble distribution provides coherence over time by advancing bubbles between time steps and updating the distribution. Our approach provides an intuitive visualization and enables direct visual comparison of simulation results to physical experiments.

Study of the effect of reactor scale on fluidization hydrodynamics using fine-grid CFD simulations based on the two-fluid model

Reliable scale-up of fluidized beds is essential to ensure that analysis and performance optimization at lab-scale can be applied to commercial scales. However, scaling fluidized beds for dynamic similarity continues to be challenging because flow hydrodynamics at lab-scale are largely influenced by bed geometry making extrapolation of conclusions to large-scales infeasible. Therefore, this study is focused on analyzing the effect of bed geometry on the fluidization hydrodynamics using large-scale CFD simulations. The two fluid model (TFM) is employed to describe the solids motion efficiently and simulations are conducted for fluidization of 1150μm LLDPE and 500μm glass beads in beds of different sizes (diameter D=15–70cm and initial bed height H0=10–75cm). The hydrodynamics are subsequently investigated qualitatively using time-resolved visualizations, bubble centroid and solids velocity maps as well as quantitatively using detailed bubble statistics and solids circulation metrics. It is shown that as the bed diameter is increased, average bubble sizes decrease although similar-sized bubbles rise faster because of lower wall resistance, both factors contributing to faster solids circulation. On the other hand, fluidization hydrodynamics in 50cm diameter bed are relatively insensitive to the choice of H0 and similarities in solids circulation patterns are observed in shallow beds as well as in the lower regions of deep beds. Finally, it is shown that the size and spatial-distribution of bubbles is crucial for maintaining dynamic similarity of bubbling beds. Specifically, the bed dimensions (D, H0) must ensure that (a) bubbles are typically much smaller than the bed diameter and (b) solids circulation patterns are similar across scales of interest. Overall, insights from this study can be used for describing the gas distribution and solids motion more accurately for better design of commercial beds.

CFD investigation of bubble effects on Taylor–Couette flow patterns in the weakly turbulent vortex regime

The effect of gas bubbles on Taylor–Couette flow patterns in the transition from laminar to weakly turbulent vortex flow is investigated by quasi-2D and 3D CFD simulation. The CFD model is based on the Eulerian–Eulerian two-fluid approach with constitutive closures for inter-phase forces and liquid turbulence and validated by comparison with previously published experimental data. The bubble spatial distribution patterns and the mechanism by which the dispersed fluid obtains a non-uniform radial distribution is discussed and the relative importance of various inter-phase forces is demonstrated. In addition, the experimental observation of vortex axial wavelength expansion due to the presence of bubbles was reproduced by the CFD simulations, and a two fluid Taylor–Couette flow regime map based on three observed regimes was constructed. Lastly, the influence of the azimuthal Reynolds number and gas superficial velocity on vortex wavelength, liquid velocity, mass transfer and turbulence are discussed in detail.

Effect of electrostatics on freely-bubbling beds of mono-sized particles

Electrostatic charges influence single bubble shape and rise velocity, as shown in our previous study. This paper extends the previous work to investigate the influence of electrostatic charges on bubbles in freely-bubbling fluidized beds, by simulating a two-dimensional column, with charged mono-sized glass beads supported by air. Simulation results are first compared with experiments in an uncharged system at two superficial gas velocities. The sensitivity of the results to frictional models is investigated, with the model chosen for the rest of the study which gives better predictions. In the second part of this work, the effect of electrostatics on hydrodynamics of freely-bubbling fluidized beds is investigated by comparing bubble diameter, distribution and time-average vertical solid velocities for uncharged and charged cases. The results predict that electrostatic charges decrease bubble size, modify bubble spatial distribution and influence time-average solid velocities.

The effects of nitroglycerin, norepinephrine and aminophylline on intrapulmonary arteriovenous anastomoses in healthy humans at rest

We have investigated the effects of the intravenous infusion of nitroglycerin (NTG), norepinephrine (NE) and aminophylline (AMP) on the opening and recruitment of intrapulmonary arteriovenous anastomoses (IPAVA) in healthy humans at rest. In ten volunteers saline contrast echocardiography was performed during administration of two doses of the NTG (3μgkg−1min−1 and 6μgkg−1min−1) and NE (0.1μgkg−1min−1 and 0.25μgkg−1min−1) as well as 30min following the administration of AMP at rate of 6mgkg−1. Echocardiography was used to assign bubble scores (0–5) based on the number and spatial distribution of bubbles in the left ventricle. Doppler ultrasound was used to estimate pulmonary artery systolic pressure. Using a Finometer the following hemodynamic parameters were assessed: heart rate, stroke volume, cardiac output, total peripheral resistance as well as systolic, diastolic and mean arterial pressure. The most important finding from the current study was that nitroglycerin, norepinephrine and aminophylline in the applied doses were not found to promote IPAVA opening in healthy humans at rest.

Numerical simulation of bubble dispersion in turbulent Taylor-Couette flow

We investigate bubble dispersion in turbulent Taylor-Couette flow. The aim of this study is to describe the main mechanisms yielding preferential bubble accumulation in near-wall structures of the flow. We first proceed to direct numerical simulation of Taylor-Couette flows for three different geometrical configurations (three radius ratios η = R1/R2: η = 0.5, η = 0.72, and η = 0.91 with the outer cylinder at rest) and Reynolds numbers corresponding to turbulent regime ranging from 3000 to 8000. The statistics of the flow are discussed using two different averaging procedures that permit to characterize the mean azimuthal velocity, the Taylor vortices contribution and the small-scale turbulent fluctuations. The simulations are compared and validated with experimental and numerical data from literature. The second part of this study is devoted to bubble dispersion. Bubble accumulation is analyzed by comparing the dispersion obtained with the full turbulent flow field to bubble dispersion occurring at lower Reynolds numbers in previous works. Several patterns of preferential accumulation of bubbles have been observed depending on bubble size and the effect of gravity. For the smaller size considered, bubbles disperse homogeneously throughout the gap, while for the larger size they accumulate along the inner wall for the large gap width (η = 0.5). Varying the intensity of buoyancy yields complex evolution of the bubble spatial distribution. For low gravity effect, bubble entrapment is strong leading to accumulation along the inner wall in outflow regions (streaks of low wall shear stress). When buoyancy effect dominates on vortex trapping, bubbles rise through the vortices, while spiral patterns stretched along the inner cylinder are clearly identified. Force balance is analyzed to identify dominating forces leading to this accumulation and accumulation patterns are compared with previous experiments.

The size range of bubbles that produce ash during explosive volcanic eruptions

Volcanic eruptions can produce ash particles with a range of sizes and morphologies. Here we morphologically distinguish two textural types: Simple (generally smaller) ash particles, where the observable surface displays a single measureable bubble because there is at most one vesicle imprint preserved on each facet of the particle; and complex ash particles, which display multiple vesicle imprints on their surfaces for measurement and may contain complete, unfragmented vesicles in their interiors. Digital elevation models from stereo-scanning electron microscopic images of complex ash particles from the 14 October 1974 sub-Plinian eruption of Volcán Fuego, Guatemala and the 18 May 1980 Plinian eruption of Mount St. Helens, Washington, U.S.A. reveal size distributions of bubbles that burst during magma fragmentation. Results were compared between these two well-characterized eruptions of different explosivities and magma compositions and indicate that bubble size distributions (BSDs) are bimodal, suggesting a minimum of two nucleation events during both eruptions. The larger size mode has a much lower bubble number density (BND) than the smaller size mode, yet these few larger bubbles represent the bulk of the total bubble volume. We infer that the larger bubbles reflect an earlier nucleation event (at depth within the conduit) with subsequent diffusive and decompressive bubble growth and possible coalescence during magma ascent, while the smaller bubbles reflect a relatively later nucleation event occurring closer in time to the point of fragmentation. Bubbles in the Mount St. Helens complex ash particles are generally smaller, but have a total number density roughly one order of magnitude higher, compared to the Fuego samples. Results demonstrate that because ash from explosive eruptions preserves the size of bubbles that nucleated in the magma, grew, and then burst during fragmentation, the analysis of the ash-sized component of tephra can provide insights into the spatial distribution of bubbles in the magma prior to fragmentation, enabling better parameterization of numerical eruption models and improved understanding of ash transport phenomena that result in pyroclastic volcanic hazards. Additionally, the fact that the ash-sized component of tephra preserves BSDs and BNDs consistent with those preserved in larger pyroclasts indicates that these values can be obtained in cases where only distal ash samples from particular eruptions are obtainable.

Microbubble behavior in an ultrasound field for high intensity focused ultrasound therapy enhancement

The enhancement of heating due to inertial cavitation has been focused to reduce the long treatment time of conventional high-intensity focused ultrasound (HIFU) therapy. The influences of the physical properties of surrounding tissues, initial void fraction, and spatial distribution of bubbles on microbubble-enhanced HIFU are examined. A bubble dynamics equation based on the Keller-Miksis equation is employed in consideration of the elasticity of surrounding tissue. The mixture phase and bubbles are coupled by the Euler-Lagrange method to take into account the interaction between ultrasound and bubbles. As a result, the temperature around the target increases with the initial void fraction. But at the high void fraction of 10(-5), ultrasound is too attenuated to heat the target, and the heating region moves to the transducer side. On the other hand, both the viscosity and shear elasticity of the surrounding media reduce the attenuation of ultrasound propagation through the bubbly mixture. Numerical results show that localized heating is induced with increasing viscosity or shear elasticity, though it depends on the pressure amplitudes. In addition, it was numerically confirmed that the localization of the microbubble distribution is important to obtain efficient localized heating.

Catecholamine-induced opening of intrapulmonary arteriovenous anastomoses in healthy humans at rest

The mechanism or mechanisms that cause intrapulmonary arteriovenous anastomoses (IPAVA) to either open during exercise in subjects breathing room air and at rest when breathing hypoxic gas mixtures, or to close during exercise while breathing 100% oxygen, remain unknown. During conditions when IPAVA are open, plasma epinephrine (EPI) and dopamine (DA) concentrations both increase, potentially representing a common mechanism. The purpose of this study was to determine whether EPI or DA infusions open IPAVA in resting subjects breathing room air and, subsequently, 100% oxygen. We hypothesized that these catecholamine infusions would open IPAVA. We performed saline-contrast echocardiography in nine subjects without a patent foramen ovale before and during serial EPI and DA infusions while breathing room air and then while breathing 100% oxygen. Bubble scores (0-5) were assigned based on the number and spatial distribution of bubbles in the left ventricle. Pulmonary artery systolic pressure (PASP) was estimated using Doppler ultrasound, while cardiac output (Q(C)) was measured using echocardiography. Bubble scores were significantly greater during EPI infusions of 80-320 ng·kg(-1)·min(-1) compared with baseline when subjects breathed room air; however, bubble scores did not increase when they breathed 100% oxygen. At comparable Q(C) and PASP, intravenous DA (16 μg·kg(-1)·min(-1)) and EPI (40 ng·kg(-1)·min(-1)) resulted in identical bubble scores. Subsequent studies revealed that β-blockade did not prevent hypoxia-induced opening of IPAVA. We suggest that increases in Q(C) or PASP (or both) secondary to EPI or DA infusions open IPAVA in normoxia. The closing mechanism associated with breathing 100% oxygen is independent from the opening mechanisms.