Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface diffusion”

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.

591 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Creep and fracture of dispersion-strengthened materials

The creep and fracture of dispersion strengthened materials is reviewed. A compilation of creep data on several alloys showed that the reported values of the stress exponent for creep varied between 3.5 and 100. The activation energy for creep exceeded that for lattice self diffusion in the matrix in the case of some materials and a threshold stress behavior was generally reported in these instances. The threshold stress is shown to be dependent on the interparticle spacing and it is significantly affected by the initial microstructure. The effect of particle size and the nature of the dispersoid on the threshold stress is not well understood at the present time. In general, most studies indicate that the microstructure after creep is similar to that before testing and very few dislocations are usually observed. It is shown that the stress acting on a dispersoid due to a rapidly moving dislocation can exceed the particle yield strength of the G sub p/1000, where G sub p is the shear modulus of the dispersoid. The case when the particle deforms is examined and it is suggested that the dislocation creep threshold stress of the alloy is equal to the yield strength of the dispersoid under these conditions. These results indicate that the possibility that the dispersoid creep threshold stress is determined by either the particle yield strength or the stress required to detach a dislocation from the dispersoid matrix interface. The conditions under which the threshold stress is influenced by one or the other mechanism are discussed and it is shown that the particle yield strength is important until the extent of dislocation core relaxation at the dispersoid matrix interface exceeds about 25 pct. depending on the nature of the particle matrix combination. Finally, the effect of grain boundaries and grain morphology on the creep and fracture behavior of dispersoid strengthened alloys is examined.

Raj, Sai V.↗
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 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 NTRS1989

Raman intensity as a probe of concentration near a crystal growing in solution

The feasibility of using Raman spectral scattering signals for measurements of concentration profiles near a crystal interface during growth or dissolution is discussed. With KH2PO4 (KDP) as a test material, optical multichannel analyzer (OMA) detection of a solute Raman vibrational band provided direct quantification of solute concentration with band intensity. The intersection of incident laser and Raman collection optics provided 3-D selective point measurements of the solution concentration field. Unlike many other techniques, the Raman band intensity is not sensitive to the typical temperature variations. Precision calibration of Raman intensity versus KDP concentration with less than 1 pct standard deviation error levels was demonstrated. A fiber optic, which sampled incident laser intensity and coupled it to the OMA, provided a fully synchronized monitor of fluctuations in laser power to correlate with observed Raman signals. With 1 W of laser power at the sample, good data statistics required eight repeated data collections at approximately 2.5 min collection. The accumulated time represents the concentration measurement time at one spatial location. Photomicroscopy documented a 30 micrometer diameter by 200 micrometer of laser Raman scattering region in the solution near the crystal surface. The laser beam was able to approach up to 25 micrometer from the crystal surface. However, a crystal surface reflected intensity contribution was weakly detectable. Nucleated microcrystals were seen in the crystal-growing solution. These microcrystals convect right up to the crystal surface and indicate no quiet diffusion region under normal gravity conditions. Translation of the solution cell with respect to the optics caused systematic intensity errors.

Wilkinson, R. Allen↗
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 NTRS2002

Smyer H-Chondrite Impact-Melt Breccia and Evidence for Sulfur Vaporization

Smyer is an H-chondrite impact-melt breccia containing approx.20 vol% 0.5- to 13-mm-thick silicate-rich melt veins surrounding unmelted subrounded chondritic clasts up to 7 cm in maximum dimension. At the interface between some of the melt veins and chondritic clasts, there are troilite-rich regions consisting of unmelted. crushed 0.2- to 140-micron-size angular silicate grains and chondrule fragments surrounded by troilite and transected by thin troilite veins. Troilite fills every available fracture in the silicates. including some as thin as 0.1 microns. Little metallic Fe-Ni is present in these regions: the FeS/Fe modal ratio ranges from -25: 1 to approx.500: 1, far higher than the eutectic weight ratio of 7.5: 1. The texture of these regions indicates that the sulfide formed from a fluid of very low viscosity. The moderately high viscosity (0.2 poise) and large surface tension of liquid FeS, its inability to wet silicate grain surfaces at low oxygen fugacities. and the supereutectic FeS/Fe ratios in the troilite-rich regions indicate that the fluid was a vapor. It seems likely that during the shock event that melted Smyer, many silicates adjacent to the melt veins were crushed. Upon release of shock pressure. some of the troilite evaporated and dissociated. Molecules of S2 were transported and condensed into fractures and around tiny silicate grains: there, they combined with Fe from small adjacent metallic Fe-Ni grains to form troilite. The Ni content at the edges of some of these metal grains increased significantly; Co from these Ni-rich grains diffused into nearby kamacite. Impact-induced S volatilization may have played a major role in depleting the surface of 433 Eros (and other chondritic asteroids) in S.

