Progress In Modeling Of Fluid Flows In Crystal Growth Processes

Modeling of fluid flows in crystal growth processes has become an important research area in theoretical and applied mechanics. Most crystal growth processes involve fluid flows, such as flows in the melt, solution or vapor. Theoretical modeling has played an important role in developing technologies used for growing semiconductor crystals for high performance electronic and optoelectronic devices. The application of devices requires large diameter crystals with a high degree of crystallographic perfection, low defect density and uniform dopant distribution. In this article, the flow models developed in modeling of the crystal growth processes such as Czochralski, ammonothermal and physical vapor transport methods are reviewed. In the Czochralski growth modeling, the flow models for thermocapillary flow, turbulent flow and MHD flow have been developed. In the ammonothermal growth modeling, the buoyancy and porous media flow models have been developed based on a single-domain and continuum approach for the composite fluid-porous layer systems. In the physical vapor transport growth modeling, the Stefan flow model has been proposed based on the flow-kinetics theory for the vapor growth. In addition, perspectives for future studies on crystal growth modeling are proposed. (c) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Oil-Water Two-Phase Flow Inside T-Junction

The oil/water two-phase flow inside T-junctions was numerically simulated with a 3-D two-fluid model, and the turbulence was described using the mixture k - epsilon model. Some experiments of oil/water flow inside a single T-junction were conducted in the laboratory. The results show that the separating performance of T-junction largely depends oil the inlet volumetric fraction and flow patterns. A reasonable agreement is reached between the numerical simulation and the experiments for both the oil fraction distribution and the separation efficiency.

Modified Generalized Beam Lattice Model Associated With Fracture Of Reinforced Fiber/Particle Composites

Lattice-type model can simulate in a straightforward manner heterogeneous brittle media. Mohr-Coulomb failure criterion has recently been involved into the generalized beam (GB) lattice model, and as a result, numerical experiments on concrete under various loading conditions can be conducted. The GB lattice model is further used to investigate the reinforced fiber/particle composites instead of only particle composites as the model did before. Numerical examples are given to show the effectiveness of the modeling procedure, and influences of inclusions (particle, fiber and rebar) on the fracture processes are also discussed. (c) 2008 Elsevier Ltd. All rights reserved.

On-Board Experimental Study Of Bubble Thermocapillary Migration In A Recoverable Satellite

An on-board space experiment of bubble thermocapillary migration was performed in the Chinese 22nd recoverable satellite in 2005. Silicone oil of nominal viscosity 5cSt was used as the continuous phase in the experiment. Air bubbles were injected into the liquid in the same direction as the constant temperature gradient in the liquid. The velocities of bubbles were obtained by recording the paths of the bubbles. The results indicate that the scaled velocity of bubbles decreases with an increase of the Marangoni number extended to 9288, which agrees with the results of previous space experiments and numerical simulation. In addition, the interaction between two bubbles was also observed in the space experiment. The trajectories and the velocities of the bubbles were obtained. The two-bubble experiment results are also consistent with the theoretical analysis.

Fractal Pattern Formation In Anodic Bonding Of Pyrex Glass/Al/Si

Anodic bonding of Pyrex glass/Al/Si is an important bonding technique in micro/nanoelectromechanical systems (MEMS/NEMS) industry. The anodic bonding of Pyrex 7740 glass/Aluminum film/Silicon is completed at the temperature from 300 degrees C to 375 degrees C with a bonding voltage between 150 V and 450 V. The fractal patterns are formed in the intermediate Al thin film. This pattern has the fractal dimension of the typical two-dimensional diffusion-limited aggregation (2D DLA) process, and the fractal dimension is around 1.7. The fractal patterns consist of Al and Si crystalline grains, and their occurrences are due to the limited diffusion, aggregation, and crystallization of Si and Al atoms in the intermediate Al layers. The formation of the fractal pattern is helpful to enhance the bonding strength between the Pyrex 7740 glass and the aluminum thin film coated on the crystal silicon substrates.

Surface Stress Induced By Adatoms Adsorption On Solid Surfaces

