Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PISTONS”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Unraveling the role of EGR olefins at advanced combustion conditions in the presence of nitric oxide: Ethylene, propene and isobutene

We report the role of EGR (exhaust gas recirculation) olefinic constituents at advanced combustion conditions in the presence of nitric oxide is unraveled in this study through experimental and modeling efforts using a twin-piston rapid compression machine operating at a stochiometric fuel loading with 20% EGR by mass, pressures of 20 and 40 bar, and temperatures from 680 to 950 K. Five different levels of olefin addition, focusing on ethylene, propene and isobutene, with a fixed amount of NO at 70 ppm are doped into test mixtures of PACE-20, a multi-component gasoline surrogate, where olefin addition effects are characterized through changes in ignition times and heat release rates. Experiments indicate that all three EGR olefins inhibit autoignition reactivity and low-temperature heat release at Tc < 850 K, with isobutene exhibiting the greatest impact, while at Tc > 850 K, low ethylene and propene additions promote reactivity. A recently updated chemical kinetic model, with detailed gasoline/NO x interacting and olefin/NO x interacting chemistry incorporated, is adopted to simulate the experiments. Simulation results are somewhat inconsistent with the experiments, where the model captures the inhibiting effects of all olefins on first-stage ignition reactivity, while consistently predicting a promoting effect on main ignition reactivity. Sensitivity and rate of production analyses reveal that adding olefins greatly alters the role of the consuming pathways for the olefins and their primary derivatives at the initial stage of the oxidation process, particularly with the presence of NO, where the olefins and their derivatives interact with both NO x species such as NO 2 and other species such as OH and HO 2 . The olefin/NO x interactions are particularly pronounced with propene and isobutene addition, where these lead to increased ignition reactivity by facilitating NO production and are mostly responsible for the disagreement between the model and experiments. Further investigations of olefin interacting chemistry, particularly those with NOx species, are needed for chemistry models to accurately predict the complicated effects of EGR.

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

Energy dissipation in polymer-bonded explosives with various levels of constituent plasticity and internal friction

The ignition of energetic materials (EM) under dynamic loading is mainly controlled by localized temperature spikes known as hotspots. Hotspots occur due to several dissipation mechanisms, including viscoplasticity, viscoelasticity, and internal friction along crack surfaces. In this work, to analyze the contributions of these mechanisms, we quantify the ignition probability, energy dissipation, damage evolution, and hotspot characteristics of polymer-bonded explosives (PBXs) with various levels of constituent plasticity of the energetic phase and internal crack face friction. Using PBX9501 consisting of HMX (Octahydro-1,3,5,7-Tetranitro-1,2,3,5-Tetrazocine) and Estane as a reference material, we analyze variants of this material with several values of the yield stress of the energetic phase and coefficients of internal crack face friction, while other parameters are kept unchanged. The impact loading involves piston velocities between 200 and 1200 m/s. The analysis uses a Lagrangian cohesive finite element framework that explicitly accounts for finite-strain elastic-viscoplastic deformation of the grains, viscoelastic deformation of the binder, arbitrary crack initiation and propagation in the grains and the binder, debonding between the grains and the binder, contact between internal surfaces, friction and frictional heating along internal surfaces, heat generation resulting from inelastic bulk deformation, and heat conduction. To determine the ignition status of the material or “go” or “no-go” state, we use a criterion based on a criticality threshold obtained from chemical kinetics calculations. For PBX with various levels of HMX plasticity and friction, the probability of ignition, the evolution of dissipation caused by plasticity and friction, the density of cracks, and the locations of cracks are quantified. Results show that samples with higher levels of constituent plasticity (lower yield strengths) or lower levels of internal friction are less likely to ignite. The relative importance of plasticity and friction depends on load intensity, with frictional heating decreasing as load intensity increases. Although the overall viscoplastic heating outweighs the overall frictional heating, friction plays a very important role in hotspot development at all load intensities analyzed, owing to the fact that frictional heating is more localized than viscoplastic heating. Finally, the predicted thresholds and ignition probabilities are expressed in a load intensity-load duration relation for PBX with different constituent properties.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Superior wear resistance of diamond and DLC coatings

