Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Accelerating Innovative Energy Solutions Using Combustion Simulations

Combustion-based transportation, electricity generation, and industrial heating in manufacturing constitute the three largest sectors of energy demand. Some of the recent technology development in these sectors are: switching to low-carbon fuels for the transportation sector, increasing energy efficiency in the power sector, and capturing carbon emissions from conventional power generators. Several teams at the National Renewable Energy Laboratory have been actively advancing research in these areas by leveraging computational modeling of combustion processes across the heavy-duty land based transportation, aviation, and power generation sectors. This article summarizes some of these efforts, demonstrating the potential of advanced computational techniques to generate technological solutions that will transform the global energy system.

08 HYDROGEN↗
DOE OSTI2020

Assessment of various tar and soot treatment methods and a priori analysis of the steady laminar flamelet model for use in coal combustion simulation

Soot formed in the fuel-rich regions of coal flames is known to have a significant impact on the radiative heat transfer characteristics of the flame. It is therefore important to account for the effect of soot radiation when simulating coal combustion systems. Several recent coal combustion studies have used acetylene as a surrogate for devolatilized tar and neglect soot formation altogether, while others have utilized empirical models for tar and soot or even chemical mechanisms which include tar and soot as species. In this work, we investigate several methods of treating coal-derived tar and soot and the impact each treatment has on predicted furnace temperature, particle mass loss, and other quantities of interest. Furthermore, we perform a priori comparisons of simulation results to calculations from a steady laminar flamelet model in order to assess the accuracy of a flamelet approach for each tar and soot treatment considered. Data used as a basis for the study is generated through One-Dimensional Turbulence (ODT) simulation of a coal combustion system in which high- fidelity models are used for gas phase chemistry and devolatilization, where char combustion is neglected.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Results of Small-scale Solid Rocket Combustion Simulator testing at Marshall Space Flight Center

The Small-scale Solid Rocket Combustion Simulator (SSRCS) program was established at the Marshall Space Flight Center (MSFC), and used a government/industry team consisting of Hercules Aerospace Corporation, Aerotherm Corporation, United Technology Chemical Systems Division, Thiokol Corporation and MSFC personnel to study the feasibility of simulating the combustion species, temperatures and flow fields of a conventional solid rocket motor (SRM) with a versatile simulator system. The SSRCS design is based on hybrid rocket motor principles. The simulator uses a solid fuel and a gaseous oxidizer. Verification of the feasibility of a SSRCS system as a test bed was completed using flow field and system analyses, as well as empirical test data. A total of 27 hot firings of a subscale SSRCS motor were conducted at MSFC. Testing of the Small-scale SSRCS program was completed in October 1992. This paper, a compilation of reports from the above team members and additional analysis of the instrumentation results, will discuss the final results of the analyses and test programs.

Goldberg, Benjamin E.↗
NASA NTRS1991

Experimental investigation of a solid rocket combustion simulator

The response of solid rocket motor materials to high-temperature corrosive gases is usually accomplished by testing the materials in a subscale solid rocket motor. While this imposes the proper thermal and chemical environment, a solid rocket motor does not provide practical features that would enhance systematic evaluations such as: the ability to throttle for margin testing, on/off capability, low test cost, and a low-hazards test article. Solid Rocket Combustion Simulators (SRCS) are being evaluated by NASA to test solid rocket nozzle materials and incorporate these essential practical features into the testing of rocket materials. The SRCS is designed to generate the thermochemical environment of a solid rocket. It uses hybrid rocket motor technology in which gaseous oxygen (Gox) is injected into a chamber containing a solid fuel grain. Specific chemicals are injected in the aft mixing chamber so that the gases entering the test section match the temperature and a non-dimensional erosion factor B' to insure similarity with a solid motor. Because the oxygen flow can be controlled, this approach allows margin testing, the ability to throttle, and an on/off capability. The fuel grains are inert which makes the test article very safe to handle. The objective of this work was to establish the baseline operating characteristics of a Labscale Solid Rocket Combustion Simulator (LSRCS). This included establishing the baseline burning rates of plexiglass fuels and the evaluation of a combustion instability for hydroxy-terminated polybutadyene (HTPB) propellants. The scope of the project included: (1) activation of MSFC Labscale Hybrid Combustion Simulator; (2) testing of plexiglass fuel at Gox ranges from 0.025 to 0.200 lb/s; (3) burning HTPB fuels at a Gox rate of 0.200 lb/s using four different mixing chamber configurations; and (4) evaluating the fuel regression and chamber pressure responses of each firing.

