SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The purpose of this project was to create a fire model of the facility. Results from the model will be used to establish standoff distances from combustible fuel packages in the facility, in order to protect material from a thermal event.
Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.
Accurate predictions of low-temperature heat release (LTHR) are critical for modeling auto-ignition processes in internal combustion engines. While LTHR is typically obscured by deflagration, extremely late ignition phasing can lead to LTHR prior to the spark, a behavior known as pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection spark-ignition engine was studied using 3-D computational fluid dynamics (CFD) and detailed chemical kinetics. The turbulent combustion was modeled via a hybrid approach that incorporates the G-equation model for tracking the turbulent flame front, and the well-stirred reactor model with detailed chemistry for assessing the low-temperature reactions in unburnt gas. Simulations were conducted using Co-Optima alkylate and E30 fuels at operating conditions characterized by different PSHR intensities. The predicted in-cylinder pressure and heat release rate were found to agree well with experiments. It was found the estimate of previous-cycle trapped residuals is of utmost importance for capturing PSHR correctly. A simulation best practice was developed which keeps the detailed chemistry solver active throughout the entire simulation, allowing to track the evolution of intermediate species from one cycle to the next. Following the validation, the dynamics of PSHR were analyzed in detail employing the pressure-temperature (P-T) trajectory framework. It was shown that PSHR correlated with the first-stage ignition delay of the fuel, hence showing close relation to the in-cylinder P-T trajectory and the chemical kinetics. Besides, it was indicated that LTHR is a self-limiting process that has the effect of attenuating the thermal stratification in the combustion chamber. Furthermore, it was observed the occurrence of PSHR caused the P-T trajectory of end-gas to overlap with the negative temperature coefficient region of the fuel’s ignition-delay maps. This effect was more significant in the fuel-rich regions where engine knock tendency would be generally higher, with potential implications on knock control and mitigation.
A perturbation analysis modeling approach based on a genetic algorithm was used to identify possible reaction pathways to explain previous experimental observations of the strong acceleration of syngas auto-ignition by trimethylsilanol (TMSO). Organosilicon reactions were taken from an existing chemical kinetic mechanism for tetramethylsilane pyrolysis which also contained TMSO reactions. The experimental targets for optimization were ignition delay times and peak OH radical concentrations in the range 1010–1070 K, 5 atm with dilute (4%) ϕ = 0.1 syngas mixtures in N 2 /Ar containing 0, 200, and 1000 ppm TMSO. Hundreds of thousands of trials resulted in three models which accurately predicted the experimental auto-ignition times and peak OH concentrations. Rate of production analyses indicated that TMSO forms dimethylsilanediol, which then reacts creating two “catalytic loops.” The net effect of the loops is to accelerate the reaction H 2 + O = OH + H at a faster effective rate than H 2 /O 2 kinetics alone, thus enhancing the reactivity of syngas. The formation of the loops may be attributable to the high Si–C bond energies in siloxanols, reducing rates of thermal decomposition and allowing dimethylsilanediol to react with H 2 for a longer duration compared with its alkane analog. Here, the results of the current work provide valuable direction for fundamental studies of these new hypothesized reaction pathways.
Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.
Catalytic Fast Pyrolysis (CFP) is a process that converts biomass into liquid intermediates suitable for transportation fuels by rapidly heating it in the presence of a catalyst, aiming to produce stable oils with reduced oxygen content. During CFP, the catalyst can become deactivated by the accumulation of coke, a carbon-rich deposit formed from the decomposition of biomass components. Unlike in petroleum refining, regenerating coked catalysts from biomass pyrolysis requires specific approaches due to the different chemical nature of the coke formed. An experimental technique, Temperature Programmed Oxidation (TPO), was used to study the de-coking process by gradually increasing temperature while monitoring the production of CO and CO2, which provides data for kinetic modeling. Utilizing data from TPO experiments, coke combustion kinetic model was developed to describe the rate of coke removal at different temperatures, allowing for simulation of regeneration processes. Then kinetic model is integrated into MFIX-Exa for the simulation of Biogenic Fluid Catalytic Cracker (BFCC) regenerator at different scales, enabling analysis of catalyst flow, temperature distribution, and regeneration efficiency under various operating conditions.
Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.
This project produced clear evidence that biofuels can reduce particulate emissions from road vehicles compared with petroleum-derived fuels. It also produced fundamental knowledge that can be used to select fuel compositions and engine designs that will maximize this benefit. This information has been disseminated through more than 20 peer-reviewed publications, more than 30 conference presentations, and 5 freely available internet databases and web tools. This project also provided educational opportunities to 17 students at the high school, undergraduate, and graduate levels.
n-Propylbenzene (PB, C9H12, 432K boiling point) is important as a constituent in surrogates for some transportation fuels but little is understood about its droplet burning characteristics. We present an experimental and numerical study of a fuel droplet burning under conditions where the sole means of gas transport is fuel evaporation at the droplet surface thereby creating a one-dimensional gas transport that is amenable to detailed numerical modeling. The one-dimensional gas transport that results is amenable to detailed numerical modeling. Initial droplet diameters (Do) ranged from 0.3 mm to about 5 mm to better understand the mechanisms that control burning over such a wide range of Do. For Do < 1 mm droplets were burned while anchored to 14 m SiC fibers within a sealed container in free-fall during which time the droplets were ignited by spark discharge and the burning process recorded by video cameras. Larger and free-floating droplets were studied in the spaced-based platform of the orbiting International Space Station for Do > 1.0 mm which provided unlimited times for experimental observation. Measurements of droplet, flame and soot shell diameters were made and are discussed in the presentation. The experimental results showed formation of optically thick soot clouds developing during burning that required a robust image analysis algorithm to enable measurement of droplet diameter. The measured droplet diameters were well predicted by a numerical model that included soot formation, a PB kinetic mechanism comprising 576 species and 27369 reactions, unsteady gas and liquid transport, variable properties and radiative transport. Flame diameters showed more variability compared to the simulations, depending on how the flame was defined (peak gas temperature and OH concentration). Droplets burned slower as Do was increased and the PB sooting propensity increased with Do. Flame extinction was observed for Do > 3.3 mm indicating a strong effect of radiation as Do was increased compared to the smaller droplet diameters examined where extinction was not experimentally observed. For the largest droplets examined the possibility that the observed extinction is indicative of a transition to cool flame burning is still under investigation.
This paper investigates the formation of nitric oxide (NO) and carbon monoxide (CO) emissions as a function of nitrogen (N 2 ) or carbon dioxide (CO 2 ) diluent content, added to a premixed reacting jet in crossflow. Here, reaction characteristics of a rich methane-air jet injected into a lean vitiated crossflow were analyzed at an elevated pressure of 5 atm. The jet was pre-heated and enriched with 0%, 15%, 30%, and 50% mass fraction to quantify the effect on pollutant emissions. Simulated results of the full chemistry Star-CCM+CFD model were verified with data taken in the experimental high-pressure combustion facility, which provides pressure, temperature and velocity profiles, as well as line-of-sight chemiluminescence and exit emission measurements. The analysis revealed the significant influence of the diluent to delay axial combustion; this mixing delay increased the axial flame lift-off and controlled any thermal ignition limitations. Increasing the diluent content increased the timescale for the flame to stabilize, which allowed for greater entrainment of crossflow oxidizer into the axial jet stream, and led to decreased pollutant emissions. Hence, crossflow entrainment is a critical driving force at high diluent content. Local nitric oxide (NO) emission formation in the axial stage was predicted numerically, showing the correlation between diluent addition, axial heat release, and the formation of nitric oxide pollution. The high diluent levels leaned out the jet mixture and delayed axial combustion, while reducing prompt and thermal NOx by mitigating flame hotspots and minimizing the timescale that the products remain at high temperatures.
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.
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.
The goal of the project is to improve the predictive accuracy and efficiency of turbulent combustion sub-models for pre-chamber turbulent jet ignition (TJI). This goal is achieved through the development of a hierarchically informed engineering model for turbulent combustion in TJI. The model development starts with the highest level of model description of turbulent combustion with direct numerical simulation (DNS) from which fundamental characteristics and scaling properties of turbulent premixed flame propagation under TJI relevant conditions are obtained.
Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.
With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.
Mixing controlled combustion of alcohol fuels has been identified as a promising technology based on their low propensity for particulate and NOx production, but the higher heats of vaporization and auto-ignition temperatures of these fuels make their direct use in diesel engine architectures a challenge. To realize the potential of alcohol-fueled combustion, a computational fluid dynamics (CFD) modeling framework is developed, validated, and exercised to identify designs that maximize engine thermal efficiency. To evaluate the use of thermal barrier coatings (TBCs), a simplified one-dimensional (1D) conjugate heat transfer (CHT) modeling framework is employed. The addition of the 1D CHT model only increases the computational expense by 15% relative to traditional approaches, yet offers more accurate heat transfer predictions over constant temperature boundary conditions. The validated model is then used to explore a range of injector orientations and piston bowl geometries. Using a design of experiments (DoE) approach, several designs were identified that improved fuel-air mixing, shortened the combustion duration, and increased thermal efficiency. The most promising design was fabricated and tested in a Caterpillar 1Y3700 single-cylinder oil test engine (SCOTE). Engine testing confirmed the findings from the CFD simulations and found that the co-optimized injector and piston bowl design yielded over 2-percentage point increase in thermal efficiency at the same equivalence ratio (0.96) and over 6-percentage point increase at the same engine load (10.1 bar indicated mean effective pressure (IMEP)), while satisfying design constraints for peak pressure and maximum pressure rise rate.