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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Experimental and Simulation Analysis of Binary Mixtures of Biomass and Inert Material

A literature review was conducted to identify experiments on various scales of biomass combustion specific to circulating fluidized bed (CFB) combustors. Typical feedstock and inert bed material characteristics for biomass combustion in a CFB were gathered. This information was used to design laboratory-scale, non-reacting experiments on biomass-inert bed material mixtures in a 1-in. cylindrical test rig in the Multiphase Flow Analysis Laboratory (MFAL) in Morgantown, WV. The biomass chosen was a hardwood. Pellets from the hardwood were milled and sieved to obtain three size ranges—small, medium, and large. Three mass loadings of each of the biomass samples were fluidized with glass beads as the inert bed material to study the mixing characteristics under varying fluidizing gas velocities. Hydrodynamic data including pressure drop measurements at select vertical locations and video recordings of the mixture behavior were generated in this effort.

09 BIOMASS FUELS↗
DOE OSTI2021

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗
DOE OSTI2021

An Experimental and Numerical Investigation to Characterize the Low-Temperature Heat Release in Stoichiometric and Lean Combustion

This work reports on an experimental and modeling study on the low-temperature heat release (LTHR) characteristics for three RON 90 binary blends (n-heptane blended with isooctane, toluene and ethanol) in a Cooperative Fuel Research (CFR) engine at lean and stoichiometric conditions that are representative of homogeneous charge compression ignition (HCCI) and spark-ignition (SI) end-gas combustion conditions, respectively. An analysis of the end-gas temperature-pressure (T-P) trajectories was performed to identify the intake conditions leading to similar T-P trajectories between the two lambdas for each fuel blend. A heat release analysis was then conducted for the identified cases, where fuel-to-fuel differences in LTHR were identified and found to be sensitive to the operating condition. Simulations were conducted for these cases using a recently updated chemical kinetic model and a 0-D engine model, where good qualitative and reasonable quantitative agreements in LTHR were obtained. Sensitivity analysis was also performed directly on the rates of LTHR, to understand the controlling chemical reactions of LTHR, providing further insights into the fuel-to-fuel differences. The results demonstrate the significant promoting effect of n-heptane on LTHR rates, while inhibiting effects were seen for ethanol and toluene. Also highlighted was the importance of H atom abstraction reactions from the chemistry of each fuel component, which could lead to contradictory fuel behavior depending on the locations of the H site of the abstraction reaction due to the different ensuing pathways for the primary fuel radicals.

LTHR, RON 90 binary blends, sensitivity analysis↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2021

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
DOE OSTI2022

Exploring the potential of hydrogen in decarbonizing China's light-duty vehicle market

The Chinese government has pledged to achieve overall carbon neutrality by 2060. Currently, the transportation sector contributes to about 10% of total greenhouse gas (GHG) emissions in China. Hence, China has created a well-defined energy vehicle development strategy to reduce GHG emissions from the transportation sector, further expanding into hydrogen vehicle technologies. In this study, the Transportation Energy Analysis Model (TEAM) investigates the potential of hydrogen internal combustion engine vehicles (H2- ICEVs) and fuel cell vehicles (FCEVs) as a reliable pathway towards the government's aspiration of carbon neutrality in the transportation sector. According to TEAM, by adopting FCEVs and H2-ICEVs in the vehicle market, hydrogen demand could reach 25% of the total light-duty transportation energy demand in 2050. Consequently, this will lead to an annual reduction of more than 35 million tons GHG compared to only counting on the electrification pathway in the decarbonization task. Besides, FCEVs would take longer to penetrate the light-duty vehicle market compared to H2-ICEVs, as the current fuel cell technology still requires much improvement to attain a competitive vehicle cost of production.

Life cycle analysis↗
DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗
DOE OSTI2021

The impact of the third O 2 addition reaction network on ignition delay times of neo-pentane

