Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NATURAL CONVECTION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

A first principles framework to predict the transient performance of latent heat thermal energy storage

Thermal energy storage (TES) is increasingly recognized as an essential component of efficient Combined Heat and Power (CHP), Concentrated Solar Power (CSP), Heating Ventilation and Air Conditioning (HVAC), and refrigeration as it reduces peak demand while helping to manage intermittent availability of energy (e.g., from solar or wind). Latent Heat Thermal Energy Storage (LHTES) is a viable option because of its high energy storage density. Parametric analysis of LHTES in terms of dimensionless numbers is highly desired as a tool to model LHTES systems. One approach is to develop a model equation so as to minimize the error between the model and data obtained from experiments or simulations. While this approach can produce an accurate correlation applicable within the range of data used for its creation, it does not provide physical understanding of the rate-limiting process controlling the transient behavior of the device. In this paper we present an alternative approach whereby the potential rate-limiting processes are identified from first principles and then the key process is determined as a function of time as a LHTES device is charged. For example, in a simple geometry, the melt-fraction can be expected to vary linearly in time if the heat transfer rate is limited by natural convection of the phase changing material and we show it scales with the PCM Grashof number as $Gr^1_p$ and PCM Prandtl number as $Pr_p^{(1/3)}$. On the other hand, if surface area of solid PCM limits the heat transfer rate, the melt fraction increases asymptotically to reach full melting. The existence of these linear and asymptotic regions and the $Gr^1_pP r^{1/3}_p$ shape of the melt fraction curve is verified using our database of 64 simulations. Of practical importance in designing LHTES devices is the melt fraction at which the heat transfer rate ceases to be limited by convection, after which the heat storage rate deteriorates. For our geometry, this is found to be about 90%. This test case of our methodology shows the value of our approach, that predicting heat storage rate based on the rate-limiting physical phenomenon as a function of time is an effective approach to modeling LHTES devices.

25 ENERGY STORAGE↗

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Thermal Model for Horizontal In-Water and In-Air Mini-Plate Experiments in the Advanced Test Reactor

Mini-Plate (MP) irradiation tests are a fueled experiment designed for irradiation in multiple test locations in the Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL). MP experiments are considered a non-instrumented drop-in test where small aluminum-clad fuel plate samples are cooled directly by the ATR Primary Coolant System (PCS) water. These fueled experiments contain aluminum-clad fuel mini plates consisting of monolithic U-10Mo. Previous thermal analyses have been documented concerning oxide growth[1] and thermal safety margin[2] for the MP-2 experiment. A thermal analysis has been performed on the MP-2 experiment using a new technique for calculating the necessary time after the experiment has been irradiated and reactor shutdown necessary to not exceed the 95% lower bound blister threshold temperature for the fuel in a horizontal in air configuration. A horizontal in water model has also been analyzed and will be discussed. After the experiment is irradiated in the ATR, it is necessary to calculate the fuel meat temperatures at 12 hours after shutdown. This is done to assure that at the 12 hour mark, if this experiment is removed from the reactor and accidentally dropped onto the working tray and is in the most severe geometrical configuration of horizontal (minimum heat transfer) that fission products will not be released. The blister threshold temperature is considered this safety limit. The horizontal in air analysis is also performed considering that the experiment is horizontal in a drained canal and filled with air and natural convection heat transfer with air on the outside. The initial mini-plate testing campaign known as MP-1 [3] has been irradiated and its purpose was to form the base fuel specimens fabricated by various processes so that a single process can be down-selected for qualification. MP-2 will then serve to populate a comprehensive fuel performance data set for the down-selected fabrication process. The MP-2 experiment will be irradiated in several locations of the ATR as shown in Figure 1. The fuel plate thickness and geometry is identical to the MP-1 experiment [3], with the only variable being the fuel meat thickness. This model of the B-12 position is treated as encompassing the I-21, I-24, and B-9 positions since it has more fuel and higher total heat than the other positions. Since the blister threshold temperature is a function of fission density, a series of fission densities is calculated for the B-12 position to cover the other three positions. The B 12 position specifically uses fuel meat that is 0.016 in thick by 3.25 in long and 0.75 in tall. The total plate thickness is 0.050 in. Aluminum cladding forms the barrier between the PCS and the fuel meat. The purpose of this paper is to demonstrate this new technique developed to calculate the necessary decay time after shutdown for a series of fission densities using the ABAQUS [4] finite element and heat transfer code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dropwindsonde observations of the mean boundary layer structure in the low-latitude Pacific during SOP-1

