Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemically 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 343 records · Page 19

The flame structure and vorticity generated by a chemically reacting transverse jet

An analytical model describing the behavior of a turbulent fuel jet injected normally into a cross flow is developed. The model places particular emphasis on the contrarotating vortex pair associated with the jet, and predicts the flame length and shape based on entrainment of the oxidizer by the fuel jet. Effects of buoyancy and density variations in the flame are neglected in order to isolate the effects of large-scale mixing. The results are compared with a simulation of the transverse reacting jet in a liquid (acid-base) system. For a wide range of ratios of the cross flow to jet velocity, the model predicts flame length quite well. In particular, the observed transitional behavior in the flame length between cross-flow velocity to jet velocity of orifice ratios of 0.0 to 0.1, yielding an approximate minimum at the ratio 0.05, is reproduced very clearly by the present model. The transformation in flow structure that accounts for this minimum arises from the differing components of vorticity dominant in the near-field and far-field regions of the jet.

Karagozian, A. R.↗

Large eddy simulation and direct numerical simulation of high speed turbulent reacting flows

The objective of this research is to make use of Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first phase of this research conducted within the past three years have been directed in several issues pertaining to intricate physics of turbulent reacting flows. In our previous 5 semi-annual reports submitted to NASA LaRC, as well as several technical papers in archival journals, the results of our investigations have been fully described. In this progress report which is different in format as compared to our previous documents, we focus only on the issue of LES. The reason for doing so is that LES is the primary issue of interest to our Technical Monitor and that our other findings were needed to support the activities conducted under this prime issue. The outcomes of our related investigations, nevertheless, are included in the appendices accompanying this report. The relevance of the materials in these appendices are, therefore, discussed only briefly within the body of the report. Here, results are presented of a priori and a posterior analyses for validity assessments of assumed Probability Density Function (PDF) methods as potential subgrid scale (SGS) closures for LES of turbulent reacting flows. Simple non-premixed reacting systems involving an isothermal reaction of the type A + B yields Products under both chemical equilibrium and non-equilibrium conditions are considered. A priori analyses are conducted of a homogeneous box flow, and a spatially developing planar mixing layer to investigate the performance of the Pearson Family of PDF's as SGS models. A posteriori analyses are conducted of the mixing layer using a hybrid one-equation Smagorinsky/PDF SGS closure. The Smagorinsky closure augmented by the solution of the subgrid turbulent kinetic energy (TKE) equation is employed to account for hydrodynamic fluctuations, and the PDF is employed for modeling the effects of scalar fluctuations. The implementation of the model requires the knowledge of the local values of the first two SGS moments. These are provided by additional modeled transport equations. In both a priori and a posteriori analyses, the predicted results are appraised by comparison with subgrid averaged results generated by DNS. Based on these results, the paths to be followed in future investigations are identified.

Adumitroaie, V.↗

A co-kurtosis PCA based dimensionality reduction with nonlinear reconstruction using neural networks

For turbulent reacting flow systems, identification of low-dimensional representations of the thermo-chemical state space is vitally important, primarily to significantly reduce the computational cost of device-scale simulations. Principal component analysis (PCA), and its variants, are a widely employed class of methods. Recently, an alternative technique that focuses on higher-order statistical interactions, co-kurtosis PCA (CoK-PCA), has been shown to effectively provide a low-dimensional representation by capturing the stiff chemical dynamics associated with spatiotemporally localized reaction zones. While its effectiveness has only been demonstrated based on a priori analyses with linear reconstruction, in this work, we employ nonlinear techniques to reconstruct the full thermo-chemical state and evaluate the efficacy of CoK-PCA compared to PCA. Specifically, we combine a CoK-PCA-/PCA-based dimensionality reduction (encoding) with an artificial neural network (ANN) based reconstruction (decoding) and examine, a priori, the reconstruction errors of the thermo-chemical state. In addition, we evaluate the errors in species production rates and heat release rates, which are nonlinear functions of the reconstructed state, as a measure of the overall accuracy of the dimensionality reduction technique. We employ four datasets to assess CoK-PCA/PCA coupled with ANN-based reconstruction: zero-dimensional (homogeneous) reactor for autoignition of an ethylene/air mixture that has conventional single-stage ignition kinetics, a dimethyl ether (DME)/air mixture which has two-stage (low and high temperature) ignition kinetics, a one-dimensional freely propagating premixed ethylene/air laminar flame, and a two-dimensional dataset representing turbulent autoignition of ethanol in a homogeneous charge compression ignition (HCCI) engine. Finally, results from the analyses demonstrate the robustness of the CoK-PCA based low-dimensional manifold with ANN reconstruction in accurately capturing the data, specifically from the reaction zones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computer user's guide for a chemically reacting viscous shock-layer program

