Experimental Study on Bénard-Marangoni Convection in Double Immiscible Liquid Layer

An experimental investigation of Bénard-Marangoni convection has been performed in double immiscible liquid layers of rectangular configuration on the ground. The two kinds of liquid are 10cst silicon oil and FC-70 respectively. The size of rectangular chamber is 100mm×40mm in horizontal cross-section. The evolution processes of convection are observed in the differential thickness ratio of two liquid layers. The critical temperature difference was measured via the detections of fluid convection by a particle image velocimetry (PIV) in the vertical cross-section of the liquid layer. The critical temperature difference or the critical Marangoni number was given. And the influence of the thickness ratio of two liquid layers on the convection instability was discussed. The evolution processes of patterns and temperature distributions on the interface are displayed by using thermal liquid crystal. The velocity distributions on the interface were also obtained. In comparison with the thermocapillary effect, the effect of buoyancy convection will relatively increase when the depth of the liquid layer increases. Because of the coupling of buoyancy and thermocapillary effect, the convection instability is much more complex than that in the microgravity environment. And the critical convection depends on the change of the thickness of liquid layers and also the change of thickness ratio of two liquid layers.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Heat Transfer of Aviation Kerosene at Supercritical Conditions

The heat transfer characteristics of China no. 3 kerosene were investigated experimentally and analytically under conditions relevant to a regenerative cooling system for scramjet applications. A test facility developed for the present study can handle kerosene in a temperature range of 300-1000 K, a pressure range of 2.6-5 MPa, and a mass How rate range of 10-100 g/s. In addition, the test section was uniquely designed such that both the wall temperature and the bulk fuel temperature were measured at the same location along the flowpath. The measured temperature distributions were then used to analytically deduce the local heat transfer characteristics. A 10-component kerosene surrogate was proposed and employed to calculate the fuel thermodynamic and transport properties that were required in the heat transfer analysis. Results revealed drastic changes in the fuel flow properties and heat transfer characteristics when kerosene approached its critical state. Convective heat transfer enhancement was also found as kerosene became supercritical. The heat transfer correlation in the relatively low-fuel-temperature region yielded a similar result to other commonly used jet fuels, such as JP-7 and JP-8, at compressed liquid states. In the high-fuel-temperature region, near and beyond the critical temperature, heat transfer enhancement was observed; hence, the associated correlation showed a more significant Reynolds number dependency.

Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge

In present study, effect of interfacial heat transfer with ambient gas on the onset of oscillatory convection in a liquid bridge of large Prandtl number on the ground is systematically investigated by the method of linear stability analyses. With both the constant and linear ambient air temperature distributions, the numerical results show that the interfacial heat transfer modifies the free-surface temperature distribution directly and then induces a steeper temperature gradient on the middle part of the free surface, which may destabilize the convection. On the other hand, the interfacial heat transfer restrains the temperature disturbances on the free surface, which may stabilize the convection. The two coupling effects result in a complex dependence of the stability property on the Biot number. Effects of melt free-surface deformation on the critical conditions of the oscillatory convection were also investigated. Moreover, to better understand the mechanism of the instabilities, rates of kinetic energy change and "thermal" energy change of the critical disturbances were investigated (C) 2009 Elsevier Ltd. All rights reserved.

Modelling Study On The Plasma Flow and Heat Transfer In A Laminar Arc Plasma Torch Operating At Atmospheric and Reduced Pressure

A modelling study is performed to investigate the characteristics of both plasma flow and heat transfer of a laminar non-transferred arc argon plasma torch operated at atmospheric and reduced pressure. It is found that the calculated flow fields and temperature distributions are quite similar for both cases at a chamber pressure of 1.0 atm and 0.1 atm. A fully developed flow regime could be achieved in the arc constrictor-tube between the cathode and the anode of the plasma torch at 1.0 atm for all the flow rates covered in this study. However the flow field could not reach the fully developed regime at 0.1 atm with a higher flow rate. The arc-root is always attached to the torch anode surface near the upstream end of the anode, i.e. the abruptly expanded part of the torch channel, which is in consistence with experimental observation. The surrounding gas would be entrained from the torch exit into the torch interior due to a comparatively large inner diameter of the anode channel compared to that of the arc constrictor-tube.

