Influence of liquid bridge volume on instability of floating half zone convection

Thermocapillary instabilities on floating half zone convection in microgravity environment were investigated by linear instability analysis method. The critical Marangoni numbers were obtained and compared with the experimental ones. The influences of the liquid bridge volume and the aspect ratio on the critical Marangoni number were analyzed. It is found that the liquid bridge volume and the aspect ratio have great influence on the critical Marangoni number. There was a gap region where the oscillatory convection will not be observed in present analyses and in experiments in the curve of the critical Marangoni number vs the liquid bridge volume for the case of large Prandtl number and small aspect ratio.

Relation between Stabilization Energy, Crystal Field Coefficient and the Magnetic Exchange Interaction for Tb3+ Ion

Based on a single ion model, Hamiltonian of the simplest form about magnetocrystalline anisotropy for Tb3+ ion was solved by using the numerical method. The relation between the stabilization energy, crystal field coefficient B-2(0) and the magnetic exchange interaction was studied as temperature approaches to 0 K. The results show that the stabilization energy contributed by Tb3+ is linear with crystal field coefficient B-2(0) approximately, but it is insensitive to the change of magnetic exchange interaction for the strong magnetic substances such as TbCo5, Tb2Co17 and Tb2Fe14B compounds.

Effects of Translational Nonequilibrium on Performance of Flowing Chemical Oxygen-Iodine Laser

The effect of the translational nonequilibrium on performance modeling of flowing chemical oxygen-iodine lasers (COIL) is emphasized in this paper. The spectral line broadening (SLB) model is a basic factor for predicting the performances of flowing COIL. The Voigt profile function is a well-known SLB model and is usually utilized. In the case of gas pressure in laser cavity less than 5 torr, a low pressure limit expression of the Voigt profile function is used. These two SLB models imply that ail lasing particles can interact with monochromatic laser radiation. Basically, the inhomogeneous broadening effects are not considered in these two SLB models and they cannot predict the spectral content. The latter requires consideration of finite translational relaxation rate. Unfortunately, it is rather difficult to solve simultaneously the Navier-Stokes (NS) equations and the conservation equations of the number of lasing particles per unit volume and per unit frequency interval. In the operating condition of flowing COIL, it is possible to obtain a perturbational solution of the conservational equations for lasing particles and deduce a new relation between the gain and the optical intensity, i.e., a new gain-saturation relation. By coupling the gain-saturation relation with other governing equations (such as the NS equations, chemical reaction equations and the optical model of gain-equal-loss), We have numerically calculated the performances of flowing COIL. The present results are compared with those obtained by the common rate-equation (RE) model, in which the Voigt profile function and its low pressure limit expression are used. The difference of different model's results is great. For instance, in the case of lasing frequency coinciding with the central frequency of line profile and very low gas pressure, the gain-saturation relation of the present model is quite different with that of the RE model.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection III. Three-dimensional model

An unsteady and three-dimensional model of the floating-half-zone convection on the ground is studied by the direct numerical simulation for the medium of 10 cSt silicon oil, and the influence of the liquid bridge volume on the critical applied temperature difference is especially discussed. The marginal curves for the onset of oscillation are separated into two branches related, respectively, to the slender liquid bridge and the fat liquid bridge. The oscillatory features of the floating-half-zone convection are also discussed.

Influence of Liquid Bridge Volume on the Floating Zone Convection

The liquid bridge volume is a critical geometrical parameter in addition to the aspect ratio for onset of oscillation in the floating zone convection. The oscillatory features are generally divided into two characteristic regions: slender liquid bridge region and fat liquid bridge region. The oscillatory modes in two regions are discussed in the present paper.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

A Modified Free Volume Model For Characterizing Of Rate Effect In Bulk Metallic Glasses

We investigate the plastic deformation and constitutive behaviour of bulk metallic glasses (BMGs). A dimensionless Deborah number De(ID) = t(r)/t(i) is proposed to characterize the rate effect in BMGs, where t(r) is the structural relaxing characteristic time of BMGs under shear load, t(i) is the macroscopic imposed characteristic time of applied stress or the characteristic time of macroscopic deformation. The results demonstrate that the modified free volume model can characterize the strain rate effect in BMGs effectively.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .1. Experiments

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was experimentally studied, and discussions were concentrated on the influence of liquid bridge volume on the onset of oscillation. Distributions of critical applied temperature difference and frequency depending on the volume of the liquid bridge were obtained, and there was a gap range of liquid volume which separated the curve of marginal stability into two parts for fixed rod diameter and aspect ratio. The results imply that the volume of the liquid bridge is a sensitive critical parameter for the onset of oscillation. The implication on the instability is also discussed in the present paper.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .2. Numerical-simulation

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was studied numerically for an unsteady and two-dimensional model, and studies were concentrated on the influence of liquid bridge volume on the onset of oscillation in comparison with the experimental results in the Paper I. The numerical results agree with the experimental ones presented in the previous paper, in which the distributions of critical applied temperature difference depending on the volume of liquid bridge and a gap range of liquid volume in marginal stability curve were obtained.

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.