Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reacting flows”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Using AMReX's Embedded Boundaries to Support MFIX-Exa’s Geometry Capabilities

MFIX-Exa is an AMReX-based application used for the simulation of chemically reacting multiphase flows inside reactors. The reactor geometry is mostly a combination of primitive shapes such as cylinders, cubes and spheres of varying dimensions connected together with potentially a lot of dead space between them. Additionally, the reactors may have particle-laden jets introduced at different parts of their geometry to enhance reactant mixing. This talk presents how MFIX-Exa uses OpenSCAD’s csg format to allow users to define the embedded boundary (EB) geometry and how it supports specifying inflow through certain sections of the EB. The talk will describe how MFIX-Exa leverages the EB features of grid-pruning, restart files and redistribution. Finally, the talk goes over some capabilities which are currently not available in MFIX-Exa but planned for the future, such as moving parts of the reactor geometry.

embedded boundaries↗
DOE OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗
DOE OSTI2022

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Fluidic Ignition Detection

Ignition sensor, utilizing fluidic oscillator device and its ability to react to flow system's molecular weight and temperature changes, detects propellant ignition. Commercial applications include fields of fire protection equipment and fluidics.

Henry, C. S.↗
NASA NTRS1973

Non-reacting and equilibrium chemically reacting turbulent boundary-layer flows.

Turbulent boundary-layer flows of nonreacting gases are predicted for both internal (nozzle) and external flows. Effects of favorable pressure gradients on two eddy viscosity models were studied in rocket and hypervelocity wind tunnel flows. Nozzle flows of equilibrium air with stagnation temperatures up to 10,000 K were computed. Predictions of equilibrium nitrogen flows through hypervelocity nozzles were compared with experimental data. A slender spherically blunted cone was studied at 70,000 ft altitude and 19,000 ft/sec in the earth's atmosphere. Comparisons with available experimental data showed good agreement. A computer program was developed and fully documented during this investigation for use by interested individuals.

Anderson, E. C.↗
NASA NTRS1974

Measurement techniques for supersonic combustion testing

Measurement techniques used during nine years of intensive research in supersonic combustion are presented. The tests employed a number of combustor and injector configurations, which were directly connected to a Mach 3.2 supersonic nozzle, and both gaseous hydrogen and storable liquid fuels. Nominal plenum conditions for the arc-heated air were 3000 to 4500 R and 30 atm. Measurements included wall surface pressures, wall temperatures, heat flux, skin friction, instream cone-static and pitot pressures, and gas samples. Measurement and analysis techniques have been developed for determining total enthalpy of an airstream containing NO concentrations in excess of equilibrium using the two-throat technique, bulk combustion efficiency using steam calorimetry, radial gradients of fluid properties in the combustor exit plane using a stream tube method, and wall shear from measurements of heat transfer in a reacting supersonic flow using Reynolds analogy.

Orth, R. C.↗
NASA NTRS1977

Prediction of hydrodynamics and chemistry of confined turbulent methane-air flames with attention to formation of oxides of nitrogen

A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a tubular combustor is presented. A numerical procedure for the solution of the governing differential equations is described, and models for chemical equilibrium and chemical kinetics calculations are presented. The chemical equilibrium model is used to characterize the hydrocarbon reactions. The chemical kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor consists of a cylindrical duct of varying cross sections with concentric streams of gaseous fuel and air entering the duct at one end. Four sample cases with specified inlet and boundary conditions are considered, and the results are discussed

Elghobashi, S.↗
NASA NTRS1982

Symposium (International) on Combustion, 19th, Technion Israel Institute of Technology, Haifa, Israel, August 8-13, 1982, Proceedings

Topics discussed are related to elementary reactions, reaction mechanisms and modeling, laminar flames, flame chemistry, turbulent reacting shear flows, turbulent premixed flames, turbulent combustion measurements, continuous combustors, detonation, detonation and explosion, heterogeneous detonation, propellant combustion, fire-ignition and thermal degradation, fire-flame spread and burning, fire-modeling, spray combustion, and droplet combustion. Coal combustion kinetics and mechanisms are considered along with coal combustion mechanisms and pyrolysis, coal combustion techniques, NOx in coal combustion, gaseous pollutants, soot and PAH, soot and inorganic pollutants, I.C. engine combustion, and ignition and extinction. Attention is given to intricate paths and simple steps in chemical kinetics and combustion, the formation of polycyclic aromatic hydrocarbons by combustion, turbulent flame structure and speed in spark ignition engines, and unresolved problems in SOx, NOx, and soot control in combustion.

Source record↗
NASA NTRS1986

Comparison of viscous shock layer and boundary layer reentry heating techniques for Orbiter nose cap

A comparison of two viscous shock layer methods and one boundary layer method for predicting the aerodynamic heating around the Orbiter nose cap during STS-5 entry is presented. The object of the study was to compare these methods with one another and with the measured Orbiter flight data for this trajectory. The nonequilibrium, chemically reacting viscous flow fields obtained by these methods are evaluated, and effects on heating rate of wall catalycity variation with time are presented. The effects of shock slip and combined wall/shock slip are considered at high altitudes (above 300,000 ft). Using the variable wall catalycity analysis, it is shown that heating rates can be predicted within a 5.7 percent flight data band for altitudes between 175,000 ft and 265,000 ft in this trajectory.

Ting, P. C.↗
NASA NTRS1987

Mechanisms by which heat release affects the flow field in a chemically reacting, turbulent mixing layer