Shroud heat transfer measurements from a rotating cavity with an axial throughflow of air

The paper discusses measurements of heat transfer obtained from the inside surface of the peripheral shroud. The experiments were carried out on a rotating cavity, comprising two 0.985-m-dia disks, separated by an axial gap of 0.065 m and bounded at the circumference by a carbon fiber shroud. Tests were conducted with a heated shroud and either unheated or heated disks. When heated, the disks had the same temperature level and surface temperature distribution. Two different temperature distributions were tested; the surface temperature either increased, or decreased with radius. The effects of disk, shroud, and air temperature levels were also studied. Tests were carried out for the range of axial throughflow rates and speeds: 0.0025 ≤ m ≤ 0.2 kg/s and 12.5 ≤ Ω ≤ 125 rad/s, respectively. Measurements were also made of the temperature of the air inside the cavity. The shroud Nusselt numbers are found to depend on a Grashof number, which is defined using the centripetal acceleration. Providing the correct reference temperature is used, the measured Nusselt numbers also show similarity to those predicted by an established correlation for a horizontal plate in air. The heat transfer from the shroud is only weakly affected by the disk surface temperature distribution and temperature level. The heat transfer from the shroud appears to be affected by the Rossby number. A significant enhancement to the rotationally induced free convection occurs in the regions 2 ≤ Ro ≤ 4 and Ro ≥ 20. The first of these corresponds to a region where vortex breakdown has been observed. In the second region, the Rossby number may be sufficiently large for the central throughflow to affect the shroud heat transfer directly. Heating the shroud does not appear to affect the heat transfer from the disks significantly.

Coupled neutronic thermo-hydraulic analysis of full PWR core with Monte-Carlo based BGCore system

BGCore reactor analysis system was recently developed at Ben-Gurion University for calculating in-core fuel composition and spent fuel emissions following discharge. It couples the Monte Carlo transport code MCNP with an independently developed burnup and decay module SARAF. Most of the existing MCNP based depletion codes (e.g. MOCUP, Monteburns, MCODE) tally directly the one-group fluxes and reaction rates in order to prepare one-group cross sections necessary for the fuel depletion analysis. BGCore, on the other hand, uses a multi-group (MG) approach for generation of one group cross-sections. This coupling approach significantly reduces the code execution time without compromising the accuracy of the results. Substantial reduction in the BGCore code execution time allows consideration of problems with much higher degree of complexity, such as introduction of thermal hydraulic (TH) feedback into the calculation scheme. Recently, a simplified TH feedback module, THERMO, was developed and integrated into the BGCore system. To demonstrate the capabilities of the upgraded BGCore system, a coupled neutronic TH analysis of a full PWR core was performed. The BGCore results were compared with those of the state of the art 3D deterministic nodal diffusion code DYN3D (Grundmann et al.; 2000). Very good agreement in major core operational parameters including k-eff eigenvalue, axial and radial power profiles, and temperature distributions between the BGCore and DYN3D results was observed. This agreement confirms the consistency of the implementation of the TH feedback module. Although the upgraded BGCore system is capable of performing both, depletion and TH analyses, the calculations in this study were performed for the beginning of cycle state with pre-generated fuel compositions. © 2011 Published by Elsevier B.V.

Thermal analysis of microstructural evolution of Al2O3 surface modified layers

Solidification behavior and microstructural evolution of surface modified layers in plasma cladding technique are studied via numerical simulations. Both the coupling effect of temperature and solid volume fraction are considered in the proposed thermal analytical model, by which the transient temperature distributions are calculated and the shape of melting pool is determined. Furthermore, we perform microscopic thermal analysis on the nucleation and growth behaviors of ceramic hardening phases and dendrites, as well as the kinetics of related two-phase flow systems. By comparing with experimental observations, the evolution mechanisms of the morphology of Al2O3 ceramic hardening layer are explained. Based on the above results, a relationship among the scanning velocity of plasma stream, dendritic growth rate and the advancing speed of solid/liquid interface is found, and an energy criterion is proposed for predicting the pushing/engulfing transition of ceramic particles by grain growth fronts. (C) 2009 Elsevier B.V. All rights reserved.

