SEARCH · Engineering Papers

Results for “Combustion modeling”

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 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Atomization modeling in a multiphase flow environment and comparison with experiments

An atomization model based on Reitz's instability wave analysis has been implemented into the ARICC3D multiphase CFD combustion code. Preliminary test runs with cold non-evaporating liquid jet and coaxial gas-liquid atomization cases appeared to have verified basic performance of the model, generating realistic-looking sprays. Furthermore, the extended liquid jet is explicitly resolved, and predicted jet lengths agree well with classical correlations. Fair agreement with test data is obtained for predicted spray tip penetrations and liquid mass flux radial distributions, with obvious room for improvement. Some numerical problems also appear to have resulted with the current implementation when low gas Mach number and high liquid velocities are involved.

Liang, P. Y.↗

NASA NTRS1989

Thermodynamic and fluid mechanic analysis of rapid pressurization in a dead-end tube

Three models have been applied to very rapid compression of oxygen in a dead-ended tube. Pressures as high as 41 MPa (6000 psi) leading to peak temperatures of 1400 K are predicted. These temperatures are well in excess of the autoignition temperature (750 K) of teflon, a frequently used material for lining hoses employed in oxygen service. These findings are in accord with experiments that have resulted in ignition and combustion of the teflon, leading to the combustion of the stainless steel braiding and catastrophic failure. The system analyzed was representative of a capped off-high-pressure oxygen line, which could be part of a larger system. Pressurization of the larger system would lead to compression in the dead-end line, and possible ignition of the teflon liner. The model consists of a large plenum containing oxygen at the desired pressure (500 to 6000 psi). The plenum is connected via a fast acting valve to a stainless steel tube 2 cm inside diameter. Opening times are on the order of 15 ms. Downstream of the valve is an orifice sized to increase filling times to around 100 ms. The total length from the valve to the dead-end is 150 cm. The distance from the valve to the orifice is 95 cm. The models describe the fluid mechanics and thermodynamics of the flow, and do not include any combustion phenomena. A purely thermodynamic model assumes filling to be complete upstream of the orifice before any gas passes through the orifice. This simplification is reasonable based on experiment and computer modeling. Results show that peak temperatures as high as 4800 K can result from recompression of the gas after expanding through the orifice. An approximate transient model without an orifice was developed assuming an isentropic compression process. An analytical solution was obtained. Results indicated that fill times can be considerably shorter than valve opening times. The third model was a finite difference, 1-D transient compressible flow model. Results from the code show the recompression effect but predict much lower peak temperatures than the thermodynamic model.

Leslie, Ian H.↗

DOE OSTI2024

A framework for strategic discovery of credible neural network surrogate models under uncertainty

The widespread integration of deep neural networks in developing data-driven surrogate models for high-fidelity simulations of complex physical systems highlights the critical necessity for robust uncertainty quantification techniques and credibility assessment methodologies, ensuring the reliable deployment of surrogate models in consequential decision-making. Here, this study presents the Occam Plausibility Algorithm for surrogate models (OPAL-surrogate), providing a systematic framework to uncover predictive neural network-based surrogate models within the large space of potential models, including various neural network classes and choices of architecture and hyperparameters. The framework is grounded in hierarchical Bayesian inferences and employs model validation tests to evaluate the credibility and prediction reliability of the surrogate models under uncertainty. Leveraging these principles, OPAL-surrogate introduces a systematic and efficient strategy for balancing the trade-off between model complexity, accuracy, and prediction uncertainty. The effectiveness of OPAL-surrogate is demonstrated through two modeling problems, including the deformation of porous materials for building insulation and turbulent combustion flow for ablation of solid fuels within hybrid rocket motors.

42 ENGINEERING↗

NASA NTRSDate not supplied

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES↗

DOE OSTI2024

Development of a rate-based ENRTL-RK process model for a water-lean solvent

