Search NASA⌕ Search

SEARCH · Search NASA

Results for “bubble dynamics”

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 37 records · Page 2

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↗

Recent progress in in-situ visualization of electrochemical reactions in electrochemical energy devices

A growing body of evidence indicates that bubble dynamics are crucial performance-limiting factors on electrochemical reactions and energy conversion efficiencies in both half-cell and full cell device levels. This review aims to summarize the recent progress in in-situ visualization characterizations of electrochemical reactions and reactant/product transport in electrochemical energy systems. The latest achievements of various new electrodes, transparent cell designs, and development for simultaneous visualization of reactions at both electrodes and electrolyte/electrode interface are highlighted. The new insights into the influence of bubbles on reaction sites and associated electrochemical performances are discussed. Lastly, some perspectives on potential strategies for mitigating the bubble-induced energy losses are proposed for future research directions and opportunities in water electrolysis and other electrochemical energy systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling multiple scattering transient of an ultrashort laser pulse by spherical particles

The multiple scattering of an ultrashort laser pulse by a turbid dispersive medium (namely a cloud of bubbles in water) is investigated by means of Monte Carlo simulations. The theory of Gouesbet and Gréhan (2000) is used to derive an energetic model of the scattering transient. It is shown that the spreading and extinction of the pulse can be decoupled from the transient of scattering, which allows to describe each phenomenon individually. The transient of scattering is modeled with the Lorenz-Mie Theory and thus is also valid for a relative refractive index lower than one, contrary to the Debye series expansion which does not converge close to the critical angle. This is made possible after the introduction of a new physical object, the Scattering Impulse Response Function (SIRF) which allows to detect the different modes of scattering transient, in time and direction. Here the present approach is more generic, as it enables to simulate clouds of air bubbles in water, which was not possible previously. Two different approaches are proposed within the Monte Carlo framework. The first is a pure Monte Carlo approach where the delay due to the scattering is randomly drawn at each event, while the second is based on the transport of the whole scattering signal. They are both embedded in the Monte Carlo code Scatter3D. Both models produce equivalent trends and are validated against published numerical results. They are then applied to the multiple scattering of ultra short pulse by a cloud of bubble in water in the forward direction. The pulse spread due to the propagation in water is computed for a wide range of traveled distances and pulse durations, and the optimal pulse duration is given to minimize the pulse spread at a given distance. The main result is that the scattered photons exit the turbid medium earlier than the ballistic photons and produce a double peak related to the refraction in the bubble. This demonstrates the possibility to develop new diagnostics to characterize dynamic bubbly flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase transitions from the fifth dimension

We study the cosmological transition of 5D warped compactifications, from the high-temperature black-brane phase to the low-temperature Randall-Sundrum I phase. The transition proceeds via percolation of bubbles of IR-brane nucleating from the black-brane horizon. The violent bubble dynamics can be a powerful source of observable stochastic gravitational waves. While bubble nucleation is non-perturbative in 5D gravity, it is amenable to semiclassical treatment in terms of a “bounce” configuration interpolating between the two phases. We demonstrate how such a bounce configuration can be smooth enough to maintain 5D effective field theory control, and how a simple ansatz for it places a rigorous lower-bound on the transition rate in the thin-wall regime, and gives plausible estimates more generally. When applied to the Hierarchy Problem, the minimal Goldberger-Wise stabilization of the warped throat leads to a slow transition with significant supercooling. We demonstrate that a simple generalization of the Goldberger-Wise potential modifies the IR-brane dynamics so that the transition completes more promptly. Supercooling determines the dilution of any (dark) matter abundances generated before the transition, potentially at odds with data, while the prompter transition resolves such tensions. We discuss the impact of the different possibilities on the strength of the gravitational wave signals. Via AdS/CFT duality the warped transition gives a theoretically tractable holographic description of the 4D Composite Higgs (de)confinement transition. Our generalization of the Goldberger-Wise mechanism is dual to, and concretely models, our earlier proposal in which the composite dynamics is governed by separate UV and IR RG fixed points. The smooth 5D bounce configuration we introduce complements the 4D dilaton/radion dominance derivation presented in our earlier work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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↗

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗

Impact of Variable Gas Mixtures on Bubble Size Distribution and Mass Transfer in Gas Fermentation Reactors

