Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional diffusion”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Anchoring Polar Magnetic Field in a Stationary Thick Accretion Disk

We investigate the properties of a hot accretion flow bathed in a poloidal magnetic field. We consider an axisymmetric viscous-resistive flow in the steady-state configuration. We assume that the dominant mechanism of energy dissipation is due to turbulence viscosity and magnetic diffusivity. A certain fraction of that energy can be advected toward the central compact object. We employ the self-similar method in the radial direction to find a system of ODEs with just one varible, θ in the spherical coordinates. For the existence and maintenance of a purely poloidal magnetic field in a rotating thick disk, we find that the necessary condition is a constant value of angular velocity along a magnetic field line. We obtain an analytical solution for the poloidal magnetic flux. We explore possible changes in the vertical structure of the disk under the influences of symmetric and asymmetric magnetic fields. Our results reveal that a polar magnetic field with even symmetry about the equatorial plane makes the disk vertically thin. Moreover, the accretion rate decreases when we consider a strong magnetic field. Finally, we notice that hot magnetized accretion flows can be fully advected even in a slim shape.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Environmental forcing does not induce diel or synoptic variation in the carbon isotope content of forest soil respiration

Recent studies have examined temporal fluctuations in the amount and carbon isotope content (δ 13 C) of CO 2 produced by the respiration of roots and soil organisms. These changes have been correlated with diel cycles of environmental forcing (e.g., sunlight and soil temperature) and with synoptic-scale atmospheric motion (e.g., rain events and pressure-induced ventilation). We used an extensive suite of measurements to examine soil respiration over 2 months in a subalpine forest in Colorado, USA (the Niwot Ridge AmeriFlux forest). Observations included automated measurements of CO 2 and δ 13 C of CO 2 in the soil efflux, the soil gas profile, and forest air. There was strong diel variability in soil efflux but no diel change in the δ 13 C of the soil efflux (δ R ) or the CO 2 produced by biological activity in the soil (δ J ). Following rain, soil efflux increased significantly, but δ R and δ J did not change. Temporal variation in the δ 13 C of the soil efflux was unrelated to measured environmental variables, and we failed to find an explanation for this unexpected result. Measurements of the δ 13 C of the soil efflux with chambers agreed closely with independent observations of the isotopic composition of soil CO 2 production derived from soil gas well measurements. Deeper in the soil profile and at the soil surface, results confirmed established theory regarding diffusive soil gas transport and isotopic fractionation. Deviation from best-fit diffusion model results at the shallower depths illuminated a pump-induced ventilation artifact that should be anticipated and avoided in future studies. There was no evidence of natural pressure-induced ventilation of the deep soil. However, higher variability in δ 13 C of the soil efflux relative to δ 13 C of production derived from soil profile measurements was likely caused by transient pressure-induced transport with small horizontal length scales.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Comparing dry aerosol size measurements

Two instruments were used to size dry aerosols for the CCN experiments: the Naval Research Laboratory (NRL) instrument which uses an electrical mobility analyzer to separate the aerosol size fractions and a diffusion chamber to grow and count the particles and the University of Wyoming (WYO) aerosol monitoring system. Measurements from these two systems were usually found to be in agreement during the CCN Workshop. Two examples are presented to compare aerosol size distribution measurements of the two instruments: experiment 8 (monodisperse NaCl) and experiment 27 (polydisperse AgI). Differential (dN/dR) and cumulative plots are shown for both instruments for experiment 8 and experiment 27; also shown are Aitken particle measurements for comparison.

Rogers, D. C.↗
NASA NTRS1991

Nitrogen abundance in Comet Halley

Data on the nitrogen-containing compounds that observed spectroscopically in the coma of Comet Halley are summarized, and the elemental abundance of nitrogen in the Comet Halley nucleus is derived. It is found that 90 percent of elemental nitrogen is in the dust fraction of the coma, while in the gas fraction, most of the nitrogen is contained in NH3 and CN. The elemental nitrogen abundance in the ice component of the nucleus was found to be deficient by a factor of about 75, relative to the solar photosphere, indicating that the chemical partitioning of N2 into NH3 and other nitrogen compounds during the evolution of the solar nebula cannot account completely for the low abundance ratio N2/NH3 = 0.1, observed in the comet. It is suggested that the low N2/NH3 ratio in Comet Halley may be explained simply by physical fractionation and/or thermal diffusion.

