Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface tension”

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 343 records · Page 19

Department of Transportation PHMSA Aerosol Spray Project: Phase-0

This report discusses the work conducted at Sandia National Laboratories during the preliminary phase (Phase-0) of the Department of Transportation (DOT) Pipeline and Hazardous Materials Safety Administration (PHMSA) project tasked to better understand what happens to a high-pressure liquid canister or vessel during a fire in which an accidental instantaneous release occurs. Phase 0 had two main objects: (1) to establish a methodology for the high temperature and pressure spray vessel release and (2) to characterize the spray droplet formation. For this test series, diagnostic equipment such as the Malvern Spraytec, high-speed photometric cameras, pencil gauges, and thermocouples were used. Additionally, the effect of surface tension on the particle size distribution was evaluated with 1-liter of distilled water and with the same including the addition of 0.5% soap/water mixture. Results from this effort indicated that the plume formed from the instantaneous release in all conditions primarily included two distinct particle size distributions including 217-312µm and 0.53µm on average. Soapy water formed a higher concentration of smaller mist (0.53µm) than distilled water, potentially due to decreased surface tension with the added surfactant. Though results showed good correlation to certain particle size distributions it is encouraged to conduct additional tests like this with a wider range of liquids with different properties, which could affect the overall particle size distribution when released.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Liquid Oxygen Liquid Acquisition Device Bubble Point Tests with High Pressure LOX at Elevated Temperatures

When transferring propellant in space, it is most efficient to transfer single phase liquid from a propellant tank to an engine. In earth s gravity field or under acceleration, propellant transfer is fairly simple. However, in low gravity, withdrawing single-phase fluid becomes a challenge. A variety of propellant management devices (PMD) are used to ensure single-phase flow. One type of PMD, a liquid acquisition device (LAD) takes advantage of capillary flow and surface tension to acquire liquid. The present work reports on testing with liquid oxygen (LOX) at elevated pressures (and thus temperatures) (maximum pressure 1724 kPa and maximum temperature 122K) as part of NASA s continuing cryogenic LAD development program. These tests evaluate LAD performance for LOX stored in higher pressure vessels that may be used in propellant systems using pressure fed engines. Test data shows a significant drop in LAD bubble point values at higher liquid temperatures, consistent with lower liquid surface tension at those temperatures. Test data also indicates that there are no first order effects of helium solubility in LOX on LAD bubble point prediction. Test results here extend the range of data for LOX fluid conditions, and provide insight into factors affecting predicting LAD bubble point pressures.

Jurns, John M.↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

Thermocapillary Motion Research

Fundamental microgravity research on interfacial surface tension driven motion by a linear temperature gradient was conducted. A vertical linear temperature gradient will be established within a test container filled with a host fluid. The thermal gradient will be positive in the upward vertical direction to avoid any free convective currents. With minimal disturbance, a single immisible fluid droplet will be introduced into the host fluid. As predicted by existing theory, the droplet will move in the direction of higher temperature due to tangential shear stresses created by surface tension variations along the periphery of the droplet.

Peterson, T. T.↗

Reduced-Gravity Measurements of the Effect of Oxygen on Properties of Zirconium

The influence of oxygen on the thermophysical properties of zirconium is being investigated using MSL-EML (Material Science Laboratory - Electromagnetic Levitator) on ISS (International Space Station) in collaboration with NASA, ESA (European Space Agency), and DLR (German Aerospace Center). Zirconium samples with different oxygen concentrations will be put into multiple melt cycles, during which the density, viscosity, surface tension, heat capacity, and electric conductivity will be measured at various undercooled temperatures. The facility check-up of MSL-EML and the first set of melting experiments have been successfully performed in 2015. The first zirconium sample will be tested near the end of 2015. As part of ground support activities, the thermophysical properties of zirconium and ZrO were measured using a ground-based electrostatic levitator located at the NASA Marshall Space Flight Center. The influence of oxygen on the measured surface tension was evaluated. The results of this research will serve as reference data for those measured in ISS.

