Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Moving interface”

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 55 records · Page 3

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics

Experiment 13: The Study of Dopant Segregation Behavior During the Growth of GaAs in Microgravity on USML-2

An investigation into the segregation behavior of selenium doped gallium arsenide (Se/GaAs) during directional solidification in the microgravity environment was conducted using the Crystal Growth Furnace (CGF) aboard the second United States Microgravity Laboratory (USML-2). Two crystals were successfully processed on USML-2, which lasted from October 20 to November 7, 1995. The first sample was processed for 67 hours, 45 minutes (MET 5/04:53:45-8/00:23:50) and included 19 hours of growth at 0.5 microns/sec which yielded 3.42 cm of sample length, and 5 hours of growth at 1.5 microns/sec which yielded 2.7 cm of sample. During the second experiment, the furnace temperature was adjusted to move the melt-solid interface position towards the hot end of the furnace. The second sample was processed for 50 hours, 10 minutes (MET 8/18:48:49-10/21:58:54) and included 11 hours of growth at 0.5 microns/sec which yielded 1.98 cm of sample, and 1 hour, 25 minutes of growth at 5.0 microns/sec which yielded 2.6 cm of sample. This sample provides an order of magnitude change in growth rate and reproduces one of the growth rates used during USML-1. In contrast to the results from USML-1, no voids were present in either crystal grown on USML-2. The absence of voids in either sample indicates that growth rate changes alone were not responsible for the formation of voids found in the crystals grown on USML-1. Sections of the ground-based and flight crystals grown on USML-2 were cut and polished. All of the interface demarcation lines expected from the current pulse interface demarcation (CPID) system have been identified. These measurements have been analyzed for interface positions, interface shapes, and growth rates. Using a newly developed technique, based on experimental and numerical results, the seeding interface reproducibility from run to run was <= 2.5 mm. The seeding interface position could be controllably moved, with respect to the furnace zones, by adjusting the control set points of the heating zones. The interface shapes flattened slightly as the interface position moved closer to the hot zone but was always an unfavorable concave into the solid shape. The growth rate was found to equal the furnace translation rate, after a 2 -hour transient, for growth rates <= 1.0 microns/sec. Segregation measurements for the ground-based crystals are indicative of complete mixing behavior, as expected. Segregation measurements of the flight crystals are still in progress.

Matthiesen, David H.

Qualitative investigation of cryogenic injected shock dissipation

The injection of fluid nitrogen into the flow discharged from a two-dimensional, ambient, static-temperature and static pressure Mach = 2.7 nitrogen flow tunnel was observed. A bow shock stands approximately one diameter off the end of the 3.2 mm o.d. x 1.6 mm i.d. injection tube placed directly into the tunnel exhaust. Cryogenic injection creates a high density region in the injected region and within the bow shock wake but the standoff distance remains unchanged. However as the temperature reaches a critical value the sharp shock interface begins to fade resembling the density gradient seen at the interface of a supercritical fluid and the shock interface begins to move into the supersonic stream. As the injection temperature decreases, the interface continues to move into the supersonic stream and becomes more diffuse; the shock interface can not be distinguished. The phenomena is completely reversible.

Hendricks, R. C.

Observation of an Aligned Gas - Solid "Eutectic" during Controlled Directional Solidification Aboard the International Space Station - Comparison with Ground-based Studies

Direct observation of the controlled melting and solidification of succinonitrile was conducted in the glovebox facility of the International Space Station (ISS). The experimental samples were prepared on ground by filling glass tubes, 1 cm ID and approximately 30 cm in length, with pure succinonitrile (SCN) in an atmosphere of nitrogen at 450 millibar pressure for eventual processing in the Pore Formation and Mobility Investigation (PFMI) apparatus in the glovebox facility (GBX) on board the ISS. Real time visualization during controlled directional melt back of the sample showed nitrogen bubbles emerging from the interface and moving through the liquid up the imposed temperature gradient. Over a period of time these bubbles disappear by dissolving into the melt. Translation is stopped after melting back of about 9 cm of the sample, with an equilibrium solid-liquid interface established. During controlled re-solidification, aligned tubes of gas were seen growing perpendicular to the planar solid/liquid interface, inferring that the nitrogen previously dissolved into the liquid SCN was now coming out at the solid/liquid interface and forming the little studied liquid = solid + gas eutectic-type reaction. The observed structure is evaluated in terms of spacing dimensions, interface undercooling, and mechanisms for spacing adjustments. Finally, the significance of processing in a microgravity environment is ascertained in view of ground-based results.