In this work, we studied the oxidation of neo-pentane by combining experiments, theoretical calculations, and mechanistic developments to elucidate the impact of the 3rd O 2 addition reaction network on ignition delay time predictions. The experiments are based on photoionization mass spectrometry in jet-stirred and time-resolved flow reactors allowing for sensitive detection of the keto-hydroperoxide (KHP) and keto-dihydroperoxide (KDHP) intermediates. With neo-pentane exhibiting a unique symmetric molecular structure, which consequently results only in single KHP and KDHP isomers, theoretical calculations of ionization and fragment appearance energies and of absolute photoionization cross sections enabled the unambiguous identification and quantification of the KHP intermediate. Its temperature and time-resolved profiles together with calculated and experimentally observed KHP-to-KDHP signal ratios were compared to simulation results based on a newly developed mechanism that describes the 3rd O 2 addition reaction network. A satisfactory agreement has been observed between the experimental data points and the simulation results, thus adding confidence to the model's overall performance. Finally, this mechanism was used to predict ignition delay times reported previously in shock tube and rapid compression machine experiments (J. Bugler et al., Combust. Flame 163 (2016) 138–156). While the model accurately reproduces the experimental data, simulations with and without the 3rd O 2 addition reaction network included reveal only a negligible effect on the predicted ignition delay times at 10 and 20 atm. According to model calculations, low temperatures and high pressures promote the importance of the 3rd O2 addition reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Jacobian-based Model Diagnostics and Application to Equation Oriented Modeling of a Carbon Capture System

Equation-oriented (EO) modeling has the potential to enable the effective design and optimization of the operation of advanced energy systems. However, advanced modeling of energy systems results in a large number of variables and non-linear equations, and it can be difficult to search through these to identify the culprit(s) responsible for convergence issues. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms so they can be rescaled. A further singular value decomposition can be per-formed to identify degenerate sets of equations and remaining scaling issues. This work presents an EO model of a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. The IDAES diagnostics tools were successfully applied to this flowsheet to identify problems to improve model robustness and enable the optimization of process design and operating conditions of a carbon capture system.

Allan, Douglas↗
DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗
DOE OSTI2025

Optimization of a Secondary Air Injector for a Rich-Quench-Lean (RQL) Ammonia Combustor Using Computational Fluid Dynamics

Ammonia has emerged as a promising medium for moving hydrogen around the globe due to its energy density, existing infrastructure for production, transportation, and storage, and multiple applications from fertilizer to power generation. However, one of the most significant challenges with ammonia combustion for power generation is the formation of nitrogen oxides (NOx) during combustion. Several combustion strategies have been developed to minimize the formation of NOx. One such strategy, known as rich-quench-lean (RQL), is a method that combusts NH3 in a fuel-rich environment, followed by a quick mix section and a lean burnout section. Rapid mixing of secondary air before lean burnout is thought to be important to minimize the formation of NOx. A genetic algorithm (GA) is used to parametrically vary the secondary air injector design, including the diameter, count, and angles over a constrained design space. The designs are evaluated using a non-reacting OpenFOAM model, with the objective function being the uniformity index of the secondary air (modeled as a scalar). Optimization campaigns show that traditional correlation-based designs might not be adequate. After running 414 OpenFOAM models, the optimal design is 32, 1 mm diameter, angled air injectors, increasing the uniformity index at the outlet of the quick mix zone by 36.8%. The most promising designs will be manufactured and tested in a small RQL burner setup, which is expected to lead to validation and insight into the practical usage of NH3 combustion for power generation.

ammonia combustion↗
DOE OSTI2021

Modeling, Design, and Testing of a Novel Biphasic Solvent-Enabled Absorption System for Post-Combustion Carbon Capture