As the hardest known material, diamond and its coatings continue to generate significant attention for stringent applications involving extreme tribological conditions. Likewise, diamond-like carbon (DLC, especially the tetragonal amorphous carbon, ta-C) coatings have also maintained a high level interest for numerous industrial applications where efficiency, performance, and reliability are of great importance. Additionally, the strong covalent bonding or sp 3 -hybridizaiton in diamond and ta-C coatings assures high mechanical hardness, stiffness, chemical and thermal stability that make them well-suited for harsh tribological conditions involving high-speeds, loads, and temperatures. In particular, unique chemical and mechanical nature of diamond and ta-C surfaces plays an important role in their unusual friction and wear behaviors. As with all other tribomaterials, both diamond and ta-C coatings strongly interact with the chemical species in their surroundings during sliding and hence produce a chemically passive top surface layer which ultimately determines the extent of friction and wear. Thick micro-crystalline diamond films are most preferred for tooling applications, while thinner nano/ultranano-crysalline diamond films are well-suited for mechanical devices ranging from nano- (such as NEMS) to micro- (MEMS and AFM tips) as well as macro-scale devices including mechanical pump seals. The ta-C coatings have lately become indispensable for a variety of automotive applications and are used in very large volumes in tappets, piston pins, rings, and a variety of gears and bearings, especially in the Asian market. This paper is intended to provide a comprehensive overview of the recent developments in tribology of super-hard diamond and DLC (ta-C) films with a special emphasis on their friction and wear mechanisms that are key to their extraordinary tribological performance under harsh tribological conditions. Finally, based on the results of recent studies, the paper will also attempt to highlight what lies ahead for these films in tribology and other demanding industrial applications.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Composition of the core from gallium metal–silicate partitioning experiments

We present gallium concentration (normalized to CI chondrites) in the mantle is at the same level as that of lithophile elements with similar volatility, implying that there must be little to no gallium in Earth's core. Metal-silicate partitioning experiments, however, have shown that gallium is a moderately siderophile element and should be therefore depleted in the mantle by core formation. Moreover, gallium concentrations in the mantle (4 ppm) are too high to be only brought by the late veneer; and neither pressure, nor temperature, nor silicate composition has a large enough effect on gallium partitioning to make it lithophile. We therefore systematically investigated the effect of core composition (light element content) on the partitioning of gallium by carrying out metal–silicate partitioning experiments in a piston–cylinder press at 2 GPa between 1673 K and 2073 K. Four light elements (Si, O, S, C) were considered, and their effect was found to be sufficiently strong to make gallium lithophile. The partitioning of gallium was then modeled and parameterized as a function of pressure, temperature, redox and core composition. A continuous core formation model was used to track the evolution of gallium partitioning during core formation, for various magma ocean depths, geotherms, core light element contents, and magma ocean composition (redox) during accretion. The only model for which the final gallium concentration in the silicate Earth matched the observed value is the one involving a light-element rich core equilibrating in a FeO-rich deep magma ocean (>1300 km) with a final pressure of at least 50 GPa. More specifically, the incorporation of S and C in the core provided successful models only for concentrations that lie far beyond their allowable cosmochemical or geophysical limits, whereas realistic O and Si amounts (less than 5 wt.%) in the core provided successful models for magma oceans deeper that 1300 km. In conclusion, these results offer a strong argument for an O- and Si-rich core, formed in a deep terrestrial magma ocean, along with oxidizing conditions.

58 GEOSCIENCES↗
DOE OSTI2025

Co-optimization of fuel properties, combustion system geometry, and injection strategy for conventional diesel fuel

