Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 145 records · Page 8

LDV measurements in dynamically separated flows

Two component, phase averaged mean velocity data have been obtained with an Laser Doppler Velocimeter (LDV) system, on the upper surface, near the leading edge of an oscillating airfoil undergoing compressible dynamic stall. In particular, the effect of oscillation amplitude has been studied. The results show that at an oscillation amplitude of 10 degrees, a separation bubble forms, that eventually bursts on the upstroke, well beyond the static stall angle. At 2 degrees amplitude, the bubble forms on the upstroke, but dynamic stall occurs on the downstroke. The results reveal new flow physics and the data sets serve as valuable quantitative information for validation of unsteady flow codes at transitional Reynolds numbers. The maximum velocity seen in the flow is about 1.6 times the free stream value and it occurs slightly downstream of the suction peak location. Some of the measurement difficulties are also discussed.

Chandrasekhara, M. S.↗

Bubble Formation at a Submerged Orifice in Reduced Gravity

The dynamic regime of gas injection through a circular plate orifice into an ideally wetting liquid is considered, when successively detached bubbles may be regarded as separate identities. In normal gravity and at relatively low gas flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole time of its evolution. If the flow rate exceeds a certain threshold value, another stage of the detachment process takes place in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical stem that lengthens as the bubble rises above the plate. The bubble shape resembles then that of a mushroom and the upper envelope continues to grow until the gas supply through the stem is completely cut off. Such a stage is always present under conditions of sufficiently low gravity, irrespective of the flow rate. Two major reasons make for bubble detachment: the buoyancy force and the force due to the momentum inflow into the bubble with the injected gas. The former force dominates the process at normal gravity whereas the second one plays a key role under negligible gravity conditions. It is precisely this fundamental factor that conditions the drastic influence on bubble growth and detachment that changes in gravity are able to cause. The frequency of bubble formation is proportional to and the volume of detached bubbles is independent of the gas flow rate in sufficiently low gravity, while at normal and moderately reduced gravity conditions the first variable slightly decreases and the second one almost linearly increases as the flow rate grows. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Buyevich, Yu A.↗

Effects of Gravity on Bubble Formation at a Plate Orifice

A model of the dynamic regime of gas injection through a submerged plate orifice into an ideally wetting liquid is developed in the circumstance when successively detached bubbles may be regarded as independent objects. Two major factors favor bubble detachment: the buoyancy force and a force due to the momentum inflow into the bubble with injected gas. In normal and moderately reduced gravity, the first factor dominates. At relatively low flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole period of its evolution till detachment. If the flow rate exceeds a critical value, another stage of bubble evolution occurs in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical gaseous stem that lengthens as the bubble rises above the plate. The bubble continues to grow until the gas supply through the stem is completely cut off, after which back flow of gas into the stem from the bubble becomes possible. In microgravity, the second factor prevails, and the latter stage is always present irrespective of the flow rate. However, the gas envelope rises and the stem lengthens very slowly. This difference in the underlying physical mechanisms provides for key properties of bubble growth and detachment being drastically different in appreciable and sufficiently reduced gravity. The frequency of bubble formation slightly decreases with and the detachment volume is almost proportional to the gas flow rate in the first case, in accordance with familiar relations. In the second case, the first variable is proportional to the flow rate whereas the second one is independent of it. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Webbon, Bruce W.↗

Transition effects on compressible dynamic stall of transiently pitching airfoils

Experimental results and analysis of the effects of boundary layer tripping on dynamic stall of a transiently pitching airfoil are presented. At low Mach numbers, the tripped airfoil exhibits qualitative similarity with the behavior of the untripped airfoil. However, the local supersonic flow at Mach numbers greater than 0.3 is significantly modified by the trip leading to vastly different shock/boundary layer interactions, dynamic stall onset and vortex formation angles. The formation of the laminar separation bubble is found to have a favorable influence in delaying dynamic stall on the untripped airfoil flow. In both Mach number regimes, the tripped flow actually stalls at slightly lower angles of attack. Further experimentation with three trips on an oscillating airfoil showed that the dynamic stall process is very sensitive to the state of the turbulence in the boundary layer. This sensitivity points to a need for finer turbulence modeling techniques for use in dynamic stall computations.

Wilder, M. C.↗

Interferometric investigations of compressible dynamic stall over a transiently pitching airfoil

The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Chandrasekhara, M. S.↗

