SEARCH · Search NASA

Results for “Compressible 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗

NASA NTRS1992

Model free simulations of a high speed reacting mixing layer

The effects of compressibility, chemical reaction exothermicity and non-equilibrium chemical modeling in a combusting plane mixing layer were investigated by means of two-dimensional model free numerical simulations. It was shown that increased compressibility generally had a stabilizing effect, resulting in reduced mixing and chemical reaction conversion rate. The appearance of 'eddy shocklets' in the flow was observed at high convective Mach numbers. Reaction exothermicity was found to enhance mixing at the initial stages of the layer's growth, but had a stabilizing effect at later times. Calculations were performed for a constant rate chemical rate kinetics model and an Arrhenius type kinetics prototype. The Arrhenius model was found to cause a greater temperature increase due to reaction than the constant kinetics model. This had the same stabilizing effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗

NASA NTRS2011

Density Weighted FDF Equations for Simulations of Turbulent Reacting Flows

In this report, we briefly revisit the formulation of density weighted filtered density function (DW-FDF) for large eddy simulation (LES) of turbulent reacting flows, which was proposed by Jaberi et al. (Jaberi, F.A., Colucci, P.J., James, S., Givi, P. and Pope, S.B., Filtered mass density function for Large-eddy simulation of turbulent reacting flows, J. Fluid Mech., vol. 401, pp. 85-121, 1999). At first, we proceed the traditional derivation of the DW-FDF equations by using the fine grained probability density function (FG-PDF), then we explore another way of constructing the DW-FDF equations by starting directly from the compressible Navier-Stokes equations. We observe that the terms which are unclosed in the traditional DW-FDF equations are now closed in the newly constructed DW-FDF equations. This significant difference and its practical impact on the computational simulations may deserve further studies.

Shih, Tsan-Hsing↗

NASA NTRS2013

Very Large Eddy Simulations of a Jet-A Spray Reacting Flow in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar DWFDF Method

This paper presents the very large eddy simulations (VLES) of a Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar DWFDF method, in which DWFDF is defined as the density weighted time filtered fine grained probability density function. The flow field is calculated by using the time filtered compressible Navier-Stokes equations (TFNS) with nonlinear subscale turbulence models, and when the Eulerian scalar DWFDF method is invoked, the energy and species mass fractions are calculated by solving the equation of DWFDF. A nonlinear subscale model for closing the convection term of the Eulerian scalar DWFDF equation is used and will be briefly described in this paper. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar DWFDF method in both improving the simulation quality and maintaining economic computing cost are observed.

Shih, Tsan-Hsing↗

NASA NTRS1991

Progress of simulations for reacting shear layers

An attempt was made to develop a high speed, chemically reactive shear layer test rig. The purpose of the experiment was to study the mixing of oxidizer and fuel streams in reacting shear layers for various density, velocity, and Mach number. The primary goal was to understand the effects of the compressibility upon mixing and combustion in a fundamental way. Therefore, a two-dimensional shear layer is highly desirable for its simplicity to quantify the compressibility effects. The RPLUS 2D code is used to calculate the flow fields of different sections of the test rig. The emphasis was on the supersonic nozzle design, the vitiation process for the hot air stream and the overall thermodynamic conditions of the test matrix. The k-epsilon turbulence model with wall function was successfully implemented in the RPLUS code. The k and epsilon equations are solved simultaneously and the LU scheme is used to make it compatible with the flow solver.

Yu, Sheng-Tao↗

NASA NTRS1988

Effect of heat release on the spatial stability of a supersonic reacting mixing layer

A numerical study of the stability of compressible mixing layers in which a diffusion flame is embedded is described. The mean velocity profile has been approximated by a hyperbolic tangent profile and the limit of infinite activation energy taken, which reduces the diffusion flame to a flame sheet. The addition of combustion in the form of a flame sheet was found to have important, and complex, effects on the flow stability.

Jackson, T. L.↗

NASA NTRS1988

Analysis of homogeneous turbulent reacting flows

Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Leonard, A. D.↗

DOE OSTI2023

PeleLMeX: an AMR Low Mach Number Reactive Flow Simulation Code without level sub-cycling

PeleLMeX simulates chemically reacting low Mach number flows with block-structured adaptive mesh refinement (AMR). The code is built upon the AMReX library, which provides the underlying data structures and tools to manage and operate on them across massively parallel computing architectures. PeleLMeX algorithmic features are inherited from its predecessor PeleLM but key improvements allow representation of more complex physical processes. Together with its compressible flow counterpart PeleC, the thermo-chemistry library PelePhysics and the multi-physics library PeleMP, it forms the Pele suite of open-source reactive flow simulation codes.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1989

Two-dimensional analysis of two-phase reacting flow in a firing direct-injection diesel engine

The flow field, spray penetration, and combustion in two-stroke diesel engines are described. Fuel injection begins at 345 degrees after top dead center (ATDC) and n-dodecane is used as the liquid fuel. Arrhenius kinetics is used to calculate the reaction rate term in the quasi-global combustion model. When the temperature, fuel, and oxygen mass fraction are within suitable flammability limits, combustion begins spontaneously. No spark is necessary to ignite a localized high temperature region. Compression is sufficient to increase the gaseous phase temperature to a point where spontaneous chemical reactions occur. Results are described for a swirl angle of 22.5 degrees.