Here, studies have shown that fuel properties can impact an engine’s operation in several ways, including ignition delay, sooting tendency, mixture formation, and combustion temperature. In mixing-controlled compression ignition (MCCI) engines, the fuel system design and piston bowl geometry significantly affect combustion performance and emissions. Based on current information, it is difficult to draw conclusions about fuel property effects and sensitivities. The central fuel hypothesis approach used in the US Department of Energy Co-Optima program has worked well for spark ignition fuels: identifying critical fuel property ranges is sufficient to screen fuel blends that are expected to maximize efficiency and reduce pollutant emissions. However, for MCCI-relevant fuels, the information gained from past studies is not sufficient to build such a merit function or to allow for performing a similar screening of fuel blends. It is hypothesized that a co-optimization of a fuel’s physical and chemical properties, combustion system geometry, and injection strategy could leverage synergies between the effects of the fuel properties and geometries, resulting in improved performance over state-of-the-art. A machine learning–assisted unconstrained global optimization algorithm was used to explore a design space comprising 23 independent variables. The results show that physical property effects were minimal even for large variations in fuel properties, and the only interaction effect that was observed was the effect of varied fuel density parameters on fuel/air mixture formation. Nevertheless, these interactions were not sufficient in magnitude to significantly affect optimization results. Therefore, analysis of the results suggests that fuel physical properties cannot be leveraged in a co-optimization context to increase engine efficiency.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

Assessment of shallow subsea hydrocarbons as a proxy for leakage at offshore geologic CO 2 storage sites

We report this study is part of a multi-phase effort to identify and characterize offshore geological carbon storage (GCS) potential along the Texas Gulf Coast. Previous efforts acquired a high-resolution 3D seismic dataset (P-cable™) and interpreted a seismically discontinuous zone as shallow gas pockets (<100 m below seafloor) associated with a vertical chimney structure. Our approach was to measure hydrocarbon concentrations and stable carbon isotopes near the seafloor to assess if gas migrated vertically from seismic anomalies. Deep-sourced thermogenic hydrocarbons would indicate that the structure is transmissive which could indicate an unacceptable risk for GCS at the site. Alternatively, hydrocarbons formed in-situ from biogenic processes would not preclude transmissivity of the structure but would add information to the risk assessment. Gases were extracted from 23 piston core samples recovered between 2.56 and 3.50 m subseafloor. Light hydrocarbons and stable carbon isotopes of methane were used for source attribution. The result was that geochemical signatures were consistent with typical background values observed within the first few meters of subsea sediment and therefore did not indicate leakage from depth. From our analysis, we offer insights into the use of hydrocarbon molecular and isotopic compositions at the seafloor for signal attribution as part of environmental assessments at geological CO 2 storage sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An energy model of droplet impingement on an inclined wall under isothermal and non-isothermal environments

The observation of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of different droplet-wall impingement conditions is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process. These conditions are grouped as: a single diesel droplet impinging on the wall with different inclined angles (α); cold wall - heated droplet and heated wall - cold droplet impingement when inclined angle of the wall is 0°, respectively. The analytical model is built up based on the energy conservation that considers kinetic energy, gravitation energy, and surface energy before impingement, as well as viscous dissipation, gravitation energy, adhesion energy, deformation energy, and heat energy after impingement. The experimental work of diesel droplet impinging on an inclined wall is performed at a certain range of the Weber number (We of 33 to 420) with various inclined angles (α of 0° to 45°), while for inclined angle is 0°, droplet and wall temperature are varied from 25°C to 150°C to study the effects of the inclined angle and temperature on the temporal evolution of the post-impingement characteristics (i.e. droplet spreading length, dynamic contact angle). The analytical model is validated and evaluated at the aforementioned experimental operating points. The validated model can be employed to predict maximum spreading length of the droplet impinged on the wall. It is further utilized to determine the transition from capillary regime to kinetic regime, then to viscous regime at different inclined angle of the wall.

42 ENGINEERING↗
DOE OSTI2018

Rapid high density cryogenic pressure vessel filling to 345 bar with a liquid hydrogen pump

