SEARCH · Search NASA

Results for “coupling coefficient”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Surrogate Models for TREAT Transient Calculations

The Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) is currently preparing for resumption of operation after more than 20 years in a long-term standby state. In parallel with preparation for operations and a new transient testing program, the DOE-NE Advanced Modeling and Simulation program (NEAMS) has been supporting the development of three-dimensional transient multi-physics modeling capabilities within the reactor physics application MAMMOTH [1]. Recently, the MAMMOTH reactor physics team has begun development of uncertainty and sensitivity analysis for the TREAT core. This research will be used to support the analysis of reactor operation and experiment planning. Traditional approaches, i.e., Monte Carlo simulations, for such tasks require a large number of model evaluations under both steady state and transient conditions. However, exhaustive sets of sampling runs are simply not practical here since evaluations of each model can take a substantial amount of time. A solution to this problem is to build a surrogate model employing a smaller number of model evaluations. The surrogate model is designed to provide an accurate and fast approximation to the detailed TREAT model. In this manner, the surrogate model can be used in place of the TREAT transient model for uncertainty and sensitivity analysis at a fraction of the computational cost. The high dimensional model representation (HDMR) [2] method is employed here for approximating the responses of interest in terms of functions of lower dimensions. This method has been implemented in RAVEN [3] and all TREAT transient calculations are performed using the multiphysics code MAMMOTH. In addition, a code interface has also been developed in order to couple RAVEN with MAMMOTH. The surrogate model is then tested and compared for the MAMMOTH model of the TREAT startup testing program Transient Test #15 core [4], and its accuracy is measured by its ability to estimate the mean, coefficient of variance, skewness, excess kurtosis and probability density function of the outputs of Transient Test #15 model, as well as root mean square error metric.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Structure and tunable temperature coefficient of magnetization of Mn 4-x Ga x C alloys prepared by induction melting method

The magnetization of most magnetic materials decreases monotonically with increasing temperature. In this work, we found that the temperature coefficient of magnetization of Mn 4-x Ga x C alloys can be tuned from negative values to positive values by controlling the composition x. The antiperovskite type Mn 4-x Ga x C (0.05 ≤ x ≤ 0.75) alloys were prepared by using induction melting method, which is more efficient in large-scale production and obtaining full-density alloys in comparison with the traditional solid-state-reaction method. The values of the temperature coefficient of magnetization of Mn 4-x Ga x C change continuously from negative to positive with decreasing x. The Mn 4-x Ga x C alloys with highly thermal-stable magnetization is expected to present in the composition range of 0.15 < x < 0.25. The saturation magnetization of Mn 4-x Ga x C increases with increasing x, owing to the reduced number of antiferromagnetically coupled Mn atoms at the cubic corner with the face centered Mn atoms. Most Mn 4-x Ga x C alloys with varying x display near-zero remanent magnetization and coercivity at room temperature. The Currie temperature of Mn 4-x Ga x C decreases with increasing x. The x-ray photoelectron spectra of Mn 2p, Ga 2p, and C 1 s reveal distinct splitting due to the diverse chemical states of these atoms at different lattice positions and/or phases. Our work has developed a class of alloys capable of offering a desired temperature coefficient of magnetization across a broad temperature range, thereby offering a method to manipulate the thermodynamics of magnetization.

Magnetic properties↗

DOE OSTI2019

Impact of phonons and spin-orbit coupling on Auger recombination in InAs

Here we extend the methodology for calculating Auger recombination from first principles to include spin-orbit interaction, which allows us to study materials in which the spin-orbit splitting and the band gap are of similar magnitude. We use this methodology to compute the direct and indirect Auger coefficients in InAs and related alloys. The direct process involves only Coulomb interaction, while the indirect process is mediated by phonons. We show how the spin-orbit split-off states allow for efficient direct Auger recombination, greatly enhancing the excitation rate of Auger holes due to a near resonance between the spin-orbit splitting and the band gap of the material. We find that the direct hole-hole-electron Auger recombination coefficient decays exponentially with increasing band gap, while the electron-electron-hole Auger recombination coefficient decays faster than exponentially. Additionally, we modeled the phonon-assisted process in InAs and show that the direct Auger process is dominant.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Evaluation of the Uncertainty in Satellite-Based Crop State Variable Retrievals Due to Site and Growth Stage Specific Factors and Their Potential in Coupling with Crop Growth Models

