Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Correlation coefficient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗
DOE OSTI2022

Mechanical and tribological properties of anodic Al coatings as a function of anodizing conditions

Enhanced wear and corrosion resistance properties are vital for numerous functional and decorative applications. Anodizing is a process frequently used to generate stable oxide coatings on metal surfaces, including aluminium (Al). In this work, we present an evaluation of process variables for optimized anodized coating properties using existing methods and as a baseline for novel coating development. A representative relatively pure Al alloy was chosen as a model system with a goal of finding conditions to achieve high wear coatings with retained ductility. AA5052 alloy coupons were anodized in two different electrolytes (sulfuric and phosphoric acid) to study the effects of anodizing conditions on the mechanical and tribological properties of the anodic coatings. Porosity was evaluated with microscopy, coating hardness was determined via nanoindentation, and tribological properties were characterized using a pin-on-disk tribometer. These measurements were correlated and a loading evaluation was performed to assess the optimum nanoindentation parameters for the use of representative coupons without the need for special processing prior to anodizing. The porosity and coating thickness were found to depend on the anodizing time and the type of electrolyte. Thicker and harder anodic coatings were generated in sulfuric acid compared to phosphoric acid. Porosity studies of the anodic coatings generated in phosphoric acid revealed that coating porosity decreased as the anodizing time increased. Both electrolytes generated anodic coatings with improved tribological properties (friction coefficient and dimensional wear rate), and an optimum processing time was identified for best wear properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dark Energy Survey Year 3 results: magnification modelling and impact on cosmological constraints from galaxy clustering and galaxy–galaxy lensing

We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy–galaxy lensing, using two different lens samples: a sample of luminous red galaxies, redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects using the Balrog image simulations. We estimate the impact of magnification on the galaxy clustering and galaxy–galaxy lensing cosmology analysis, finding it to be a significant systematic for the MagLim sample. We show cosmological constraints from the galaxy clustering autocorrelation and galaxy–galaxy lensing signal with different magnifications priors, finding broad consistency in cosmological parameters in ΛCDM and wCDM. However, when magnification bias amplitude is allowed to be free, we find the two-point correlation functions prefer a different amplitude to the fiducial input derived from the image simulations. We validate the magnification analysis by comparing the cross-clustering between lens bins with the prediction from the baseline analysis, which uses only the autocorrelation of the lens bins, indicating that systematics other than magnification may be the cause of the discrepancy. We show that adding the cross-clustering between lens redshift bins to the fit significantly improves the constraints on lens magnification parameters and allows uninformative priors to be used on magnification coefficients, without any loss of constraining power or prior volume concerns.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Low-temperature spin Seebeck effect in nonmagnetic vanadium dioxide

The spin Seebeck effect (SSE) is sensitive to thermally driven magnetic excitations in magnetic insulators. Vanadium dioxide in its insulating low-temperature phase is expected to lack magnetic degrees of freedom, as vanadium atoms are thought to form singlets upon dimerization of the vanadium chains. Instead, we find a paramagnetic SSE response in V⁢O 2 films that grows as the temperature decreases below 50 K. The field and temperature-dependent SSE voltage is qualitatively consistent with a general model of paramagnetic SSE response and inconsistent with triplet spin transport. Quantitative estimates find a spin Seebeck coefficient comparable in magnitude to that observed in strongly magnetic materials. Finally, the microscopic nature of the magnetic excitations in V⁢O 2 requires further examination.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗
DOE OSTI2021

A semiparametric latent factor model for large scale temporal data with heteroscedasticity