Dynamical behavior of surface tension on rotating fluids in low and microgravity environments

Consideration is given to the time-dependent evolutions of the free surface profile (bubble shapes) of a cylindrical container, partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry in low and microgravity environments. The dynamics of the bubble shapes are calculated for four cases: linear time-dependent functions of spin-up and spin-down in low and microgravity, linear time-dependent functions of increasing and decreasing gravity at high and low rotating cylinder speeds, time-dependent step functions of spin-up and spin-down in low gravity, and sinusoidal function oscillation of the gravity environment in high and low rotating cylinder speeds. It is shown that the computer algorithms developed by Hung et al. (1988) may be used to simulate the profile of time-dependent bubble shapes under variations of centrifugal, capillary, and gravity forces.

Hung, R. J.↗

Hydrodynamics of bubble flow through a porous medium with applications to packed bed reactors

Gas-liquid flows through packed bed reactors (PBRs) are challenging to predict due to the tortuous flow paths that fluid interfaces must traverse. Experiments at the International Space Station showed that bubble and pulse flows are predominately observed under microgravity conditions, while the trickle and spray flows observed under terrestrial conditions are not present in microgravity. To understand the physics behind the former experiments, we simulate bubble flow through a PBR for different packing-particle-diameter-based Weber numbers and under different gravity conditions. We demonstrate different pore-scale mechanisms, such as capillary entrapment, buoyancy entrapment, and inertia-induced bubble displacement. Then, we perform a quantitative analysis by introducing new dynamic scales, dependent upon the evolving gas-liquid interfacial area, to understand the dynamic trade-offs between the inertia, capillary, and buoyancy forces on a bubble passing through a PBR. This analysis leads us to define new dimensionless Weber-like numbers that delineate bubble entrapment from bubble displacement.

Pranay P Nagrani↗

Numerical Modeling of Three-Dimensional Fluid Flow with Phase Change

We present a numerical method to compute phase change dynamics of three-dimensional deformable bubbles. The full Navier-Stokes and energy equations are solved for both phases by a front tracking/finite difference technique. The fluid boundary is explicitly tracked by discrete points that are connected by triangular elements to form a front that is used to keep the stratification of material properties sharp and to calculate the interfacial source terms. Two simulations are presented to show robustness of the method in handling complex phase boundaries. In the first case, growth of a vapor bubble in zero gravity is studied where large volume increase of the bubble is managed by adaptively increasing the front resolution. In the second case, growth of a bubble under high gravity is studied where indentation at the rear of the bubble results in a region of large curvature which challenges the front tracking in three dimensions.

Esmaeeli, Asghar↗

Analytical model of foreshock ion interaction with a discontinuity: A statistical study

When a solar wind discontinuity interacts with foreshock ions, foreshock transients such as hot flow anomalies and foreshock bubbles can form. These create significant dynamic pressure perturbations disturbing the bow shock, magnetopause, and magnetosphere-ionosphere system. However, presently these phenomena are not predictable. In the accompanying paper, we derived analytical equations of foreshock ion partial gyration around a discontinuity and the resultant current density. In this study, we utilize the derived current density strength to model the energy conversion from the foreshock ions, which drives the outward motion or expansion of the solar wind plasma away from the discontinuity. We show that the model expansion speeds match those from local hybrid simulations for varying foreshock ion parameters. Using MMS, we conduct a statistical study showing that the model expansion speeds are moderately correlated with the magnetic field strength variations and the dynamic pressure decreases around discontinuities with correlation coefficients larger than 0.5. We use conjunctions between ARTEMIS and MMS to show that the model expansion speeds are typically large for those already-formed foreshock transients. Our results show that our model can be reasonably successful in predicting significant dynamic pressure disturbances caused by foreshock ion-discontinuity interactions. We discuss ways to improve the model in the future.

Terry Z. Liu↗

Transition Analysis for the Pazy Wing