Frederick, Robert A., Jr.↗
DOE OSTI2023

SootLib: A soot model library for combustion simulation

Soot formation in combustion is an important process that affects radiative heat transfer, flame temperatures, and emissions with health and environmental impacts. Soot formation involves complex chemistry for nucleation, growth, oxidation, and coagulation processes. The soot particles vary widely in size and accurate modeling requires representation of the particle size distribution (PSD). Modeling soot is not trivial, and is only one of several physical processes active in combustion systems. This paper presents a software package called SootLib, which is an open-source library for modeling soot formation and other aerosol systems. SootLib is written in C++, is documented with Doxygen, and is available on GitHub. The library includes several models for soot chemistry and coagulation, and it represents the PSD using either a sectional model or the method of moments (MOM). Four closure approaches for the MOM are implemented allowing up to eight moments: monodispersed, an assumed-shape lognormal distribution, the quadrature method of moments, and the method of moments with interpolative closure. SootLib provides an interface for inclusion in other combustion packages including CFD or reacting flow solvers. The range of models allows comparisons and sensitivity studies, and the modularity facilitates extension to other soot models.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗
DOE OSTI2021

Data-assisted combustion simulations with dynamic submodel assignment using random forests

This investigation outlines a data-assisted approach that employs random forest classifiers for local and dynamic submodel assignment in turbulent-combustion simulations. This method is demonstrated in simulations of a single-element GOX/GCH4 rocket combustor; a priori as well as a posteriori assessments are conducted to (i) evaluate the accuracy and adjustability of the classifier for targeting different quantities of interest (QoIs), and (ii) assess improvements, resulting from the data-assisted combustion model assignment, in predicting target QoIs during simulation runtime. Results from the a priori study show that random forests, trained with local flow properties as input variables and combustion model errors as training labels, assign three different combustion models – finite-rate chemistry (FRC), flamelet progress variable (FPV) model, and inert mixing (IM) – with reasonable classification performance even when targeting multiple QoIs. Applications in a posteriori studies demonstrate improved predictions from data-assisted simulations, in temperature and CO mass fraction, when compared with monolithic FPV calculations. An additional a posteriori data-assisted simulation of a modified configuration demonstrates that the present approach can be successfully applied to different configurations, as long as thermophysical behavior can be represented by the training data. Furthermore, these results demonstrate that this data-driven framework holds promise for dynamic combustion submodel assignments in reacting flow simulations.

42 ENGINEERING↗
NASA NTRS1990

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-2

Jet-A spray combustion has been evaluated in gas turbine combustion with the use of propane chemical kinetics as the first approximation for the chemical reactions. Here, the numerical solutions are obtained by using the KIVA-2 computer code. The KIVA-2 code is the most developed of the available multidimensional combustion computer programs for application of the in-cylinder combustion dynamics of internal combustion engines. The released version of KIVA-2 assumes that 12 chemical species are present; the code uses an Arrhenius kinetic-controlled combustion model governed by a four-step global chemical reaction and six equilibrium reactions. Researchers efforts involve the addition of Jet-A thermophysical properties and the implementation of detailed reaction mechanisms for propane oxidation. Three different detailed reaction mechanism models are considered. The first model consists of 131 reactions and 45 species. This is considered as the full mechanism which is developed through the study of chemical kinetics of propane combustion in an enclosed chamber. The full mechanism is evaluated by comparing calculated ignition delay times with available shock tube data. However, these detailed reactions occupy too much computer memory and CPU time for the computation. Therefore, it only serves as a benchmark case by which to evaluate other simplified models. Two possible simplified models were tested in the existing computer code KIVA-2 for the same conditions as used with the full mechanism. One model is obtained through a sensitivity analysis using LSENS, the general kinetics and sensitivity analysis program code of D. A. Bittker and K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. The other model is based on the work published by C. K. Westbrook and F. L. Dryer.

Nguyen, H. L.↗
DOE OSTI2025

In situ multi-tier auto-ignition detection applied to dual-fuel combustion simulations

