Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster characterization”

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

Characterization of Precipitation-Induced Radon Progeny Deposition Events Using a City-Scale Sensor Network

Networks of radiation detectors provide a platform for real-time radioactive source detection and identification in urban environments. Detection algorithms in these systems must adapt to naturally-occurring changes in background, which requires well-characterized relationships between precipitation events and their corresponding radiological signature. Here, we present a quantitative and qualitative description of rain-induced radon progeny deposition events occurring in Chicago from September 2023 to February 2024. We measure ambient gamma radiation levels, precipitation rate, temperature, pressure, and relative humidity in a network of sensor nodes. For each identified precipitation period, we decompose spectra into static- and radon-associated components as defined by a non-negative matrix factorization (NMF) algorithm. We find a consistent power-law relationship between a precipitation-dependent peak of the radon progeny proxy (RPP) and the peak strength of the radon-associated NMF component for most precipitation events. We conduct a case study of a rainfall period with abnormally high levels of implied radon progeny concentration and describe its temporal and spatial evolution. We hypothesize that this phenomenon is due to the air mass path that intersects a uranium-rich region of Wyoming. Finally, we cluster precipitation events into three distinct categories. One category roughly corresponds to events with deep low-pressure systems and high relative radon concentration, while another is characteristic of light stratiform rain with slightly higher temperatures and intermediate relative radon concentration. The third category appears to contain weak-gradient or lake breeze convection showers with intermittent precipitation and low relative radon concentration. These findings suggest that radiological anomaly detection could be improved by training unique background models corresponding to each category of meteorological event.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Structural and functional characterization of IgG- and non-IgG-based T-cell-engaging bispecific antibodies

Bispecific T-cell-engaging antibodies are a growing class of therapeutics with numerous molecules being tested in clinical trials and, currently, seven of them have received market approval. They are structurally complex and function as adaptors to redirect the cytotoxicity of T cells to kill tumor cells. T-cell-engaging bispecific antibodies can be generally divided into two categories: IgG/IgG-like and non-IgG-like formats. Different formats may have different intrinsic potencies and physiochemical properties, and comprehensive studies are needed to gain a better understanding of how the differences in formats impact on structural and functional characteristics. In this study, we designed and generated bispecific T-cell-engaging antibodies with IgG-like (DVD-Ig) and non-IgG (BiTE) formats. Both target the same pair of antigens (EGFR and CD3) to minimize the possible influence of targets on functional characterization. We performed a side-by-side comparison to assess differences in the physiochemical and biological properties of these two bispecific T-cell-engaging antibodies using a variety of breast and ovarian cancer cell-based functional assays to delineate the structural–functional relationships and anti-tumor activities/potency. We found that the Fc portion of T-cell-engaging bispecific antibodies can significantly impact antigen binding activity, potency, and stability in addition to eliciting different mechanisms of action that contribute the killing of cancer cells.

cluster of differentiation 3 (CD3)↗

Analysis of the Ion Conversion Mechanisms in the Effluent of Atmospheric Pressure Plasma Jets in Ar with Admixtures of O 2 , H 2 O and Air

Ionic species in atmospheric pressure plasma jets (APPJs) play an important role in plasmasurface and plasma-liquid interactions, nonetheless, they have not received the same attention as their neutral reactive species counterparts. In this work, a molecular beam mass spectrometer (MBMS) was used to characterize the ion compositions in the effluent of an APPJ operating in ambient air for different feed gases including Ar + O 2 , Ar + air and Ar + H 2 O mixtures inspired by gas compositions used for biomedical applications. Changes in compositions of positive and negative ions as a function of nozzle-to-substrate distance along the plasma plume were analyzed and compared with a pseudo-1D plug flow model. Positive and negative ions were detected up to distances of 12 mm from the visible plasma plume tip. The measurements enable to follow the ion conversion pathways in the effluent of the APPJs as a function of distance from the nozzle. The trends in ion yield obtained from a pseudo-1D plug flow model showed generally a good agreement with the experimentally observed trends after addition of ionic reactions to the previously reported reaction set but also some distinctive differences were observed. The dominant positive ions in the far effluent are water ion clusters, the most stable ion for all gas mixtures investigated, while a large variety of negative ions was found for different gas mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncertainty in Synthetic Tropical Cyclone Hazard and Risk Estimates: Insights from RAFT, CHAZ, MIT, STORM, and CLIMADA

