Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffused interface model”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Particle Engulfment and Pushing By Solidifying Interfaces

The study of particle behavior at solid/liquid interfaces (SLI s) is at the center of the Particle Engulfment and Pushing (PEP) research program. Interactions of particles with SLI s have been of interest since the 1960 s, starting with geological observations, i.e., frost heaving. Ever since, this field of research has become significant to such diverse areas as metal matrix composite materials, fabrication of superconductors, and inclusion control in steels. The PEP research effort is geared towards understanding the fundamental physics of the interaction between particles and a planar SLI. Experimental work including 1-g and mu-g experiments accompany the development of analytical and numerical models. The experimental work comprised of substantial groundwork with aluminum (Al) and zinc (Zn) matrices containing spherical zirconia particles, mu-g experiments with metallic Al matrices and the use of transparent organic metal-analogue materials. The modeling efforts have grown from the initial steady-state analytical model to dynamic models, accounting for the initial acceleration of a particle at rest by an advancing SLI. To gain a more comprehensive understanding, numerical models were developed to account for the influence of the thermal and solutal field. Current efforts are geared towards coupling the diffusive 2-D front tracking model with a fluid flow model to account for differences in the physics of interaction between 1-g and -g environments. A significant amount of this theoretical investigation has been and is being performed by co-investigators at NASA MSFC.

Source record↗
NASA NTRS2014

Behavior of Particle Depots in Molten Silicon During Float-Zone Growth in Strong Static Magnetic Fields

In the photovoltaics industry, the largest market share is represented by solar cells made from multicrystalline silicon, which is grown by directional solidification. During the growth process, the silicon melt is in contact with the silicon nitride coated crucible walls and the furnace atmosphere which contains carbon monoxide. The dissolution of the crucible coating, the carbon bearing gas, and the carbon already present in the feedstock, lead to the precipitation of silicon carbide, and silicon nitride, at later stages of the growth process. The precipitation of Si3N4 and SiC particles of up to several hundred micrometers in diameter leads to severe problems during the wire sawing process for wafering the ingots. Furthermore the growth of the silicon grains can be negatively influenced by the presence of particles, which act as nucleation sources and lead to a grit structure of small grains and are sources for dislocations. If doped with Nitrogen from the dissolved crucible coating, SiC is a semi conductive material, and can act as a shunt, short circuiting parts of the solar cell. For these reasons, the incorporation of such particles needs to be avoided. In this contribution we performed model experiments in which the transport of intentionally added SiC particles and their interaction with the solid-liquid interface during float zone growth of silicon in strong steady magnetic fields was investigated. SiC particles of 7μm and 60μm size are placed in single crystal silicon [100] and [111] rods of 8mm diameter. This is achieved by drilling a hole of 2mm diameter, filling in the particles and closing the hole by melting the surface of the rod until a film of silicon covers the hole. The samples are processed under a vacuum of 1x10(exp -5) mbar or better, to prevent gas inclusions. An oxide layer to suppress Marangoni convection is applied by wet oxidation. Experiments without and with static magnetic field are carried out to investigate the influence of melt convection on the distribution of particles and their incorporation into the crystal. The field strengths applied by a superconducting magnet are 1T, 3T, 4.5T, and 5T. The increase in field strength dampens the melt flow, and so this study provides comparative data to the crystal growth experiment to be carried out onboard the sounding rocket mission TEXUS 51, where purely diffusive growth condition will be achieved under microgravity conditions.

Jauss, T.↗
NASA NTRS2001

Effect of Marangoni Convection Generated by Voids on Segregation During Low-G and 1-G Solidification

