Search NASASearch

Engineering topics

Alexander, J. I. D.

Publications and source records attributed to Alexander, J. I. D..

Investigation of Vibrational Control of the Bridgman Crystal Growth Technique

The character of natural buoyant convection in rigidly contained inhomogeneous fluids can be drastically altered by vibrating the container. Vibrations are expected to play a crucial influence on heat and mass transfer onboard the International Space Station (ISS). It is becoming evident that substantial vibrations will exist on the ISS in the wide frequency spectrum. In general, vibrational flows are very complex and governed by many parameters. In many terrestrial crystal growth situations, convective transport of heat and constituent components is dominated by buoyancy driven convection arising from compositional and thermal gradients. Thus, it may be concluded that vibro-convective flow can potentially be used to influence and even control transport in some crystal growth situations.

Fedoseyev, Alexandre I.

Orbital Processing of High-Quality Zn-Alloyed CdTe Compound Semiconductors

The objective of this research is to investigate the influences of gravitationally-dependent phenomena (hydrostatic and buoyant) on the growth and quality of doped and alloyed Cadmium-Zinc-Telluride (CdZnTe) crystals grown by the modified seeded Bridgman-Stockbarger technique. It is hypothesized that the damping of the gravitationally-dependent buoyancy convection will substantially enhance chemical homogeneity and the near-elimination of hydrostatic pressure will enable significant reduction in defect (dislocations and twins) density.

Larson, David J., Jr.

Experiment 1: Orbital Processing of High-Quality Zn-Alloyed CdTe Compound Semiconductors

Zinc-alloyed cadmium telluride (Zn:CdTe) crystals were grown in unit gravity and in microgravity for comparative analysis. Two crystals were grown on USML-1/STS-50 in 1992, and two additional crystals were grown on USML-2/STS-73 in 1995. The Crystal Growth Furnace (CGF) in the seeded Bridgman-Stockbarger crystal growth geometry was utilized on both missions. Crystals grown on USML-1/STS-50 were found to have solidified with partial wall contact due to the near-absence of the hydrostatic pressure in microgravity, a residual g-vector that was not axial, and the nonwetting sample/ampoule wetting conditions. Crystals grown on USML-2/STS-73 included: a sample/ampoule identical to the USML-1/STS-73 sample/ampoule, with the addition of a restraining spring to simulate hydrostatic pressure internally, and a tapered ampoule which accomplished 2.2 cm of crystal growth without wall contact. Infrared transmission of all ground and flight samples was found to be 63-66 percent, very close to the theoretical 66 percent, suggesting good stoichiometric control. Infrared microscopy confirmed that the primary precipitates were Te and their size (1-10 microns) and density suggested that the flight- and ground-base samples experienced similar thermal histories. Longitudinal macrosegregation, calculated using scaling analysis, was predicted to be low. Nearly diffusion controlled growth was achieved even in unit gravity and macrosegregation data could be fit with a diffusion controlled model. Radial segregation was monitored and was found to vary with fraction solidified, particularly through the shoulder region, where the sample cross section was varying significantly. It was also disturbed in the flight samples in regions where asymmetric wall contact was noted. In regions where a steady state was established, the radial segregation was invariant within our experimental measurement error. Flight samples were found to be much higher in structural perfection than samples processed in unit gravity under identical growth conditions. In regions where solidification had occurred without wall contact, the free surfaces evidenced virtually no twinning, although twins appeared in the flight samples in regions of wall contact and were pervasive in the ground samples. These results were confirmed using optical microscopy and synchrotron x-ray white beam topography. Full-width half-maximum rocking curve widths (FWHM), recorded in arc-seconds, were significantly reduced from 20 a-s (1 g) to 9 a-s (micro-g) for the best regions of the crystals. The 9 a-s (FWHM) rocking curve value in the unconfined flight samples equals the best value reported terrestrially for this material. The ground samples exhibited a fully developed (111)[110] dislocation mosaic structure, whereas dislocations within the flight samples were discrete and no mosaic structure was evident. The defect density was quantitatively reduced from 75,000 (1 g) to 800 (micro-g) +/-50 percent. Dislocation etch pit density results were confirmed using transmission synchrotron white beam and monochromated beam topography. The low defect density is thought to have resulted from the near absence of hydrostatic pressure, which allowed the molten boule to solidify with little or no wall contact. This minimized the transfer of hoop stresses during solidification and post-solidification processing.