Coupling crop growth models and remote sensing provides the potential to improve our understanding of the genotype x environment x management (G × E × M) variability of crop growth on a global scale. Unfortunately, the uncertainty in the relationship between the satellite measurements and the crop state variables across different sites and growth stages makes it difficult to perform the coupling. In this study, we evaluate the effects of this uncertainty with MODIS data at the Mead, Nebraska Ameriflux sites (US-Ne1, US-Ne2, and US-Ne3) and accurate, collocated Hybrid-Maize (HM) simulations of leaf area index (LAI) and canopy light use efficiency (LUE Canopy ). The simulations are used to both explore the sensitivity of the satellite-estimated genotype × management (G × M) parameters to the satellite retrieval regression coefficients and to quantify the amount of uncertainty attributable to site and growth stage specific factors. Additional ground-truth datasets of LAI and LUE Canopy are used to validate the analysis. The results show that uncertainty in the LAI/satellite measurement regression coefficients lead to large uncertainty in the G × M parameters retrievable from satellites. In addition to traditional leave-one-site-out regression analysis, the regression coefficient uncertainty is assessed by evaluating the retrieval performance of the temporal change in LAI and LUE Canopy . The weekly change in LAI is shown to be retrievable with a correlation coefficient absolute value (|r|) of 0.70 and root-mean square error (RMSE) value of 0.4, which is significantly better than the performance expected if the uncertainty was caused by random error rather than secondary effects caused by site and growth stage specific factors (an expected |r| value of 0.36 and RMSE value of 1.46 assuming random error). As a result, this study highlights the importance of accounting for site and growth stage specific factors in remote sensing retrievals for future work developing methods coupling remote sensing with crop growth models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Magnetotransport in semiconductors and two-dimensional materials from first principles

Here we demonstrate a first-principles method to study magnetotransport in materials by solving the Boltzmann transport equation (BTE) in the presence of an external magnetic field. Our approach employs ab initio electron-phonon interactions and takes spin-orbit coupling into account. We apply our method to various semiconductors (Si and GaAs) and two-dimensional (2D) materials (graphene) as representative case studies. The magnetoresistance, Hall mobility, and Hall factor in Si and GaAs are in very good agreement with experiments. In graphene, our method predicts a large magnetoresistance, consistent with experiments. Analysis of the steady-state electron occupations in graphene shows the dominant role of optical phonon scattering and the breaking of the relaxation time approximation. Our paper provides a detailed understanding of the microscopic mechanisms governing magnetotransport coefficients, establishing the BTE in a magnetic field as a broadly applicable first-principles tool to investigate transport in semiconductors and 2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Performance evaluation of CMIP6 models on the Arctic-Siberian Plain teleconnection affecting the East Asian heat waves

The frequency and intensity of summer heat waves in East Asia have increased sharply in recent decades, significantly impacting public health and the economy. The Arctic-Siberian Plain (ASP) teleconnection pattern has been identified as a key driver, with ASP warming amplifying atmospheric circulation patterns conducive to extreme temperatures. This study evaluates the ability of Coupled Model Inter-comparison Project phase 6 models to simulate the ASP pattern across interannual variability (IAV) and intra-seasonal variability (ISV) timescales using the Common Basis Function method. The multi-model mean shows statistically significant pattern correlations with ERA5 reanalysis, with correlation coefficients of 0.90 and 0.99 for IAV and ISV, respectively. While the ASP pattern is generally well captured, models exhibit substantial inter-model diversity in the intensity and position of anticyclonic anomalies over the ASP and East Asia. Models with ASP pattern variability similar to reanalysis better reproduce extreme East Asian temperatures, whereas those over- or underestimating ASP variability exhibit lower skill. These performance differences are related to differences in simulating key variables associated with the development of the ASP pattern. Our findings highlight the role of the ASP pattern in modulating extreme heat events, as models with improved ASP simulations align more closely with observed temperature extremes. Refining ASP representations in models could enhance seasonal heat wave predictions, improving climate adaptation strategies.

Arctic-Siberian Plain (ASP)↗

