Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Experimental Study on Bénard-Marangoni Convection in Double Immiscible Liquid Layer

An experimental investigation of Bénard-Marangoni convection has been performed in double immiscible liquid layers of rectangular configuration on the ground. The two kinds of liquid are 10cst silicon oil and FC-70 respectively. The size of rectangular chamber is 100mm×40mm in horizontal cross-section. The evolution processes of convection are observed in the differential thickness ratio of two liquid layers. The critical temperature difference was measured via the detections of fluid convection by a particle image velocimetry (PIV) in the vertical cross-section of the liquid layer. The critical temperature difference or the critical Marangoni number was given. And the influence of the thickness ratio of two liquid layers on the convection instability was discussed. The evolution processes of patterns and temperature distributions on the interface are displayed by using thermal liquid crystal. The velocity distributions on the interface were also obtained. In comparison with the thermocapillary effect, the effect of buoyancy convection will relatively increase when the depth of the liquid layer increases. Because of the coupling of buoyancy and thermocapillary effect, the convection instability is much more complex than that in the microgravity environment. And the critical convection depends on the change of the thickness of liquid layers and also the change of thickness ratio of two liquid layers.

Two Bubbles Thermocapillary Migration Experimental System

An experimental investigation will be performed on the thermocapillary motion of two bubbles in Chinese return-satellite. The experiment will study the migration process of bubble caused by thermocapillary effect in microgravity environment, and their interaction between two bubbles. The bubble is driven by the thermocapillary stress on the surface on account on the variation of the surface tension with temperature. The interaction between two bubbles becomes significant as the separation distance between them is reduced drastically so that the bubble interaction has to be considered. Recently, the problem has been discussed on the method of successive reflections, and accurate migration velocities of two arbitrarily oriented bubbles were derived for the limit of small Marangoni and Reynolds numbers. Numerical results for the migration of the two bubbles show that the interaction between two bubbles has significant influence on their thermocapillary migration velocities with a bubble approaching another. However, there is a lack of experimental validate for the theoretic results. Now the experimental facility is designed for experimenting time after time. A cone-shaped top cover is used to expel bubble from the cell after experiment. But, the cone-shaped top cover can cause temperature uniformity on horizontal plane in whole cell. Therefore, a metal board with multi-holes is fixed under the top cover. The board is able to let the temperature distribution on the board uniform because of their high heat conductivity, and the bubble can pass through it. In the system two bubbles are injected into the test cell respectively by two sets of cylinder. And the bubbles sizes are controlled by two sets of step-by-step motor. It is very important problem that bubble can be divorced from the injecting mouth in microgravity environment. Thus, other two sets of device for injecting mother liquid were used to push bubble. The working principle of injecting mother liquid is to utilize pressure difference directly between test cell and reservoir

Influence of Evaporation on Marangoni-Benard Instability

Evaporative convection and instability give rise to both scientific and technological interests. Practically, a number of the industrial applications such as thin-film evaporators, boiling technologies and heat pipes concern with the evaporation process of which through the vapor-liquid interface the heat and mass transfer occur. From a physical viewpoint, one of interesting questions is the mechanisms of convection instability in thin-liquid layers induced by the coupling of evaporation phenomenon and Marangoni effect at the mass exchanged interface. Classical theories, including Rayleigh’s and Pearson’s, have only successfully explained convection in a liquid layer heated from below without evaporation. However these theories are unable to explain the convection in an evaporating thin layer, especially liquid layer is cooled from below. In present paper, a new two-sided model is put forward rather than the one-sided model in previous works. In previous works, the vapor is treated as passive gas and dynamics of vapor has been ignored. In this case, the vapor liquid system can be described by one-sided model. In our two-sided model, the dynamics of vapor should be considered. Linear instability analysis of the Marangoni-Bénard convection in the two-layer system with an evaporation interface is performed. We define a new evaporating Biot number which is different from the Biot number in one-sided model and obtain the curves of critical Marangoni number versus wave number. In our theoretical results, the Biot number and the evaporating velocity play a major role in the stability of the vapor-liquid system.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Oscillatory Instabilities of Two-Layer Rayleigh-Marangoni-Benard Convection

The Rayleigh-Marangoni-Benard convective instability (R-M-B instability) in the two-layer systems such as Silicone oil (10cSt)/Fluorinert (FC70) and Silicone oil (2cSt)/water liquids are studied. Both linear instability analysis and nonlinear instability analysis (2D numerical simulation) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. The results show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. The secondary instability phenomenon found in the real two-layer system of Silicone oil over water could explain the difference in the comparison of the Degen's experimental observation with the previous linear stability analysis results of Renardy et al.

Two bifurcation processes for onset of oscillatory thermocapillary convection in a floating half zone