Rubin, Alan E.↗
NASA NTRS2001

Growth and Morphology of Supercritical Fluids, a Fluid Physics Experiment Conducted on Mir, Complete

The Growth and Morphology of Supercritical Fluids (GMSF) is an international experiment facilitated by the NASA Glenn Research Center and under the guidance of U.S. principal investor Professor Hegseth of the University of New Orleans and three French coinvestigators: Daniel Beysens, Yves Garrabos, and Carole Chabot. The GMSF experiments were concluded in early 1999 on the Russian space station Mir. The experiments spanned the three science themes of near-critical phase separation rates, interface dynamics in near-critical boiling, and measurement of the spectrum of density fluctuation length scales very close to the critical point. The fluids used were pure CO2 or SF6. Three of the five thermostats used could adjust the sample volume with the scheduled crew time. Such a volume adjustment enabled variable sample densities around the critical density as well as pressure steps (as distinct from the usual temperature steps) applied to the sample. The French-built ALICE II facility was used for these experiments. It allows tightly thermostated (left photograph) samples (right photograph) to be controlled and viewed/measured. Its diagnostics include interferometry, shadowgraph, high-speed pressure measurements, and microscopy. Data were logged on DAT tapes, and PCMCIA cards and were returned to Earth only after the mission was over. The ground-breaking near critical boiling experiment has yielded the most results with a paper published in Physical Review Letters (ref. 1). The boiling work also received press in Science Magazine (ref. 2). This work showed that, in very compressible near-critical two-phase pure fluids, a vapor bubble was induced to temporarily overheat during a rapid heating of the sample wall. The temperature rise in the vapor was 23-percent higher than the rise in the driving container wall. The effect is due to adiabatic compression of the vapor bubble by the rapid expansion of fluid near the boundary during heatup. Thermal diffusivity is low near the critical point, so getting heat out of the compressed bubble is observably slow. This gives the appearance of a backward heat flow, or heat flow from a cold surface to a warm fluid.

Wilkinson, R. Allen↗
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 NTRS1993

Immiscible phase incorporation during directional solidification of hypermonotectics

Solidification processes in immiscible samples were investigated by directly observing the events taking place at the solid-liquid interface during directional solidification. Visualization of these events was made possible through the use of a transparent metal analog system and a temperature gradient stage assembly fitted to an optical microscope. The immiscible transparent analog system utilized was the succinonitrile-glycerol system. This system has been shown to exhibit the same morphological transitions as observed in metallic alloys of monotectic composition. Both monotectic and hypermonotectic composition samples were directionally solidified in order to gain an improved understanding of the manner in which the excess hypermonotectic liquid is incorporated into the solidifying structure. The processing conditions utilized prevented sedimentation of the excess hypermonotectic liquid by directionally solidifying the samples in very thin (13 microns), horizontally oriented cells. High thermal gradient to growth rate ratios (G/R) were used in an effort to prevent constitutional supercooling and the subsequent formation of L(sub 2) droplets in advance of the solidification front during the growth of fibrous composite structures. Results demonstrated that hypermonotectic composites could be produced in samples up to two weight percent off of the monotectic composition by using a G/R ratio greater than or equal to 4.6 x 10(exp 4) C(s)/mm(sup 2) to avoid constitutional supercooling. For hypermonotectic samples processed with G/R ratios below 4.6 x 10(exp 4) C(s)/mm(sup 2), constitutional supercooling occurred and resulted in slight interfacial instability. For these samples, two methods of incorporation of the hypermonotectic liquid were observed and are reported. The correlation between the phase spacing, lambda, and the growth rate, R, was examined and was found to obey a relationship generally associated with a diffusion controlled coupled growth process. For samples with compositions ranging from the monotectic composition up to 2 percent off of the monotectic composition, data indicated that the square of the phase spacing (lambda) varied linearly with the inverse of the growth rate (R).