Grugel, R. N.

A Real Time Investigation of Morphological Evolution During Solidification of Different Alloy Systems

Solidification phenomenon which occur at the solid/liquid (s/I) interface play a major role in the determination of structure and hence the technologically important properties of a casting. However, metals being opaque, conclusions related to several important phenomenon such as boundary layer thickness, morphological evolution, and eutectic and cell spacing are deduced from quenching experiments and subsequent post solidification metallographic analysis. Consequently, limited information is obtained about the dynamics of the process. This paper will discuss the recent efforts at the Space Science Laboratory, NASA Marshall Space Flight Center, to view and quantify in-situ and in real time the dynamics of the solidification process and to measure interfacial undercooling. First, a high resolution x-ray transmission microscope (XTM) has been developed to monitor fundamental interfacial phenomena during directional solidification of metals and alloys. The XTM operates in the range of 10-100 KeV and through projection is capable of achieving magnification of up to 16OX. Secondly, an innovative collapsible furnace has been designed to quantify interfacial undercooling by measuring the temperature of a moving s/I interface in reference to a fixed s/l interface. This measurement technique is non-intrusive in nature and is based on the Seebeck principle. In this paper real time results obtained to characterize the dynamics of irregular eutectic spacing will be presented. As an example fiber to lamella or plate transition in the Al-Al2Au eutectic system will be discussed. Further, a resolution limit of 25 micron has permitted viewing in real time morphological instability and cellular growth in Al-Au and Al-Ag systems. Simultaneously, a systematic investigation has been carried out to measure interfacial undercooling for Pb-1 wt.% Sn at and near the marginal stability regime. In conjunction with the XTM observations this study attempts to validate existing relationships between undercooling and growth velocity during plane front growth, marginal stability regime, and stable cellular growth.

Sen, S.

Long Life Technology Work at Rockwell International Space Division

This paper presents highlights of long-life technology oriented work performed at the Space Division of Rockwell International Corporation under contract to NASA. This effort included evaluation of Saturn V launch vehicle mechanical and electromechanical components for potential extended life capabilities, endurance tests, and accelerated aging experiments. A major aspect was evaluation of the components at the subassembly level (i.e., at the interface between moving surfaces) through in-depth wear analyses and assessments. Although some of this work is still in progress, preliminary conclusions are drawn and presented, together with the rationale for each. The paper concludes with a summary of the effort still remaining.

D K Huzel

Finite element simulation of temperature dependent free surface flows

The method of Engelman and Sani (1984) for a finite-element simulation of incompressible surface flows with a free and/or moving fluid interface, such as encountered in crystal growth and coating and polymer technology, is extended to temperature-dependent flows, including the effect of temperature-dependent surface tension. The basic algorithm of Saito and Scriven (1981) and Ruschak (1980) has been generalized and implemented in a robust and versatile finite-element code that can be employed with relative ease for the simulation of free-surface problems in complex geometries. As a result, the costly dependence on the Newton-Raphson algorithm has been eliminated by replacing it with a quasi-Newton iterative method, which nearly retains the superior convergence properties of the Newton-Raphson method.

Engelman, M. S.

Numerical study of 3-D inducer and impeller for pump model development

Current design of high-performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study is to develop a robust and effective CFD pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, two key components of the turbopump, the inducer and impeller, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne.

Cheng, G. C.