An overview on the onset of thermocapillary oscillatory convection in a floating half zone is provided, and it is a typical subject in the microgravity sciences related to the space materials science, especially the floating zone processing, and also to the microgravity fluid physics. The main interests are focused around the process for onset of oscillatory thermocapillary convection, which is known also as the bifurcation transition from quasi-steady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, such as the Marangoni number, Prandtl number, geometrical parameters, and heat transfer parameters. Recent studies show that, there exists the bifurcation transition from steady and axial symmetric convection to the steady and axial non-symmetric convection before the onset of oscillation in cases of small Prandtl number fluids and in cases of larger Prandtl number fluids of fat liquid bridge with small aspect ratio. The transition process is a strong non-linear process because the velocity deviation has the same order of magnitude as that of an average flow after the onset of oscillation, and unsteady 3-D numerical simulation is suitable to do in depth analysis on strong non-linear process, and leads generally to a better comparison with the experimental results.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

Instability Analysis Of Torsional Mems/Nems Actuators Under Capillary Force

The static and dynamic instabilities of a torsional MEMS/NEMS actuator caused by capillary effects are studied, respectively. An instability number, eta, is defined, and the critical gap distance, g(cr), between the mainplate and the substrate is derived. According to the values of eta and g, the instability criteria of the actuator are presented. The dimensionless motion equation of the MEMS/NEMS torsional actuator is derived when it makes nonlinear oscillation under capillary force. The qualitative analysis of the nonlinear equation is made, and the phase portraits are presented on the phase plane. In addition, the bifurcation phenomena in the system are also analyzed. (C) 2008 Elsevier Inc. All rights reserved.

On the origin of shear banding instability in metallic glasses

To uncover the physical origin of shear-banding instability in metallic glass (MG), a theoretical description of thermo-mechanical deformation of MG undergoing one-dimensional simple shearing is presented. The coupled thermo-mechanical model takes into account the momentum balance, the energy balance and the dynamics of free volume. The interplay between free-volume production and temperature increase being two potential causes for shear-banding instability is examined on the basis of the homogeneous solution. It is found that the free-volume production facilitates the sudden increase in the temperature before instability and vice versa. A rigorous linear perturbation analysis is used to examine the inhomogeneous deformation, during which the onset criteria and the internal length and time scales for three types of instabilities, namely free-volume softening, thermal softening and coupling softening, are clearly revealed. The shear-banding instability originating from sole free-volume softening takes place easier and faster than that due to sole thermal softening, and dominates in the coupling softening. Furthermore, the coupled thermo-mechanical shear-band analysis does show that an initial slight distribution of local free volume can incur significant strain localization, producing a shear band. During such a localization process, the local free-volume creation occurs indeed prior to the increase in local temperature, indicating that the former is the cause of shear localization, whereas the latter is its consequence. Finally, extension of the above model to include the shear-induced dilatation shows that such dilatation facilitates the shear instability in metallic glasses.

Self-instability and bending behaviors of nano plates

This paper aims at investigating the size-dependent self-buckling and bending behaviors of nano plates through incorporating surface elasticity into the elasticity with residual stress fields. In the absence of external loading, positive surface tension induces a compressive residual stress field in the bulk of the nano plate and there may be self-equilibrium states corresponding to the plate self-buckling. The self-instability of nano plates is investigated and the critical self-instability size of simply supported rectangular nano plates is determined. In addition, the residual stress field in the bulk of the nano plate is usually neglected in the existing literatures, where the elastic response of the bulk is often described by the classical Hooke’s law. The present paper considered the effect of the residual stress in the bulk induced by surface tension and adopted the elasticity with residual stress fields to study the bending behaviors of nano plates without buckling. The present results show that the surface effects only modify the coefficients in corresponding equations of the classical Kirchhoff plate theory.

Influence of Rayleigh-Taylor Instability on Liquid Propellant Reorientation in a Low-Gravity Environment

A computational simulation is conducted to investigate the influence of Rayleigh-Taylor instability on liquid propellant reorientation flow dynamics for the tank of CZ-3A launch vehicle series fuel tanks in a low-gravity environment. The volume-of-fluid (VOF) method is used to simulate the free surface flow of gas-liquid. The process of the liquid propellant reorientation started from initially flat and curved interfaces are numerically studied. These two different initial conditions of the gas-liquid interface result in two modes of liquid flow. It is found that the Rayleigh-Taylor instability can be reduced evidently at the initial gas-liquid interface with a high curve during the process of liquid reorientation in a low-gravity environment.

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.

Self-focusing, Raman, and modulation instability of shaped relativistic laser pulse in plasmas

The interaction of shaped laser pulses with plasmas is studied in a strict theoretical framework without adopting the slow-varying envelope approximation (SVEA). Any physical quantities involved in the interaction are denoted as a summation of different real quantities of respective phases. The relationships among the phases of those real quantities and their moduli are strictly analyzed. Such strict analyses lead to a more exact equation set for the three-dimensional envelope of the laser pulse, which is not based on SVEA. Based on this equation set, self-focusing, Raman, and modulation instabilities could be discussed in a unified framework. The solutions of this equation set for the laser envelope reveal many possible multicolor laser modes in plasmas. The energy and the shape of a pulse determine its propagation through plasmas in a multicolor mode or in a monochromic mode. A global growth rate is introduced to measure the speed of the transition from the monochromic mode in vacuum to a possible mode in plasmas. (c) 2006 American Institute of Physics.