A novel absorption process enabled by a new class of biphasic solvents for post-combustion carbon capture (BiCAP) is presently being developed at the University of Illinois at Urbana-Champaign. The new biphasic solvents are water-lean solvent blends that develop dual liquid phases with the absorbed CO2 highly enriched in one of the phases. These solvents are superior in CO2 capacity and have high thermal and oxidative stability compared to MEA. The BiCAP technology features a unique process configuration of multi-stage CO2 absorption and liquid–liquid phase separation during CO2 absorption, allowing continual separation and removal of the CO2-enriched liquid phase for maintaining low solvent viscosity and a fast absorption rate. A 40 kWe-scale skid has recently been constructed at the Abbott Power Plant in Champaign, IL, to test the BiCAP solvent in continuous operation with a slipstream of real post-combustion flue gas. Along with improvements in the stripping configuration including cold feed bypass, the process is expected to achieve a minimum reboiler duty of 2,210 kJ/kg CO2 captured at a stripper operating pressure of 6 bar.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Collaborative investigation of the internal flow and near-nozzle flow of an eight-hole gasoline injector (Engine Combustion Network Spray G)

The internal details of fuel injectors have a profound impact on the emissions from gasoline direct injection engines. However, the impact of injector design features is not currently understood, due to the difficulty in observing and modeling internal injector flows. Gasoline direct injection flows involve moving geometry, flash boiling, and high levels of turbulent two-phase mixing. In order to better simulate these injectors, five different modeling approaches have been employed to study the engine combustion network Spray G injector. Here these simulation results have been compared to experimental measurements obtained, among other techniques, with X-ray diagnostics, allowing the predictions to be evaluated and critiqued. The ability of the models to predict mass flow rate through the injector is confirmed, but other features of the predictions vary in their accuracy. The prediction of plume width and fuel mass distribution varies widely, with volume-of-fluid tending to overly concentrate the fuel. All the simulations, however, seem to struggle with predicting fuel dispersion and by inference, jet velocity. This shortcoming of the predictions suggests a need to improve Eulerian modeling of dense fuel jets.

ECN↗
DOE OSTI2022

Effect of fuel composition and EGR on spark-ignited engine combustion with LPG fueling: Experimental and numerical investigation

This paper presents an experimental and numerical investigation of a spark-ignited (SI) cooperative fuel research (CFR) engine fueled with different liquefied petroleum gas (LPG) fuels and exhaust gas recirculation (EGR). Here, the effects of LPG fuel composition on engine combustion characteristics are initially evaluated at two different compression ratios (CR). Results show normal combustion at CR 7 and heavy knocking combustion at CR 10 for all the tested fuels, with a more substantial impact for the LPG fuel with high proportions of n-butane species. The Livengood-Wu (LW) integral method is then used to analyze the knock occurrence risk of individual fuel based on the reactivity of the tested fuels. The introduction of EGR then demonstrates the potential of knock intensity reduction below the borderline knock limit. A zonal-based kinetic interactions study is also performed to understand the knock mitigation effectiveness of EGR over the pressure–temperature domain relevant to SI engine operation. Finally, a multidimensional, computational fluid dynamics (CFD) simulation model is shown to predict the LPG combustion characteristics and presents the evolution of in-cylinder temperature and chemical species to demonstrate the development of end-gas autoignition events without and with EGR.

02 PETROLEUM↗
DOE OSTI2020

A detailed chemical kinetic modeling and experimental investigation of the low- and high-temperature chemistry of n-butylcyclohexane

