Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67

Study of high speed combustion flows by laser velocimetry

The feasibility of laser velocimetry in a high temperature jet was assessed in a model of an aircraft engine combustor. Experiments show that the problems encountered in measuring combustion flow can flow can be overcome by a carefully designed optical set-up and an appropriate signal processing and data acquisition system. Laser Doppler velocimetry provides useful information about coherent structures in hot free jets. The measurements agree with measurements in an isothermal jet.

Schaefer, H. J.↗

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NO x emissions under low-load gasoline compression ignition conditions

Here, in the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NO x ), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NO x sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NO x emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NO x and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NO x emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NO x emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NO x decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NO x emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Computational fluid dynamics↗

Prediction of the Thrust Performance and the Flowfield of Liquid Rocket Engines

In an effort to improve the current solutions in the design and analysis of liquid propulsive engines, a computational fluid dynamics (CFD) model capable of calculating the reacting flows from the combustion chamber, through the nozzle to the external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea level, was investigated as a sample case. The CFD model, FDNS, is a pressure based, non-staggered grid, viscous/inviscid, ideal gas/real gas, reactive code. An adaptive upwinding differencing scheme is employed for the spatial discretization. The upwind scheme is based on fourth order central differencing with fourth order damping for smooth regions, and second order central differencing with second order damping for shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a PARASOL finite rate chemistry algorithm using the point implicit method. The computed flow results and performance compared well with those of other standard codes and engine hot fire test data. In addition, the transient nozzle flowfield calculation was also performed to demonstrate the ability of FDNS in capturing the flow separation during the startup process.

Wang, T.-S.↗

Biochemical Conversion of Lignocellulosic Biomass to Hydrocarbon Fuels and Products: 2019 State of Technology and Future Research

The purpose of this report is to benchmark the latest experimental developments across a number of potential bioconversion pathways as quantified by modeled minimum fuel selling prices (MFSPs), as a measure of current status relative to those final targets. For this state of technology, TEA models were run for two separate biological conversion pathways to fuels, based on available data for integrated biomass deconstruction and hydrolysate processing; namely carboxylic acids (primarily butyric acid) and diols (2,3-butanediol [BDO]), reflecting NREL’s recently-published 2018 biochemical design report focused on those two pathways. The models were run across three scenarios for lignin utilization, namely combustion, conversion to coproducts based on “base case” performance with biomass hydrolysate, and conversion to coproducts based on “high” performance demonstrated with model lignin monomer components.

09 BIOMASS FUELS↗

Biochemical Conversion of Lignocellulosic Biomass to Hydrocarbon Fuels and Products: 2020 State of Technology and Future Research

The annual State of Technology (SOT) assessment is an essential activity for biochemical platform research. It allows the impact of research progress to be quantified in terms of economic improvements in the overall cellulosic biofuel production process for a particular conversion pathway. As such, initial benchmarks can be established for currently demonstrated performance and progress can be tracked towards out-year goals to ultimately demonstrate cost-competitive cellulosic biofuel technology. The purpose of this report is to benchmark the latest experimental developments across a number of potential bioconversion pathways as quantified by modeled minimum fuel selling prices (MFSPs), as a measure of current status relative to those final targets. For this state of technology, TEA models were run for two separate biological conversion pathways to fuels, based on available data for integrated biomass deconstruction and hydrolysate processing; namely carboxylic acids (primarily butyric acid) and diols (2,3-butanediol [BDO]), reflecting NREL's recently-published 2018 biochemical design report focused on those two pathways. The models were run across three scenarios for lignin utilization, namely combustion, conversion to coproducts based on "base case" performance with biomass hydrolysate, and conversion to coproducts based on "high" performance demonstrated with model lignin monomer components. Relative to prior 2019 benchmarks, the 2020 SOT achieved further MFSP reductions of approximately $1/gallon gasoline equivalent, primarily attributed to improvements in pretreatment that led to both cost reductions and yield improvements.

09 BIOMASS FUELS↗

Biochemical Conversion of Lignocellulosic Biomass to Hydrocarbon Fuels and Products (2021 State of Technology and Future Research)