Advanced water-lean solvents (WLS) for post-combustion CO2 capture offer several advantages over the aqueous amine solvents . WLS have lower parasitic energy penalty, lower corrosion, lower temperature and high-pressure CO2 regeneration leading to lower cost of CO2 capture. RTI International, with funding from the US Department of Energy, has been developing its novel water-lean solvent, that has shown specific reboiler duty of 2.3 GJ/t-CO2 at the 60-kWe pilot testing unit (Tiller Plant, SINTEF, Norway) and 2.6 GJ/t-CO2 at the engineering scale testing system (12 MWe) at the Technology Centre Mongstad (TCM) in Norway. All heat duties, including the one from TCM testing, were consistent with Aspen Plus modeling of the specific configuration of each test plant. This work focuses on the development of a detailed process model using in-house laboratory measurements and process data at pilot scale. The eNTRL-RK model used in this work is based on an unsymmetric activity coefficient model with the reference states chosen to be pure liquids for solvents and ideal dilute solution at unit solute molality (resulting in activity coefficient of unity at infinite dilution) for electrolytes. It uses the Redlich-Kwong equation of state for vapor phase properties and Henry’s law for solubility of supercritical gases. The model was validated using process data from the pilot-scale campaign at the Tiller plant, and the engineering scale test campaign at TCM. Data on CO2 capture rate, absorber, and regenerator temperature profiles and specific reboiler duties from two different test campaigns at Tiller and TCM, were used to further refine and validate the model and the model compares favorably to experimental data. The validation results against TCM campaign will be presented in this work.

CO2 capture↗

NASA NTRS1996

CFD Analysis of the 24-inch JIRAD Hybrid Rocket Motor

A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.

Liang, Pak-Yan↗

DOE OSTI2021

Advanced Light-Duty Spark Ignition Engine Research: Co-Optimization of Fuels and Engines and Partnership to Advance Combustion Engines (FY2020 Annual Progress Report)

This report covers recent progress on research tasks that support both the Co-Optimization of Fuels and Engines (Co-Optima) initiative and the Partnership to Advance Combustion Engines (PACE) consortium. The Co-Optima tasks further the science-base needed by industry stakeholders to co-evolve the next generation of highly efficient direct injection spark ignition (DISI) engines and new gasoline-type fuels. The research emphasis is on fuel effects on multimode spark ignition (SI) engine operation, which uses traditional non-dilute stoichiometric operation for peak load and power but reverts to lean operation at lower loads to provide higher fuel economy. This work focuses on determining desirable fuel specifications in terms of well-established metrics like research octane number (RON) and motor octane number, but it also involves the assessment of new fuel metrics, including fuel sooting propensity and phi-sensitivity. The PACE task supports the development of predictive computational fluid dynamics (CFD) modeling, which promises to unlock new strategies for high-efficiency combustion while minimizing tailpipe emissions. Here, the primary fuel is a regular E10 gasoline (i.e., a regular gasoline blend containing 10% ethanol), and focus is on fuel-spray dynamics and soot emissions. Soot-formation pathways are studied to determine how the pathways change with injection strategies and the thermal state of the engine (i.e., cold-starting vs. fully warmed-up operation). This PACE task also contributed to the development of an optimal E10 gasoline surrogate fuel, as reported in detail elsewhere

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1995

The effects of detailed chemistry and transport on microgravity droplet combustion

A brief overview of recent advances in the theoretical study of microgravity droplet combustion is presented. Much of this work has centered on the development and utilization of sphero-symmetric transient numerical models which consider detailed gas phase chemistry and transport as well as energy and/or species transport within a regressing condensed phase. Numerical results for microgravity combustion and vaporization of methanol, methanol/water, heptane, and heptane/hexadecane droplets are summarized along with refinements in chemical kinetics and the development of a new two-dimensional axi-symmetric model.

Marchese, A. J.↗

DOE OSTI2022

Development of A Multi-Spectral Pyrometry Sensor for High-Speed Transient Surface-Temperature Measurements in Combustion-Relevant Harsh Environments

Accurate and high-speed transient surface-temperature measurements of combustion devices including internal combustion (IC) engines, gas turbines, etc., provide validation targets and boundary conditions for computational fluid dynamics models, and are broadly relevant to technology advancements such as performance improvement and emissions reduction. Development and demonstration of a multi-infrared-channel pyrometry-based optical instrument for high-speed surface-temperature measurement is described. The measurement principle is based on multi-spectral radiation thermometry (MRT) and uses surface thermal radiation at four discrete spectral regions and a corresponding emissivity model to obtain surface temperature via non-linear least squares (NLLS) optimization. Rules of thumb for specifying the spectral regions and considerations to avoid interference with common combustion products are developed; the impact of these along with linear and non-linear MRT analysis are assessed as a function of temperature and signal-to-noise ratio. A multi-start method to determine the MRT-solution global optimum is described and demonstrated. The resulting multi-channel transient pyrometry instrument is described along with practical considerations including optical-alignment drift, matching intra-channel transient response, and solution-confidence indicators. The instrument demonstrated excellent >97% accuracy and >99% 2-sigma precision over the 400–800 °C range, with ~20 µs (50 kHz, equivalent to 0.2 cad at 2000 RPM IC-engine operation) transient response in the bench validation.

