Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “benchmark correlations”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fundamental gap of fluorographene by many-body GW and fixed-node diffusion Monte Carlo methods

Fluorographene (FG) is a promising graphene-derived material with a large bandgap. Currently existing predictions of its fundamental gap (Δ f ) and optical gap (Δ opt ) significantly vary when compared with experiment. We provide here an ultimate benchmark of Δ f for FG by many-body GW and fixed-node diffusion Monte Carlo (FNDMC) methods. Both approaches independently arrive at Δf ≈ 7.1 ± 0.1 eV. In addition, the Bethe–Salpeter equation enabled us to determine the first exciton binding energy, Eb = 1.92 eV. We also point to the possible misinterpretation problem of the results obtained for gaps of solids by FNDMC with single-reference trial wave functions of Bloch orbitals. We argue why instead of Δ opt , in the thermodynamic limit,

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Explicitly correlated coupled cluster method for accurate treatment of open-shell molecules with hundreds of atoms

We present a near-linear scaling formulation of the explicitly correlated coupled-cluster singles and doubles with the perturbative triples method [CCSD(T)F12¯] for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g., the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the C160H322 n-alkane, with the measured size complexity of ≈1.2). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)F12¯ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580–6602 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)], the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)F12¯ yielded a RMSD of ∼0.4 kcal/mol with respect to the CBS CCSD(T) estimate.

Kumar, Ashutosh (ORCID:0000000175896030)↗
NASA NTRS1993

Numerical Study of Base Pressure Characteristic Curve for a Four-Engine Clustered Nozzle Configuration

Excessive base heating has been a problem for many launch vehicles. For certain designs such as the direct dump of turbine exhaust in the nozzle section and at the nozzle lip of the Space Transportation Systems Engine (STME), the potential burning of the turbine exhaust in the base region has caused tremendous concern. Two conventional approaches have been considered for predicting the base environment: (1) empirical approach, and (2) experimental approach. The empirical approach uses a combination of data correlations and semi-theoretical calculations. It works best for linear problems, simple physics and geometry. However, it is highly suspicious when complex geometry and flow physics are involved, especially when the subject is out of historical database. The experimental approach is often used to establish database for engineering analysis. However, it is qualitative at best for base flow problems. Other criticisms include the inability to simulate forebody boundary layer correctly, the interference effect from tunnel walls, and the inability to scale all pertinent parameters. Furthermore, there is a contention that the information extrapolated from subscale tests with combustion is not conservative. One potential alternative to the conventional methods is computational fluid dynamics (CFD), which has none of the above restrictions and is becoming more feasible due to maturing algorithms and advancing computer technology. It provides more details of the flowfield and is only limited by computer resources. However, it has its share of criticisms as a predictive tool for base environment. One major concern is that CFD has not been extensively tested for base flow problems. It is therefore imperative that CFD be assessed and benchmarked satisfactorily for base flows. In this study, the turbulent base flowfield of a experimental investigation for a four-engine clustered nozzle is numerically benchmarked using a pressure based CFD method. Since the cold air was the medium, accurate prediction of the base pressure distributions at high altitudes is the primary goal. Other factors which may influence the numerical results such as the effects of grid density, turbulence model, differencing scheme, and boundary conditions are also being addressed. Preliminary results of the computed base pressure agreed reasonably well with that of the measurement. Basic base flow features such as the reverse jet, wall jet, recompression shock, and static pressure field in plane of impingement have been captured.

Wang, Ten-See↗
NASA NTRS1993

Fuel injector: Air swirl characterization aerothermal modeling, phase 2, volume 1

A well integrated experimental/analytical investigation was conducted to provide benchmark quality relevant to a prefilming type airblast fuel nozzle and its interaction with the combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM), and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗
NASA NTRS1993

Fuel Injector: Air swirl characterization aerothermal modeling, phase 2, volume 2

A well integrated experimental/analytical investigation was conducted to provide benchmark quality data relevant to prefilming type airblast fuel nozzle and its interaction with combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) equipment was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM) and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗
NASA NTRS2012

The Violent Early Solar System, as Told by Lunar Sample Geochronology

One of the legacies of the samples collected by the Apollo and Luna missions is the link forged between radiometric ages of rocks and relative ages according to stratigraphic relationships and impact crater size-frequency distributions. Our current understanding of the history of the inner solar system is based on the relative chronology of individual planets, tied to the absolute geochronology of the Moon via these important samples. Samples from these nearside locations reveal a preponderance of impact-disturbed or recrystallized ages between 3.75 and 3.95 billion years. Argon and lead loss (and correlated disturbances in the Rb-Sr system) have been attributed to metamorphism of the lunar crust by an enormous number of impacts in a brief pulse of time, called the Lunar Cataclysm or Late Heavy Bombardment. Subsequent high-precision geochronometric analyses of Apollo samples and lunar highlands meteorites show a wider range of ages, but very few older than 4 Ga. The paucity of ancient impact melt rocks has been interpreted to mean that either that most impact basins formed at this time, or that ejecta from the large, near-side, young basins dominates the Apollo samples. Selenochronology is getting more complicated: new results question meaning of sample ages, crater counts, crater production functions, and the solar system itself. Improved geological mapping of lunar geologic units and boundaries using multiple remote sensing datasets. High-resolution image-based crater counting of discrete geologic units and relating them to location. Improved understanding of the regolith thickness and its global variation (GRAIL). Tying the sampling of impact-melt rocks to the lunar impact flux. Using improved techniques (magnetic fields, diffusion studies, isotopic analysis) on existing samples. New sample return from benchmark craters, particularly SPA, which appears in 2013 Decadal Survey.

