Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gas dynamics”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Development of Flow and Heat Transfer Models for the Carbon Fiber Rope in Nozzle Joints of the Space Shuttle Reusable Solid Rocket Motor

Methodologies have been developed for modeling both gas dynamics and heat transfer inside the carbon fiber rope (CFR) for applications in the space shuttle reusable solid rocket motor joints. Specifically, the CFR is modeled using an equivalent rectangular duct with a cross-section area, friction factor and heat transfer coefficient such that this duct has the same amount of mass flow rate, pressure drop, and heat transfer rate as the CFR. An equation for the friction factor is derived based on the Darcy-Forschheimer law and the heat transfer coefficient is obtained from pipe flow correlations. The pressure, temperature and velocity of the gas inside the CFR are calculated using the one-dimensional Navier-Stokes equations. Various subscale tests, both cold flow and hot flow, have been carried out to validate and refine this CFR model. In particular, the following three types of testing were used: (1) cold flow in a RSRM nozzle-to-case joint geometry, (2) cold flow in a RSRM nozzle joint No. 2 geometry, and (3) hot flow in a RSRM nozzle joint environment simulator. The predicted pressure and temperature history are compared with experimental measurements. The effects of various input parameters for the model are discussed in detail.

Wang, Q.↗
NASA NTRS1987

Development and application of unified algorithms for problems in computational science

A framework is presented for developing computationally unified numerical algorithms for solving nonlinear equations that arise in modeling various problems in mathematical physics. The concept of computational unification is an attempt to encompass efficient solution procedures for computing various nonlinear phenomena that may occur in a given problem. For example, in Computational Fluid Dynamics (CFD), a unified algorithm will be one that allows for solutions to subsonic (elliptic), transonic (mixed elliptic-hyperbolic), and supersonic (hyperbolic) flows for both steady and unsteady problems. The objectives are: development of superior unified algorithms emphasizing accuracy and efficiency aspects; development of codes based on selected algorithms leading to validation; application of mature codes to realistic problems; and extension/application of CFD-based algorithms to problems in other areas of mathematical physics. The ultimate objective is to achieve integration of multidisciplinary technologies to enhance synergism in the design process through computational simulation. Specific unified algorithms for a hierarchy of gas dynamics equations and their applications to two other areas: electromagnetic scattering, and laser-materials interaction accounting for melting.

Shankar, Vijaya↗
NASA NTRS1986

Detonation wave compression in gas turbines

A study was made of the concept of augmenting the performance of low pressure ratio gas turbines by detonation wave compression of part of the flow. The concept exploits the constant volume heat release of detonation waves to increase the efficiency of the Brayton cycle. In the models studied, a fraction of the compressor output was channeled into detonation ducts where it was processed by transient transverse detonation waves. Gas dynamic studies determined the maximum cycling frequency of detonation ducts, proved that upstream propagation of pressure pulses represented no problems and determined the variations of detonation duct output with time. Mixing and wave compression were used to recombine the combustor and detonation duct flows and a concept for a spiral collector to further smooth the pressure and temperature pulses was presented as an optional component. The best performance was obtained with a single firing of the ducts so that the flow could be re-established before the next detonation was initiated. At the optimum conditions of maximum frequency of the detonation ducts, the gas turbine efficiency was found to be 45 percent while that of a corresponding pressure ratio 5 conventional gas turbine was only 26%. Comparable improvements in specific fuel consumption data were found for gas turbines operating as jet engines, turbofans, and shaft output machines. Direct use of the detonation duct output for jet propulsion proved unsatisfactory. Careful analysis of the models of the fluid flow phenomena led to the conclusion that even more elaborate calculations would not diminish the uncertainties in the analysis of the system. Feasibility of the concept to work as an engine now requires validation in an engineering laboratory experiment.

Wortman, A.↗
DOE OSTI2020

A Novel use of Direct Simulation Monte-Carlo to Model Dynamics of COVID-19 Pandemic Spread