We synthesize five complementary tropical cyclone (TC) hazard frameworks—RAFT (physics-based machine learning), CHAZ and MIT (statistical–dynamical), STORM (fully statistical), and CLIMADA (observation-driven resampling)—to characterize uncertainty in wind-related TC metrics relevant to energy applications. All datasets and the IBTrACS observational record are harmonized to a common 6-hourly, 2.5° grid. We compare basin-wide and coastal properties using consistent definitions for TC frequency, mean and maximum intensity, 24-hour intensification, and 6-hour translation speed, and quantify agreement with Pearson r, RMSE, and Kling–Gupta efficiency (KGE) alongside resampling-based confidence intervals. CLIMADA is included for basin context but excluded from coastal skill scoring because it resamples historical IBTrACS; if supplied with projected future tracks from an external hazard model, CLIMADA can be used to simulate future TC scenarios. Results show robust, cross-model signals: (i) a corridor of activity from the tropical Atlantic through the Caribbean into the Bahamas and western subtropical Atlantic; (ii) a meridional dipole in 24-hour intensification (low-latitude strengthening, subtropical weakening); and (iii) a transition from slower tropical motion to faster midlatitude translation. Coastal winds (mean and maximum) consistently cluster from the eastern Gulf into the Bahamas–western Atlantic transition. The largest structural spread occurs in the amplitude and footprint of lifetime maximum intensity and, secondarily, in translation speed; intensification exhibits similar central behavior across frameworks with variability in extremes. Translation speed shows the most uniform coastal agreement. These findings provide a decision envelope for wind-focused risk screening and clarify where uncertainty should be carried forward; wind-only results represent a lower bound on total hazard, motivating integration of surge and rainfall modules and a companion, asset-level damage analysis.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Characterization of local parameters of inclined upward air–water bubbly two-phase flows in a round pipe

The present work seeks to develop a detailed experimental database and characterize the local effects of inclination on horizontal and upward bubbly air–water two-phase flows. Experimental conditions are selected to investigate bubbly flows at a variety of orientations while reducing the influence of flow regime transition. Detailed time-averaged data is collected using state-of-the-art four-sensor conductivity probes for measurements across the flow area for five flow conditions at five inclination orientations in a 25.4 mm inner diameter pipe. Using the newly established database, effects of inclination on bubbly flows are analyzed, as are changes in liquid and gas flow rate at different orientations. The observations include (1) bubble migration and agglomeration at the top pipe wall as orientation changes from vertical to horizontal; (2) reducing the angle from vertical reduces the local gas velocity, and correspondingly reduces the global slip from greater than unity at 90° to less than unity for 60° and lower; (3) changes in the bubble size distribution are primarily at the interface between the bubble cluster and the liquid region of the pipe, not coalescence within the bubble cluster. In conclusion, physical mechanisms for these observations are proposed, and the present trends at different orientations support previous experiments in horizontal two-phase flows.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Introducing the DREAMS Project: DaRk mattEr and Astrophysics with Machine Learning and Simulations

We introduce the DaRk mattEr and Astrophysics with Machine learning and Simulations (DREAMS) project, an innovative approach to understanding the astrophysical implications of alternative dark matter (DM) models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over DM physics, astrophysics, and cosmology in modeling a range of systems—from galaxy clusters to ultra-faint satellites. Such extensive simulation suites can provide adequate training sets for machine-learning-based analyses. This paper introduces two new cosmological hydrodynamical suites of warm dark matter (WDM), each comprising 1024 simulations generated using the arepo code. One suite consists of uniform-box simulations covering a ${(25\,{h}^{-1}\,{\rm{Mpc}})}^{3}$ volume, while the other consists of Milky Way zoom-ins with sufficient resolution to capture the properties of classical satellites. For each simulation, the WDM particle mass is varied along with the initial density field and several parameters controlling the strength of baryonic feedback within the IllustrisTNG model. We provide two examples, separately utilizing emulators and convolutional neural networks, to demonstrate how such simulation suites can be used to disentangle the effects of DM and baryonic physics on galactic properties. The DREAMS project can be extended further to include different DM models, galaxy formation physics, and astrophysical targets. In this way, it will provide an unparalleled opportunity to characterize uncertainties on predictions for small-scale observables, leading to robust predictions for testing the particle physics nature of DM on these scales.

cosmological parameters↗

Molecular interactions in short-chain perfluoroalkyl carboxylic acids and aqueous solutions

The presence of short-chain per- and polyfluoroalkyl substances in water poses a major health and environmental challenge. Here, we have performed high-energy small- and wide-angle X-ray scattering measurements on CF 3 [CF 2 ] n COOH (where n = 1, 2, 3 represents the chain length) and their aqueous solutions at 10% mole concentrations to characterize their molecular interactions at the atomic and nanometer length scales. The experimental wide-angle structure factors have been modelled using Empirical Potential Structural Refinement. The oxygen–oxygen partial X-ray pair distribution functions show that the coordination number between the hydroxyl oxygen on the acid and surrounding oxygen water molecules increases significantly with acid chain length, rising from 3.2 for n = 1 to 4.1 for n = 3. The small-angle scattering is dominated by a sharp, high-intensity peak at Q 1 ~ 0.2 Å -1 and a smaller peak at Q 2 = 1.2 Å -1 for n = 3, both of which decrease with decreasing chain length. The Q 2 peak is attributed to groups of adjacent non-bonded acid molecules, and Q 1 has contributions from both correlations between acid molecules and water–water interactions. In all cases, the models show nanoscale aggregation occurs in the form of denser channels of winding hydrogen-bonded chains, approximately 20 water molecules in length, surrounding clusters of acid molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗

