Search NASA⌕ Search

SEARCH · Search NASA

Results for “effective volume”

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 271 records · Page 15

Effective pair potential for random fcc CoCrFeMnNi alloys

The single-phase equiatomic CoCrFeMnNi alloy is a random solid solution of five elements on the face-centered cubic lattice, whose pure constituents crystallize in very different structures and exhibit diverse magnetic properties. Due to the randomness of the alloy, 80% of nearest neighbor bonds are between unlike elements and thus the details of bonding in pure structures are less important. The elastic moduli of this alloy give rise to small Cauchy pressure C 12 - C 44 , which suggests that the dominant part of bonding may be described by a simple pair potential. We test this hypothesis by developing a long-range Lennard-Jones potential in which the equilibrium crystal structures of pure constituents are taken as reference. The standard mixing rules for regular solutions are then adopted to obtain parameters for bonds between unlike elements in the quinary system. The transferability of this potential to quaternary CoCrFeNi, ternary CoCrNi, and binary FeNi alloys is investigated and the predictions compared with experiments and density functional theory calculations. By sampling over a large number of random configurations, we investigate the effect of compositional randomness on misfit volumes, energies of point defects and stacking faults, and the dislocation friction stresses experienced by moving edge and screw dislocations.

Materials Science↗

Impact of property covariance on cluster weak lensing scaling relations

ABSTRACT We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, ΔΣ, and the ‘true’ cluster galaxy number count, Ngal, as measured within a spherical volume that is void of projection effects. By quantifying the impact of this covariance on mass calibration, this work reveals a significant source of systematic uncertainty. Using the MDPL2 simulation with galaxies traced by the SAGE semi-analytic model, we measure the intrinsic property covariance between these observables within the three-dimensional vicinity of the cluster, spanning a range of dynamical mass and redshift values relevant for optical cluster surveys. Our results reveal a negative covariance at small radial scales (R ≲ R200c) and a null covariance at large scales (R ≳ R200c) across most mass and redshift bins. We also find that this covariance results in a $2{\!-\!}3~{{\ \rm per\ cent}}$ bias in the halo mass estimates in most bins. Furthermore, by modelling Ngal and ΔΣ as multi-(log)-linear equations of secondary halo properties, we provide a quantitative explanation for the physical origin of the negative covariance at small scales. Specifically, we demonstrate that the Ngal–ΔΣ covariance can be explained by the secondary properties of haloes that probe their formation history. We attribute the difference between our results and the positive bias seen in other works with (mock)-cluster finders to projection effects. These findings highlight the importance of accounting for the covariance between observables in cluster mass estimation, which is crucial for obtaining accurate constraints on cosmological parameters.

Astronomy & Astrophysics↗

Halide segregation to boost all-solid-state lithium-chalcogen batteries

Mixing electroactive materials, solid-state electrolytes, and conductive carbon to fabricate composite electrodes is the most practiced but least understood process in all-solid-state batteries, which strongly dictates interfacial stability and charge transport. Here, we report on universal halide segregation at interfaces across various halogen-containing solid-state electrolytes and a family of high-energy chalcogen cathodes enabled by mechanochemical reaction during ultrahigh-speed mixing. Bulk and interface characterizations by multimodal synchrotron x-ray probes and cryo–transmission electron microscopy show that the in situ segregated lithium halide interfacial layers substantially boost effective ion transport and suppress the volume change of bulk chalcogen cathodes. Various all-solid-state lithium-chalcogen cells demonstrate utilization close to 100% and extraordinary cycling stability at commercial-level areal capacities.

36 MATERIALS SCIENCE↗

Collective Summary of sCO2 Materials Development (Supercritical Transformational Electric Power Generation (STEP) Level 2 Milestone Report) (Parts I - II)

