Drag coefficients for free molecule flow in the velocity range 7-37 km/sec
Monoenergetic nitrogen free molecule beam impingement on solid surface, calculating satellite drag coefficients from momentum transfer measurements
SEARCH · Engineering Papers
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.
Monoenergetic nitrogen free molecule beam impingement on solid surface, calculating satellite drag coefficients from momentum transfer measurements
This project is aimed at using a rotating magnetic field (RMF) to control fluid flow and transport during directional solidification of elemental and compound melts. Microgravity experiments have demonstrated that small amounts of residual acceleration of less than a micro-g can initiate and prolong fluid flow, particularly when there is a static component of the field perpendicular to the liquid solid interface. Thus a true diffusion boundary layer is not formed, and it becomes difficult to verify theories of solidification or to achieve diffusion controlled solidification. The RMF superimposes a stirring effect on an electrically conducting liquid, and with appropriate field strengths and frequencies, controlled transport of material through a liquid column can be obtained. As diffusion conditions are precluded and complete mixing conditions prevail, the technique is appropriate for traveling solvent zone or float zone growth methods in which the overall composition of the liquid can be maintained throughout the growth experiment. Crystals grown by RMF techniques in microgravity in previous, unrelated missions have shown exceptional properties. The objective of the project is two-fold, namely (1) using numerical modeling to simulate the behavior of a solvent zone with applied thermal boundary conditions and demonstrate the effects of decreasing gravity levels, or an increasing applied RMF, or both, and (2) to grow elements and II-VI compounds from traveling solvent zones both with and without applied RMFs, and to determine objectively how well the modeling predicts solidification parameters. Numerical modeling has demonstrated that, in the growth of CdTe from a tellurium solution, a rotating magnetic field can advantageously modify the shape of the liquid solid interface such that the interface is convex as seen from the liquid. Under such circumstances, the defect structure is reduced as any defects which are formed tend to grow out and not propagate. The flow of liquid, however, is complex due to the competing flow induced by the rotating magnetic field and the buoyancy driven convection. When the acceleration forces are reduced to one thousandth of gravity, the flow pattern is much simplified and well controlled material transport through the solvent zone can be readily achieved. Triple axis diffractometry and x-ray synchrotron topography have demonstrated that there is no significant improvement in crystal quality for HgCdTe grown on earth from a tellurium solution when a rotating magnetic field is applied. However, modeling shows that the flow in microgravity with a rotating magnetic field would produce a superior product.
Several studies have shown a significant increase in drag on a distribution of solid spherical particles within a fluid with increasing particle volume fraction. As a result, many empirical drag laws accounting for the dependence on the Reynolds number and volume fraction can be found in the literature. This study investigates the possibility of a similar effect of the particle volume fraction on the mean hydrodynamic lift force on randomly distributed spherical particles in a linear shear flow. Particle-resolved direct numerical simulations are performed to evaluate the mean lift force, and the results are compared with the case of an isolated particle in a linear shear flow for the same Reynolds number and shear rate. The mean lift force acting on the particles appears to remain nearly the same as that on an isolated particle. However, due to the influence of neighboring particles, there is a substantial force variation in transverse directions on each individual particle, whose magnitude is comparable to the mean drag force. The distribution of drag force in a linear shear flow is shown to be nearly the same as in a uniform flow at the same volume fraction and Reynolds number. A simple stochastic model based on a Gaussian distribution is presented for the lift force variation, and its performance is compared to the prediction of the deterministic pairwise interaction extended point-particle model.
This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.
As the first step of the study of the combustion control of solid propellants by electrical discharges, the effects of an arc discharge, which flows along the burning surface, on the burning rate and on the increase of enthalpy of the combustion product were investigated. For specially devised composite propellants, which are composed of Al and Teflon powders, it was shown that the combination can be controlled by an arc discharge; the combustion continues when the arc discharge is applied and is interrupted when the arc discharge breaks. In the present investigation, it was also shown that an arc discharge coupled with a high-frequency electrical discharge has potential as an effective ignition method for solid propellants. For the application of this type of combustion control to an ignitor for a solid propellant rocket motor or to a control rocket motor, this method lacks flexibility in the configuration scale and needs relatively high electric power at the present stage.
A stability analysis based on Arnol'd-Liapunov arguments in a gauge-variable formalism is applied to a generalized form of propagating modons. The method was previously applied to the special case of Stern's modon in a quiescent background flow. The analysis presented here shows that propagating-modon solutions in the shallow-water equations are stable to small-amplitude perturbations, regardless of the sign of their translational speed, as long as this is much smaller in magnitude than the solid-rotation speed of the earth. Sufficient conditions for stability are that the given flow be quasigeostrophic, i.e., that the Rossby number be small, and that the associated Froude number be one order of magnitude smaller than the Rossby number.
Planar flow visualizations were obtained in a wind tunnel test in the NASA Langley Research Center’s Unitary Plan Wind Tunnel using the laser-light-sheet method. This method uses a laser to illuminate fine particles generated in the wind tunnel to visualize flow structures. The test article was designed to simulate the separation of the two solid rocket boosters (SRBs) from the core stage of the NASA Space Launch System (SLS) at Mach 4 using a scale model. The test was run on of the SLS Block 1B Cargo (27005) configuration and the SLS Block 1B Crew (28005) configuration. Planar flow visualization was obtained only on the crew configuration. Air at pressures up to 1500 psi was used to simulate plumes from the booster separation motors (BSMs) located at the nose, and aft skirt of the two boosters. The facility free stream was seeded with water vapor, which condensed and froze into small ice crystals in the tunnel nozzle expansion. A continuous wave green (532 nm) laser sheet was used to illuminate the ice crystals, and the resulting Mie-scattered light was collected with a camera. The resulting images clearly identify shock waves and other flow features including BSM plume shapes. Measurements were acquired for different BSM pressures and booster separation locations.
A previously developed computer program was used to model two transient hybrid combustion processes involving tubes of solid Plexiglas. In the first study, representing combustion of a hybrid rocket, the oxidizing gas was oxygen, and calculations were continued sufficiently long to obtain steady-state values. Systematic variations were made in reaction rate constant, mass flow rate, and pressure, alternatively using constant and temperature dependent regression rate models for the fuel surface. Consistent results were obtained, as is evidenced by the values for the mass function of the reaction product and the flame temperature, for which plots are supplied. In the second study, fire initiation in a duct was studied, with an air mixture as the oxidizing gas. It was demonstrated that a satisfactory flame spread mechanism could be reproduced on the computer. In both of the above applications, the general, transient, two-dimensional conservation equations were represented, together with chemical reactions, solid-fuel interface conditions, and heat conduction in the solid fuel.
The Liquid Sheet Radiator (LSR) is an external flow radiator that uses a triangular-shaped flowing liquid sheet as the radiating surface. It has potentially much lower mass than solid wall radiators such as pumped loop and heat pipe radiators, along with being nearly immune to micrometeoroid penetration. The LSR has an added advantage of simplicity. Surface tension causes a thin (100-300 microns) liquid sheet to coalesce to a point, causing the sheet flow to have a triangular shape. Such a triangular sheet is desirable since it allows for simple collection of the flow at a single point. A major problem for all external flow radiators is the requirement that the working fluid be of very low (approx. 10(sup -8) torr) vapor pressure to keep evaporative losses low. As a result, working fluids are limited to certain oils (such as used in diffusion pumps) for low temperatures (300-400 K) and liquid metals for higher temperatures. Previous research on the LSR has been directed at understanding the fluid mechanics of thin sheet flows and assessing the stability of such flows, especially with regard to the formation of holes in the sheet. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. The latest research has been directed at determining the emittance of thin sheet flows. The emittance was calculated from spectral transmittance data for the Dow Corning 705 silicone oil. By experimentally setting up a sheet flow, the emittance was also determined as a function of measurable quantities, most importantly, the temperature drop between the top of the sheet and the temperature at the coalescence point of the sheet. Temperature fluctuations upstream of the liquid sheet were a potential problem in the analysis and were investigated.
The development of aerodynamic theories of streaming motions around bodies with unsteady vortical and entropic disturbances is reviewed. The basic concepts associated with such motions, their interaction with solid boundaries and their noise generating mechanisms are described. The theory was first developed in the approximation wherein the unsteady flow is linearized about a uniform mean lfow. This approach has been extensively developed and used in aeroelastic and aeroacoustic calculations. The theory was recently extended to account for the effect of distortion of the incident disturbances by the nonuniform mean flow around the body. This effect is found to have a significant influence on the unsteady aerodynamic force along the body surface and the sound radiated in the far field. Finally, the nonlinear characteristics of unsteady transonic flows are reviewed and recent results of linear and nonlinear computations are presented.
Future NASA architectures have baselined cryogenic propulsion systems as well as cryogenic fluid management to support lunar missions and ultimately to support future missions to Mars. These missions will require chilling hardware down prior to engine restart as well as chilling lines and tanks prior to transferring and refueling these propulsion elements in orbit. In lieu of expensive tests conducted on-orbit, accurate predictive computational models of these chilldown processes can be used to reduce system and propellant mass as well as mission risk. To gain confidence in these computational models, appropriate anchoring and validation to experimental data in a relevant environment needs to be performed. Recent ground and sub-orbital flight experiments conducted by the Japan Aerospace Exploration Agency (JAXA) investigated chilldown of a complex channel resembling a turbopump bearing cavity at low flow rates. This work presents Computational Fluid Dynamics (CFD) model development of the chilldown experiment employing two-phase flow boiling models available in commercial CFD software STAR-CCM+ using Volume of Fluid (VOF) and the traditional Euler-Euler multiphase flow solvers. Comparisons of the numerical and experimental results under normal and low-gravity conditions are presented. An assessment of solid wall temperatures and phase distribution yielded important insights into multiphase solver choice, dependence on gravity environment, and challenges associated with cryogenic flow boiling prediction and validation.
The traditional Zn/MnO 2 battery has attracted great interest due to its low cost, high safety, high output voltage, and environmental friendliness. However, it remains a big challenge to achieve long-term stability, mainly owing to the poor reversibility of the cathode reaction. Different from previous studies where the cathode redox reaction of MnO 2 /MnOOH is in solid state with limited reversibility, in this work a new aqueous rechargeable Zn/MnO 2 flow battery is constructed with dissolution–precipitation reactions in both cathodes (Mn 2+ /MnO 2 ) and anodes (Zn 2+ /Zn), which allow mixing of anolyte and catholyte into only one electrolyte and remove the requirement for an ion selective membrane for cost reduction. Impressively, this new battery exhibits a high discharge voltage of ≈1.78 V, good rate capability (10C discharge), and excellent cycling stability (1000 cycles without decay) at the areal capacity ranging from 0.5 to 2 mAh cm -2 . More importantly, this battery can be readily enlarged to a bench scale flow cell of 1.2 Ah with good capacity retention of 89.7% at the 500th cycle, displaying great potential for large-scale energy storage.
A novel method for treating certain troublesome boundary conditions in the numerical solution of time-dependent incompressible viscous flow problems is presented. This new method is developed on the basis of an integral representation for the velocity vector which contains the entire kinematics of the problem, including the boundary conditions of concern. It is shown that for the exterior flow problem the free-stream condition is satisfied at infinity exactly and the need to treat a far-field condition is removed by the use of the integral representation. The distribution of a nonvelocity variable on the solid boundary - i.e., the 'extraneous' boundary condition needed for both the exterior and the interior flows - is shown to be governed by the kinematics of the problem.
Sublimating chemicals offer accurate, low-cost way of indicating laminarto-turbulent flow transisions on surfaces of aircraft. Aerodynamic surfaces coated with thin film of such volatile chemical solids as naphthalene, diphenyl, acenaphthene, or fluorene. Film sublimes rapidly because of high local shear stress and heat transfer in boundary layer. Coating appears white in regions where chemical remained on surface indicating laminar flow; regions where chemical disappeared indicate turbulent flow.
In this work, we numerically investigate the transient combustion of a premixed gas mixture in a narrow, perfectly-insulated, catalytic channel which can represent an interior channel of a catalytic monolith. The model assumes a quasi-steady gas-phase and a transient, thermally thin solid phase. The gas phase is one-dimensional, but it does account for heat and mass transfer in a direction perpendicular to the flow via appropriate heat and mass transfer coefficients. The model neglects axial conduction in both the gas and in the solid. The model includes both detailed gas-phase reactions and catalytic surface reactions. The reactants modeled so far include lean mixtures of dry CO and CO/H2 mixtures, with pure oxygen as the oxidizer. The results include transient computations of light-off and system response to inlet condition variations. In some cases, the model predicts two different steady-state solutions depending on whether the channel is initially hot or cold. Additionally, the model suggests that the catalytic ignition of CO/O2 mixtures is extremely sensitive to small variations of inlet equivalence ratios and parts per million levels of H2.
We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.
The effects of surface and gas-phase radiation on the rate and the structure of laminar flame spread over thin fuels are investigated using a flame-spread model which consists of the continuity, momentum, species, and energy equations in the gas and the continuity and energy equations in the solid. Numerical calculations, complemented by scaling arguments, show that, at high velocities of the oxidizer flow, radiation effects are unimportant; the spread rate decreases with increasing opposing velocity due to finite-rate gas-phase kinetics. However, radiation becomes progressively important when the opposing velocity is below a certain value: the flame cools, shrinks in size, and its spread rate falls sharply with decreasing opposing velocity.
The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.