Influence Of Rayleigh Effect Combined With Marangoni Effect On The Onset Of Convection In A Liquid Layer Overlying A Porous Layer

The coupling mechanism of Rayleigh effect and Marangoni effect in a liquid-porous system is investigated using a linear stability analysis. The eigenvalue problem is solved by means of a Chebyshev tau method. Results indicate that there are three coupling modes between the Rayleigh effect and the Marangoni effect for different depth ratios. (C) 2008 Elsevier Ltd. All rights reserved.

Hydrothermal Instability Of Thermocapillary Convection In Large-Prandtl-Number Liquid Bridges Under Microgravity

Linear stability analysis was performed to study the mechanism of transition of thermocapillary convection in liquid bridges with liquid volume ratios ranging from 0.4 to 1.2, aspect ratio of 0.75 and Prandtl number of 100. 2-D governing equations were solved to obtain the steady axi-symmetric basic flow and temperature distributions. 3-D perturbation equations were discretized at the collocation grid points using the Chebyshev-collocation method. Eigenvalues and eigenfunctions were obtained by using the Q-R. method. The predicted critical Marangoni numbers and critical frequencies were compared with data from space experiments. The disturbance of the temperature distribution on the free surface causes the onset of oscillatory convection. It is shown that the origin of instability is related to the hydrothermal origin for convections in large-Prandtl-number liquid bridges. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

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.

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.

The Convective Instabilities in a Liquid-Vapor System with a Non-equilibrium Evaporation Interface

Rayleigh-Marangoni-B,nard instability in a system consisting of a horizontal liquid layer and its own vapor has been investigated. The two layers are separated by a deformable evaporation interface. A linear stability analysis is carried out to study the convective instability during evaporation. In previous works, the interface is assumed to be under equilibrium state. In contrast with previous works, we give up the equilibrium assumption and use Hertz-Knudsen's relation to describe the phase change under non-equilibrium state. The influence of Marangoni effect, gravitational effect, degree of non-equilibrium and the dynamics of the vapor on the instability are discussed.

Spatiotemporal instabilities in nonlinear Kerr media in the presence of arbitrary higher-order dispersions

Spatiotemporal instabilities in nonlinear Kerr media with arbitrary higher-order dispersions are studied by use of standard linear-stability analysis. A generic expression for instability growth rate that unifies and expands on previous results for temporal, spatial, and spatiotemporal instabilities is obtained. It is shown that all odd-order dispersions contribute nothing to instability, whereas all even-order dispersions not only affect the conventional instability regions but may also lead to the appearance of new instability regions. The role of fourth-order dispersion in spatiotemporal instabilities is studied exemplificatively to demonstrate the generic results. Numerical simulations confirm the obtained analytic results. (C) 2002 Optical Society of America.

Linear spatio-temporal instability analysis of ice growth under a falling water film

A linear spatio-temporal stability analysis is conducted for the ice growth under a falling water film along an inclined ice plane. The full system of linear stability equations is solved by using the Chebyshev collocation method. By plotting the boundary curve between the linear absolute and convective instabilities (AI/CI) of the ice mode in the parameter plane of the Reynolds number and incline angle, it is found that the linear absolute instability exists and occurs above a minimum Reynolds number and below a maximum inclined angle. Furthermore, by plotting the critical Reynolds number curves with respect to the inclined angle for the downstream and upstream branches, the convectively unstable region is determined and divided into three parts, one of which has both downstream and upstream convectively unstable wavepackets and the other two have only downstream or upstream convectively unstable wavepacket. Finally, the effect of the Stefan number and the thickness of the ice layer on the AI/CI boundary curve is investigated.

Effects of wall temperature on boundary layer stability over a blunt cone at Mach 7.99

Effects of wall temperature on stabilities of hypersonic boundary layer over a 7-degree half-cone-angle blunt cone are studied by using both direct numerical simulation (DNS) and linear stability theory (LST) analysis. Four isothermal wall cases with Tw/T0= 0.5, 0.7, 0.8 and 0.9, as well as an adiabatic wall case are considered. Results of both DNS and LST indicate that wall temperature has significant effects on the growth of disturbance waves. Cooling the surface accelerates unstable Mack II mode waves and decelerates the first mode (Tollmien–Schlichting mode) waves. LST results show that growth rate of the most unstable Mack II mode waves for the cases of cold wall Tw/T0=0.5 and 0.7 are about 45% and 25% larger than that for the adiabatic wall, respectively. Numerical results show that surface cooling modifies the profiles of rdut/dyn and temperature in the boundary layers, and thus changes the stability haracteristic of the boundary layers, and then effects on the growth of unstable waves. The results of DNS indicate that the disturbances with the frequency range from about 119.4 to 179.1 kHz, including the most unstable Mack modes, produce strong mode competition in the downstream region from about 11 to 100 nose radii. And adiabatic wall enhances the amplitudes of disturbance according to the results of DNS, although the LST indicates that the growth rate of the disturbance of cold wall is larger. That because the growth of the disturbance does not only depend on the development of the second unstable mode.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Oscillatory instability of Rayleigh-Marangoni-Bénard convection in two-layer liquid systems