Andrews, J. Barry↗
NASA NTRS1999

Microgravity Transport Phenomena Experiment (MTPE) Overview

The Microgravity Transport Phenomena Experiment (MTPE) is a fluids experiment supported by the Fundamentals in Biotechnology program in association with the Human Exploration and Development of Space (BEDS) initiative. The MTP Experiment will investigate fluid transport phenomena both in ground based experiments and in the microgravity environment. Many fluid transport processes are affected by gravity. Osmotic flux kinetics in planar membrane systems have been shown to be influenced by gravimetric orientation, either through convective mixing caused by unstably stratified fluid layers, or through a stable fluid boundary layer structure that forms in association with the membrane. Coupled transport phenomena also show gravity related effects. Coefficients associated with coupled transport processes are defined in terms of a steady state condition. Buoyancy (gravity) driven convection interferes with the attainment of steady state, and the measurement of coupled processes. The MTP Experiment measures the kinetics of molecular migration that occurs in fluids, in response to the application of various driving potentials. Three separate driving potentials may be applied to the MTP Experiment fluids, either singly or in combination. The driving potentials include chemical potential, thermal potential, and electrical potential. Two separate fluid arrangements are used to study membrane mediated and bulk fluid transport phenomena. Transport processes of interest in membrane mediated systems include diffusion, osmosis, and streaming potential. Bulk fluid processes of interest include coupled phenomena such as the Soret Effect, Dufour Effect, Donnan Effect, and thermal diffusion potential. MTP Experiments are performed in the Microgravity Transport Apparatus (MTA), an instrument that has been developed specifically for precision measurement of transport processes. Experiment fluids are contained within the MTA fluid cells, designed to create a one dimensional flow geometry of constant cross sectional area, and to facilitate fluid filling and draining operations in microgravity. The fluid cells may be used singly for bulk solutions, or in a Stokes diaphragm configuration to investigate membrane mediated phenomena. Thermal and electrical driving potentials are applied to the experiment fluids through boundary plates located at the ends of the fluid cells. In the ground based instrument, two constant temperature baths circulate through reservoirs adjacent to the boundary plates, and establish the thermal environment within the fluid cells. The boundary plates also serve as electrodes for measurement and application of electrical potentials. The Fluid Manipulation System associated with the MTA is a computer controlled system that enables storage and transfer of experiment fluids during on orbit operations. The system is used to automatically initiate experiments and manipulate fluids by orchestrating pump and valve operations through scripted sequences. Unique technologies are incorporated in the MTA for measurement of fluid properties. Volumetric Flow Sensors have been developed for precision measurement of total fluid volume contained within the fluid cells over time. This data is most useful for measuring the kinetics of osmosis, where fluid is transported from one fluid cell to another through a semipermeable membrane. The MicroSensor Array has been designed to perform in situ measurement of several important fluid parameters, providing simultaneous measurement of solution composition at multiple locations within the experiment fluids. Micromachined sensors and interface electronics have been developed to measure temperature, electrical conductivity, pH, cation activity, and anion activity. The Profile Refractometer uses a laser optical system to directly image the fluid Index of Refraction profile that exists along the MTA fluid cell axis. A video system acquires images of the RI profile over time, and records the transport kinetics that occur upon application of chemical, thermal, or electrical driving potentials. Image processing algorithms have been developed to analyze the refractometer images on a pixel by pixel basis, calibrating and scaling the measured Index of Refraction profile to correlated solution properties of interest such as density, concentration, and temperature. Additional software has been developed to compile the processed images into a three dimensional matrix that contains fluid composition data as a function of experiment time and position in the fluid cell. These data are combined with data from the other sensor systems, and analyzed in the context of transport coefficients associated with the various transport phenomena. Analysis protocols have been developed to measure the transient kinetics, and steady state distribution of fluid components that occur in response to the applied driving potentials. The results are expressed in terms of effective transport coefficients. Experiments have been performed using a variety of solutes, and results generated are that are in agreement with published transport coefficient values.

