Monte Carlo Simulation Of Thermal Fluctuations Below The Onset Of Rayleigh-Benard Convection

The density fluctuations below the onset of convection in the Rayleigh-Benard problem are studied with the direct simulation Monte Carlo method. The particle simulation results clearly show the connection between the static correlation functions of fluctuations below the critical Rayleigh number and the flow patterns above the onset of convection for small Knudsen number flows (Kn=0.01 and Kn=0.005). Furthermore, the physical nature for no convection in the Rayleigh-Benard problem under large Knudsen number conditions (Kn>0.028) is explained based on the dynamics of fluctuations.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

Study On Buoyancy Convection Phenomenon In The Crystal Growth Process

Real-time phase shift Mach-Zehnder interference technique, imaging technique, and computer image processing technique were combined to perform a real-time diagnosis of NaClO3 crystal, which described both the dissolution process and the crystallization process of the NaClO3 crystal in real-time condition. The dissolution fringes and the growth fringes in the process were obtained. Moreover, a distribution of concentration field in this process was obtained by inversion calculation. Finally, the buoyancy convection phenomenon caused by gravity in the crystal growth process was analyzed. The results showed that this convection phenomenon directly influences the growth rate of each crystal face in the crystal.

Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge

In present study, effect of interfacial heat transfer with ambient gas on the onset of oscillatory convection in a liquid bridge of large Prandtl number on the ground is systematically investigated by the method of linear stability analyses. With both the constant and linear ambient air temperature distributions, the numerical results show that the interfacial heat transfer modifies the free-surface temperature distribution directly and then induces a steeper temperature gradient on the middle part of the free surface, which may destabilize the convection. On the other hand, the interfacial heat transfer restrains the temperature disturbances on the free surface, which may stabilize the convection. The two coupling effects result in a complex dependence of the stability property on the Biot number. Effects of melt free-surface deformation on the critical conditions of the oscillatory convection were also investigated. Moreover, to better understand the mechanism of the instabilities, rates of kinetic energy change and "thermal" energy change of the critical disturbances were investigated (C) 2009 Elsevier Ltd. All rights reserved.

Experimental Study on the Combined Evaporation Effect and Thermocapillary Convection in a Thin Liquid Layer

The coupling mechanisms and flow characteristics of thermocapillary convection in a thin liquid layer with evaporating interface were studied. The planar liquid layer, with the upper surface open to air, was imposed externally horizontal temperature differences. The measured average evaporating rates and interfacial temperature profiles indicated the relative importance of evaporation effect and thermocapillary convection under different temperature gradients. A temperature jump was found at the interface, which was thought to be related to the influence of evaporation effect. All above mentioned results were repeated in a rarely evaporating liquid to compare the influence of evaporation effect.

Thermocapillary migrations of drops and bubbles

The transport phenomenon of drops or bubbles is a very important topic in fundamental hydrodynamics research and practical applications such as material processing and the chemical engineering. In microgravity environment, if drops or bubbles stay in a continuous phase with non-uniform temperature ¯eld, they will start to move as a result of the variance of the interface tension. This kind of movement is called the Marangoni migration. This review tries to sum up the main results in this ¯eld on theoretical analysis, numerical simulations and experiments. So far the theoretical analysis is still limited to the linear or weak nonlinear steady questions, while the current numerical simulations can already obtain the time- dependent process of the bubble/drop migration when the e®ect of heat convection is small. For strong heat convection problem, or when the Marangoni number is bigger than 100, no numerical result is in consistence with those of experiments so far. Some of the lastest numerical results are shown when heat convection is strong, and the main di®erence between strong and weak heat convection is analyzed. Finally, we also discuss the main unresolved problems in this ¯eld and some possible directions in the future.

Effects of convection and solid wall on the diffusion in microscale convection flows

The diffusive transport properties in microscale convection flows are studied by using the direct simulation Monte Carlo method. The effective diffusion coefficient D is computed from the mean square displacements of simulated molecules based on the Einstein diffusion equation D = x2 t /2t. Two typical convection flows, namely, thermal creep convection and Rayleigh– Bénard convection, are investigated. The thermal creep convection in our simulation is in the noncontinuum regime, with the characteristic scale of the vortex varying from 1 to 100 molecular mean free paths. The diffusion is shown to be enhanced only when the vortex scale exceeds a certain critical value, while the diffusion is reduced when the vortex scale is less than the critical value. The reason for phenomenon of diffusion reduction in the noncontinuum regime is that the reduction effect due to solid wall is dominant while the enhancement effect due to convection is negligible. A molecule will lose its memory of macroscopic velocity when it collides with the walls, and thus molecules are hard to diffuse away if they are confined between very close walls. The Rayleigh– Bénard convection in our simulation is in the continuum regime, with the characteristic length of 1000 molecular mean free paths. Under such condition, the effect of solid wall on diffusion is negligible. The diffusion enhancement due to convection is shown to scale as the square root of the Péclet number in the steady convection regime, which is in agreement with previous theoretical and experimental results. In the oscillation convection regime, the diffusion is more strongly enhanced because the molecules can easily advect from one roll to its neighbor due to an oscillation mechanism. © 2010 American Institute of Physics. doi:10.1063/1.3528310