Large scale temporal data have flourished in a vast array of applications, and their sophisticated structures, especially the heteroscedasticity among subjects with inter- and intra-temporal dependence, have fueled a great demand for new statistical models. In this paper, with covariate information, we consider a flexible model for large scale temporal data with subject-specific heteroscedasticity. Formally, the model employs latent semiparametric factors to simultaneously account for the subject-specific heteroscedasticity and the contemporaneous and/or serial correlations. The subject-specific heteroscedasticity is modeled as the product of the unobserved factor process and subject’s covariate effect, which is further characterized via additive models. For estimation, we propose a two-step procedure. First, the latent factor process and nonparametric loading are recovered through projection-based methods, and following, we estimate the regression components by approaches motivated from the generalized least squares. By scrupulously examining the non-asymptotic rates for recovering the factor process and its loading, we show the consistency and efficiency of estimated regression coefficients in the absence of prior knowledge of latent factor process and subject’s covariate effect. Here, the statistical guarantees remain valid even for finite time points that makes our method particularly appealing when the subjects significantly outnumber the observation time points. Using comprehensive simulations, we demonstrate the finite sample performance of our method, which corroborates the theoretical findings. Finally, we apply our method to a data set of air quality and energy consumption collected at 129 monitoring sites in the United States in 2015.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Quantitative measurement of pebble fuel anisotropy using generalized ellipsometry

The renewed development of high-temperature gas-cooled reactors has created a need to modernize pebble fuel quality control methods. One important specification in these methods is the coefficient of thermal expansion (CTE) anisotropy, which reflects graphite grain alignment and affects fuel performance. This study investigates the feasibility of using two-modulator generalized ellipsometry microscopy (2-MGEM) as a potential surrogate metric for CTE anisotropy. Graphite matrix compacts with controlled levels of anisotropy were fabricated by varying pre-pressures during forming and were subsequently characterized using both optical dilatometry and 2-MGEM. Measured responses were compared to an isotropic reference pebble produced using an industrial manufacturing process to provide a practical baseline. The resulting measurements revealed a strong quantitative correlation between the CTE anisotropy ratio and the corresponding diattenuation ratio. However, propagated measurement uncertainties in the optical method were found to limit definitive pass/fail discrimination, highlighting the need for improved precision before qualification as a quality control tool. Still, the results establish a quantitative relationship between optical and CTE anisotropy and support further development of 2-MGEM as a potential surrogate method for evaluating graphite matrix anisotropy in pebble fuel manufacturing.

Raftery, Alicia [ORNL] (ORCID:0000000259076427)↗
DOE OSTI2025

Order-by-order uncertainties of nucleon-nucleon Wolfenstein amplitudes in chiral effective field theory

Quantum mechanical invariance principles dictate the most general operator structure that can be present in the nucleon-nucleon (NN) interaction. Five independent operators appear in the on-shell NN amplitude together with five corresponding coefficient functions. The usual choice for these coefficient functions is known as the NN Wolfenstein amplitudes. We analyze the order-by-order convergence of each of the five NN Wolfenstein amplitudes predicted by a semilocal coordinate space potential implementation of chiral effective field theory (𝜒⁢EFT). We do this at laboratory kinetic energies between 25 and 200 MeV for both neutron-proton and proton-proton scattering. Our analysis uses the Gaussian-process methods developed by the BUQEYE Collaboration to describe the contributions of each 𝜒⁢EFT order, and so yields truncation uncertainties for each Wolfenstein amplitude that are correlated across scattering angles. We combine information on the size of different orders in the EFT to infer the 𝜒⁢EFT breakdown scale for each amplitude, finding, on average, Λ 𝑏 between 750 and 800 MeV. Furthermore, with this choice of Λ 𝑏 , the EFT truncation uncertainties cover both higher-order results and empirical Wolfenstein amplitudes well for all orders other than the leading order.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
DOE OSTI2020

Sulfur isotope fractionation between aqueous and carbonate-associated sulfate in abiotic calcite and aragonite

