-
High-power laser experiment on developing supercritical shock propagating in homogeneously magnetized plasma of ambient gas origin
Authors:
S. Matsukiyo,
R. Yamazaki,
T. Morita,
K. Tomita,
Y. Kuramitsu,
S. J. Tanaka,
T. Takezaki,
S. Isayama,
T. Higuchi,
H. Murakami,
Y. Horie,
N. Katsuki,
R. Hatsuyama,
M. Edamoto,
H. Nishioka,
M. Takagi,
T. Kojima,
S. Tomita,
N. Ishizaka,
S. Kakuchi,
S. Sei,
K. Sugiyama,
K. Aihara,
S. Kambayashi,
M. Ota
, et al. (8 additional authors not shown)
Abstract:
A developing supercritical collisionless shock propagating in a homogeneously magnetized plasma of ambient gas origin having higher uniformity than the previous experiments is formed by using high-power laser experiment. The ambient plasma is not contaminated by the plasma produced in the early time after the laser shot. While the observed developing shock does not have stationary downstream struc…
▽ More
A developing supercritical collisionless shock propagating in a homogeneously magnetized plasma of ambient gas origin having higher uniformity than the previous experiments is formed by using high-power laser experiment. The ambient plasma is not contaminated by the plasma produced in the early time after the laser shot. While the observed developing shock does not have stationary downstream structure, it possesses some characteristics of a magnetized supercritical shock, which are supported by a one-dimensional full particle-in-cell simulation taking the effect of finite time of laser-target interaction into account.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
High-power laser experiment forming a supercritical collisionless shock in a magnetized uniform plasma at rest
Authors:
Ryo Yamazaki,
S. Matsukiyo,
T. Morita,
S. J. Tanaka,
T. Umeda,
K. Aihara,
M. Edamoto,
S. Egashira,
R. Hatsuyama,
T. Higuchi,
T. Hihara,
Y. Horie,
M. Hoshino,
A. Ishii,
N. Ishizaka,
Y. Itadani,
T. Izumi,
S. Kambayashi,
S. Kakuchi,
N. Katsuki,
R. Kawamura,
Y. Kawamura,
S. Kisaka,
T. Kojima,
A. Konuma
, et al. (29 additional authors not shown)
Abstract:
We present a new experimental method to generate quasi-perpendicular supercritical magnetized collisionless shocks. In our experiment, ambient nitrogen (N) plasma is at rest and well-magnetized, and it has uniform mass density. The plasma is pushed by laser-driven ablation aluminum (Al) plasma. Streaked optical pyrometry and spatially resolved laser collective Thomson scattering clarify structures…
▽ More
We present a new experimental method to generate quasi-perpendicular supercritical magnetized collisionless shocks. In our experiment, ambient nitrogen (N) plasma is at rest and well-magnetized, and it has uniform mass density. The plasma is pushed by laser-driven ablation aluminum (Al) plasma. Streaked optical pyrometry and spatially resolved laser collective Thomson scattering clarify structures of plasma density and temperatures, which are compared with one-dimensional particle-in-cell simulations. It is indicated that just after the laser irradiation, the Al plasma is magnetized by a self-generated Biermann battery field, and the plasma slaps the incident N plasma. The compressed external field in the N plasma reflects N ions, leading to counter-streaming magnetized N flows. Namely we identify the edge of the reflected N ions. Such interacting plasmas form a magnetized collisionless shock.
△ Less
Submitted 7 February, 2022; v1 submitted 19 January, 2022;
originally announced January 2022.
-
On viscid-inviscid interactions of a pair of bubbles rising near the wall
Authors:
Kazuki Maeda,
Masanobu Date,
Kazuyasu Sugiyama,
Shu Takagi,
Yoichiro Matsumoto
Abstract:
Series of experiments on turbulent bubbly channel flows observed bubble clusters near the wall which can change large-scale flow structures. To gain insights into clustering mechanisms, we study the interaction of a pair of spherical bubbles rising in a vertical channel through combined experiments and modeling. Experimental imaging identifies that pairwise bubbles of 1.0 mm diameter take two pref…
▽ More
Series of experiments on turbulent bubbly channel flows observed bubble clusters near the wall which can change large-scale flow structures. To gain insights into clustering mechanisms, we study the interaction of a pair of spherical bubbles rising in a vertical channel through combined experiments and modeling. Experimental imaging identifies that pairwise bubbles of 1.0 mm diameter take two preferred configurations depending on their mutual distance: side-by-side positions for a short distance ($S<5$) and nearly inline, oblique positions for a long distance ($S>5$), where $S$ is the mutual distance normalized by the bubble radius. In the model, we formulate the motions of pairwise bubbles rising at $Re=O(100)$. Analytical drag and lift, and semi-empirical, spatio-temporal stochastic forcing are employed to represent the mean acceleration and the fluctuation due to turbulent agitation, respectively. The model is validated against the experiment through comparing Lagrangian statistics of the bubbles. Simulations using this model identify two distinct timescales of interaction dynamics which elucidate the preferred configurations. For pairs initially in-line, the trailing bubble rapidly escapes from the viscous wake of the leading bubble to take the oblique position. Outside of the wake, the trailing bubble travels on a curve-line path with a slower velocity driven by potential interaction and horizontally approaches the leading bubble to become side-by-side. Moreover, statistical analysis identifies that the combination of the wake and the agitation can significantly accelerate the side-by-side clustering of in-line pairs. These results indicate positive contributions of liquid viscosity and turbulence to the formation of bubble clusters.
△ Less
Submitted 16 October, 2020;
originally announced October 2020.
-
Direct Phase Locking of a Ti:Sapphire Optical Frequency Comb to a Remote Optical Frequency Standard
Authors:
Eunmi Chae,
Kota Nakashima,
Takuya Ikeda,
Kei Sugiyama,
Kosuke Yoshioka
Abstract:
We report on an ultralow noise optical frequency transfer from a remotely located Sr optical lattice clock laser to a Ti:Sapphire optical frequency comb through telecom-wavelength optical fiber networks. The inherent narrow linewidth of the Ti:Sapphire optical frequency comb eliminates the need for a local reference high-finesse cavity. The relative fractional frequency instability of the optical…
▽ More
We report on an ultralow noise optical frequency transfer from a remotely located Sr optical lattice clock laser to a Ti:Sapphire optical frequency comb through telecom-wavelength optical fiber networks. The inherent narrow linewidth of the Ti:Sapphire optical frequency comb eliminates the need for a local reference high-finesse cavity. The relative fractional frequency instability of the optical frequency comb with respect to the remote optical reference was $6.7(1) \times 10^{-18}$ at 1 s and $1.05(3) \times 10^{-19}$ at 1,000 s including a 2.9 km-long fiber network. This ensured the optical frequency comb had the same precision as the optical standard. Our result paves the way for ultrahigh-precision spectroscopy and conversion of the highly precise optical frequency to radio frequencies in a simpler setup.
△ Less
Submitted 2 April, 2019; v1 submitted 19 July, 2018;
originally announced July 2018.
-
Equatorial jet in the lower to middle cloud layer of Venus revealed by Akatsuki
Authors:
T. Horinouchi,
S. Murakami,
T. Satoh,
J. Peralta,
K. Ogohara,
T. Kouyama,
T. Imamura,
H. Kashimura,
S. S. Limaye,
K. McGouldrick,
M. Nakamura,
T. M. Sato,
K. Sugiyama,
M. Takagi,
S. Watanabe,
M. Yamada,
A. Yamazaki,
E. F. Young
Abstract:
The Venusian atmosphere is in a state of superrotation where prevailing westward winds move much faster than the planet's rotation. Venus is covered with thick clouds that extend from about 45 to 70 km altitude, but thermal radiation emitted from the lower atmosphere and the surface on the planet's nightside escapes to space at narrow spectral windows of the near-infrared. The radiation can be use…
▽ More
The Venusian atmosphere is in a state of superrotation where prevailing westward winds move much faster than the planet's rotation. Venus is covered with thick clouds that extend from about 45 to 70 km altitude, but thermal radiation emitted from the lower atmosphere and the surface on the planet's nightside escapes to space at narrow spectral windows of the near-infrared. The radiation can be used to estimate winds by tracking the silhouettes of clouds in the lower and middle cloud regions below about 57 km in altitude. Estimates of wind speeds have ranged from 50 to 70 m/s at low to mid-latitudes, either nearly constant across latitudes or with winds peaking at mid-latitudes. Here we report the detection of winds at low latitude exceeding 80 m/s using IR2 camera images from the Akatsuki orbiter taken during July and August 2016. The angular speed around the planetary rotation axis peaks near the equator, which we suggest is consistent with an equatorial jet, a feature that has not been observed previously in the Venusian atmosphere. The mechanism producing the jet remains unclear. Our observations reveal variability in the zonal flow in the lower and middle cloud region that may provide clues to the dynamics of Venus's atmospheric superrotation.
△ Less
Submitted 7 September, 2017;
originally announced September 2017.
-
Lattice Boltzmann Simulations of Droplet formation in confined Channels with Thermocapillary flows
Authors:
A. Gupta,
M. Sbragaglia,
D. Belardinelli,
K. Sugiyama
Abstract:
Based on mesoscale lattice Boltzmann simulations with the "Shan-Chen" model, we explore the influence of thermocapillarity on the break-up properties of fluid threads in a microfluidic T-junction, where a dispersed phase is injected perpendicularly into a main channel containing a continuous phase, and the latter induces periodic break-up of droplets due to the cross-flowing. Temperature effects a…
▽ More
Based on mesoscale lattice Boltzmann simulations with the "Shan-Chen" model, we explore the influence of thermocapillarity on the break-up properties of fluid threads in a microfluidic T-junction, where a dispersed phase is injected perpendicularly into a main channel containing a continuous phase, and the latter induces periodic break-up of droplets due to the cross-flowing. Temperature effects are investigated by switching on/off both positive/negative temperature gradients along the main channel direction, thus promoting a different thread dynamics with anticipated/delayed break-up. Numerical simulations are performed at changing the flow-rates of both the continuous and dispersed phases, as well as the relative importance of viscous forces, surface tension forces and thermocapillary stresses. The range of parameters is broad enough to characterize the effects of thermocapillarity on different mechanisms of break-up in the confined T-junction, including the so-called "squeezing" and "dripping" regimes, previously identified in the literature. Some simple scaling arguments are proposed to rationalize the observed behaviour, and to provide quantitative guidelines on how to predict the droplet size after break-up.
△ Less
Submitted 14 November, 2016; v1 submitted 19 July, 2016;
originally announced July 2016.
-
Flow states in two-dimensional Rayleigh-Bénard convection as a function of aspect-ratio and Rayleigh number
Authors:
Erwin P. van der Poel,
Richard J. A. M. Stevens,
Kazuyasu Sugiyama,
Detlef Lohse
Abstract:
In this numerical study on two-dimensional Rayleigh-Bénard convection we consider $10^7 \leq Ra \leq 10^{12}$ in aspect ratio $0.23 \leq Γ\leq 13$ samples. We focus on several cases. First we consider small aspect ratio cells, where at high Ra number we find a sharp transition from a low Ra number branch towards a high Ra number branch, due to changes in the flow structure. Subsequently, we show t…
▽ More
In this numerical study on two-dimensional Rayleigh-Bénard convection we consider $10^7 \leq Ra \leq 10^{12}$ in aspect ratio $0.23 \leq Γ\leq 13$ samples. We focus on several cases. First we consider small aspect ratio cells, where at high Ra number we find a sharp transition from a low Ra number branch towards a high Ra number branch, due to changes in the flow structure. Subsequently, we show that the influence of the aspect ratio on the heat transport decreases with increasing aspect ratio, although even at very large aspect ratio of $Γ\approx10$ variations up to 2.5% in the heat transport as a function of Γare observed. Finally, we observe long-lived transients up to at least $Ra=10^9$, as in certain aspect ratio cells we observe different flow states that are stable for thousands of turnover times.
△ Less
Submitted 3 July, 2012; v1 submitted 18 June, 2012;
originally announced June 2012.
-
Horizontal Structures of Velocity and Temperature Boundary Layers in 2D Numerical Turbulent Rayleigh-Bénard Convection
Authors:
Quan Zhou,
Kazuyasu Sugiyama,
Richard J. A. M. Stevens,
Siegfried Grossmann,
Detlef Lohse,
Ke-Qing Xia
Abstract:
We investigate the structures of the near-plate velocity and temperature profiles at different horizontal positions along the conducting bottom (and top) plate of a Rayleigh-Bénard convection cell, using two-dimensional (2D) numerical data obtained at the Rayleigh number Ra=10^8 and the Prandtl number Pr=4.4 of an Oberbeck-Boussinesq flow with constant material parameters. The results show that mo…
▽ More
We investigate the structures of the near-plate velocity and temperature profiles at different horizontal positions along the conducting bottom (and top) plate of a Rayleigh-Bénard convection cell, using two-dimensional (2D) numerical data obtained at the Rayleigh number Ra=10^8 and the Prandtl number Pr=4.4 of an Oberbeck-Boussinesq flow with constant material parameters. The results show that most of the time, and for both velocity and temperature, the instantaneous profiles scaled by the dynamical frame method [Q. Zhou and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010) agree well with the classical Prandtl-Blasius laminar boundary layer (BL) profiles. Therefore, when averaging in the dynamical reference frames, which fluctuate with the respective instantaneous kinematic and thermal BL thicknesses, the obtained mean velocity and temperature profiles are also of Prandtl-Blasius type for nearly all horizontal positions. We further show that in certain situations the traditional definitions based on the time-averaged profiles can lead to unphysical BL thicknesses, while the dynamical method also in such cases can provide a well-defined BL thickness for both the kinematic and the thermal BLs.
△ Less
Submitted 30 September, 2011;
originally announced September 2011.
-
Flow reversals in thermally driven turbulence
Authors:
Kazuyasu Sugiyama,
Rui Ni,
Richard J. A. M. Stevens,
Tak Shing Chan,
Sheng-Qi Zhou,
Heng-Dong Xi,
Chao Sun,
Siegfried Grossmann,
Ke-Qing Xia,
Detlef Lohse
Abstract:
We analyze the reversals of the large scale flow in Rayleigh-Bénard convection both through particle image velocimetry flow visualization and direct numerical simulations (DNS) of the underlying Boussinesq equations in a (quasi) two-dimensional, rectangular geometry of aspect ratio 1. For medium Prandtl number there is a diagonal large scale convection roll and two smaller secondary rolls in the t…
▽ More
We analyze the reversals of the large scale flow in Rayleigh-Bénard convection both through particle image velocimetry flow visualization and direct numerical simulations (DNS) of the underlying Boussinesq equations in a (quasi) two-dimensional, rectangular geometry of aspect ratio 1. For medium Prandtl number there is a diagonal large scale convection roll and two smaller secondary rolls in the two remaining corners diagonally opposing each other. These corner flow rolls play a crucial role for the large scale wind reversal: They grow in kinetic energy and thus also in size thanks to plume detachments from the boundary layers up to the time that they take over the main, large scale diagonal flow, thus leading to reversal. Based on this mechanism we identify a typical time scale for the reversals. We map out the Rayleigh number vs Prandtl number phase space and find that the occurrence of reversals very sensitively depends on these parameters.
△ Less
Submitted 30 September, 2011;
originally announced September 2011.
-
Electromagnetic response of a metamaterial with field-gradient-induced transparency
Authors:
Yasuhiro Tamayama,
Toshihiro Nakanishi,
Yasuhiro Wakasa,
Tetsuo Kanazawa,
Kazuhiko Sugiyama,
Masao Kitano
Abstract:
We investigate a dynamically controllable electromagnetically induced transparency-like metamaterial. The unit structure of the metamaterial consists of a low-quality-factor resonator and a high-quality-factor resonator. The field gradient of the incident electromagnetic wave in the transverse direction induces coupling between these two types of resonators and causes a transparency phenomenon. We…
▽ More
We investigate a dynamically controllable electromagnetically induced transparency-like metamaterial. The unit structure of the metamaterial consists of a low-quality-factor resonator and a high-quality-factor resonator. The field gradient of the incident electromagnetic wave in the transverse direction induces coupling between these two types of resonators and causes a transparency phenomenon. We present the simulation and experimental results for the dynamic control of the electromagnetic response.
△ Less
Submitted 30 July, 2011;
originally announced August 2011.
-
Numerical simulations of compressible Rayleigh-Taylor turbulence in stratified fluids
Authors:
A. Scagliarini,
L. Biferale,
M. Sbragaglia,
K. Sugiyama,
F. Toschi
Abstract:
We present results from numerical simulations of Rayleigh-Taylor turbulence, performed using a recently proposed lattice Boltzmann method able to describe consistently a thermal compressible flow subject to an external forcing. The method allowed us to study the system both in the nearly-Boussinesq and strongly compressible regimes. Moreover, we show that when the stratification is important, the…
▽ More
We present results from numerical simulations of Rayleigh-Taylor turbulence, performed using a recently proposed lattice Boltzmann method able to describe consistently a thermal compressible flow subject to an external forcing. The method allowed us to study the system both in the nearly-Boussinesq and strongly compressible regimes. Moreover, we show that when the stratification is important, the presence of the adiabatic gradient causes the arrest of the mixing process.
△ Less
Submitted 29 September, 2010;
originally announced September 2010.
-
A full Eulerian finite difference approach for solving fluid-structure coupling problems
Authors:
Kazuyasu Sugiyama,
Satoshi Ii,
Shintaro Takeuchi,
Shu Takagi,
Yoichiro Matsumoto
Abstract:
A new simulation method for solving fluid-structure coupling problems has been developed. All the basic equations are numerically solved on a fixed Cartesian grid using a finite difference scheme. A volume-of-fluid formulation (Hirt and Nichols (1981, J. Comput. Phys., 39, 201)), which has been widely used for multiphase flow simulations, is applied to describing the multi-component geometry. The…
▽ More
A new simulation method for solving fluid-structure coupling problems has been developed. All the basic equations are numerically solved on a fixed Cartesian grid using a finite difference scheme. A volume-of-fluid formulation (Hirt and Nichols (1981, J. Comput. Phys., 39, 201)), which has been widely used for multiphase flow simulations, is applied to describing the multi-component geometry. The temporal change in the solid deformation is described in the Eulerian frame by updating a left Cauchy-Green deformation tensor, which is used to express constitutive equations for nonlinear Mooney-Rivlin materials. In this paper, various verifications and validations of the present full Eulerian method, which solves the fluid and solid motions on a fixed grid, are demonstrated, and the numerical accuracy involved in the fluid-structure coupling problems is examined.
△ Less
Submitted 19 September, 2010;
originally announced September 2010.
-
On the lateral migration of a slightly deformed bubble rising near a vertical plane wall
Authors:
Kazuyasu Sugiyama,
Fumio Takemura
Abstract:
Deformation-induced lateral migration of a bubble slowly rising near a vertical plane wall in a stagnant liquid is numerically and theoretically investigated. In particular, our focus is set on a situation with a short clearance $c$ between the bubble interface and the wall. Motivated by the fact that numerically and experimentally measured migration velocities are considerably higher than the vel…
▽ More
Deformation-induced lateral migration of a bubble slowly rising near a vertical plane wall in a stagnant liquid is numerically and theoretically investigated. In particular, our focus is set on a situation with a short clearance $c$ between the bubble interface and the wall. Motivated by the fact that numerically and experimentally measured migration velocities are considerably higher than the velocity estimated by the available analytical solution using the Faxén mirror image technique for $a/(a+c)\ll 1$ (here $a$ is the bubble radius), when the clearance parameter $\varepsilon(= c/a)$ is comparable to or smaller than unity, the numerical analysis based on the boundary-fitted finite-difference approach solving the Stokes equation is performed to complement the experiment. The migration velocity is found to be more affected by the high-order deformation modes with decreasing $\varepsilon$. The numerical simulations are compared with a theoretical migration velocity obtained from a lubrication study of a nearly spherical drop, which describes the role of the squeezing flow within the bubble-wall gap. The numerical and lubrication analyses consistently demonstrate that when $\varepsilon\leq 1$, the lubrication effect makes the migration velocity asymptotically $μV_{B1}^2/(25\varepsilon γ)$ (here, $V_{B1}$, $μ$, and $γ$ denote the rising velocity, the dynamic viscosity of liquid, and the surface tension, respectively).
△ Less
Submitted 11 July, 2010;
originally announced July 2010.
-
Prandtl-Blasius temperature and velocity boundary layer profiles in turbulent Rayleigh-Bénard convection
Authors:
Quan Zhou,
Richard J. A. M. Stevens,
Kazuyasu Sugiyama,
Siegfried Grossmann,
Detlef Lohse,
Ke-Qing Xia
Abstract:
The shape of velocity and temperature profiles near the horizontal conducting plates in turbulent Rayleigh-Bénard convection are studied numerically and experimentally over the Rayleigh number range $10^8\lesssim Ra\lesssim3\times10^{11}$ and the Prandtl number range $0.7\lesssim Pr\lesssim5.4$. The results show that both the temperature and velocity profiles well agree with the classical Prandt…
▽ More
The shape of velocity and temperature profiles near the horizontal conducting plates in turbulent Rayleigh-Bénard convection are studied numerically and experimentally over the Rayleigh number range $10^8\lesssim Ra\lesssim3\times10^{11}$ and the Prandtl number range $0.7\lesssim Pr\lesssim5.4$. The results show that both the temperature and velocity profiles well agree with the classical Prandtl-Blasius laminar boundary-layer profiles, if they are re-sampled in the respective dynamical reference frames that fluctuate with the instantaneous thermal and velocity boundary-layer thicknesses.
△ Less
Submitted 5 February, 2010;
originally announced February 2010.
-
Flow organization in non-Oberbeck-Boussinesq Rayleigh-Benard convection in water
Authors:
Kazuyasu Sugiyama,
Enrico Calzavarini,
Siegfried Grossmann,
Detlef Lohse
Abstract:
Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; the…
▽ More
Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; these together describe the various features of the rather complex flow organization. The results are presented both as functions of the Rayleigh number Ra (with Ra up to 10^8) for fixed temperature difference (Delta) between top and bottom plates and as functions of Delta ("non-Oberbeck-Boussinesqness") for fixed Ra with Delta up to 60 K. All results are consistent with the available experimental NOB data for the center temperature Tc and the Nusselt number ratio Nu_{NOB}/Nu_{OB} (the label OB meaning that the Oberbeck-Boussinesq conditions are valid).
Beyond Ra ~ 10^6 the flow consists of a large diagonal center convection roll and two smaller rolls in the upper and lower corners. In the NOB case the center convection roll is still characterized by only one velocity scale.
△ Less
Submitted 20 December, 2008;
originally announced December 2008.
-
Wetting failure and contact line dynamics in a Couette flow
Authors:
M. Sbragaglia,
K. Sugiyama,
L. Biferale
Abstract:
Liquid-liquid wetting failure is investigated in a two-dimensional Couette system with two immiscible fluids of arbitrary viscosity. The problem is solved exactly using a sharp interface treatment of hydrodynamics (lubrication theory) as a function of the capillary number, viscous ratio and separation of scale, i.e. slip length versus macroscopic scale of the system. The existence of critical ve…
▽ More
Liquid-liquid wetting failure is investigated in a two-dimensional Couette system with two immiscible fluids of arbitrary viscosity. The problem is solved exactly using a sharp interface treatment of hydrodynamics (lubrication theory) as a function of the capillary number, viscous ratio and separation of scale, i.e. slip length versus macroscopic scale of the system. The existence of critical velocities, above which no stationary solutions are found, is analyzed in detail in terms of the relevant parameters of the system. Comparisons with existing analysis for other geometries are also carried out. A numerical method of analysis is also presented, based on diffuse interface models obtained from multiphase extensions of the lattice Boltzmann equation (LBE). Sharp interface and diffuse interface models are quantitatively compared face to face indicating the correct limit of applicability of the diffuse interface models.
△ Less
Submitted 17 December, 2008;
originally announced December 2008.
-
Soft-lubrication effect on the lateral migration of a slightly deformed bubble rising near a vertical plane wall
Authors:
Kazuyasu Sugiyama,
Fumio Takemura
Abstract:
Deformation-induced lateral migration of a bubble slowly rising near a vertical plane wall in a stagnant liquid is numerically and theoretically investigated. In particular, our focus is set on a situation with a small clearance $c$ between the bubble interface and the wall. Motivated by the fact that experimentally measured migration velocity (Takemura et al. (2002, J. Fluid Mech. {\bf 461}, 27…
▽ More
Deformation-induced lateral migration of a bubble slowly rising near a vertical plane wall in a stagnant liquid is numerically and theoretically investigated. In particular, our focus is set on a situation with a small clearance $c$ between the bubble interface and the wall. Motivated by the fact that experimentally measured migration velocity (Takemura et al. (2002, J. Fluid Mech. {\bf 461}, 277)) is higher than the velocity estimated by the available analytical solution (Magnaudet et al. (2003, J. Fluid Mech. {\bf 476}, 115)) using the Faxén mirror image technique for $κ(=a/(a+c))\ll 1$ (here $a$ is the bubble radius), when the clearance parameter $ε(=c/a)$ is comparable to or smaller than unit, the numerical analysis based on the boundary-fitted finite-difference approach by solving the Stokes equation is performed to complement the experiment. To improve the understandings of a role of the squeezing flow within the bubble-wall gap, the theoretical analysis based on a soft-lubrication approach (Skotheim & Mahadevan (2004, Phys. Rev. Lett. {\bf 92}, 245509)) is also performed. The present analyses demonstrate the migration velocity scales $\propto{\rm Ca}\ ε^{-1}V_{B1}$ (here, $V_{B1}$ and ${\rm Ca}$ denote the rising velocity and the capillary number, respectively) in the limit of $ε\to 0$.
△ Less
Submitted 19 September, 2010; v1 submitted 21 November, 2008;
originally announced November 2008.
-
An invisible medium for circularly polarized electromagnetic waves
Authors:
Y. Tamayama,
T. Nakanishi,
K. Sugiyama,
M. Kitano
Abstract:
We study the no reflection condition for a planar boundary between vacuum and an isotropic chiral medium. In general chiral media, elliptically polarized waves incident at a particular angle satisfy the no reflection condition. When the wave impedance and wavenumber of the chiral medium are equal to the corresponding parameters of vacuum, one of the circularly polarized waves is transmitted to t…
▽ More
We study the no reflection condition for a planar boundary between vacuum and an isotropic chiral medium. In general chiral media, elliptically polarized waves incident at a particular angle satisfy the no reflection condition. When the wave impedance and wavenumber of the chiral medium are equal to the corresponding parameters of vacuum, one of the circularly polarized waves is transmitted to the medium without reflection or refraction for all angles of incidence. We propose a circular polarizing beam splitter as a simple application of the no reflection effect.
△ Less
Submitted 8 December, 2008; v1 submitted 4 November, 2008;
originally announced November 2008.
-
Non-Oberbeck-Boussinesq effects in turbulent thermal convection in ethane close to the critical point
Authors:
Guenter Ahlers,
Enrico Calzavarini,
Francisco Fontenele Araujo,
Denis Funfschilling,
Siegfried Grossmann,
Detlef Lohse,
Kazuyasu Sugiyama
Abstract:
As shown in earlier work (Ahlers et al., J. Fluid Mech. 569, p.409 (2006)), non-Oberbeck Boussinesq (NOB) corrections to the center temperature in turbulent Rayleigh-Benard convection in water and also in glycerol are governed by the temperature dependences of the kinematic viscosity and the thermal diffusion coefficient. If the working fluid is ethane close to the critical point the origin of n…
▽ More
As shown in earlier work (Ahlers et al., J. Fluid Mech. 569, p.409 (2006)), non-Oberbeck Boussinesq (NOB) corrections to the center temperature in turbulent Rayleigh-Benard convection in water and also in glycerol are governed by the temperature dependences of the kinematic viscosity and the thermal diffusion coefficient. If the working fluid is ethane close to the critical point the origin of non-Oberbeck-Boussinesq corrections is very different, as will be shown in the present paper. Namely, the main origin of NOB corrections then lies in the strong temperature dependence of the isobaric thermal expansion coefficient β(T). More precisely, it is the nonlinear T-dependence of the density ρ(T) in the buoyancy force which causes another type of NOB effect. We demonstrate that through a combination of experimental, numerical, and theoretical work, the latter in the framework of the extended Prandtl-Blasius boundary layer theory developed in Ahlers et al., J. Fluid Mech. 569, p.409 (2006). The latter comes to its limits, if the temperature dependence of the thermal expension coefficient β(T) is significant.
△ Less
Submitted 22 February, 2008; v1 submitted 2 December, 2007;
originally announced December 2007.
-
Non-Oberbeck-Boussinesq effects in two-dimensional Rayleigh-Benard convection in glycerol
Authors:
Kazuyasu Sugiyama,
Enrico Calzavarini,
Siegfried Grossmann,
Detlef Lohse
Abstract:
We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of "non-Oberbeck-Boussinesqness'' or "NOBness''…
▽ More
We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of "non-Oberbeck-Boussinesqness'' or "NOBness'' ($Δ$) up to 50 K (for fixed Ra). For this large NOBness the center temperature $T_c$ is more than 5 K larger than the arithmetic mean temperature $T_m$ between top and bottom plate and only weakly depends on Ra. To physically account for the NOB deviations of the Nusselt numbers from its Oberbeck-Boussinesq values, we apply the decomposition of $Nu_{NOB}/Nu_{OB}$ into the product of two effects, namely first the change in the sum of the top and bottom thermal BL thicknesses, and second the shift of the center temperature $T_c$ as compared to $T_m$. While for water the origin of the $Nu$ deviation is totally dominated by the second effect (cf. Ahlers et al., J. Fluid Mech. 569, pp. 409 (2006)) for glycerol the first effect is dominating, in spite of the large increase of $T_c$ as compared to $T_m$.
△ Less
Submitted 2 October, 2007; v1 submitted 1 July, 2007;
originally announced July 2007.
-
Microbubbly drag reduction in Taylor-Couette flow in the wavy vortex regime
Authors:
Kazuyasu Sugiyama,
Enrico Calzavarini,
Detlef Lohse
Abstract:
We investigate the effect of microbubbles on Taylor-Couette flow by means of direct numerical simulations. We employ an Eulerian-Lagrangian approach with a gas-fluid coupling based on the point-force approximation. Added mass, drag, lift, and gravity are taken into account in the modeling of the motion of the individual bubble. We find that very dilute suspensions of small non-deformable bubbles…
▽ More
We investigate the effect of microbubbles on Taylor-Couette flow by means of direct numerical simulations. We employ an Eulerian-Lagrangian approach with a gas-fluid coupling based on the point-force approximation. Added mass, drag, lift, and gravity are taken into account in the modeling of the motion of the individual bubble. We find that very dilute suspensions of small non-deformable bubbles (volume void fraction below 1%, zero Weber number and bubble Reynolds number <10) induce a robust statistically steady drag reduction (up to 20%) in the so called wavy vortex flow regime (Re = 600-2500). The Reynolds number dependence of the normalized torque (the so-called Torque Reduction Ratio (TRR) which corresponds to the drag reduction) is consistent with a recent series of experimental measurements performed by Murai et al. (J. Phys. 14, 143 (2005)). Our analysis suggests that the physical mechanism for the torque reduction in this regime is due to the local axial forcing, induced by rising bubbles, that is able to break the highly dissipative Taylor wavy vortices in the system. We finally show that the lift force acting on the bubble is crucial in this process. When neglecting it, the bubbles preferentially accumulate near the inner cylinder and the bulk flow is less efficiently modified.
△ Less
Submitted 18 July, 2007; v1 submitted 8 March, 2007;
originally announced March 2007.
-
Linear shear flow past a hemispherical droplet adhering to a solid surface
Authors:
K. Sugiyama,
M. Sbragaglia
Abstract:
This paper investigates the properties of a three dimensional shear flow overpassing a hemispherical droplet resting on a plane wall. The exact solution is computed as a function of the viscosity ratio between the droplet and the surrounding fluid and generalizes the solution for the hemispherical no-slip bump given in an earlier paper by Price (1985). Several expressions including the torque an…
▽ More
This paper investigates the properties of a three dimensional shear flow overpassing a hemispherical droplet resting on a plane wall. The exact solution is computed as a function of the viscosity ratio between the droplet and the surrounding fluid and generalizes the solution for the hemispherical no-slip bump given in an earlier paper by Price (1985). Several expressions including the torque and the force acting on the drop will be considered as well as the importance of the deformations on the surface for small Capillary numbers.
△ Less
Submitted 20 August, 2007; v1 submitted 19 March, 2007;
originally announced March 2007.
-
Generalized Lattice Boltzmann Method with multi-range pseudo-potential
Authors:
M. Sbragaglia,
R. Benzi,
L. Biferale,
S. Succi,
K. Sugiyama,
F. Toschi
Abstract:
The physical behaviour of a class of mesoscopic models for multiphase flows is analyzed in details near interfaces. In particular, an extended pseudo-potential method is developed, which permits to tune the equation of state and surface tension independently of each other. The spurious velocity contributions of this extended model are shown to vanish in the limit of high grid refinement and/or h…
▽ More
The physical behaviour of a class of mesoscopic models for multiphase flows is analyzed in details near interfaces. In particular, an extended pseudo-potential method is developed, which permits to tune the equation of state and surface tension independently of each other. The spurious velocity contributions of this extended model are shown to vanish in the limit of high grid refinement and/or high order isotropy. Higher order schemes to implement self-consistent forcings are rigorously computed for 2d and 3d models. The extended scenario developed in this work clarifies the theoretical foundations of the Shan-Chen methodology for the lattice Boltzmann method and enhances its applicability and flexibility to the simulation of multiphase flows to density ratios up to O(100).
△ Less
Submitted 7 November, 2006;
originally announced November 2006.
-
Boundary induced non linearities at small Reynolds Numbers
Authors:
M. Sbragaglia,
K. Sugiyama
Abstract:
We investigate the influence of boundary slip velocity in Newtonian fluids at finite Reynolds numbers. Numerical simulations with Lattice Boltzmann method (LBM) and Finite Differences method (FDM) are performed to quantify the effect of heterogeneous boundary conditions on the integral and local properties of the flow. Non linear effects are induced by the non homogeneity of the boundary conditi…
▽ More
We investigate the influence of boundary slip velocity in Newtonian fluids at finite Reynolds numbers. Numerical simulations with Lattice Boltzmann method (LBM) and Finite Differences method (FDM) are performed to quantify the effect of heterogeneous boundary conditions on the integral and local properties of the flow. Non linear effects are induced by the non homogeneity of the boundary condition and change the symmetry properties of the flow inducing an overall mean flow reduction. To explain the observed drag modification, reciprocal relations for stationary ensembles are used, predicting a reduction of the mean flow rate from the creeping flow to be proportional to the fourth power of the friction Reynolds number. Both numerical schemes are then validated within the theoretical predictions and reveal a pronounced numerical efficiency of the LBM with respect to FDM.
△ Less
Submitted 20 March, 2006;
originally announced March 2006.
-
Observation of Brewster's effect for transverse-electric electromagnetic waves in metamaterials: Experiment and theory
Authors:
Y. Tamayama,
T. Nakanishi,
K. Sugiyama,
M. Kitano
Abstract:
We have experimentally realized Brewster's effect for transverse-electric waves with metamaterials. In dielectric media, Brewster's no-reflection effect arises only for transverse-magnetic waves. However, it has been theoretically predicted that Brewster's effect arises for TE waves under the condition that the relative permeability r is not equal to unity. We have designed an array of split-rin…
▽ More
We have experimentally realized Brewster's effect for transverse-electric waves with metamaterials. In dielectric media, Brewster's no-reflection effect arises only for transverse-magnetic waves. However, it has been theoretically predicted that Brewster's effect arises for TE waves under the condition that the relative permeability r is not equal to unity. We have designed an array of split-ring resonators as a metamaterial with mu_r 1 using a finite-difference time-domain method. The reflection measurements were carried out in a 3-GHz region and the disappearance of reflected waves at a particular incident angle was confirmed.
△ Less
Submitted 19 May, 2006; v1 submitted 31 October, 2005;
originally announced November 2005.