Here we use an anomaly detection methodology that is centered on analyzing fourth-order joint moments (co-kurtosis), particularly focusing on its application in auto-ignition of combustion problems with large numbers of species. Unsupervised anomaly detection is challenging to generalize across problem types and domains. A recent technique, centered on analyzing information in the fourth-order joint moment co-kurtosis, has shown promise, especially for high-dimensional scientific data. In this work we present developments to the co-kurtosis based anomaly detection method needed to make it effective and scalable for large-scale distributed scientific data, such as those generated by massively parallel simulations. An in situ co-kurtosis algorithm is employed as the anomaly detection method for identifying ignition kernels in simulations of turbulent combustion. Here, we extend an existing methodology which identifies regions of the domain where anomalies are present, and add another tier of anomaly detection where the individual samples contributing to the anomaly are identified. We apply this algorithm on-the-fly to a variety of turbulent reacting flow problems and compare it to the widely used (but significantly more expensive) chemical explosive mode analysis (CEMA). We demonstrate the ability of the method to detect and identify the onset of low and high temperature ignition which can be used for computational steering, as chemical and combustion anomalies occur intermittently at spatio-temporal locations unknown a priori. Finally, we apply our lightweight in situ algorithm to an exascale high-fidelity simulation with a total of 2.4 Trillion degrees of freedom, performed using an adaptive mesh refinement solver. Furthermore, through a scalability analysis, we show that the relative computational cost of this in-situ anomaly detection algorithm compared to an iteration of the reacting flow solver is negligible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
NASA NTRS1999

A Sectored-One-Dimensional Model for Simulating Combustion Instabilities in Premix Combustors

A one-dimensional, CFD based combustor simulation has been developed that exhibits self-excited, thermoacoustic oscillations in premixed combustor geometries that typically have large, abrupt changes in cross sectional area. The combustor geometry is approximated by dividing it into a finite number of one-dimensional sectors. Within each sector, the equations of motion are integrated numerically, along with a species transport and a reaction equation. Across the sectors, mass and energy are conserved, and momentum loss is prescribed using appropriately compatible boundary conditions that account for the area change. The resulting simulation and associated boundary conditions essentially represent a one-dimensional, multi-block technique. Details of the simulation code are presented herein. Results are then shown comparing experimentally observed and simulated operation of a particular combustor rig that exhibited different instabilities at different operating points. It will be shown that the simulation closely matched the rig data in oscillation amplitudes, frequencies, and operating points at which the instabilities occurred. Finally, advantages and limitations of the simulation technique are discussed.

Paxson, Daniel E.↗
NASA NTRS1990

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-II

Numerical solutions of the Jet-A spray combustion were obtained by means of the KIVA-II computer code after Jet-A properties were added to the 12 chemical species the program had initially contained. Three different reaction mechanism models are considered. The first model consists of 131 reactions and 45 species; it is evaluated by comparing calculated ignition delay times with available shock tube data, and it is used in the evaluation of the other two simplified models. The simplified mechanisms consider 45 reactions and 27 species and 5 reactions and 12 species, respectively. In the prediction of pollutants NOx and CO, the full mechanism of 131 reactions is considered to be more reliable. The numerical results indicate that the variation of the maximum flame temperature is within 20 percent as compared with that of the full mechanism of 131 reactions. The chemical compositions of major components such as C3H8, H2O, O2, CO2, and N2 are of the same order of magnitude. However, the concentrations of pollutants are quite different.

Nguyen, H. L.↗
DOE OSTI2025

Collaborative: in situ visual analytics technologies for extreme scale combustion simulations

This project aims to drastically enhance the usability of in situ analysis and visualization for extreme-scale scientific simulations. Current exascale computing capabilities promise to offer greater predictive ability of simulations and to further push the frontiers of science and technology. However, to validate the simulation output at extreme scale, examine the modeled phenomena, and discover previously unknowns from the output data, the output must be reduced or transformed in situ as it is being generated during the simulation such that the amount of data to examine and store is kept to a minimum. Such in situ approaches allow us to process and analyze the data and any embedded geometry to an extent that would be prohibitively expensive, if not impossible, to perform as a post hoc task. While in situ processing has been demonstrated to be a feasible and promising approach, its full potential has not yet been leveraged. In this project, we have developed comprehensive enhancements to in situ technology based on probability distributions in data. Our research focuses on jointly developing new ways of interacting with massive statistical samples while creatively utilizing new state-of-the-art computational resources to push the boundaries of in situ exploration. Moreover, we have developed new time-dependent techniques to enable previously unattainable capabilities in areas such as intelligent simulation steering and precise feature identification. We have experimentally studied our design and implementation at NERSC and OLCF, and are able to leverage existing in situ infrastructures whenever possible. While the exemplar in this project is combustion, many other fields for which turbulent transport is important, e.g., fusion, climate, astrophysics among others, encounter similar issues as simulations scale up to the exascale. This project shows its potential to generate high impact on DOE missions since the resulting technology promises to improve scientists’ ability to rapidly and correctly interpret and tune extreme-scale simulations, leading to new scientific understanding and advancements.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