Nguyen, H. Lee↗

NASA NTRS1985

Implicit/explicit analysis of interactive phenomena in: supersonic, chemically-reacting, mixing and boundary layer problems

An efficient implicit/explicit procedure for solving the parabolized Navier-Stokes (PNS) equations to analyze supersonic, viscous, chemically-reacting flows is presented. The approach combines a fully-implicit algorithm for solving the parabolic mixing/boundary layer equations with an explicit algorithm for determining the hyperbolic pressure field. The unified model, SPLITP, contains provisions for treating bounding/embedded subsonic regions using a pressure-splitting procedure which suppresses elliptic (upstream influence) effects without recourse to sublayer approximations. Several two-equation turbulence models with compressibility and curvature-correction terms are incorporated into SPLITP. Applications to a variety of flow problems are presented which include: ducted supersonic mixing and combustion flowfields, free jets and plumes, and, curved boundary layers and wall jets.

Dash, S. M.↗

NASA NTRS1991

Radiative interactions in chemically reacting supersonic internal flows

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with finite-rate chemistry and radiation for hydrogen-air systems. The chemistry source terms in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed flow in a channel with a ten-degree compression ramp. Three different chemistry models are used, accounting for increasing number of reactions and participating species. Two chemistry models assure nitrogen as inert, while the third model accounts for nitrogen reactions and NO(x) formation. The tangent slab approximation is used in the radiative flux formulation. A pseudo-gray model is used to represent the absorption-emission characteristics of the participating species. Results obtained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species and that this can have a significant influence on the flowfield.

Tiwari, S. N.↗

NASA NTRS1991

Simulation of mixing in the quick quench region of a rich burn-quick quench mix-lean burn combustor

A computer program was developed to study the mixing process in the quick quench region of a rich burn-quick quench mix-lean burn combustor. The computer program developed was based on the density-weighted, ensemble-averaged conservation equations of mass, momentum (full compressible Navier-Stokes), total energy, and species, closed by a k-epsilon turbulence model with wall functions. The combustion process was modeled by a two-step global reaction mechanism, and NO(x) formation was modeled by the Zeldovich mechanism. The formulation employed in the computer program and the essence of the numerical method of solution are described. Some results obtained for nonreacting and reacting flows with different main-flow to dilution-jet momentum flux ratios are also presented.

Shih, Tom I.-P.↗

NASA NTRS1989

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.↗

NASA NTRS1993

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗

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.↗

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.↗

NASA NTRS1993

Combustion chamber analysis code

A three-dimensional, time dependent, Favre averaged, finite volume Navier-Stokes code has been developed to model compressible and incompressible flows (with and without chemical reactions) in liquid rocket engines. The code has a non-staggered formulation with generalized body-fitted-coordinates (BFC) capability. Higher order differencing methodologies such as MUSCL and Osher-Chakravarthy schemes are available. Turbulent flows can be modeled using any of the five turbulent models present in the code. A two-phase, two-liquid, Lagrangian spray model has been incorporated into the code. Chemical equilibrium and finite rate reaction models are available to model chemically reacting flows. The discrete ordinate method is used to model effects of thermal radiation. The code has been validated extensively against benchmark experimental data and has been applied to model flows in several propulsion system components of the SSME and the STME.

Przekwas, A. J.↗

NASA NTRS1989

Structural development of laminar flow control aircraft chordwise wing joint designs

For laminar flow to be achieved, any protuberances on the surface must be small enough to avoid transition to turbulent flow. However, the surface must have joints between the structural components to allow assembly or replacement of damaged parts, although large continuous surfaces can be utilized to minimize the number the number of joints. Aircraft structural joints usually have many countersunk bolts or rivets on the outer surface. To maintain no mismatch on outer surfaces, it is desirable to attach the components from the inner surface. It is also desirable for the panels to be interchangeable, without the need for shims at the joint, to avoid surface discontinuities that could cause turbulence. Fabricating components while pressing their outer surfaces against an accurate mold helps to ensure surface smoothness and continuity at joints. These items were considered in evaluating the advantages and disadvantages of the joint design concepts. After evaluating six design concepts, two of the leading candidates were fabricated and tested using many small test panels. One joint concept was also built and tested using large panels. The small and large test panel deflections for the leading candidate designs at load factors up to +1.5 g's were well within the step and waviness requirements for avoiding transition.The small panels were designed and tested for compression and tension at -65 F, at ambient conditions, and at 160 F. The small panel results for the three-rib and the sliding-joint concepts indicated that they were both acceptable. The three-rib concept, with tapered splice plates, was considered to be the most practical. A modified three-rib joint that combined the best attributes of previous candidates was designed, developed, and tested. This improved joint met all of the structural strength, surface smoothness, and waviness criteria for laminar flow control (LFC). The design eliminated all disadvantages of the initial three-rib concept except for unavoidable eccentricity, which was reduced and reacted satisfactorily by the rib supports. It should also result in a relatively simple low-cost installation, and makes it easy to replace any panels damaged in the field.

Fischler, J. E.↗