3D Finite Volume Scheme for Czochralski Crystal Growth

A general three-dimensional model is developed for simulation of the growth process of silicon single crystals by Czochralski technique. The numerical scheme is based on the curvilinear non-orthogonal finite volume discretization. Numerical solutions show that the flow and temperature fields in the melt are asymmetric and unsteady for 8’’ silicon growth. The effects of rotation of crystal on the flow structure are studied. The rotation of crystal forms the Ekman layer in which the temperature gradient along solid/melt surface is small.

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.

Integration property of a wedge-elliptic cone waverider with scramjet engine

The hypersonic waverider forebody is designed in this paper. For the present waverider, the undersurface is carved out as a stream surface of a hypersonic inviscid flow field around wedge-elliptic cone, and the upper surface is assumed to be a freestream surface. A finite-volume code is used to generate the three-dimensional flow field. The leading edge is determined by satisfying the condition that the lip is situated at the intersection line of shocks.

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.

EQUIVALENT MODEL OF EXPANSION OF CEMENT MORTAR UNDER SULPHATE EROSION

The expansion property of cement mortar under the attack of sulfate ions is studied by experimental and theoretical methods. First, cement mortars are fabricated with the ratio of water to cement of 0.4, 0.6, and 0.8. Secondly, the expansion of specimen immerged in sulphate solution is measured at different times. Thirdly, a theoretical model of expansion of cement mortar under sulphate erosion is suggested by virtue of represent volume element method. In this model, the damage evolution due to the interaction between delayed ettringite and cement mortar is taken into account. Finally, the numerical calculation is performed. The numerical and experimental results indicate that the model perfectly describes the expansion of the cement mortar.

A new method to analyse the Peierls-Nabarro model of an edge dislocation in a rectangular plate

A new method is presented here to analyse the Peierls-Nabarro model of an edge dislocation in a rectangular plate. The analysis is based on the superposition scheme and series expansions of complex potentials. The stress field and dislocation density field on the slip plane can be expressed as the first and the second Chebyshev polynomial series respectively. Two sets of governing equations are obtained on the slip plane and outer boundary of the rectangular plate respectively. Three numerical methods are used to solve the governing equations.

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

A special method to evaluate the strong adhesion between brittle coating and ductile substrate

The evaluation of the interfacial adhesion of coating system has always been a rough task. In this paper, a special testing method of cross-sectional indentation is applied on a model coating system, i.e. electroplated chromium on a steel substrate which is generally regarded as an example of materials pair with strong adhesion. Based on fractography analysis with SEM and interfacial stress simulation with FEM, it is found that interfacial shear stress may induce coating spalling. More interestingly, spalling location is sensitive to substrate pretreatment process. This shows the feasibility of cross-sectional indentation to distinguish interfacial strength at a high level.

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.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Vacuum heating in the interaction of ultrashort, relativistically strong laser pulses with solid targets

An analytical fluid model for vacuum heating during the oblique incidence by an ultrashort ultraintense p-polarized laser on a solid-density plasma is proposed. The steepening of an originally smooth electron density profile as the electrons are pushed inward by the laser is included self-consistently. It is shown that the electrons being pulled out and then returned to the plasma at the interface layer by the wave field can lead to a phenomenon like wave breaking since the front part of the returning electrons always move slower than the trailing part. This can lead to heating of the plasma at the expense of the wave energy. An estimate for the efficiency of laser energy absorption by the vacuum heating is given. It is also found that for the incident laser intensity parameter, a(L)> 0.5, the absorption rate peaks at an incident angle 45 degrees-52 degrees and it reaches a maximum of 30% at a(L)approximate to 1.5.

Anisotropic Bragg diffraction of finite-sized volume holographic grating in photorefractive crystals