Mason, Larry W.↗
NASA NTRS1995

First spectral observations of the diffuse background with the Extreme Ultraviolet Explorer

We present the first results from the analysis of the spectroscopic observations of diffuse extreme ultraviolet (EUV) emission taken with the Extreme Ultraviolet Explorer (EUVE) spectrometers in the wavelength range 160-740 A. Although not designed or optimized for diffuse observation, the EUVE spectrometers are the most sensitive diffuse EUV spectrometer in orbit. The spectral resolution for diffuse emission of the medium and long-wavelength spectrometers are 17 and 34 A FWHM, respectively. During the period from 1992 July 25 to 1992 August 19, the spectrometers surveyed a 2 x 20 deg field scanned from (l, b) = (24 deg, -28 deg) to (44 deg, -74 deg) with a total effective exposure time of 575,232 s. The only emission lines detected were those of He I and He II (584, 537, 304 A) with intensities consistent with local geocoronal and/or interplanetary scattering of solar radiation (584 A = 1.30 rayleighs; 537 A = 0.040 R; and 304 A = 0.029 R). Models of the soft X-ray background, which results from a 10(exp 6) K plasma (Local Bubble) surrounding the neutral gas near the Sun (Local Cloud), predict that most of the flux from the hot plasma appears as emission lines in the EUV. We have compared these spectral predictions with our observations to place limits on the emission measure versus temperature of the proposed hot plasma. Using the same plasma model, we derived emissions measures for our data and the C and B soft X-ray bands of the Wisconsin rocket survey. We find that our limits for the plasma emission measure are a factor of 5-10 below the C- and B-band emission measures over the temperature range from 10(exp 5.7) to 10(exp 6.4) K. We explore possible scenarios that could reconcile our results with the X-ray surveys and conclude that depletion or a nonequilibrium plasma state rather than absorption are the more likely explanations of the discrepancy. We also show that our spectrum is inconsistent with the spectrum from the approximately 10(exp 5) K gas at the conductive interface between the hot Local Bubble and the cooler Local Cloud given by Slavin (1989). In addition, we place new limits on the helium ionization parameter in the Sun's vicinity caused by the 10(exp 6) K plasma in the Local Bubble.

Jelinksy, P.↗
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.↗
NASA NTRS2000

Over-Expression, Purification and Crystallization of Human Dihydrolipoamide Dehydrogenase

Dehydrolipoamide dehydrogenase (E3; dihydrolipoan-tide:NAD+ oxidoreductase, EC 1.8.1.4) is a common catalytic component found in pyruvate dehydrogenase complex, alpha-ketoglutarate dehydrogenase complex, and branched-chain cc-keto acid dehydrogenase complex. E3 is also a component (referred to as L protein) of the glycine cleavage system in bacterial metabolism (2). Active E3 forms a homodimer with four distinctive subdomain structures (FAD binding, NAD+ binding, central and interface domains) with non-covalently but tightly bound FAD in the holoenzyme. Deduced amino acids from cloned full-length human E3 gene showed a total of 509 amino acids with a leader sequence (N-terminal 35 amino acids) that is excised (mature form) during transportation of expressed E3 into mitochondria membrane. So far, three-dimensional structure of human E3 has not been reported. Our effort to achieve the elucidation of the X-ray crystal structure of human E3 will be presented. Recombinant pPROEX-1 expression vector (from GIBCO BRL Life Technologies) having the human E3 gene without leader sequence was constructed by Polymerase Chain Reaction (PCR) and subsequent ligation, and cloned in E.coli XL1-Blue by transformation. Since pPROEX-1 vector has an internal His-tag (six histidine peptide) located at the upstream region of a multicloning site, one-step affinity purification of E3 using nickelnitriloacetic acid (Ni-NTA) agarose resin, which has a strong affinity to His-tag, was feasible. Also a seven-amino-acid spacer peptide and a recombinant tobacco etch virus protease recognition site (seven amino acids peptide) found between His-tag and first amino acid of expressed E3 facilitated the cleavage of His-tag from E3 after the affinity purification. By IPTG induction, ca. 15 mg of human E3 (mature form) was obtained from 1L LB culture with overnight incubation at 25C. Over 98% of purity of E3 from one-step Ni-NTA agarose affinity purification was confirmed by SDS-PAGE analysis. For crystallization, E3 samples were prepared with and without His-tag. To minimize the aggregation of E3, apo- and holo- forms of E3s were tested, as well as a mutated E3. Dynamic light scattering measurements revealed that the E3 preparations without His-tag and substrate are highly monodispersive with regard to homodimers. Consequent crystallization trials of this E3 preparation led to single crystals of E3 grown by the vapor diffusion method. Crystals were obtained within a few days from solution containing poly (ethylene glycol) monomethyl ether 5000 as a precipitant. Autoindexing and integration of the X-ray diffraction data showed that E3 crystals belong to an orthorhombic system with unit cell parameters a-- 123. 1, b= 165.3 and c=214.3A. Further optimization of protein preparation and crystallization experiments for the structural determination will be discussed.

