Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “combustion kinetics”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
NASA NTRS2012

Foundational Performance Analyses of Pressure Gain Combustion Thermodynamic Benefits for Gas Turbines

A methodology is described whereby the work extracted by a turbine exposed to the fundamentally nonuniform flowfield from a representative pressure gain combustor (PGC) may be assessed. The method uses an idealized constant volume cycle, often referred to as an Atkinson or Humphrey cycle, to model the PGC. Output from this model is used as input to a scalable turbine efficiency function (i.e., a map), which in turn allows for the calculation of useful work throughout the cycle. Integration over the entire cycle yields mass-averaged work extraction. The unsteady turbine work extraction is compared to steady work extraction calculations based on various averaging techniques for characterizing the combustor exit pressure and temperature. It is found that averages associated with momentum flux (as opposed to entropy or kinetic energy) provide the best match. This result suggests that momentum-based averaging is the most appropriate figure-of-merit to use as a PGC performance metric. Using the mass-averaged work extraction methodology, it is also found that the design turbine pressure ratio for maximum work extraction is significantly higher than that for a turbine fed by a constant pressure combustor with similar inlet conditions and equivalence ratio. Limited results are presented whereby the constant volume cycle is replaced by output from a detonation-based PGC simulation. The results in terms of averaging techniques and design pressure ratio are similar.

Paxson, Daniel E.↗
NASA NTRS1976

Composite propellant combustion modeling with a porous plate burner

A burner is designed to model on a large scale (millimeters) the complex vapor phase processes that take place on awkwardly small scales (about a hundred microns) in the combustion of practical AP/ composite propellants at conventional pressures. Binder vapor evolution is modeled with gaseous fuel (ethane in the experiments reported) flow through a porous plate and oxidizer vapor with the flow of a gaseous oxidizer (air and enriched air in the experiments reported) through discrete holes in the porous plate. Measured flame standoff distance and surface temperature variations are consistent with theoretical predictions at atmospheric pressure. Data obtained at several atmospheres are discussed in the light of the current theories of propellant burning that place varying emphasis on the roles of chemical kinetics and fluid dynamic diffusion/mixing in the vapor phase of a burning composite propellant. The potential and proposed future applications of the burner are indicated.

Kumar, R. N.↗
NASA NTRS1997

Turbulent Premixed Flame Propagation in Microgravity

A facility in which turbulent Couette flow could be generated in a microgravity environment was designed and built. To fit into the NASA Lewis drop tower the device had to be very compact. This means that edge effects and flow re-circulation were expected to affect the flow. The flow was thoroughly investigated using LDV and was found to be largely two dimensional away from the edges with constant turbulence intensities in the core. Slight flow asymmetries are introduced by the non symmetric re-circulation of the fluid outside the test region. Belt flutter problems were remedied by adding a pair of guide plates to the belt. In general, the flow field was found to be quite similar to previously investigated Couette flows. However, turbulence levels and associated shear stresses were higher. This is probably due to the confined re-circulation zone reintroducing turbulence into the test section. An estimate of the length scales in the flow showed that the measurements were able to resolve nearly all the length scales of interest. Using a new LES method for subgrid combustion it has been demonstrated that the new procedure is computational feasible even on workstation type environment. It is found that this model is capable of capturing the propagation of the premixed names by resolving the flame in the LES grid within 2-3 grid points. In contrast, conventional LES results in numerical smearing of the flame and completely inaccurate estimate of the turbulent propagation speed. Preliminary study suggests that there is observable effect of buoyancy in the 1g environment suggesting the need for microgravity experiments of the upcoming experimental combustion studies. With the cold flow properties characterized, an identical hot flow facility is under construction. It is assumed that the turbulence properties ahead of the flame in this new device will closely match the results obtained here. This is required since the hot facility will not enable LDV measurements. The reacting flow facility is also being constructed with planned drop tower experiments in mind. Therefore, issues related to safety and structural integrity are being take into account. Further development of the numerical model will also be carried out to include finite-rate kinetics for representative premixed cases. More detail analysis of the flame structure and propagation nature will be investigated. Simulations will also be compared to the flame properties observed in the experiments.

