Modeling of Ammonothermal Growth of Nitrides

Bulk single crystals of GaN and AlN 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 nitride is used to attack a bulk nitride feedstock at temperatures from 200°C to 500°C and pressures from 1 to 4 kbar. 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 fluid flow is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. The effects of particle size on flow pattern and temperature distribution in an autoclave are analyzed.

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.

Analysis of Temperature and Pressure Changes in Liquefied Natural Gas (LNG) Cryogenic Tanks

Liquefied natural gas (LNG) is being developed as a transportation fuel for heavy vehicles such as trucks and transit buses, to lessen the dependency on oil and to reduce greenhouse gas emissions. The LNG stations are properly designed to prevent the venting of natural gas (NG) from LNG tanks, which can cause evaporative greenhouse gas emissions and result in fluctuations of fuel flow and changes of fuel composition. Boil-off is caused by the heat added into the LNG fuel during the storage and fueling. Heat can leak into the LNG fuel through the shell of tank during the storage and through hoses and dispensers during the fueling. Gas from tanks onboard vehicles, when returned to LNG tanks, can add additional heat into the LNG fuel. A thermodynamic and heat transfer model has been developed to analyze different mechanisms of heat leak into the LNG fuel. The evolving of properties and compositions of LNG fuel inside LNG tanks is simulated. The effect of a number of buses fueled each day on the possible total fuel loss rate has been analyzed. It is found that by increasing the number of buses, fueled each day, the total fuel loss rate can be reduced significantly. It is proposed that an electric generator be used to consume the boil-off gas or a liquefier be used to re-liquefy the boiloff gas to reduce the tank pressure and eliminate fuel losses. These approaches can prevent boil-off of natural gas emissions, and reduce the costs of LNG as transportation fuel.

Porous materials for thermal management under extreme conditions

A brief analysis is presented of how heat transfer takes place in porous materials of various types. The emphasis is on materials able to withstand extremes of temperature, gas pressure, irradiation, etc., i.e. metals and ceramics, rather than polymers. A primary aim is commonly to maximize either the thermal resistance (i.e. provide insulation) or the rate of thermal equilibration between the material and a fluid passing through it (i.e. to facilitate heat exchange). The main structural characteristics concern porosity (void content), anisotropy, pore connectivity and scale. The effect of scale is complex, since the permeability decreases as the structure is refined, but the interfacial area for fluid-solid heat exchange is, thereby, raised. The durability of the pore structure may also be an issue, with a possible disadvantage of finer scale structures being poor microstructural stability under service conditions. Finally, good mechanical properties may be required, since the development of thermal gradients, high fluid fluxes, etc. can generate substantial levels of stress. There are, thus, some complex interplays between service conditions, pore architecture/scale, fluid permeation characteristics, convective heat flow, thermal conduction and radiative heat transfer. Such interplays are illustrated with reference to three examples: (i) a thermal barrier coating in a gas turbine engine; (ii) a Space Shuttle tile; and (iii) a Stirling engine heat exchanger. Highly porous, permeable materials are often made by bonding fibres together into a network structure and much of the analysis presented here is oriented towards such materials. © 2005 The Royal Society.

Thermal transport in high porosity cellular metal foams

The heat dissipation capability of highly porous cellular metal foams with open cells subject to forced air convection is studied using a combined experimental and analytical approach. The cellular morphologies of six FeCrAlY (an iron-based alloy) foams and six copper alloy foams with a range of pore sizes and porosities are quantified with the scanning electronic microscope and image analysis. Experimental measurements on pressure drop and heat transfer for copper foams are carried out. A numerical model for forced convection across open-celled metal foams is subsequently developed, and the predictions are compared with those measured. Reasonably good agreement with test data is obtained, given the complexity of the cellular foam morphology and the associated momentum/energy transport. The results show that cell size has a more significant effect on the overall heat transfer than porosity. An optimal porosity is obtained based on the balance between pressure drop and overall heat transfer, which decreases as the Reynolds number is increased.

Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0.similar to 4.7) x 10(7)/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30 degrees for a sweepback angle of 67.6 degrees. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10 degrees and secondary separation were detected at deflection of theta greater than or equal to 20 degrees. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements at M greater than or equal to 6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers.

Pool Boiling In Microgravity: Recent Results And Perspectives For The Project Depa-Sj10

Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments have also been performed both in normal gravity and in short-term microgravity in the Drop Tower Beijing. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. In the first case, slight enhancement of heat transfer is observed in microgravity, while diminution is evident for high heat flux in the second one. Lateral motions of bubbles on the heaters are observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drives it to detach from the heaters. The Marangoni effect on the bubble behavior is also discussed. The perspectives for a new project DEPA-SJ10, which has been planned to be flown aboard the Chinese recoverable satellite SJ-10 in the future, are also presented.

Progress In Modeling Of Fluid Flows In Crystal Growth Processes