In this report, we evaluate a novel method for modeling the spread of COVID-19 pandemic. In this new approach we leverage methods and algorithms developed for fully-kinetic plasma physics simulations using Particle-In-Cell (PIC) Direct Simulation Monte-Carlo (DSMC) models. This approach then leverages Sandia-unique simulation capabilities, and High-Performance Computer (HPC) resources and expertise in particle-particle interactions using stochastic processes. Our hypothesis is that this approach would provide a more efficient platform with assumptions based on physical data that would then enable the user to assess the impact of mitigation strategies and forecast different phases of infection. This work addresses key scientific questions related to the assumptions this new approach must make to model the interactions of people using algorithms typically used for modeling particle interactions in physics codes (kinetic plasma, gas dynamics). The model developed uses rational/physical inputs while also providing critical insight; the results could serve as inputs to, or alternatives for, existing models. The model work presented was developed over a four-week time frame, thus far showing promising results and many ways in which this model/approach could be improved. This work is aimed at providing a proof-of-concept for this new pandemic modeling approach, which could have an immediate impact on the COVID-19 pandemic modeling, while laying a basis to model future pandemic scenarios in a manner that is timely and efficient. Additionally, this new approach provides new visualization tools to help epidemiologists comprehend and articulate the spread of this and other pandemics as well as a more general tool to determine key parameters needed in order to better predict pandemic modeling in the future. In the report we describe our model for pandemic modeling, apply this model to COVID-19 data for New York City (NYC), assess model sensitivities to different inputs and parameters and , finally, propagate the model forward under different conditions to assess the effects of mitigation and associated timing. Finally, our approach will help understand the role of asymptomatic cases, and could be extended to elucidate the role of recovered individuals in the second round of the infection, which is currently being ignored.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1972

Interrelated structures of the transport shock and collisional relaxation layer in a multitemperature, multilevel ionized gas

The gas dynamic structures of the transport shock and the downstream collisional relaxation layer are evaluated for partially ionized monatomic gases. Elastic and inelastic collisional nonequilibrium effects are taken into consideration. Three electronic levels are accounted for in the microscopic model of the atom. Nonequilibrium processes with respect to population of levels and species plus temperature are considered. By using an asymptotic technique the shock morphology is found on a continuum flow basis. The asymptotic procedure gives two distinct layers in which the nonequilibrium effects to be considered are different. A transport shock appears as the inner solution to an outer collisional relaxation layer in which the gas reaches local equilibrium. A family of numerical examples is displayed for different flow regimes. Argon and helium models are used in these examples.

Vinolo, A. R.↗
DOE OSTI2026

Gauging the Impact of Cosmic-ray Feedback on the Stellar Initial Mass Function

Cosmic rays (CRs) drive ionization and influence gas dynamics in molecular clouds (MCs), potentially impacting the resulting star formation outcomes. Although previous simulations of individual star formation have included methods for CR transport (CRT), none have been large enough to resolve the stellar initial mass function (IMF). We conduct numerical simulations following the collapse of a 20,000 M ⊙ MC and the subsequent star formation including CRT, both with and without CRs accelerated by winds from the young massive stars, and compare against a non-CRT simulation. We show that after the first massive stars form, the cavity produced by feedback is more pronounced in the CRT simulations because the external CRs are able to propagate inward and compress the gas into higher-density structures. This increases the subsequent star formation in the cloud; by the end of the simulation, the star formation efficiency (SFE) in the CRT simulation including stellar wind CRs is 43% higher than the non-CRT simulation. The IMF is also top-heavy in comparison, with a slope above 1 M ⊙ that is shallower by ∼20%. These effects are also present in the simulation without wind-accelerated CRs, but they are not as pronounced; the SFE is only 16% higher than the non-CRT simulation, and the IMF high-mass slope is shallower by ∼10%. These results may explain some of the observed top-heavy IMFs, which typically occur in high-CR environments such as the Galactic center.