Solidification experiments, especially microgravity solidification experiments, are often compromised by the evolution of unwanted voids or bubbles in the melt. Although these voids and/or bubbles are highly undesirable, there is currently no effective means of preventing their formation or of eliminating their adverse effects, particularly during microgravity experiments. Marangoni convection caused by these voids can drastically change the transport processes in the melt. Recent microgravity experiments by Matthiesen (1) Andrews (2) and Fripp (3) are perfect examples of how voids and bubbles can affect the outcome of costly space experiments and significantly increase the level of difficulty in interpreting their results. Formation of bubbles have caused problems in microgravity experiments for a long time. Even in the early Skylab mission an unexpectedly large number of bubbles were detected in the four materials processing experiments reported by Papazian and Wilcox (4). They demonstrated that while during ground-based tests bubbles were seen to detach from the interface easily and float to the top of the melt, in low-gravity tests no detachment from the interface occurred and large voids were grown in the crystal. More recently, the lead-tin-telluride crystal growth experiment of Fripp et al.(3) flown aboard the USMP-3 mission has provided very interesting results. The purpose of the study was to investigate the effect of natural convection on the solidification process by growing the samples at different orientations with respect to the gravitational field. Large pores and voids were found in the three solid crystal samples processed in space. Post-growth characterization of the compositional profiles of the cells indicated considerable levels of mixing even in the sample grown in the hot-on-top stable configuration. The mixing was attributed to thermocapillary convection caused by the voids and bubbles which evolved during growth. Since the thermocapillary convection is orientation-independent, diffusion-controlled growth was not possible in any of the samples, even the top-heated one. These results are consistent with recent studies of thermocapillary convection generated by a bubble on a heated surface undertaken by Kassemi and Rashidnia (5-7) where it is numerically and experimentally shown that the thermocapillary flow generated by a bubble in a model fluid (silicone oil) can drastically modify the temperature field through vigorous mixing of the fluid around it, especially under microgravity conditions.

Kassemi, M.↗
NASA NTRS2010

NASA Tech Briefs, September 2010

Topics covered include: Instrument for Measuring Thermal Conductivity of Materials at Low Temperatures; Multi-Axis Accelerometer Calibration System; Pupil Alignment Measuring Technique and Alignment Reference for Instruments or Optical Systems; Autonomous System for Monitoring the Integrity of Composite Fan Housings; A Safe, Self-Calibrating, Wireless System for Measuring Volume of Any Fuel at Non-Horizontal Orientation; Adaptation of the Camera Link Interface for Flight-Instrument Applications; High-Performance CCSDS Encapsulation Service Implementation in FPGA; High-Performance CCSDS AOS Protocol Implementation in FPGA; Advanced Flip Chips in Extreme Temperature Environments; Diffuse-Illumination Systems for Growing Plants; Microwave Plasma Hydrogen Recovery System; Producing Hydrogen by Plasma Pyrolysis of Methane; Self-Deployable Membrane Structures; Reactivation of a Tin-Oxide-Containing Catalys; Functionalization of Single-Wall Carbon Nanotubes by Photo-Oxidation; Miniature Piezoelectric Macro-Mass Balance; Acoustic Liner for Turbomachinery Applications; Metering Gas Strut for Separating Rocket Stages; Large-Flow-Area Flow-Selective Liquid/Gas Separator; Counterflowing Jet Subsystem Design; Water Tank with Capillary Air/Liquid Separation; True Shear Parallel Plate Viscometer; Focusing Diffraction Grating Element with Aberration Control; Universal Millimeter-Wave Radar Front End; Mode Selection for a Single-Frequency Fiber Laser; Qualification and Selection of Flight Diode Lasers for Space Applications; Plenoptic Imager for Automated Surface Navigation; Maglev Facility for Simulating Variable Gravity; Hybrid AlGaN-SiC Avalanche Photodiode for Deep-UV Photon Detection; High-Speed Operation of Interband Cascade Lasers; 3D GeoWall Analysis System for Shuttle External Tank Foreign Object Debris Events; Charge-Spot Model for Electrostatic Forces in Simulation of Fine Particulates; Hidden Statistics Approach to Quantum Simulations; Reconstituted Three-Dimensional Interactive Imaging; Determining Atmospheric-Density Profile of Titan; Digital Microfluidics Sample Analyzer; Radiation Protection Using Carbon Nanotube Derivatives; Process to Selectively Distinguish Viable from Non-Viable Bacterial Cells; and TEAMS Model Analyzer.

