Numerical analysis of LEC growth of GaAs with an axial magnetic field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Experimental and numerical investigations on dissolution and recrystallization processes of GaSb/InSb/GaSb under microgravity and terrestrial conditions

The effects of gravity and crystal orientation on the dissolution of GaSb into InSb melt and the recrystallization of InGaSb were investigated under microgravity condition using a Chinese recoverable satellite and under normal gravity condition on earth. To investigate the effect of gravity on the solid/liquid interface and compositional profiles. a numerical simulation was carried out. The InSb crystal melted at 525 degrees C and then a part of GaSb dissolved into the InSb melt during heating to 706 degrees C and this process led to the formation of InGaSb solution. InGaSb solidified during the cooling process. The experimental and calculation results clearly show that the shape of the solid/liquid interface and compositional profiles in the solution were significantly affected by gravity. Under microgravity, as the Ga compositional profiles were uniform in the radial direction. the interfaces were almost parallel. On the contrary, for normal gravity condition, as large amounts of Ga moved up in the upper region due to buoyancy, the dissolved zone broadened towards gravitational direction. Also. during the cooling process, needle crystals of InGaSb started appearing and the value of x of InxGa1-xSb crystals increased with the decrease of temperature. The GaSb with the (111)B plane dissolved into the InSb melt much more than that of the (111)A plane. (C) 2000 Elsevier Science B.V. All rights reserved.

Simulation on Gain Recovery of Quantum Dot Semiconductor Optical Amplifiers by Rate Equation

The gain recoveries in quantum dot semiconductor optical amplifiers are numerically studied by rate equation models. Similar to the optical pump-probe experiment, the injection of double optical pulses is used to simulate the gain recovery of a weak continuous signal for the QD SOAs. The gain recoveries are fitted by a response function with multiple exponential terms. For the pulses duration of 10 ps, the gain recovery can be described by three exponential terms with the time constants, and for the pulse with the width of 150 fs, the gain recovery can be described by two exponential terms, the reason is that the short pulse does not consume lot of carriers.

Efficient Simulation of Freak Waves in Random Oceanic Sea States

Numerical simulations of freak wave generation are studied in random oceanic sea states described by JONSWAP spectrum. The evolution of initial random wave trains is numerically carried out within the framework of the modified four-order nonlinear Schroedinger equation (mNLSE), and some involved influence factors are also discussed. Results show that if the sideband instability is satisfied, a random wave train may evolve into a freak wave train, and simultaneously the setting of the Phillips parameter and enhancement coefficient of JONSWAP spectrum and initial random phases is very important for the formation of freak waves. The way to increase the generation efficiency of freak waves though changing the involved parameters is also presented.

Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

Numerical analysis of four-wave-mixing based multichannel wavelength conversion techniques in fibers

We numerically investigate four-wave-mixing (FWM) based multichannel wavelength conversion for amplitude-modulated signals, phase-modulated signals, together with mixed amplitude and phase modulated signals. This paper also discusses the influence of stimulated Brillouin scattering (SBS) effects on high-efficiency FWM-based wavelength conversion applications. Our simulation results show that DPSK signals are more suitable for FWM-based multichannel wavelength conversion because the OOK signals will suffer from the inevitable datapattern-dependent pump depletion. In future applications, when the modulation format is partially upgraded from OOK to DPSK, the influence of OOK signals on the updated DPSK signals must be considered when using multichannel wavelength conversion. This influence becomes severe with the increase of OOK channel number. It can be concluded that DPSK signals are more appropriate for both transmission and multichannel wavelength conversion,especially in long haul and high bit-rate system.

Effects and numerical analysis of argon gas flow on the oxygen concentration in Czochralski silicon single crystal growth

Argon gas, as a protective environment and carrier of latent heat, has an important effect on the temperature distribution in crystals and melts. Numeric simulation is a potent tool for solving engineering problems. In this paper, the relationship between argon gas flow and oxygen concentration in silicon crystals was studied systematically. A flowing stream of argon gas is described by numeric simulation for the first time. Therefore, the results of experiments can be explained, and the optimum argon flow with the lowest oxygen concentration can be achieved. (C) 2002 Elsevier Science B.V. All rights reserved.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Numerical Study on the Spatial and Temporal Characteristics of Water Quality in the Taihu Lake of China