The low-level mean thermodynamic and wind profiles in the transition between the trade winds and the ICZ are described using the FGGE level II-b dropwindsonde data set collected during the period of SOP-1 (January 15-February 20, 1979). The results illustrate the influence of decaying midlatitude frontal systems in the northwesternmost region, the general tradelike nature of the eastern regions, and the highly convective nature of a region near the equator. Contrary to the results obtained with the traditional approach to the tropical boundary layer, the FGGE data indicate a trade inversion that remains quasi-horizontal in the north-south direction despite an equatorward decrease in its frequency of occurrence, indicating that subsidence may be more important in equatorial regions than was previously thought.

Firestone, James K.↗

Direct numerical simulation of a combusting droplet with convection

The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.

Liang, Pak-Yan↗

Observations of soot during droplet combustion at low gravity - Heptane and heptane/monochloroalkane mixtures

Experimental observations of the combustion of sooting fuel droplets, performed in a drop tower to create a low gravity environment, are reported. Free n-heptane droplets and suspended droplets of heptane, monochloroalkanes, and mixtures of monochloro-octane and heptane were studied. Initial droplet diameters ranged from 0.4 to 1.1. mm. The results suggest that soot may influence droplet vaporization rates. Spherical symmetry of the flame allowed for extended observation of soot agglomerates inside the droplet flame. Effects of slight convective flows were also observed, both through variations of natural convection around the suspended droplets and through variations in the drift velocities of the unsupported droplets. Slight convective flows around the suspended droplets reduced flame luminosity as well as soot accumulation inside the flame and increased droplet vaporization rates. Mixing monochloro-octane with n-heptane demonstrated the effectiveness of n-heptane in reducing soot emissions from the flames of the chlorinated fuels. Finally, trends of initial droplet diameter with burning rate were observed and may be linked to the effect of droplet size upon soot formation inside the flame.

Jackson, G. S.↗

Thermosolutal convection during dendritic solidification of a binary alloy

The FEM model of Heinrich (1988) is presently extended to the case of dendritic vertical solidification under steady-state conditions. An analysis is conducted of the nonlinear governing equations for linear stability, and FEM calculations are performed in cases where instability has been predicted for Pb-Sn melt systems; the results obtained show that the model is more stable than the plane front liquid-solid model, and shows similar, finger-like covection when unstable. It is found that convection in the 'mush' zone is confined to the top of that region; convective motion is driven by convection in the bulk fluid, in the absence of natural convection.

Heinrich, Juan C.↗

Electric fields in the magnetosphere

Two techniques, tracking the motions of Ba(+) clouds and measuring the differences in floating potential between symmetric double probes, were successful in: (1) demonstrating the basic convective nature of magnetospheric electric fields, (2) mapping global patterns of convection at upper ionosphere levels, and (3) revealing the physics of electric currents in the ionosphere and the importance of magnetosphere-ionosphere feedback in altering the imposed convection.

Heppner, J. P.↗

Secondary arm coarsening and microsegregation in superalloy PWA-1480 single crystals: Effect of low gravity

Single crystal specimens of nickel base superalloy PWA-1480 were directionally solidified on ground and during low gravity (20 sec) and high gravity (90 sec) parabolic maneuver of KC-135 aircraft. Thermal profiles were measured during solidification by two in-situ thermocouples positioned along the sample length. The samples were quenched during either high or low gravity cycles so as to freeze the structures of the mushy zone developing under different gravity levels. Microsegregation was measured by examining the solutal profiles on several transverse cross-sections across primary dendrites along their length in the quenched mushy zone. Effect of gravity level on secondary arm coarsening kinetics and microsegregation have been investigated. The results indicate that there is no appreciable difference in the microsegregation and coarsening kinetics behavior in the specimens grown under high or low gravity. This suggests that short duration changes in gravity/levels (0.02 to 1.7 g) do not influence convection in the interdendritic region. Examination of the role of natural convection, in the melt near the primary dendrite tips, on secondary arm spacings requires low gravity periods longer than presently available on KC-135. Secondary arm coarsening kinetics show a reasonable fit with the predictions from a simple analytical model proposed by Kirkwood for a binary alloy.

Vijayakumar, M.↗