Modeling study to compare the flow and heat transfer characteristics of low-power hydrogen, nitrogen and argon arc-heated thrusters

A modelling study is performed to compare the plasma °ow and heat transfer char- acteristics of low-power arc-heated thrusters (arcjets) for three di®erent propellants: hydrogen, nitrogen and argon. The all-speed SIMPLE algorithm is employed to solve the governing equa- tions, which take into account the e®ects of compressibility, Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. The temperature, veloc- ity and Mach number distributions calculated within the thruster nozzle obtained with di®erent propellant gases are compared for the same thruster structure, dimensions, inlet-gas stagnant pressure and arc currents. The temperature distributions in the solid region of the anode-nozzle wall are also given. It is found that the °ow and energy conversion processes in the thruster nozzle show many similar features for all three propellants. For example, the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appear- ing near the cathode tip; the °ow transition from the subsonic to supersonic regime occurs within the constrictor region; the highest axial velocity appears inside the nozzle; and most of the input propellant °ows towards the thruster exit through the cooler gas region near the anode-nozzle wall. However, since the properties of hydrogen, nitrogen and argon, especially their molecular weights, speci¯c enthalpies and thermal conductivities, are di®erent, there are appreciable di®er- ences in arcjet performance. For example, compared to the other two propellants, the hydrogen arcjet thruster shows a higher plasma temperature in the arc region, and higher axial velocity but lower temperature at the thruster exit. Correspondingly, the hydrogen arcjet thruster has the highest speci¯c impulse and arc voltage for the same inlet stagnant pressure and arc current. The predictions of the modelling are compared favourably with available experimental results.

Effective thermal conductivity of graded nonlinear composites with heat contact resistance

The perturbation expansion method is used to find the effective thermal conductivity of graded nonlinear composites having thermal contact resistance on the inclusion surface. As an example, we have studied the graded composites with cylindrical inclusions immersed in a homogeneous matrix. The thermal conductivity of the cylindrical inclusion is assumed to have a power-law profile of the radial distance r measured from its origin. For weakly nonlinear constitutive relations between the heat flow density q and the temperature field T, namely, q = -mu del T - chi vertical bar del T vertical bar(2) del T, in both the inclusion and the matrix regions, we have derived the temperature distributions using the perturbation expansion method. A nonlinear effective medium approximation of graded composites is proposed to estimate the effective linear and nonlinear thermal conductivities. by considering the temperature singularity on the inclusion surface due to the heat contact resistance. (c) 2006 Elsevier B.V. All rights reserved.

Modeling of a compact plate-fin reformer for methanol steam reforming in fuel cell systems

The characteristics of a compact plate-fin reformer (PFR) which integrates endothermic and exothermic reactions into one unit have been investigated by experiment as well as by numerical simulation. One reforming chamber was integrated with two vaporization chambers and two combustion chambers to constitute a single unit of PFR. In the PFR, which is based on a plate-fin beat exchanger, catalytic combustion of the reforming gas is used to simulate the fuel cell anode off gas (AOG) which supplies the necessary heat for the methanol steam reforming. Temperature distributions in all chambers and composition distribution in reforming chamber have been studied, and the effect of the ratio of H2O/CH3OH on the performance of the PFR has also been investigated. A model of the PFR was derived using a three-dimensional numerical model for a cross-current flow arrangement. Theoretical predictions of the temperature distributions in the PFR were in good agreement with experimental values. In addition, the numerical model was able to accurately predict the methanol conversion and the reformate composition in reforming chamber. © 2005 Elsevier B.V. All rights reserved.

Methanol steam reforming in a compact plate-fin reformer for fuel-cell systems