Cryogenic pressurized hydrogen (H 2 ) vessels promise maximum storage density with potential to enabling practical H 2 vehicles with maximum driving autonomy and minimum cost of ownership. This paper adds a more complete evaluation of the benefits of cryogenic vessel technology by establishing a methodology for evaluating fill density for any initial vessel thermodynamic state. This is accomplished by assessing 24 cryogenic pressure vessel fill experiments with a liquid hydrogen (LH 2 ) pump manufactured by Linde and installed at the Lawrence Livermore National Laboratory (Livermore, CA) campus. The LH 2 piston pump takes LH 2 from the station Dewar at near ambient pressure (3 bar) and very low temperature (24.6 K) and pressurizes it to the vessel pressure in two stages of compression, up to 875 bar, although the rating of the cryogenic vessel used for these experiments limits fill pressure to 345 bar. Experiments spanned initial vessel temperatures from ambient to 22 K, enabling pump testing over a broad range of conditions. Experimental results validate many of the virtues that make LH 2 pumping a promising technology for H 2 dispensing, both for filling cryogenic as well as ambient temperature vessels: high throughput (1.67 kg/min, 100 kg/h), unlimited capacity for back to back refueling, low electricity consumption at the station due to high density of LH 2 minimizing compression work, and highest fill density, up to 80 g/L (estimated) when dispensing cryogenic hydrogen at 700 + bar. Analysis of the experiments shows that fill density can be predicted with reasonable accuracy (±0.7 g/L) by assuming 10 kJ/kg K cryogenic vessel inlet entropy (pump delivery hose outlet entropy). Fill lines are finally generated on a H 2 phase diagram, producing charts that can be used for rapidly determining fill density for any vessel condition at the moment of fill initiation.

08 HYDROGEN↗
DOE OSTI2021

The simulation of 3D hypervelocity spallation using a hydrocode PAGOSA with FLIP+MPM

n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.

42 ENGINEERING↗
DOE OSTI2025

Crystallographic stability of the BCC phase during quenching of metastable beta titanium alloy Ti-5553 and comparison of structural criteria for the predictive capability of α” martensite

This work evaluates the compositional standard currently developed for Ti-5553 powder and explores the metastable β region for sensitivity to martensitic formation during rapid quenching from the melt. Ti-5553 is a beta-stabilized titanium alloy that is increasingly being used in additive manufacturing applications. A series of alloys within the Ti-5553 compositional space was processed using two-piston spat quenching to perform rapid solidification and quenching for each of these alloy compositions. Typical empirical models such as molybdenum equivalency are not able to fully separate the retained BCC and martensitic compositions. Further, the use of thermodynamic data to estimate the transformation energy needed to form martensite can differentiate the alloys and provide a metric to further develop compositional limits for metastable beta titanium alloy development.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Effect of thermally stable Cu- and Mg-rich aluminides on the high temperature strength of an AlSi12CuMgNi alloy

The internal architecture of an AlSi12CuMgNi piston alloy, revealed by synchrotron tomography, consists of three dimensional interconnected hybrid networks of Cu-rich aluminides, Mg-rich aluminides and eutectic/primary Si embedded in an α-Al matrix. The strength at room temperature and at 300°C is studied as a function of solution treatment time at 490°C and compared with results previously reported for an AlSi12Ni alloy. The addition of 1 wt% Cu and 1 wt% Mg to AlSi12CuMgNi increases the room temperature strength by precipitation hardening while the strength at 300°C is similar for both alloys in as-cast condition. The strength of AlSi12CuMgNi decreases with solution treatment time and stabilizes at 4 h solution treatment. The effect of solution treatment time on the strength of the AlSi12CuMgNi alloy is less pronounced than for the AlSi12Ni alloy both at room temperature and at 300°C. - Highlights: • The 3D microstructure of AlSi12CuMgNi is revealed by synchrotron tomography. • An imaging analysis procedure to segment phases with similar contrasts is presented. • 1 wt% Cu and Mg results in the formation of 3D networks of rigid phases. • AlSi12CuMgNi is stronger than AlSi12Ni owing to the stability of the 3D networks.

36 MATERIALS SCIENCE↗
DOE OSTI2015

An improved understanding of the natural resonances of moonpools contained within floating rigid-bodies: Theory and application to oscillating water column devices

The fundamental interactions between waves, a floating rigid-body, and a moonpool that is selectively open to atmosphere or enclosed to purposefully induce pressure fluctuations are investigated. The moonpool hydrodynamic characteristics and the hydrodynamic coupling to the rigid-body are derived implicitly through reciprocity relations on an array of field points. By modeling the free surface of the moonpool in this manner, an explicit hydrodynamic coupling term is included in the equations of motion. This coupling results in the migration of the moonpool's natural resonance frequency from the piston frequency to a new frequency when enclosed in a floating rigid-body. Two geometries that highlight distinct aspects of marine vessels and oscillating water column (OWC) renewable energy devices are analyzed to reveal the coupled natural resonance migration. The power performance of these two OWCs in regular waves is also investigated. The air chamber is enclosed and a three-dimensional, linear, frequency domain performance model that links the rigid-body to the moonpool through a linear resistive control strategy is detailed. Furthermore, an analytic expression for the optimal linear resistive control values in regular waves is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

