Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

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 613 records · Page 34

Absence of first-order unbinding transitions of fluid and polymerized membranes

Unbinding transitions of fluid and polymerized membranes are studied by renormalization-group (RG) methods. Two different RG schemes are used and found to give rather consistent results. The fixed-point structure of both RG's exhibits a complex behavior as a function of the decay exponent tau for the fluctuation-induced interaction of the membranes. For tau greater than tau(S2) interacting membranes can undergo first-order transitions even in the strong-fluctuation regime. These estimates for tau(S2) imply, however, that both fluid and polymerized membranes unbind in a continuous way in the absence of lateral tension.

Grotehans, Stefan↗

An Experimental Investigation of Incompressible Richtmyer-Meshkov Instability

Richtmyer-Meshkov (RM) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and eventually become a turbulent flow. It is closely related to Rayleigh-Taylor instability, which is the instability of a planar interface undergoing constant acceleration, such as caused by the suspension of a heavy fluid over a lighter one in the earth's gravitational field. Like the well-known Kelvin-Helmholtz instability, RM instability is a fundamental hydrodynamic instability which exhibits many of the nonlinear complexities that transform simple initial conditions into a complex turbulent flow. Furthermore, the simplicity of RM instability (in that it requires very few defining parameters), and the fact that it can be generated in a closed container, makes it an excellent test bed to study nonlinear stability theory as well as turbulent transport in a heterogeneous system. However, the fact that RM instability involves fluids of unequal densities which experience negligible gravitational force, except during the impulsive acceleration, requires RM instability experiments to be carried out under conditions of microgravity. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The impulsive acceleration is generated by bouncing a rectangular tank containing two different density liquids off a retractable vertical spring. The initial perturbation is produced prior to release by oscillating the tank in the horizontal direction to produce a standing wave. The instability evolves in microgravity as the tank travels up and then down the vertical rails of a drop tower until hitting a shock absorber at the bottom. Planar Laser Induced Fluorescence (PLIF) is employed to visualize the flow. PLIF images are captured by a video camera that travels with the tank. Figure 1 is as sequence of images showing the development of the instability from the initial sinusoidal disturbance far into the nonlinear regime which is characterized by the appearance of mushroom structures resulting from the coalescence of baroclinic vorticity produced by the impulsive acceleration. At later times in this sequence the vortex cores are observed to become unstable showing the beginnings of the transition to turbulence in this flow. The amplitude of the growing disturbance after the impulsive acceleration is measured and found to agree well with theoretical predictions. The effects of Reynolds number (based on circulation) on the development of the vortices and the transition to turbulence are also determined.

Jacobs, J. W.↗

Experimental evaluation of a pilot multinozzle-duct apparatus

A pilot multinozzle and duct were tested at ambient enthalpy to evaluate the suitability of such apparatus for testing thermal protection system panels mounted in the sidewalls of the duct downstream of the nozzle array. The flow field in the duct was complex: effects of wakes and shock waves from the nozzle dominated the flow field; the wakes continually mixed with the surrounding fluid; the boundary layer on the sidewalls of the duct was nonuniform; and near the exit of the duct the sidewall pressure variation was as much as 8.5 percent about the mean wall pressure. Starting loads on the duct walls were higher than those of a similar conventional nozzle and duct. It was concluded that the multinozzle-duct apparatus was not suitable for testing TPS panels, although the design and flow-field information should be of interest to designers of high-energy gasdynamic lasers.

Puster, R. L.↗

Thermal environment around the Space Shuttle with hot-gas jets for ice suppression

