• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter iResearch App Facebook

Search Issue | RSS Feeds RSS
Next Issue

Jan 2011

Volume 23, Issue 1, Articles (01xxxx)

Issue Cover Spotlight Figure

Phys. Fluids 23, 011702 (2011); http://dx.doi.org/10.1063/1.3541844 (4 pages)

Hongjie Zhong, Shiyi Chen, and Cunbiao Lee
back to top
RSS Feeds
back to top Interfacial Flows

Numerical simulation of a liquid bridge in a coaxial gas flow

Miguel A. Herrada, José M. López-Herrera, Emilio J. Vega, and José M. Montanero

Phys. Fluids 23, 012101 (2011); http://dx.doi.org/10.1063/1.3534076 (11 pages) | Cited 7 times

Online Publication Date: 11 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The dynamical response of an isothermal liquid bridge to a coaxial gas stream is examined from axisymmetric numerical simulations of the Navier–Stokes equations. The simulation method is previously validated by calculating the temporal evolution of the first oscillation mode in both cylindrical and axisymmetric liquid bridges. The comparison with other theoretical approaches and experiments shows good agreement in most cases, although significant discrepancies are found between the simulation and the experimental values of the damping rate for hexadecane. The simulation of a liquid bridge in a coaxial gas stream shows that a recirculation cell always appears in the liquid driven by the gas viscous stress on the free surface. The recirculation cell speed depends quasilinearly on the gas velocity for the range of gas flow rates considered. If the gas stream and gravity have the same direction, then the speed of the recirculation cell increases considerably due to the free surface deformation of the liquid bridge at equilibrium. This effect does not occur when gravity has the opposite direction because viscous dissipation in the liquid increases in this case. If the gas stream and gravity point downward, the liquid bridge shrinks at the upper part and bulges at the lower owing to the accumulation of momentum there. The same occurs for zero gravity, but noncylindrical liquid bridges deform more than cylindrical shapes with the same slenderness. If one inverts the direction of the gravity force, the interface deformation caused by the gas stream is the opposite, and its magnitude is smaller. The magnitude of the free surface deformation depends almost linearly on the gas stream velocity for both zero and normal gravity conditions.
Show PACS
47.10.ad Navier-Stokes equations
47.55.nb Capillary and thermocapillary flows
02.60.Lj Ordinary and partial differential equations; boundary value problems

Parametric stability and dynamic buckling of an encapsulated microbubble subject to acoustic disturbances

Kostas Tsiglifis and Nikos A. Pelekasis

Phys. Fluids 23, 012102 (2011); http://dx.doi.org/10.1063/1.3536646 (28 pages) | Cited 4 times

Online Publication Date: 11 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Stability analysis of the radial pulsations of a gas microbubble that is encapsulated by a thin viscoelastic shell and surrounded by an ideal incompressible liquid is carried out. Small axisymmetric disturbances in the microbubble shape are imposed and their long and short term stability is examined depending on the initial bubble radius, the shell properties, and the parameters, i.e., frequency and amplitude, of the external acoustic excitation. Owing to the anisotropy of the membrane that is forming the encapsulating shell, two different types of elastic energy are accounted for, namely, the membrane and bending energy per unit of initial area. They are used to describe the tensions that develop on the shell due to shell stretching and bending, respectively. In addition, two different constitutive laws are used in order to relate the tensions that develop on the membrane as a result of stretching, i.e., the Mooney–Rivlin law describing materials that soften as deformation increases and the Skalak law describing materials that harden as deformation increases. The limit for static buckling is obtained when the external overpressure exerted upon the membrane surpasses a critical value that depends on the membrane bending resistance. The stability equations describing the evolution of axisymmetric disturbances, in the presence of an external acoustic field, reveal that static buckling becomes relevant when the forcing frequency is much smaller than the resonance frequency of the microbubble, corresponding to the case of slow compression. The resonance frequencies for shape oscillations of the microbubble are also obtained as a function of the shell parameters. Floquet analysis shows that parametric instability, similar to the case of an oscillating free bubble, is possible for the case of a pulsating encapsulated microbubble leading to shape oscillations as a result of subharmonic or harmonic resonance. These effects take place for acoustic amplitude values that lie above a certain threshold but below those required for static buckling to occur. They are quite useful in providing estimates for the shell elasticity and bending resistance based on a frequency/amplitude sweep that monitors the onset of shape oscillations when the forcing frequency resonates with the radial pulsation, ωf = ω0, or with a certain shape mode, ωf = 2ωn. An acceleration based instability, identified herein as dynamic buckling, is observed during the compression phase of the pulsation, evolving over a small number of periods of the forcing, when the amplitude of the acoustic excitation is further increased. It corresponds to the Rayleigh–Taylor instability observed for free bubbles, and has been observed with contrast agents as well, e.g., BR-14. Finally, phase diagrams for contrast agent BR-14 are constructed and juxtaposed with available experimental data, illustrating the relevance and range of the above instabilities.
Show PACS
47.55.dd Bubble dynamics
47.20.Ma Interfacial instabilities (e.g., Rayleigh-Taylor)
47.50.Gj Instabilities
83.60.Df Nonlinear viscoelasticity
83.50.-v Deformation and flow