Source record↗
NASA NTRS2004

Notes on Interface Growth Kinetics 50 Years after Burton, Cabrera and Frank

This overview is devoted to some unresolved basic problems in crystal growth kinetics. The density wave approach to propagation of a spatially diffuse interface between a growing crystal and its simple (e.g., metallic) melt is discussed is Section 2. This approach allows for the calculation of kinetic coefficients and is an alternative to the localized interface concept in which each atom belongs to either a solid or a liquid. Sections 3 and 4 deals mainly with layer growth from solution. Mutual retardation of the growth steps via their bulk and d a c e diffusion fields is the major subject. The influence of solution flow on step bunching (Section 4) suggests the essential influence of bulk diffusion on the surface morphology. The flow within the solution boundary layer enhances step-step interaction, influences the step bunching process and the resulting step pattern morphology on the growing surface. Recent experiments on the rates at which strongly polygonized steps on protein and small molecule crystals propagate during growth from solution are analyzed in Section 5 . We have shown that the step segments may be "singular" and that "one-dimensional nucleation" may be the rate limiting stage for the segments that are shorter or comparable in length to the thermodynamically equilibrium interlink distance. In this case, the reciprocal dependence of the segment propagation rate on the segment length that follow from the Gibbs-Thompson law, should be replaced by an abrupt switch from zero to a finite constant velocity. Until recently, the Kossel crystal remained the only model used in crystal growth theory. In such Kossel Gibbs-Thomson law, should be replaced by an abrupt switch &om zero to a finite constant velocity. crystals, all kinks at the steps are identical and the kink rate is a linear function of the supersaturation. In the non-Kossel crystals, there may be several kink configurations characterized by different geometries and energies. These configurations should appear in a specific sequence when each new lattice unit cell is filled. As a result of such a cooperative interaction within the unit cell, a non-linear dependence of the kink rate on the vapor pressure or solution concentration in excess over the equilibrium value should be expected.

Chernov, A. A.↗
NASA NTRS2004

The formation of spikes in the displacement of miscible fluids

We report on experiments in which a more viscous fluid displaces a less viscous one in a vertical cylindrical tube. These experiments were performed using silicone oils in a vertical pipette of small diameter. The more viscous fluid also had a slightly larger density than the less viscous fluid. In the initial configuration, the fluids were at rest, and the interface was nominally flat. A dye was added to the more viscous fluid for ease of observation of the interface between the fluids. The flow was initiated by pumping the more viscous fluid into the less viscous one. The displacement velocity was such that the Reynolds number was smaller than unity and the Peclet number for mass transfer between the fluids was large compared to unity. For upward displacement of the more viscous fluid from an initially stable configuration, an axisymmetric finger was observed under all conditions. However, a needle-shaped spike was seen to propagate from the main finger in many cases, similar to that observed by Petitjeans and Maxworthy for the displacement of a more viscous fluid by a less viscous one.

Rheology/methods↗
NASA NTRS1999

Crystal Growth Using MEPHISTO

The shuttle flight experiment "In Situ Monitoring of Crystal Growth Using MEPHISTO" was accomplished during STS-87 as part of the fourth flight of the United States Microgravity Payload (USMP-4), which was flown from November 19 to December 5, 1997. The data returned from that flight are just now beginning to yield quantitative results. This project is an international collaboration: the furnace system known as MEPHISTO was built in France by CNES (French National Space Agency) and CEA (French Atomic Energy Commission); the principal investigator, Prof. Reza Abbaschian, is from the University of Florida at Gainesville; and numerical and analytical modeling support includes collaborators from the University of New South Wales, Australia, the University of Wisconsin at Milwaukee, the National Institute of Standards and Technology, and the NASA Lewis Research Center. MEPHISTO is a French acronym that translates into English as Materials for the Study of Interesting Phenomena of Solidification on Earth and in Orbit. Since this was the fourth flight of the MEPHISTO furnace, the experiment is referred to as MEPHISTO-4. MEPHISTO-4 was a directional solidification experiment that studied the liquid-to-solid transformation of bismuth alloyed with tin. Directional solidification is a freezing technique common to the processing of the electronic materials used in integrated circuits and detectors, such as silicon and germanium. When liquids are frozen on Earth, they must be cooled. The cooling causes stirring because of density variations in the liquid. This stirring, known as natural convection, influences the quality of the resulting solid. During freezing, regions of high and low concentrations of tin are created. This introduces another important phenomenon: diffusion, or the movement by molecular action of matter from regions of high concentration to regions of lower concentration. In MEPHISTO-4, it is tin that diffuses from the high-concentration region in front of the solid-liquid interface to more distant low-concentration regions.