Taihu Lake is the third largest fresh water lake in China. With the fast economic development, abundant industrial and agricultural waste water has been discharged into Taihu Lake, causing the eutrophication of the water quality, which greatly affected the water utility. In the past decades, the treatment of Taihu Lake has witnessed limited success. Therefore, it is practically and theoretically significant to study the eutrophication of Taihu Lake. This research has focused on the issue of water quality including the characteristics of spatial and temporal distributions, and the rules of nutrient diffusion in the Taihu lake area. Based on the monitoring data, the basis distribution characteristics of water quality in Taihu Lake are analyzed. Comparing Taihu Lake with other Lakes shows that one important reason for Taihu eutrophication is the long period of water retention. A transporting and diffusing model of Taihu nutrient is developed by combining with the hydrodynamics model. Using the model, the concentration field of the total phosphorus (TP) and the influence of wind-driven current are numerically investigated, which leads to the conclusion that the flow field has a great influence on the spatial and temporal distributions of TP in Taihu Lake. Furthermore, the effect for improving the water quality by the project of water diversion from the Yangtze River to Taihu Lake was analyzed by simulation. The results demonstrate that short-term water diversion cannot improve the water quality of the heavily-polluted Meiliang Bay and the western bank areas of Taihu Lake.

Experimental and numerical study on power consumptions in a double-tube-socket pneumatic conveying system

A new methodology is proposed in this paper to predict the lowest power consumption for a double-tube-socket (DTS) pneumatic conveying system. This methodology is established on both experimental work and numerical simulation. After parametric studies by numerical simulation, the desired conveying cases which have the lowest power consumption were obtained. Finally those cases were carried out in our experimental system. The measured power consumption was close to that predicted. In this paper the experimental work is discussed and the numerical simulation introduced. (c) 2010 Elsevier B.V. All rights reserved.

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Simulation of Banded Spherulite Patterns by Coupled Logistic Map Lattice Model in Three Dimensions

Banded spherulite patterns are simulated in three dimensions by means of a Coupled Logistic map lattice model. The patterns obtained by numerical calculation are consistent with those in experiments. The simulation results also indicate that the hand spacing is decreased with the increase of parameter mu in the Logistic map and increased with the increase of the coupling parameter e for cube lattices, and increased with the increase of the thickness of the lattice for polymer film, which is quite similar to the results in some experiments. Spiral pattern in three dimensions is also shown in this paper, which helps us understand the form of banded spherulite in polymers.

Computer simulation of reactive extrusion processes for free radical polymerization

The reactive extrusion for polymerization is an integrated polymer processing technology. A new semi-implicit iterative algorithm was proposed to deal with the complicated relationships among the chemical reaction, the macromolecular structure and the chemorheological property. Then the numerical computation expressions of the average molecular weight, the monomer conversion, and the initiator concentration were deduced, and the computer simulation of the reactive extrusion process for free radical polymerization was carried out, on basis of which reactive processing conditions can be optimized.

Digital simulation in a thin layer spectroelectrochemical cell with minigrid platinum electrode by microregion approximation explicit finite difference method

The microregion approximation explicit finite difference method is used to simulate cyclic voltammetry of an electrochemical reversible system in a three-dimensional thin layer cell with minigrid platinum electrode. The simulated CV curve and potential scan-absorbance curve were in very good accordance with the experimental results, which differed from those at a plate electrode. The influences of sweep rate, thickness of the thin layer, and mesh size on the peak current and peak separation were also studied by numerical analysis, which give some instruction for choosing experimental conditions or designing a thin layer cell. The critical ratio (1.33) of the diffusion path inside the mesh hole and across the thin layer was also obtained. If the ratio is greater than 1.33 by means of reducing the thickness of a thin layer, the electrochemical property will be far away from the thin layer property.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone). V: Estimation of temperature distribution during cure process

A numerical method to estimate temperature distribution during the cure of epoxy-terminated poly(phenylene ether ketone) (E-PEK)-based composite is suggested. The effect of the temperature distribution on the selection of cure cycle is evaluated using a suggested alternation criterion. The effect of varying heating rate and thickness on the temperature distribution, viscosity distribution and distribution of the extent of cure reaction are discussed based on the combination of the here-established temperature distribution model and the previously established curing kinetics model and chemorheological model. It is found that, for a thin composite (<=10mm) and low heating rate (<=2.5K/min), the effect of temperature distribution on cure cycle and on the processing window for pressure application can be neglected. Low heating rate is of benefit to reduce the temperature gradient. The processing window for pressure application becomes narrower with increasing thicknesses of composite sheets. The validity of the temperature distribution model and the modified processing window is evaluated through the characterization of mechanical and physical properties of E-PEK-based composite fabricated according to different temperature distribution conditions.