Search NASA⌕ Search

SEARCH · Search NASA

Results for “cold pools”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Forest Breeze–Cold Pool Interactions Drive Convective Organization over Heterogeneous Vegetation

Heterogeneous landscapes can influence the development of convection through the generation of thermally driven mesoscale circulations. To assess the impacts of these circulations and their interaction with sea breezes, we simulated convection in an idealized coastal environment using the Regional Atmospheric Modeling System (RAMS). We compared simulations with striped patterns of surface vegetation to those of uniform vegetation to identify the importance of vegetation heterogeneity in impacting convective development. Under dry soil conditions representative of those during the Tracking Aerosol Convection Interactions Experiment (TRACER) and Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE) campaigns in June 2022, we found that these vegetation-induced circulations, referred to in the literature as “forest breezes,” are more important than the sea breeze in determining the location of convection initiation. Convection and precipitation are also found to be favored over forests and suppressed over pasture and suburban landscapes as a result of greater surface sensible heat flux over the forest. Our findings also indicate that forest breezes are important for initiating convection along the boundaries of the forest, but that cold pools may play a key role in propagating the forest breezes toward the center of the forest stripe. In our simulations, the collisions of these breezes in the center of the forest stripe lead to uplift and strong convection there; however, a different width of the forest stripe would alter when the forest breezes collide or whether they collide at all. The presence of these cold pools may therefore impact the “ideal stripe width,” the width of each vegetation stripe which maximizes domain-wide precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of Monin‐Obukhov Similarity Theory Wind Profiles in Convective Storm Environments and Cold Pools at the ARM Southern Great Plains Atmospheric Observatory

Monin-Obukhov similarity theory (MOST) is widely used in numerical weather prediction to model surface fluxes of momentum, heat, and water vapor. However, MOST is based on assumptions of steady state and horizontally homogeneous turbulence that can lead to prediction errors in and around convective storms. To understand the nature of these errors, we used wind and eddy covariance flux measurements from the Atmospheric Radiation Measurement Southern Great Plains Atmospheric Observatory to evaluate MOST wind profiles in fair-weather and convective storm environments, specifically those of mesoscale convective systems and ordinary thunderstorms. MOST wind profiles agreed well with observations in fair-weather cases, while in convective storm cases the theory systematically overestimated shear in cold pools after gust front passage. Surface layer stability was found to be important in assessing MOST within convective storm environments. The overestimation of wind shear in cold pools suggests the role of non-local fluxes in transferring momentum downward. We discuss reasons for differences and agreement with past studies, and conclude with recommendations to improve prediction of surface winds and fluxes in convective storm simulations.

58 GEOSCIENCES↗
DOE OSTI2025

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗
DOE OSTI2025

Investigation of the Sensitivity of Tropical Cyclogenesis to Aerosol Intervention

Abstract As risks from tropical cyclones (TCs) are fueled by climate change escalation, there is an urgent need for transformational solutions to complement traditional approaches. Seeding TCs using aerosols can be a promising method to reduce cyclone intensity, supported by theoretical understanding of the microphysical effects of aerosols on TC clouds. The ideal time to intervene effectively in TCs is likely during their initial stage, before TC wind speeds reach their peak. However, studies exploring potential aerosol effects on TC formation remain scarce. This study investigates how a TC embryo responds to the addition of aerosols of varying sizes using the Weather Research & Forecasting (WRF) model coupled with a spectral‐bin microphysics model. We found that aerosols of different sizes and concentrations distinctively affect the pre‐TC vortex's microstructure and dynamics. Fine and ultrafine aerosols enhance the latent heat of condensation, freezing, deposition, and riming, initially intensifying the vortex. However, this results in enhancement of the cold pool, thereby reducing inflow and surface fluxes, subsequently weakening the vortex. Coarse aerosols produce the opposite effect to that of fine and ultrafine aerosols. Coarse aerosols lead to a slower initial acceleration owing to enhanced warm rain. However, the resulting weaker cold pool is insufficient to effectively reduce the strength of the vortex at the later stage. This study provides critical insights into how aerosols of varying sizes and concentrations modulate the energy cascade and impact the evolution of a TC embryo, laying the groundwork for further research on TC risk management through aerosol intervention.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Urban Impacts on Convective Squall Lines over Chicago in the Warm Season—Part I: Observations of Multi-Scale Convective Evolution

