Transport Phenomena in Ammonothermal Growth of Gallium Nitride

GaN can be used to fabricate blue/green/UV LEDs and high temperature, high power electronic devices. Ideal substrates are needed for high quality III-nitride epitaxy, which is an essential step for the manufacture of LEDs. GaN substrates are ideal to be lattice matched and isomorphic to nitride-based films. Bulk single crystals of GaN can be grown from supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. In this process, a mineralizer such as amide, imide or azide is used to attack a bulk nitride feedstock at temperatures from 200 - 500癈 and pressures from 1 - 4 kbar. Baffle design is essential for successful growth of GaN crystals. Baffle is used to separate the dissolving zone from the growth zone, and to maintain a temperature difference between the two zones. For solubility curve with a positive coefficient with respect to temperature, the growth zone is maintained at a lower temperature than that in the dissolving zone, thus the nutrient becomes supersaturated in the growth zone. The baffle opening is used to control the mixing of nutrients in the two zones, thus the transfer of nutrient from the lower part to the upper part. Ammonothermal systems have been modeled here using fluid dynamics, thermodynamics and heat transfer models. The nutrient is considered as a porous media bed and the flow is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. We investigated the effects of baffle opening and position on the transport phenomena of nutrient from dissolving zone to the growth zone. Simulation data have been compared qualitatively with experimental data.

Determination of plastic properties by instrumented spherical indentation: Expanding cavity model and similarity solution approach

The present paper aims to develop a robust spherical indentation-based method to extract material plastic properties. For this purpose, a new consideration of-piling-up effect is incorporated into the expanding cavity model; an extensive numerical study on the similarity Solution has also been performed. As a consequence, two semi-theoretical relations between the indentation response and material plastic properties are derived, with which plastic properties of materials can be identified from a single instrumented spherical indentation curve, the advantage being that this approach no longer needs estimations of contact radius with given elastic modulus. Moreover, the inconvenience in using multiple indenters with different tip angles can be avoided. Comprehensive sensitivity analyses show that the present algorithm is reliable. Also, by experimental verification performed oil three typical materials, good agreement of the material properties between those obtained from the reverse algorithm and experimental data is obtained.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

A Breakdown Criterion of Free Molecular Flows and an Optimum Analysis of EBPVD

Two important issues in electron beam physical vapor deposition (EBPVD) are addressed. The first issue is a validity condition of the classical cosine law widely used in the engineering context. This requires a breakdown criterion of the free molecular assumption on which the cosine law is established. Using the analytical solution of free molecular effusion flow, the number of collisions (N-c) for a particle moving from an evaporative source to a substrate is estimated that is proven inversely proportional to the local Knudsen number at the evaporation surface. N-c = 1 is adopted as a breakdown criterion of the free molecular assumption, and it is verified by experimental data and DSMC results. The second issue is how to realize the uniform distributions of thickness and component over a large-area thin film. Our analysis shows that at relatively low evaporation rates the goal is easy achieved through arranging the evaporative source positions properly and rotating the substrate.

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.

Nonlinear Fluid Damping In Structure-Wake Oscillators In Modeling Vortex-Induced Vibrations

A Nonlinear Fluid Damping (NFD) in the form of the square-velocity is applied in the response analysis of Vortex-induced Vibrations (VIV). Its nonlinear hydrodynamic effects oil the coupled wake and structure oscillators are investigated. A comparison between the coupled systems with the linear and nonlinear fluid dampings and experiments shows that the NFD model can well describe response characteristics, such as the amplification of body displacement at lock-in and frequency lock-ill, both at high and low mass ratios. Particularly, the predicted peak amplitude of the body in the Griffin plot is ill good agreement with experimental data and empirical equation, indicating the significant effect of the NFD on the structure motion.

Determination Of Plastic Properties By Instrumented Spherical Indentation: Expanding Cavity Model And Similarity Solution Approach