Polymers such as PTFE (polytetrafluorethylene or Teflon), EPDM (ethylene propylene diene monomer) rubber, FKM fluoroelastomer (Viton), Nylon 11, Nitrile butadiene (NBR) rubber, hydrogenated nitrile rubber (HNBR) and perfluoroelastomers (FF_202) are commonly employed in super critical CO 2 (sCO2) energy conversion systems. O-rings and gaskets made from these polymers face stringent performance conditions such as elevated temperatures, high pressures, pollutants, and corrosive humid environments. In FY 2019, we conducted experiments at high temperatures (100°C and 120°C) under isobaric conditions (20 MPa). Findings showed that elevated temperatures accelerated degradation of polymers in sCO2, and that certain polymer microstructures are more susceptible to degradation over others. In FY 2020, the focus was to understand the effect of sCO2 on polymers at low (10 MPa) and high pressures (40 MPa) under isothermal conditions (100°C). It was clear that the same selectivity was observed in these experiments wherein certain polymeric functionalities showed more propensity to failure over others. Fast diffusion, supported by higher pressures and long exposure times (1000 hours) at the test temperature, caused increased damage in sCO2 environments to even the most robust polymers. We also looked at polymers under compression in sCO2 at 100°C and 20 MPa pressure to imitate actual sealing performance required of these materials in sCO2 systems. Compression worsened the physical damage that resulted from chemical attack of the polymers under these test conditions. In FY 2021, the effect of cycling temperature (from 50°C to 150°C to 50°C) for polymers under a steady sCO2 pressure of 20 MPa was studied. The aim was to understand the influence of cycling temperatures of sCO2 for typical polymers under isobaric (20 MPa) conditions. Thermoplastic polymers (Nylon, and PTFE) and elastomers (EPDM, Viton, Buna N, Neoprene, FF202, and HNBR) were subjected to 20 MPa sCO2 pressure for 50 cycles and 100 cycles in separate experiments. Samples were extracted for ex-situ characterization at 50 cycles and upon the completion of 100 cycles. Each cycle constituted of 175 minutes of cycling from 50°C to 150°C. The polymer samples were examined for physical and chemical changes by Dynamic Mechanical and Thermal Analysis (DMTA), Fourier Transform Infrared (FTIR) spectroscopy, and compression set. Density and mass changes immediately after removal from test were measured for degree of swell comparisons. Optical microscopy techniques and micro computer tomography (micro CT) images were collected on select specimens. Evaluations conducted showed that exposures to super-critical CO2 environments resulted in combinations of physical and/or chemical changes. For each polymer, the dominance of cycling temperatures under sCO2 pressures, were evaluated. Attempts were made to qualitatively link the permanent sCO2 effects to polymer micro- structure, free volume, backbone substitutions, presence of polar groups, and degree of crystallinity differences. This study has established that soft polymeric materials are conducive to failure in sCO2 through mechanisms of failure that are dependent on polymer microstructure and chemistry. Polar pendant groups, large atom substitutions on the backbone are some of the factors that are influential structural factors.

36 MATERIALS SCIENCE↗

Rational optimization of substituted α–MnO 2 cathode for aqueous zinc–ion battery

Density functional theory (DFT) is utilized to explore the effects of increasing concentrations of vanadium (V) and chromium (Cr) substitution on the discharge of α-MnO 2 cathode in hydrated zinc-ion batteries (ZIBs). During H + -intercalation and Zn 2+ -intercalation, Cr-substitution proves to be more effective than V-substitution in promoting discharge behaviors. Transitioning from Mn 0.875 Cr 0.125 O 2 , Mn 0.75 Cr 0.25 O 2 to Mn 0.625 Cr 0.375 O 2 is found to consistently enhance the discharge voltage, along with improved tunnel structure retention and volume expansion suppression. In comparison, the promoting effect of increasing V-substitution is relatively small at initial discharge stages and leads to degradation at later stages, primarily due to an increased concentration of unstable Mn 2+ ion. The superior effect of Cr-substitution is attributed to the unique atomic and electronic structures of substituted Cr 4+ and reduced Cr 3+ ions during discharge. These ions serve as active electron acceptors to limit the formation of Mn 3+ and Mn 2+ ions, and as anchors to stabilize the α-MnO 2 framework and intercalated H + /Zn 2+ ions, respectively. Finally, our study highlights fine-tuning through substitution to enhance the performance of α-MnO 2 -based cathode materials in ZIBs.

25 ENERGY STORAGE↗

Asymmetric Miktoarm Star Polymers as Polyester Thermoplastic Elastomers