The annual State of Technology (SOT) assessment is an essential activity for biochemical platform research. It allows the impact of research progress to be quantified in terms of economic improvements in the overall cellulosic biofuel production process for a particular conversion pathway. As such, initial benchmarks can be established for currently demonstrated performance and progress can be tracked towards out-year goals to ultimately demonstrate cost-competitive cellulosic biofuel technology. The purpose of this report is to benchmark the latest experimental developments across a number of potential bioconversion pathways as quantified by modeled minimum fuel selling prices (MFSPs), as a measure of current status relative to those final targets. For this state of technology, TEA models were run for two separate biological conversion pathways to fuels, based on available data for integrated biomass deconstruction and hydrolysate processing; namely carboxylic acids (primarily butyric acid) and diols (2,3-butanediol [BDO]), reflecting NREL's recently-published 2018 biochemical design report focused on those two pathways. The models were run across three scenarios for lignin utilization, namely combustion, conversion to coproducts based on "base case" performance with biomass hydrolysate, and conversion to coproducts based on "high" performance demonstrated with model lignin monomer components. A key improvement reflected in the 2021 SOT is centered around making use of the latest lignin conversion data, which over the past year focused primarily on production of ß-ketoadipate (BKA) as a more optimal molecule compared to the closely-related adipic acid coproduct of prior recent focus, both in terms of superior product properties and biology, as well as reduced processing complexity (reducing two steps for sequential production of muconate followed by hydrogenation to adipic acid down to a single step for direct production of BKA). This update translated to a roughly 17% increase in mass yield of final coproduct output at a nearly four-fold increase in fermentation productivity on lignin monomers relative to prior 2020 SOT benchmarks for muconic/adipic acid production.

09 BIOMASS FUELS↗

Unsteady loads measurements in a generic high speed engine model by means of recessed transducers

Results are presented from measurements of unsteady loads during performance tests of a generic high-speed engine model, which were made using high-frequency pressure gages installed in existing calorimeter ports of the engine model and recessed into the interior wall surface in order to reduce thermal flux to the gage diaphragm. In was found that the boundary layer pressure spectra at the model wall start to deviate from their flat plate counterpart at a short distance into the model inlet, which suggests the contributions to the spectra from the shock/boundary layer interaction. It was also found that significant levels of combustion noise propagate up through the subsonic portion of the boundary layer well into the inlet region. At the condition of an unstart, the combustion noise apparently couples with the acoustic modes of the model to cause acoustic 'hot spots' well upstream of the combustor.

Parrott, Tony L.↗

Interannual Variability and Trends of Combustion Aerosol and Dust in Major Continental Outflows Revealed by MODIS Retrievals and CAM5 Simulations During 2003-2017

Emissions and long-range transport of mineral dust and combustion-related aerosol from burning fossil fuels and biomass vary from year to year, driven by the evolution of the economy and changes in meteorological conditions and environmental regulations. This study offers both satellite and model perspectives of interannual variability and possible trend of combustion aerosol and dust in major continental outflow regions over the past 15 years (2003-2017). The decade-long record of aerosol optical depth (AOD, denoted as t), separately for combustion aerosol (τ(sub c)) and dust (τ(sub d)), over global oceans is derived from the Collection 6 aerosol products of the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard both Terra and Aqua. These MODIS/Aqua datasets, complemented by aerosol source-tagged simulations using the Community Atmospheric Model Version 5 (CAM5), are then analyzed to understand the interannual variability and potential trend of τ(sub c) and τ(sub d) in the major continental outflows. Both MODIS and CAM5 consistently yield a similar decreasing trend of -0.017 to - 0.020 decade(exp -1) for τ(sub c) over the North Atlantic Ocean and the Mediterranean Sea that is attributable to reduced emissions from North America and Europe, respectively. On the contrary, both MODIS and CAM5 display an increasing trend of +0.017 to +0.036 decade(exp -1) for τ(sub c) over the tropical Indian Ocean, the Bay of Bengal, and the Arabian Sea, which reflects the influence of increased anthropogenic emissions from South Asia and Middle East in the last two decades. Over the northwestern Pacific Ocean that is often affected by East Asian emissions of pollution and dust, the MODIS retrievals show a decreasing trend of - 0.021 decade(exp -1) for τ(sub c) and -0.012 decade(exp -1) for τ(sub d), which is however not reproduced by the CAM5 model. In other outflow regions strongly influenced by biomass burning smoke or dust, both MODIS retrievals and CAM5 simulations show no statistically significant trends; and the MODIS observed interannual variability is usually larger than that of the CAM5 simulation.

Hongbin Yu↗

Modeling the Uptake of Hydrogen Chloride onto Interior Spacecraft Materials