Menon, S.↗
NASA NTRS2009

Wind-US Code Physical Modeling Improvements to Complement Hypersonic Testing and Evaluation

This report gives an overview of physical modeling enhancements to the Wind-US flow solver which were made to improve the capabilities for simulation of hypersonic flows and the reliability of computations to complement hypersonic testing. The improvements include advanced turbulence models, a bypass transition model, a conjugate (or closely coupled to vehicle structure) conduction-convection heat transfer capability, and an upgraded high-speed combustion solver. A Mach 5 shock-wave boundary layer interaction problem is used to investigate the benefits of k- s and k-w based explicit algebraic stress turbulence models relative to linear two-equation models. The bypass transition model is validated using data from experiments for incompressible boundary layers and a Mach 7.9 cone flow. The conjugate heat transfer method is validated for a test case involving reacting H2-O2 rocket exhaust over cooled calorimeter panels. A dual-mode scramjet configuration is investigated using both a simplified 1-step kinetics mechanism and an 8-step mechanism. Additionally, variations in the turbulent Prandtl and Schmidt numbers are considered for this scramjet configuration.

Georgiadis, Nicholas J.↗
DOE OSTI2025

CRN Modeling of Ammonia RQL Combustion using a Partially-Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However low flammability and a propensity for high nitrogen oxide (NOx) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NOx emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel, while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NOx concentrations below that of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing non-adiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and post-flame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors’ past work, this arrangement has demonstrated the significant sensitivity of exit NOx to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and non-premixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NOx emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have shown to create conditions which are conducive to increased N2O formation in the lean stage. The results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia..

advanced gas turbines↗
NASA NTRS2000

A Model For The Burning of Teflon

Teflon has been identified as suitable material for use in oxygen-enriched atmospheres because of its low specific enthalpy of combustion that is, less than 1500 cal/gram. However, once ignited, Teflon burns in a heterogeneous reaction until total consumption or depletion of oxygen occurs. A model is developed for the burning of Teflon in pure, high-pressure oxygen, 3.4 to 68.9 MPa. The Teflon polymer chain dissociates via monomer units due to pyrolysis. These monomer units diffuse to the surface due to free convection. The model consists of coupled mass and energy balances for the polymer and an energy balance for the free convection of product gases. The model is used to obtain appropriate kinetic parameters for the dissociation and surface reactions. The model is validated against experimental measurements of Teflon discs supported on aluminum rods (2216) at oxygen pressures of 3.4, 6.8, 34.4, and 68.9 MPa. The model simulates the temperature distribution in the product gases above burning Teflon.

Waller, Jess M.↗
NASA NTRS1990

Carbon deposition model for oxygen-hydrocarbon combustion. Task 6: Data analysis and formulation of an empirical model

The formation and deposition of carbon (soot) was studied in the Carbon Deposition Model for Oxygen-Hydrocarbon Combustion Program. An empirical, 1-D model for predicting soot formation and deposition in LO2/hydrocarbon gas generators/preburners was derived. The experimental data required to anchor the model were identified and a test program to obtain the data was defined. In support of the model development, cold flow mixing experiments using a high injection density injector were performed. The purpose of this investigation was to advance the state-of-the-art in LO2/hydrocarbon gas generator design by developing a reliable engineering model of gas generator operation. The model was formulated to account for the influences of fluid dynamics, chemical kinetics, and gas generator hardware design on soot formation and deposition.

Makel, Darby B.↗
DOE OSTI2021

Oxidation Behavior of Candidate NiCr Alloys for Engine Exhaust Valves: Part I—Effect of Minor Alloying Elements