Cohen, Barbara↗
DOE OSTI2025

Temporal Forecasting of Distributed Temperature Sensing in a Thermal Hydraulic System With Machine Learning and Statistical Models

We benchmark performance of long-short term memory (LSTM) network machine learning model and autoregressive integrated moving average (ARIMA) statistical model in temporal forecasting of distributed temperature sensing (DTS). Data in this study consists of fluid temperature transient measured with two co-located Rayleigh scattering fiber optic sensors (FOS) in a forced convection mixing zone of a thermal tee. We treat each gauge of a FOS as an independent temperature sensor. We first study prediction of DTS time series using Vanilla LSTM and ARIMA models trained on prior history of the same FOS that is used for testing. The results yield maximum absolute percentage error (MaxAPE) and root mean squared percentage error (RMSPE) of 1.58% and 0.06% for ARIMA, and 3.14% and 0.44% for LSTM, respectively. Next, we investigate zero-shot forecasting (ZSF) with LSTM and ARIMA trained on history of the co-located FOS only, which is advantageous when limited training data is available. The ZSF MaxAPE and RMSPE values for ARIMA are comparable to those of the Vanilla use case, while the error values for LSTM increase. We show that in ZSF, performance of LSTM network can be improved by training on most correlated gauges between the two FOS, which are identified by calculating the Pearson correlation coefficient. The improved ZSF MaxAPE and RMSPE for LSTM are 4.4% and 0.33%, respectively. Performance of ZSF LSTM can be further enhanced through transfer learning (TL), where LSTM is re-trained on a subset of the FOS that is the target of forecasting. We show that LSTM pre-trained on correlated dataset and re-trained on 30% of testing target dataset achieves MaxAPE and RMSPE values of 2.32% and 0.28%, respectively.

ARIMA↗
DOE OSTI2023

Benchmarks and results of the two-band Hubbard model from the Gutzwiller conjugate gradient minimization theory

Ground-state properties, such as energies and double occupancies, of a one-dimensional two-band Hubbard model are calculated using a first-principles Gutzwiller conjugate gradient minimization theory. The favorable agreement with the results from the density matrix renormalization group theory demonstrates the accuracy of our method. A rotationally invariant approach is further incorporated into the method to greatly reduce the computational complexity with a speedup of approximately 50 times. Moreover, we investigate the Mott transition between a metal and a Mott insulator by evaluating the charge gap. In conclusion, with greatly reduced computational effort, our method reproduces the phase diagram in reasonable agreement with the density matrix renormalization group theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effects of error, chimera, bias, and GC content on the accuracy of amplicon sequencing

ABSTRACT Targeted amplicon sequencing is widely used in microbial ecology studies. However, sequencing artifacts and amplification biases are of great concern. To identify sources of these artifacts, a systematic analysis was performed using mock communities comprised of 16S rRNA genes from 33 bacterial strains. Our results indicated that while sequencing errors were generally isolated to low-abundance operational taxonomic units, chimeric sequences were a major source of artifacts. Singleton and doubleton sequences were primarily chimeras. Formation of chimeric sequences was significantly correlated with the GC content of the targeted sequences. Low-GC-content mock community members exhibited lower rates of chimeric sequence formation. GC content also had a large impact on sequence recovery. The quantitative capacity was notably limited, with substantial recovery variations and weak correlation between anticipated and observed strain abundances. The mock community strains with higher GC content had higher recovery rates than strains with lower GC content. Amplification bias was also observed due to the differences in primer affinity. A two-step PCR strategy reduced the number of chimeric sequences by half. In addition, comparative analyses based on the mock communities showed that several widely used sequence processing pipelines/methods, including DADA2, Deblur, UCLUST, UNOISE, and UPARSE, had different advantages and disadvantages in artifact removal and rare species detection. These results are important for improving sequencing quality and reliability and developing new algorithms to process targeted amplicon sequences. IMPORTANCE Amplicon sequencing of targeted genes is the predominant approach to estimate the membership and structure of microbial communities. However, accurate reconstruction of community composition is difficult due to sequencing errors, and other methodological biases and effective approaches to overcome these challenges are essential. Using a mock community of 33 phylogenetically diverse strains, this study evaluated the effect of GC content on sequencing results and tested different approaches to improve overall sequencing accuracy while characterizing the pros and cons of popular amplicon sequence data processing approaches. The sequencing results from this study can serve as a benchmarking data set for future algorithmic improvements. Furthermore, the new insights on sequencing error, chimera formation, and GC bias from this study will help enhance the quality of amplicon sequencing studies and support the development of new data analysis approaches.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Improvements of Nuclear Data Evaluations for Lead Isotopes in Support of Next Generation Lead-Cooled Fast Systems

