3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Modeling of Ammonothermal Growth of Gallium Nitride Single Crystals

Ammonothermal growth of GaN crystals with a retrograde solubility has been modeled and simulated here using fluid dynamics, thermodynamics and heat transfer models. The nutrient is considered as a porous media bed and the flow in the porous charge is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. For the case of retrograde solubility, the charge is put above the baffle. The temperature difference between the dissolving zone and growth zone is found smaller than that applied on the sidewall of autoclave. The baffle opening has a strong effect on the nutrient transport and supersaturation of GaN species in the growth zone.

Single curved fiber sedimentation under gravity

Dynamics of single curved fiber sedimentation under gravity are simulated by using the lattice Boltzmann method. The results of migration and rotation of the curved fiber at different Reynolds numbers are reported. The results show that the rotation and migration processes are sensitive to the curvature of the fiber. (c) 2007 Elsevier Ltd. All rights reserved.

Hydroelectric Voltage Generation Based On Water-Filled Single-Walled Carbon Nanotubes

A DFT/MD mutual iterative method was employed to give insights into the mechanism of voltage generation based on water-fitted single-walled carbon nanotubes (SWCNTs). Our calculations showed that a constant voltage difference of several mV would generate between the two ends of a carbon nanotube, due to interactions between the water dipole chains and charge carriers in the tube. Our work validates this structure of a water-fitted SWCNT as a promising candidate for a synthetic nanoscale power cell, as well as a practical nanopower harvesting device at the atomic level. 

Transport Properties And Induced Voltage In The Structure Of Water-Filled Single-Walled Boron-Nitrogen Nanotubes

Density functional theory/molecular dynamics simulations were employed to give insights into the mechanism of voltage generation based on a water-filled single-walled boron-nitrogen nanotube (SWBNNT). Our calculations showed that (1) the transport properties of confined water in a SWBNNT are different from those of bulk water in view of configuration, the diffusion coefficient, the dipole orientation, and the density distribution, and (2) a voltage difference of several millivolts would generate between the two ends of a SWBNNT due to interactions between the water dipole chains and charge carriers in the tube. Therefore, this structure of a water-filled SWBNNT can be a promising candidate for a synthetic nanoscale power cell as well as a practical nanopower harvesting device.

Formation of Single and Multiple Deformation Twins in Nanocrystalline fcc Metals

Deformation twins are often observed to meet each other to form multi-fold twins in nanostructured face-centered cubic (fcc) metals.Here we propose two types of mechanism for the nucleation and growth of four different single and multiple twins. These mechanisms provide continuous generation of twinning partials for the growth of the twins after ucleation. A relatively high stress or high strain rate is needed to activate these mechanisms, making them more prevalent in nanocrystalline materials than in their coarse-grained counterparts.Experimental observations that support the proposed mechanisms are presented.

A comparative study of Young's modulus of single-walled carbon nanotube by CPMD, MD, and first principle simulations

Carbon nanotubes (CNTs), due to their exceptional magnetic, electrical and mechanical properties, are promising candidates for several technical applications ranging from nanoelectronic devices to composites. Young's modulus holds the special status in material properties and micro/nano-electromechanical systems (MEMS/NEMS) design. The excellently regular structures of CNTs facilitate accurate simulation of CNTs' behavior by applying a variety of theoretical methods. Here, three representative numerical methods, i.e., Car-Parrinello molecular dynamics (CPMD), density functional theory (DFT) and molecular dynamics (MD), were applied to calculate Young's modulus of single-walled carbon nanotube (SWCNT) with chirality (3,3). The comparative studies showed that the most accurate result is offered by time consuming DFT simulation. MID simulation produced a less accurate result due to neglecting electronic motions. Compared to the two preceding methods the best performance, with a balance between efficiency and precision, was deduced by CPMD.

A surface energy model and application to mechanical behavior analysis of single crystals at sub-micron scale

Size effects of mechanical behaviors of materials are referred to the variation of the mechanical behavior due to the sample sizes changing from macroscale to micro-/nanoscales. At the micro-/nanoscale, since sample has a relatively high specific surface area (SSA) (ratio of surface area to volume), the surface although it is often neglected at the macroscale, becomes prominent in governing the energy effect, although it is often neglected at the macroscale, becomes prominent in governing the mechanical behavior. In the present research, a continuum model considering the surface energy effect is developed through introducing the surface energy to total potential energy. Simultaneously, a corresponding finite element method is developed. The model is used to analyze the axial equilibrium strain problem for a Cu nanowire at the external loading-free state. As another application of the model, from dimensional analysis, the size effects of uniform compression tests on the microscale cylinder specimens for Ni and Au single crystals are analyzed and compared with experiments in literatures. (C) 2009 Elsevier B.V. All rights reserved.