deGroh, Henry C., III↗
NASA NTRS1999

Chain Dynamics in Magnetorheological Suspensions

Magnetorheological (MR) suspensions are composed of colloidal particles which acquire dipole moments when subjected to an external magnetic field. At sufficient field strengths and concentrations, the dipolar particles rapidly aggregate to form long chains. Subsequent lateral cross-linking of the dipolar chains is responsible for a rapid liquid-to-solid-like rheological transition. The unique, magnetically-activated rheological properties of MR suspensions make them ideal for interfacing mechanical systems to electronic controls. Additionally, the ability to experimentally probe colloidal suspensions interacting through tunable anisotropic potentials is of fundamental interest. Our current experimental work has focused on understanding the fluctuations of dipolar chains. It has been proposed by Halsey and Toor (HT) that the strong Landau-Peierls thermal fluctuations of dipolar chains could be responsible for long-range attractions between chains. Such interactions will govern the long-time relaxation of MR suspensions. We have synthesized monodisperse neutrally buoyant MR suspensions by density matching stabilized ferrofluid emulsion droplets with D2O. This allows us to probe the dynamics of the dipolar chains using light scattering without gravitational, interfacial, and polydispersity effects to resolve the short-wavelength dynamics of the dipolar chains. We used diffusing wave spectroscopy to measure these dynamics. The particle displacements at short times that show an independence to the field strength, but at long times exhibit a constrained, sub-diffusive motion that slows as the dipole strength is increased. The experiments are in good qualitative agreement with Brownian dynamics simulations of dipolar chains. Although there have been several important and detailed studies of the structure and interactions in MR suspensions, there has not been conclusive evidence that supports or contradicts the HT model prediction that long-range interactions exist between fluctuating chains of dipolar particles. Resolving this issue would contribute greatly to the understanding of these interesting and important materials. We have begun to test the predictions of the HT model by both examining the dynamics of individual chains and by measuring the forces between dipolar chains directly to accurately and quantitatively assess the interactions that they experience. To do so, we employ optical trapping techniques and video-microscopy to manipulate and observe our samples on the microscopic level. With these techniques, it is possible to observe chains that are fluctuating freely in three-dimensions, independent of interfacial effects. More importantly, we are able to controllably observe the interactions of two chains at various separations to measure the force-distance profile. The techniques also allow us to study the mechanical properties of individual chains and chain clusters. Our work to this point has focused on reversibly-formed dipolar chains due to field induced dipoles where the combination of this chaining, the dipolar forces, and the hydrodynamic interactions that dictate the rheology of the suspensions. One can envision, however, many situations where optical, electronic, or rheological behavior may be optimized with magneto-responsive anisotropic particles. Chains of polarizable particles may have the best properties as they can coil and flex in the absence of a field and stiffen and orient when a field is applied. We have recently demonstrated a synthesis of stable, permanent paramagnetic chains by both covalently and physically linking paramagnetic colloidal particles. The method employed allows us to create monodisperse chains of controlled length. We observed the stability, field-alignment, and rigidity of this new class of materials. The chains may exhibit unique rheological properties in an applied magnetic field over isotropic suspensions of paramagnetic particles. They are also useful rheological models as bead-spring systems. These chains form the basis for our current experiments with optical traps.

Gast, A. P.↗
NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗
NASA NTRS1999

Role of Dynamic Nucleation at Moving Boundaries in Phase and Microstructure Selection