47 OTHER INSTRUMENTATION↗

NASA NTRS2003

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗

DOE OSTI2020

Learning reference governor for cycle-to-cycle combustion control with misfire avoidance in spark-ignition engines at high exhaust gas recirculation–diluted conditions

Cycle-to-cycle feedback control is employed to achieve optimal combustion phasing while maintaining high levels of exhaust gas recirculation by adjusting the spark advance and the exhaust gas recirculation valve position. The control development is based on a control-oriented model that captures the effects of throttle position, exhaust gas recirculation valve position, and spark timing on the combustion phasing. Under the assumption that in-cylinder pressure information is available, an adaptive extended Kalman filter approach is used to estimate the exhaust gas recirculation rate into the intake manifold based on combustion phasing measurements. The estimation algorithm is adaptive since the cycle-to-cycle combustion variability (output covariance) is not known a priori and changes with operating conditions. A linear quadratic regulator controller is designed to maintain optimal combustion phasing while maximizing exhaust gas recirculation levels during load transients coming from throttle tip-in and tip-out commands from the driver. During throttle tip-outs, however, a combination of a high exhaust gas recirculation rate and an overly advanced spark, product of the dynamic response of the system, generates a sequence of misfire events. In this work, an explicit reference governor is used as an add-on scheme to the closed-loop system in order to avoid the violation of the misfire limit. The reference governor is enhanced with model-free learning which enables it to avoid misfires after a learning phase. Experimental results are reported which illustrate the potential of the proposed control strategy for achieving an optimal combustion process during highly diluted conditions for improving fuel efficiency.

42 ENGINEERING↗

DOE OSTI2021

Role of Inlet Boundary Conditions on Fuel-Air Mixing at Supercritical Conditions

Advanced engine design and alternative fuels present the possibility of fuel injection at purely supercritical conditions in diesel engines and gas turbines. The complex interactions that govern this phenomenon still need significant research, particularly the boundary conditions for fuel injection are critical for accurate simulation. However, the flow inside the injector itself is often omitted to reduce the computational efforts, and thus, velocity, mass flux, or total pressure is specified at the injector exit (or domain inlet), often with simplified velocity profiles and turbulence levels. This simplified inlet boundary treatment has minimal effects on results for conventional fuel injection conditions, however, the validity of this approach at supercritical conditions has not been assessed. Comprehensive real-gas and binary fluid mixing models have been implemented for computational fluid dynamics (CFD) analysis of fuel-air mixing at supercritical conditions. The model is verified using prior CFD results from the literature. The model is used to investigate the effects of the shape of axial velocity and mass fraction profiles at the inlet boundary with the goal to improve the comparison of predictions to experimental data. Results show that the boundary conditions have a significant effect on the predictions, and none of the cases match precisely with experimental data. Furthermore, the study reveals that the physical location of the inlet boundary might be difficult to infer correctly from the experiments and highlights the need for high-quality, repeatable measurements at supercritical conditions to support the development of relevant high-fidelity models for fuel-air mixing.

02 PETROLEUM↗

NASA NTRS2007

Combustion

An overview of the emissions related research being conducted as part of the Fundamental Aeronautics Subsonics Fixed Wing Project is presented. The overview includes project metrics, milestones, and descriptions of major research areas. The overview also includes information on some of the emissions research being conducted under NASA Research Announcements. Objective: Development of comprehensive detailed and reduced kinetic mechanisms of jet fuels for chemically-reacting flow modeling. Scientific Challenges: 1) Developing experimental facilities capable of handling higher hydrocarbons and providing benchmark combustion data. 2) Determining and understanding ignition and combustion characteristics, such as laminar flame speeds, extinction stretch rates, and autoignition delays, of jet fuels and hydrocarbons relevant to jet surrogates. 3) Developing comprehensive kinetic models for jet fuels.

Bulzan, Dan↗

NASA NTRS1994

Sources and distribution of NO(x) in the upper troposphere at northern midlatitudes