The neutron evaluation of the isotopes that comprise natural lead were undertaken as a part of DOE-NEUP Project #19-16739. The goal of the project was to update the neutron cross sections to account for new differential measurements and incorporate the most up-to-date physics. Shortcomings in the lead cross sections was made known by several independent reports. The work performed here repeated the simulation of all the “benchmark” validation systems and concluded that the major issue is the scattering cross sections in the major lead isotopes above 100 keV. Re-evaluation of 206,207,208 Pb included both the resolved resonance region and fast region evaluations of the cross sections. Combined these regions cover energies from thermal to 20 MeV. The most drastic improvement is the resolved resonance region evaluation of 208 Pb which is now extended from 1.0 to 1.5 MeV. Parameterization of these resonances in the R-matrix provides a superior reconstruction of not only the experimental cross section but also the scattering distributions via the Blatt-Biedenharn formalism. Fast region evaluations of the three major isotopes were done to include new inelastic experimental data from the neutron Time-of-Flight facility at CERN. The culmination of all the changes to the cross section is a drastic improvement in the scattering kernel as shown in Rensselaer Polytechnic Institute (RPI) Quasi-Differential scattering measurements and improved prediction of keff for fast integral experiments. Alongside the new cross sections, new nuclear data covariance (uncertainties) have been computed and are included in the evaluation. Little is changed in the magnitude of the uncertainties but the correlations within the covariance display non-trivial changes. The new evaluations of 206,207,208 Pb have been submitted to the National Nuclear Data Center to be included in the ENDF/B-VIII.1 library. While the cross sections have been updated extensively, knowledge and modeling of the cross sections between 1.0 and 3.0 MeV for the isotopes remain a challenge. Most notably is the double differential elastic cross section for 208 Pb and 206 Pb. New quasi-differential measurements for pure samples of these two nuclei would go a long way in reducing compensating errors in the evaluations. Qualitatively, the project has produced a prosperous collaboration between the nuclear data group at RPI with staff scientists at Brookhaven National Laboratory, Naval Nuclear Laboratory, Oak Ridge National Laboratory, Los Alamos National Laboratory, Lawrence Livermore National Laboratory, and Sandia National Laboratories. Quantitatively this is reflected in seven conference presentations, at least one journal submission, and one doctoral thesis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
NASA NTRS2023

Characterization of Oscillation Modes in Levitated Droplets Using Image and Non-Image Based Techniques

The dynamics of levitated liquid droplets can be used to measure their thermophysical properties by correlating the frequencies at which normal modes of oscillation most strongly resonate when subject to an external oscillatory force. In two preliminary works, it was shown via electrostatic levitation and processing of various metals and alloys that 1) the resonance of the first principal mode of oscillation (mode n = 2) can be used to accurately measure surface tension and 2) that so-called ‘higher-order resonance’ of n = 3 is observable at a predictable frequency. It was also shown, in the context of future space-based experimentation on the Electrostatic Levitation Furnace (ELF), a setup on the International Space Station (ISS) operated by Japan Aerospace Exploration Agency (JAXA), that while the shadow array method in which droplet behavior is visualized would be challenging to identify the n = 3 resonance, the normal mode n = 4 was predicted to be more easily identifiable. In this short communication, experimental evidence of the first three principal modes of oscillation is provided using molten samples of Tin and Indium and it is subsequently shown that, as predicted, an “image-less” approach can be used to identify both n = 2 and n = 4 resonances in levitated liquid droplets. This suggests that the shadow array method may be satisfactorily used to obtain a self-consistent benchmark of thermophysical properties by comparing results from two successive even-mode natural frequencies.

Nevin Brosius↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2021

Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in i-band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established redMaGiC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016)MNRAA40035-871110.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies i<4zphot+18 and has ∼30% wider redshift distributions but ∼3.5 times more galaxies than redMaGiC. Assuming a wCDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of Ωm, w, and σ8, and gains of 16% in σ8, 10% in Ωm, and 12% in w. Similarly, in ΛCDM, we find an improvement of 19% and 27% on σ8 and Ωm, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i<22.2, has ∼7 times more galaxies and ∼20% wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Minimizing Cr-Evaporation from Balance Of Plant Components By Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al 2 O 3 -forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Minimizing Cr-Evaporation from Balance of Plant Components by Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al2O3-forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

36 MATERIALS SCIENCE↗