Larson, D. J., Jr.

Orbital processing of high-quality CdTe compound semiconductors

CdZnTe crystals were grown in one-g and in micro-g for comparative analysis. The two micro-g crystals were grown in the Crystal Growth Furnace during the First United States Microgravity Laboratory mission (USML-1). The samples were analyzed for chemical homogeneity, structural perfection, and optoelectronic performance (infrared transmission). Fourier Transform Infrared (FTIR) transmission of both ground and flight materials showed that the infrared transmission was close to theoretical, 63% versus 66%, suggesting that the material was close to the stochiometric composition during both the ground and flight experiments. Infrared microscopy confirmed that the principal precipitates were Te and their size (1-10 microns) and density suggested that the primary flight and ground base samples experienced similar cooling rates. Macrosegregation was predicted, using scaling analysis, to be low even in one-g crystals and this was confirmed experimentally, with nearly diffusion controlled growth achieved even in the partial mixing regime on the ground. Radial segregation was monitored in the flight samples and was found to vary with fraction solidified, but was disturbed due to the asymmetric grvitational and thermal fields experienced by the flight samples. The flight samples, however, were found to be much higher in structural perfection than the ground samples produced in the same furnace under identical growth conditions except for the gravitational level. Rocking curve widths were found to be substantially reduced, from 20/35 (one-g) to 9/20 (micro-g) for the best regions of the crystals. The full width at half maximum (FWHM) of 9 arc seconds is as good as the best reported terrestrially for this material. The ground samples were found to have a fully developed mosaic structure consisting of subgrains, whereas the flight sample dislocations were discrete and no mosaic substructure was evident. The defect density was reduced from 50-100,000 (one-g) to 500-25000 EPD (micro-g). These results were confirmed using rocking curve analysis, synchrotron topography, and etch pit analysis. The low dislocation density is thought to have resulted from the near-absence of hydrostatic pressure which allowed the melt to solidify with minimum or no wall contact, resulting in very low stress being exerted on the crystal during growth or during post-solidification cooling.

Larson, David J., Jr.

A Chebyshev Collocation Method for Moving Boundaries, Heat Transfer, and Convection During Directional Solidification

Free and moving boundary problems require the simultaneous solution of unknown field variables and the boundaries of the domains on which these variables are defined. There are many technologically important processes that lead to moving boundary problems associated with fluid surfaces and solid-fluid boundaries. These include crystal growth, metal alloy and glass solidification, melting and name propagation. The directional solidification of semi-conductor crystals by the Bridgman-Stockbarger method is a typical example of such a complex process. A numerical model of this growth method must solve the appropriate heat, mass and momentum transfer equations and determine the location of the melt-solid interface. In this work, a Chebyshev pseudospectra collocation method is adapted to the problem of directional solidification. Implementation involves a solution algorithm that combines domain decomposition, finite-difference preconditioned conjugate minimum residual method and a Picard type iterative scheme.

Zhang, Yiqiang

Detailed analysis of Honeywell In-Space Accelerometer data - STS-32

The Honeywell In-Space Accelerometer (HISA) system collected data in the mid-deck area of the Shuttle Columbia during the flight of STS-32, January 1990. The resulting data were to be used to investigate the response of crystal microstructure to different types of residual acceleration. The HISA is designed to detect and record transient and oscillatory accelerations. The sampling and electronics package stored averaged accelerations over two sampling periods; two sampling rates were available: 1 Hz and 50 Hz. Analysis of the HISA data followed the CMMR Acceleration Data Processing Guide, considering in-house computer modelling of a float-zone indium crystal growth experiment. Characteristic examples of HISA data showing the response to the primary reaction control system, Orbiter Maneuvering System operations, and crew treadmill activity are presented. Various orbiter structural modes are excited by these and other activities.