DOE OSTI2016

Dynamic simulation of CO 2 -injection-induced fault rupture with slip-rate dependent friction coefficient

It is known that poro-elastic stress and effective stress reduction associated with deep underground fluid injection can potentially trigger shear rupture along pre-existing faults. We modeled an idealized CO 2 injection scenario, to assess the effects on faults in the first phase of a generic CO 2 aquifer storage operation. We used coupled multiphase fluid flow and geomechanical numerical modeling to evaluate the stress and pressure perturbations induced by fluid injection and the response of a nearby normal fault. Slip-rate dependent friction and inertial effects have been taken into account during rupture. Contact elements have been used to take into account the frictional behavior of the rupture plane. We investigated different scenarios of injection rate to induce rupture on the fault, employing various fault rheologies. Published laboratory data on CO 2 -saturated intact and crushed rock samples, representative of a potential target aquifer, sealing formation and fault gouge, have been used to define a scenario where different fault rheologies apply at different depths. Nucleation of fault rupture takes place at the bottom of the reservoir, in agreement with analytical poro-elastic stress calculations, depending on injection-induced reservoir inflation and the tectonic stress scenario. For the stress state considered here, the first triggered rupture always produces the largest rupture length and slip magnitude, both of which correlate with the fault rheology. Velocity weakening produces larger ruptures and generates larger magnitude seismic events. Heterogeneous faults have been considered including velocity-weakening or velocity strengthening sections inside and below the aquifer, with the upper sections being velocity-neutral. Nucleation of rupture in a velocity-strengthening section results in a limited rupture extension, both in terms of maximum slip and rupture length. For a heterogeneous fault with nucleation in a velocity-weakening section, the rupture may propagate into the overlying velocity-neutral section, if the extent of velocity-weakening and associated friction drop are large enough.

58 GEOSCIENCES↗

DOE OSTI2019

Pressure-induced magnetovolume effect in CoCrFeAl high-entropy alloy

High-entropy alloys (HEAs) composed of multiple-principal elements with (nearly) equimolar ratio establish a new conceptual framework for alloy design and hold a promise for extensive applications in industry, akin to the controlled expansion alloys (CEAs), such as Invar alloys. Spontaneously, one question emerges - would it be possible to synthesize a novel class of alloys combining the virtues of both CEAs and HEAs? Here, we report the pressure-induced magnetovolume effect in the body-centered-cubic CoCrFeAl HEA coupled with magnetic phase transitions from ferromagnetic to paramagnetic, and to non-magnetic states, originating from the successive collapses of local magnetic moments of Co and Fe. The observed magnetovolume anomalies, occurring in a progressive way, tailor appreciably the coefficient of thermal expansion of CoCrFeAl. These results further strengthen HEAs’ anticipated potential for designing multifunctional materials in virtue of their multiple outstanding properties, and reveal possible routes for their future synthesis.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Rotational relaxation in molecular hydrogen and deuterium: Theory versus acoustic experiments

An explicit formulation of the rotational relaxation time in terms of state-to-state rate coefficients associated to inelastic collisions is reported. The state-to-state rates needed for the detailed interpretation of relaxation in H{sub 2} and D{sub 2}, including isotopic variant mixtures, have been calculated by solving the close-coupling Schrödinger equations using the H{sub 2}–H{sub 2} potential energy surface by Diep and Johnson [J. Chem. Phys. 112, 4465 (2000)]. Relaxation related quantities (rotational effective cross section, bulk viscosity, relaxation time, and collision number) calculated from first principles agree reasonably well with acoustic absorption experimental data on H{sub 2} and D{sub 2} between 30 and 293 K. This result confirms at once the proposed formulation, and the validation of the H{sub 2}–H{sub 2} potential energy surface employed, since no approximations have been introduced in the dynamics. Accordingly, the state-to-state rates derived from Diep and Johnson potential energy surface appear to be overestimated by up to 10% for H{sub 2}, and up to 30% for D{sub 2} at T = 300 K, showing a better agreement at lower temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Mixed quantum/classical theory for rotationally and vibrationally inelastic scattering of open-shell molecules and its application to the NH(X 3 Σ – ) + He collisional system

