Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Acoustic waves superimposed on incompressible flows

The use of incompressible approximations in deriving solutions to the Lighthill wave equation was investigated for problems where an analytical solution could be found. A particular model problem involves the determination of the sound field of a spherical oscillating bubble in an ideal fluid. It is found that use of incompressible boundary conditions leads to good approximations in the important region of high acoustic wave number.

Hodge, Steve↗
NASA NTRS2001

Effect of Marangoni Convection Generated by Voids on Segregation During Low-G and 1-G Solidification

Solidification experiments, especially microgravity solidification experiments, are often compromised by the evolution of unwanted voids or bubbles in the melt. Although these voids and/or bubbles are highly undesirable, there is currently no effective means of preventing their formation or of eliminating their adverse effects, particularly during microgravity experiments. Marangoni convection caused by these voids can drastically change the transport processes in the melt. Recent microgravity experiments by Matthiesen (1) Andrews (2) and Fripp (3) are perfect examples of how voids and bubbles can affect the outcome of costly space experiments and significantly increase the level of difficulty in interpreting their results. Formation of bubbles have caused problems in microgravity experiments for a long time. Even in the early Skylab mission an unexpectedly large number of bubbles were detected in the four materials processing experiments reported by Papazian and Wilcox (4). They demonstrated that while during ground-based tests bubbles were seen to detach from the interface easily and float to the top of the melt, in low-gravity tests no detachment from the interface occurred and large voids were grown in the crystal. More recently, the lead-tin-telluride crystal growth experiment of Fripp et al.(3) flown aboard the USMP-3 mission has provided very interesting results. The purpose of the study was to investigate the effect of natural convection on the solidification process by growing the samples at different orientations with respect to the gravitational field. Large pores and voids were found in the three solid crystal samples processed in space. Post-growth characterization of the compositional profiles of the cells indicated considerable levels of mixing even in the sample grown in the hot-on-top stable configuration. The mixing was attributed to thermocapillary convection caused by the voids and bubbles which evolved during growth. Since the thermocapillary convection is orientation-independent, diffusion-controlled growth was not possible in any of the samples, even the top-heated one. These results are consistent with recent studies of thermocapillary convection generated by a bubble on a heated surface undertaken by Kassemi and Rashidnia (5-7) where it is numerically and experimentally shown that the thermocapillary flow generated by a bubble in a model fluid (silicone oil) can drastically modify the temperature field through vigorous mixing of the fluid around it, especially under microgravity conditions.

Kassemi, M.↗
NASA NTRS2008

Embossed Teflon AF Laminate Membrane Microfluidic Diaphragm Valves

A microfluidic system has been designed to survive spaceflight and to function autonomously on the Martian surface. It manipulates microscopic quantities of liquid water and performs chemical analyses on these samples to assay for the presence of molecules associated with past or present living processes. This technology lies at the core of the Urey Instrument, which is scheduled for inclusion on the Pasteur Payload of the ESA ExoMars rover mission in 2013. Fabrication processes have been developed to make the microfabricated Teflon-AF microfluidic diaphragm pumps capable of surviving extreme temperature excursions before and after exposure to liquid water. Two glass wafers are etched with features and a continuous Teflon membrane is sandwiched between them (see figure). Single valves are constructed using this geometry. The microfabricated devices are then post processed by heating the assembled device while applying pneumatic pressure to force the Teflon diaphragm against the valve seat while it is softened. After cooling the device, the embossed membrane retains this new shape. This solves previous problems with bubble introduction into the fluid flow where deformations of the membrane at the valve seat occurred during device bonding at elevated temperatures (100-150 C). The use of laminated membranes containing commercial Teflon AF 2400 sheet sandwiched between spun Teflon AF 1600 layers performed best, and were less gas permeable than Teflon AF 1600 membranes on their own. Spinning Teflon AF 1600 solution (6 percent in FLOURINERT(Registered TradeMark) FC40 solvent, 3M Company) at 500 rpm for 1.5 seconds, followed by 1,000 rpm for 3 seconds onto Borofloat glass wafers, results in a 10-micron-thick film of extremely smooth Teflon AF. This spinning process is repeated several times on flat, blank, glass wafers in order to gradually build a thick, smooth membrane. After running this process at least five times, the wafer and Teflon coating are heated under vacuum at 220 C for one hour in order to drive off any residual solvent present in the composite film. After this, a second blank, glass wafer is brought down from above and the stack is held under vacuum at 3 atm mechanical pressure for ten 10 hours.

Willis, Peter↗
NASA NTRS2008

Optical Measurement of Mass Flow of a Two-Phase Fluid