Solidification microstructures that form under steady-state growth conditions (cells, dendrites, regular eutectics, etc.) are reasonably well understood in comparison to other, more complex microstructures, which form under intrinsically non-steady-state growth conditions due to the competition between the nucleation and growth of several phases. Some important practical examples in this latter class include microstructures forming in peritectic systems in highly undercooled droplets, and in strip cast stainless steels. Prediction of phase and microstructure selection in these systems has been traditionally based on (1) heterogeneous nucleation on a static interface, and (2) comparing the relative growth rate of different phase/microstructures under steady-state growth conditions. The formation of new phases, however, occurs via nucleation on, or ahead of, a moving boundary. In addition, the actual selection process is controlled by a complex interaction between the nucleation process and the growth competition between the nuclei and the pre-existing phase under non-steady-state conditions. As a result, it is often difficult to predict which microstructure will form and which phases will be selected under prescribed processing conditions. This research addresses this critical role of nucleation at moving boundaries in the selection of phases and solidification microstructures through quantitative experiments and numerical modeling in peritectic systems. In order to create a well characterized system in which to study this problem, we focus on the directional solidification of hypo- and hyper-peritectic alloys in the two-phase region, imposing a large enough ratio of temperature gradient/growth rate (G/V(sub p)) to suppress the morphological instability of both the parent (alpha) and peritectic (Beta) phases, i.e. each phase alone would grow as a planar front. Our combined experimental and theoretical results show that, already in this simplified case, the growth competition of these two phases leads to a rich variety of microstructures that depend sensitively upon the relative importance of nucleation, diffusion, and convection.

Karma, Alain↗
NASA NTRS1999

Magnetic Damping of Solid Solution Semiconductor Alloys

The objective of this study is to: (1) experimentally test the validity of the modeling predictions applicable to the magnetic damping of convective flows in electrically conductive melts as this applies to the bulk growth of solid solution semiconducting materials; and (2) assess the effectiveness of steady magnetic fields in reducing the fluid flows occurring in these materials during processing. To achieve the objectives of this investigation, we are carrying out a comprehensive program in the Bridgman and floating-zone configurations using the solid solution alloy system Ge-Si. This alloy system has been studied extensively in environments that have not simultaneously included both low gravity and an applied magnetic field. Also, all compositions have a high electrical conductivity, and the materials parameters permit reasonable growth rates. An important supporting investigation is determining the role, if any, that thermoelectromagnetic convection (TEMC) plays during growth of these materials in a magnetic field. TEMC has significant implications for the deployment of a Magnetic Damping Furnace in space. This effect will be especially important in solid solutions where the growth interface is, in general, neither isothermal nor isoconcentrational. It could be important in single melting point materials, also, if faceting takes place producing a non-isothermal interface. In conclusion, magnetic fields up to 5 Tesla are sufficient to eliminate time-dependent convection in silicon floating zones and possibly Bridgman growth of Ge-Si alloys. In both cases, steady convection appears to be more significant for mass transport than diffusion, even at 5 Tesla in the geometries used here. These results are corroborated in both growth configurations by calculations.

Szofran, Frank R.↗
NASA NTRS2003

An Alternative Hypothesis for How Microgravity Improves Macromolecular Crystal Quality

There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests that for many proteins crystal nucleation and growth is by addition of associated species that are preformed by reversible concentration-driven self association processes in the bulk solution. We have developed a self-association model for the crystal nucleation and growth of the protein chicken egg lysozyme. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small an octamer in the 43 helix configuration (the proposed average sized growth unit) would have a M.W. approx. 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm (a 24-mer) commonly attach to the crystal. AFM results from Weichmann et al. (Ultramicroscopy 86, 159-166, 2001) suggest that associated species of up to 40-mers in size add to the (101) faces. These measurements reflect the sizes of units that both added and desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. On Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species is more likely to be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.

Pusey, Marc↗
NASA NTRS2001

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗
NASA NTRS2002

Tetragonal Lysozyme Nucleation and Crystal Growth: The Role of the Solution Phase