Gas fermentation has emerged as a promising new technology for the generation of fuels and chemicals from mixtures of greenhouse and energy rich gas streams (CO2/CH4/H2/CO) via microbial bioreaction. Example pathways include biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Effective Gas-liquid mass-transfer is an important physical phenomenon that determines the design and scale-up of these systems. There is currently a knowledge-gap regarding bubble-size distributions when using a mixture of gases with vastly different properties, which can have a significant impact on overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a detailed multiphase computational-fluid-dynamics (CFD) model to study the impact of gas-mixtures on overall mass-transfer in bubble column and air-lift reactors. The CFD tool previously developed by the authors (1) for simulating aerobic fermentation reactors at scale is used in this study. The Reynolds-averaged mass, momentum, energy, and species transport equations are solved for interpenetrating gas and liquid phase in this model. We use a population balance-based bubble-size-distribution model that is validated against small-scale experiments in our solver. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient is analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗

Multiphysics Computational-Fluid-Dynamics (CFD) for Design and Scale-Up of Gas Bioreactors that Utilize CO2

Gas fermentation technology offers sustainable and high carbon-efficiency routes to fuels and chemicals. However, efficiently delivering low-solubility substrate gases to production-scale bioreactors is a significant technical challenge that could limit productivity. We perform fundamental and applied multiphysics computational fluid dynamics (CFD) research targeting specific challenges of gas delivery in fermentation systems that utilize varied gas-compositions, for e.g. CO2/CO/CH4/H2 substrates: bubble dynamics, coupled mass-transfer and reaction, that can lead to optimal reactor designs at scale. This work will support BETO goals towards greenhouse gas reduction, derisking technologies toward commercialization/deployment and sustainable-aviation-fuel synthesis targets.

BIOMASS FUELS↗

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗

Molecular-Resolution Electron Imaging of Defects and Dynamics at the Ice-Water Interface

Water crystallization into hexagonal ice (type I h ) is one of the most critical processes relevant to the Earth’s environment and human activities. However, despite recent breakthroughs in imaging non-equilibrium condensed ice structures, the ice-water interface has never been imaged at a molecular resolution. This is primarily due to the low stability of the hydrogen bonds in ice under high-resolution microscopy conditions and a lack of methods to prepare compatible samples. Here, this presentation describes the first molecular-resolution imaging of ice crystallized from liquid water and the ice-water interface using high-resolution transmission electron microscopy (HRTEM). By encapsulating deionized (DI) water between two amorphous carbon (a-C) TEM grids and subsequently freezing it with liquid N 2 on a cryo sample stage, we generated two types of ice: non-equilibrium, condensed ice from the atmosphere and encapsulated ice from the DI water (Fig. 1A). Condensed ice usually shows irregular, spherulitic shapes (Fig. 1B). Selected area electron diffraction (SAED) shows that they are a mixture of cubic and hexagonal crystals (Fig. 1C). On the contrary, encapsulated ice forms thin films that contain large-area single-crystalline regions of hexagonal ice oriented along the [0001] zone axis (Fig. 1D). Differential electron energy-loss spectroscopy (EELS) confirmed the high purity of the encapsulated ice free from organic contaminations that are common in other encapsulation methods for HRTEM such as graphene liquid cells. These single-crystalline areas are robust under the electron beam up to ~100 e/Å 2 s. Aberration-corrected HRTEM imaging in these areas achieved a line resolution of ~1.3 Å (Fig. 1E and F). This platform allows us to study near-equilibrium ice structures and dynamics at an unprecedented spatial resolution (Fig. 2). For example, we discovered subdomain-rich regions near the defective crystal edges despite the structure appearing single-crystalline according to diffraction criteria. These subdomains connect via low-angle grain boundaries with flat energy landscapes as a function of tilt angles (according to simulations), showing the high tolerance of ice to defect structures. When we tuned the sample temperature and electron flux rate, we observed radiolysis-controlled bubble generation and dissolution in ice single crystals near a steady state of bubble dynamics. Furthermore, rich beam-induced melting and recrystallization dynamics were observed at the ice-water interface with lattice resolution. These data represent the first observation of the ice-water phase transformation at the sub-nanometer level. In summary, the methods developed in this work enabled molecular-resolution observations of ice and the ice-water interface and shed light on the microstructures and phase transformation pathways. Finer control on the temperature, electron irradiation profile, and imaging detector could eventually lead to real-time observation of ice nucleation in water and address long-standing questions in the nucleation pathways.

74 ATOMIC AND MOLECULAR PHYSICS↗

Tracking Local pH Dynamics during Water Electrolysis via In-Line Continuous Flow Raman Spectroscopy

The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.