We developed the mixed quantum/classical theory (MQCT) for rotationally and vibrationally inelastic scattering of a diatomic molecule in a 2S+1 Σ electronic state with a closed-shell atom. In this approach, the rotational (and vibrational) fine-structure states of the molecule are treated quantum-mechanically, whereas the relative motion of collisional partners is treated classically. This theory is benchmarked against full quantum (close-coupling) calculations for the NH(X 3 Σ – ) + He system. Good agreement is found in all cases, and MQCT reproduces all main features of energy dependence of the cross section, except narrow scattering resonances at very low energies. Our method recovers propensity rules that govern values of close-coupling cross sections. Particularly, MQCT correctly predicts that the magnitudes of F -conserving transitions are significantly larger than those of F -changing ones. Furthermore, this new development makes MQCT a viable candidate for obtaining fine-structure resolved rotational rate coefficients at higher temperatures and for complex polyatomic systems where the standard full quantum treatment is computationally infeasible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Enhancement of the effective mass at high magnetic fields in CeRhIn 5

The Kondo-lattice compound CeRhIn 5 displays a field-induced Fermi surface reconstruction at B*≈ 30 T, which occurs within the antiferromagnetic state, prior to the quantum critical point at B c0 ≈ 50 T. Here, in order to investigate the nature of the Fermi surface change, we measured the magnetostriction, specific heat, and magnetic torque of CeRhIn 5 across a wide range of magnetic fields. Our observations uncover the field-induced itineracy of the 4 f electrons, where above B onset ≈ 17 T there is a significant enhancement of the Sommerfeld coefficient, and spin-dependent effective cyclotron masses determined from quantum oscillations. Upon crossing B onset , the temperature dependence of the specific heat also shows distinctly different behavior from that at low fields. Finally, our results indicate that the Kondo coupling is remarkably robust upon increasing the magnetic field. This is ascribed to the delocalization of the 4 f electrons at the Fermi surface reconstruction at B*.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Requirements Description of DASSH-F

This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Induced spin-accumulation and spin-polarization in a quantum-dot ring by using magnetic quantum dots and Rashba spin-orbit effect

The effect of magnetic contacts on spin-dependent electron transport and spin-accumulation in a quantum ring, which is threaded by a magnetic flux, is studied. The quantum ring is made up of four quantum dots, where two of them possess magnetic structure and other ones are subjected to the Rashba spin-orbit coupling. The magnetic quantum dots, referred to as magnetic quantum contacts, are connected to two external leads. Two different configurations of magnetic moments of the quantum contacts are considered; the parallel and the anti-parallel ones. When the magnetic moments are parallel, the degeneracy between the transmission coefficients of spin-up and spin-down electrons is lifted and the system can be adjusted to operate as a spin-filter. In addition, the accumulation of spin-up and spin-down electrons in non-magnetic quantum dots are different in the case of parallel magnetic moments. When the intra-dot Coulomb interaction is taken into account, we find that the electron interactions participate in separation between the accumulations of electrons with different spin directions in non-magnetic quantum dots. Furthermore, the spin-accumulation in non-magnetic quantum dots can be tuned in the both parallel and anti-parallel magnetic moments by adjusting the Rashba spin-orbit strength and the magnetic flux. Thus, the quantum ring with magnetic quantum contacts could be utilized to create tunable local magnetic moments which can be used in designing optimized nanodevices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗

DOE OSTI2022

Process-based modeling of soil nitrous oxide emissions from United States corn fields under different management and climate scenarios coupled with evaluation using regional estimates