A library of polyester-based A(BA') n asymmetric miktoarm star polymers was synthesized with A, A' = poly(l-lactide) (PLLA) as the semicrystalline hard blocks and B = poly(4-methylcaprolactone) (PMCL) as the soft segment using a grafting-through platform known as μSTAR. The synthetic versatility of μSTAR enabled a systematic investigation of architectural design parameters, in particular the number of BA' arms (n), while holding the A, A', and B block lengths constant. The value of n has a pronounced impact on the mechanical properties of these high-molecular-weight miktoarm materials. Tensile toughness increases with n, an effect likely related to bridging, as the modulus drops because the hard-block volume fraction decreases. Furthermore, these insights expand our understanding of architecture effects in sustainable thermoplastic elastomers.

36 MATERIALS SCIENCE↗

Validation and Mechanistic Studies of the Headspace Effect in MOF Synthesis

Abstract The influence of headspace volume and composition on metal–organic framework (MOF) synthesis has been largely overlooked with no mechanisms currently invoking these factors. Herein, we reveal that an increased volume of oxygen in the headspace during the synthesis of zinc carboxylate‐based MOFs (MOF‐5, MOF‐177, and IRMOF‐8) as well as UiO‐66 and ZJU‐28 accelerates crystal formation. Investigation into the underlying mechanism of this phenomenon elucidates the critical roles played by all reagents, and some of their byproducts, in MOF synthesis. These insights are the key to developing more efficient and sustainable synthetic strategies for MOF production.

Donovan, Julia C. [Department of Chemistry Institu↗

Validation and Mechanistic Studies of the Headspace Effect in MOF Synthesis

Abstract The influence of headspace volume and composition on metal–organic framework (MOF) synthesis has been largely overlooked with no mechanisms currently invoking these factors. Herein, we reveal that an increased volume of oxygen in the headspace during the synthesis of zinc carboxylate‐based MOFs (MOF‐5, MOF‐177, and IRMOF‐8) as well as UiO‐66 and ZJU‐28 accelerates crystal formation. Investigation into the underlying mechanism of this phenomenon elucidates the critical roles played by all reagents, and some of their byproducts, in MOF synthesis. These insights are the key to developing more efficient and sustainable synthetic strategies for MOF production.

Donovan, Julia C. [Department of Chemistry Institu↗

First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics

Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.

36 MATERIALS SCIENCE↗

Volatile aldehyde emissions from “sub-ohm” vaping devices

"Sub-ohm" atomizers with reduced resistance can deliver more power than conventional electronic cigarettes. Typical battery outputs are 100 watts or more. These devices are particularly popular among young users, and can be a significant source of volatile carbonyls in the indoor environment. Emissions from next-generation sub-ohm vaping products were characterized by determining e-liquid consumption and volatile aldehydes emissions for several combinations of popular high-power configurations. Tests explored the effect of dilution air flow (air vent opening), puffing volume, and coil assembly configuration. The mass of liquid consumed per puff increased as the puff volume increased from 50 to 100 mL, then remained relatively constant for larger puff volumes up to 500 mL. This is likely due to mass transfer limitations at the wick and coil assembly, which reduced the vaporization rate at higher puff volumes. Carbonyl emission rates were systematically evaluated using a 0.15 Ω dual coil atomizer as a function of the puffing volume and dilution air flow, adjusted by setting the air vents to either 100% (fully open), 50%, 25%, or 0% (closed). The highest formaldehyde emissions were observed for the lowest puff volume (50 mL) when the vents 31 were closed (48 ng mg-1), opened at 25% (39 ng mg -1 ) and at 50% (32 ng mg -1 ).. By contrast, 50-mL puffs with 100% open vents, and puff volumes >100 mL for any vent aperture, generated formaldehyde yields of 20 ng mg -1 or lower, suggesting that a significant cooling effect resulted in limited carbonyl formation. Considering the effect of the coil resistance when operated at a voltage of 3.8 V, the amount of liquid evaporated per puff decreased as the resistance increased, in the order of 0.15 Ω > 0.25 Ω > 0.6 Ω consistent with decreasing aerosol temperatures measured at the mouthpiece. Three different configurations of 0.15 Ω coils (dual, quadruple and octuple) were evaluated, observing significant variability. No clear trend was found between carbonyl emission rates and coil resistance or configuration, with highest emissions corresponding to a 0.25 Ω dual coil atomizer. Carbonyl emission rates were compared with those determined using the same methodology for conventional e-cigarettes (lower power tank systems), observing overall lower yields for the sub-ohm devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potential-induced wetting and dewetting in pH-responsive block copolymer membranes for mass transport control

Wetting and dewetting behavior in channel-confined hydrophobic volumes is used in biological membranes to effect selective ion/molecular transport. Artificial biomimetic hydrophobic nanopores have been devised utilizing wetting and dewetting, however, tunable mass transport control utilizing multiple transport modes is required for applications such as controllable release/transport, water separation/purification and energy conversion. Here, we investigate the potential-induced wetting and dewetting behavior in a pH-responsive membrane composed of a polystyrene-b-poly(4-vinylpyridine) (PS-b-P4VP) block copolymer (BCP) when fabricated as a hierarchically-organized sandwich structure on a nanopore electrode array (NEA), i.e. BCP@NEA. At pH < pK a (P4VP) (pK a ~ 4.8), the BCP acts as an anion-exchange membrane due to the hydrophilic, protonated P4VP cylindrical nanodomains, but at pH > pK a (P4VP), the P4VP domains exhibit charge-neutral, hydrophobic and collapsed structures, blocking mass transport via the hydrophobic membrane. However, when originally prepared in a dewetted condition, mass transport in the BCP membrane may be switched on if sufficiently negative potentials are applied to the BCP@NEA architecture. When the hydrophobic BCP membrane is introduced on top of 2-electrode-embedded nanopore arrays, electrolyte solution in the nanopores is introduced, then isolated, by exploiting the potential-induced wetting and dewetting transitions in the BCP membrane. The potential-induced wetting/dewetting transition and the effect on cyclic voltammetry in the BCP@NEA structures is characterized as a function of the potential, pH and ionic strength. Additionally chronoamperometry and redox cycling experiments are used to further characterize the potential response. The multi-modal mass transport system proposed in this work will be useful for ultrasensitive sensing and single-molecule studies, which require long-time monitoring to explore reaction dynamics as well as molecular heterogeneity in nanoconfined volumes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗

Toward Water Volume Fraction Calculation in Multiphase Flows Using Electrical Capacitance Tomography Sensors

Real-time monitoring of water volume fraction in multiphase flows is an important problem for a number of industrial applications. The water phase in the multiphase flows may correspond to either the dispersed phase or the continuous phase. In the past, several lowcost and nonintrusive techniques based on the electrical capacitance tomography (ECT) has been developed to image and monitor in real-time multiphase flows containing water. Furthermore, such monitoring becomes increasingly challenging for high salinity levels, and no reliable ECTbased method is presently available which could work for obtaining water volume fraction in multiphase flows for all water salinity levels. In this paper, we propose a new approach based on the Hanai’s formula for complex dielectric constant and taking advantage of the Maxwell- Wagner-Sillars effect to obtain, to a good approximation, water volume fractions in multiphase flows containing water as either dispersed or continuous phase.

02 PETROLEUM↗

Evaluation of the economic implications of varied pressure drawdown strategies generated using a real-time, rapid predictive, multi-fidelity model for unconventional oil and gas wells

Experience has suggested that pressure maintenance in hydraulically fractured reservoirs via lower, more sustained production drawdowns may offer improved cumulative recovery and overall resource extraction efficiency compared to more rapid drawdown approaches aimed at generating high initial production. However, given the inherent variability of oil and natural gas markets, operators pursue production strategies that maximize profitability over resource extraction efficiency. This study focuses on evaluating the implications of contrasting pressure drawdown strategies on the long-term production and resulting economics for a real, producing unconventional gas well in the Marcellus Shale of the Appalachian Basin using a techno-economic analysis approach. Our research combines elements of well-specific horizontal well design, production forecasting, equipment sizing and capital cost estimation, operating cost estimation, and revenue and tax calculations. Gas production forecast outlook scenarios were generated under varying pressure drawdowns using two approaches: 1) a novel physics-informed machine learning workflow and 2) traditional reservoir simulation. A discounted cash flow model was used to evaluate the resulting economic implications for each drawdown scenario—generating output for exploring the coupled effect of factors like the timing and volume of gas production, prevailing economic and market conditions for natural gas, and overall estimated ultimate recovery on profitability metrics such as internal rate of return and net present value. Results show that there is potential to maximize the cumulative gas produced in the specific case study well by employing a lower pressure drawdown. Conversely, the greatest profitability is achieved using rapid drawdown as signified by a small, specific subset of our outlook scenarios. On an averaging basis, we find that the combinations of highest cumulative producing and most profitable scenarios occur under lower drawdowns with long (>40 years) producing timeframes, but require higher relative gas price and lower discounting considerations. Further, the machine learning predictive outlooking capability proved effective for enabling rapid generation of a multitude of scenario forecasts. As a result, a variety of prominent example cases could be generated to strike the balance of greater productivity and economic return given their associated producing features and economic conditions when compared to similar producing scenarios—critical insight that offers improved decision support for unconventional oil and gas operations.

42 ENGINEERING↗

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗

Dose Rate Analysis of the WCS Consolidated Interim Storage Facility

This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗

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)↗