Hong, Y. S.↗
NASA NTRS1993

CFD Parametric Study of Consortium Impeller

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design impact in a productive manner. The main goal of this study is to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A Finite Difference Navier-Stokes flow solver, FDNS, which includes the extended k-epsilon turbulence model and appropriate moving interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. A second-order central difference scheme plus adaptive dissipation terms was employed in the FDNS code, along with a predictor plus multi-corrector pressure-based solution procedure. The multi-zone, multi-block capability allows the FDNS code to efficiently solve flow fields with complicated geometry. The FDNS code has been benchmarked by analyzing the pump consortium inducer, and it provided satisfactory results. In the present study, a CFD parametric study of the pump consortium impeller was conducted using the FDNS code. The pump consortium impeller, with partial blades, is a new design concept of the advanced rocket engines. The parametric study was to analyze the baseline design of the consortium impeller and its modification which utilizes TANDEM blades. In the present study, the TANDEM blade configuration of the consortium impeller considers cut full blades for about one quarter chord length from the leading edge and clocks the leading edge portion with an angle of 7.5 or 22.5 degrees. The purpose of the present study is to investigate the effect and trend of the TANDEM blade modification and provide the result as a design guideline. A 3-D flow analysis, with a 103 x 23 x 30 mesh grid system and with the inlet flow conditions measured by Rocketdyne, was performed for the baseline consortium impeller. The numerical result shows that the mass flow rate splits through various blade passages are relatively uniform. Due to the complexity of blade geometries, the TANDEM blade configurations were analyzed with the multi-zone grid structure. Both the 7.5 deg- and the 22.5 deg-clocking TANDEM blade cases utilized a 80K mesh system. The numerical result of two TANDEM blade modifications indicates the efficiency and the head are worse than those of the baseline case due to larger flow distortion. The gap between the TANDEM blade and the full blade allows the flow passes through and heavily loads the pressure side of the partial blade such that flow reversal occurs near the suction side of the splitter. The flow split at the exit of impeller blades is very non-uniform for TANDEM blade cases, and this will greatly induce the side load on the diffuser. Therefore, the TANDEM blade modification in the present CFD analysis does not improve the performance of the consortium impeller.

Cheng, Gary C.↗
NASA NTRS2000

A Collaborative Analysis Tool for Integrated Hypersonic Aerodynamics, Thermal Protection Systems, and RBCC Engine Performance for Single Stage to Orbit Vehicles