An optoelectronic system utilizes wavelength-dependent scattering of light for measuring the density and mass flow of a two-phase fluid in a pipe. The apparatus was invented for original use in measuring the mass flow of a two-phase cryogenic fluid (e.g., liquid hydrogen containing bubbles of hydrogen gas), but underlying principles of operation can readily be adapted to non-cryogenic two-phase fluids. The system (see figure) includes a laser module, which contains two or more laser diodes, each operating at a different wavelength. The laser module also contains beam splitters that combine the beams at the various wavelengths so as to produce two output beams, each containing all of the wavelengths. One of the multiwavelength output beams is sent, via a multimode fiberoptic cable, to a transmitting optical coupler. The other multiwavelength output beam is sent, via another multimode fiber-optic cable, to a reference detector module, wherein fiber-optic splitters split the light into several multiwavelength beams, each going to a photodiode having a spectral response that is known and that differs from the spectral responses of the other photodiodes. The outputs of these photodiodes are digitized and fed to a processor, which executes an algorithm that utilizes the known spectral responses to convert the photodiode outputs to obtain reference laser-power levels for the various wavelengths. The transmitting optical coupler is mounted in (and sealed to) a hole in the pipe and is oriented at a slant with respect to the axis of the pipe. The transmitting optical coupler contains a collimating lens and a cylindrical lens that form the light emerging from the end of the fiber-optic cable into a fan-shaped beam in a meridional plane of the pipe. Receiving optical couplers similar to the transmitting optical couplers are mounted in the same meridional plane at various longitudinal positions on the opposite side of the pipe, approximately facing the transmitting optical coupler along the same slant. Light collected by each receiving optical coupler is sent, via a multimode fiber-optic cable, to a detector module similar to the reference detector module. The outputs of the photodiodes in each detector module are digitized and processed, similarly to those of the reference detector module, to obtain indications of the amounts of light of each wavelength scattered to the corresponding receiving position. The value for each wavelength at each position is also normalized to the reference laser-power level for that wavelength. From these normalized values, the density and the mass flow rate of the fluid are estimated.

Wiley, John↗
NASA NTRS2005

A New Unsteady Model for Dense Cloud Cavitation in Cryogenic Fluids

A new unsteady, cavitation model is presented wherein the phase change process (bubble growth/collapse) is coupled to the acoustic field in a cryogenic fluid. It predicts the number density and radius of bubbles in vapor clouds by tracking both the aggregate surface area and volume fraction of the cloud. Hence, formulations for the dynamics of individual bubbles (e.g. Rayleigh-Plesset equation) may be integrated within the macroscopic context of a dense vapor cloud i.e. a cloud that occupies a significant fraction of available volume and contains numerous bubbles. This formulation has been implemented within the CRUNCH CFD, which has a compressible real fluid formulation, a multi-element, unstructured grid framework, and has been validated extensively for liquid rocket turbopump inducers. Detailed unsteady simulations of a cavitating ogive in liquid nitrogen are presented where time-averaged mean cavity pressure and temperature depressions due to cavitation are compared with experimental data. The model also provides the spatial and temporal history of the bubble size distribution in the vapor clouds that are shed, an important physical parameter that is difficult to measure experimentally and is a significant advancement in the modeling of dense cloud cavitation.

Hosangadi, A.↗
DOE OSTI2020

DEMONSTRATION OF A DATA-DRIVEN PHYSICS-BASED APPROACH FOR COMPUTATIONALLY EFFICIENT CFD PREDICTION OF TWO-PHASE BUBBLY FLOW

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale physics-guided data-driven approach, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. By defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

42 ENGINEERING↗
NASA NTRS2020

Formation and Topology of Foreshock Bubbles

We use global and local hybrid (kinetic ions and fluid electrons) simulations to investigate the conditions under which foreshock bubbles (FBs) form and how their topology changes with solar wind conditions. FBs form as a result of the interaction between solar wind discontinuities and backstreaming ion beams in the foreshock. They consist of an outer shock and its associated sheath plasma and a low density high temperature core with low magnetic field strength. The structure of FBs is determined by the angle between the interplanetary magnetic field and the normal to the solar wind discontinuity. We show that interaction of rotational discontinuities with the foreshock during small angles between the interplanetary magnetic field and discontinuity normal results in the formation of a nearly spherical bubble with a radius that scales with the width of the foreshock. As this angle increases, FBs become more elongated and eventually become nearly planar structures with dimensions that scale with the length of the foreshock. Despite this transformation, the signatures of FBs in spacecraft time series data remain the same in agreement with the observations. Global simulation results show that FBs form when the solar wind flow speed corresponds to high or intermediate Alfvén Mach numbers (approximately >7 MA). In general, this is tied to the relative speed between the solar wind and ion beams and drop in density of the back streaming ions.

Nick Omidi↗
NASA NTRS2000

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗
NASA NTRS1996

Effects of Gravity and Shear on the Dynamics and Stability of Particulate and Multiphase Flows

The main objectives of this project are to understand the differing particulate and multiphase flow behaviors that will occur in space and in Earth's gravity. More specifically, the project is concerned with understanding the effect of shear and gravity on two relatively ideal suspensions with significant inertial effects. The first is a gas-solid suspension at small Reynolds numbers and finite Stokes numbers. In this type of suspensions the inertia of the particle phase is significant while the hydrodynamic interactions are dominated by viscous forces in the suspending fluid. The other is a bubble suspension at small Weber and large Reynolds numbers. The hydrodynamic interactions in such suspensions are dominated by the inertial effects in the suspending fluid, but these inertial interactions can be described using potential flow theory. Our main objective is to examine the effects of shear and gravity on the average properties and stability of these two suspensions.