Electric fields in the magnetosphere.

Two techniques, tracking the motions of Ba(+) clouds and measuring the differences in floating potential between symmetric double probes, have been highly successful in: (1) demonstrating the basic convective nature of magnetospheric electric fields, (2) mapping the global patterns of convection at upper ionosphere levels, and (3) revealing the physics of electric currents in the ionosphere and the importance of magnetosphere-ionosphere feedback in altering the imposed convection. The basic pattern of anti-solar convection across the polar cap and night toward day convection in both the evening and morning sectors at auroral belt latitudes persists at all levels of activity. The dawn-dusk potential drop across the polar cap (anti-solar convection) ranges from 20 to 100 kilovolts with the most typical values in the center of this range. The sum of morning and evening (night toward day convection) potential drops in the adjacent auroral belts roughly equals the polar cap drop in the opposite sense as expected.

Heppner, J. P.↗

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗

Tropical Cyclone Diurnal Cycle as Observed by TRMM

Using infrared satellite data, previous work has shown a consistent diurnal cycle in the pattern of cold cloud tops around mature tropical cyclones. In particular, an increase in the coverage by cold cloud tops often occurs in the inner core of the storm around the time of sunset and subsequently propagates outward to several hundred kilometers over the course of the following day. This consistent cycle may have important implications for structure and intensity changes of tropical cyclones and the forecasting of such changes. Because infrared satellite measurements are primarily sensitive to cloud top, the goal of this study is to use passive and active microwave measurements from the Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) and Precipitation Radar (PR), respectively, to examine and better understand the tropical cyclone diurnal cycle throughout a larger depth of the storm's clouds. The National Hurricane Center's best track dataset was used to extract all PR and TMI pixels within 1000 km of each tropical cyclone that occurred in the Atlantic basin between 1998-2011. Then the data was composited according to radius (100-km bins from 0-1000 km) and local standard time (LST; 3-hr bins). Specifically, PR composites involved finding the percentage of pixels with reflectivity greater than or equal to 20 dBZ at various heights (i.e., 2-14 km in increments of 2 km) as a function of radius and time. The 37- and 85- GHz TMI channels are especially sensitive to scattering by precipitation-sized ice in the mid to upper portions of clouds. Hence, the percentage of 37- and 85-GHz polarization corrected temperatures less than various thresholds were calculated using data from all storms as a function of radius and time. For 37 GHz, thresholds of 260 K, 265 K, 270 K, and 275 K were used, and for 85 GHz, thresholds of 200-270 K in increments of 10 K were utilized. Note that convection forced by the interactions of a tropical cyclone with land (e.g., due to frictional convergence) may disrupt the natural convective cycle of a cyclone. Hence, only data pertaining to storms whose centers were greater than 300 km from land were included in the composites. Early results suggest the presence of a diurnal cycle in the PR composites of all Atlantic basin tropical cyclones from a height of 2-12 km from approximately 0-400 km radius, but the cycle is most apparent above 6 km. At a height of 8 km, there is a peak (minimum) in the percentage of PR pixels greater than or equal to 20 dBZ near 0 (21) LST in the inner core with some indication that this signal propagates outward with time. In contrast, the 37- and 85-GHz composites show little indication of a diurnal cycle at any radii, regardless of the threshold used. Ongoing work with this project will involve sub-setting the composites according to storm intensity to see if the diurnal cycle varies with storm strength. Moderate to strong vertical wind shear often leads to asymmetries in tropical cyclone convection and may disrupt the cyclone's natural diurnal cycle. Therefore, wind shear thresholds will be applied to the composites to determine if the diurnal cycle becomes more apparent in a low shear environment. Finally, other work to be completed will involve developing composites for other tropical cyclone basins, including the East Pacific, Northwest Pacific, South Pacific, and Indian Ocean.

tropical cyclones↗

Experimental investigation of Multi-Mode heat transfer to a Free-Falling dilute particle cloud in a heated vertical tube