Nonlinear development of oscillatory instability in a three-layer system under the joint action of buoyancy and thermocapillary effect

Ilya B. Simanovskii, Antonio Viviani, Frank Dubois, and Jean-Claude Legros

Phys. Fluids 23, 012103 (2011); http://dx.doi.org/10.1063/1.3536655 (12 pages) | Cited 4 times

Online Publication Date: 11 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The nonlinear development of oscillatory instability under the joint action of buoyant and thermocapillary effects in multilayer system is investigated. The nonlinear convective regimes are studied by the finite difference method. The calculations have been performed for two-dimensional flows. The interfaces are assumed to be nondeforming. Rigid heat-insulated lateral walls are considered. Transitions between the flows with different spatial structures are studied. Specific types of nonlinear flows—symmetric and asymmetric oscillations—have been found. It is shown that the oscillatory flow takes place in an interval of Grashof number values bounded both from below by the quiescent mechanical equilibrium, and from above by a convecting steady state. Cavities with different lengths are considered.
Show PACS
47.20.Bp Buoyancy-driven instabilities (e.g., Rayleigh-Benard)
47.55.nb Capillary and thermocapillary flows
47.11.Bc Finite difference methods
02.70.Bf Finite-difference methods
47.55.Hd Stratified flows

Droplet charging regimes for ultrasonic atomization of a liquid electrolyte in an external electric field

Thomas P. Forbes, F. Levent Degertekin, and Andrei G. Fedorov

Phys. Fluids 23, 012104 (2011); http://dx.doi.org/10.1063/1.3541818 (10 pages)

Online Publication Date: 11 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Distinct regimes of droplet charging, determined by the dominant charge transport process, are identified for an ultrasonic droplet ejector using electrohydrodynamic computational simulations, a fundamental scale analysis, and experimental measurements. The regimes of droplet charging are determined by the relative magnitudes of the dimensionless Strouhal and electric Reynolds numbers, which are a function of the process (pressure forcing), advection, and charge relaxation time scales for charge transport. Optimal (net maximum) droplet charging has been identified to exist for conditions in which the electric Reynolds number is of the order of the inverse Strouhal number, i.e., the charge relaxation time is on the order of the pressure forcing (droplet formation) time scale. The conditions necessary for optimal droplet charging have been identified as a function of the dimensionless Debye number (i.e., liquid conductivity), external electric field (magnitude and duration), and atomization drive signal (frequency and amplitude). The specific regime of droplet charging also determines the functional relationship between droplet charge and charging electric field strength. The commonly expected linear relationship between droplet charge and external electric field strength is only found when either the inverse of the Strouhal number is less than the electric Reynolds number, i.e., the charge relaxation is slower than both the advection and external pressure forcing, or in the electrostatic limit, i.e., when charge relaxation is much faster than all other processes. The analysis provides a basic understanding of the dominant physics of droplet charging with implications to many important applications, such as electrospray mass spectrometry, ink jet printing, and drop-on-demand manufacturing.
Show PACS
47.55.db Drop and bubble formation
47.65.-d Magnetohydrodynamics and electrohydrodynamics
47.11.-j Computational methods in fluid dynamics
47.80.-v Instrumentation and measurement methods in fluid dynamics

Analysis of time-dependent nonlinear dynamics of the axisymmetric liquid film on a vertical circular cylinder: Energy integral model

E. Novbari and A. Oron

Phys. Fluids 23, 012105 (2011); http://dx.doi.org/10.1063/1.3541856 (12 pages) | Cited 2 times

Online Publication Date: 14 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The nonlinear dynamics of an axisymmetric liquid film falling on the outer surface of a vertical cylinder is investigated in the framework of the set of two coupled evolution equations derived recently using the energy integral method (EIM). This set of EIM evolution equations is solved numerically and its solutions are compared with the traveling wave solutions derived from it using AUTO. We find that traveling wave solutions of EIM equations can bifurcate either supercritically or subcritically from the base state. The type of bifurcation depends on the parameter set of the problem. The set of EIM equations studied here admits both traveling wave and nonstationary wave flows. We demonstrate that in the case of subcritical primary bifurcation the film dynamics is sensitive to the choice of the initial condition and coexistence of up to five different flows is possible for the same parameter set in the domain of a given periodicity. The case of supercritical primary bifurcation exhibits much lesser dependence on the initial condition, though coexistence of two different flows for the same parameter set is possible. The synergetic approach based on both direct numerical solution of the governing evolution equations and search of traveling wave solutions using AUTO facilitate a discovery of a large variety of flows and help to conclude about stability of the traveling wave flows found using AUTO.
Show PACS
68.15.+e Liquid thin films
47.35.-i Hydrodynamic waves
47.20.Ky Nonlinearity, bifurcation, and symmetry breaking
FREE

