Harunori Yoshikawa

Dielectrophoretic force-driven thermal convection in annular geometry

The thermal convection driven by the dielectrophoretic force is investigated in annular geometry under microgravity conditions. A radial temperature gradient and a radial alternating electric field are imposed on a dielectric fluid that fills the gap of two concentric infinite-length cylinders. The resulting dielectric force is regarded as thermal buoyancy with a radial effective gravity. This electric gravity varies in space and may change its sign depending on the temperature gradient and the cylinder radius ratio. The linear stability problem is solved by a spectral-collocation method. The critical mode is stationary and non-axisymmetric. The critical Rayleigh number and wavenumbers depend sensitively on the electric gravity and the radius ratio. The mechanism behind the instability is examined from an energetic viewpoint. The instability in wide gap annuli is an exact analogue to the gravity-driven thermal instability.

Instability of the vertical annular flow with a radial heating and rotating inner cylinder

A linear stability analysis of the flow confined in a differentially rotating cylindrical annulus with a radial temperature gradient has been performed. Depending on values of control parameters (the Taylor number, the Grashof number, and the Froude number), it has shown flow destabilization to axisymmetric or non-axisymmetric modes. Analysis of different terms involved in the evolution rate of the perturbation kinetic energy has allowed us to isolate the dominant terms (centrifugal force or buoyancy force) in the destabilization process. We have shown that the centrifugal buoyancy can induce the asymmetry of the temperature gradient on critical states.

Forces on a boiling bubble in a developing boundary layer, in microgravity with g-jitter and in terrestrial conditions

Terrestrial and microgravity flow boiling experiments were carried out with the same test rig, comprising a locally heated artificial cavity in the center of a channel near the frontal edge of an intrusive glass bubble generator. Bubble shapes were in microgravity generally not far from those of truncated spheres, which permitted the computation of inertial lift and drag from potential flow theory for truncated spheres approximating the actual shape. For these bubbles, inertial lift is counteracted by drag and both forces are of the same order of magnitude as g-jitter. A generalization of the Laplace equation is found which applies to a deforming bubble attached to a plane wall and yields the pressure difference between the hydrostatic pressures in the bubble and at the wall, Δp. A fully independent way to determine the overpressure Δp is given by a second Euler-Lagrange equation. Relative differences have been found to be about 5% for both terrestrial and microgravity bubbles. A way is found to determine...

Bubble splitting in oscillatory flows on ground and in reduced gravity

The stability of centimeter scale air bubbles is studied in quiescent suspending liquid under an imposed oscillatory acceleration field. Experiments were performed in reduced- and normal-gravity environments. A strong acceleration resulted in an instability leading to the breakups of the bubbles in both gravity environments. The breakup onset was investigated and found to be characterized by a critical acceleration acr . The influence of the liquid viscosity and the gravitational environment was studied. Empirical correlations for the onset are presented and discussed with the intention to reveal splitting mechanism. The inertial mechanism often deemed to cause the breakup of drops subjected to a rapid gas stream is shown to give explanations consistent with the experiments. A breakup criterion for both gravitational environments is proposed through discussions from an energetic point of view.

Oscillatory Kelvin-Helmholtz instability. Part 2. An experiment in fluids with a large viscosity contrast

The stability of two-layer oscillatory flows was studied experimentally in a cylindrical container with a vertical axis. Two superposed immiscible liquids, differing greatly in viscosity, were set in relative oscillatory motion by alternating container rotation. Waves arising beyond a threshold were observed in detail for small oscillation frequencies ranging from 0.1 to 6 Hz. Measurements were performed on the growth rate and the wavenumber of these waves. The instability threshold was determined from the growth rate data. It was found that the threshold and the wavenumber varied with the frequency. In particular, significantly lower thresholds and longer waves were found than those predicted by the inviscid theory of the oscillatory Kelvin-Helmholtz instability. Favourable agreement with the predictions of an existing viscous theory for small oscillation amplitude flows indicates the important role of viscosity, even at the highest frequency, and suggests a similar mechanism behind the instability as that for the short wave instability in steady Couette flows. A semi-numerical stability determination for finite amplitude flows was also performed to improve the prediction in experiments with a frequency lower than 1 Hz.

Oscillatory Kelvin-Helmholtz instability. Part 1. A viscous theory