In this study, our aim is to diagnose how two quasi-linear convective systems (QLCS) are organized so one can determine the possible role of the city of Chicago, IL, USA, in modifying convective precipitation systems. In this Part I of a two-part study, we employ large-scale analyses, radiosonde soundings, surface observations, and Doppler radar data to diagnose the precursor atmospheric circulations that organize the evolution of two mesoscale convective systems and compare those circulations to radar and precipitation. Several multi-scale processes are found that organize and modify convection over the Chicago metroplex. Two sequential quasi-linear convective systems (QLCS #1 and #2) were organized that propagated over Chicago, IL, USA, during an eight-hour period on 5–6 July 2018. The first squall line (QLCS #1) built from the southwest to the northeast while strengthening as it propagated over the city, and the second (QLCS #2) propagated southeastwards and weakened as it passed over the city in association with a polar cold front. The weak upper-level divergence associated with a diffluent flow poleward of an expansive ridge built over and strengthened a low-level trough and confluence zone, triggering QLCS #1. Convective downdrafts from QLCS #1 produced a cold pool that interacted with multiple confluent low-level jets surrounding and focused on the metroplex urban heat island, thus advecting the convection poleward over the metroplex. The heaviest precipitation occurred just south-southeast of Midway Airport, Chicago. Subsequently, a polar cold front propagated into the metroplex, which triggered QLCS #2. However, the descending air above it under the polar jet and residual cold pool from QLCS #1 rapidly dissipated the cold frontal convection. This represents a case study where very active convection built over the metroplex and was likely modified by it, as evidenced in numerical simulations to be described in Part II.

Kaplan, Michael L.↗
DOE OSTI2025

Interactions Between Clouds and Wind-Driven Surface Heat Exchanges over Land

Earth system model experiments show that increasing horizontal resolution fundamentally alters the simulated soil-moisture-precipitation feedback. Kilometer-scale simulations often produce weaker or even negative feedback compared to coarse-resolution models. A key difference of kilometer-scale models is that they resolve mesoscale secondary circulations, including boundary layer horizontal rolls and cellular structures, in addition to cold pools and downdrafts associated with convective precipitation. However, because the relevant processes occur on yet-smaller scales, these circulations are often poorly resolved. This project demonstrated that boundary layer secondary circulations significantly affect surface heat exchanges and wind gusts, and that current model parameterizations can misrepresent these processes at kilometer-scale resolution. Using DOE Atmospheric Radiation Measurement (ARM) observations and targeted experiments with the DOE Energy Exascale Earth System Model (E3SM), we identified physically unrealistic wind gust and surface flux responses to secondary circulations, diagnosed a systematic overestimation of wind shear in convective cold pools, and uncovered a multivariate relationship between land surface fluxes and the scales of updrafts that form shallow cumulus clouds. These findings provide observation-based recommendations for improving parameterizations of surface fluxes and wind gusts in high-resolution Earth system models, thereby reducing uncertainty in convective storm prediction and land-atmosphere feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Convectively Induced Secondary Circulations and Wind‐Driven Heat Fluxes in the Surface Energy Balance Over Land

Increased resolution has enabled kilometer-scale weather and climate models to partially resolve secondary circulations, including horizontal convective rolls (HCRs) and cold pool gust fronts. Although these circulations are ubiquitous in convective boundary layers over land, their impacts on the surface energy balance are largely unknown. Doppler lidar and surface observations were combined with DOE E3SM land model experiments, revealing increased surface winds (5 m/s) and heat fluxes (50 W/m 2 ) in convergent branches of HCRs. Larger wind-driven flux responses (up to 150 W/m 2 ) were found along gust fronts. Surface energy balance shifts to accommodate wind-driven fluxes, reducing ground heat conduction and longwave cooling. Our findings from the US Southern Great Plains are broadly relevant to modeling convective boundary layers. In particular, widely used subgrid wind gust parameterizations were found to be physically inconsistent with resolved secondary circulations and could worsen climate prediction biases at kilometer-scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Factors Controlling Precipitation Onset and Maintenance: Inferences From Large Eddy Simulations of Two Afternoon Deep Convective Regimes Over Amazon

To investigate the environmental factors controlling the onset and maintenance of afternoon precipitation over tropical rainforests, two contrasting cases are created for large eddy simulations (LES) using GoAmazon observations: one features a single-pulse rain dissipating quickly in earlier afternoon, while the other shows double pulses lasting until evening. For these two specific cases under consideration, LES confirms that early-morning relative humidity dominates afternoon rain patterns. The attributed impacts are distinct: moisture determines the rain onset timing while temperature affects the peak rain intensity. Single-pulse day observes one round of strong precipitation and cold pools, which further suppresses convection. On double-pulse day, the first precipitation peak results from an intermediate development of congestus clouds, whose detrainment leads to a gradual moistening of the lower-to-middle troposphere. This behavior of convection favors a second pulse of stronger precipitation with more convective organization, whose development decouples from surface fluxes and sustains until evening.

Amazon↗
DOE OSTI2025

Investigating the Relative Roles of INPs and CCN in a Simulated Thunderstorm Using a New Immersion Freezing Algorithm

Microphysical processes in deep convective clouds are sensitive to the number concentrations of cloud condensation nuclei (CCN) and ice nucleating particles (INPs), but the effects of INPs are less studied. Modeling studies investigating the effects of INPs and/or CCN on deep convection typically retain a volume‐dependent raindrop freezing relation. The resulting neglect of aerosol accumulation in raindrops via drop collisions has likely produced unrealistic storm responses to INPs in past studies. To address this deficiency, a new immersion freezing algorithm was developed and embedded in a bulk microphysics scheme that freezes both cloud drops and raindrops using the same immersion freezing INP (IF‐INP) activity spectrum based on measurements. Multiple idealized simulations of a single case of deep convection observed during the Clouds, Aerosols, and Complex Terrain Interactions (CACTI) field campaign were conducted, with microphysical differences produced by independently altering IF‐INP temperature dependencies and CCN number concentrations from their observed values. Surface precipitation in all simulations resulted almost exclusively from riming graupel that melted upon descending to the surface. Rainfall and cold pools were substantially and systematically weakened with increased CCN due to decreased graupel riming rates but were relatively insensitive to variations in the magnitude and slope of IF‐INP spectra due to compensating depletion of supercooled liquid water. These compensating processes were a consequence of the accumulation of IF‐INPs in raindrops, encouraging caution in studying IF‐INP effects upon thunderstorms using traditional volume‐dependent drop freezing relationships.

CCN↗
DOE OSTI2024

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2024

Effect of cold-cap-bottom shear stress on primary foam stability

Here, we analyze the heat transfer between the melt pool and the layer of reacting materials called cold cap, which is the determining factor for the glass production rate in Joule-heated air-bubbled nuclear waste glass melters. We consider the effects of bubbling rate, melt viscosity, melt thermal conductivity, conversion kinetics, and shear rate at the cold-cap bottom, and develop simple relationships for the cold-cap bottom temperature and glass production rate. We show that the melt convection enhances the heat transfer to the cold cap by suppressing the thermal boundary layer and primary foam, while the conversion kinetics opposes these effects. Additionally, we mention how dissolving silica particles affects the primary foam stability and discuss the limitations of the presented approach for vigorously foaming feeds.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Compositional Variation Tests on DuraMelter 100 with LAW Sub-Envelope A1 Feed (LAWA44 Glass) in Support of the LAW Pilot Melter, VSL-02R62N0-4, Rev. 0 (Nov 2002)

The two tests described in this report utilized blended feed (glass formers plus waste simulant) prepared by Optima Chemicals according to VSL specifications generating about 1. 7 metric tons of glass. Sugar was added (at VSL) to the nominal feed at a "sugar ratio" of 0.5 for each of the two variation tests1; however, since the sugar addition is assumed to be "blind" to the variations, the actual sugar ratios were 0.57 and 0.44. The DMl00-WV melter was used in order to provide a direct comparison with the LAW tests previously conducted on the same melter. Two 100-hour melter tests were conducted: one with a 15% deficiency in simulant and one with 15% excess in simulant. Key operating parameters including cold cap coverage, feed rate, and glass pool temperature were held constant to investigate the effects of the glass compositional changes on processing characteristics (including salt formation) and the product glass. The bubbling rate was adjusted to provide the desired glass production rate with a near complete cold cap (90-100% of melt surface covered with feed). Quantitative measurements of glass production rates, melter operating conditions (temperatures, pressures, power, flows, etc.), and off-gas characteristics (NOx, SO 2, CO, particulate load and composition, and acid gases) were made for each test.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗
DOE OSTI2025

Investigation of fuel film formation and soot emissions in a GDI engine during cold-start with Split-injection strategies

This study investigates the impact of fuel-film formation on engine-out soot emissions in a gasoline direct injection (GDI) engine under cold-start conditions. Split-injection strategies were applied by varying the number of injections, injection duration, and total fuel quantity to affect wall wetting and control the average in-cylinder equivalence ratio. A combined experimental and numerical approach was employed to analyze fuel-film deposition, combustion efficiency, and engine-out soot and unburnt hydrocarbons (UHC) emissions. In particular, fuel film distribution estimated by means of non-reacting, 3-D, computational fluid dynamics (CFD) simulations, together with experimentally measured soot data, were used to investigate fuel film formation and its role in soot generation. In the experiments, a skip-fired engine control strategy was applied to mimic the transient nature of engine cold-start operation. The results indicate that, under the same total number of injections, increasing the average in-cylinder equivalence ratio through longer injection durations improves combustion stability, as indicated by the decrease of the coefficient of variation of IMEP n (Net Indicated Mean Effective Pressure) from 6.0 % to 0.6 %. However, this strategy leads to higher soot emissions, which increased by nearly an order of magnitude, primarily due to enhanced wall-film formation. In contrast, increasing the number of injections while maintaining a constant equivalence ratio significantly impacts fuel-film deposition and, consequently, soot emissions, with a fivefold reduction of the measured engine-out soot, decreasing from 3.5 mg to 0.7 mg. A soot-film correlation was developed and achieved a high coefficient of determination (r 2 = 0.95) and was further extended to account for spark timing effects. These findings confirm the effectiveness of split-injection for avoiding wall film formation and soot emissions, and the critical role of fuel film in soot generation, supporting the hypothesis that pool fires play a crucial role in contributing to soot formation under these cold-start conditions. In conclusion, the study also indicates the value of a predictive soot-film correlation for developing cold-start emission control strategies.

Engine Cold-start↗