The Pazy wing, designed as a benchmark for highly flexible wings that experience large deformations, presents intriguing nonlinear phenomena due to its unique geometry with spanwise varying static deformation and low Reynolds number flow conditions. These phenomena include limit cycle oscillation related to laminar separation bubbles, boundary-layer transition, and dynamic stall. The primary focus of this study is to investigate transition behavior on the sagged Pazy wing under static loading by using the shear-stress-transport (SST)-based Langtry-Menter 𝜸 − 𝑹𝒆𝜽𝒕 transition model in the NASA OVERFLOW and FUN3D solvers. Based on the steady-state results obtained via different variants of the Langtry-Menter model, we investigate the likely transition behavior over a range of incidence angles. The predicted flowfield on the sagged Pazy wing is strongly three dimensional and the application of turbulence index as a viable indicator of the transition front in such flows is also explored. The study demonstrates the limitations of the Langtry-Menter model when applied to low Reynolds number flows. We introduce a spanwise periodic configuration to simulate a narrow-span section of the Pazy wing, specifically centered on one of the 14 nearly equidistant ribs spanning the wing. The computational results demonstrate that the spanwise periodic configuration provides a straightforward environment for studying transition characteristics as well as demonstrating systematic grid convergence and successful code-to-code comparisons. The research offers useful insights into the transition characteristics of flexible wings under various loading conditions, contributing to the field of aeroelastic simulations.

CFD Modeling↗

Transition Analysis for the Pazy Wing

The Pazy wing, designed as a benchmark for highly flexible wings that experience large deformations, presents intriguing nonlinear phenomena due to its unique geometry with spanwise varying static deformation and low Reynolds number flow conditions. These phenomena include limit cycle oscillation related to laminar separation bubbles, boundary-layer transition, and dynamic stall. The primary focus of this study is to investigate transition behavior on the sagged Pazy wing under static loading by using the shear-stress-transport (SST)-based Langtry-Menter 𝜸 − 𝑹𝒆𝜽𝒕 transition model in the NASA OVERFLOW and FUN3D solvers. Based on the steady-state results obtained via different variants of the Langtry-Menter model, we investigate the likely transition behavior over a range of incidence angles. The predicted flowfield on the sagged Pazy wing is strongly three dimensional and the application of turbulence index as a viable indicator of the transition front in such flows is also explored. The study demonstrates the limitations of the Langtry-Menter model when applied to low Reynolds number flows. We introduce a spanwise periodic configuration to simulate a narrow-span section of the Pazy wing, specifically centered on one of the 14 nearly equidistant ribs spanning the wing. The computational results demonstrate that the spanwise periodic configuration provides a straightforward environment for studying transition characteristics as well as demonstrating systematic grid convergence and successful code-to-code comparisons. The research offers useful insights into the transition characteristics of flexible wings under various loading conditions, contributing to the field of aeroelastic simulations.

CFD Modeling↗

Potential uses of vacuum bubbles in noise and vibration control

Vacuum bubbles are new acoustic elements which are dynamically more compliant than the gas volume they replace, but which are statically robust. They are made of a thin metallic shell with vacuum in their cavity. Consequently, they pose no danger in terms of contamination or fire hazard. The potential of the vacuum bubble concept for noise and vibration control was assessed with special emphases on spacecraft and aircraft applications. The following potential uses were identified: (1) as a cladding, to reduce sound radiation of vibrating surfaces and the sound excitation of structures, (2) as a screen, to reflect or absorb an incident sound wave, and (3) as a liner, to increase low frequency sound transmission loss of double walls and to increase the low frequency sound attenuation of muffler baffles. It was found that geometric and material parameters must be controlled to a very high accuracy to obtain optimal performance and that performance is highly sensitive to variations in static pressure. Consequently, it was concluded that vacuum bubbles have more potential in spacecraft applications where static pressure is controlled more than in aircraft applications where large fluctuations in static pressure are common.

Ver, Istvan L.↗

Laser velocimetry measurements of oscillating airfoil dynamic stall flow field

Ensemble-averaged two-component velocity measurements over an airfoil experiencing oscillatory dynamic stall under compressibility conditions were obtained. The measurements show the formation of a separation bubble over the airfoil that persists till angles of attack close to when the dynamic stall vortex forms and convects. The fluid attains mean velocities as large as 1.6 times the free stream velocity with instantaneous values of 1.8 times the free stream velocity. The airfoil motion induces these large velocities in regions that are far removed from the surface.

Chandrasekhara, M. S.↗

Controlling the Mobility of the Fluid Interface of Moving Gas Bubbles or Liquid Drops by Using Micellar Solutions of Surfactants