To prevent prelaunch ice formation on the external tank of the Space Shuttle, the final selected approach is to heat the surrounding air with vertical hot-gas jets located at the launch pad. This approach was considerably more cost-effective than other ice suppression methods considered, although its feasibility was not easily discernible due to the complex flowfield interactions. This paper describes how the use of vertical jets was first evaluated with the aid of a computational fluid dynamics (CFD) technique. An existing general-purpose CFD program (PHOENICS) was used to predict the thermal environment around the Shuttle under various jet configurations and wind conditions. The program accounts for effects of buoyancy, turbulence, and structural obstructions in the flowfield. The computed results showed physically plausible and consistent trends. The wind wake effects were found to be significant, and normally resulted in higher temperatures on the leeward side of the tank. High wind conditions were found to be more severe than calm wind conditions. The use of four jets with two different temperatures was identified as a promising option in which the air temperatures were raised sufficiently to prevent ice formation on the external tank, without excessive increase in Orbiter surface temperature. The use of the numerical model also facilitated the selection of test configurations and a test matrix for verifying the approach. Selected results of an experimental verification (by using a 2 percent scale model in a wind tunnel) are also presented.

Singhal, A. K.↗

Numerical simulation of single and multiple jet interaction within a combustor

Single and multiple jet interaction in a combustor is simulated numerically to determine the effects of various k-epsilon models and boundary conditions as part of an attempt to understand the physics of the flow field inside the preburner of the Space Shuttle Main Engine. Due to the complexity of the physical and chemical phenomena involved, the model is simplified to two dimensions and a single-component nonreacting fluid. The computer code used is the modified TEACH code, which is a steady state isothermal fluid code using the SIMPLE algorithm. The physical validity and numerical accuracy of the model are established through comparisons with experimentally determined isothermal flow fields in a combustor.

Bai, S. Don↗

Hypersonic propulsion: Status and challenge

Scientists in the U.S. are again focusing on the challenge of hypersonic flight with the proposed National Aerospace Plane (NASP). This renewed interest has led to an expansion of research related to high speed airbreathing propulsion, in particular, the supersonic combustion ramjet, or scramjet. The history is briefly traced of scramjet research in the U.S., with emphasis on NASA sponsored efforts, from the Hypersonic Research Engine (HRE) to the current status of today's airframe integrated scramjets. The challenges of scramjet technology development from takeover to orbital speeds are outlined. Existing scramjet test facilities such as NASA Langley's Scramjet Test Complex as well as new high Mach number pulse facilities are discussed. The important partnership role of experimental methods and computational fluid dynamics is emphasized for the successful design of single stage to orbit vehicles.

Guy, R. Wayne↗

Gravity Effects in Diffusive Coarsening of Bubble Lattices: von Neumann's Law

von Neumann modelled the evolution of two-dimensional soap froths as a purely diffusive phenomenon; the area growth of a given cell was found to depend only on the geometry of the bubble lattice. In the model, hexagons are stable, pentagons shrink and heptagons grow. The simplest equivalent to the area growth law is / approximately t(sub beta). The result depends on assuming (1) an incompressible gas; (2) bubble walls which meet at 120 deg and (3) constant wall thickness and curvature. Each assumption is borne out in experiments except the last one: bubble wall thickness between connecting cells varies in unit gravity because of gravity drainage. The bottom part of the soap membrane is thickened, the top part is thinned, such that gas diffusion across the membrane shows a complex dependence on gravity. As a result, experimental tests of von Neumann's law have been influenced by effects of gravity; fluid behavior along cell borders can give non-uniform wall thicknesses and thus alter the effective area and gas diffusion rates between adjacent bubbles. For area plotted as a function of time, Glazier (J.A. Glazier, S.P. Gross, and I. Stavans, Phys. Rev. A. 36, 306 (1987); J. Stavans, J.A, Glazier, Phys. Rev. Lett. 62, 1318 (1989).) suggest that in some cases their failure to observe von Neumann's predicted growth exponent ((sup beta)theor(sup =1; beta)exp(sup =0.70 + 0.10)) may have been the result of such "fluid drainage onto the lower glass plate". Additional experiments which varied plate spacing gave different beta exponents in a fashion consistent with this suggestion. During preliminary long duration experiments (approximately 100 h) aboard Spacelab-J, a low-gravity test of froth coarsening has examined (1) power law scaling of von Neumann's law (beta values) in the appropriate diffusive limits; (2) new bubble lattice dynamics such as greater fluid wetting behavior on froth membranes in low gravity; and (3) explicit relations for the gravity dependence of the second moment (or disorder parameter) governing the geometric spread in cell-sidedness around the mean of perfect hexagonal filling. By reducing the gravity-induced distortion in lattice wall thickness, the diffusion-limited regime of bubble coarsening becomes available for performing critical tests of network dynamics.