Inducer analysis/pump model development

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study was 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 an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, three key components of the turbopump, the inducer, impeller, and diffuser, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne. For the numerical calculation of inducer flows with tip clearance, the turbulence model and grid spacing are very important. Meanwhile, the development of the cross-stream secondary flow, generated by curved blade passage and the flow through tip leakage, has a strong effect on the inducer flow. Hence, the prediction of the inducer performance critically depends on whether the numerical scheme of the pump model can simulate the secondary flow pattern accurately or not. The impeller and diffuser, however, are dominated by pressure-driven flows such that the effects of turbulence model and grid spacing (except near leading and trailing edges of blades) are less sensitive. The present CFD pump model has been proved to be an efficient and robust analytical tool for pump design due to its very compact numerical structure (requiring small memory), fast turnaround computing time, and versatility for different geometries.

Cheng, Gary C.

Analytical studies of hypersonic viscous dissociated flows

This project primarily dealt with integral boundary-layer solution techniques that are directly applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles like X-33 in the vicinity of stagnation points, windward centerlines, and swept-wing leading edges. The analyses include effects of finite-rate gas chemistry across the boundary layer and finite-rate catalysis of atom recombination at the surface. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (and expensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations was developed. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. The post-processed data base allows a designer at a workstation to ask and answer the following questions: (1) How much does the heating change if one uses a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change when one moves the interface of two different TPS materials with different catalytic efficiencies for the purpose of reducing vehicle weight and expense? The answer to the second question is particularly critical, because abrupt changes from low catalytic efficiency to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. A secondary issue that was addressed involves the prediction of heating levels in the vicinity of sharp corners that are transverse to or aligned with the flow. An example of the first case is heating at the edge of the COMET reentry module. An example of the second case is heating along the side edge of a deflected body flap on an SSV. The difficulty of putting grids in the vicinity of such corners with continuously varying metric coefficients causes problems in CFD predictions. A preliminary theory for prediction that says the heating at the corner is X percent of the heating N boundary-layer thicknesses inboard was developed. This will prove useful to analytically evaluate the possible benefits of rounding the edges of these configurations and defining how much rounding is sufficient.

Inger, George R.

Numerical Calculation of the Drag Force Acting on an Insoluble Particle Moving in Front of a Solidifying Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. Whether a particle will be engulfed or continuously pushed by SLI is essentially determined by the balance of forces acting on the particle. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that not only the steady-state approach is inappropriate to model the interaction process but also that even at steady-state the theoretical solution underestimates the drag force. A regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gave a relationship of the form F(sub d) (exp num) = (square root of 3) pi eta V(sub p)((R(sub p) / d) (exp 10 gamma/3)) as compared to the theoretical solution F(sub d) (exp theor) = 3 X (square root of 2) pi eta V(sub p)((R(sub p) / d) (exp 3/2)) where F(sub d) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma is approximately equal to 0.577), and d is the width of the gap between the particle and the SLI. Model validation as well as the influence of the interface shape on the value of F(sub d) will also be discussed.

Catalina, Adrian V.

Advanced Payload Transfer Measurement System (APT-MS) Mechanical Features

During the last two feet or so of transfer, for a large or heavy payload to its restraining fixture, the consequences of unplanned contact (or impact) between payload and support structure can range from merely annoying to something approaching disaster. Current transfer methods employ technicians with meter sticks stationed at the various hold-down locations to estimate the distances to contact. This information is communicated orally to the crane operator. It is understandable that this work proceeds carefully and therefore slowly. The objectives of this study are (1) to explore mechanical sources of measurement error; and (2) to develop an improved prototype design that is compact, inexpensive, and addresses the measurement error issues. APT is a relatively inexpensive electro-mechanical device that reduces both time and manpower required to make shuttle payload transfers. APT measurement system would provide a GUI for the 'move-conductor' (crane operator) so he could see the displacements of all hold-down interfaces as they move together and thus have a more accurate, comprehensive and 'real-time' picture of the engagement activity. An error model attempts to include all estimatable sources of mechanical error and design features were introduced to reduce or eliminate major sources of error.

Brewer, William V.

The interaction of plume heads with compositional discontinuities in the Earth's mantle