Zhao, J.↗

Macrovoid Defect Growth during Evaporative Casting of Polymeric Membranes

Macrovoid (MV) formation is a significant problem in evaporatively cast polymeric membranes. MVs are large, elongated or teardrop-shaped pores (~10-50 micron) that can impair membrane structural integrity. Although MVs have been extensively studied, there is no general agreement on the mechanisms governing MV growth. Recently, our research group has formulated the solutocapillary convection (SC) hypothesis, which contends that MV growth involves three principal forces: a Marangoni force generated by surface tension gradients within the MV interface, a viscous drag force, and a gravitationally induced body force. Two sets of complementary experiments were conducted to test the SC hypothesis. Ground-based videomicroscopy flow-visualization (VMFV) was utilized to measure the flow velocities at the MV-casting solution interface and deep within the casting solution. The measurements were performed with casting solutions containing 10 wt% cellulose acetate (CA), 30 wt% H2O, 60 wt% acetone, and 200- ppm TiO2 particles for flow visualization, and the surface tension was controlled by surfactant addition. Qualitatively, the experiments indicated that MV growth occurs in three distinct phases: (1) a very rapid initial growth period, (2) a much slower growth phase, and (3) absorption of selected MVs into the expanding demixed region. The presence of tracer particles inside the MVs suggests the presence of a convective flow, which transfers the particles from the bulk solution to the MV interior. Although the VMFV experiments did not establish any surfactant effect on the interfacial velocities, a statistically significant effect on the MV number density was observed. In the second set of experiments, membranes were cast aboard a KC-135 aircraft under 0-g and 2-g conditions. Despite careful attention to the design and fabrication of the membrane casting apparatus (MCA), several problems were encountered, the most significant of which was the contamination of the casting solution by the activated carbon particles used for solvent absorption.

Greenberg, A. R.↗

Thermocapillary stabilization of the capillary breakup of an annular film of liquid

It is known that the breakup by surface tension of a cylindrical interface containing a viscous liquid can be damped by axial motion of the underlying liquid and that for an annular film the capillary instability can be completely suppressed (disturbances of all wavelengths decay) by certain axial velocity profiles. Here, using a linear stability analysis, it is shown that complete stabilization can also occur for thermocapillary-driven axial motions. However, the influence of thermocapillary instabilities typically shrinks the window in parameter space where stabilization is found, relative to the isothermal case. The influence of Reynolds, surface tension, Prandtl, and Biot parameters on limits of stabilization is calculated using continuation techniques. It is observed that windows of stabilization first open with topological changes of the neutral curves in parameter space.

Dijkstra, Henk A.↗

Solvability condition for needle crystals at large undercooling in a nonlocal model of solidification

It is explicitly shown that, in a realistic model of diffusion-controlled dendritic solidification, Ivantsov's continuous family of steady-state needle crystals is destroyed by the addition of surface tension. The starting point is in the exact integro-differential equation for the one-sided model, in two dimensions, in a moving frame of reference. In the limit of large undercooling, where the range of the diffusion field is much smaller than the radius of curvature of the tip of the needle, this problem is reduced to a linear, inhomogeneous differential equation of infinite order. A solvability condition for this equation is derived and it is shown that solutions cease to exist for arbitrarily small but finite isotropic surface tension.

Caroli, B.↗

Measurements of Thermophysical Properties of Molten Silicon by a High-Temperature Electrostatic Levitator