Rogers, Melissa J. B.

Numerical simulation of thermocapillary flow under zero and low gravity conditions

This paper discusses the numerical solution methods and results of steady and unsteady thermocapillary (surface-tension) and buoyancy driven flows in 2D cavities and liquid columns. The 2D cavity was assumed to be square with one free surface with a zero Capillary number (i.e., the free surface was constrained to be flat). A pseudospectral method was used to solve steady and unsteady surface tension-driven and mixed buoyancy-surface tension flows in a square cavity. For the liquid column a finite-difference scheme based on a Picard iteration was used to solve for the flow, temperature and free surface shape. The surface of the liquid column was allowed to deform and, as for the 2D cavity, the surface tension was assumed to depend on temperature.

Alexander, J. I. D.

Effects of buoyancy-driven flow and thermal boundary conditions on physical vapor transport

A 2D numerical model was developed in order to ascertain if reduced gravity conditions are beneficial to physical vapor transport (PVT) and to determine its tolerance limits to residual accelerations. This was solved using the PHOENICS finite-volume code. Reduction of gravitational accelerations to less than 0.1 g0 was found to be sufficient to suppress buoyancy-driven convection to an extent that diffusion was the dominant transport mode, whence a greater uniformity in the growth rate could be obtained. It is shown that a uniform temperature gradient on the ampoule walls causes the vapor to be supersaturated throughout the ampoule, potentially resulting in undesirable nucleation at the walls. A 'hump' in the wall temperature profile can be used to avoid this. The prevailing transport conditions determine the size of the hump needed.

Nadarajah, Arunan

Residual acceleration data analysis for Spacelab missions

A two-level reduction plan is developed analytically with references to existing experimental data that can be used to develop accelerometer databases for postflight experiment analyses. Residual-acceleration data from previous experiments are studied with attention given to the differences in magnitude among the databases. First-level processing is based on the identification of instances in which tolerable acceleration levels are exceeded so that the time intervals requiring further scrutiny are found. Acceleration time histories including magnitude and orientation are utilized with frequency data for the second level of processing, and the data required for a given time window are reduced by means of data decimation. The identification and processing of limited windows of acceleration data facilitates the characterization of the acceleration environment. These characterizations are of use in the design, location, and utilization of low-gravity laboratory equipment.

Rogers, Melissa J. B.

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang

Experiment specific processing of residual acceleration data

To date, most Spacelab residual acceleration data collection projects have resulted in data bases that are overwhelming to the investigator of low-gravity experiments. This paper introduces a simple passive accelerometer system to measure low-frequency accelerations. Model responses for experiments using actual acceleration data are produced and correlations are made between experiment response and the accelerometer time history in order to test the idea that recorded acceleration data and experimental responses can be usefully correlated. Spacelab 3 accelerometer data are used as input to a variety of experiment models, and sensitivity limits are obtained for particular experiment classes. The modeling results are being used to create experiment-specific residual acceleration data processing schemes for interested investigators.

Rogers, Melissa J. B.

The sensitivity of a non-isothermal liquid bridge to residual acceleration

An examination is conducted of the steady and unsteady behavior of isothermal and nonisothermal liquid bridge systems under both normal and low-gravity conditions, recasting the full nonlinear equations in terms of a stream-function vorticity formulation in conjunction a nonorthogonal coordinate transformation. For the cases examined, the shape of the bridge is more sensitive to typical spacecraft accelerations than the buoyancy-driven flow. The effect of thermocapillary flow on the surface shape is found to be small, for the range of capillary and Reynolds numbers considered.