Lysozyme, and most particularly the tetragonal form of the protein, has become the default standard protein for use in macromolecule crystal nucleation and growth studies. There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests this proteins crystal nucleation and growth is by addition of associated species that are preformed by standard reversible concentration-driven self association processes in the bulk solution. The evidence includes high resolution AFM studies of the surface packing and of growth unit size at incorporation, fluorescence resonance energy transfer measurements of intermolecular distances in dilute solution, dialysis kinetics, and modeling of the growth rate data. We have developed a selfassociation model for the proteins crystal nucleation and growth. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. Further, it indicates that nucleation and crystal growth are not distinct mechanistically, but identical, with the primary difference being the probability that the particle will continue to grow or dissolve. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small (M.W. = 14,400), a structured octamer in the 4(sub 3) helix configuration (the proposed average sized growth unit) would have a M.W. = 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm commonly attach to the crystal. These measurements were made at approximately saturation conditions, and they reflect the sizes of species that both added or desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. Flow-driven effects resulting in misorientation upon addition to and incorporation into the crystal need only be a small fraction of a percentage to significantly affect the resulting crystal. One Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species will be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.

Pusey, Marc L.↗
NASA NTRS1993

Limited subsolidus diffusion in type B1 CAI: Evidence from Ti distribution in spinel

Most models of calcium aluminum-rich inclusions (CAI) have focused on early stages of formation by equilibrium crystallization of a homogeneous liquid. Less is known about the subsolidus cooling history of CAI. Chemical and isotopic heterogeneties on a scale of tens to hundreds of micrometers (e.g. MacPherson et al. (1989) and Podosek, et al. (1991)) suggest fairly rapid cooling with a minimum of subsolidus diffusion. However, transmission electron microscopy indicates that solid state diffusion may have been an important process at a smaller scale (Barber et al. 1984). If so, chemical evidence for diffusion could provide constraints on cooling times and temperatures. With this in mind, we have begun an investigation of the Ti distribution in spinels from two type B1 CAI from Allende to determine if post-crystallization diffusion was a significant process. The type B1 CAIs, 3529Z and 5241 have been described by Podosek et al. (1991) and by El Goresy et al. (1985) and MacPherson et al. (1989). We have analyzed spinels in these inclusions using the electron microprobe. These spinels are generally euhedral, range in size from less than 10 to 15 micron and are poikilitically enclosed by millimeter-sized pyroxene, melilite, and anorthite. Analyses were obtained from both the mantles and cores of the inclusions. Compositions of pyroxene in the vicinity of individual spinel grains were obtained by analyzing at least two points on opposite sides of the spinel and averaging the compositions. The pyroxene analyses were obtained within 15 microns of the spinel-pyroxene interface. No compositional gradients were observed within single spinel crystals. Ti concentrations in spinels included within pyroxene, melilite, and anorthite are presented.

Meeker, G. P.↗
NASA NTRS1999

Control of Meridional Flow in Circular Cylinders by a Travelling Axial Magnetic Field