Swift heavy ion irradiation effects in zirconium and hafnium carbides

The behavior of microcrystalline zirconium carbide (ZrC) and hafnium carbide (HfC) was studied under highly ionizing irradiation conditions at room temperature. The induced structural modifications were characterized via synchrotron-based X-ray diffraction experiments. Unit-cell expansion and buildup of microstrain were determined across a wide fluence range and linked to chemical compositions of the target material. The observed swelling resulting from irradiation with 198 MeV Xe ions in both carbide materials is characterized by two distinct mechanisms that operate within different fluence regimes. Unit-cell expansion initially proceeds by a direct-impact behavior that reaches saturation, followed at higher fluences by a second, linear swelling regime. The overall behavior, particularly the direct-impact regime, is similar for ZrC and HfC, with a more pronounced second defect accumulation process in HfC. Swelling in ZrC shows the same two distinct mechanisms upon irradiation with 198 MeV Xe ions and 946 MeV Au ions, but swelling induced by the lower-energy ions is greater across the entire fluence series. Accounting for the difference in energy deposition density between the two irradiation conditions reveals that the first swelling mechanism (direct-impact behavior) is likely related to the formation of more simple defects. In contrast, the second damaging mechanism at higher fluences (linear increase) cannot be fully explained by the induced energy density, and swelling remains somewhat higher for the low-velocity Xe irradiation. This may suggest that more complex defects and defect clusters are responsible for this swelling regime, with either their size and/or morphology modified at different energy densities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Detecting Dwarf Galaxies Around NGC 55

In order to improve our understanding of the satellite populations of low-mass galaxies, the DECam Local Volume Exploration (DELVE) – DEEP Survey is performing a search for satellites around the isolated, low-mass galaxy, NGC 55. As part of this search, we characterize our dwarf galaxy detection sensitivity by injecting artificial dwarf galaxies into our data set and measuring the fraction that we can recover. Additionally, we inject and recover artificial dwarfs to test how the blending of resolved stars affects the detection of dwarf galaxies. We inject 650 artificial dwarfs across 26 DELVE-DEEP coadded images of the NGC 55 halo and use the software tool, Source Extractor, to perform source recovery. A dwarf-search pipeline is then run on our dwarf-injected catalogs, utilizing the photutils function, find peaks, and detections are recorded as a function of dwarf size and luminosity. Our results demonstrate that we are most sensitive to large, bright dwarf galaxies, and our sensitivity begins to decrease as we approach DELVE-DEEP’s detection limit or blending causes significant reductions in recovered magnitude. We also find that our results are hindered by the limitations of star-galaxy separation applied to distant, faint sources, which is possibly a consequence of the blending of small star clusters within dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Tuning the Molecular Structure and Reaction Mechanism of Olefin Metathesis by Model Bilayered Supported MoO x /AlO x /SiO 2 Catalysts

The molecular structure and activity of supported MoO x olefin metathesis catalysts are heavily impacted by the choice of catalyst support. In this study, surface modification of the SiO 2 support with AlO x and selective anchoring of the MoO x on the surface AlO x sites were used to tune the structure, activation, and reactivity of the resulting surface MoO x sites. Extensive in situ molecular characterization, chemical probe studies, and density functional theory (DFT) calculations reveal that the enhanced activity of the supported MoO x /AlO x /SiO 2 catalyst over the MoO x / SiO 2 catalyst is associated with more favorable activation and kinetics of surface MoO x anchored at AlO x sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reactivity of Pd–MO 2 encapsulated catalytic systems for CO oxidation

In this work, we present an investigation aimed at characterizing and understanding the synergistic interactions in encapsulated catalytic structures between the metal core (i.e., Pd) and oxide shell (i.e., TiO 2 , ZrO 2 , and CeO 2 ). Encapsulated catalysts were synthesized using a two-step procedure involving the initial colloidal synthesis of Pd nanoparticles (NPs) capped by various ligands and subsequent sol–gel encapsulation of the NPs with porous MO 2 (M = Ti, Zr, Ce) shells. The encapsulated catalytic systems displayed higher activity than the Pd/MO 2 supported structures due to unique physicochemical properties at the Pd–MO 2 interface. Pd@ZrO 2 exhibited the highest catalytic activity for CO oxidation. Results also suggested that the active sites in Pd encapsulated by an amorphous ZrO 2 shell structure were significantly more active than the crystalline oxide encapsulated structures at low temperatures. Furthermore, CO DRIFTS studies showed that Pd redispersion occurred under CO oxidation reaction conditions and as a function of the oxide shell composition, being observed in Pd@TiO 2 systems only, with potential formation of smaller NPs and oxide-supported Pd clusters after reaction. This investigation demonstrated that metal oxide composition and (in some cases) crystallinity play major roles in catalyst activity for encapsulated catalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning Assisted HPC Workload Trace Generation for Leadership Scale Storage Systems