Thermocapillary Migration of Deformable Bubbles at Moderate to Large Marangoni Number in Microgravity

Using the level-set method and the continuum interface model, the axisymmetric thermocapillary migration of gas bubbles in an immiscible bulk liquid with a temperature gradient at moderate to large Marangoni number is simulated numerically. Constant material properties of the two phases are assumed. Steady state of the motion can always be reached. The terminal migration velocity decreases monotonously with the increase of the Marangoni number due to the wrapping of isotherms around the front surface of the bubble. Good agreements with space experimental data and previous theoretical and numerical studies in the literature are evident. Slight deformation of bubble is observed, but no distinct influence on the motion occurs. It is also found that the influence of the convective transport of heat inside bubbles cannot be neglected at finite Marangoni number, while the influence of the convective transport of momentum inside bubbles may be actually negligible.

Period-doubling bifurcations of the thermocapillary convection in a floating half zone

This study experimentally explored the fine structures of the successive period-doubling bifurcations of the time-dependent thermocapillary convection in a floating half zone of 10 cSt silicone oil with the diameter d (0)=3.00 mm and the aspect ratio A=l/d (0)=0.72 in terrestrial conditions. The onset of time-dependent thermocapillary convection predominated in this experimental configuration and its subsequent evolution were experimentally detected through the local temperature measurements. The experimental results revealed a sequence of period-doubling bifurcations of the time-dependent thermocapillary convection, similar in some way to one of the routes to chaos for buoyant natural convection. The critical frequencies and the corresponding fractal frequencies were extracted through the real-time analysis of the frequency spectra by Fast-Fourier-Transformation (FFT). The projections of the trajectory onto the reconstructed phase-space were also provided. Furthermore, the experimentally predicted Feigenbaum constants were quite close to the theoretical asymptotic value of 4.669 [Feigenbaum M J. Phys Lett A, 1979, 74: 375-378].

Effect of volume ratio on thermocapillary flow in liquid bridges of high-Prandtl-number fluids

In present study, the transition of thermocapillary convection from the axisymmetric stationary flow to oscillatory flow in liquid bridges of 5cst silicon oil (aspect ratio 1.0 and 1.6) is investigated in microgravity conditions by the linear instability analysis. The corresponding marginal instability boundary is closely related to the gas/liquid configuration of the liquid bridge noted as volume ratio. With the increasing volume ratio, the marginal instability boundary consists of the increasing branch and the decreasing branch. A gap region exists between the branches where the critical Marangoni number of the corresponding axisymmetric stationary flow increases drastically. Particularly, a unique axisymmetric oscillatory flow (the critical azimuthal wave number is m=0) in the gap region is reported for the liquid bridge of aspect ratio 1.6. Moreover, the energy transfer between the basic state and the disturbance fields of the thermocapillary convection is analyzed at the corresponding critical Marangoni number, which reveals different major sources of the energy transfer for the development of the disturbances in regimes of the increasing branch, the gap region and the decreasing branch, respectively.

Coupling of thermocapillary convection and evaporation effect in a liquid layer when the evaporating interface is open to air

The coupling mechanism of thermocapillary convection and evaporation effect in evaporating liquids was studied experimentally. The experiments were carried out to study a thin evaporating liquid layer in a rectangular test cell when the upper surface was open to air. By altering the imposed horizontal temperature differences and heights of liquid layers, the average evaporating rate and interfacial temperature profiles were measured. The flow fields were also visualized by PIV method. For comparison, the experiments were repeated by use of another two non-evaporating liquids to study the influence of evaporation effect. The results reveal evidently the role that evaporation effect plays in the coupling with thermocapillary convection.

气泡/液滴Marangoni迁移数值模拟研究

本文利用Level Set方法,数值模拟了微重力情况下气泡/液滴的Marangoni迁移现象,分析了Marangoni数对迁移速度的影响.数值模拟结果表明,随着Marangoni数的增大,非线性热对流效应的影响会逐渐增大,导致沿相界面温度分布趋于均匀,从而降低迁移运动的驱动力,使气泡/液滴的迁移速度随Marangoni数的增加而逐渐减小.

Pool Boiling on Wire in Microgravity - What Have We Learned

In the past years, steady pool boiling of degassed R113 on thin platinum wires has been studied systematically in our lab, including experiments in long-term microgravity aboard RS-22, in short-term microgravity in the Drop Tower Beijing / NMLC, and in normal gravity on the ground. Slight enhancement of nucleate boiling heat transfer is observed in microgravity, while dramatic changes of bubble behaviors are much evident. The value of CHF in microgravity is lower than that in normal gravity, but it can be predicted well by the Lienhard-Dhir correlation, although the dimensionless radius in the present case is far beyond its initial application range. The scaling of CHF with gravity is thus much different from the traditional viewpoint. Considering the influence of the Marangoni effects, the different characteristics of bubble behaviors in microgravity have been explained. A new bubble departure model has also been proposed, which can predict the whole observation both in microgravity and in normal gravity.