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.
Computer program contains mathematical model which provides relationship between engine gas dynamics and combustion processes. Mathematically simulated explosions initiate gas disturbances. Design methods for damping disturbances can be studied to prevent future engine shutdown or destruction.
Deep learning models are being widely used in the field of combustion. Given the black-box nature of typical neural network based models, uncertainty quantification (UQ) is critical to ensure the reliability of predictions as well as the training datasets, and for a principled quantification of noise and its various sources. Deep learning surrogate models for predicting properties of chemical compounds and mixtures have been recently shown to be promising for enabling data-driven fuel design and optimization, with the ultimate goal of improving efficiency and lowering emissions from combustion engines. In this study, UQ is performed for a multi-task deep learning model that simultaneously predicts the research octane number (RON), Motor Octane Number (MON), and Yield Sooting Index (YSI) of pure components and multicomponent blends. The deep learning model is comprised of three smaller networks: Extractor 1, Extractor 2, and Predictor, and a mixing operator. The molecular fingerprints of individual components are encoded via Extractor 1 and Extractor 2, the mixing operator generates fingerprints for mixtures/blends based on linear mixing operation, and the predictor maps the fingerprint to the target properties. Two different classes of UQ methods, Monte Carlo ensemble methods and Bayesian neural networks (BNNs), are employed for quantifying the epistemic uncertainty. Combinations of Bernoulli and Gaussian distributions with DropConnect and DropOut techniques are explored as ensemble methods. All the DropConnect, DropOut and Bayesian layers are applied to the predictor network. Aleatoric uncertainty is modeled by assuming that each data point has an independent uncertainty associated with it. The results of the UQ study are further analyzed to compare the performance of BNN and ensemble methods. Although this study is confined to UQ of fuel property prediction, the methodologies are applicable to other deep learning frameworks that are being widely used in the combustion community.
The characteristics of propellant mixing near the injector have a profound effect on the liquid rocket engine performance. However, the flow features near the injector of liquid rocket engines are extremely complicated, for example supercritical-pressure spray, turbulent mixing, and chemical reactions are present. Previously, a homogeneous spray approach with a real-fluid property model was developed to account for the compressibility and evaporation effects such that thermodynamics properties of a mixture at a wide range of pressures and temperatures can be properly calculated, including liquid-phase, gas- phase, two-phase, and dense fluid regions. The developed homogeneous spray model demonstrated a good success in simulating uni- element shear coaxial injector spray combustion flows. However, the real-fluid model suffered a computational deficiency when applied to a pressure-based computational fluid dynamics (CFD) code. The deficiency is caused by the pressure and enthalpy being the independent variables in the solution procedure of a pressure-based code, whereas the real-fluid model utilizes density and temperature as independent variables. The objective of the present research work is to improve the computational efficiency of the real-fluid property model in computing thermal properties. The proposed approach is called an efficient real-fluid model, and the improvement of computational efficiency is achieved by using a combination of a liquid species and a gaseous species to represent a real-fluid species.
The paper describes the results of an experimental study of the quantitative determination of the capabilities of the combustion processes associated with coaxial gaseous propellant rocket injectors to drive combustor pressure oscillations. The data, obtained by employing the modified impedance tube technique with compressed air as the oxidizer and acetylene gas as the fuel, describe the frequency dependence of the admittance of the combined injector-combustion process. The measured data are compared with the predictions of the Feiler and Heidmann analytical model utilizing different values for the characteristic combustion time tau sub b. The values of tau sub b which result in a best fit between the measured and predicted data are indicated for different equivalence ratios. It is shown that for the coaxial injector investigated in this study the tau sub b varies between 0.7 and 1.2 msec for equivalence ratios in the range of 0.57 to 1.31. In addition, the experimental data indicate that the tested injector system could drive combustion instabilities over a frequency range that is in qualitative agreement with the predictions of the Feiler and Heidmann model.
We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.
An effort was undertaken to analyze the performance of a model Lean-Direct Injection (LDI) combustor designed to meet emissions and performance goals for NASA's N+3 program. Computational predictions of Emissions Index (EINOx) and combustor exit temperature were obtained for operation at typical power conditions expected of a small-core, high pressure-ratio (greater than 50), high T3 inlet temperature (greater than 950K) N+3 combustor. Reacting-flow computations were performed with the National Combustion Code (NCC) for a model N+3 LDI combustor, which consisted of a nine-element LDI flame-tube derived from a previous generation (N+2) thirteen-element LDI design. A consistent approach to mesh-optimization, spray-modeling and kinetics-modeling was used, in order to leverage the lessons learned from previous N+2 flame-tube analysis with the NCC. The NCC predictions for the current, non-optimized N+3 combustor operating indicated a 74% increase in NOx emissions as compared to that of the emissions-optimized, parent N+2 LDI combustor.
An effort was undertaken to analyze the performance of a model Lean-Direct Injection (LDI) combustor designed to meet emissions and performance goals for NASA's N+3 program. Computational predictions of Emissions Index (EINOx) and combustor exit temperature were obtained for operation at typical power conditions expected of a small-core, high pressure-ratio (greater than 50), high T3 inlet temperature (greater than 950K) N+3 combustor. Reacting-flow computations were performed with the National Combustion Code (NCC) for a model N+3 LDI combustor, which consisted of a nine-element LDI flame-tube derived from a previous generation (N+2) thirteen-element LDI design. A consistent approach to mesh-optimization, spraymodeling and kinetics-modeling was used, in order to leverage the lessons learned from previous N+2 flame-tube analysis with the NCC. The NCC predictions for the current, non-optimized N+3 combustor operating indicated a 74% increase in NOx emissions as compared to that of the emissions-optimized, parent N+2 LDI combustor.
Ethanol is an attractive oxygenate increasingly used for blending with petroleum-derived gasoline yielding beneficial combustion and emissions behavior for a range of internal combustion engine schemes, including stoichiometric spark-ignition and low temperature combustion (LTC). As such, it is important to fundamentally understand the autoignition behavior of gasoline/ethanol blends. This work utilizes a rapid compression machine (RCM) and a homogeneous charge compression ignition (HCCI) engine to experimentally quantify changes in fuel reactivity, through ignition delay times and preliminary heat release, for blends of 0 to 30% vol./vol. into a full boiling range research gasoline (FACE-F). Diluted/stoichiometric and undiluted/fuel-lean conditions are explored covering a wide range of compressed temperatures and pressures relevant to conventional and advanced, gasoline combustion engines. Detailed chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model in conjunction with a five-component surrogate to model the RCM experiments and provide chemical insight into the perturbative effects of ethanol on the autoignition process. The diluted/stoichiometric RCM measurements reveal that within the low-temperature regime ethanol retards first-stage and main ignition delay times, and suppresses both the rates and extents of low-temperature heat release (LTHR), while within the intermediate-temperature regime ethanol only causes slight changes. Good agreement of ignition delay time and preliminary heat release prediction is found between model and experimental results. Sensitivity and flux analyses further show that ethanol blending effects are dominated by the competition between the H-atom abstraction from ethanol and other fuel components by OH radical at low temperatures and by HO 2 radical at intermediate temperatures. These findings are consistent across both fuel loading conditions explored in this study. In addition, when HCCI engine experiments are mapped onto undiluted/lean RCM measurements under a constant combustion phasing scenario, good correspondence between the two apparatuses is observed for LTHR and start of high-temperature heat release. Finally, the current study highlights the importance of characterizing LTHR in predicting fuel behaviors in high-boost/low-temperature engines, and demonstrates that RCM experiments can provide an alternative, and more-efficient avenue for such characterization.
Pulsed combustion is receiving renewed interest as a potential route to higher performance in air breathing propulsion systems. Pulsejets offer a simple experimental device with which to study unsteady combustion phenomena and validate simulations. Previous computational fluid dynamic (CFD) simulation work focused primarily on the pulsejet combustion and exhaust processes. This paper describes a new inlet sub-model which simulates the fluidic and mechanical operation of a valved pulsejet head. The governing equations for this sub-model are described. Sub-model validation is provided through comparisons of simulated and experimentally measured reed valve motion, and time averaged inlet mass flow rate. The updated pulsejet simulation, with the inlet sub-model implemented, is validated through comparison with experimentally measured combustion chamber pressure, inlet mass flow rate, operational frequency, and thrust. Additionally, the simulated pulsejet exhaust flowfield, which is dominated by a starting vortex ring, is compared with particle imaging velocimetry (PIV) measurements on the bases of velocity, vorticity, and vortex location. The results show good agreement between simulated and experimental data. The inlet sub-model is shown to be critical for the successful modeling of pulsejet operation. This sub-model correctly predicts both the inlet mass flow rate and its phase relationship with the combustion chamber pressure. As a result, the predicted pulsejet thrust agrees very well with experimental data.
Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Pressure strongly impacts thermal radiation, leading to a substantial difference in heat transfer between POC and conventional atmospheric pressure combustion. A good understanding of thermal radiation in POC is needed to aid new combustor development. The present work is a step in this direction, initiating a systematic analysis of thermal radiation and heat transfer in a pilot-scale POC combustor, which has been developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Few researchers have investigated the thermal radiation of pressurized oxy-combustion. This work adopted the discrete ordinate model to simulate the radiation of a 1-D cylinder with predetermined conditions (temperature, pressure, gas compositions) to predict the influence of radiative properties of flue gas on the radiation in the combustor.
A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.
A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.
A multiphase computational fluid dynamics code (ARICC-3D) is presented and results of two simulations are discussed. The numerical framework of the CFD code is reviewed as well as some of the two-phase physical submodels. The simulations performed include a single coaxial element injector and a multielement injector using LOX/Hydrogen reactants. The single element injector simulation verified the interaction among the code's submodels. The multielement injector simulation transient results include the chamber response to a transverse pressure wave with and without a chamber baffle. The results of these simulations demonstrate the current capabilities and their limitations to model complex two-phase combustion phenomena. Possible ways to exceed these limitations are suggested.
The main objective of this research was to develop an efficient three-dimensional computer code for chemically reacting flows. The main computer code developed is ALLSPD-3D. The ALLSPD-3D computer program is developed for the calculation of three-dimensional, chemically reacting flows with sprays. The ALL-SPD code employs a coupled, strongly implicit solution procedure for turbulent spray combustion flows. A stochastic droplet model and an efficient method for treatment of the spray source terms in the gas-phase equations are used to calculate the evaporating liquid sprays. The chemistry treatment in the code is general enough that an arbitrary number of reaction and species can be defined by the users. Also, it is written in generalized curvilinear coordinates with both multi-block and flexible internal blockage capabilities to handle complex geometries. In addition, for general industrial combustion applications, the code provides both dilution and transpiration cooling capabilities. The ALLSPD algorithm, which employs the preconditioning and eigenvalue rescaling techniques, is capable of providing efficient solution for flows with a wide range of Mach numbers. Although written for three-dimensional flows in general, the code can be used for two-dimensional and axisymmetric flow computations as well. The code is written in such a way that it can be run in various computer platforms (supercomputers, workstations and parallel processors) and the GUI (Graphical User Interface) should provide a user-friendly tool in setting up and running the code.
Two fuel-rich and two oxidizer-rich preburner injectors are tested with LOX/RP-1 in an investigation of performance, stability and gas temperature uniformity over a chamber pressure range from 1292 to 2540 psia. Fuel-rich mixture ratios range from 0.238 to 0.367 and oxidizer-rich mixture ratios range from 27 to 48, and carbon deposition data are collected by measuring the pressure drop across a turbine simulator flow device. The oxidizer-rich testing demonstrates the feasibility of oxidizer-rich preburners, indicating equilibrium combustion as predicted, and the measured fuel-rich gas composition and C-asterisk performance are in excellent agreement with kinetic model predictions indicating kinetically-limited combustion.
Internal combustion engine design is increasingly driven by computational models used to predict change in performance due to change of design. The design process previously depended on limited intuition and expensive and time-consuming physical testing. Improved model capabilities shorten design cycles and enable the production of cleaner and more efficient engines. This project focuses on advancing the state of the art in internal combustion engine simulations. The overarching goal is to enable predictive models and reduced time to solution for simulations that impact combustion engine design.
Coming into the Combustion Branch of the Turbomachinery and Propulsion Systems Division, there was not any set project planned out for me to work on. This was understandable, considering I am only at my sophmore year in college. Also, my mentor was a division chief and it was expected that I would be passed down the line. It took about a week for me to be placed with somebody who could use me. My first project was to write a macro for TecPlot. Commonly, a person would have a 3D contour volume modeling something such as a combustion engine. This 3D volume needed to have slices extracted from it and made into 2D scientific plots with all of the appropriate axis and titles. This was very tedious to do by hand. My macro needed to automate the process. There was some education I needed before I could start, however. First, TecPlot ran on Unix and Linux, like a growing majority of scientific applications. I knew a little about Linux, but I would need to know more to use the software at hand. I took two classes at the Learning Center on Unix and am now comfortable with Linux and Unix. I already had taken Computer Science I and II, and had undergone the transformation from Computer Programmer to Procedural Epistemologist. I knew how to design efficient algorithms, I just needed to learn the macro language. After a little less than a week, I had learned the basics of the language. Like most languages, the best way to learn more of it was by using it. It was decided that it was best that I do the macro in layers, starting simple and adding features as I went. The macro started out slicing with respect to only one axis, and did not make 2D plots out of the slices. Instead, it lined them up inside the solid. Next, I allowed for more than one axis and placed each slice in a separate frame. After this, I added code that transformed each individual slice-frame into a scientific plot. I also made frames for composite volumes, which showed all of the slices in the same XYZ space. I then designed an addition companion macro that exported each frame into its own image file. I then distributed the macros to a test group, and am awaiting feedback. In the meantime, a am researching the possible applications of distributed computing on the National Combustor Code. Many of our Linux boxes were idle for most of the day. The department thinks that it would be wonderful if we could get all of these idle processors to work on a problem under the NCC code. The client software would have to be easily distributed, such as in screensaver format or as a program that only ran when the computer was not in use. This project proves to be an interesting challenge.
Analytical and computational techniques were developed to predict the stability behavior of liquid propellant rocket combustors using damping devices such as acoustic liners, slot absorbers, and injector face baffles. Models were developed to determine the frequency and decay rate of combustor oscillations, the spatial and temporal pressure waveforms, and the stability limits in terms of combustion response model parameters.