Monitoring and analyzing a wide range of I/O activities in an HPC cluster is important in maintaining mission-critical performance in a large-scale, multi-user, parallel storage system. Center-wide I/O traces can provide high-level information and fine-grained activities per application or per user running in the system. Studying such large-scale traces can provide helpful insights into the system. It can be used to develop predictive methods for making predictive decisions, adjusting scheduling policies, or providing decisions for the design of next-generation systems. However, sharing real-world I/O traces to expedite such research efforts leaves a few concerns; i) the cost of sharing the large traces is expensive due to this large size, and ii) privacy concern is an issue.We address such issues by building an end-to-end machine learn- ing (ML) workflow that can generate I/O traces for large-scale HPC applications. We leverage ML based feature selection and gener- ative models for I/O trace generation. The generative models are trained on I/O traces collected by the darshan I/O characterization tool over a period of one year. We present a two-step generation process consisting of two deep-learning models, called the feature generator and the trace generator. The combination of two-step generative models provides robustness by reducing the bias of the model and accounting for the stochastic nature of the I/O traces across different runs of an application. We evaluate the performance of the generative models and show that the two-step model can generate time-series I/O traces with less than 20% root mean square error.

Paul, Arnab↗

Design, Processing, and Properties of WTaCrV-Hf Multi-principal Element Alloys

Refractory multi-principal element alloys are candidates for high-temperature structural components due, in part, to their high strength and high melting points. Single-phase materials are initially preferred for isotropic material properties as a function of time and temperature in service conditions. This work outlines a computational rank-ordering and experimental validation methodology for single-phase body-centered-cubic phase stability in WTaCrV-Hf alloys using order–disorder transition temperature. Eight compositions were fabricated by arc-melting and heat-treated at 1400 °C for 24 hrs. X-ray diffraction, energy-dispersive x-ray spectroscopy, and Vickers hardness testing showed alloys with order–disorder transition temperatures below 600 °C formed a single-phase body-centered-cubic structure during solidification and remained single-phase after heat-treatment. The sample possessing the lowest order–disorder transition temperature exhibited slip traces suggestive of room-temperature plastic deformation under Vickers indentation, with both heat-treated single-phase samples exhibiting hardnesses over 800 HV with little cracking compared to tungsten. These results establish order–disorder transition temperature as a viable predictive parameter for multi-principal element alloy phase stability. The methodology outlined in this work provides a framework for future design, fabrication, and characterization of high-temperature structural multi-principal element alloys.

CALPHAD↗

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗

A comprehensive numerical investigation on spray models for Direct-Injection Spark-Ignition engines

Gasoline direct-injection spark-ignition (DISI) engines generate a large portion of their unburned hydrocarbon (UHC) and soot emissions during the cold-start phase. A predictive computational fluid dynamics (CFD) modeling framework can be used to understand the physical processes that characterize fuel spray evolution and fuel-film formation at cold start conditions, which can help to reduce engine-out particulate emissions. This study systematically evaluated spray submodels and developed a set of simulation best practices for physical-numerical submodels with the goal of enabling accurate simulations of liquid spray behavior in a DISI engine. Three comprehensive experimental datasets containing free-spray projected liquid volume (PLV), liquid volume fraction (LVF), and near-field X-ray radiography data were used to validate the simulation results and evaluate the spray submodels. Systematic analysis delved into injected parcel distribution, droplet collision, spray breakup, and evaporation via a detailed assessment of the relevant spray submodels. Moreover, the effects of turbulence models and the initial turbulent flow properties on the liquid spray evolution were examined. Based on extensive calibration efforts, a set of simulation best practices for the free spray was developed and validated against the PLV/LVF data. Simulation results indicated that the uniform distribution for parcel initialization, coupled with appropriate droplet collision submodels, provides an improved spray morphology compared to the cluster distribution. The findings also underscored the importance of calibrating the Kelvin-Helmholtz Rayleigh-Taylor (KH-RT) breakup model constants and droplet heat transfer coefficient scaling factor to achieve favorable agreement regarding measured liquid penetration and spray widths. In conclusion, this study marks a substantial stride towards accurately predicting fuel film evolution and soot formation within DISI engine performance.

ECN Spray G↗