Oxidation-induced degradation of structural materials employed as exhaust valves within internal combustion engines (ICEs) will be a relevant life-limiting mechanism, in addition to creep and mechanical fatigue, due to ever-increasing severity of operating temperatures and pressures. Ni–Cr-based alloys, which form external chromia-based scales at the relevant operating temperatures are being considered as suitable candidate materials. Thermal cycling of these alloys in water vapor-containing atmospheres, such as those present during hydrocarbon fuel combustion within ICEs, can considerably influence their oxidation behavior. Here, the role of typical alloying additions such as Mn, Si, Al and Ti on the cyclic oxidation behavior of model NiCr–X (X = Mn,Si,Al,Ti) alloys exposed in dry air and air + 10% H 2 O at 800°C and 950°C was investigated. Combined additions of Mn and Si reduced scaling rates compared to binary Ni–22Cr alloys. The presence of water vapor possibly suppressed formation of NiMnCr spinel and thereby the Cr depletion in the alloy. Combined Al and Ti additions mainly resulted in accelerated oxidation kinetics due to the Ti doping of chromia scales. More porous external scales were observed in water vapor leading to a much deeper depth of nitridation in the Ni–22Cr–Al–Ti alloys.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Interactive computer modeling of combustion chemistry and coalescence-dispersion modeling of turbulent combustion

An interactive computer code for simulation of a high-intensity turbulent combustor as a single point inhomogeneous stirred reactor was developed from an existing batch processing computer code CDPSR. The interactive CDPSR code was used as a guide for interpretation and direction of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser diagnostic techniques to experimentally determine the appropriate mixing frequency, and for validation of CDPSR as a mixing-chemistry model for a laboratory jet-stirred reactor. The coalescence-dispersion model for finite rate mixing was incorporated into an existing interactive code AVCO-MARK I, to enable simulation of a combustor as a modular array of stirred flow and plug flow elements, each having a prescribed finite mixing frequency, or axial distribution of mixing frequency, as appropriate. Further increase the speed and reliability of the batch kinetics integrator code CREKID was increased by rewriting in vectorized form for execution on a vector or parallel processor, and by incorporating numerical techniques which enhance execution speed by permitting specification of a very low accuracy tolerance.

Pratt, D. T.↗
DOE OSTI2022

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗
NASA NTRS2013

Advanced Materials and Manufacturing for Low-Cost, High-Performance Liquid Rocket Combustion Chambers

A document describes the low-cost manufacturing of C103 niobium alloy combustion chambers, and the use of a high-temperature, oxidation-resistant coating that is superior to the standard silicide coating. The manufacturing process involved low-temperature spray deposition of C103 on removable plastic mandrels produced by rapid prototyping. Thin, vapor-deposited platinum-indium coatings were shown to substantially improve oxidation resistance relative to the standard silicide coating. Development of different low-cost plastic thrust chamber mandrel materials and prototyping processes (selective laser sintering and stereolithography) yielded mandrels with good dimensional accuracy (within a couple of mils) for this stage of development. The feasibility of using the kinetic metallization cold-spray process for fabrication of free-standing C1O3 thrusters on removable plastic mandrels was also demonstrated. The ambient and elevated temperature mechanical properties of the material were shown to be reasonably good relative to conventionally processed C103, but the greatest potential benefit is that coldsprayed chambers require minimal post-process machining, resulting in substantially lower machining and material costs. The platinum-iridium coating was shown to provide greatly increased oxidation resistance over the silicide when evaluated through oxyacetylene torch testing to as high as 300 F (= 150 C). The iridium component minimizes reaction with the niobium alloy chamber at high temperatures, and provides the high-temperature oxidation resistance needed at the throat.

Williams, Brian E.↗
DOE OSTI2023

Designing Electrolytes With Controlled Solvation Structure for Fast‐Charging Lithium‐Ion Batteries