Alexander, J. I. D.

Numerical analysis of the sensitivity of crystal growth experiments to spacecraft residual acceleration

An analysis is conducted of the sensitivity of the Bridgman-Stockbarger crystal growth method, using an idealized model for a range of operating and boundary conditions over a variety of accelerations. Attention is given to the dopant nonuniformity at the melt-crystal interface. The largest compositional nonuniformities are found to occur for disturbances whose amplitudes are greater than 10 exp 6 g, and frequencies below 0.1 Hz.

Alexander, J. I. D.

Analysis of Spacelab 3 residual acceleration data

A data reduction plan is being developed to efficiently process residual acceleration data, making such data more accessible to principal investigators of low-gravity experiments. Accelerometer data collected during the Spacelab 3 mission is being processed using a prototype version of this plan. The plan addresses various aspects of acceleration data analysis: the identification of disturbances that are intolerable to experiments, the investigation of acceleration orientations, the definition and characterization of the background acceleration corresponding to a given experiment time line, the isolation and examination of particular significant disturbances, and the identification of disturbance sources. Acceleration magnitude, frequency, and orientation are discussed: transient accelerations can have magnitudes as large as 0.01 g with frequency components no greater than 0.001 g. These accelerations fluctuate rapidly in orientation. The occurrence of disturbance sources in an orbiter is tentatively identified as a random process, whereas the response of the orbiter to given accelerations is considered deterministic. The need to continue monitoring of residual accelerations in orbiting space laboratories is stressed.

Rogers, Melissa J. B.

A strategy for residual acceleration data reduction and dissemination

A data reduction plan is being developed to efficiently process residual acceleration data from orbiting space laboratories. Implementation of the reduction plan will result in a useful, manageable accelerometer data base which can be readily employed by principal investigators during post-flight analysis of experimental results. The data reduction plan will also assist in the characterization of the acceleration environment of orbiters, which is important for the planning of future experimental missions. Prototype versions of the reduction plan are being tested using accelerometer data recorded during the Spacelab 3 (SL3) mission. Transient disturbances caused by shuttle and experiment operations and crew activities are being cataloged. The magnitudes of such disturbances can be as much as 0.01 g, but are rarely sustained for more than a fraction of a second and tend to vary greatly in orientation. It is found that the recorded frequency components of the SL3 disturbances do not exceed tolerance limits determined from modeling for a selected experiment set.

Rogers, Melissa J. B.

Influence of residual gravity on crystal growth processes

The effects of steady and time-dependent residual accelerations on representative melt and solution crystal growth systems have been investigated using numerical models. For the melt growth system, it was found that the lateral nonuniformity in crystal composition, depending on the particular operating conditions, can be very sensitive to the orientation of the steady component of the residual gravity vector down to 10 to the -6th g. Oscillatory residual accelerations (g-jitter) were found to lead to transients in the dopant distribution. For the solution growth, longtime simulations showed that at 10 to the -5th g, and the cooling rates used in Spacelab 3 experiments, diffusive transport dominates only in the beginning. Convective transport becomes increasingly significant after 3 to 5 hours, depending on the orientation of the crystal in the residual gravity field.

Rosenberger, F.

Bridgman crystal growth in low gravity - A scaling analysis

The results of an order-of-magnitude or scaling analysis are compared with those of numerical simulations of the effects of steady low gravity on compositional nonuniformity in crystals grown by the Bridgman-Stockbarger technique. In particular, the results are examined of numerical simulations of the effect of steady residual acceleration on the transport of solute in a gallium-doped germanium melt during directional solidification under low-gravity conditions. The results are interpreted in terms of the relevant dimensionless groups associated with the process, and scaling techniques are evaluated by comparing their predictions with the numerical results. It is demonstrated that, when convective transport is comparable with diffusive transport, some specific knowledge of the behavior of the system is required before scaling arguments can be used to make reasonable predictions.

Alexander, J. I. D.