Marquez, Raul A. [Univ. of Texas, Austin, TX (Unit↗

Modelling the ejection of primary aerosols during the fast pyrolysis of biomass anisotropic particles

A model for the fast pyrolysis of anisotropic biomass particles is presented which considers bubbling dynamics within the liquid intermediate phase (metaplast) and aerosol ejection from this phase. The model employs the population balance equation and the method of moments to estimate the production rate and resultant size distribution of aerosol ejections, incorporating a detailed CRECK reaction mechanism, and considers the effect of anisotropic biomass microstructure on the intraparticle transport of mass and energy. Here, this study investigates the impact of particle size, heating rate (heat transfer coefficient), and lignocellulosic composition on aerosol ejection. The model predicts that, at high heating rates (convective heat transfer coefficient of 359 W/m 2 .K), aerosols can contribute over 20% to the heavy fraction yield in bio-oil for small particles (1 mm diameter, 4 mm length). The model can predict aerosol size distribution and surface area, indicating an average size of 20 μm for bubbles and 5 μm for aerosols during increased bubble production and aerosol ejection rates. These findings are consistent with prior experimental results and provide essential information for future modeling of extra-particle reactions of the aerosols as they progress through the reactor.

09 BIOMASS FUELS↗

Favorable morphology and electronic conductivity of functional sublayers for highly efficient water splitting electrodes

Low electronic conductivities and improper morphologies of anode electrodes greatly limit the reaction area, catalyst utilization and efficiency in proton exchange membrane water electrolyzers. In this study, conductive sublayers with different conductivities and morphologies were introduced into anode electrodes in membrane-based water electrolyzers. In-situ and ex-situ investigation results showed that conductive sublayers (Au mesh and carbon nanotube (CNT) film) augmented the sheet conductivity of anode electrodes by up to 4000 times (from 2000 to 0.5 ohm square -1 ), and the ohmic resistance of water electrolyzers was reduced to 1/3 when inserting conductive sublayers. In addition, CNT film provided a higher electrochemical active area than Au mesh, because of favorable morphologies (large porosity and surface area) of CNT fibers on CNT films. Therefore, the current density of water splitting was increased by 3 times (from 4.55 mA cm -2 to 14.83 mA cm -2 ) at 2.5 V compared to a conventional anode electrode. Visualizations on bubble dynamics showed improved performances with conductive sublayers; this was mainly due to greatly increased number of reaction sites, highly spread reaction area (from 50 to 1000 um), and reduced activation overpotential. Therefore, a balance between high electronic conductivity and nanoporous morphology is essential to the anode electrode for larger reaction sites and areas in highly efficient water electrolyzers.

25 ENERGY STORAGE↗

Computational Analysis of Different Sparging Systems and their Influence in the Fluid-Dynamic Behavior of Bubble Column Reactors

Bubble column bioreactors are being actively considered for gas fermentation applications, specifically for CO2 utilization, and sugars to fuels conversion. Their main advantages include good mass transfer without any moving parts and low-cost of operation and maintenance. However, the design and scale-up of such reactors is challenging specifically for carbon capture applications where a mixture of gases (e.g. CO2/CO/H2) with variable solubilities is used. The overall performance of scaled-up bioreactors (e.g., mass transfer rate) is largely affected by gas holdup, bubble size distribution (BSD), and multiphase hydrodynamics. We investigate the effect of gas sparger designs on the performance of these large-scale bioreactors using computational fluid dynamics simulations in this work, so as to improve CO2 conversion at scale. The gas distribution systems in bubble column reactors not only determines operational regime, but also affects the evolution of the BSD, which in turn influences interfacial mass transfer and ultimately the efficiency of the gas-liquid exchange process. In addition to the BSD, uniformity in gas sparging affects gas holdup and bubble residence time which constitute important metrics of performance in gas-liquid systems. In this work, we use computational models to simulate high fidelity representations of different sparger designs and their effect on the operation of a bubble column reactor. Four different types of spargers have been selected for the computational study (Fig. 1): ladder, multi-ring, single-ring and toroidal. Their effect on superficial velocity, gas holdup mixing efficiency, and BSD will be evaluated in this work. The model uses a multiphase Eulerian framework similar to [1] and include a composition of mixtures of H2/CO/CO2 gases, common in fermentation applications.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

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

Dynamics of Gas Bubbles From a Submarine Hydrocarbon Seep Within the Hydrate Stability Zone

We validate a new model for mass transfer and bubble transport for natural seeps on the continental margins using an integrated observation of a seep at 883 m in the Northern Gulf of Mexico. In the model, mass transfer is assumed to transition from clean to dirty bubble mass transfer rates following a characteristic hydrate formation time that depends on the initial bubble surface area and the hydrate subcooling. We show that buoyancy-induced upwelling is negligible for the bubble stream. We initialize the model using precise data observed at the seafloor and validate the model results to optical and acoustic measurements within the bubble stream. The model accurately predicts the acoustic attenuation of the bubble flare by lateral spreading and bubble shrinkage throughout the water column. In conclusion, the model shows that up to 99.4% by mass of the released gases from this seep dissolves into the ocean within the hydrate stability zone.

03 NATURAL GAS↗