The present paper aims to develop a robust spherical indentation-based method to extract material plastic properties. For this purpose, a new consideration of-piling-up effect is incorporated into the expanding cavity model; an extensive numerical study on the similarity Solution has also been performed. As a consequence, two semi-theoretical relations between the indentation response and material plastic properties are derived, with which plastic properties of materials can be identified from a single instrumented spherical indentation curve, the advantage being that this approach no longer needs estimations of contact radius with given elastic modulus. Moreover, the inconvenience in using multiple indenters with different tip angles can be avoided. Comprehensive sensitivity analyses show that the present algorithm is reliable. Also, by experimental verification performed oil three typical materials, good agreement of the material properties between those obtained from the reverse algorithm and experimental data is obtained.

Measurement of velocity profiles in a rectangular microchannel with aspect ratio α = 0.35

In this work, we measured 14 horizontal velocity profiles along the vertical direction of a rectangular microchannel with aspect ratio alpha = h/w = 0.35 (h is the height of the channel and w is the width of the channel) using microPIV at Re = 1.8 and 3.6. The experimental velocity profiles are compared with the full 3D theoretical solution, and also with a Poiseuille parabolic profile. It is shown that the experimental velocity profiles in the horizontal and vertical planes are in agreement with the theoretical profiles, except for the planes close to the wall. The discrepancies between the experimental data and 3D theoretical results in the center vertical plane are less than 3.6%. But the deviations between experimental data and Poiseuille's results approaches 5%. It indicates that 2D Poiseuille profile is no longer a perfect theoretical approximation since a = 0.35. The experiments also reveal that, very near the hydrophilic wall (z = 0.5-1 mu m), the measured velocities are significantly larger than the theoretical velocity based on the no-slip assumption. A proper discussion on some physical effects influencing the near wall velocity measurement is given.

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 Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

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.

Parallel Fully-Implicit Computation of MHD Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

Subcooled pool boiling on thin wire in microgravity

A new set of experimental data of subcooled pool boiling on a thin wire in rnicrogravity aboard the 22nd Chinese recoverable satellite is reported in the present paper. The temperature-control led heating method is used. The results of the experiments in normal gravity before and after the flight experiment are also presented, and compared with those in microgravity. The working fluid is degassed R113 at 0.1 MPa and subcooled by 26 degrees C nominally. A thin platinum wire of 60 mu m in diameter and 30 mm in length is simultaneously used as heater and thermometer. It is found that the heat transfer of nucleate pool boiling is slightly enhanced in microgravity comparing with those in normal gravity. It is also found that the correlation of Lienhard and Dhir can predict the CHF with good agreement, although the range of the dimensionless radius is extended by three or more decades above the originally set limit. Three critical bubble diameters are observed in microgravity, which divide the observed vapor bubbles into four regimes with different sizes. Considering the Marangoni effect, a qualitative model is proposed to reveal the mechanism underlying the bubble departure processes, and a quantitative agreement can also be acquired.

Laminar Burning Velocities and Markstein Lengths of Premixed Methane/Air Flames near the Lean Flammability Limit in Microgravity

Effects of flame stretch on the laminar burning velocities of near-limit fuel-lean methane/air flames have been studied experimentally using a microgravity environment to minimize the complications of buoyancy. Outwardly propagating spherical flames were employed to assess the sensitivities of the laminar burning velocity to flame stretch, represented by Markstein lengths, and the fundamental laminar burning velocities of unstretched flames. Resulting data were reported for methane/air mixtures at ambient temperature and pressure, over the specific range of equivalence ratio that extended from 0.512 (the microgravity flammability limit found in the combustion chamber) to 0.601. Present measurements of unstretched laminar burning velocities were in good agreement with the unique existing microgravity data set at all measured equivalence ratios. Most of previous 1-g experiments using a variety of experimental techniques, however, appeared to give significantly higher burning velocities than the microgravity results. Furthermore, the burning velocities predicted by three chemical reaction mechanisms, which have been tuned primarily under off-limit conditions, were also considerably higher than the present experimental data. Additional results of the present investigation were derived for the overall activation energy and corresponding Zeldovich numbers, and the variation of the global flame Lewis numbers with equivalence ratio. The implications of these results were discussed. 2010 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.