Sangani, Ashor S.↗
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Interpretation of high-speed flows in the plasma sheet

Pursuing an idea suggested by Pontius and Wolf (1990), we propose that the `bursty bulk flows' observed by Baumjohann et al. (1990) and Angelopoulos et al. (1992) are `bubbles' in the Earth's plasma sheet. Specifically, they are flux tubes that have lower values of pV(exp 5/3) than their neighbors, where p is the thermal pressure of the particles and V is the volume of a tube containing one unit of magnetic flux. Whether they are created by reconnection or some other mechanism, the bubbles are propelled earthward by a magnetic buoyancy force, which is related to the interchange instability. Most of the major observed characteristics of the bursty bulk flows can be interpreted naturally in terms of the bubble picture. We propose a new `stratified fluid' picture of the plasma sheet, based on the idea that bubbles constitute the crucial transport mechanism. Results from simple mathematical models of plasma sheet transport support the idea that bubbles can resolve the pressure balance inconsistency, particularly in cases where plasma sheet ions are lost by gradient/curvature drift out the sides of the tail or bubbles are generated by reconnection in the middle of plasma sheet.

Chen, C. X.↗
NASA NTRSDate not supplied

Ground-based Characterization of Plant Water Management (PWM) Hydroponic Root Modules for Spaceflight

Hydroponic crop production in space is crucial for long-term space travel but faces numerous challenges – one of which is providing sufficient dissolved oxygen (DO) in nutrient solution. In microgravity environments, surface tension is the primary force acting on liquids, causing water to form into suspended spherical droplets. This can suffocate plants as the liquid clings onto plant roots and the lack of aeration deprives the plant of oxygen needed for growth. To overcome these challenges, plant water management (PWM) systems explore options of growing plants in space autonomously and passively through capillary forces. From 2018-2023, there have been 6 PWM experiments conducted on the ISS. In past experiments conducted in microgravity, bubbles formed in test cells have disrupted fluid dynamics and may adversely affect plant growth. The accumulation of bubbles may lead to inconsistent nutrient delivery, break prime in tubing, and suspend plant roots in air, leading to plant stress and eventually death. Current research efforts focus on the oxygenation capabilities of the PWM system along with a comprehensive sensor array that will improve nutrient and DO monitoring capabilities.

L Wang↗
NASA NTRS1999

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗
NASA NTRS2014

A Predictive Bubble Point Pressure Model for Porous Liquid Acquisition Device Screens

This article presents a simplified model for porous screen channel liquid acquisition devices based on a maximum bubble point pressure method from Adamson and Gast (1997). To validate the model, three 304 stainless steel (325 × 2300, 450 × 2750, and 510 × 3600) mesh samples were tested in methanol, acetone, isopropyl alcohol, and water. Screen pores are estimated based on analysis from scanning electron microscopy, historical data, and current test data. Results show that the bubble point pressure is proportional to the surface tension of the fluid only when accounting for nonzero contact angles. The previous assumption that bubble point pressure scales inversely with effective pore diameter is shown to be invalid, as the second finest 450 × 2750 produced the highest bubble point of the three screens. The simplified bubble point model can be used to make predictions for any pure fluid when pore diameters are based on bubble point tests and not SEM analysis.

Liquid Acquisition Device↗
NASA NTRS1998

Investigation of a Method to Reduce Cavitation in Diesel Engine Bearings

Sonoluminescence is the effect of producing light from sound and occurs when a gas bubble is trapped in a fluid filled cavity and is forced to collapse under a barrage of sound waves. Frenzel and Schultes discovered this phenomenon in 1934 while exposing acoustic waves to photographic plates. This effect was not well understood until 1988 when Crum and Gaitan discovered the necessary conditions for producing single bubble sonoluminescence in the laboratory. The luminescence is a result of the bubble violently collapsing from sound waves and this shares a close association with vibratory cavitation. Cavitation erosion is known to cause damage to rotational machinery when the collapse is near to surfaces due to the high pressures associated with bubble collapse. With these high pressures and temperatures there is a considerable amount of damage to the outside layer of a bearing, thereby, reducing its useful life. An experiment was constructed to generate sonoluminescence in the laboratory in order to obtain a greater understanding of this phenomenon and its association with bubble cavitation. Most of the research was done to investigate how to obtain single bubble sonoluminescence under different conditions and to determine how to detect it. Success in this has inspired several theories on how to use the methods for generating sonoluminescence to control cavitation in fluids under industrial conditions.

Keith, Theo G., Jr.↗
NASA NTRS1989

Analytic theory for the determination of velocity and stability of bubbles in a Hele-Shaw cell. I - Velocity selection. II - Stability

An asymptotic theory is presented for the determination of velocity and linear stability of a steady symmetric bubble in a Hele-Shaw cell for small surface tension. First the bubble velocity relative to the fluid velocity at infinity is determined for small surface tension by means of a transcendentally small correction to the asymptotic series solution. In addition, a linear stability analysis shows that only the solution branch corresponding to the largest possible bubble velocity for given surface tension is stable, while all the others are unstable.

Tanveer, S.↗