Numerical Simulation of Heat Transfer and Kinetics in the Bulk Growth of Silicon Carbide

The growth process of 2-inch silicon carbide (SiC) single crystals by the physical vapor transport method (or modified Lely method) has been modeled and simulated. The comprehensive process model incorporates the calculations of radio frequency (RF) induction heating, heat and mass transfer and growth kinetics. The transport equations for electromagnetic field, heat transfer, and species transport are solved using a finite volume-based numerical scheme called MASTRAPP (Multizone Adaptive Scheme for Transport and Phase Change Process). Temperature distribution for a 2-inch growth system is calculated, and the effects of induction heating frequency and current on the temperature distribution and growth rate are investigated. The predicted results have been compared with the experimental data.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Processing Simulation of Virtual Manufacture in Laser Surface Treatmenthang

Processing simulation is at the bottom of the coral technology of VM and is also difficult due to the complexity of mechanism and diversity of parameters. Previously much research has been mainly carried out on the geometrical simulation or physical simulation respectively. The aim of this paper is to study the processing simulation in laser surface treatment based on the mechanism, put forward the architecture of the whole processing simulation and give the models of the processing. As a result the data structure layers in the whole simulation is presented.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

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.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

Design of super-resolution filters with a Gaussian beam in optical data storage systems

Super-resolution filters based on a Gaussian beam are proposed to reduce the focusing spot in optical data storage systems. Both of amplitude filters and pure-phase filters are designed respectively to gain the desired intensity distributions. Their performances are analysed and compared with those based on plane wave in detail. The energy utilizations are presented. The simulation results show that our designed super-resolution filters are favourable for use in optical data storage systems in terms of performance and energy utilization.

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Fiber Dispersion and Nonlinearity Influences on Transmissions of AM and FM Data Modulation Signals in Radio-Over-Fiber System

Transmission properties of data amplitude modulation (AM) and frequency modulation (FM) in radio-over-fiber (RoF) system are studied numerically. The influences of fiber dispersion and nonlinearity on different microwave modulation schemes, including double side band (DSB), single side band (SSB) and optical carrier suppression (OCS), are investigated and compared. The power penalties at the base station (BS) and the eye opening penalties of the recovered data at the end users are both calculated and analyzed. Numerical simulation results reveal that the power penalty of FM can be drastically decreased due to the larger modulation depth it can achieve than that of AM. The local spectrum broadening around subcarrier microwave frequency of AM due to fiber nonlinearity can also be eliminated with FM. It is demonstrated for the first time that the eye openings of the FM recovered data can be controlled by its modulation depths and the coding formats. Negative voltage encoding format was used to further decrease the RF frequency thus increase the fluctuation period considering their inverse relationship.

First charge collection and position-precision data on the medium-resistivity silicon strip detectors before and after neutron irradiation up to 2x10(14) n/cm(2)

Test strip detectors of 125 mu m, 500 mu m, and 1 mm pitches with about 1 cm(2) areas have been made on medium-resistivity silicon wafers (1.3 and 2.7 k Ohm cm). Detectors of 500 mu m pitch have been tested for charge collection and position precision before and after neutron irradiation (up to 2 x 10(14) n/cm(2)) using 820 and 1030 nm laser lights with different beam-spot sizes. It has been found that for a bias of 250 V a strip detector made of 1.3 k Ohm cm (300 mu m thick) can be fully depleted before and after an irradiation of 2 x 10(14) n/cm(2). For a 500 mu m pitch strip detector made of 2.7 k Ohm cm tested with an 1030 nm laser light with 200 mu m spot size, the position reconstruction error is about 14 mu m before irradiation, and 17 mu m after about 1.7 x 10(13) n/cm(2) irradiation. We demonstrated in this work that medium resistivity silicon strip detectors can work just as well as the traditional high-resistivity ones, but with higher radiation tolerance. We also tested charge sharing and position reconstruction using a 1030 nm wavelength (300 mu m absorption length in Si at RT) laser, which provides a simulation of MIP particles in high-physics experiments in terms of charge collection and position reconstruction, (C) 1999 Elsevier Science B.V. All rights reserved.

First charge collection and position-precision data on the medium-resistivity silicon strip detectors before and after neutron irradiation up to 2x10(14) n/cm(2)

Test strip detectors of 125 mu m, 500 mu m, and 1 mm pitches with about 1 cm(2) areas have been made on medium-resistivity silicon wafers (1.3 and 2.7 k Ohm cm). Detectors of 500 mu m pitch have been tested for charge collection and position precision before and after neutron irradiation (up to 2 x 10(14) n/cm(2)) using 820 and 1030 nm laser lights with different beam-spot sizes. It has been found that for a bias of 250 V a strip detector made of 1.3 k Ohm cm (300 mu m thick) can be fully depleted before and after an irradiation of 2 x 10(14) n/cm(2). For a 500 mu m pitch strip detector made of 2.7 k Ohm cm tested with an 1030 nm laser light with 200 mu m spot size, the position reconstruction error is about 14 mu m before irradiation, and 17 mu m after about 1.7 x 10(13) n/cm(2) irradiation. We demonstrated in this work that medium resistivity silicon strip detectors can work just as well as the traditional high-resistivity ones, but with higher radiation tolerance. We also tested charge sharing and position reconstruction using a 1030 nm wavelength (300 mu m absorption length in Si at RT) laser, which provides a simulation of MIP particles in high-physics experiments in terms of charge collection and position reconstruction, (C) 1999 Elsevier Science B.V. All rights reserved.

Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored.

Direct Numerical Simulation of Shock/Turbulent Boundary Layer Interaction in a Supersonic Compression Ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9. The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition. Both the mean wall pressure and the velocity profiles agree with those of the experimental data, which validates the simulation. The turbulent kinetic energy budget in the separation region is analyzed. Results show that the turbulent production term increases fast in the separation region, while the turbulent dissipation term reaches its peak in the near-wall region. The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation. Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble, the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance, but rather the instability of separation bubble.

Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.