Based on the statistical thermodynamics theory, a theoretical model of adsorbate induced surface stress of adatoms adsorption on solid surface is presented. For the low coverage, the interaction between the adsorbed molecules is entirely negligible and the adsorption induced surface stress is found to be the function of the coverage and the adsorption energy change with strain. For the high coverage, the adsorbate-adsorbate interaction contributes to the adsorption-induced surface stress effectively. In the case of carbon adsorption on the Ni(100) surface, the value of 0.5 is obtained as a characteristic coverage to decide whether to take the interaction between the adsorabtes into consideration and the results also show that the adsorption induces a compressive surface stress.

Residual Stress And Deformation Analysis Of Anodic Bonded Multi-Layer Of Glass And Aluminum

The bonding of glass wafer to aluminum foils in multi-layer assemblies was made by the common anodic bonding process. The bonding was performed at temperatures in the range 350-450 degrees C and with an applied voltage in the range 400-700 V under a pressure of 0.05 MPa. Residual stress and deformation in samples of two-layer (aluminum/glass) and three-layer (glass/aluminum/glass) were analyzed by nonlinear finite element simulation software MARC. The stress and strain varying with cooling time were obtained. The analyzed results show that deformation of the three-layer sample is significantly smaller than that of the two-layer sample, because of the symmetric structure of the three-layer sample. This has an important advantage in MEMS fabrication. The maximum equivalent stresses locate in the transition layer in both samples, which will become weakness in bonded sample.

Receptivity To Free-Stream Disturbance Waves For Hypersonic Flow Over A Blunt Cone

A high-order shock-fitting finite difference scheme is studied and used to do direction numerical simulation (DNS) of hypersonic unsteady flow over a blunt cone with fast acoustic waves in the free stream, and the receptivity problem in the blunt cone hypersonic boundary layers is studied. The results show that the acoustic waves are the strongest disturbance in the blunt cone hypersonic boundary layers. The wave modes of disturbance in the blunt cone boundary layers are first, second, and third modes which are generated and propagated downstream along the wall. The results also show that as the frequency decreases, the amplitudes of wave modes of disturbance increase, but there is a critical value. When frequency is over the critial value, the amplitudes decrease. Because of the discontinuity of curvature along the blunt cone body, the maximum amplitudes as a function of frequencies are not monotone.

Thermocapillary Migration Of Nondeformable Drops

In this paper, we present a numerical study on the thermocapillary migration of drops. The Navier-Stokes equations coupled with the energy conservation equation are solved by the finite-difference front-tracking scheme. The axisymmetric model is adopted in Our simulations, and the drops are assumed to be perfectly spherical and nondeformable. The benchmark simulation starts from the classical initial condition with a uniform temperature gradient. The detailed discussions and physical explanations of migration phenomena are presented for the different values of (1) the Marangoni numbers and Reynolds numbers of continuous phases and drops and (2) the ratios of drop densities and specific heats to those of continuous phases. It is found that fairly large Marangoni numbers may lead to fluctuations in drop velocities at the beginning part of simulations. Finally, we also discuss the influence of initial conditions on the thermocapillary migrations. (C) 2008 American Institute of Physics.

Thermocapillary Flows In Liquid Bridges Of Molten Tin With Small Aspect Ratios

Numerical simulations were conducted to study thermocapillary flows in short half-zone liquid bridges of molten tin with Prandtl number Pr = 0.009, under ramped temperature difference. The spatio-temporal structures in the thermocapillary flows in short half-zone liquid bridges with aspect ratios As = 0.6, 0.8, and 1.0 were investigated. The first critical Marangoni numbers were compared with those predicted by linear stability analyses (LSA). The second critical Marangoni numbers for As = 0.6 and 0.8 were found to be larger than that for As = 1.0. The time evolutions of the thermocapillary flows exhibited unusual features such as a change in the azimuthal wave number during the three-dimensional stationary (non-oscillating) flow regime, a change in the oscillation mode during the three-dimensional oscillatory flow regime, and the decreasing and then increasing of amplitudes in a single oscillation mode. The effects of the ramping rate of the temperature difference on the flow modes and critical conditions were studied as well. In this paper, the experimental observability of the critical conditions was also discussed. (C) 2008 Elsevier Inc. All rights reserved.

Role Of Vertical Component Of Surface Tension Of The Droplet On The Elastic Deformation Of Pdms Membrane