The oscillatory behaviour of the Rayleigh-Marangoni-Bénard convective instability (R-M-B instability) regarding two combinations of two-layer fluid systems has been investigated theoretically and numerically. For the two-layer system of Silicone oil (10cSt) over Fluorinert (FC70), both linear instability analysis and 2D numerical simulation show that the instability of the system depends strongly on the depth ratio Hr = H1/H2 of the two-layer liquid. The oscillatory regime at the onset of R-M-B convection enlarges with reducing Γ = Ra/Ma values. In the two-layer system of Silicone oil (2cSt) over water, it loses its stability and onsets to steady convection at first, then the steady convection bifurcates to oscillatory convection with increasing Rayleigh number Ra. This behaviour was found through numerical simulation above the onset of steady convection in the case of r = 2.9, ε=(Ra-Ruc)/Rac = 1.0, and Hr = 0.5. Our findings are different from the previous study of the Rayleigh-Benard instability and show the strong effects of the thermocapillary force at the interface on the time-dependent oscillations at or after the onset of convection. We propose a secondary oscillatory instability mechanism to explain the experimental observation of Degen et al. [Phys. Rev. E, 57 (1998), 6647-6659].

Analysis of critical excavation depth for a jointed rock slope using a face-to-face discrete element method

The critical excavation depth of a jointed rock slope is an important problem in rock engineering. This paper studies the critical excavation depth for two idealized jointed rock slopes by employing a face-to-face discrete element method (DEM). The DEM is based on the discontinuity analysis which can consider anisotropic and discontinuous deformations due to joints and their orientations. It uses four lump-points at each surface of rock blocks to describe their interactions. The relationship between the critical excavation depth D-s and the natural slope angle alpha, the joint inclination angle theta as well as the strength parameters of the joints c(r) ,phi(r) is analyzed, and the critical excavation depth obtained with this DEM and the limit equilibrium method (LEM) is compared. Furthermore, effects of joints on the failure modes are compared between DEM simulations and experimental observations. It is found that the DEM predicts a lower critical excavation depth than the LEM if the joint structures in the rock mass are not ignored.

Oscillatory Instability of Rayleigh-Marangoni-Benard Convection in Two-Layer System

The convective instabilities in two or more superposed layers heated from below were studied extensively by many scientists due to several interfacial phenomena in nature and crystal growth application. Most works of them were performed mainly on the instability behaviors induced only by buoyancy force, especially on the oscillatory behavior at onset of convection (see Gershuni et. Al.(1982), Renardy et. Al. (1985,2000), Rasenat et. Al. (1989), and Colinet et. Al.(1994)) . But the unstable situations of multi-layer liquid convection will become more complicated and interesting while considering at the same time the buoyancy effect combined with thermocapillary effect. This is the case in the gravity reduced field or thin liquid layer where the thermocapillary effect is as important as buoyancy effect. The objective of this study was to investigate theoretically the interaction between Rayleigh-Bénard instability and pure Marangoni instability in a two-layer system, and more attention focus on the oscillatory instability both at the onset of convection and with increasing supercriticality. Oscillatory behavious of Rayleigh-Marangoni-Bénard convective instability (R-M-B instability) and flow patterns are presented in the two-layer system of Silicon Oil (10cSt) over Fluorinert (FC70) for a larger various range of two-layer depth ratios (Hr=Hupper/Hdown) from 0.2 to 5.0. Both linear instability analysis and 2D numerical simulation (A=L/H=10) show that the instability of the system depends strongly on the depth ratio of two-layer liquids. The oscillatory instability regime at the onset of R-M-B convection are found theoretically in different regions of layer thickness ratio for different two-layer depth H=12,6,4,3mm. The neutral stability curve of the system displaces to right while we consider the Marangoni effect at the interface in comparison with the Rayleigh-Bénard instability of the system without the Marangoni effect (Ma=0). The numerical results show different regimes of the developing of convection in the two-layer system for different thickness ratios and some differences at the onset of pure Marangoni convection and the onset of Rayleigh-Bénard convections in two-layer liquids. Both traveling wave and standing wave were detected in the oscillatory instability regime due to the competition between Rayleigh-Bénard instability and Marangoni effect. The mechanism of the standing wave formation in the system is presented numerically in this paper. The oscillating standing wave results in the competition of the intermediate Marangoni cell and the Rayleigh convective rolls. In the two-layer system of 47v2 silicone oil over water, a transition form the steady instability to the oscillatory instability of the Rayleigh-Marangoni-Bénard Convection was found numerically above the onset of convection for ε=0.9 and Hr=0.5. We propose that this oscillatory mechanism is possible to explain the experimental observation of Degen et. Al.(1998). Experimental work in comparison with our theoretical findings on the two-layer Rayleigh-Marangoni-Bénard convection with thinner depth for H<6mm will be carried out in the near future, and more attention will be paid to new oscillatory instability regimes possible in the influence of thermocapillary effects on the competition of two-layer liquids

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.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.