A description is given of the computer code for predicting viscous shock-layer flows over nonanalytic blunt bodies (Program VISLNABB) for hypersonic, low Reynolds number flows. Four specific and one general body geometries are considered. In addition to sphere-cones, cylinder wedges and geometries defined in tabular form, options for hyperboloids and paraboloids are included. Details of the theory and results are included in a separate engineering report. The program, subroutines, variables in common, and input and output data are described. Listings of the program code, output data for a sample case, and the input data for this sample case are included.

Miner, E. W.↗

The Effects of Buoyancy and Dilution on the Structure and Lift-off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of the factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. The goal of our characterizations of coflow laminar diffusion flames is to bring to microgravity the multidimensional diagnostic tools available in normal gravity, and in so doing provide a broader understanding of the successes and limitations of current combustion models. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in both buoyant and nonbuoyant environments. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2)delta, denoted CH*), and excited-state OH (A(exp 2)Sigma, denoted OH*) are measured in mu-g and normal gravity. Two-dimensional CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC-135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. Although CH* and OH* kinetics are not well understood, the CH*, OH*, and ground-state CH distributions are spatially coincident in the flame anchoring region. Therefore, the ground-state CH distribution, which is easily computed, and the readily measured CH*/OH* distributions can be used to provide a consistent and convenient way of measuring lift-off height and flame shape in the diffusion flame under investigation. Given that the fuel composition affects flame chemistry and that buoyancy influences the velocity profile of the flow, we have the opportunity to computationally and experimentally study the roles of fluids and chemistry. In performing this microgravity study, improvements to the computational model have been made and new calculations performed for a range of gravity and flow conditions. Furthermore, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented in the following sections.

Walsh, Kevin T.↗