High temperature materials for heavy duty diesel engines: Historical and future trends

We report that heavy duty (HD) vehicles are projected to be the largest fuel-use subsector in transportation, with current demand for diesel fuel projected to grow 30% by 2040. Historically, a primary strategy for increasing diesel engine efficiency has been to increase peak cylinder pressure (PCP). However, increasing PCP imparts greater mechanical and thermal loads on engine components and materials. In recent years, the material property limits for many components have been reached and further increases in PCP above ~20 MPa have been difficult, while still maintaining the necessary affordability and longevity of on-road HD diesel engines. This paper reviews the historical evolution and major metallurgical advancements of high temperature materials in HD on road diesel engines (10–15 L displacement) up to the current state of the art, focusing on materials in the engine block, cylinder heads, pistons, valves, and exhaust components. These components cover a wide range of material classes, including cast iron, ferritic steel, austenitic steel, titanium alloys, nickel based super-alloys, and high temperature coatings. Finally, the microstructural degradation and failure mechanisms of the materials associated with the complex mechanical and thermal loading during service are discussed and key areas for future materials research are suggested that overcome technical barriers.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A review of micromirror arrays

The aim of this article is to provide a review of micromirror array (MMA) technologies (2631 MMA research papers and patents were reviewed for this effort). The performance capabilities of 277 MMA designs from 49 companies and 23 academic research groups are categorized and compared. The designs are categorized according to (i) their array’s dimension (e.g., 2D arrays consisting of mirrors that cover a surface, 1D arrays consisting of mirrors in a row, and 0D arrays consisting of only a single mirror), (ii) the nature of the surface of their mirrors (e.g., continuous or discrete), (iii) what combination of tip, tilt, and/or piston degrees of freedom (DOFs) they achieve, and (iv) how they are actuated. Standardized performance metrics that can be systematically applied to every MMA design (e.g., mirror area, fill factor, pitch, range of motion, maximum acceleration, actuator energy density, and number of uncontrolled DOFs) are defined and plotted for existing designs to enable their fair comparison. Theoretical bounds on what is physically possible for MMAs to achieve are also derived and depicted in these plots to highlight the amount of performance improvement that remains to be achieved by future designs and guidelines are provided to aid in the development of these future designs.

42 ENGINEERING↗
DOE OSTI2018

A simplified CFD model for spectral radiative heat transfer in high-pressure hydrocarbon–air combustion systems

Detailed radiation modeling in piston engines has received relatively little attention to date. Recently, it is being revisited in light of current trends towards higher operating pressures and higher levels of exhaust-gas recirculation (EGR), both of which enhance molecular gas radiation. Advanced high-efficiency engines also are expected to function closer to the limits of stable operation, where even small perturbations to the energy balance can have a large influence on system behavior. Detailed radiation modeling using sophisticated tools like photon Monte Carlo/line-by-line (PMC/LBL) is computationally expensive. Here, guided by results from PMC/LBL, a simplified stepwise-gray spectral model in combination with a first-order spherical harmonics (P1 method) radiative transfer equation (RTE) solver is proposed and tested for engine-relevant conditions. Radiative emission, reabsorption and radiation reaching the walls are computed for a heavy-duty compression-ignition engine at part-load and full-load operating conditions with different levels of EGR and soot. The results are compared with those from PMC/LBL, P1/FSK (P1 with a full-spectrum k-distribution spectral model) and P1/Gray radiation models to assess the proposed model’s accuracy and computational cost. The results show that the proposed P1/StepwiseGray model can calculate reabsorption locally and globally with less than 10% error (with respect to PMC/LBL) at a small fraction of the computational cost of PMC/LBL (a factor of 30) and P1/FSK (a factor of 15). In contrast, error in computed reabsorption by the P1/Gray model is as high as 60%. Finally, it is expected that the simplified model should be broadly applicable to high-pressure hydrocarbon–air combustion systems, with or without soot.

42 ENGINEERING↗