Study of random wave phenomena in hypergolic propellant combustion
Analytical model for describing combustion chamber disturbances during start transient and steady state operations of hypergolic liquid rocket engines
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Analytical model for describing combustion chamber disturbances during start transient and steady state operations of hypergolic liquid rocket engines
A transient combustion model of nitramine propellants is combined with an isentropic compression shock formation model to determine the role of nitramine propellant combustion in DDT, excluding effects associated with propellant structural properties or mechanical behavior. The model is derived to represent the closed pipe experiment that is widely used to characterize explosives, except that the combustible material is a monolithic charge rather than compressed powder. Computations reveal that the transient combustion process cannot by itself produce DDT by this model. Compressibility of the solid at high pressure is the key factor limiting pressure buildups created by the combustion. On the other hand, combustion mechanisms which promote pressure buildups are identified and related to propellant formulation variables. Additional combustion instability data for nitramine propellants are presented. Although measured combustion response continues to be low, more data are required to distinguish HMX and active binder component contributions. A design for a closed vessel apparatus for experimental studies of high pressure combustion is discussed.
Several efficient pseudo time techniques have been developed for calculating steady state chemically reacting flows. The techniques include the implicit treatment of the chemical source term, point implicit multiple grid accelerator and a constant CFL condition. It turns out that these methods can be viewed as ways of rescaling the equations in time such that all chemical and convective phenomena evolve at comparable pseudo time scales. Consequently the number of iterations needed to solve reacting problems is approximately the same as for non-reacting problems. The techniques are demonstrated for a simple dissociation model and a nontrivial H2 - Air combustion model.
An a priori investigation of the applicability of flamelet-based combustion models to dual-mode scramjet combustion was performed utilizing Reynolds-averaged simulations (RAS). For this purpose, the HIFiRE Direct Connect Rig (HDCR) flowpath, fueled with a JP-7 fuel surrogate and operating in dual- and scram-mode was considered. The chemistry of the JP-7 fuel surrogate was modeled using a 22 species, 18-step chemical reaction mechanism. Simulation results were compared to experimentally-obtained, time-averaged, wall pressure measurements to validate the RAS solutions. The analysis of the dual-mode operation of this flowpath showed regions of predominately non-premixed, high-Damkohler number, combustion. Regions of premixed combustion were also present but associated with only a small fraction of the total heat-release in the flow. This is in contrast to the scram-mode operation, where a comparable amount of heat is released from non-premixed and premixed combustion modes. Representative flamelet boundary conditions were estimated by analyzing probability density functions for temperature and pressure for pure fuel and oxidizer conditions. The results of the present study reveal the potential for a flamelet model to accurately model the combustion processes in the HDCR and likely other high-speed flowpaths of engineering interest.
A flame sheet model for heat release is incorporated into a multi-dimensional fluid mechanical simulation for gas turbine application. The model assumes that the chemical reaction takes place in thin sheets compared to the length scale of mixing, which is valid for the primary combustion zone in a gas turbine combustor. In this paper, the details of the model are described and computational results are discussed.
A flexible and computationally economical model of the internal combustion engine was developed for use on large digital computer systems. It is based on a system of ordinary differential equations for cylinder-averaged properties. The computer program is capable of multicycle calculations, with some parameters varying from cycle to cycle, and has restart capabilities. It can accommodate a broad spectrum of reactants, permits changes in physical properties, and offers a wide selection of alternative modeling functions without any reprogramming. It readily adapts to the amount of information available in a particular case because the model is in fact a hierarchy of five models. The models range from a simple model requiring only thermodynamic properties to a complex model demanding full combustion kinetics, transport properties, and poppet valve flow characteristics. Among its many features the model includes heat transfer, valve timing, supercharging, motoring, finite burning rates, cycle-to-cycle variations in air-fuel ratio, humid air, residual and recirculated exhaust gas, and full combustion kinetics.
A steady-state lean combustion model for monolithic catalytic combustors is given. The model, consisting of several semi-global chemical reaction steps in the gas-phase and on the surface, is capable of analyzing CO and THC emissions. In the model computation presented, the influence of operating and design parameters on the minimum combustor length is studied. Special attention is given to the effect of after-bed gas-phase reaction space. Comparison with experimental data indicates good agreement in the range of parameters covered.
A finite element algorithm for the computation of confined, axisymmetric, turbulent diffusion flames is developed. The mean mixture properties were obtained by three methods based on diffusion flame concept: without using a probability density function (PDF), with a double-delta PDF, and with a beta PDF. A comparison is made for the combustion models, and the effect of turbulence on combustion are discussed.
The one point probability density function (PDF) method is examined in light of its use in actual engineering problems. The PDF method, although relatively complicated, appears to be the only format available to handle the nonlinear stochastic difficulties caused by typical reaction kinetics. Turbulence modeling, if it is to play a central role in combustion modeling, has to be integrated with the chemistry in a way which produces accurate numerical solutions to combustion problems. It is questionable whether the development of turbulent models in isolation from the peculiar statistics of reactant concentrations is a fruitful line of development as far as propulsion is concerned. There are three issues for which additional viewgraphs are prepared: the one point pdf method; the amplitude mapping closure; and a hybrid strategy for replacing a full two point pdf treatment of reacting flows by a single point pdf and correlation functions. An appeal is made for the establishment of an adequate data base for compressible flow with reactions for Mach numbers of unity or higher.
The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.
This viewgraph presentation discusses key features of current combustion system CFD modeling capabilities at GE Aircraft Engines provided by the CONCERT code; CONCERT development history; modeling applied for designing engine combustion systems; modeling applied to improve fundamental understanding; CONCERT3D results for current production combustors; CONCERT3D model of NASA/GE E3 combustor; HYBRID CONCERT CFD/Monte-Carlo modeling approach; and future modeling directions.
The phenomenon of gas dissolution in burning metals was observed in recent metal combustion studies, but it could not be adequately explained by the traditional metal combustion models. The research reported here addresses heterogeneous metal combustion with emphasis on the processes of oxygen penetration inside burning metal and its influence on the metal combustion rate, temperature history, and disruptive burning. The unique feature of this work is the combination of the microgravity environment with a novel micro-arc generator of monodispersed metal droplets, ensuring repeatable formation and ignition of uniform metal droplets with a controllable initial temperature and velocity. Burning droplet temperature is measured in real time with a three wavelength pyrometer. In addition, particles are rapidly quenched at different combustion times, cross-sectioned, and examined using SEM-based techniques to retrieve the internal composition history of burning metal particles. When the initial velocity of a spherical particle is nearly zero, the microgravity environment makes it possible to study the flame structure, the development of flame nonsymmetry, and correlation of the flame shape with the heterogeneous combustion processes.
The analyses, designs, fabrication, and cold-flow acceptance testing of LOX/RP-1 preburner components required for a high-pressure staged-combustion rocket engine are discussed. Separate designs of injectors, combustion chambers, turbine simulators, and hot-gas mixing devices are provided for fuel-rich and oxidizer-rich operation. The fuel-rich design addresses the problem of non-equilibrium LOX/RP-1 combustion. The development and use of a pseudo-kinetic combustion model for predicting operating efficiency, physical properties of the combustion products, and the potential for generating solid carbon is presented. The oxygen-rich design addresses the design criteria for the prevention of metal ignition. This is accomplished by the selection of materials and the generation of well-mixed gases. The combining of unique propellant injector element designs with secondary mixing devices is predicted to be the best approach.
This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.
Single-drop and droplet array studies have become common methods to isolate and investigate the effects of any of the complexities that enter into the drop combustion process. Microgravity environments are required to allow larger drops to be studied while minimizing or eliminating the confounding effects of buoyancy. Based on the results from current isolated drop, drop array, and spray studies funded through the Microgravity Science and Applications Division, it has become clear that even with the effects of buoyancy removed, the extrapolation of results from droplet array studies to spray flames is difficult. The problem occurs because even the simplest spray systems introduce complexities of multi-disperse drop sizes and drop-drop interactions, coupled with more complicated fluid dynamics. Not only do these features make the interpretation of experimental data difficult, they also make the problem very difficult to analyze computationally. Group combustion models, in which the interaction between droplets is treated on a statistical manner, have become a popular method to investigate the behavior of large numbers of interacting droplets, particularly through the work of Ryan et al. and Bellan and co-workers. While these models idealize the actual spray systems to a point where they can be treated computationally, the experimental analogy to these models is difficult to achieve because it requires the formation and Combustion of drop clusters without the effects of buoyancy. Therefore, even though these models have provided useful and insightful information, the verification of the results by direct comparison with experimental data is still lacking.
The performance of a scalar PDF hydrogen-air combustion model in predicting a complex reacting flow is evaluated. In addition the results are compared to those obtained by running the same case with the so-called laminar chemistry model and also a new model based on the concept of mapping partially stirred reactor data onto perfectly stirred reactor data. The results show that the scalar PDF model produces significantly different results from the other two models, and at a significantly higher computational cost.
A predictive tool for hypergolic bipropellant spray combustion and flow evolution in small rocket combustion chambers is described. It encompasses a computational technique for the gas-phase governing equations, a discrete particle method for liquid bipropellant sprays, and constitutive models for combustion chemistry, interphase exchanges, and unlike impinging hypergolic spray interactions. Emphasis is placed on the phenomenological modeling of the hypergolic liquid bipropellant gasification processes. Sample computations with the N2H4-N2O4 propellant system are given in order to show some of the capabilities and inadequacies of this tool.
A transient ballistics and combustion model is derived to represent the closed vessel experiment that is widely used to characterize the high pressure burning rates of solid propellants. The model is applied to explain why burning rates deduced from the closed vessel are in basic agreement with those measured from an equilibrium strand burner in the case of homogeneous propellants, but differ significantly in the case of nitramine composite propellants. Thermal profile time lag effects become small at high pressure because the burning rates become high. However, the development of the burning surface structure of those nitramine propellants which exhibit shifts in pressure exponent causes the mass burning rate to lag and then exceed the equilibrium value. It is necessary to consider this mechanism in applications dealing with high pressures and pressurization rates.