Microgravity processes must rely on mechanisms other than buoyancy to move bubbles or droplets from one region to another in a continuous liquid phase. One suggested method is thermocapillary migration in which a temperature gradient is applied to the continuous phase. A significant and as yet unresolved impediment to the use of thermocapillary migration to direct bubble or drop motion is that these migrations can be significantly retarded by the adsorption onto the fluid particle surface of surface active impurities unavoidably present in the continuous or (if the particle is a liquid) droplet phases. The focus of our research was to develop a theory for remobilizing fluid particle interfaces retarded by a surfactant impurity in an effort to make more viable the use of thermocapillary migrations for the management of bubbles and drops in microgravity. We postulated that a surfactant at high bulk concentration which kinetically exchanges rapidly with the surface can restore interface mobility. The scaling arguments along with a discussion of the previous literature is reviewed in the context of the scaling framework. The specific objectives of the research were twofold. The first was to prove the remobilization theory by studying a model problem. As the mechanism for remobilization is independent of the force which drives the particle, the fluid particle shape and the presence of fluid inertia, we chose the simplest model consisting of a spherical bubble rising steadily by buoyancy in creeping flow. We solved the hydrodynamic and surfactant transport equations for rapid kinetic exchange to demonstrate that as the concentration increases, the Marangoni retardation at first increases (the low k behavior) and then decreases (the high k behavior). The second objective was to develop a method to determine the kinetic rate constants of a surfactant molecule, since this information is necessary to select surfactants which will exchange rapidly enough relative to the convective rate in the thermocapillary process of interest. To measure the kinetic rate, we measure the dynamic tension change accompanying adsorption onto an initially clean interface, or the re-equilibration in tension when an equilibrium interface is compressed. The dynamic tension measurements are made by a pendant bubble method, in which surfactant adsorbs onto a pendant bubble, and the tension is measured by analyzing the shape change in the bubble. We conclude this report by detailing the publications, presentations and doctoral thesis completed under the auspices of this grant.

Maldarelli, Charles↗

Modeling Low-G Slosh Using Negative Mass

The modeling of the behavior of fluids in low- or microgravity is a field with little validated data but of great interest to space agencies and companies in the face of humanity’s expanding presence beyond low Earth orbit. Efforts to refuel and repair satellites in geostationary orbit require an understanding of the behavior of sloshing fuel in microgravity. The stability of remote sensing equipment and satellites depends on robust fluid slosh dampening or very accurate modeling, and the Artemis program includes the refueling of large landing systems in its architecture. However, current methodologies are quite binary: either use a relatively low-accuracy equivalent mechanical model or use computational fluid dynamics which are not practical for onboard systems due to their relatively high computational cost. In this work, a novel method of tracking the bubble rather than tracking the fuel itself is presented. Computational fluid dynamics simulations are used to validate the model based on available experimental data, and the dynamics of a refueler spacecraft based on Gateway’s Logistics Module are simulated both for the sloshing case and without sloshing considered.

Matthew M. Wittal↗

Bubble mass center and fluid feedback force fluctuations activated by constant lateral impulse with variable thrust

Sloshing dynamics within a partially filled rotating dewar of superfluid helium 2 are investigated in response to constant lateral impulse with variable thrust. The study, including how the rotating bubble of superfluid helium 2 reacts to the constant impulse with variable time period of thrust action in microgravity, how amplitudes of bubble mass center fluctuates with growth and decay of disturbances, and how fluid feedback forces fluctuates in activating on the rotating dewar through the dynamics of sloshing waves are investigated. The numerical computation of sloshing dynamics is based on the non-inertial frame spacecraft bound coordinate with lateral impulses actuating on the rotating dewar in both inertial and non-inertial frames of thrust. Results of the simulations are illustrated.

Hung, R. J.↗

Leading edge velocity field of an oscillating airfoil in compressible dynamic stall

Phase-averaged mean-velocity and turbulence data are obtained and analyzed for the leading-edge region of an oscillating airfoil under compressibility conditions. A two-component laser-Doppler velocimetry system was used to make the measurements. Results are compared for the two Mach numbers 0.3 and 0.4 at a reduced frequency of 0.05 with varying airfoil angles of attack. For a Mach number of 0.3, a separation bubble is present on the airfoil throughout the oscillation cycle and no dynamic stall occurs as the peak angle of attack is below the static stall angle. However, a slight imprint of vortical structures is seen in the shear layer enveloping the bubble at the top of the cycle, a result confirmed also by the vorticity contours and in agreement with the earlier stroboscopic schlieren studies. When the Mach number is 0.4, dynamic stall occurs with its origin in the break-up of the separation bubble. Turbulence intensities in the bubble were found to be very large.

Vandyken, R. D.↗