A compact plate-fin reformer (PFR) consisting of closely spaced plate-fins, in which endothermic and exothermic reactions take place in alternate chambers, has been studied. In the PFR, which was based on a plate-fin heat exchanger, catalytic combustion of the reforming gas, as a simulation of the fuel cell anode off gas (AOG), supplied the necessary heat for the reforming reaction. One reforming chamber, which was for hydrogen production, was integrated with two vaporization chambers and two combustion chambers to constitute a single unit of PFR. The PFR is very compact, easy to be placed and scaled up. The effect of the ratio of H2O/CH3OH on the performance of the PFR has been investigated, and temperature distributions in different chambers were studied. Besides, the stationary behavior of the PFR was also investigated. Heat transfer of the reformer was enhanced by internal plate-fins as well as by external catalytic combustion, which offer both high methanol conversion ratio and low CO concentration. In addition, the fully integrated reformer exhibited good test stability. Based on the PFR, a scale-up reformer was designed and operated continuously for 1000 h, with high methanol conversion ratio and low CO concentration. (c) 2004 International Association for Hydrogen Energy. Published by Elsevier Ltd. All rights reserved.

Simulation of 2-D metal cutting by means of a distributed algorithm

Temperature distributions involved in some metal-cutting or surface-milling processes may be obtained by solving a non-linear inverse problem. A two-level concept on parallelism is introduced to compute such temperature distribution. The primary level is based on a problem-partitioning concept driven by the nature and properties of the non-linear inverse problem. Such partitioning results to a coarse-grained parallel algorithm. A simplified 2-D metal-cutting process is used as an example to illustrate the concept. A secondary level exploitation of further parallel properties based on the concept of domain-data parallelism is explained and implemented using MPI. Some experiments were performed on a network of loosely coupled machines consist of SUN Sparc Classic workstations and a network of tightly coupled processors, namely the Origin 2000.

Thermal analysis of LEDS for liquid crystal display's backlighting

This paper investigated the thermal design of the light emitting diode (LED)onto the board and its packaging. The LED was a 6-lead MultiLED with three chips designed for LCD backlighting and other lighting purposes. A 3D finite element model of this LED was built up and thermal analysis was carried out using the multi physics software package PHYSICA. The modeling results were presented as temperature distributions in each LED, and the predicted junction temperature was used for thermal resistance calculation. The results for the board structure indicated that (1) removing the foil attach decreased the thermal resistance, (2) Increasing the copper foil thickness reduced the thermal resistance. package design indicated that the SMT designed LED with integrated slug gave lower thermal resistance. Pb-free solder material gave lower thermal resistance and junction temperature when compared with conductive adhesive

Theoretical analysis on microwave heating of oil–water emulsions supported on ceramic, metallic or composite plates

A detailed theoretical analysis has been carried out to study efficient heating due to microwaves for one-dimensional (1D) oil–water emulsion samples placed on various ceramic, metallic (reflective) and ceramic–metallic composite supports. Two typical emulsion systems are considered such as oil-in-water (o/w) and water-in-oil (w/o). A preliminary study has been carried out via average power vs emulsion thickness diagram to estimate microwave power absorption within emulsion samples for various cases. The maxima in average power, also termed as ‘resonances’, are observed for specific emulsion thicknesses and the two consecutive resonances of significant magnitudes are termed as R1 and R2 modes. For both o/w and w/o emulsions, it is observed that microwave power absorption is enhanced in presence of metallic and composite supports during both R1 and R2 modes. The efficient heating strategies characterized by ‘large heating rates’ with ‘minimal thermal runaway’ i.e. uniform temperature distributions within the sample have been assessed for each type of emulsion. Based on the detailed spatial distributions of power and temperature for various cases, SiC-metallic composite support may be recommended as an optimal heating strategy for o/w samples with higher oil fractions (0.45) whereas metallic and Alumina-metallic composite supports may be favored for samples with smaller oil fractions (=0.3) during R1 mode. For w/o samples, SiC-metallic composite support may be suitable heating strategy for all ranges of water fractions during R1 mode. During R2 mode, metallic and Alumina-metallic composite supports are favored for both o/w and w/o emulsion samples. Current study recommends the efficient way to use microwaves in a single mode waveguide and the heating strategy can be suitably extended for heating of any other emulsions for which dielectric properties are easily measurable or available in the literature.