Fitz Axen, Margot [Oak Ridge National Laboratory (↗
NASA NTRS2004

History of Thermally Processed Solids in the Protoplanetary Disk: Reconciling Theoretical Models and Meteoritical Evidence

In this talk we assess theoretical models of the radial, temporal, and thermal evolution of nebula solids, and their ultimate accretion into planetesimals such as we see today, using meteorite evidence as a guide. Each class of chondrites contains a characteristic suite of chondrules and CAIs that may have formed over a period of several Myr during which planetesimals were accreting in the disk. Details of the various models for transient melting of chondrules and igneous CAIs will be left to others. However, high-temperature processes of different kinds evaporation, alteration, etc did affect these constituents and their environment over this time span. Here we describe evolutionary scenarios consistent with a large time gap between CAI and chondrule formation and the presence of distinctive suites of chondrules and CAIs in each chondrite class. Particle-gas dynamical processes transport particles of all relevant sizes (microns to many meters) within the nebula and affect their evolution in a variety of important ways. Turbulent radial diffusion spreads particles radially down their concentration gradients - as one example, it can prevent CAIs from being lost into the sun on several Myr timescales [1]. Vertical diffusion spreads the dense midplane particle layer, determining its volume density, which in turn affects the particle growth rate and even the dominant growth process [2-4]. Turbulent concentration selects aerodynamically sorted particles for orders-of-magnitude density enhancement, and is applicable to porous, fluffy particles of appropriate size as well as to solid chondrules [5]. Inward radial drift under gas drag brings a surprisingly large amount of material to regions where it evaporates; these evaporation fronts cause significant chemical modification of the nebula gas over a wide range of radii [6]. Radial transport by stellar winds can be important for small particles [7].

Cuzzi, J. N.↗
DOE OSTI2019

Detonation propagation for shock-driven, subsonic and supersonic confiner flow

In this work, we study the compressible flow dynamics of two-dimensional, steady detonation wave propagation in a high explosive (HE) confined by aluminium (Al) or stainless steel (SS), outside of which is an air layer. We examine how the thickness of the confinement affects the subsonic detonation driving zone structure (DDZ) and the detonation speed (D 0 ), demonstrating a strong dependence on whether the oblique shock-driven flow in the confiner is supersonic, as for SS, or subsonic, as for Al. A characteristic path analysis is used to examine the information flow from the material boundaries through the supersonic flow regions in both the HE and confiner that can impact the sonic surfaces bounding the subsonic flow regions. It is shown that the nature of gas-dynamic wave reflection off the SS–air or Al–air boundary can significantly influence the DDZ and D 0 .

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Macroscopic balance model for wave rotors

A mathematical model for multi-port wave rotors is described. The wave processes that effect energy exchange within the rotor passage are modeled using one-dimensional gas dynamics. Macroscopic mass and energy balances relate volume-averaged thermodynamic properties in the rotor passage control volume to the mass, momentum, and energy fluxes at the ports. Loss models account for entropy production in boundary layers and in separating flows caused by blade-blockage, incidence, and gradual opening and closing of rotor passages. The mathematical model provides a basis for predicting design-point wave rotor performance, port timing, and machine size. Model predictions are evaluated through comparisons with CFD calculations and three-port wave rotor experimental data. A four-port wave rotor design example is provided to demonstrate model applicability. The modeling approach is amenable to wave rotor optimization studies and rapid assessment of the trade-offs associated with integrating wave rotors into gas turbine engine systems.

Welch, Gerard E.↗
NASA NTRS1983

Numerical simulations of loops heated to solar flare temperatures. I - Gasdynamics. II - X-ray and UV spectroscopy

The NRL's Dynamic Flux Tube Model is used to numerically simulate the dynamic response of a coronal magnetic loop to an energy input of the order encountered in solar flares. The coronal plasma is heated by the deposition of flare energy at the top of the loop to more than 10 million K, yielding a conduction front that moves toward the chromosphere, where the plasma is heated by the large downward conductive flux and ablates upward to the coronal part of the loop at velocities of a few hundred km/sec. The conduction front simultaneously produces chromospheric ablation and compresses the material ahead of it. With the aid of compressional instabilities, the compressed plasma grows throughout the flare heating phase, presenting a possible source of the flare optical continuum emission which is correlated with soft X-ray radiation. The observational consequences of rapidly heated loop gas dynamic processes are discussed. In the second part of this presentation, the dynamical calculation results previously obtained are used to predict the spectral line intensities, profiles and wavelengths of several X-ray lines and the UV line of Fe XXI at 1354.1 A. Three different viewing orientations of the loop are considered.

Cheng, C.-C.↗
NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗
NASA NTRS2011

Effect of Inductive Coil Geometry and Current Sheet Trajectory of a Conical Theta Pinch Pulsed Inductive Plasma Accelerator

Results are presented demonstrating the e ect of inductive coil geometry and current sheet trajectory on the exhaust velocity of propellant in conical theta pinch pulsed induc- tive plasma accelerators. The electromagnetic coupling between the inductive coil of the accelerator and a plasma current sheet is simulated, substituting a conical copper frustum for the plasma. The variation of system inductance as a function of plasma position is obtained by displacing the simulated current sheet from the coil while measuring the total inductance of the coil. Four coils of differing geometries were employed, and the total inductance of each coil was measured as a function of the axial displacement of two sep- arate copper frusta both having the same cone angle and length as the coil but with one compressed to a smaller size relative to the coil. The measured relationship between total coil inductance and current sheet position closes a dynamical circuit model that is used to calculate the resulting current sheet velocity for various coil and current sheet con gura- tions. The results of this model, which neglects the pinching contribution to thrust, radial propellant con nement, and plume divergence, indicate that in a conical theta pinch ge- ometry current sheet pinching is detrimental to thruster performance, reducing the kinetic energy of the exhausting propellant by up to 50% (at the upper bound for the parameter range of the study). The decrease in exhaust velocity was larger for coils and simulated current sheets of smaller half cone angles. An upper bound for the pinching contribution to thrust is estimated for typical operating parameters. Measurements of coil inductance for three di erent current sheet pinching conditions are used to estimate the magnetic pressure as a function of current sheet radial compression. The gas-dynamic contribution to axial acceleration is also estimated and shown to not compensate for the decrease in axial electromagnetic acceleration that accompanies the radial compression of the plasma in conical theta pinches.

Hallock, Ashley K.↗
NASA NTRS1993

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗
NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗
NASA NTRS1959

Measurements of Local Heat Transfer and Pressure on Six 2-Inch-Diameter Blunt Bodies at a Mach Number of 4.95 and at Reynolds Numbers Per Foot up to 81 x 10(exp 6)

Measurements of the local heat transfer and pressure distribution have been made on six 2-inch-diameter, blunt, axially symmetric bodies in the Langley gas dynamics laboratory at a Mach number of 4.95 and at Reynolds numbers per foot up to 81 x 10(exp 6). During the investigation laminar flow was observed over a hemispherical-nosed body having a surface finish from 10 to 20 microinches at the highest test Reynolds number per foot (for this configuration) of 77.4 x 10(exp 6). Though it was repeatedly possible to measure completely laminar flow at this Reynolds number for the hemisphere, it was not possible to observe completely laminar flow on the flat-nosed body for similar conditions. The significance of this phenomenon is obscured by the observation that the effects of particle impacts on the surface in causing roughness were more pronounced on the flat-nosed body. For engineering purposes, a method developed by M. Richard Dennison while employed by Lockheed Aircraft Corporation appears to be a reasonable procedure for estimating turbulent heat transfer provided transition occurs at a forward location on the body. For rearward-transition locations, the method is much poorer for the hemispherical nose than for the flat nose. The pressures measured on the hemisphere agreed very well with those of the modified Newtonian theory, whereas the pressures on all other bodies, except on the flat-nosed body, were bracketed by modified Newtonian theory both with and without centrifugal forces. For the hemisphere, the stagnation-point velocity gradient agreed very well with Newtonian theory. The stagnation-point velocity gradient for the flat- nosed model was 0.31 of the value for the hemispherical-nosed model. If a Newtonian type of flow is assumed, the ratio 0.31 will be independent of Much number and real-gas effects.

Cooper, Morton↗
NASA NTRS1996

Wave Engine Topping Cycle Assessment

The performance benefits derived by topping a gas turbine engine with a wave engine are assessed. The wave engine is a wave rotor that produces shaft power by exploiting gas dynamic energy exchange and flow turning. The wave engine is added to the baseline turboshaft engine while keeping high-pressure-turbine inlet conditions, compressor pressure ratio, engine mass flow rate, and cooling flow fractions fixed. Related work has focused on topping with pressure-exchangers (i.e., wave rotors that provide pressure gain with zero net shaft power output); however, more energy can be added to a wave-engine-topped cycle leading to greater engine specific-power-enhancement The energy addition occurs at a lower pressure in the wave-engine-topped cycle; thus the specific-fuel-consumption-enhancement effected by ideal wave engine topping is slightly lower than that effected by ideal pressure-exchanger topping. At a component level, however, flow turning affords the wave engine a degree-of-freedom relative to the pressure-exchanger that enables a more efficient match with the baseline engine. In some cases, therefore, the SFC-enhancement by wave engine topping is greater than that by pressure-exchanger topping. An ideal wave-rotor-characteristic is used to identify key wave engine design parameters and to contrast the wave engine and pressure-exchanger topping approaches. An aerodynamic design procedure is described in which wave engine design-point performance levels are computed using a one-dimensional wave rotor model. Wave engines using various wave cycles are considered including two-port cycles with on-rotor combustion (valved-combustors) and reverse-flow and through-flow four-port cycles with heat addition in conventional burners. A through-flow wave cycle design with symmetric blading is used to assess engine performance benefits. The wave-engine-topped turboshaft engine produces 16% more power than does a pressure-exchanger-topped engine under the specified topping constraints. Positive and negative aspects of wave engine topping in gas turbine engines are identified.

Welch, Gerard E.↗
NASA NTRS1974

Effect of cyclic conditions on the dynamic oxidation of gas turbine superalloys

The effects of operating parameters of a dynamic apparatus used to study oxidation and thermal fatigue of gas turbine materials were studied. IN-100, TD-NiCr, and WI-52 were tested at a maximum temperature of 1,090 deg C. Heating time per cycle was varied from 1/20 hr to 10 hr. Minimum temperatures between heating cycles were room temperature, 430 deg, and 650 deg C. Cooling air velocities were zero, Mach 0.7, and Mach 1. Increasing the number of cycles for a given time at temperature increased weight loss. Thermal fatigue was related to number of cycles more than to time at temperature.

Johnston, J. R.↗