Poly(dimethylsiloxane) (PDMS) has been widely used in lab-on-a-chip and micro- total analysis systems (mu-TAS), thus wetting and electrowetting behaviors of PDMS are of great importance in these devices. PDMS is a kind of soft polymer material, so the elastic deformation of PDMS membrane by a droplet cannot be neglected due to the vertical component of the interfacial tension between the liquid and vapor, and this vertical component of liquid-vapor surface tension is also balanced by the stress distribution within the PDMS membrane. Such elastic deformation and stress distribution not only affect the exact measurement of contact angle, but also have influence on the micro-fluidic behavior of the devices. Using ANSYS code, we simulated numerically the elastic deformation and stress distribution of PDMS membrane on a rigid substrate due to the liquid-vapor surface tension. It is found that the vertical elastic deformation of the PDMS membrane is on the order of several tens of nanometers due to the application of a droplet with a diameter of 2.31 mm, which is no longer negligible for lab-on-a-chip and mu-TAS. The vertical elastic deformation increases with the thickness of the PDMS membrane, and there exists a saturated membrane thickness, regarded as a semi-infinite membrane thickness, and the vertical elastic deformation reaches a limiting value when the membrane thickness is equal to or thicker than such saturated thickness. (C) Koninklijke Brill NV, Leiden, 2008.

Pull-In Instability Of Micro-Switch Actuators: Model Review

The existing three widely used pull-in theoretical models (i.e., one-dimensional lumped model, linear supposition model and planar model) are compared with the nonlinear beam mode in this paper by considering both cantilever and fixed-fixed type micro and nano-switches. It is found that the error of the pull-in parameters between one-dimensional lumped model and the nonlinear beam model is large because the denominator of the electrostatic force is minimal when the electrostatic force is computed at the maximum deflection along the beam. Since both the linear superposition model and the slender planar model consider the variation of electrostatic force with the beam's deflection, these two models not only are of the same type but also own little error of the pull-in parameters with the nonlinear beam model, the error brought by these two models attributes to that the boundary conditions are not completely satisfied when computing the numerical integration of the deflection.

Space Experimental Investigation On Thermocapillary Migration Of Bubbles

Results from a space experiment on bubble thermocapillary migration conducted on board the Chinese 22nd recoverable satellite were presented. Considering the temperature field in the cell was disturbed by the accumulated bubbles, the temperature gradient was corrected firstly with the help of the temperature measurement data at six points and numerical simulation. Marangoni number (Ma) of single bubble migrating in the space experiment ranged from 98.04 to 9288, exceeding that in the previous experiment data. The experiment data including the track and the velocity of two bubble thermocapillary migration showed that a smaller bubble would move slower as it was passed by a larger one, and the smaller one would even rest in a short time when the size ratio was large enough.

Thermocapillary Motion Of Droplets At Large Marangoni Numbers

In this paper, the thermocapillary motion problem of drops is investigated using the axisymmetric model. The front-tracking method is employed to capture the drop interface. We find that the migration velocity of the drop is greatly influenced by the temperature field in the drop when Ma is fairly large (>100), which leads to an increase-decrease migration velocity at the beginning of our simulations. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Experimental Study Of Vortex-Induced Vibrations Of A Pipeline Near An Erodible Sandy Seabed

Based on similarity analyses, a series of experiments have been conducted with a newly established hydro-elastic facility to investigate the transverse vortex-induced vibrations (VIVs) of a submarine pipeline near an erodible sandy seabed under the influence of ocean currents. Typical characteristics of coupling processes between pipe vibration and soil scour in the currents have been summarized for Case 1: pipe is laid above seabed and Case 11: pipe is partially embedded in seabed on the basis of the experimental observations. Pipe vibration and the corresponding local scour are usually two coupled physical processes leading to an equilibrium state. The influence of initial gap-to-diameter ratio (e(0)/D) on the interaction between pipe vibration and local scour has been studied. Experimental results show that the critical values of V-r for the initiation of VIVs of the pipe near an erodible sand bed get bigger with decreasing initial gap-to-diameter ratio within the examined range of e(0)/D (-0.25 < e(0)/D < 0.75). The comparison of the pipe vibrations near an erodible soil with those near a rigid boundary and under wall-free conditions indicates that the vibration amplitudes of the pipe near an erodible sand bed are close to the curve fit under wall-free conditions; nevertheless, for the same stability parameter, the maximum amplitudes for the VIV coupled with local scour increase with the increase of initial gap-to-diameter ratio. (c) 2007 Elsevier Ltd. All rights reserved.