The effects of compositional discontinuities of density and viscosity in the Earth's mantle on the ascent of mantle plume heads is studied using a boundary integral numerical technique. Three specific problems are considered: (1) a plume head rising away from a deformable interface, (2) a plume head passing through an interface, and (3) a plume head approaching the surface of the Earth. For the case of a plume attached to a free-surface, the calculated time-dependent plume shapesare compared with experimental results. Two principle modes of plume head deformation are observed: plume head elingation or the formation of a cavity inside the plume head. The inferred structure of mantle plumes, namely, a large plume head with a long tail, is characteristic of plumes attached to their source region, and also of buoyant material moving away from an interface and of buoyant material moving through an interface from a high- to low-viscosity region. As a rising plume head approaches the upper mantle, most of the lower mantle will quickly drain from the gap between the plume head and the upper mantle if the plume head enters the upper mantle. If the plume head moves from a high- to low-viscosity region, the plume head becomes significantly elongated and, for the viscosity contrasts thought to exist in the Earth, could extend from the 670 km discontinuity to the surface. Plume heads that are extended owing to a viscosity decrease in the upper mantle have a cylindrical geometry. The dynamic surface topography induced by plume heads is bell-shaped when the top of the plume head is at depths greater than about 0.1 plume head radii. As the plume head approaches the surface and spreads, the dynamic topography becomes plateau-shaped. The largest stresses are produced in the early stages of plume spreading when the plume head is still nearly spherical, and the surface expression of these stresses is likely to be dominated by radial extension. As the plume spreads, compressional stresses on the surface are produced beyond the edges of the plume; consequently, extensional features will be produced above the plume head and may be surrounded by a ring of compressional features.

Manga, Michael

Mathematical model for the Bridgman-Stockbarger crystal growing system

In a major technical breakthrough, a computer model for Bridgman-Stockbarger crystal growth was developed. The model includes melt convection, solute effects, thermal conduction in the ampule, melt, and crystal, and the determination of the curved moving crystal-melt interface. The key to the numerical method is the use of a nonuniform computational mesh which moves with the interface, so that the interface is a mesh surface. In addition, implicit methods are used for advection and diffusion of heat, concentration, and vorticity, for interface movement, and for internal gracity waves. This allows large time-steps without loss of stability or accuracy. Numerical results are presented for the interface shape, temperature distribution, and concentration distribution, in steady-state crystl growth. Solutions are presented for two test cases using water, with two different salts in solution. The two diffusivities differ by a factor of ten, and the concentrations differ by a factor of twenty.

Roberts, G. O.

Investigation of compositional segregation during unidirectional solidification of solid solution semiconducting alloys

Compositional segregation of solid solution semiconducting alloys in the radial direction during unidirectional solidification was investigated by calculating the effect of a curved solid liquid interface on solute concentration at the interface on the solid. The formulation is similar to that given by Coriell, Boisvert, Rehm, and Sekerka except that a more realistic cylindrical coordinate system which is moving with the interface is used. Analytical results were obtained for very small and very large values of beta with beta = VR/D, where V is the velocity of solidification, R the radius of the specimen, and D the diffusivity of solute in the liquid. For both very small and very large beta, the solute concentration at the interface in the solid C(si) approaches C(o) (original solute concentration) i.e., the deviation is minimal. The maximum deviation of C(si) from C(o) occurs for some intermediate value of beta.

Wang, J. C.

Instability Mechanisms of Thermally-Driven Interfacial Flows in Liquid-Encapsulated Crystal Growth

During the past year, a great deal of effort was focused on the enhancement and refinement of the computational tools developed as part of our previous NASA grant. In particular, the interface mollification algorithm developed earlier was extended to incorporate the effects of surface-rheological properties in order to allow the study of thermocapillary flows in the presence of surface contamination. These tools will be used in the computational component of the proposed research in the remaining years of this grant. A detailed description of the progress made in this area is provided elsewhere. Briefly, the method developed allows for the convection and diffusion of bulk-insoluble surfactants on a moving and deforming interface. The novelty of the method is its grid independence: there is no need for front tracking, surface reconstruction, body-fitted grid generation, or metric evaluations; these are all very expensive computational tasks in three dimensions. For small local radii of curvature there is a need for local grid adaption so that the smearing thickness remains a small fraction of the radius of curvature. A special Neumann boundary condition was devised and applied so that the calculated surfactant concentration has no variations normal to the interface, and it is hence truly a surface-defined quantity. The discretized governing equations are solved subsequently using a time-split integration scheme which updates the concentration and the shape successively. Results demonstrate excellent agreement between the computed and exact solutions.