Modeling of fluid flows in crystal growth processes has become an important research area in theoretical and applied mechanics. Most crystal growth processes involve fluid flows, such as flows in the melt, solution or vapor. Theoretical modeling has played an important role in developing technologies used for growing semiconductor crystals for high performance electronic and optoelectronic devices. The application of devices requires large diameter crystals with a high degree of crystallographic perfection, low defect density and uniform dopant distribution. In this article, the flow models developed in modeling of the crystal growth processes such as Czochralski, ammonothermal and physical vapor transport methods are reviewed. In the Czochralski growth modeling, the flow models for thermocapillary flow, turbulent flow and MHD flow have been developed. In the ammonothermal growth modeling, the buoyancy and porous media flow models have been developed based on a single-domain and continuum approach for the composite fluid-porous layer systems. In the physical vapor transport growth modeling, the Stefan flow model has been proposed based on the flow-kinetics theory for the vapor growth. In addition, perspectives for future studies on crystal growth modeling are proposed. (c) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Bubble Dynamics In Nucleate Pool Boiling On Thin Wires In Microgravity

A temperature-controlled pool boiling (TCPB) device has been developed to study the bubble behavior and heat transfer in pool boiling phenomenon both in normal gravity and in microgravity. A thin platinum wire of 60 mu m in diameter and 30 mm in length is simultaneously used as heater and thermometer. The fluid is R113 at 0.1 MPa and subcooled by 26 degrees C nominally for all cases. Three modes of heat transfer, namely single-phase natural convection, nucleate boiling, and two-mode transition boiling, are observed in the experiment both in microgravity aboard the 22nd Chinese recoverable satellite and in normal gravity on the ground before and after the space flight. Dynamic behaviors of vapor bubbles observed in these experiments are reported and analyzed in the present paper. In the regime of fully developed nucleate boiling, the interface oscillation due to coalescence of adjacent tiny bubbles is the primary reason of the departure of bubbles in microgravity. On the contrary, in the discrete bubble regime, it's observed that there exist three critical bubble diameters in microgravity, dividing the whole range of the observed bubbles into four regimes. Firstly, tiny bubbles are continually forming and growing on the heating surface before departing slowly from the wire when their sizes exceed some value of the order of 10(-1) mm. The bigger bubbles with about several millimeters in diameter stay on the wire, oscillate along the wire, and coalesce with adjacent bubbles. The biggest bubble with diameter of the order of 10 mm, which was formed immediately after the onset of boiling, stays continuously

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.

Effects of Mach number on turbulent separation behaviours induced by blunt fin

An experimental study of the interaction between shock wave and turbulent boundary layer induced by blunt fin has been carried out at M-infinity = 7.8 using oil flow visualization and simultaneous measurements of fluctuating wall pressure and heat transfer. This paper presents the effects of Mach number on turbulent separation behaviours induced by blunt fin.

Electrothermal stress in conductive body with collinear cracks

The temperature and stress field in a thin plate with collinear cracks interrupting an electric current field are determined. This is accomplished by using a complex function method that allows a direct means of finding the distribution of the electric current, the temperature and stress field. Temperature dependency for the heat-transfer coefficient, coefficient of linear expansion and the elastic modulus are considered. As an example, temperature distribution is calculated for an alloy (No. GH2132) plate with two collinear cracks under high temperature. Relationships between the stress, temperature, electric density and crack length are obtained. Crack trajectories emanating from existing crack are predicted by application of the strain energy density criterion which can also be used for finding the load carrying capacity of the cracked plate. (C) 2003 Elsevier Ltd. All rights reserved.

Observation and modeling for terrestrial processes in alpine meadow

The water-heat transfer process between land and atmosphere in Haibei alpine meadow area has been systematically observed. A multi-layer coupling model for land-atmosphere interaction was presented with special attention paid to the moisture transfer in leaf stomata under unsaturated condition. A profound investigation on the physical process of turbulent transfer inside the vegetation has been performed with a revised formula of water absorption for root system. The present model facilitates the study of vertically distributed physical variables in detail. Numerical simulation was conducted according to the transfer process of Kinesia humility meadow in the area of Haibei Alpine Meadow Ecosystem Station, CAS. The calculated results agree well with observation.

Numerical Study on Flow Field and Temperature Distribution in Growth Process of 200 mm Czochralski Silicon Crystals

The melt flow and temperature distribution in a 200 mm silicon Czochralski furnace with a cusp magnetic field was modeled and simulated by using a finite-volume based FLUTRAPP ( Fluid Flow and Transport Phenomena Program) code. The melt flow in the crucible was focused, which is a result of the competition of buoyancy, the centrifugal forces caused by the rotations of the crucible and crystal, the thermocapillary force on the free surfaces and the Lorentz force induced by the cusp magnetic field. The zonal method for radiative heat transfer was used in the growth chamber, which was confined by the crystal surface, melt surface, crucible, heat shield, and pull chamber. It was found that the cusp magnetic field could strength the dominant counter-rotating swirling flow cell in the crucible and reduce the flow oscillation and the pulling-rate fluctuation. The fluctuation of dopant and oxygen concentration in the growing crystal could thus be smoothed.

Drag on an ellipsoid particle in free-molecular plasma flow

Based on the analysis of molecular gas dynamics, the drag and moment acting on an ellipsoid particle of revolution X-2/a(2) + Y-2/a(2) + Z(2)/c(2) = 1, as an example of nonspherical particles, are studied under the condition of free-molecular plasma flow with thin plasma sheaths. A nonzero moment which causes nonspherical particle self-oscillation and self-rotation around its own axis in the plasma flow-similar to the pitching moment in aerodynamics-is discovered for the first time. When the ratio of axis length c/a is unity, the moment is zero and the drag formula are reduced to the well-known results of spherical particles. The effects of the particle-plasma relative velocity, the plasma temperature, and the particle materials on the drag and moment are also investigated.