Transport properties and induced voltage in the structure of water-filled single-walled boron-nitrogen nanotubes

Density functional theory/molecular dynamics simulations were employed to give insights into the mechanism of voltage generation based on a water-filled single-walled boron-nitrogen nanotube (SWBNNT). Our calculations showed that (1) the transport properties of confined water in a SWBNNT are different from those of bulk water in view of configuration the diffusion coefficient the dipole orientation and the density distribution and (2) a voltage difference of several millivolts would generate between the two ends of a SWBNNT due to interactions between the water dipole chains and charge carriers in the tube. Therefore this structure of a water-filled SWBNNT can be a promising candidate for a synthetic nanoscale power cell as well as a practical nanopower harvesting device.

Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

Effect of humidity on microstructure and properties of YBCO film prepared by TFA-MOD method

Epitaxial YBCO superconducting films were deposited on the single crystal LaAlO3 (001) substrate by metal organic deposition method. All YBCO films were fired at 820 degrees C in humidity range of 2.6%-19.7% atmosphere. Microstructure of YBCO thin films was analyzed by means of X-ray diffraction (XRD) and scanning electron microscopy (SEM). Superconducting properties of YBCO films were measured by four-probe method. XRD results showed that the second phase (such as BaF2)and a-axis-oriented grains existed in the films prepared at 2.6% humidity condition; a-axis-oriented grains increased in the film prepared at higher than 4.2% humidity condition; almost pure c-axias-oriented grains existed in the films fired at 4.2% humidity condition. Morphologies of the YBCO films showed that all films had a smooth and crack-free surface. YBCO film prepared at 4.2% humidity condition showed J(c) value of 3.3 MA/cm(2) at 77 K in self-field.

On optimal message vector length for block single parallel partition algorithm in a three-dimensional ADI solver

It has long been recognized that many direct parallel tridiagonal solvers are only efficient for solving a single tridiagonal equation of large sizes, and they become inefficient when naively used in a three-dimensional ADI solver. In order to improve the parallel efficiency of an ADI solver using a direct parallel solver, we implement the single parallel partition (SPP) algorithm in conjunction with message vectorization, which aggregates several communication messages into one to reduce the communication costs. The measured performances show that the longest allowable message vector length (MVL) is not necessarily the best choice. To understand this observation and optimize the performance, we propose an improved model that takes the cache effect into consideration. The optimal MVL for achieving the best performance is shown to depend on number of processors and grid sizes. Similar dependence of the optimal MVL is also found for the popular block pipelined method.

Membrane Deformability and Membrane Tension of Single Isolated Mitochondria

Mitochondria dynamics is crucial to many biological processes such as mitochondria fusion and fission, which is highly correlated to the mechanics of single mitochondria. However, the mechanobiological coupling of mitochondria has been poorly understood. Here membrane deformability and membrane tension of individual mitochondria isolated from MtDsRed labeled human embryonic T-Rex-293 kidney cells were measured using a micropipette aspiration assay. The results demonstrated that membrane deformation of isolated mitochondria exhibited an elastic transition phase followed by an equilibrium phase, and mitochondrial membrane tension was proportional to the area compressibility. It was also indicated that mitochondrial membrane deformability was significantly affected by physical chemical factors such as osmotic pressure or pH value, and was further correlated to mitochondrial functionality in different respiratory states and Ca2+ regulation. These findings provide a new insight into understanding the mechanical regulation of mitochondrial physiology.

Experimental Study on the Ignition Process of Single Coal Particles at Microgravity

In this paper, the first Chinese microgravity (μ-g) experimental study on coal combustion was introduced. An experimental system used to study the ignition process of single coal particles was built up, complying with the requirements of the 3.5 s drop tower in the National Microgravity Laboratory of China (NMLC). High volatile bituminous and lignite coal particles with diameter of 1.5 and 2.0 mm were tested. The ignition and combustion process was recorded by a color CCD and the particle surface temperature before and at the ignition was determined by the RGB colorimetric method. Comparative experiments were conducted at normal gravity (1-g). The experiments revealed that at different gravity levels, the ignition of all tested coal particles commenced in homogeneous phase, while the shape, structure, brightness and development of the flames, as well as the volatile matter release during the ignition process are different. At μ-g, the part of volatile was released as a jet, while such a phenomenon was barely observed at 1-g. Also, after ignition, flames were more spherical, thicker, laminated and dimmer at μ-g. It was confirmed that ignition temperature decreased as the particle size or volatile content increased. However, contradicted to existing experimental results, provided other experimental conditions except gravity level were the same, ignition temperature of coal particles was about 50–80 K lower at μ-g than that at 1-g.