Wall energy relaxation in the Cahn–Hilliard model for moving contact lines

Pengtao Yue and James J. Feng

Phys. Fluids 23, 012106 (2011); http://dx.doi.org/10.1063/1.3541806 (8 pages) | Cited 13 times

Online Publication Date: 25 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The Cahn–Hilliard model uses diffusion between fluid components to regularize the stress singularity at a moving contact line. In addition, it represents the dynamics of the near-wall layer by the relaxation of a wall energy. The first part of the paper elucidates the role of the wall relaxation in a flowing system, with two main results. First, we show that wall energy relaxation produces a dynamic contact angle that deviates from the static one, and derive an analytical formula for the deviation. Second, we demonstrate that wall relaxation competes with Cahn–Hilliard diffusion in defining the apparent contact angle, the former tending to “rotate” the interface at the contact line while the latter to “bend” it in the bulk. Thus, varying the two in coordination may compensate each other to produce the same macroscopic solution that is insensitive to the microscopic dynamics of the contact line. The second part of the paper exploits this competition to develop a computational strategy for simulating realistic flows with microscopic slip length at a reduced cost. This consists in computing a moving contact line with a diffusion length larger than the real slip length, but using the wall relaxation to correct the solution to that corresponding to the small slip length. We derive an analytical criterion for the required amount of wall relaxation, and validate it by numerical results on dynamic wetting in capillary tubes and drop spreading.
Show PACS
47.60.Dx Flows in ducts and channels
47.32.Ef Rotating and swirling flows
47.45.Gx Slip flows and accommodation

Dynamics of nearly unstable axisymmetric liquid bridges

José M. Perales and José M. Vega

Phys. Fluids 23, 012107 (2011); http://dx.doi.org/10.1063/1.3541814 (11 pages) | Cited 1 time

Online Publication Date: 25 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The dynamics of a noncylindrical, axisymmetric, marginally unstable liquid bridge between two equal disks is analyzed in the inviscid limit. The resulting model allows for the weakly nonlinear description of both the (first stage of) breakage for unstable configurations and the (slow) dynamics for stable configurations. The analysis is made for both slender and short liquid brides. In the former range, the dynamics breaks reflection symmetry on the midplane between the supporting disks and can be described by a standard Duffing equation, while for short bridges reflection symmetry is preserved and the equation is still Duffing-like but exhibiting a quadratic nonlinearity. The asymptotic results compare well with existing experiments.
Show PACS
47.20.Cq Inviscid instability
47.60.Dx Flows in ducts and channels

The Rayleigh–Taylor instability of a surface of arbitrary cross section with pinned edges

L. E. Johns and R. Narayanan

Phys. Fluids 23, 012108 (2011); http://dx.doi.org/10.1063/1.3541819 (3 pages) | Cited 1 time

Online Publication Date: 25 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We determine the critical points of the Rayleigh–Taylor instability of a surface of arbitrary cross section having pinned edges. Often these points coincide with the diffusion eigenvalues but sometimes they do not.
Show PACS
47.20.Ma Interfacial instabilities (e.g., Rayleigh-Taylor)
02.10.Ud Linear algebra

First drop dissimilarity in drop-on-demand inkjet devices

Amin Famili, Saurabh A. Palkar, and William J. Baldy, Jr.

Phys. Fluids 23, 012109 (2011); http://dx.doi.org/10.1063/1.3543758 (6 pages) | Cited 2 times

Online Publication Date: 25 January 2011

Full Text: Read Online (HTML) | Download PDF

Show Abstract
As inkjet printing technology is increasingly applied in a broader array of applications, careful characterization of its method of use is critical due to its inherent sensitivity. A common operational mode in inkjet technology known as drop-on-demand ejection is used as a way to deliver a controlled quantity of material to a precise location on a target. This method of operation allows ejection of individual or a sequence (burst) of drops based on a timed trigger event. This work presents an examination of sequences of drops as they are ejected, indicating a number of phenomena that must be considered when designing a drop-on-demand inkjet system. These phenomena appear to be driven by differences between the first ejected drop in a burst and those that follow it and result in a break-down of the linear relationship expected between driving amplitude and drop mass. This first drop, as quantified by high-speed videography and subsequent image analysis, can be different in morphology, trajectory, velocity, and volume from subsequent drops within a burst. These findings were confirmed orthogonally by both volume and mass measurement techniques which allowed quantitation down to single drops.
Show PACS
47.55.df Breakup and coalescence
47.80.Jk Flow visualization and imaging
06.30.Dr Mass and density
Close
Google Calendar
ADVERTISEMENT

close