Noever, David A.↗

The SHIVA Project: Spaceflight Holography in a Virtual Apparatus

The poster summarizes theoretical and experimental concepts used in the development of the proposed NASA flight experiment SHIVA (Spaceflight Holography in a Virtual Apparatus). SHIVA exploits a unique holography-based, diagnostics tools to understand the behavior of small particles subjected to transient accelerations. Flight experiment protocols and apparatus will test model equations, characterize the acceleration environment and other microgravity phenomena. The primary objective of SHIVA is to enhance the current understanding of complex dynamics of small particles subjected to transient microgravity conditions. Existing theory did not fully describe the movement of particles in fluids in the microgravity environment. Results from recent ground-based experiments and comparison with the model recently developed by Rangel and Coimbra are presented.

Jim Trolinger↗

Monte Carlo Methodology Serves Up a Software Success

Widely used for the modeling of gas flows through the computation of the motion and collisions of representative molecules, the Direct Simulation Monte Carlo method has become the gold standard for producing research and engineering predictions in the field of rarefied gas dynamics. Direct Simulation Monte Carlo was first introduced in the early 1960s by Dr. Graeme Bird, a professor at the University of Sydney, Australia. It has since proved to be a valuable tool to the aerospace and defense industries in providing design and operational support data, as well as flight data analysis. In 2002, NASA brought to the forefront a software product that maintains the same basic physics formulation of Dr. Bird's method, but provides effective modeling of complex, three-dimensional, real vehicle simulations and parallel processing capabilities to handle additional computational requirements, especially in areas where computational fluid dynamics (CFD) is not applicable. NASA's Direct Simulation Monte Carlo Analysis Code (DAC) software package is now considered the Agency s premier high-fidelity simulation tool for predicting vehicle aerodynamics and aerothermodynamic environments in rarified, or low-density, gas flows.

Source record↗

AirShow 1.0 CFD Software Users' Guide

AirShow is visualization post-processing software for Computational Fluid Dynamics (CFD). Upon reading binary PLOT3D grid and solution files into AirShow, the engineer can quickly see how hundreds of complex 3-D structured blocks are arranged and numbered. Additionally, chosen grid planes can be displayed and colored according to various aerodynamic flow quantities such as Mach number and pressure. The user may interactively rotate and translate the graphical objects using the mouse. The software source code was written in cross-platform Java, C++, and OpenGL, and runs on Unix, Linux, and Windows. The graphical user interface (GUI) was written using Java Swing. Java also provides multiple synchronized threads. The Java Native Interface (JNI) provides a bridge between the Java code and the C++ code where the PLOT3D files are read, the OpenGL graphics are rendered, and numerical calculations are performed. AirShow is easy to learn and simple to use. The source code is available for free from the NASA Technology Transfer and Partnership Office.

Mohler, Stanley R., Jr.↗

Adjustable Autonomy Testbed

The Adjustable Autonomy Testbed (AAT) is a simulation-based testbed located in the Intelligent Systems Laboratory in the Automation, Robotics and Simulation Division at NASA Johnson Space Center. The purpose of the testbed is to support evaluation and validation of prototypes of adjustable autonomous agent software for control and fault management for complex systems. The AA T project has developed prototype adjustable autonomous agent software and human interfaces for cooperative fault management. This software builds on current autonomous agent technology by altering the architecture, components and interfaces for effective teamwork between autonomous systems and human experts. Autonomous agents include a planner, flexible executive, low level control and deductive model-based fault isolation. Adjustable autonomy is intended to increase the flexibility and effectiveness of fault management with an autonomous system. The test domain for this work is control of advanced life support systems for habitats for planetary exploration. The CONFIG hybrid discrete event simulation environment provides flexible and dynamically reconfigurable models of the behavior of components and fluids in the life support systems. Both discrete event and continuous (discrete time) simulation are supported, and flows and pressures are computed globally. This provides fast dynamic simulations of interacting hardware systems in closed loops that can be reconfigured during operations scenarios, producing complex cascading effects of operations and failures. Current object-oriented model libraries support modeling of fluid systems, and models have been developed of physico-chemical and biological subsystems for processing advanced life support gases. In FY01, water recovery system models will be developed.