Wyckoff, Susan↗
NASA NTRS2012

A Miniaturized Laser Heterodyne Radiometer for a Global Ground-Based Column Carbon Monitoring Network

We present progress in the development of a passive, miniaturized Laser Heterodyne Radiometer (mini-LHR) that will measure key greenhouse gases (C02, CH4, CO) in the atmospheric column as well as their respective altitude profiles, and O2 for a measure of atmospheric pressure. Laser heterodyne radiometry is a spectroscopic method that borrows from radio receiver technology. In this technique, a weak incoming signal containing information of interest is mixed with a stronger signal (local oscillator) at a nearby frequency. In this case, the weak signal is sunlight that has undergone absorption by a trace gas of interest and the local oscillator is a distributive feedback (DFB) laser that is tuned to a wavelength near the absorption feature of the trace gas. Mixing the sunlight with the laser light, in a fast photoreceiver, results in a beat signal in the RF. The amplitude of the beat signal tracks the concentration of the trace gas in the atmospheric column. The mini-LHR operates in tandem with AERONET, a global network of more than 450 aerosol sensing instruments. This partnership simplifies the instrument design and provides an established global network into which the mini-LHR can rapidly expand. This network offers coverage in key arctic regions (not covered by OCO-2) where accelerated warming due to the release of CO2 and CH4 from thawing tundra and permafrost is a concern as well as an uninterrupted data record that will both bridge gaps in data sets and offer validation for key flight missions such as OCO-2, OCO-3, and ASCENDS. Currently, the only ground global network that routinely measures multiple greenhouse gases in the atmospheric column is TCCON (Total Column Carbon Observing Network) with 18 operational sites worldwide and two in the US. Cost and size of TCCON installations will limit the potential for expansion, We offer a low-cost $30Klunit) solution to supplement these measurements with the added benefit of an established aerosol optical depth measurement. Aerosols induce a radiative effect that is an important modulator of regional carbon cycles. Changes in the diffuse radiative flux fraction (DRF) due to aerosol loading have the potential to alter the terrestrial carbon exchange.

Wilson, Emily L.↗
DOE OSTI2022

Influence of the Soret effect on binary-species boundary layers at high pressure

Direct numerical simulations of binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate the influence of the Soret effect on flow physics of binary-species boundary layers where the fluid has a uniform composition. The working fluid is a mixture of 25 % methane and 75 % nitrogen in mass fraction. Although the fluid composition is uniform at the initial condition, the mass fraction of methane increases near the wall when the wall temperature is hotter than the free stream temperature, whereas it decreases when the wall temperature is colder. The non-uniform mass fraction indicates that the uphill diffusion occurs near the wall. Investigation of fluctuations of the mass fraction reveals that the mass fraction fluctuates in the whole boundary layer, indicating that the uphill diffusion occurs even far from the wall. Examination of the species-mass diffusion balance for mean flow fields clarified that the Soret effect flux becomes large near the wall, and the large flux causes the non-uniform profile of the mass fraction near the wall.

Binary Diffusion Coefficient↗
DOE OSTI2019

Porosity Prediction in A356 Wheel Casting