The mechanisms by which heat release affects the fluid dynamics in a turbulent reacting mixing layer are studied by direct numerical simulation. In agreement with previous laboratory experiments, the heat release is observed to lower the rate at which the mixing layer grows and to reduce the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer are shown to produce changes in the flame vortex structure that act to produce more diffuse vortices than in the constant density case, resulting in lower rotation rates of fluid elements. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers are shown to result from vorticity generation by baroclinic torques. The density reductions also lower the generation rates of turbulent kinetic energy and the turbulent shear stresses, resulting in less turbulent mixing of fluid elements. Calculations of the energy in the various wave numbers show that the heat release has a stabilizing effect on the growth rates of individual modes. A linear stability analysis of a simlified model problem confirms this, showing that low density fluid in the mixing region will result in a shift of the frequency of the unstable modes to lower wave numbers (longer wavelengths). The growth rates of the unstable modes decrease, contributing to the slower growth of the mixing layer.

Riley, J. J.↗
NASA NTRS1991

New approach in direct-simulation of gas mixtures

Results are reported for an investigation of a new direct-simulation Monte Carlo method by which energy transfer and chemical reactions are calculated. The new method, which reduces to the variable cross-section hard sphere model as a special case, allows different viscosity-temperature exponents for each species in a gas mixture when combined with a modified Larsen-Borgnakke phenomenological model. This removes the most serious limitation of the usefulness of the model for engineering simulations. The necessary kinetic theory for the application of the new method to mixtures of monatomic or polyatomic gases is presented, including gas mixtures involving chemical reactions. Calculations are made for the relaxation of a diatomic gas mixture, a plane shock wave in a gas mixture, and a chemically reacting gas flow along the stagnation streamline in front of a hypersonic vehicle. Calculated results show that the introduction of different molecular interactions for each species in a gas mixture produces significant differences in comparison with a common molecular interaction for all species in the mixture. This effect should not be neglected for accurate DSMC simulations in an engineering context.

Chung, Chan-Hong↗
NASA NTRS1991

A new technique for low-to-high altitude predictions of ablative hypersonic flowfields

A parabolized Navier-Stokes (PNS) technique for predicting chemically reacting reentry flows over a wide range of altitudes extending from zero to 259 kft is presented. This coupled-chemistry PNS scheme can use either a strongly implicit solution approach or a two-step solution approach and is used to study a hypersonic flow over various configurations over a Mach-number range of 6.89 to 20. It is demonstrated that the two-step solution provides a combination of accuracy, stability, and computational efficiency. Nonablating low-altitude calculations show that the flowfield is in chemical equilibrium only in the near-wall region, while finite-rate chemistry effects are significant over a major portion of the shock layer, and the Blottner-type linearization of the production terms is unstable for low-altitude high-Reynolds-number conditions.

Bhutta, Bilal A.↗
NASA NTRS1992

Computation of losses in a scramjet combustor

The losses in a conceptual scramjet combustor at flight Mach numbers of 8, 10, 12, 16 and 20 are computed. These losses are extracted from three-dimensional parabolized Navier-Stokes solutions of the turbulent, reacting combustor flow field. A combustor performance index was defined based on the rationale that an efficient scramjet combustor should add heat to the fluid in such a manner as to maximize the stream thrust at the combustor exit while minimizing the losses. This index showed a decrease of more than 40 percent as the flight Mach number increased from 8 to 20, indicative of a drop in the thrust-producing potential of the scramjet at the upper end of the speed regime studied. A breakdown of the losses showed that dissipation, nonequilibrium chemistry and heat diffusion contributed roughly 15 percent, 35 percent, and 50 percent to the irreversible increase in entropy at Mach 8 and 22 percent, 13 and 65 percent at Mach 20.

Kamath, Pradeep S.↗
NASA NTRS1992

Lunar oxygen - The reduction of glass by hydrogen

The direct reduction of volcanic glass by hydrogen has been proposed as a method of extracting oxygen from the lunar soil. Experiments using lunar simulant glasses reacted with flowing hydrogen gas have demonstrated reduction at temperatures from 1000 to 1200 C. For melted samples ferrous iron was reduced to the metal, which formed large crystals at the expense of the glass. Samples held below the melting point rapidly devitrified, and iron was formed from submicrometer crystals of ilmenite and pyroxene. Weight losses of 3.6 - 4.5 percent, depending on glass composition, were achieved in 3 hours at 1100 C. A lunar oxygen plant operating at this efficiency and utilizing Apollo 17 orange glass as a feedstock could produce 50 kg of oxygen per ton of soil. The processes of reduction and sintering of lunar soil are synergistic, and could be combined to produce both oxygen and construction material at a moon base.

Allen, Carlton C.↗
NASA NTRS1993

Modeling of turbulent supersonic H2-air combustion with a multivariate beta PDF

Recent calculations of turbulent supersonic reacting shear flows using an assumed multivariate beta PDF (probability density function) resulted in reduced production rates and a delay in the onset of combustion. This result is not consistent with available measurements. The present research explores two possible reasons for this behavior: use of PDF's that do not yield Favre averaged quantities, and the gradient diffusion assumption. A new multivariate beta PDF involving species densities is introduced which makes it possible to compute Favre averaged mass fractions. However, using this PDF did not improve comparisons with experiment. A countergradient diffusion model is then introduced. Preliminary calculations suggest this to be the cause of the discrepancy.

Baurle, R. A.↗
NASA NTRS1994

Numerical simulation of shock-induced combustion past blunt bodies using shock-fitting technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier-Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model than before.

Ahuja, J. K.↗