Sulfate (SO 4 2– ) incorporated into calcium carbonate minerals enables measurements of sulfur (S) isotope ratios in carbonate rocks. This Carbonate Associated Sulfate (CAS) in marine carbonate minerals is thought to faithfully represent the S isotope composition of the seawater sulfate incorporated into the mineral, with little or no S isotope fractionation in the process. However, comparison between different calcifying species reveals both positive and negative S isotope fractionation between CAS and seawater sulfate, and a large range of S isotope ratios can be found within a single rock sample, depending on the component measured. To better understand the isotopic effects associated with sulfate incorporation into carbonate minerals, we precipitated inorganic calcite and aragonite over a range covering more than two orders of magnitude of sulfate concentration and precipitation rate. Coupled measurements of CAS concentration, S isotope composition and X-ray absorption near-edge spectra (XANES) permit characterization and explanation of the observed dependence of S isotope fractionation between CAS and aqueous sulfate (CAS-SO 4 2– isotope fractionation) on sulfate concentration and precipitation rate. In aragonite, the CAS-SO 4 2– isotope fractionation is 1.0 ± 0.3‰ and independent of the sulfate (and CAS) concentration. In contrast, we note the CAS-SO 4 2– isotope fractionation in calcite covaries strongly with the sulfate concentration and weakly with the precipitation rate, between values of 1.3 ± 0.1 and 3.1 ± 0.6‰. We suggest that the correlation between aqueous sulfate concentration and CAS-SO 4 2– isotope fractionation in calcite reflects a dependence of the equilibrium S isotope fractionation on the concentration of CAS, through the effect of the sulfate impurity on the carbonate mineral’s energetic state.

58 GEOSCIENCES↗
DOE OSTI2025

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗
DOE OSTI2024

Study of Magnetic Field and Turbulence in the TeV Halo around the Monogem Pulsar

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants and diffuse halos around pulsars. Its turbulent characteristics govern the diffusion of cosmic rays and the multiwavelength emission from pulsar wind nebulae (PWNe). However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this study, we present the first direct observational evidence of the orientation of the mean magnetic field and turbulent characteristics by employing a recently developed statistical parameter "Y turb " in the extended halo around the Monogem pulsar. Our study points to two possible scenarios: nearly aligned toward the line of sight (LOS) with compressible modes dominance or high inclination angle toward the LOS and characterized by Alfvénic turbulence. The first scenario appears consistent with other observational signatures. Furthermore, we report that the magnetic field has an observed correlation length of approximately 3 ± 0.6 pc in the Monogem halo. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Forward and Off-Forward Parton Distributions from Lattice QCD

The interpretation of (semi-)inclusive and certain exclusive scattering processes relies on the factorization of hard parton level cross sections from long-range and non-perturbative parton correlations. The familiar Parton Distribution Functions (PDFs) and Generalized Parton Distributions quantify the non-perturbative dynamics in these situations and address a number of key questions surrounding the structure of hadrons. A certain class of matrix elements accessible in lattice QCD, so called Lattice Cross Sections, have been shown to factorize into these collinear distributions in a manner akin to the factorization of hadronic cross sections. In the short-distance regime, matrix elements of space-like separated two-current operators and parton bilinears can be expressed as the convolution of perturbative coefficient functions and the PDFs. Matrix elements of this type are isolated in the pion and nucleon, each offering a glimpse of the unpolarized valence quark content of these phenomenologically important hadronic states. The calculations within the nucleon represent the first application of the distillation spatial smearing paradigm to the collinear structure of hadrons, and is found to offer higher precision data compared to similar calculations in the literature. A novel method to obtain PDFs from these lattice data, while simultaneously controlling systematic effects, is developed and applied to the nucleon dataset. The coordinate space factorization of space-like separated parton bilinears has also recently been extended to include Generalized Parton Distributions. Preliminary results in off-forward nucleon matrix elements using distillation are explored.

Egerer, Colin↗
DOE OSTI2024

Correlation of Fabrication Methods and Enhanced Wear Performance in Nanoporous Anodic Aluminum Oxide with Incorporated Molybdenum Disulfide (MoS 2 ) Nanomaterials