Validation of a Pressure-Based Combustion Simulation Tool Using a Single Element Injector Test Problem

The traditional design and analysis practice for advanced propulsion systems, particularly chemical rocket engines, relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment by non-CFD specialists. A computational tool, called Loci-STREAM is being developed for this purpose. It is a pressure-based, Reynolds-averaged Navier-Stokes (RANS) solver for generalized unstructured grids, which is designed to handle all-speed flows (incompressible to hypersonic) and is particularly suitable for solving multi-species flow in fixed-frame combustion devices. Loci-STREAM integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective of the ongoing work is to develop a robust simulation capability for combustion problems in rocket engines. As an initial step towards validating this capability, a model problem is investigated in the present study which involves a gaseous oxygen/gaseous hydrogen (GO2/GH2) shear coaxial single element injector, for which experimental data are available. The sensitivity of the computed solutions to grid density, grid distribution, different turbulence models, and different near-wall treatments is investigated. A refined grid, which is clustered in the vicinity of the solid walls as well as the flame, is used to obtain a steady state solution which may be considered as the best solution attainable with the steady-state RANS methodology. From a design point of view, quick turnaround times are desirable; with this in mind, coarser grids are also employed and the resulting solutions are evaluated with respect to the fine grid solution.

Thakur, Siddarth↗
DOE OSTI2021

Transforming Energy through Computational Excellence. Exascale Computing: Combustion; Simulating Effects of Fuel Injection Location in Supersonic Jet Engines

Computational tractable simulations using an adaptive-mesh-refinement solverfor compressible reacting flows help researchers understand how variations in fuel injection location within the supersonic flow cavity impacts combustion efficiency. By identifying the important physical determinants of the combustion processes, this study shows a promising pathway to improving flame stability and combustion efficiency, as well as reducing emissions.

adaptive mesh refinement↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
NASA NTRS1994

Solid rocket combustion simulator technology using the hybrid rocket for simulation

The hybrid rocket is reexamined in light of several important unanswered questions regarding its performance. The well-known heat transfer limited burning rate equation is quoted, and its limitations are pointed out. Several inconsistencies in the burning rate determination through fuel depolymerization are explicitly discussed. The resolution appears to be through the postulate of (surface) oxidative degradation of the fuel. Experiments are initiated to study the fuel degradation in mixtures of nitrogen/oxygen in the 99.9 percent/0.1 percent to 98 percent/2 percent range. The overall hybrid combustion behavior is studied in a 2 in-diameter rocket motor, where a PMMA tube is used as the fuel. The results include detailed, real-time infrared video images of the combustion zone. Space- and time-averaged images give a broad indication of the temperature reached in the gases. A brief outline is shown of future work, which will specifically concentrate on the exploration of the role of the oxidizer transport to the fuel surface, and the role of the unburned fuel that is reported to escape below the classical time-averaged boundary layer flame.

Ramohalli, Kumar↗
DOE OSTI2020

Simulating combustion of a seven-component surrogate for a gasoline/ethanol blend including soot formation and comparison with experiments

Combustion of a seven-component surrogate for a research grade 87 octane gasoline mixed with 10% ethanol is investigated experimentally and numerically from the perspective of an isolated droplet burning under conditions that promote one-dimensional gas transport. The numerical analysis included a kinetic mechanism comprised of 398 species and 24,814 reactions and a soot model that accounted for nucleation, surface growth, coalescence/aggregation of soot particles, and luminous flame radiation. Measurements of droplet and flame diameters were made for an initial droplet diameter (D o ) of approximately 0.63 mm. The simulations agreed well with the measurements including the location of the soot shell. Preferential vaporization was revealed by simulations of the liquid concentrations in the droplet. Predicted peak soot volume fractions coincided with temperatures between 1300 K and 1400 K as a soot inception temperature. Simulations were also carried out for D o between 0.25 mm and 5 mm to explore the effect of radiation and D o on burning. Below 0.25 mm radiation was negligible and burning rates and flame temperatures converged to a single value. Increasing D o up to 1.8 mm lowered the burning rate with luminous radiation having a strong effect. When radiation was entirely removed from the model the burning rate was nearly constant. Above D o = 2 mm droplets extinguished almost immediately after ignition. The flame temperature decreased with increasing D o while it increased when radiation was omitted. As a result, the simulations show that soot precursors including polyaromatic hydrocarbons were concentrated around the soot shell.

42 ENGINEERING↗