A FLUENT™-based thermal-fluid-compositional model has been developed and applied to a low-pressure, die-cast, A356 aluminum alloy wheel to explore the effects of macro-segregation of hydrogen and late-stage liquid encapsulation on pore formation during solidification. Based on the evolution of temperature, pressure, and hydrogen concentration in the liquid output from FLUENT, the porosity size distribution has been predicted at a number of locations in a commercially produced wheel using an in-house micro-porosity model. The porosity model incorporates pore nucleation kinetics based on hydrogen supersaturation, and two modes of pore growth: the first, occurring at low solid fractions is based on hydrogen diffusion; and the second, occurring at high solid fractions is based on conservation of volume. Samples from a commercially cast wheel have been analyzed using X-ray micro-tomography to provide basic validation of the micro-porosity model. The predicted results showed that the evolution in pressure has the dominant effect on pore growth, but only under conditions where pores have nucleated prior to late-stage liquid encapsulation. Otherwise, the cooling rate appears to have the dominant effect. The accuracy of the model is discussed in terms of the predicted pore size distribution and number density of pores. Areas for continued development are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Diffusion and perfusion MR parameters to assess preoperative short-course radiotherapy response in locally advanced rectal cancer: a comparative explorative study among Standardized Index of Shape by DCE-MRI, intravoxel incoherent motion- and diffusion kurtosis imaging-derived parameters

Purpose: To assess preoperative short-course radiotherapy (SCR) tumor response in locally advanced rectal cancer (LARC) by means of Standardized Index of Shape (SIS) by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), apparent diffusion coefficient (ADC), intravoxel incoherent motion (IVIM) and diffusion kurtosis imaging (DKI) parameters derived from diffusion-weighted MRI (DW-MRI). Materials and methods: Thirty-four patients with LARC who underwent MRI scans before and after SCR followed by delayed surgery, retrospectively, were enrolled. SIS, ADC, IVIM parameters [tissue diffusion (D{sub t}), pseudo-diffusion (D{sub p}), perfusion fraction (f{sub p})] and DKI parameters [mean diffusivity (MD), mean of diffusional kurtosis (MK)] were calculated for each patient. IVIM parameters were estimated using two methods, namely conventional biexponential fitting (CBFM) and variable projection (VARPRO). After surgery, the pathological TNM and tumor regression grade (TRG) were estimated. For each parameter, percentage changes between before and after SCR were evaluated. Furthermore, an artificial neural network was trained for outcome prediction. Nonparametric sample tests and receiver operating characteristic curve (ROC) analysis were performed. Results: Fifteen patients were classified as responders (TRG ≤ 2) and 19 as not responders (TRG > 3). Seven patients had TRG 1 (pathological complete response, pCR). Mean and standard deviation values of pre-treatment CBFM D{sub p} and mean value of VARPRO D{sub p} pre-treatment showed statistically significant differences to predict pCR. (p value at Mann–Whitney test was 0.05, 0.03 and 0.008, respectively.) Exclusively SIS percentage change showed significant differences between responder and non-responder patients after SCR (p value << 0.001) and to assess pCR after SCR (p value << 0.001). The best results to predict pCR were obtained by VARPRO Fp mean value pre-treatment with area under ROC of 0.84, a sensitivity of 96.4%, a specificity of 71.4%, a positive predictive value (PPV) of 92.9%, a negative predictive value (NPV) of 83.3% and an accuracy of 91.2%. The best results to assess after treatment complete pathological response were obtained by SIS with an area under ROC of 0.89, a sensitivity of 85.7%, a specificity of 92.6%, a PPV of 75.0%, a NPV of 96.1% and an accuracy of 91.2%. Moreover, the best results to differentiate after treatment responders vs. non-responders were obtained by SIS with an area under ROC of 0.94, a sensitivity of 93.3%, a specificity of 84.2%, a PPV of 82.4%, a NPV of 94.1% and an accuracy of 88.2%. Promising initial results were obtained using a decision tree tested with all ADC, IVIM and DKI extracted parameter: we reached high accuracy to assess pathological complete response after SCR in LARC (an accuracy of 85.3% to assess pathological complete response after SCR using VARPRO D{sub p} mean value post-treatment, ADC standard deviation value pre-treatment, MD standard deviation value post-treatment). Conclusion: SIS is a hopeful DCE-MRI angiogenic biomarker to assess preoperative treatment response after SCR with delayed surgery. Furthermore, an important prognostic role was obtained by VARPRO F{sub p} mean value pre-treatment and by a decision tree composed by diffusion parameters derived by DWI and DKI to assess pathological complete response.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Effects of surface diffusion in electrocatalytic CO 2 reduction on Cu revealed by kinetic Monte Carlo simulations