Hydrogen chloride (HCl) is a major combustion product from the pyrolysis of polyvinyl chloride (PVC) insulated electrical wires, a common spacecraft fire safety concern. Models at two different scales were developed to predict HCl uptake on anodized, chromate conversion coated (Iridite), and bare aluminum surfaces, as well as on Nomex fabric: a macroscopic one-step global surface reaction model where all the active sites are on the exterior surface, and a pore model where the interior active sites deeper into the oxide layer can also be accessed by HCl. Experiments were performed to calibrate kinetic and diffusion constants in the models. A cast acrylic test cell was used to measure the differences between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. For the materials with a thin (< 200 µm) or no oxide layer, the macroscopic surface reaction model adequately predicts the experimental measurements. For the anodized aluminum with a thicker oxide layer, the pore model provided a better match to experimental results. The results will be discussed with respect to the spacecraft fire safety project (Saffire).

modeling↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗

Analysis of the effects of fuel spray characteristics on NO(x) formation

Results from experimental and analytical studies on NO(x) formation in a one-dimensional monodisperse spray combustion system are reported. In a previous investigation, an optimum droplet diameter was reported for minimizing NO(x) emissions. In this paper, the effects of enhanced prevaporization and droplet interactions in addition to droplet size effects on NO(x) formation are investigated by air preheating experiments. These phenomena have been incorporated into an analytical model which predicts NO(x) formed from the combustion of monodisperse fuel sprays. Both, experimental and computed results indicate the importance of the extent of prevaporization and droplet interactions on NO(x) production. Preflame vaporization controls the gas phase stoichiometry which has a significant effect on the volume of the hot gases surrounding a fuel droplet where NO(x) is formed. On the other hand, the release of fuel vapor from the droplet to the reaction zone is influenced by droplet interactions which reduce the volume of the hot gases and NO(x) production by bringing the flame closer to the droplet surface.

Sarv, H.↗

System Modeling and Diagnostics for Liquefying-Fuel Hybrid Rockets

A Hybrid Combustion Facility (HCF) was recently built at NASA Ames Research Center to study the combustion properties of a new fuel formulation that burns approximately three times faster than conventional hybrid fuels. Researchers at Ames working in the area of Integrated Vehicle Health Management recognized a good opportunity to apply IVHM techniques to a candidate technology for next generation launch systems. Five tools were selected to examine various IVHM techniques for the HCF. Three of the tools, TEAMS (Testability Engineering and Maintenance System), L2 (Livingstone2), and RODON, are model-based reasoning (or diagnostic) systems. Two other tools in this study, ICS (Interval Constraint Simulator) and IMS (Inductive Monitoring System) do not attempt to isolate the cause of the failure but may be used for fault detection. Models of varying scope and completeness were created, both qualitative and quantitative. In each of the models, the structure and behavior of the physical system are captured. In the qualitative models, the temporal aspects of the system behavior and the abstraction of sensor data are handled outside of the model and require the development of additional code. In the quantitative model, less extensive processing code is also necessary. Examples of fault diagnoses are given.

Poll, Scott↗

A generalized analytical energy balance model for evaluating agglomeration from a binary collision of wet particles

Agglomeration of wet particles, i.e., particles coated with a thin liquid layer, is a common phenomenon in many processes like fluidized bed combustion of low rank fuels. The availability of an agglomeration model that can evaluate the outcome of a binary collision between wet particles differing in solid particle properties, liquid layer thicknesses, and initial collision (impact) speeds is essential for obtaining a comprehensive understanding on the existing processes experiencing wet particle agglomeration or for a successful development of new processes with high chances of wet particle agglomeration. This study presents a generalized agglomeration model on the basis of energy conservation before and after collision when colliding wet particles may differ in solid particle properties, liquid layer thicknesses, and impact speeds. The model was established based on the approximate values of energy losses that may happen during the collision. It incorporates body forces, solid-solid contacting, liquid capillary, and viscous contributions, as well as the liquid bridge volume effect. Predictions of the new model for collision outcomes of identical wet particles were like those from an analytical energy balance model developed recently by the group for identical wet particles. We also validated the new model by experimental data from literature. The results of a collision direction analysis indicated that the direction often has a minimal effect on the collision outcome in many practical scenarios. The results of Monte Carlo uncertainty analyses with the new model revealed that proper estimations of impact speed, under capillary limiting conditions, and thickness of coating layers and asperity heights, under viscous limiting conditions, are critical for the realistic prediction of collision outcomes at impact speeds close to critical impact speed, i.e., the minimum particle speed required for the particles to rebound.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stability analysis of a liquid fuel annular combustion chamber