Recharging battery-powered electric vehicles (EVs) in a similar timeframe as those used for refueling gas-powered internal combustion vehicles is highly desirable for rapid penetration of the EV market. It is well known that the electrolyte in a battery plays a critical role in fast-charging capability of the battery because it determines the rate of ion transport together with its derived electrode/electrolyte interphases on both cathode and anode of the battery. In this study, the effects of contents of salt, coordinating solvent, and noncoordinating diluent on salt dissociation degree and electrolyte ionic conductivity are investigated, and a controlled solvation structure electrolyte is developed to improve the lithium ion mobility and conductivity in the electrolyte and to enhance the kinetics and stability of the electrode/electrolyte interphases in the battery. This electrolyte enables fast-charging capability of high energy density lithium-ion batteries (LIBs) at up to 5 C rate (12-min charging), which significantly outperforms the state-of-the-art electrolyte. The controlled solvation structure sheds light on the future electrolyte design for fast-charging LIBs.

25 ENERGY STORAGE↗
DOE OSTI2023

Effect of water droplet growth dynamics on electrode current in fuel-cell catalyst layers

Fuel cells are a promising next-generation energy-conversion technology designed to replace internal combustion engines in transportation applications. However, much work remains to optimize them. Operation at high humidities causes liquid water droplet formation on Pt catalyst particles during oxygen reduction, potentially impeding reactant arrival to the reactive electrode. In this work, four different cases of water droplet growth in fuel-cell catalyst layers are considered: pinned or advancing droplets on a bare Pt surface, advancing droplets on a Nafion film, and water-layer growth in carbon nanopores. Transient drop growth is captured with a combination of mass, species mass, and momentum balances, and the subsequent limiting current is determined via oxygen diffusion and Tafel kinetics. Further, water droplets are found not to be mass-transfer limiting due to the relatively large liquid-gas area compared to the Pt nanoparticle. Mass-transfer-limited behavior is calculated in carbon nanopores.

25 ENERGY STORAGE↗
DOE OSTI2020

PAH formation from jet stirred reactor pyrolysis of gasoline surrogates

Soot particles and their precursor polycyclic aromatic hydrocarbon (PAH) species, formed during combustion, are responsible for particulate emissions in gasoline direct injection (GDI) engines. To better understand the effects of fuel composition on formation of soot in GDI engines, the pyrolysis of several gasoline surrogates was studied in a jet-stirred reactor across a broad temperature range at atmospheric pressure and 1 s residence time. Fuel and intermediate species, including aromatics up to naphthalene, were measured using gas chromatography (GC). PAH concentrations from pyrolysis of surrogate fuels were compared to gain insight into the effects of fuel composition on PAH formation. In addition, synergistic effects were observed in pyrolysis experiments of binary blends. A detailed kinetic model, recently developed at Lawrence Livermore National Laboratory (LLNL), successfully captured the effects of blending and the concentration of major PAHs. Major reaction pathways are discussed, as well as the role of important intermediate species, such as acetylene, and resonantly stabilized radicals such as allyl, propargyl, cyclopentadienyl, and benzyl in the formation of PAH.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2024

Computational Analysis of the Effect of Structured Packing Design on Absorption Column Hydrodynamics for Post-Combustion Carbon Capture Applications

Solvent based post-combustion carbon capture technologies have a potential for reducing carbon emissions from fossil-fuel-fired power plants and industrial sources where CO2 emissions are inherently harder to mitigate, such as steel or cement industries. A prominent technology to achieve this is by retrofitting absorption columns to the existing infrastructure. While these systems have been among the less costly alternatives for carbon capture, they still impose a considerable energy penalty to the operation of power plants or industrial facilities. The optimization of CO2 capture rate in solvent-based absorption process is complex as it depends on several factors including CO2 solubility, solvent reaction kinetics and temperature effects on the solubility, reaction rates, surface tension, and thermophysical properties of the solvent and the flue gas. The overall heat and mass transfer also depends on the hydrodynamics, which in turn, is affected by the packing geometry. In the current work, we systematically quantify the effects the design of the structured packing has on the column hydrodynamics, by performing detailed CFD simulations for different geometrical configurations and operating conditions. We then obtain relationships between the key hydrodynamic metrics, such as liquid holdup, interfacial area, wetted area, and pressure drop to the parameters defining the packing geometries and identify new more effective packing designs for the given operating conditions.

Shah, Yash Girish↗