Anisotropic diffraction of uniform plane wave by finite-sized volume holographic grating in photorefractive crystals is considered. It is found that the anisotropic diffraction can take place when some special conditions are satisfied. The diffracted image is obtained in experiment for the anisotropic Bragg diffraction in Fe:LiNbO3 crystals. A coupled wave analysis is presented to study the properties of anisotropic diffraction. An analytical integral solution for the amplitudes of the diffracted beams is submitted. A trade off between high diffraction efficiency and the deterioration of reconstruction fidelity is analyzed. Numerical evaluations also show that the finite-sized anisotropic volume grating exhibits strong angular and wavelength selectivity. All the results are useful for the optimizing design of VHOE based on finite-sized volume grating structures. (c) 2006 Elsevier GmbH. All rights reserved.

Prolonged cardiac xenograft survival in guinea pig-to-rat model by a highly active cobra venom factor

A highly active cobra venom factor (CVF) was isolated from the venom of Naja kaouthia by sequential column chromatography. It displays strong anticomplementary activity, and has 1515 U of anti complementary activity per mg protein. A single dose of 0.1 mg/kg CVF given i.v. to rats completely abrogated complement activity for nearly 5 days. Given 0.02 mg/kg of CVF. the complement activity of rats was reduced by more than 96.5% in 6 It. In guinea pig-to-rat heart transplant model, rats treated with a single dose of 0.05 mg/kg CVF had significantly prolonged xenograft survival (56.12 +/- 6.27 h in CVF-treated rats vs. 0.19 +/- 0.07 h in control rats, P < 0.001). (C) 2003 Elsevier Ltd. All rights reserved.

Quantitative structure-property relationships for octanol-air partition coefficients of polychlorinated naphthalenes, chlorobenzenes and p,p '-DDT

The octanol-air partition coefficient (K-OA) is a key descriptor of chemicals partitioning between the atmosphere and environmental organic phases. Quantitative structure-property relationships (QSPR) are necessary to model and predict KOA from molecular structures. Based on 12 quantum chemical descriptors computed by the PM3 Hamiltonian, using partial least squares (PLS) analysis, a QSPR model for logarithms of K-OA to base 10 (log K-OA) for polychlorinated naphthalenes (PCNs), chlorobenzenes and p,p'-DDT was obtained. The cross-validated Q(cum)(2) value of the model is 0.973, indicating a good predictive ability of the model. The main factors governing log K-OA of the PCNs, chlorobenzenes, and p,p'-DDT are, in order of decreasing importance, molecular size and molecular ability of donating/accepting electrons to participate in intermolecular interactions. The intermolecular dispersive interactions play a leading role in governing log K-OA. The more chlorines in PCN and chlorobenzene molecules, the greater the log K-OA values. Increasing E-LUMO (the energy of the lowest unoccupied molecular orbital) of the molecules leads to decreasing log K-OA values, implying possible intermolecular interactions between the molecules under study and octanol molecules. (C) 2002 Elsevier Science Ltd. All rights reserved.

Quantitative structure-property relationships for octanol-air partition coefficients of polychlorinated biphenyls

Based on nine quantum chemical descriptors computed by PM3 Hamiltonian, using partial least squares analysis, a significant quantitative structure-property relationship for the logarithm of octanol-air partition coefficients (log K-OA) of polychlorinated biphenyls (PCBs) was obtained. The cross-validated Q(cum)(2) value of the model is 0.962, indicating a good predictive ability. The intermolecular dispersive interactions and thus the size of the PCB molecules play a key role in governing log K-OA. The greater the size of PCB molecules, the greater the log K-OA values. Increasing E-LUMO (the energy of the lowest unoccupied molecular orbital) values of the PCBs leads to decreasing log K-OA values, indicating possible interactions between PCB and octanol molecules. Increasing Q(Cl)(+) (the most positive net atomic charges on a chlorine atom) and Q(C)(-) (the largest negative net atomic charge on a carbon atom) values of PCBs results in decreasing log K-OA values, implying possible intermolecular electrostatic interactions between octanol and PCB molecules. (C) 2002 Elsevier Science Ltd. All rights reserved.