The problems of combustion instability in an annular combustion chamber are investigated. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations. By evaluating these effects, parameters which cause instabilities to occur in the combustion chamber can be ascertained. It is assumed that in the annular combustion chamber, the liquid propellants are injected uniformly across the injector face, the combustion processes are distributed throughout the combustion chamber, and that no time delay occurs in the combustion processes.

Mcdonald, G. H.↗

Numerical study on spray collapse process of ECN spray G injector under flash boiling conditions

This work features a numerical study of the spray collapse process that occurs under flash-boiling, gasoline-direct-injection conditions. The Engine Combustion Network's Spray G injector was modeled, and iso-octane, n-hexane, and propane were used as fuels with the goal to investigate fuel property effects. The development of individual jets and jet-to-jet interaction were analyzed to reveal the physical mechanisms of spray collapse. The results agreed qualitatively well with the experimental data of propane available from the literature and the predictions revealed the model's ability to capture the near-nozzle shock structures witnessed in the experiments. Simulations showed that the flow development of flashing sprays was strongly tied with the fuel properties. Specifically, under-expanded jets were not observed with iso-octane due to the fuel's low vapor pressure. In the case of n-hexane, the jets were under-expanded, and a low-pressure core was observed inside the jets due to the flow expansion. With propane, the very strong expansion caused the flow to become transonic and generate a unique shock structure near the nozzle exit. The shock structure revealed the formation of primary shock cells within the jets and secondary shock cells induced by jet-to-jet interaction. Based on the results, the spray collapse under flash boiling conditions was mainly attributed to two effects, i.e., the decrease in pressure at the spray center induced by the low-pressure cores inside and between the jets, and the isolation effect due to the jet-to-jet interaction which contributed to maintain the pressure differential between the spray center and the ambient.

42 ENGINEERING↗

Multicycle large-eddy simulations of a direct-injection hydrogen-fueled optical engine

Hydrogen (H 2 ) is a carbon-free chemical energy carrier and one promising solution for achieving effective decarbonization of the transportation sector, particularly for internal combustion engines (ICEs). With a focus on ICEs, and compared to port-fuel injection, direct injection (DI) of gaseous H 2 during the compression stroke offers potential advantages, which include backfire avoidance and reduction of preignition occurrence. In these last two decades, much research, experimental and numerical, has been devoted to understanding H 2 's mixing and combustion processes in ICEs. Computational fluid dynamics modeling efforts commonly rely on unsteady Reynolds-averaged Navier Stokes (URANS) turbulence frameworks, mostly due to their computational affordability. However, many authors have pointed out the opportunity to perform large-eddy simulations (LESs) to investigate the cyclic variability of H 2 engines and assess potential advantages of using LES in place of URANS, especially for lean operation. This study addresses this knowledge gap and presents a computational fluid dynamics (CFD) study of the H 2 DI process in an optical engine operating at relatively low tumble conditions, using multicycle LESs. In conclusion, the manuscript presents a thorough validation of the results against experimental data available from the literature as well as direct comparison with URANS, demonstrating the feasibility of multicycle LESs for CFD modeling of DI H 2 -fueled ICEs.

Direct injection↗

AMReX: Block-structured adaptive mesh refinement for multiphysics applications

Block-structured adaptive mesh refinement (AMR) provides the basis for the temporal and spatial discretization strategy for a number of Exascale Computing Project applications in the areas of accelerator design, additive manufacturing, astrophysics, combustion, cosmology, multiphase flow, and wind plant modeling. AMReX is a software framework that provides a unified infrastructure with the functionality needed for these and other AMR applications to be able to effectively and efficiently utilize machines from laptops to exascale architectures. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving accurate descriptions of different physical processes in complex multiphysics algorithms. AMReX supports algorithms that solve systems of partial differential equations in simple or complex geometries and those that use particles and/or particle–mesh operations to represent component physical processes. In this article, we will discuss the core elements of the AMReX framework such as data containers and iterators as well as several specialized operations to meet the needs of the application projects. In addition, we will highlight the strategy that the AMReX team is pursuing to achieve highly performant code across a range of accelerator-based architectures for a variety of different applications.

Zhang, Weiqun↗