Presented is a computer-based tool that connects several disciplines that are needed in the complex and integrated design of high performance reusable single stage to orbit (SSTO) vehicles. Every system is linked to every other system, as is the case of SSTO vehicles with air breathing propulsion, which is currently being studied by NASA. An RBCC propulsion system integrates airbreathing and rocket propulsion into a single engine assembly enclosed within a cowl or duct. A typical RBCC propulsion system operates as a ducted rocket up to approximately Mach 3. Then there is a transition to a ramjet mode for supersonic-to-hypersonic acceleration. Around Mach 8 the engine transitions to a scramjet mode. During the ramjet and scramjet modes, the integral rockets operate as fuel injectors. Around Mach 10-12 (the actual value depends on vehicle and mission requirements), the inlet is physically closed and the engine transitions to an integral rocket mode for orbit insertion. A common feature of RBCC propelled vehicles is the high degree of integration between the propulsion system and airframe. At high speeds the vehicle forebody is fundamentally part of the engine inlet, providing a compression surface for air flowing into the engine. The compressed air is mixed with fuel and burned. The combusted mixture must be expanded to an area larger than the incoming stream to provide thrust. Since a conventional nozzle would be too large, the entire lower after body of the vehicle is used as an expansion surface. Because of the high external temperatures seen during atmospheric flight, the design of an airbreathing SSTO vehicle requires delicate tradeoffs between engine design, vehicle shape, and thermal protection system (TPS) sizing in order to produce an optimum system in terms of weight (and cost) and maximum performance. To adequately determine the performance of the engine/vehicle, the Hypersonic Flight Inlet Model (HYFIM) module was designed to interface with the RBCC engine model. HYFIM performs the aerodynamic analysis of forebodies and inlet characteristics of RBCC powered SSTO launch vehicles. HYFIM is applicable to the analysis of the ramjet/scramjet engine operations modes (Mach 3-12), and provides estimates of parameters such as air capture area, shock-on-lip Mach number, design Mach number, compression ratio, etc., based on a basic geometry routine for modeling axisymmetric cones, 2-D wedge geometries. HYFIM also estimates the variation of shock layer properties normal to the forebody surface. The thermal protection system (TPS) is directly linked to determination of the vehicle moldline and the shaping of the trajectory. Thermal protection systems to maintain the structural integrity of the vehicle must be able to mitigate the heat transfer to the structure and be lightweight. Herein lies the interdependency, in that as the vehicle's speed increases, the TPS requirements are increased. And as TPS masses increase the effect on the propulsion system and all other systems is compounded. The need to analyze vehicle forebody and engine inlet is critical to be able to design the RBCC vehicle. To adequately determine insulation masses for an RBCC vehicle, the hypersonic aerodynamic environment and aeroheating loads must be calculated and the TPS thicknesses must be calculated for the entire vehicle. To accomplish this an ascent or reentry trajectory is obtained using the computer code Program to Optimize Simulated Trajectories (POST). The trajectory is then used to calculate the convective heat rates on several locations on the vehicles using the Miniature Version of the JA70 Aerodynamic Heating Computer Program (MINIVER). Once the heat rates are defined for each body point on the vehicle, then insulation thicknesses that are required to maintain the vehicle within structural limits are calculated using Systems Improved Numerical Differencing Analyzer (SINDA) models. If the TPS masses are too heavy for the performance of the vehicle the process may be repeated altering the trajectory or some other input to reduce the TPS mass. E-PSURBCC is an "engine performance" model and requires the specification of inlet air static temperature and pressure as well as Mach number (which it pulls from the HYFIM and POST trajectory files), and calculates the corresponding stagnation properties. The engine air flow path geometry includes inlet, a constant area section where the rocket is positioned, a subsonic diffuser, a constant area afterburner, and either a converging nozzle or a converging-diverging nozzle. The current capabilities of E-PSURBCC ejector and ramjet mode treatment indicated that various complex flow phenomena including multiple choking and internal shocks can occur for combinations of geometry/flow conditions. For a given input deck defining geometry/flow conditions, the program first goes through a series of checks to establish whether the input parameters are sound in terms of a solution path. If the vehicle/engine performance fails mission goals, the engineer is able to collaboratively alter the vehicle moldline to change aerodynamics, or trajectory, or some other input to achieve orbit. The problem described is an example of the need for collaborative design and analysis. RECIPE is a cross-platform application capable of hosting a number of engineers and designers across the Internet for distributed and collaborative engineering environments. Such integrated system design environments allow for collaborative team design analysis for performing individual or reduced team studies. To facilitate the larger number of potential runs that may need to be made, RECIPE connects the computer codes that calculate the trajectory data, aerodynamic data based on vehicle geometry, heat rate data, TPS masses, and vehicle and engine performance, so that the output from each tool is easily transferred to the model input files that need it.

Stanley, Thomas Troy↗