We report Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO 2 reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of 12 CO 2 and 13 CO on Cu; differences in the 13 C content of C2 products ethylene and ethanol (Δ 13 C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates 12 CO* and 13 CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ 13 C that can be observed. Δ 13 C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO 2 reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

The faint galaxy contribution to the diffuse extragalactic background light

Models of the faint galaxy contribution to the diffuse extragalactic background light are presented, which are consistent with current data on faint galaxy number counts and redshifts. The autocorrelation function of surface brightness fluctuations in the extragalactic diffuse light is predicted, and the way in which these predictions depend on the cosmological model and assumptions of biasing is determined. It is confirmed that the recent deep infrared number counts are most compatible with a high density universe (Omega-0 is approximately equal to 1) and that the steep blue counts then require an extra population of rapidly evolving blue galaxies. The faintest presently detectable galaxies produce an interesting contribution to the extragalactic diffuse light, and still fainter galaxies may also produce a significant contribution. These faint galaxies still only produce a small fraction of the total optical diffuse background light, but on scales of a few arcminutes to a few degrees, they produce a substantial fraction of the fluctuations in the diffuse light.

Cole, Shaun↗
DOE OSTI2021

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Gas Phase Spectroscopy of Cold PAH Ions: Contribution to the Interstellar Extinction and the Diffuse Interstellar Bands

Polycyclic Aromatic Hydrocarbon molecules (PAHs) are ubiquitous in the interstellar medium (ISM) and constitute the building blocks of interstellar dust grains. Despite their inferred important role in mediating the energetic and chemical processes in thc ISM, their exact contribution to the interstellar extinction, and in particular to the diffuse interstellar bands (DIBs) remains unclear. The DIBs are spectral absorption features observed in the line of sight of stars that are obscured by diffuse interstellar clouds. More than 200 bands have been reported to date spanning from the near UV to the near IR with bandwidths ranging from 0.4 to 40 Angstroms (Tielens & Snow 1995). The present consensus is that the DIBs arise from free flying, gas-phase, organic molecules and/or ions that are abundant under the typical conditions reigning in the diffuse ISM. PAHs have been proposed as possible carriers (Allamandola et al. 1985; Leger & DHendecourt 1985). The PAH hypothesis is consistent with the cosmic abundance of Carbon and Hydrogen and with the required photostability of the DIB carriers against the strong VUV radiation field in the diffuse interstellar clouds. A significant fraction of PAHs is expected to be ionized in the diffuse ISM.

Biennier, L.↗
NASA NTRS1992

Thermal and ultrasonic evaluation of porosity in composite laminates

The effects of porosity on damage incurred by low-velocity impact are investigated. Specimens of graphite/epoxy composite were fabricated with various volume fractions of voids. The void fraction was independently determined using optical examination and acid resin digestion methods. Thermal diffusivity and ultrasonic attenuation were measured, and these results were related to the void volume fraction. The relationship between diffusivity and fiber volume fraction was also considered. The slope of the ultrasonic attenuation coefficient was found to increase linearly with void content, and the diffusivity decreased linearly with void volume fraction, after compensation for an approximately linear dependence on the fiber volume fraction.

Johnston, Patrick H.↗
DOE OSTI2024

Toward a better understanding of Ni coarsening in solid oxide cells: NiH on Ni (111) examined at the level of density-functional theory (and KMC)

Ni coarsening in the fuel electrode of solid oxide cells (SOCs) is known to be significantly faster under a humid atmosphere. In this talk, density-functional theory and kinetic Monte Carlo methods are used to explore the hypothesis that the surface diffusion of NiH on Ni may promote Ni coarsening in the SOC operated in electrolysis cell mode. Defining the surface diffusivity as the product of the fractional surface coverage and single-molecule surface diffusivity, the surface diffusivity of NiH on Ni (111) is found to be large enough under a significant overpotential to support the above hypothesis. However, as NiH could dissociate on Ni (111), more work is needed to show that NiH may promote Ni coarsening.

Mantz, Yves↗
DOE OSTI2017

Importance of turbulence-chemistry interactions at low temperature engine conditions

The role of turbulence-chemistry interaction in autoignition and flame stabilization is investigated for spray flames at low temperature combustion (LTC) conditions by performing high-fidelity three-dimensional computational fluid dynamics (CFD) simulations. A recently developed Tabulated Flamelet Model (TFM) is coupled with a large eddy simulation (LES) framework and validated across a range of Engine Combustion Network (ECN) ambient temperature conditions for n-dodecane fuel. High resolution grids with 0.0625 mm minimum cell size and 25 million total cell count are implemented using adaptive mesh refinement over the spray and combustion regions. Simulations with these grids and multiple LES realizations, with a 103 species n-dodecane mechanism show good agreement with experimental data for all the ambient conditions investigated. This modeling approach with the computational cost advantage of tabulated chemistry is then extended towards understanding the auto-ignition and flame stabilization at an ambient temperature of 750 K. These low temperature conditions lead to substantially higher ignition delays and flame liftoff lengths, and significantly leaner combustion compared to conventional high temperature diesel combustion. These conditions also require the simulations to span significantly larger temporal and spatial dimensions thereby increasing the computational cost. The TFM approach is able to capture autoignition and flame liftoff length at the low temperature conditions. Significant differences with respect to mixing, species formation and flame stabilization are observed under low temperature compared to conventional diesel combustion. At higher ambient temperatures, formation of formaldehyde is observed in the rich region (phi > 1) followed by the formation of OH in the stoichiometric regions. Under low temperature conditions, formaldehyde is observed to form at leaner regions followed by the onset of OH formation in significantly lean regions of the flame. Qualitative differences between species formation and transient flame development for the high and low temperature conditions are presented. The two stage ignition process is further investigated by studying the species formation in mixture fraction space by solving 1D flamelet equations for different scalar dissipation rates and homogeneous reactor assumption. Results show that scalar dissipation causes these radicals to diffuse within the mixture fraction space. As a result, this significantly enhances ignition and plays a dominant role at such low temperature conditions which cannot be captured by the homogeneous reaction assumption based model.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Absolute age determination of impact cratering events remains difficult and often controversial; a challenge that has resulted in < 50% of known impact craters regarded as accurately and precisely dated. Besides conventional 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He (ZHe) dating of impactites has been applied to large- to medium-sized impact structures. ZHe dates can be fully reset in minutes at 1000°C, which is commonly reached in central sections of the melt sheet, whereas resetting of ZHe at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-6 years. There is a critical need to quantify the effects of shock-induced microstructures and impact metamorphism on helium diffusion kinetics in well-characterized, variably shocked zircon to further establish the reliability of (U-Th)/He for dating impacts. For this purpose, we investigated suevite and impact melt samples from two impact structures, the Chicxulub multi-ring basin and the Ries complex crater, which enables us to compare zircon helium diffusion kinetics from impact structures with differing sizes, ages, and hydrothermal system longevities. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incremental 10°C steps from 250°C to 600°C. Afterward, we characterize the diffusion domain sizes and their interconnectivity within the shocked zircon grains using electron backscatter diffraction (EBSD). We find that zircon with few shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features (PDFs) and granular textures classified by BSE and EBSD are characterized by a dramatic decrease in helium retentivity, similar to radiation damaged grains, due to a reduction in the effective domain size and the introduction of interconnected fast diffusion pathways created by shock microstructures. A subset of grains were dated by ZHe after the external morphology of the grains was determined by BSE imaging. The euhedral grains yielded a weighted mean age within the uncertainty of the accepted impact ages, whereas the grains with PDFs or granular textures exhibited younger ages. Thus, these diffusion experiments and ZHe dates suggest that the dramatic decrease in domain size likely renders shocked grains more susceptible to impact-induced hydrothermal resetting and subsequent overprinting. Hence, characterization of shock microstructures is critical for determining accurate impact ages using ZHe methods especially when applied to previously unconstrained craters. The thermochronometer also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

zircon↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