The stability of oscillatory two-layer flows is investigated with a linear perturbation analysis. An asymptotic case is considered where the oscillation amplitude is small when compared to the perturbation wavelength. The focus of the analysis is on the influence of viscosity and its contrast at the interface. The flows are unstable when the relative velocity of the layers is larger than a critical value. Depending on the oscillation frequency, the flows are in different dynamical regimes, which are characterized by the relative importance of the capillary wavelength and the thicknesses of the Stokes boundary layers developed on the interface. A particular regime is found in which instability occurs at a substantially lower critical velocity. The mechanism behind the instability is studied by identifying the velocity- and shear-induced components in the disturbance growth rate. They interchange dominance depending on the frequency and the viscosity contrast. Results of the analysis are compared with the experiments in the literature. Good agreement is found with the experiments that have a small oscillation amplitude. The validity condition of the asymptotic theory is estimated.

Dynamics of transient eddy above rolling-grain ripples

This study deals with the flow motion over the so-called rolling-grain ripples which are generated by water oscillations above a sand bed. We focus our efforts on quantifying by means of laboratory experiments and numerical calculations the morphology and the dynamics of transient flow patterns. We report, for the first time, on the formation of an unsteady pattern with closed streamlines ~we call it ‘‘eddy’’ ! above rolling-grain ripples using flow visualizations and particle image velocimetry measurements. This structure appears in the ripple trough during flow reversal and scales with the ripple wavelength. The experimental results are in qualitative agreement with the perturbative flow solution calculated by Vittori in 1989. Even if the relative ripple amplitude is not small in the experiment the perturbative expansion at the first order gives an accurate description of the flow dynamics.

Heat transfer in the thermo-electro-hydrodynamic convection under microgravity conditions

This article deals with the thermal convection in a dielectric fluid confined in a finite-length plane capacitor with a temperature gradient under microgravity conditions. The dielectrophoretic force resulting from differential polarization of the fluid plays the role of buoyancy force associated with an electric effective gravity. It induces the convection when the Rayleigh number based on this electric gravity exceeds a critical value. Two-dimensional numerical simulation for a geometry with a large aspect ratio is used to determine the convective flow in the saturated state. The Nusselt number Nu is computed for a wide range of Prandtl number (0.01 ≤ Pr ≤ 103) and its dependence on the distance from the critical condition is determined. A correlation between Nu and Pr in the vicinity of criticality is obtained and compared with that of the Rayleigh-Bénard convection. The behavior of the convection is analyzed in detail from an energetic viewpoint: electrostatic energy, power inputs by different components of the electric gravity and viscous and thermal dissipations are computed.Graphical abstract

Application of two-color LIF thermometry to nucleate boiling

The laser-induced fluorescence (LIF) thermometry is applied to measure the temperature field surrounding a single vapor bubble growing at an artificial nucleation site. In order to correct measurement errors due to the non-uniformity of the incident laser intensity, the two-color LIF thermometry technique is used in this nucleate boiling experiment. This technique is based on the use of two fluorescent dyes: the temperature sensitive dye Rhodamine B and the temperature insensitive dye Sulforhodamine-101. The concentration of the dyes is optimized by analyzing the behavior of fluorescence intensities. The mapping between the two images is determined through a geometrical calibration procedure. This technique presents a success in correcting the non uniformities due to the reflection of the light at the bubble surface and to the temperature gradient. The obtained temperature fields show that the two-color LIF is a promising technique in the investigation of nucleate boiling.

Pattern formation in bubbles emerging periodically from a liquid free surface

Abstract.Patterns formed by centimeter scale bubbles on the free surface of a viscous liquid are investigated in a cylindrical container. These bubbles emerge periodically at the surface and interact with each other in the central zone. Their radial emission, due to interaction and radial surface flow, leads to the formation of a variety of patterns. Different star-like and spiral patterns appear spontaneously by increasing the bubble emergence frequency. It is found that these patterns are due to a constant angular shift in the bubble emission direction. Measurements of this angular shift show a supercritical bifurcation accompanied by a transition from a pattern of two opposed straight arms to spiral patterns. By applying the tools and concepts from the study of leaf arrangement in botany (phyllotaxis), the recognized patterns and the mechanism of the pattern formation are discussed. Close similarities to the leaf arrangement are found in the behavior of the angular shift and the patterns. These findings suggest that the observed patterns are formed by a packing mechanism of successively appearing elements (bubbles), which is similar to that of the leaves at the earliest stage of phyllotaxis.