Malin, Jane T.↗

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗

Domain decomposition for aerodynamic and aeroacoustic analyses, and optimization

The overarching theme was the domain decomposition, which intended to improve the numerical solution technique for the partial differential equations at hand; in the present study, those that governed either the fluid flow, or the aeroacoustic wave propagation, or the sensitivity analysis for a gradient-based optimization. The role of the domain decomposition extended beyond the original impetus of discretizing geometrical complex regions or writing modular software for distributed-hardware computers. It induced function-space decompositions and operator decompositions that offered the valuable property of near independence of operator evaluation tasks. The objectives have gravitated about the extensions and implementations of either the previously developed or concurrently being developed methodologies: (1) aerodynamic sensitivity analysis with domain decomposition (SADD); (2) computational aeroacoustics of cavities; and (3) dynamic, multibody computational fluid dynamics using unstructured meshes.

Baysal, Oktay↗

Computational fluid dynamics - The coming revolution

The development of aerodynamic theory is traced from the days of Aristotle to the present, with the next stage in computational fluid dynamics dependent on superspeed computers for flow calculations. Additional attention is given to the history of numerical methods inherent in writing computer codes applicable to viscous and inviscid analyses for complex configurations. The advent of the superconducting Josephson junction is noted to place configurational demands on computer design to avoid limitations imposed by the speed of light, and a Japanese projection of a computer capable of several hundred billion operations/sec is mentioned. The NASA Numerical Aerodynamic Simulator is described, showing capabilities of a billion operations/sec with a memory of 240 million words using existing technology. Near-term advances in fluid dynamics are discussed.

Graves, R. A., Jr.↗

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

Liquid rocket combustor computer code development

The Advanced Rocket Injector/Combustor Code (ARICC) that has been developed to model the complete chemical/fluid/thermal processes occurring inside rocket combustion chambers are highlighted. The code, derived from the CONCHAS-SPRAY code originally developed at Los Alamos National Laboratory incorporates powerful features such as the ability to model complex injector combustion chamber geometries, Lagrangian tracking of droplets, full chemical equilibrium and kinetic reactions for multiple species, a fractional volume of fluid (VOF) description of liquid jet injection in addition to the gaseous phase fluid dynamics, and turbulent mass, energy, and momentum transport. Atomization and droplet dynamic models from earlier generation codes are transplated into the present code. Currently, ARICC is specialized for liquid oxygen/hydrogen propellants, although other fuel/oxidizer pairs can be easily substituted.

Liang, P. Y.↗

Invariance of Hypersonic Normal Force Coefficients with Reynolds Number and Determination of Inviscid Wave Drag from Laminar Experimental Results

Observations have been made and reported that the experimental normal force coefficients at a constant angle of attack were constant with a variation of more than 2 orders of magnitude of Reynolds number at a free-stream Mach number M(sub infinity) of 8.00 and more than 1 order of magnitude variation at M(sub infinity) = 6.00 on the same body-wing hypersonic cruise configuration. These data were recorded under laminar, transitional, and turbulent boundary layer conditions with both hot-wall and cold-wall models. This report presents experimental data on 25 configurations of 17 models of both simple and complex geometry taken at M(sub infinity) = 6.00, 6.86, and 8.00 in 4 different hypersonic facilities. Aerodynamic calculations were made by computational fluid dynamics (CID) and engineering methods to analyze these data. The conclusions were that the normal force coefficients at a given altitude are constant with Reynolds numbers at hypersonic speeds and that the axial force coefficients recorded under laminar boundary-layer conditions at several Reynolds numbers may be plotted against the laminar parameter (the reciprocal of the Reynolds number to the one-half power) and extrapolated to the ordinate axis to determine the inviscid-wave-drag coefficient at the intercept.

Hawkins, Richard↗

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