Convective flow in a Bridgman or float zone configuration significantly affects the interface shape and segregation phenomena. While the primary causative factor for this flow is buoyancy induced convection in an enclosed Bridgman melt, the presence of a free surface gives rise to surface tension driven flows in the floating zone processing of melts. It is of interest to curtail these flows in order to realize near quiescent growth conditions that have shown to result in crystals with good longitudinal and radial homogeneity and thereby of better overall quality. While buoyancy effects can be reduced by careful processing in a low gravity (space) environment, the reduction of Marangoni flows due to surface tension variations is not that straight forward. Attempts have been made with some limited success with the use of external fields to affect the melt thermo-fluid behavior. The use of a static magnetic field that reduces convective contamination through the effects of a non-intrusively induced, dissipative Lorentz force in an electrically conducting melt is one such approach. Experiments have shown that axial fields of the order of 5 Tesla can significantly eliminate convection and yield close to diffusion limited crystal growth conditions. The generation and use of such high magnetic fields require substantial hardware and incur significant costs for its operation. Lately, the use of rotating magnetic fields has been tested in semiconductor crystal growth. The method is fairly well known and commonly used in metal processing but its adaptation to crystal growth of semiconductors is fairly recent. The elegance of the technique rests in its low power requirement (typically 10-20 milli-Tesla at 50-400 Hz) and its efficacy in curtailing deleterious temperature fluctuations in the melt. A rotating magnetic field imposes a rotational force and thereby induces a circulation within the melt that tends to dominate other sporadic convective effects. Thus a known low level of convective flow is introduced into the system. A new novel variation of the Lorentz force mechanism is proposed and investigated in this study. Since one of the desired process conditions in melt crystal growth is the minimization of convective effects, this investigation examines the use of an external field of magnetic origin to counteract existing convective flow within the melt. This is accomplished by utilizing a running or traveling axial magnetic wave in the system. The concept is similar to the use of vibrational means in order to induce streaming flows that oppose buoyant or surface tension driven convection in the system. The rotation direction as well as the magnitude (strength) of this circulation can be easily controlled by external inputs thus affording a direct means of controlling the developing shape of the crystallizing front (interface). The theoretical model of this technique is fully developed and presented in this paper. Results from the solution of the developed governing equations and boundary conditions are also presented. An experimental demonstration of the concept is presented through the suppression of natural convective flow in a mercury column. Implications to crystal growth systems will be fully explored in the final manuscript.

Mazuruk, K.↗
NASA NTRS2002

Development of a New Membrane Casting Apparatus for Studying Macrovoid Defects in Low-G

A new membrane-casting apparatus is developed for studying macrovoid defects in polymeric membranes made by the wet- and dry-casting process in low-gravity. Macrovoids are large (10-50 micron), open cavities interspersed among the smaller pores in the substructure under the gelled skin surface layer of the cast membrane. Although their occurrence is considered endemic to the wet- and dry-casting process since they can lead to compaction or skin rupture in the membrane process, recent studies suggest several useful applications such as transdermal and osmotic drug delivery systems, miniature bioreactors, etc. However, lack of knowledge about the macrovoid formation mechanism is an obstacle to further development of applications using them. An on-going debate is the role of the surface-tension-driven solutocapillary convection during macrovoid formation. The rapid growth of macrovoids within 1-5 seconds and the high polymer concentration in and near macrovoids make it difficult to explain the mechanism of macrovoid growth by diffusion alone, which is the widely accepted hypothesis proposed by Reuvers et al. The hypothesis advanced by our research group can explain this rapid growth via a mechanism that involves diffusion from the casting solution in the meta-stable region to the macrovoid enhanced by solutocapillary convection induced by the steep nonsolvent concentration gradient in the vicinity of the macrovoid. Since macrovoid growth is hypothesized to be the interplay of a solutocapillary-induced driving force counteracted by viscous drag and buoyancy, eliminate the latter provides a means for testing this hypothesis. Moreover, free convection mass transfer in the nonsolvent immersion bath used to cause phase-separation in membrane casting complicates developing a model for both the wet-casting process and macrovoid growth. The low-g environment minimizes gravitationally induced free convection thereby permitting a tractable solution to the ternary diffusion equations that characterize membrane formation. NASA's Parabolic Flight Research Aircraft provides a small window of low-g (approximately 25 s) that can be used to study macrovoid development in both wet- and dry-cast membranes if an appropriate casting apparatus is used. This casting apparatus should be able to cast the membrane in both low- and high-g in a manner so that essential one-dimensional mass transfer conditions are achieved to insure lateral uniformity in the membrane. The apparatus used in previous research on membrane casting in low-gravity was operated with the plunger driven mechanism. The spring-loaded plunger pushes the bottom block containing the polymer casting solution well directly under the absorbent chamber located in the upper stationary block. However, membranes made via this casting apparatus often displayed lateral nonuniformities that precluded obtaining quantitative information on the macrovoid growth process. Thus, it was necessary to determine the reason for these structural irregularities observed in the low-g casting apparatus. Both experimental as well as computer simulation studies of the low-g casting apparatus established that the impulsive action of the plunger caused the undesired structural nonuniformities. The simulation results showed that the width-to-depth aspect ratio of the shallow well that contains the casting solution in this apparatus was not an important factor in minimizing this problem. Even for a 40:1 (width : depth) aspect ratio, any convection induced by the horizontal motion of the interface of the casting solution will be damped out within 6.25x10(exp 4) seconds. However, the experimental studies revealed that the impulsive motion of the plunger caused a 'sloshing' of the casting solution that had to be eliminated. Therefore, the plungerdriven mechanism was changed to a cam-driven mechanism that did not cause any impulsive motion of the casting solution. Other refinements to this new membrane-casting apparatus include provision for removing the membranes from the casting wells in a less destructive manner. This was accomplished by using a slit geometry for the casting well that permitted disassembly for removal of the cast membrane. The materials used in the construction of this casting apparatus were chosen to insure wetting at the side walls and to maintain precise control of the thickness of the polymer solution in the casting well. An additional provision in this new casting apparatus is the ability to carry out both wet- as well as dry-casting. As such, this apparatus permitted the first studies of the wet-casting of polymeric membranes in low-g. Both wet- and dry-casting experiments on NASA's KC-135 research aircraft employing this new membrane-casting apparatus are scheduled in July 2002. The morphology of the resulting membranes will be characterized using an environmental scanning electron microscope (ESEM). The results of these low-g studies will be reported later.