Direct emissions of soil nitrous oxide during a growing season (N 2 O gs ) can be quantified with process-based models considering interactions between management, climate, and soil moisture when key data are available. We used an adapted “parameterized CENTURY/DAYCENT-model” ( p CENTURY) calibrated with crop growth and soil organic matter decay coefficients at the county-level for the estimation of N 2 O gs in the United States Corn Belt. Model estimated N 2 O-emissions from corn-based biofuels scenarios considering crop rotation, fertilizer inputs, tillage, and weather were compared against meta-summary of field observations from 55 studies. Both model and meta-summary ranked N 2 O gs -emissions to be corn > wheat > soybean phase while model likely underestimated cover crop N 2 O gs -emissions. The N 2 O gs -emissions and the associated emission factors (EFs) were modeled and summarized to be greater after anhydrous ammonia than urea application and from conventional tilled than non-tilled fields. Modeled and observed N 2 O gs -emissions after organic and inorganic fertilizer amendment did not differ due to high variability associated with the treatments. However, the organic fertilizer associated EFs were greater according to meta-summary data because of N input rates. Regionalized weather scenarios indicate hotspots for N 2 O gs -emissions can occur where crop N uptake is limited during dry years and in eastern states also during normal or wet seasons. The p CENTURY-derived N 2 O gs EFs (0.91 ± 0.19%) for counties investigated were only slightly lower than literature (1.07 ± 0.57%) or Tier-1 (1%) values. Our preliminary evaluation of regional soil moisture estimates showed reasonable agreement between monthly soil moisture estimates and the North American Soil Moisture Dataset during the growing season, but overestimation of soil moisture in winter-spring can influence the estimates of annual N 2 O emissions so future work is needed to calibrate soil moisture-associated model parameters. Our work provided scenario-based estimates of climate and management impacts on soil N 2 O gs -emissions together with valuable spatial insights into EFs that will be improved by more accurate information of fertilizer inputs and more temporally refined model evaluation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Combination of the W boson polarization measurements in top quark decays using ATLAS and CMS data at $\sqrt{s} =$ 8 TeV

The combination of measurements of the W boson polarization in top quark decays performed by the ATLAS and CMS collaborations is presented. The measurements are based on proton-proton collision data produced at the LHC at a centre-of-mass energy of 8 TeV, and corresponding to an integrated luminosity of about 20 fb$^{−1}$ for each experiment. The measurements used events containing one lepton and having different jet multiplicities in the final state. The results are quoted as fractions of W bosons with longitudinal (F$_{0}$), left-handed (F$_{L}$), or right-handed (F$_{R}$) polarizations. The resulting combined measurements of the polarization fractions are F$_{0}$ = 0.693 ± 0.014 and F$_{L}$ = 0.315 ± 0.011. The fraction F$_{R}$ is calculated from the unitarity constraint to be F$_{R}$ = −0.008 ± 0.007. These results are in agreement with the standard model predictions at next-to-next-to-leading order in perturbative quantum chromodynamics and represent an improvement in precision of 25 (29)% for F$_{0}$ (F$_{L}$) with respect to the most precise single measurement. A limit on anomalous right-handed vector (V$_{R}$), and left- and right-handed tensor (g$_{L}$, g$_{R}$) tWb couplings is set while fixing all others to their standard model values. The allowed regions are [−0.11, 0.16] for V$_{R}$, [−0.08, 0.05] for g$_{L}$, and [−0.04, 0.02] for g$_{R}$, at 95% confidence level. Limits on the corresponding Wilson coefficients are also derived.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

A continuous analogue of the tensor-train decomposition

We develop new approximation algorithms and data structures for representing and computing with multivariate functions using the functional tensor-train (FT), a continuous extension of the tensor-train (TT) decomposition. The FT represents functions using a tensor-train ansatz by replacing the three-dimensional TT cores with univariate matrix-valued functions. The main contribution of this paper is a framework to compute the FT that employs adaptive approximations of univariate fibers, and that is not connected to any tensorized discretization. The algorithm can be coupled with any univariate linear or nonlinear approximation procedure. We demonstrate that this approach can generate multivariate function approximations that are several orders of magnitude more accurate, for the same cost, than those based on the conventional approach of compressing the coefficient tensor of a tensor-product basis. Our approach is in the spirit of other continuous computation packages such as Chebfun, and yields an algorithm which requires the computation of “continuous” matrix factorizations such as the LU and QR decompositions of vector-valued functions. To support these developments, we describe continuous versions of an approximate maximum-volume cross approximation algorithm and of a rounding algorithm that re-approximates an FT by one of lower ranks. We demonstrate that our technique improves accuracy and robustness, compared to TT and quantics-TT approaches with fixed parameterizations, of high-dimensional integration, differentiation, and approximation of functions with local features such as discontinuities and other nonlinearities.

42 ENGINEERING↗

DOE OSTI2024

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↗