The development of dilute particle heat exchangers and reactors for advanced energy systems requires an understanding of the multi-mode heat transfer from a heated wall to falling particles. This study presents experimental results of the overall heat transfer coefficient for a free-falling, dilute flow of particles with solid volume fraction from 0.0005 to 0.006 corresponding to feed rates from 3.7 kg s -1 m -2 to 44 kg s -1 m -2 in a vertical, heated tube containing quiescent air at atmospheric pressure. Tube wall temperatures are varied between 300°C to 900°C while keeping the particle inlet temperature constant. The experimental results show that the overall heat transfer coefficient is a strong function of particle feed rate and surface temperature. Good agreement was obtained with prior studies conducted at comparable temperatures but lower particle feed rates (< 4 kg m -2 s -1 ). The established correlations for particle-to-wall radiation and particle-to-gas convection were used to estimate the wall-to-gas convective contribution from the measured overall heat transfer coefficient. The experimental results indicated a 4 to 6 times improvement in the wall convection in the solid-gas mixture compared to that expected from natural convection in a single-phase gas. Furthermore, the data presented here are applicable to characterize heat transfer in dilute particle heat exchangers, furnaces, and solar receivers.

14 SOLAR ENERGY↗

Laminar free convection over two-dimensional bodies with uniform surface heat flux

A general analysis is presented of the steady nondissipative constant-property laminar boundary layer flow over a two-dimensional body of uniform surface heat flux situated in an infinite ambient fluid of undisturbed temperature. The analysis is then applied to a long horizontal circular cylinder. Numerical solutions to the universal functions associated with the first two terms in the derived series are given for Prandtl numbers 0.7, 1, 3, 5, 7, and 10. The results are compared with those obtained by Koh (1964) whose method is patterned after the Blasius-Frossling procedure for forced convection flow. The study reveals that Wilks' (1972) analysis concerning the external natural convection about two-dimensional bodies with constant heat flux is in error.

Lin, F. N.↗

Effect of Pigmentation in Particulate Formation from Fluoropolymer Thermodegradation in Microgravity

Fires aboard spacecraft have occurred as a result of overheated electrical wires and thermodegradation of their insulation, which is composed of fluoropolymers. The particulate products of polymer thermodegradation are only 20-50 run in diameter and are thought to play a role in "polymer fume fever". Therefore an experimental study of the particulates produced by intense ohmic heating of various fluoropolymer-insulated 20 AWG copper wire (representative of spacecraft materials) was undertaken in normal gravity and in microgravity. The 2.2 s drop facility at NASA LeRC and 1.5 s drop facility at the Colorado School of Mines were used to achieve low gravity conditions. Thermophoretic sampling was used for particulate collection. Transmission electron microscopy (TEM) and scanning transmission electron microscopy (STEM) were used to characterize the smoke particulates. It was found that the color of PTFE (Polytetrafluoroethylene) insulation has an overwhelming effect on the size, shape, morphology and, composition of the particulates. Size distributions and shape analyses using computerized image analysis showed that particle size distributions were also dependent on the pigment of the fluoropolymer insulation. The influence of pigment was observed in experiments under both normal and microgravity. Under microgravity conditions, owing to the lack of natural convective transport of particulates, much more particle aggregation was observed, and the nature of the aggregates was dependent on the color of the insulation.

Srivastava, Rajiv↗

Drop evaporation in a single-axis acoustic levitator

A 20 kHz single-axis acoustic positioner is used to levitate aqueous-solution drops (volumes less than or approximately equal to 100 micro-liters). Drop evaporation rates are measured under ambient, isothermal conditions for different relative humidities. Acoustic convection around the levitated sample enhances the mass loss over that due to natural convection and diffusion. A theoretical treatment of the mass flow is developed in analogy to previous studies of the heat transfer from a sphere in an acoustic field. Predictions of the enhanced mass loss, in the form of Nusselt (Sherwood) numbers, are compared with observed rages of drop shrinking. The work is part of an ESA crystal growth from levitated solution drops.

Lierke, E. G.↗

A half-zone study of Marangoni convection in floating-zone crystal growth under microgravity

Floating-zone crystal growth under microgravity, though free from natural convection, can still be susceptible to Marangoni convection. The use of a ring heater in contact with most of the surface of the melt zone, has recently been considered as one way of reducing Marangoni convection in the melt zone. In order to study Marangoni convection in the melt zone in this modified floating-zone crystal growth process, a half-zone system was adopted. Computer simulation and direct observation of Marangoni convection in the half-zone system were carried out for two different cases. In the first case, the surface of the melt zone was completely free, while in the second it was mostly in contact with a quartz ring. The results of computer simulation and direct observation both indicated that in the second case Marangoni convection was significantly reduced near the melt/solid interface as well as in the bulk melt.

Lan, C. W.↗