Wear performance is integral to component longevity, minimizing industrial waste and excess energy costs in a wide variety of applications. Anodized aluminum oxide (AAO) has many beneficial properties leading to its wide use across industries as a surface treatment for many aluminum components, but the wear properties of the coating could be improved significantly. Here, we used an electrochemical method to incorporate molybdenum disulfide (MoS 2 ), a nanomaterial used as a dry lubricant, to modify alloys of aluminum during AAO preparation. Using Raman spectroscopy and tribological scratch measurements, we thoroughly characterized the structure and wear behavior of the films. The MoS 2 deposition procedure was optimal on aluminum 5052 anodized in higher acid concentrations, with friction coefficients at around 0.05 (~10× better than unmodified AAO). Changing anodization conditions to produce harder films with smaller pores led to worsened wear properties, likely because of lower MoS 2 content. Studying a commercial MoS 2 /AAO film of a different Al alloy (7075) showed that a heat treatment step intended to fully convert all deposited MoS x species to MoS 2 can adversely affect wear in some alloys. While Al 6061 and 1100 produced films with worse wear performance compared to Al 5052 or 7075, our results show evidence that acid cleaning after initial anodization likely removes residual alloying elements, affecting MoS 2 incorporation. This study demonstrates a nanomaterial modified AAO film with superior wear characteristics to unmodified AAO and relates fabrication procedure, film structure, and practical performance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A formalism for extracting track functions from jet measurements

The continued success of the jet substructure program will require widespread use of tracking information to enable increasingly precise measurements of a broader class of observables. The recent reformulation of jet substructure in terms of energy correlators has simplified the incorporation of universal non-perturbative matrix elements, so called “track functions”, in jet substructure calculations. These advances make it timely to understand how these universal non-perturbative functions can be extracted from hadron collider data, which is complicated by the use jet algorithms. In this paper we introduce a new class of jet functions, which we call (semi-inclusive) track jet functions, which describe measurements of the track energy fraction in identified jets. These track jet functions can be matched onto the universal track functions, with perturbatively calculable matching coefficients that incorporate the jet algorithm dependence. We perform this matching, and present phenomenological results for the charged energy fraction in jets at the LHC and EIC/HERA at collinear next-to-leading logarithmic accuracy. We show that higher moments of the charged energy fraction directly exhibit non-linear Lorentzian renormalization group flows, allowing the study of these flows with collider data. Our factorization theorem enables the extraction of universal track functions from jet measurements, opening the door to their use for a precision jet substructure program.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Impact of new experimental data on the C2HDM: the strong interdependence between LHC Higgs data and the electron EDM

The complex two-Higgs doublet model (C2HDM) is one of the simplest extensions of the Standard Model with a source of CP-violation in the scalar sector. It has a $\mathbb{Z}$ 2 symmetry, softly broken by a complex coefficient. There are four ways to implement this symmetry in the fermion sector, leading to models known as Type-I, Type-II, Lepton Specific and Flipped. In the latter three models, there is a priori the surprising possibility that the 125 GeV Higgs boson couples mostly as a scalar to top quarks, while it couples mostly as a pseudoscalar to bottom quarks. This “maximal” scenario was still possible with the data available in 2017. Since then, there have been more data on the 125 GeV Higgs boson, direct searches for CP-violation in angular correlations of $τ$-leptons produced in Higgs boson decays, new results on the electron electric dipole moment, new constraints from LHC searches for additional Higgs bosons and new results on $b$ → $sγ$ transitions. Highlighting the crucial importance of the physics results of LHC’s Run 2, we combine all these experiments and show that the “maximal” scenario is now excluded in all models. Still, one can have a pseudoscalar component in $hτ\overline{τ}$ couplings in the Lepton-Specific case as large as 87% of the scalar component for all mass orderings of the neutral scalar bosons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Medium-Scale Methanol Pool Fire Model Validation

In this work, medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, non-sooting flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