The Effects of Buoyancy and Dilution on the Structure and Lift-Off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of those factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. A goal of this work is to bring to microgravity flame studies the detailed experimental and numerical tools that have been used to study ground-based systems. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in a nonbuoyant environment. To better understand these phenomena, experimental and computational studies of a coflow laminar diffusion flame have been carried out. To date, these studies have focused on a single set of flow conditions, in which a nitrogen-diluted methane fuel stream (65% methane by volume) was surrounded by an air coflow, with exit velocities matched at 35 cm/s. Of particular interest is the change in flame shape due to the absence of buoyant forces, as well as the amount of diluent in the fuel stream and the coflow velocity. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2 delta) denoted CH*), and excited-state OH (A(exp 2 sigma, denoted OH*) are measured. CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. In extending the study to microgravity conditions, improvements to the computational model have been made and new calculations performed for a range of gravity conditions. In addition, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented.

Walsh, Kevin T.↗

Nonequilibrium effects on the aerothermodynamics of transatmospheric and aerobraking vehicles

A 3D CFD algorithm is used to study the effect of thermal and chemical nonequilibrium on slender and blunt body aerothermodynamics. Both perfect gas and reacting gas air models are used to compute the flow over a generic transatmospheric vehicle and a proposed lunar transfer vehicle. The reacting air is characterized by a translational-rotational temperature and a vibrational-electron-electronic temperature and includes eight chemical species. The effects of chemical reaction, vibrational excitation, and ionization on lift-to-drag ratio and trim angle are investigated. Results for the NASA Ames All-body Configuration show a significant difference in center of gravity location for a reacting gas flight case when compared to a perfect gas wind tunnel case at the same Mach number, Reynolds number, and angle of attack. For the same center of gravity location, the wind tunnel model trims at lower angle of attack than the full-scale flight case. Nonionized and ionized results for a proposed lunar transfer vehicle compare well to computational results obtained from a previously validated reacting gas algorithm. Under the conditions investigated, effects of weak ionization on the heat transfer and aerodynamic coefficients were minimal.

Hassan, Basil↗

Numerical Investigation of Film Coefficient Engineering Methodology for Dissociated, Chemically Reacting Boundary Layers

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flow field for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for their application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. The sensitivity of the engineering model is assessed based on the impact of the CFD assumptions as an effort to make it available to design engineers across industry.

numerical boundary layer↗

Physical Simulation of Rocket Exhaust Aerodynamics Using Heated Ethane: Conceptual Foundations

Heated ethane (C2H6) has been proposed as an alternative to inert gases for use as a motive fluid in the experimental simulation of rocket exhaust plumes. By adjusting stagnation temperature, the isentropic exponent of ethane can be tuned to approximate those produced by common rocket propellants including hydrogen, hypergols, alcohols, and hydrocarbons. As a result, ethane can be made to follow a nozzle expansion process which is nearly identical to realistic rocket engine flow fields. Additionally, its high auto-ignition temperature and resistance to condensation enable the testing of expansion ratios much larger than conventional inertgas testing. NASA SSC has performed quasi-one-dimensional analyses using the Chemical Equilibrium with Applications (CEA) code as a preliminary means to compare flow fields produced by non-reacting ethane to those of reacting combustion products. A LO2/LH2 rocket engine operating at a chamber pressure of 5.0 MPa and a mixture ratio of 6.1 was used as an example case to demonstrate ethane’s efficacy as a simulant. Errors for key similarity parameters were compared to legacy cold-flow test methods. Additional errors induced by machining tolerances and chemical impurities were also examined. Results suggest that at a 3% geometric scale and ~500 K ethane stagnation temperature, an error of less than 2.5% throughout the flow field is realistically achievable along the dimensions of Mach number, Reynolds number, pressure ratio, and isentropic exponent. The development of an experimental test bed for validation of this configuration is currently underway.

Jones, Daniel R.↗

PeleMP: The Multiphysics Solver for the Combustion Pele Adaptive Mesh Refinement Code Suite

Combustion encompasses multiscale, multiphase reacting flow physics spanning a wide range of scales from the molecular scales, where chemical reactions occur, to the device scales, where the turbulent flow is affected by the geometry of the combustor. This scale disparity and the limited measurement capabilities from experiments make modeling combustion a significant challenge. Recent advancements in high-performance computing (HPC), particularly with the Department of Energy's Exascale Computing Project (ECP), have enabled high-fidelity simulations of practical applications to be performed. The major physics submodels, including chemical reactions, turbulence, sprays, soot, and thermal radiation, exhibit distinctive computational characteristics that need to be examined separately to ensure efficient utilization of computational resources. This paper presents the multiphysics solver for the Pele code suite, called PeleMP, which consists of models for spray, soot, and thermal radiation. Here, the mathematical and algorithmic aspects of the model implementations are described in detail as well as the verification process. The computational performance of these models is benchmarked on multiple supercomputers, including Frontier, an exascale machine. Results are presented from production simulations of a turbulent sooting ethylene flame and a bluff-body swirl stabilized spray flame with sustainable aviation fuels to demonstrate the capability of the Pele codes for modeling practical combustion problems with multiphysics. This work is an important step toward the exascale computing era for high-fidelity combustion simulations providing physical insights and data for predictive modeling of real-world devices.

42 ENGINEERING↗

The impact of time step definition on code convergence and robustness

We have implemented preconditioning for multi-species reacting flows in two independent codes, an implicit (ADI) code developed in-house and the RPLUS code (developed at LeRC). The RPLUS code was modified to work on a four-stage Runge-Kutta scheme. The performance of both the codes was tested, and it was shown that preconditioning can improve convergence by a factor of two to a hundred depending on the problem. Our efforts are currently focused on evaluating the effect of chemical sources and on assessing how preconditioning may be applied to improve convergence and robustness in the calculation of reacting flows.

Venkateswaran, S.↗

Computational reacting gas dynamics

In the study of high speed flows at high altitudes, such as that encountered by re-entry spacecrafts, the interaction of chemical reactions and other non-equilibrium processes in the flow field with the gas dynamics is crucial. Generally speaking, problems of this level of complexity must resort to numerical methods for solutions, using sophisticated computational fluid dynamics (CFD) codes. The difficulties introduced by reacting gas dynamics can be classified into three distinct headings: (1) the usually inadequate knowledge of the reaction rate coefficients in the non-equilibrium reaction system; (2) the vastly larger number of unknowns involved in the computation and the expected stiffness of the equations; and (3) the interpretation of the detailed reacting CFD numerical results. The research performed accepts the premise that reacting flows of practical interest in the future will in general be too complex or 'untractable' for traditional analytical developments. The power of modern computers must be exploited. However, instead of focusing solely on the construction of numerical solutions of full-model equations, attention is also directed to the 'derivation' of the simplified model from the given full-model. In other words, the present research aims to utilize computations to do tasks which have traditionally been done by skilled theoreticians: to reduce an originally complex full-model system into an approximate but otherwise equivalent simplified model system. The tacit assumption is that once the appropriate simplified model is derived, the interpretation of the detailed numerical reacting CFD numerical results will become much easier. The approach of the research is called computational singular perturbation (CSP).

Lam, S. H.↗

Time accurate simulations of compressible shear flows

The objectives of this research are to employ direct numerical simulation (DNS) to study the phenomenon of mixing (or lack thereof) in compressible free shear flows and to suggest new means of enhancing mixing in such flows. The shear flow configurations under investigation are those of parallel mixing layers and planar jets under both non-reacting and reacting nonpremixed conditions. During the three-years of this research program, several important issues regarding mixing and chemical reactions in compressible shear flows were investigated.

Givi, Peyman↗

Numerical solutions of the complete Navier-Stokes equations

The physical phenomena within supersonic flows that sustain chemical reactions are investigated. An earlier study to develop accurate physical models for supersonic reacting flowfields focused on 2-D laminar shear layers. The objective is to examine the mixing and subsequent combustion within turbulent reacting shear layers. To conduct this study, a computer program has been written to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method uses a cell-centered finite volume approach and a Runge Kutta time stepping scheme. The Reynolds averaged equations are closed using the eddy viscosity concept. Several zero-equation models have been tested by making calculations for an H2-air nonreacting coaxial jet flow. Comparisons made with experimental data show that Cohen's eddy viscosity model provides best agreement. The finite rate chemistry model used in the study of 2-D laminar shear layers is incorporated into the computer program and data is compared from a recent experiment performed at NASA Langley.

Hassan, H. A.↗

Nitric oxide formation by meteoroids in the upper atmosphere

The process of nitric oxide formation during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source in a uniform flow with a continuum regime evolving in its wake. The amount of nitric oxide produced by high-temperature reactions of air in the continuum regime is calculated by numerical integration of chemical-rate equations. This is accomplished by assuming that flow properties are constant across the reacting region, its radius being determined from considerations of shock-wave formation and molecular diffusion. The results, when summed over the observed mass, velocity, and entry-angle distributions of meteoroids, provide annual global production rates of nitric oxide as a function of altitude. The peak production is found to occur between 90- and 100-km altitude, the total annual rate being around 40 million kg.

Menees, G. P.↗

Calculations of condensation and chemistry in an aircraft contrail

The flow field, chemistry, and condensation nucleation behind a transport airplane are calculated in two regimes using two separate reacting flow codes: first the axisymmetric plume, then the three dimensional vortex wake. The included chemical kinetics equations follow the evolution of the NO(y) and SO(x) chemical families. In the plume regime, the chemistry is coupled with the binary homogeneous formation of sulfate condensation nuclei, where the calculated nucleation rates predict that copious quantities of H2SO4/H2O nuclei are produced in subnanometer sizes. These sulfate aerosols could play a major role in the subsequent condensation of water vapor and the formation of contrails under favorable atmospheric conditions.

Miake-Lye, Richard C.↗

Nitric Oxide Formation by Meteoroids in the Upper Atmosphere

The process of nitric oxide formation during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source in a uniform flow with a continuum regime evolving in its wake. The amount of nitric oxide produced by high-temperature reactions of air in the continuum regime is calculated by numerical integration of chemical-rate equations. This is accomplished by assuming that flow properties are constant across the reacting region, the radius of the region being determined from considerations of shock-wave formation and molecular diffusion. The results, when summed over the observed mass, velocity, and entry-angle distributions of meteoroids, provide annual global production rates of nitric oxide as a function of altitude. The peak production of nitric oxide is found to occur at altitudes between 9 x 10(exp 4) and 10(exp 5) m, the total annual rate being about 4 x 10(exp 7) kg. The present results suggest that the large concentration of nitric oxide observed below 9.5 x 10(exp 4) m could be attributed to meteoroids instead of photodissociation of nitrogen into metastable, 2D-state atoms, as has been previously hypothesized.

Menees, Gene P.↗

Sphere based fluid systems

Systems are described for using multiple closely-packed spheres. In one system for passing fluid, a multiplicity of spheres lie within a container, with all of the spheres having the same outside diameter and with the spheres being closely nested in one another to create multiple interstitial passages of a known size and configuration and smooth walls. The container has an inlet and outlet for passing fluid through the interstitial passages formed between the nested spheres. The small interstitial passages can be used to filter out material, especially biological material such as cells in a fluid, where the cells can be easily destroyed if passed across sharp edges. The outer surface of the spheres can contain a material that absorbs a constitutent in the flowing fluid, such as a particular contamination gas, or can contain a catalyst to chemically react the fluid passing therethrough, the use of multiple small spheres assuring a large area of contact of these surfaces of the spheres with the fluid. In a system for storing and releasing a fluid such as hydrogen as a fuel, the spheres can include a hollow shell containing the fluid to be stored, and located within a compressable container that can be compressed to break the shells and release the stored fluid.

Elleman, Daniel D.↗