Several thermophysical properties of molten silicon measured by the high-temperature electrostatic levitator at JPL are presented. They are density, constant-pressure specific heat capacity, hemispherical total emissivity, and surface tension. Over the temperature range investigated (1350 less than T(sub m) less than 1825 K), the measured liquid density (in g/cc) can be expressed by a quadratic function, pi(P) = p(sub m) - 1.69 x 10(exp -4)(T - T(sub m)) - 1.75 x 10(exp -7)((T - T(sub m))(exp 2) with T(sub m) and p(sub m) being 1687 K and 2.56 g/cc, respectively. The hemispherical total emissivity of molten silicon at the melting temperature was determined to be 0.18, and the constant-pressure specific heat was evaluated as a function of temperature. The surface tension (in 10(exp -3) N/m) of molten silicon over a similar temperature range can be expressed by sigma(T) = 875 - 0.22( T - T(sub m)).

Rhim, W. K.↗

Bubble Formation on a Wall in Cross-Flowing Liquid and Surrounding Fluid Motion,With and Without Heating

Application of gas-liquid two-phase flow systems for space-based thermal management and for the HEDS program demands a precise control of bubble size distribution in liquid. The necessity of bulk liquid motion for controlling bubble size and frequency in the space environment has been suggested by recent studies on pool, forced convection boiling and bubble formation in flowing liquid. The present work, consisting of two parts, explores bubble generation at wall in a cross-flowing liquid, i.e., in a forced convection boiling configuration. A schematic is shown. The first part looks into the bubble formation process under isothermal conditions in a reduced gravity environment, by injecting gas through a hole in the wall of a flowing liquid channel. In the latter part with channel wall heating, flow and temperature fields near a single bubble are studied under normal (1-g) and micro-gravity (mu-g) conditions. The objective of the isothermal experiments is to experimentally investigate the effects of liquid cross-flow velocity, gas flow rate, and orifice diameter on bubble formation. Data were taken mainly under reduced gravity conditions but some data were taken in normal gravity for comparison. The reduced gravity experiment was conducted aboard the NASA DC-9 Reduced Gravity Aircraft. The results show that the process of bubble formation and detachment depends on gravity, the orifice diameter (D(sub N)), the gas flow rate (Q(sub g)), and the liquid cross-flow velocity (U(sub L)). The reduced gravity data are shown. The data are analyzed based on a force balance, and two different detachment mechanisms are identified. When the gas momentum is large, the bubble detaches from the injection orifice as the gas momentum overcomes the attaching effects of liquid drag and inertia. The surface tension force is much reduced because a large part of the bubble pinning edge at the orifice is lost as the bubble axis is tilted by the liquid flow. When the gas momentum is small, the force balance in the liquid flow direction is important, and the bubble detaches when the bubble axis inclination exceeds a certain angle. With wall heating, liquid motion around an air bubble in cross-flowing 2cs silicone oil is experimentally investigated in 1-g. A spectral element based steady 2D numerical model is also developed. The traces of particles from experimental flow visualization and the corresponding computed streamlines are shown. At the upstream side of the bubble facing the cross-flow, thermocapillary and forced convection create liquid motion away from the wall, up along the surface. At the downstream side, a competing interaction between the two creates a recirculation cell, causing the bulk liquid to stagnate on the surface and separate thereafter. The important dimensionless parameters are - Surface tension and local cross-flow Reynolds numbers R(sub sigma)(U(sub ref)a/v) and Re(sub loc)(U(sub L)a/v), respectively based on reference thermocapillary U(sub ref)(sigma(sub T)Delta.T/mu, Delta T=T(sub wall)-T(sub liquid)) and local cross-flow velocity U(sub L), Prandtl number P(sub r) and Grashoff number Gr(rho.g.beta.DeltaTa(sup 3)/mu.v). Variation of the stagnation point with R(sub sigma) and Re(sub loc) is shown. Also shown are good agreement between experimental and numerical results in 1-g. The computational model is extended to mu-g condition to investigate temperature and velocity on the bubble surface, stagnation and reattachment points of the recirculation cell and wall heat transfer. It is observed that wall heating significantly alters the flow field around the bubble and thus the forces acting on the bubble, which control its detachment. Thus a combination of heating and liquid cross-flow can be utilized to precisely control bubble formation in a mu-g environment.

Bhunia, Avijit↗

Pinching solutions of slender cylindrical jets

Simplified equations for slender jets are derived for a circular jet of one fluid flowing into an ambient second fluid, the flow being confined in a circular tank. Inviscid flows are studied which include both surface tension effects and Kelvin-Helmholtz instability. For slender jets a coupled nonlinear system of equations is found for the jet shape and the axial velocity jump across it. The equations can break down after a finite time and similarity solutions are constructed, and studied analytically and numerically. The break-ups found pertain to the jet pinching after a finite time, without violation of the slender jet ansatz. The system is conservative and admissible singular solutions are those which conserve the total energy, mass, and momentum. Such solutions are constructed analytically and numerically, and in the case of vortex sheets with no surface tension certain solutions are given in closed form.

Papageorgiou, Demetrios T.↗

Studies of Two-Phase Gas-Liquid Flow in Microgravity

Two-phase gas-liquid flows are expected to occur in many future space operations. Due to a lack of buoyancy in the microgravity environment, two-phase flows are known to behave differently than those in earth gravity. Despite these concerns, little research has been conducted on microgravity two-phase flow and the current understanding is poor. This dissertation describes an experimental and modeling study of the characteristics of two-phase flows in microgravity. An experiment was operated onboard NASA aircraft capable of producing short periods of microgravity. In addition to high speed photographs of the flows, electronic measurements of void fraction, liquid film thickness, bubble and wave velocity, pressure drop and wall shear stress were made for a wide range of liquid and gas flow rates. The effects of liquid viscosity, surface tension and tube diameter on the behavior of these flows were also assessed. From the data collected, maps showing the occurrence of various flow patterns as a function of gas and liquid flow rates were constructed. Earth gravity two-phase flow models were compared to the results of the microgravity experiments and in some cases modified. Models were developed to predict the transitions on the flow pattern maps. Three flow patterns, bubble, slug and annular flow, were observed in microgravity. These patterns were found to occur in distinct regions of the gas-liquid flow rate parameter space. The effect of liquid viscosity, surface tension and tube diameter on the location of the boundaries of these regions was small. Void fraction and Weber number transition criteria both produced reasonable transition models. Void fraction and bubble velocity for bubble and slug flows were found to be well described by the Drift-Flux model used to describe such flows in earth gravity. Pressure drop modeling by the homogeneous flow model was inconclusive for bubble and slug flows. Annular flows were found to be complex systems of ring-like waves and a substrate film. Pressure drop was best fitted with the Lockhart- Martinelli model. Force balances suggest that droplet entrainment may be a large component of the total pressure drop.

Bousman, William Scott↗

Similarity rules in gravity jitter-related spacecraft liquid propellant slosh waves excitation

The dynamical behavior of fluids, in particular the effect of surface tension on partially filled rotating fluids in a full-scale prototype Gravity Probe-B Spacecraft propellant tank and various 10 percent subscale containers with identical values of similarity parameters such as Bond number, dynamical capillary number, rotational Reynolds number, and Weber number, as well as imposed gravity jitters have been investigated. It is shown that the Bond number can be used to simulate the wave characteristics of slosh wave excitation, whereas the Weber number can be used to simulate the wave amplitude of slosh-mode excitation. It is shown that a dynamical capillary number can be used to simulate the induced perturbation of the fluid stress distribution exerted on the wall. This distribution is governed by the interaction between surface tension (slosh-wave excitation along the liquid-vapor interface) and viscous (fluid stress exerted on the wall) forces.

Hung, R. J.↗

Thermophysical Property Measurements of Molten Semiconductors in 1-g and Reduced-g Condition

Understanding and controlling the formation kinetics of varieties of crystal imperfections such as point defects, non uniform distribution of doping atoms, and impurity atoms in growing crystals are very important. Theoretical (numerical) modeling of the crystal growth process is an essential step to achieving these objectives. In order to obtain reliable modeling results, input parameters, i.e. various thermophysical parameters, must be accurate. The importance of accurate thermophysical properties of semiconductors in crystal growth cannot be overly emphasized. The total hemispherical emissivity, for instance, has a dramatic impact on the thermal environment. It determines the radiative emission from the surface of the melt which determines to a large extent the profile of the solidified crystal. In order to understand the convection and the turbulence in a melt, viscosity becomes an important parameter. The liquid surface tension determines the shape of the liquid-atmosphere interface near the solid-liquid-atmosphere triple point. Currently used values for these parameters are rather inaccurate, and this program intends to provide more reliable measurements of these thermophysical properties. Thus, the objective of this program is in the accurate measurements of various thermophysical properties which can be reliably used in the modeling of various crystal growth processes. In this program, thermophysical properties of molten semiconductors, such as Si, Ge, Si-Ge, and InSb will be measured as a function of temperature using the High Temperature Electrostatic Levitator at JPL. Each material will be doped by different kinds of impurities at various doping levels. Thermophysical properties which will be measured include: density, thermal expansion coefficient, surface tension, viscosity, specific heat, hemispherical total emissivity, and perhaps electrical and thermal conductivities. Many molten semiconductors are chemically reactive with crucibles. As a result, these dispersed impurities in the melts tend to substantially modify the properties of pure semiconductors. Sample levitation done in a vacuum clearly helps maintain the sample purity. However, in the 1-g environment, all gravity caused effects such as convection, sedimentation and buoyancy are still present in the sample. In addition, large forces needed to levitate a sample in the presence of the gravity can cause additional flows in the melt. The use of the High Temperature Electrostatic Levitator (HTESL) for the present research is a recent development and little is known about the flows induced by the electrostatic forces. In this ground base program, we will define the limits of HTESL technology as various thermophysical properties of molten semiconductors are measured.

Rhim, Won-Kyu↗

Measurements of the Viscosity of the Undercooled Melts Under the Conditions of Microgravity and Supporting MHD Calculations

This report covers the work to support the measurement of the surface tension and viscosity of undercooled metals with the TEMPUS as a part of the IML-2 mission. This work consisted of two parts. First was the ground-based research program, whose purpose was to establish the feasibility of the flight program and to establish the knowledge base necessary to plan a successful flight program. Analytical calculations established a fundamental understanding of the problem, and then a more rigorous program of numerical calculations was employed in the development of the experimental program. The second part of the covered research was the actual flight program, developed in close collaboration with Prof. Egry. A metal sphere was levitated and melted, and then set oscillating by 'squeezing' it with the magnetic field. The droplet's oscillations were recorded on videotape and digital image analysis was employed in reduction of the data. The experiment was successful because surface tension measurements were obtained for the gold and gold-copper alloy samples. This data provides unique experimental evidence in support of theories about electromagnetic levitation on the ground. Because of unforeseen difficulties with stability of liquid samples in TEMPUS, the viscosity measurements were not possible, but the data collection and analysis techniques were well proven in this mission. Further research is being conducted regarding the measurement of viscosity by this technique.

Szekely, Julian↗

Forced Oscillations of Supported Drops

Oscillations of supported liquid drops are the subject of wide scientific interest, with applications in areas as diverse as liquid-liquid extraction, synthesis of ceramic powders, growing of pure crystals in low gravity, and measurement of dynamic surface tension. In this research, axisymmetric forced oscillations of arbitrary amplitude of viscous liquid drops of fixed volume which are pendant from or sessile on a rod with a fixed or moving contact line and surrounded by an inviscid ambient gas are induced by moving the rod in the vertical direction sinusiodally in time. In this paper, a preliminary report is made on the computational analysis of the oscillations of supported drops that have 'clean' interfaces and whose contact lines remain fixed throughout their motions. The relative importance of forcing to damping can be increased by either increasing the amplitude of rod motion A or Reynolds number Re. It is shown that as the ratio of forcing to damping rises, for drops starting from an initial rest state a sharp increase in deformation can occur when they are forced to oscillate in the vicinity of their resonance frequencies, indicating the incipience of hysteresis. However, it is also shown that the existence of a second stable limit cycle and the occurrence of hysteresis can be observed if the drop is subjected to a so-called frequency sweep, where the forcing frequency is first increased and then decreased over a suitable range. Because the change in drop deformation response is abrupt in the vicinity of the forcing frequencies where hysteresis occurs, it should be possible to exploit the phenomenon to accurately measure the viscosity and surface tension of the drop liquid.

Wilkes, Edward D.↗

Surface deformation coupled with self-organized pattern on a liquid anode of 1 atm DC glow discharge

Intricate, self-organized plasma structures observed above the surface of a liquid anode of atmospheric DC glow discharge were found to give rise to coherent, organized surface deformation and mechanical wave formation at the plasma–liquid interface. This new phenomenon indicates that the liquid is closely coupled to the plasma by the anode sheath’s electrohydrodynamic (EHD) force. A scientific question then arises: Do surface perturbations, coupled with the nonuniform surface charge distribution, enhance the electric field and induce self-organization? Using the reflective background-oriented schlieren technique, the liquid surface profile under the plasma pattern was measured for the first time. The results show that surface distortions are driven by the repulsive Coulomb force of nonuniform net-negative surface charge acted by the anode sheath field. The impacts of various operating parameters on the patterns and surface waves were examined, revealing the significance of gas heating and liquid charge relaxation time in the pattern formation mechanism. Time-resolved dynamics of a pulsed DC discharge indicated that the surface deformation only appeared after the establishment of plasma patterns. Statistically, the surface wave under the plasma has high wave numbers (8000–16000 m −1 ) and small amplitudes ($<$ 10 µm), generally found in the capillary wave regime. Yet the motion of surface deformations is in tandem with the plasma pattern and exhibits a nondispersive nature of constant phase velocity (0.1–0.4 m s −1 ), suggesting the dominant role of EHD force over the surface tension in the observed surface wave. These results indicate that the deformed liquid surface is driven by the EHD force of the plasma sheath. Although they share a similar geometry, the deformation and wave dynamics of the liquid surface do not stimulate a plasma pattern. Importantly, the complex EHD coupling in the plasma–liquid system raises awareness and new challenges for plasma control engineering, and the quantitative characteristics of the nondispersive surface wave are informative for advancing relevant theory and modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Convective effects in float-zone and Czochralski melts

The hydrodynamics of crystal-growth melts is a relatively new research area. Numerical modeling of these processes is necessary. The work discussed herein is in two parts: numerical simulations of the flow in a Czochralski melt, and also of that in a float zone. In addition, for the float-zone case, energy stability theory will be used to determine stability bounds for the onset of oscillatory thermo-capillary flow. Convective effects in crystal-growth melts arise from a variety of mechanisms. Temperature gradients both in the direction of gravity and normal to it give rise to convection due to buoyancy effects. Rotation of the crucible and/or crystal causes a forced convection which may augment or oppose the buoyancy-driven flow. Finally, thermo-capillary forces (due to the variation of surface tension with temperature) drive surface motions which in turn generate convection in the bulk fluid. All of these mechanisms are present in either Czochralski or float-zone growth. The objective of the Czochralski modeling is to develop an accurate numerical simulation of the flow in a Czochralski silicon melt and to investigate the effects of various parameters on the flow properties. Like some earlier investigations, the intent is to simulate the effects of buoyancy, forced and thermo-capillary convection, including unsteady effects. Unlike earlier work, the aim is to include the effects of a variable free surface and freezing interface and, possibly incorporate nonaxisymmetric effects.

Neitzel, G. P.↗