Haj-Hariri, Hossein

Heat Transfer Analysis in Wire Bundles for Aerospace Vehicles

Design of wiring for aerospace vehicles relies on an understanding of "ampacity" which refers to the current carrying capacity of wires, either, individually or in wire bundles. Designers rely on standards to derate allowable current flow to prevent exceedance of wire temperature limits due to resistive heat dissipation within the wires or wire bundles. These standards often add considerable margin and are based on empirical data. Commercial providers are taking an aggressive approach to wire sizing which challenges the conventional wisdom of the established standards. Thermal modelling of wire bundles may offer significant mass reduction in a system if the technique can be generalized to produce reliable temperature predictions for arbitrary bundle configurations. Thermal analysis has been applied to the problem of wire bundles wherein any or all of the wires within the bundle may carry current. Wire bundles present analytical challenges because the heat transfer path from conductors internal to the bundle is tortuous, relying on internal radiation and thermal interface conductance to move the heat from within the bundle to the external jacket where it can be carried away by convective and radiative heat transfer. The problem is further complicated by the dependence of wire electrical resistivity on temperature. Reduced heat transfer out of the bundle leads to higher conductor temperatures and, hence, increased resistive heat dissipation. Development of a generalized wire bundle thermal model is presented and compared with test data. The steady state heat balance for a single wire is derived and extended to the bundle configuration. The generalized model includes the effects of temperature varying resistance, internal radiation and thermal interface conductance, external radiation and temperature varying convective relief from the free surface. The sensitivity of the response to uncertainties in key model parameters is explored using Monte Carlo analysis.

S L Rickman

Avionics Architectures for Exploration: Ongoing Efforts in Human Spaceflight

The field of Avionics is advancing far more rapidly in terrestrial applications than in spaceflight applications. Spaceflight Avionics are not keeping pace with expectations set by terrestrial experience, nor are they keeping pace with the need for increasingly complex automation and crew interfaces as we move beyond Low Earth Orbit. NASA must take advantage of the strides being made by both space-related and terrestrial industries to drive our development and sustaining costs down. This paper describes ongoing efforts by the Avionics Architectures for Exploration (AAE) project chartered by NASA's Advanced Exploration Systems (AES) Program to evaluate new avionic architectures and technologies, provide objective comparisons of them, and mature selected technologies for flight and for use by other AES projects. The AAE project team includes members from most NASA centers, and from industry. It is our intent to develop a common core avionic system that has standard capabilities and interfaces, and contains the basic elements and functionality needed for any spacecraft. This common core will be scalable and tailored to specific missions. It will incorporate hardware and software from multiple vendors, and be upgradeable in order to infuse incremental capabilities and new technologies. It will maximize the use of reconfigurable open source software (e.g., Goddard Space Flight Center's (GSFC's) Core Flight Software (CFS)). Our long-term focus is on improving functionality, reliability, and autonomy, while reducing size, weight, and power. Where possible, we will leverage terrestrial commercial capabilities to drive down development and sustaining costs. We will select promising technologies for evaluation, compare them in an objective manner, and mature them to be available for future programs. The remainder of this paper describes our approach, technical areas of emphasis, integrated test experience and results as of mid-2014, and future plans. As a part of the AES Program, we are encouraged to set aggressive goals and fall short if necessary, rather than to set our sights too low. We are also asked to emphasize providing our personnel with hands-on experience in development, integration, and testing. That we have embraced both of these philosophies will be evident in the descriptions below.

Goforth, Montgomery B.