Chemical kinetic models of gasoline, jet, and diesel fuels and their mixtures with bioderived fuels are needed to assess fuel property effects on efficiency, emissions, and other performance metrics in internal combustion and gas turbine engines. As these real fuels have too many fuel components to be included in a chemical kinetic model, surrogate fuels containing fewer components are used to represent them. These surrogate fuels mimic the chemical classes or molecular structures contained in the real fuel. One of the important chemical classes in gasoline, jet, and diesel fuels comprises cyclohexanes. Cyclohexanes comprise about 30% or more by weight in diesel fuel. Also, Mueller et al proposed n-butylcyclohexane (nBCH) as a component in a nine-component surrogate palette to represent the ignition properties, distillation curve, density, and molecular structures of a diesel certification fuel. Herein, experimental measurements of the ignition delay times (IDTs) of nBCH in a shock tube and in a rapid compression machine are reported over a wide range of temperature, pressure, and equivalence ratio important for enabling the validation of a chemical kinetic model for nBCH for combustion in diesel engines. The range of conditions are temperatures of 630–1420 K, pressures of 10, 30, and 50 bar, and equivalence ratios of 0.3, 0.5, 1.0, and 2.0 in ‘air’. A detailed chemical model is developed for nBCH to simulate its ignition at both low and high-temperature conditions and at relevant elevated pressures. The experimentally measured IDTs are used to improve and validate the chemical kinetic model.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

An improved detailed chemical kinetic model for C3-C4 linear and iso-alcohols and their blends with gasoline at engine-relevant conditions

Propanol and butanol isomers have received significant research attention as promising fuel additives or neat biofuels. Robust chemical kinetic models are needed that can provide accurate and efficient predictions of combustion performance across a wide range of engine relevant conditions. Herein, this study seeks to improve the understanding of ignition and combustion behavior of pure C3-C4 linear and iso-alcohols, and their blends with gasoline at engine-relevant conditions. In this work, a kinetic model with improved thermochemistry and reaction kinetics was developed based on recent theoretical calculations of H-atom abstraction and peroxy radical reaction rates. Kinetic model validations are reported, and the current model reproduces the ignition delay times of the C3 and C4 alcohols well. Variations in reactivity over a wide range of temperatures and other operating conditions are also well predicted by the current model. Recent ignition delay time measurements from a rapid compression machine of neat iso-propanol and iso-butanol [Cheng et al., Proc. Combust Inst. (2020)] and blends with a research grade gasoline [Goldsborough et al., Proc. Combust Inst. (2020)] at elevated pressure (20–40 bar) and intermediate temperatures (780–950 K) were used to demonstrate the accuracy of the current kinetic model at conditions relevant to boosted spark-ignition engines. The effects of alcohol blending with gasoline on the autoignition behavior are discussed. The current model captures the suppression of reactivity in the low-temperature and negative-temperature-coefficient (NTC) region when either isopropanol and isobutanol are added to a research grade gasoline. Sensitivity and reaction flux analysis were performed to provide insights into the relevant fuel chemistry of the C3-C4 alcohols.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

CFD unified approach under Eulerian–Lagrangian framework for methanol and gasoline direct injection sprays in evaporative and flash boiling conditions

Innovative synthetic fuels for advanced propulsion systems, such as methanol and ammonia, and synthetic blended fuels (E00, E10, and E30), known for their high volatility, are often injected directly into combustion chambers. It follows that Eulerian–Lagrangian spray models need to accurately capture the spray collapse as a consequence of flash boiling onset and be capable of proficiently handling the preferential evaporation of multi-component fuels in evaporative scenarios. So, we performed the assessment of an Eulerian–Lagrangian CFD code for simulating methanol and E00 gasoline blend sprays in both early and late injection conditions involving flash boiling conditions and preferential evaporation. The adoption of an effervescent breakup model and of a non-equilibrium phase transition model for the discrete phase allows the adoption of a setup that is almost completely free from specific constant tuning, especially for what concerns the breakup model. We validated the simulations using experimental PLV maps of methanol and E00 sprays issued from the ECN Spray M injector. The results highlight a significantly different morphology of the methanol spray compared to the E00 one under late injection conditions. Under stratified combustion, low-volatile fuels are likely to be ignited first, and the flame propagates toward the high-volatile fuels. In conclusion, the spray collapse was also correctly reproduced, inducing the presence of a low-pressure zone and modifying the spray morphology.

E00↗