Pareto-efficient combustion modeling for improved CO-emission prediction in LES of a piloted turbulent dimethyl ether jet flame
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
Model for prediction of rocket engine combustion with performance parameters of gas phase turbulent diffusion, and mixing of oxidizer and fuel vapors
An interactive computer code for simulation of a high-intensity turbulent combustor as a single point inhomogeneous stirred reactor was developed from an existing batch processing computer code CDPSR. The interactive CDPSR code was used as a guide for interpretation and direction of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser diagnostic techniques to experimentally determine the appropriate mixing frequency, and for validation of CDPSR as a mixing-chemistry model for a laboratory jet-stirred reactor. The coalescence-dispersion model for finite rate mixing was incorporated into an existing interactive code AVCO-MARK I, to enable simulation of a combustor as a modular array of stirred flow and plug flow elements, each having a prescribed finite mixing frequency, or axial distribution of mixing frequency, as appropriate. Further increase the speed and reliability of the batch kinetics integrator code CREKID was increased by rewriting in vectorized form for execution on a vector or parallel processor, and by incorporating numerical techniques which enhance execution speed by permitting specification of a very low accuracy tolerance.
Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.
Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.
Combustion instability prediction using nonlinear bipropellant vaporization model
In this project a transient, fully three-dimensional computer simulation code was developed to simulate the effects of realistic g-jitter on a number of combustion systems. The simulation code is capable of simulating flame spread on a solid and nonpremixed or premixed gaseous combustion in nonturbulent flow with simple combustion models. Simple combustion models were used to preserve computational efficiency since this is meant to be an engineering code. Also, the use of sophisticated turbulence models was not pursued (a simple Smagorinsky type model can be implemented if deemed appropriate) because if flow velocities are large enough for turbulence to develop in a reduced gravity combustion scenario it is unlikely that g-jitter disturbances (in NASA's reduced gravity facilities) will play an important role in the flame dynamics. Acceleration disturbances of realistic orientation, magnitude, and time dependence can be easily included in the simulation. The simulation algorithm was based on techniques used in an existing large eddy simulation code which has successfully simulated fire dynamics in complex domains. A series of simulations with measured and predicted acceleration disturbances on the International Space Station (ISS) are presented. The results of this series of simulations suggested a passive isolation system and appropriate scheduling of crew activity would provide a sufficiently "quiet" acceleration environment for spherical diffusion flames.
Major research advances in heat transfer and fluid dynamics are outlined, with particular reference to relevant energy problems. Of significant importance are such topics as synthetic fuels in combustion, turbulence models, combustion modeling, numerical methods for interacting boundary layers, and light-scattering diagnostics for gases. The discussion covers thermal convection, two-phase flow and boiling heat transfer, turbulent flows, combustion, and aerospace heat transfer problems. Other areas discussed include compressible flows, fluid mechanics and drag, and heat exchangers. Featured topics comprise heat and salt transfer in double-diffusive systems, limits of boiling heat transfer in a liquid-filled enclosure, investigation of buoyancy-induced flow stratification in a cylindrical plenum, and digital algorithms for dynamic analysis of a heat exchanger. Individual items are announced in this issue.
A computerized mathematical model of the combustion response function of composite solid propellants was developed with particular attention to the contributions of the solid phase heterogeneity. The one-dimensional model treats the solid phase as alternating layers of ammonium perchlorate and binder, with an exothermic melt layer at the surface. Solution of the Fourier heat equation in the solid provides temperature and heat flux distributions with space and time. The problem is solved by conserving the heat flux at the surface from that produced by a suitable model of the gas phase. An approximation of the BDP flame model is utilized to represent the gas phase. By the use of several reasonable assumptions, it is found that a significant portion of the problem can be solved in closed form. A method is presented by which the model can be applied to tetramodal particle size distributions. A computerized steady-state version of the model was completed, which served to validate the various approximations and lay a foundation for the combustion response modeling. The combustion response modeling was completed in a form which does not require an iterative solution, and some preliminary results were acquired.
Combustion instabilities can lead to increased development time and cost for aeroengine gas turbines. This problem has been evident in the development of very-low emissions stationary gas turbines, and will likely be encountered in the newer, more aggressive aeroengine designs. In order to minimize development time and cost, it is imperative that potential combustion dynamics issues be resolved using analyses and smaller-scale experimentation. This paper discusses a methodology through which a problem in a full-scale engine was replicated in a single-nozzle laboratory combustor. Specifically, this approach is valid for longitudinal and "bulk" mode combustion instabilities. An explanation and partial validation of the acoustic analyses that were used to achieve this replication are also included. This approach yields a testbed for the diagnosis of combustion dynamics problems and for their solution through passive and active control techniques.
Combustion instabilities in gas turbine engines are most frequently encountered during the late phases of engine development, at which point they are difficult and expensive to fix. The ability to replicate an engine-traceable combustion instability in a laboratory-scale experiment offers the opportunity to economically diagnose the problem more completely (to determine the root cause), and to investigate solutions to the problem, such as active control. The development and validation of active combustion instability control requires that the casual dynamic processes be reproduced in experimental test facilities which can be used as a test bed for control system evaluation. This paper discusses the process through which a laboratory-scale experiment and be designed to replicate an instability observed in a developmental engine. The scaling process used physically-based analyses to preserve the relevant geometric, acoustic, and thermo-fluid features, ensuring that results achieved in the single-nozzle experiment will be scalable to the engine.
A three-dimensional computer code (KIVA) is being modified to include state-of-the-art submodels for diesel engine flow and combustion: spray atomization, drop breakup/coalescence, multi-component fuel vaporization, spray/wall interaction, ignition and combustion, wall heat transfer, unburned HC and NOx formation, soot and radiation, and the intake flow process. Improved and/or new submodels which were completed are: wall heat transfer with unsteadiness and compressibility, laminar-turbulent characteristic time combustion with unburned HC and Zeldo'vich NOx, and spray/wall impingement with rebounding and sliding drops. Results to date show that adding the effects of unsteadiness and compressibility improves the accuracy of heat transfer predictions; spray drop rebound can occur from walls at low impingement velocities (e.g., in cold-starting); larger spray drops are formed at the nozzle due to the influence of vaporization on the atomization process; a laminar-and-turbulent characteristic time combustion model has the flexibility to match measured engine combustion data over a wide range of operating conditions; and finally, the characteristic time combustion model can also be extended to allow predictions of ignition. The accuracy of the predictions is being assessed by comparisons with available measurements. Additional supporting experiments are also described briefly. To date, comparisons with measured engine cylinder pressure and heat flux data were made for homogeneous charge, spark-ignited and compression-ignited engines. The model results are in good agreement with the experiments.
A model of spray combustion which employs a second-order turbulence model was developed. The assumption of locally homogeneous flow is made, implying infinitely fast transport rates between the phase. Measurements to test the model were completed for a gaseous n-propane flame and an air atomized n-pentane spray flame, burning in stagnant air at atmospheric pressure. Profiles of mean velocity and temperature, as well as velocity fluctuations and Reynolds stress, were measured in the flames. The predictions for the gas flame were in excellent agreement with the measurements. The predictions for the spray were qualitatively correct, but effects of finite rate interphase transport were evident, resulting in a overstimation of the rate development of the flow. Predictions of spray penetration length at high pressures, including supercritical combustion conditions, were also completed for comparison with earlier measurements. Test conditions involved a pressure atomized n-pentane spray, burning in stagnant air at pressures of 3, 5, and 9 MPa. The comparison between predictions and measurements was fair. This is not a very sensitive test of the model, however, and further high pressure experimental and theoretical results are needed before a satisfactory assessment of the locally homogeneous flow approximation can be made.
Liquid rocket engines using oxygen and methane propellants are being considered by the National Aeronautics and Space Administration (NASA) for future in-space vehicles. This propellant combination has not been previously used in flight-qualified engine systems developed by NASA, so limited test data and analysis results are available at this stage of early development. As part of activities for the Propulsion and Cryogenic Advanced Development (PCAD) project funded under the Exploration Technology Development Program, the NASA Marshall Space Flight Center (MSFC) has been evaluating capability to model combustion performance and stability for oxygen and methane propellants. This activity has been proceeding for about two years and this paper is a summary of results to date. Hot-fire test results of oxygen/methane propellant rocket engine combustion devices for the modeling investigations have come from several sources, including multi-element injector tests with gaseous methane from the 1980s, single element tests with gaseous methane funded through the Constellation University Institutes Program, and multi-element injector tests with both gaseous and liquid methane conducted at the NASA MSFC funded by PCAD. For the latter, test results of both impinging and coaxial element injectors using liquid oxygen and liquid methane propellants are included. Configurations were modeled with two one-dimensional liquid rocket combustion analysis codes, the Rocket Combustor Interactive Design and Analysis code and the Coaxial Injector Combustion Model. Special effort was focused on how these codes can be used to model combustion and performance with oxygen/methane propellants a priori, and what anchoring or calibrating features need to be applied, improved or developed in the future. Low frequency combustion instability (chug) occurred, with frequencies ranging from 150 to 250 Hz, with several multi-element injectors with liquid/liquid propellants, and was modeled using techniques from Wenzel and Szuch. High-frequency combustion instability also occurred at the first tangential (1T) mode, at about 4500 Hz, with several multi-element injectors with liquid/liquid propellants. Analyses of the transverse mode instability were conducted by evaluating injector resonances and empirical methods developed by Hewitt.
Modeling plays a vital role in providing physical insights into behavior revealed by experiment. The program at the University of Illinois is designed to improve our understanding of basic combustion phenomena through the analytical and numerical modeling of a variety of configurations undergoing experimental study in NASA's microgravity combustion program. Significant progress has been made in two areas: (1) flame-balls, studied experimentally by Ronney and his co-workers; (2) particle-cloud flames studied by Berlad and his collaborators. Additional work is mentioned below. NASA funding for the U. of Illinois program commenced in February 1991 but work was initiated prior to that date and the program can only be understood with this foundation exposed. Accordingly, we start with a brief description of some key results obtained in the pre - 2/91 work.
The status of modelling the coherent anti-Stokes Raman spectroscopy (CARS) spectra of molecules important in combustion, such as N2, H2O, and CO2, is reviewed. It is shown that accurate modelling generally requires highly precise knowledge of line positions and reasonable estimates of Raman linewidths, and the sources of these data are discussed. CARS technique and theory is reviewed, and the status of modelling the phenomenon of collisional narrowing at pressures well above atmospheric for N2, H2O, and CO2 is described. It is shown that good agreement with experiment can be achieved using either the Gordon rotational diffusion model or phenomenological models for inelastic energy transfer rates.
Combustion in turbulent flow - models
Volume 1, the first of three volumes is a compilation of 47 unclassified/unlimited-distribution technical papers presented at the Joint Army-Navy-NASA-Air Force (JANNAF) 36th Combustion Subcommittee held jointly with the 24th Airbreathing Propulsion Subcommittee and 18th Propulsion Systems Hazards Subcommittee. The meeting was held on 18-21 October 1999 at NASA Kennedy Space Center and The DoubleTree Oceanfront Hotel, Cocoa Beach, Florida. Solid phase propellant combustion topics covered in this volume include cookoff phenomena in the pre- and post-ignition phases, solid rocket motor and gun propellant combustion, aluminized composite propellant combustion, combustion modeling and combustion instability and instability measurement techniques.