A simple quasi 2-D model is used to study the zonal distribution of NO(x). The model includes vertical transport in form of eddy diffusion and deep convection, zonal transport by a vertically uniform wind, and a simplified chemistry of NO, NO2 and HNO3. The NO(x) sources considered are surface emissions (mostly from the combustion of fossil fuel), lightning, aircraft emissions, and downward transport from the stratosphere. The model is applied to the latitude band of 40 deg N to 50 deg N during the month of June; the contributions to the zonal NO(x) distribution from the individual sources and transport processes are investigated. The model predicted NO(x) concentration in the upper troposphere is dominated by air lofted from the polluted planetary boundary layer over the large industrial areas of Eastern North America and Europe. Aircraft emissions are also important and contribute on average 30 percent. Stratospheric input is minor about 10 percent, less even than that by lightning. The model provides a clear indication of intercontinental transport of NO(x) and HNO3 in the upper troposphere. Comparison of the modelled NO profiles over the Western Atlantic with those measured during STRATOZ 3 in 1984 shows good agreement at all altitudes.

Rohrer, Franz↗

DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗

NASA NTRS1973

Fuel droplet burning rates at high pressures.

Combustion of methanol, ethanol, propanol-1, n-pentane, n-heptane, and n-decane was observed in air under natural convection conditions, at pressures up to 100 atm. The droplets were simulated by porous spheres, with diameters in the range from 0.63 to 1.90 cm. The pressure levels of the tests were high enough so that near-critical combustion was observed for methanol and ethanol. Due to the high pressures, the phase-equilibrium models of the analysis included both the conventional low-pressure approach as well as high-pressure versions, allowing for real gas effects and the solubility of combustion-product gases in the liquid phase. The burning-rate predictions of the various theories were similar, and in fair agreement with the data. The high-pressure theory gave the best prediction for the liquid-surface temperatures of ethanol and propanol-1 at high pressure. The experiments indicated the approach of critical burning conditions for methanol and ethanol at pressures on the order of 80 to 100 atm, which was in good agreement with the predictions of both the low- and high-pressure analysis.

Canada, G. S.↗

DOE OSTI2023

Soot Formation and Ignition Characteristics of Ethanol/Gasoline Blends in a Rapid Compression Machine

With the ever-increasing demand for sustainable energy, alcohol fuels have garnered interest for use in heavy duty engines. The significant infrastructure for ethanol production and blending of ethanol with gasoline make these fuels/fuel blends desirable candidates. However, development of heavy duty engine technology that is capable of burning alcohol fuels while retaining the advantages of traditional diesel combustion requires an improved understanding of the soot formation for these fuels under conditions relevant to mixing-controlled combustion. This work uses an extinction diagnostic to study the sooting tendency of ethanol and gasoline/ethanol blends ranging from E10 to E98 during ignition in a homogeneous environment. Experiments were conducted in a rapid compression machine (RCM) for compressed conditions of 20 ± 1 bar and an approximately constant temperature (± 10K) which was unique for each fuel. For a given soot volume fraction, a linear relationship was observed between ethanol content and the equivalence ratio in which that soot volume fraction was formed. Accounting for the oxygenated nature of ethanol, E85 and E98 fuels produced similar amounts of soot at a given Φ ox, , suggesting other factors outside of fuel oxygen content, such as fuel morphology, impact soot formation. Ignition delay data is reported for compressed pressures of 20 ± 1 bar and compressed temperatures ranging from 633 – 670 K for E10 and 771 – 789 K for E98. Varying pressures for E10 and E98 at conditions producing similar soot volume fractions demonstrated a linear dependence of soot formation on pressure, regardless of if the pressure considered was at top dead center or peak combustion pressure. Furthermore, the data gleaned from this work will be used to select soot models and chemical kinetic mechanisms for RCM simulations to ultimately model heavy duty engine technology with the studied fuels.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Classification and computation of extreme events in turbulent combustion

In the design of practical combustion systems, ensuring safety and reliability is an important requirement. For instance, reliably avoiding lean blowout, flame flashback or inlet unstart is critical for ensuring safe operation. Currently, the science of predicting such events is based on prior experience, limited modeling or diagnostic tools and purely statistical approaches. Even though computational and experimental tools for studying combustion devices have vastly advanced in the last three decades, the analysis of such failure events has not been pursued widely. While the use of data for model development and calibration is being widely accepted, the extension to failure events introduces numerous challenges. In particular, the focus here is on so-called data-poor problems, where the cost of generating data is extremely high and is not easily amenable to existing computational and experimental approaches. Data-poor problems are particularly relevant when related to extreme events (also called anomalous events) that can lead to catastrophic failure of the system. It is argued that transient events that describe such failure can have different causal mechanisms. To develop the scientific inference process, a classification of such problems is used to determine specific modeling paths as well as computational tools needed. Research opportunities in the emerging field of extreme event prediction are highlighted in order to identify critical and immediate needs.

97 MATHEMATICS AND COMPUTING↗