Lee, Hanyong↗
NASA NTRS1998

Profile Evolution Simulation in Etching Systems Using Level Set Methods

Semiconductor device profiles are determined by the characteristics of both etching and deposition processes. In particular, a highly anisotropic etch is required to achieve vertical sidewalls. However, etching is comprised of both anisotropic and isotropic components, due to ion and neutral fluxes, respectively. In Ar/Cl2 plasmas, for example, neutral chlorine reacts with the Si surfaces to form silicon chlorides. These compounds are then removed by the impinging ion fluxes. Hence the directionality of the ions (and thus the ion angular distribution function, or IAD), as well as the relative fluxes of neutrals and ions determines the amount of undercutting. One method of modeling device profile evolution is to simulate the moving solid-gas interface between the semiconductor and the plasma as a string of nodes. The velocity of each node is calculated and then the nodes are advanced accordingly. Although this technique appears to be relatively straightforward, extensive looping schemes are required at the profile corners. An alternate method is to use level set theory, which involves embedding the location of the interface in a field variable. The normal speed is calculated at each mesh point, and the field variable is updated. The profile comers are more accurately modeled as the need for looping algorithms is eliminated. The model we have developed is a 2-D Level Set Profile Evolution Simulation (LSPES). The LSPES calculates etch rates of a substrate in low pressure plasmas due to the incident ion and neutral fluxes. For a Si substrate in an Ar/C12 gas mixture, for example, the predictions of the LSPES are identical to those from a string evolution model for high neutral fluxes and two different ion angular distributions.(2) In the figure shown, the relative neutral to ion flux in the bulk plasma is 100 to 1. For a moderately isotropic ion angular distribution function as shown in the cases in the left hand column, both the LSPES (top row) and rude's string method (bottom row) predict tapered profiles. The LSPES uses an AD with a FWHM = 13.5 degrees, and rude's model uses a ratio of sheath voltage to ion temperature of 50. The more anisotropic IADs produce profiles with more vertical sidewalls and a wider bottom surface, as shown in the right hand column. Here, the LSPES has an AD with a FWHM = 2 degrees, and rude's model uses a ratio of 500 of sheath voltage to ion temperature. The agreement between the LSPES and rude's model is excellent in both cases. We will present etching profiles generated by the LSPES, including calculations of the re-emitted fluxes of both neutrals and ions off of the profile walls. In addition, we will show the effect of geometric structures (overhangs, etc.) on the etching profiles. Other physical aspects, such as surface diffusion, will also be included in the model.

Hwang, Helen H.↗