Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Yield function”

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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Transient Multiphysics Simulations with Pin Power Reconstruction in the Griffin Reactor Physics Code

This work introduces the pin power reconstruction capability available in the Griffin reactor physics code. This capability is implemented in an unstructured mesh framework, and the methods introduced are applied to the 2D SIMBA reactor core, which has assemblies and pins arranged in a hexagonal lattice. Since this reactor has a non-Cartesian geometry and also operates in the thermal spectrum, a general approach to pin power reconstruction is adopted, where SPH-based equivalence is leveraged to preserve assembly-wise reaction rates, while computing full-core form functions to preserve pin-wise fission production rates within the fuel pins of the reactor core. In a 2D microreactor benchmark problem, this pin power reconstruction approach was shown to reproduce pin powers compared to the Serpent2 Monte Carlo code for fixed temperature conditions and control drum rotation angles, yielding a core-wide RMS error level of 0.6\% and a maximum absolute pin error of 2.3\%. In addition, a tabulated library of multigroup cross sections, SPH factors, and form functions was generated to demonstrate the applicability of pin power reconstruction to a thermal feedback problem. Finally, a control drum transient was successfully simulated, showcasing the application of pin power reconstruction in a transient multiphysics feedback problem.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fission fragment decay simulations with the CGMF code

The CGMF code implements the Hauser-Feshbach statistical nuclear reaction model to follow the de-excitation of fission fragments by successive emissions of prompt neutrons and γ rays. The Monte Carlo technique is used to facilitate the analysis of complex distributions and correlations among the prompt fission observables. Starting from initial configurations for the fission fragments in mass, charge, kinetic energy, excitation energy, spin, and parity, Y(A,Z,KE,U,J,π), CGMF samples neutron and γ-ray probability distributions at each stage of the decay process, conserving energy, spin and parity. Nuclear structure and reaction input data from the RIPL3 library are used to describe fission fragment properties and decay probabilities. Characteristics of prompt fission neutrons, prompt fission γ rays, and independent fission yields can be studied consistently. Correlations in energy, angle and multiplicity among the emitted neutrons and γ rays can be easily analyzed as a function of the emitting fragments. Finally, in its current version, CGMF is limited to certain spontaneous fission reactions on Pu and Cf isotopes, as well as neutron-induced fission reactions up to 20 MeV on U, Np, and Pu isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Laboratory transferability of optimally shaped laser pulses for quantum control

Optimal control experiments can readily identify effective shaped laser pulses, or “photonic reagents,” that achieve a wide variety of objectives. An important additional practical desire is for photonic reagent prescriptions to produce good, if not optimal, objective yields when transferred to a different system or laboratory. Building on general experience in chemistry, the hope is that transferred photonic reagent prescriptions may remain functional even though all features of a shaped pulse profile at the sample typically cannot be reproduced exactly. As a specific example, we assess the potential for transferring optimal photonic reagents for the objective of optimizing a ratio of photoproduct ions from a family of halomethanes through three related experiments. First, applying the same set of photonic reagents with systematically varying second- and third-order chirp on both laser systems generated similar shapes of the associated control landscape (i.e., relation between the objective yield and the variables describing the photonic reagents). Second, optimal photonic reagents obtained from the first laser system were found to still produce near optimal yields on the second laser system. Third, transferring a collection of photonic reagents optimized on the first laser system to the second laser system reproduced systematic trends in photoproduct yields upon interaction with the homologous chemical family. These three transfers of photonic reagents are demonstrated to be successful upon paying reasonable attention to overall laser system characteristics. The ability to transfer photonic reagents from one laser system to another is analogous to well-established utilitarian operating procedures with traditional chemical reagents. The practical implications of the present results for experimental quantum control are discussed.

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

Detecting element specific electrons from a single cobalt nanocluster with synchrotron x-ray scanning tunneling microscopy

Here, we use a nanofabricated scanning tunneling microscope tip as a detector to investigate local X-ray induced tunneling and electron emission from a single cobalt nanocluster on a Au(111) surface. The tip-detector is positioned a few angstroms above the nanocluster, and ramping the incident X-ray energy across the Co photoabsorption K-edge enables the detection of element specific electrons. Atomic-scale spatial dependent changes in the X-ray absorption cross section are directly measured by taking the X-ray induced current as a function of X-ray energy. From the measured sample and tip currents, element specific X-ray induced current components can be separated and thereby the corresponding yields for the X-ray induced processes of the single cobalt nanocluster can be determined. The detection of element specific synchrotron X-ray induced electrons of a single nanocluster opens a new avenue for materials characterization on a one particle at-a-time basis.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Toward Molecular Catalysts by Computer

Rational design of molecular catalysts requires a systematic approach to designing ligands with specific functionality and precisely tailored electronic and steric properties. It then becomes possible to devise computer protocols to predict accurately the required properties and ultimately to design catalysts by computer. In this account we first review how thermodynamic properties such as oxidation-reduction potentials (E0), acidities (pKa), and hydride donor abilities (ΔGH-) form the basis for a systematic design of molecular catalysts for reactions that are critical for a secure energy future (hydrogen evolution and oxidation, oxygen and nitrogen reduction, and carbon dioxide reduction). We highlight how density functional theory allows us to determine and predict these properties within “chemical” accuracy (~ 0.06 eV for redox potentials, ~ 1 pKa unit for pKa values, and ~ 1.5 kcal/mol for hydricities). These quantities determine free energy maps and profiles associated with catalytic cycles, i.e. the relative energies of intermediates, and help us distinguish between desirable and high-energy pathways and mechanisms. Good catalysts have flat profiles that avoid high activation barriers due to low and high energy intermediates. We illustrate how the criterion of a flat energy profile lends itself to the prediction of design points by computer for optimum catalysts. This research was carried out in the Center for Molecular Electro-catalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences. Pacific Northwest National Laboratory (PNNL) is operated for the DOE by Battelle.

acidity↗
DOE OSTI2025

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗
DOE OSTI2013

Comparison of sugar content for ionic liquid pretreated Douglas-fir woodchips and forestry residues

The development of affordable woody biomass feedstocks represents a significant opportunity in the development of cellulosic biofuels. Primary woodchips produced by forest mills are considered an ideal feedstock, but the prices they command on the market are currently too expensive for biorefineries. In comparison, forestry residues represent a potential low-cost input but are considered a more challenging feedstock for sugar production due to complexities in composition and potential contamination arising from soil that may be present. We compare the sugar yields, changes in composition in Douglas-fir woodchips and forestry residues after pretreatment using ionic liquids and enzymatic saccharification in order to determine if this approach can efficiently liberate fermentable sugars. : These samples were either mechanically milled through a 2 mm mesh or pretreated as received with the ionic liquid (IL) 1-ethyl-3-methylimidazolium acetate [C 2 mim][OAc] at 120°C and 160°C. IL pretreatment of Douglas-fir woodchips and forestry residues resulted in approximately 71-92% glucose yields after enzymatic saccharification. X-ray diffraction (XRD) showed that the pretreated cellulose was less crystalline after IL pretreatment as compared to untreated control samples. Two-dimensional nuclear magnetic resonance spectroscopy (2D-NMR) revealed changes in lignin and hemicellulose structure and composition as a function of pretreatment. Mass balances of sugar and lignin streams for both the Douglas-fir woodchips and forestry residues throughout the pretreatment and enzymatic saccharification processes are presented. While the highest sugar yields were observed with the Douglas-fir woodchips, reasonably high sugar yields were obtained from forestry residues after ionic liquid pretreatment. Structural changes to lignin, cellulose and hemicellulose in the woodchips and forestry residues of Douglas-fir after [C 2 mim][OAc] pretreatment are analyzed by XRD and 2D-NMR, and indicate that significant changes occurred. Irrespective of the particle sizes used in this study, ionic liquid pretreatment successfully allowed high glucose yields after enzymatic saccharification. These results indicate that forestry residues may be a more viable feedstock than previously thought for the production of biofuels.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

KEPLER-LIKE MULTI-PLEXING FOR MASS PRODUCTION OF MICROLENS PARALLAXES

We show that a wide-field Kepler-like satellite in solar orbit could obtain microlens parallaxes for several thousand events per year that are identified from the ground, yielding masses and distances for several dozen planetary events. This is roughly an order of magnitude larger than previously considered narrow-angle designs. Such a satellite would, in addition, roughly double the number of planet detections (and mass/distance determinations). It would also yield a trove of brown-dwarf binaries with masses, distances, and (frequently) full orbits, enable new probes of the stellar mass function, and identify isolated black-hole candidates. We show that the actual Kepler satellite, even with degraded pointing, can demonstrate these capabilities and make substantial initial inroads into the science potential. We discuss several ''Deltas'' to the Kepler satellite aimed at optimizing microlens parallax capabilities. Most of these would reduce costs. The wide-angle approach advocated here has only recently become superior to the old narrow-angle approach, due to the much larger number of ground-based microlensing events now being discovered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Near-Barrier Fusion of the {sup 8}B+{sup 58}Ni Proton-Halo System

Fusion cross sections were measured for the exotic proton-halo nucleus {sup 8}B incident on a {sup 58}Ni target at several energies near the Coulomb barrier. This is the first experiment to report on the fusion of a proton-halo nucleus. The resulting excitation function shows a striking enhancement with respect to expectations for normal projectiles. Evidence is presented that the sum of the fusion and breakup yields saturates the total reaction cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model

We address the inverse problem of inferring the basal geothermal heat flux from surface velocity observations using a steady-state thermomechanically coupled nonlinear Stokes ice flow model. This is a challenging inverse problem since the map from basal heat flux to surface velocity observables is indirect: the heat flux is a boundary condition for the thermal advection–diffusion equation, which couples to the nonlinear Stokes ice flow equations; together they determine the surface ice flow velocity. This multiphysics inverse problem is formulated as a nonlinear least-squares optimization problem with a cost functional that includes the data misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to render the problem well posed. We derive adjoint-based gradient and Hessian expressions for the resulting partial differential equation (PDE)-constrained optimization problem and propose an inexact Newton method for its solution. As a consequence of the Petrov–Galerkin discretization of the energy equation, we show that discretization and differentiation do not commute; that is, the order in which we discretize the cost functional and differentiate it affects the correctness of the gradient. Using two- and three-dimensional model problems, we study the prospects for and limitations of the inference of the geothermal heat flux field from surface velocity observations. The results show that the reconstruction improves as the noise level in the observations decreases and that short-wavelength variations in the geothermal heat flux are difficult to recover. We analyze the ill-posedness of the inverse problem as a function of the number of observations by examining the spectrum of the Hessian of the cost functional. Motivated by the popularity of operator-split or staggered solvers for forward multiphysics problems – i.e., those that drop two-way coupling terms to yield a one-way coupled forward Jacobian – we study the effect on the inversion of a one-way coupling of the adjoint energy and Stokes equations. Here, we show that taking such a one-way coupled approach for the adjoint equations can lead to an incorrect gradient and premature termination of optimization iterations. This is due to loss of a descent direction stemming from inconsistency of the gradient with the contours of the cost functional. Nevertheless, one may still obtain a reasonable approximate inverse solution particularly if important features of the reconstructed solution emerge early in optimization iterations, before the premature termination.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Compositional phase stability of correlated electron materials within DFT+DMFT

Predicting the compositional phase stability of strongly correlated electron materials is an outstanding challenge in condensed matter physics, requiring precise computations of total energies. In this work, we employ the density functional theory plus dynamical mean-field theory (DFT+DMFT) formalism to address local correlations due to transition metal $\textit{d}$ electrons on compositional phase stability in the prototype rechargeable battery cathode material Li x CoO 2 , and detailed comparisons are made with the simpler DFT+$\textit{U}$ approach (i.e., the Hartree-Fock solution of the DMFT impurity problem). Local interactions are found to strongly impact the energetics of the band insulator LiCoO 2 , most significantly via the $E_g$ orbitals, which are partially occupied via hybridization with $\textit{O p}$ states. We find CoO 2 and Li 1/2 CoO 2 to be moderately correlated Fermi liquids with quasiparticle weights of 0.6–0.8 for the $T_{2g}$ states, which are most impacted by the interactions. Further, as compared to DFT+$\textit{U}$, DFT+DMFT considerably dampens the increase in total energy as $\textit{U}$ is increased, which indicates that dynamical correlations are important to describe this class of materials despite the relatively modest quasiparticle weights. Unlike DFT+$\textit{U}$, which can incorrectly drive Li x CoO 2 toward spurious phase separating or charge-ordered states, DFT+DMFT correctly captures the system's phase stability and does not exhibit a strong charge-ordering tendency. Most importantly, the error within DFT+$\textit{U}$ varies strongly as the composition changes, challenging the common practice of artificially tuning $\textit{U}$ within DFT+$\textit{U}$ to compensate the errors of Hartree-Fock. DFT+DMFT predicts the average intercalation voltage decreases relative to DFT, opposite to the result of DFT+$\textit{U}$, which would yield favorable agreement with experiment in conjunction with the overprediction of the voltage by the strongly constrained and appropriately normed DFT functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Long- and short-range correlations and their event-scale dependence in high-multiplicity pp collisions at $ \sqrt{s} $ = 13 TeV

Two-particle angular correlations are measured in high-multiplicity proton-proton collisions at √s = 13 TeV by the ALICE Collaboration. The yields of particle pairs at short-(Δη ~ 0) and long-range (1.6 < |Δη| < 1.8) in pseudorapidity are extracted on the near-side (Δ~ 0). They are reported as a function of transverse momentum (p T ) in the range 1 < p T < 4 GeV/c. Furthermore, the event-scale dependence is studied for the first time by requiring the presence of high-p T leading particles or jets for varying p T thresholds. The results demonstrate that the long-range “ridge” yield, possibly related to the collective behavior of the system, is present in events with high-p T processes as well. The magnitudes of the short- and long-range yields are found to grow with the event scale. The results are compared to EPOS LHC and PYTHIA 8 calculations, with and without string-shoving interactions. It is found that while both models describe the qualitative trends in the data, calculations from EPOS LHC show a better quantitative agreement for the p T dependency, while overestimating the event-scale dependency.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

THE PANCHROMATIC HUBBLE ANDROMEDA TREASURY. IV. A PROBABILISTIC APPROACH TO INFERRING THE HIGH-MASS STELLAR INITIAL MASS FUNCTION AND OTHER POWER-LAW FUNCTIONS

We present a probabilistic approach for inferring the parameters of the present-day power-law stellar mass function (MF) of a resolved young star cluster. This technique (1) fully exploits the information content of a given data set; (2) can account for observational uncertainties in a straightforward way; (3) assigns meaningful uncertainties to the inferred parameters; (4) avoids the pitfalls associated with binning data; and (5) can be applied to virtually any resolved young cluster, laying the groundwork for a systematic study of the high-mass stellar MF (M {approx}> 1 M {sub Sun }). Using simulated clusters and Markov Chain Monte Carlo sampling of the probability distribution functions, we show that estimates of the MF slope, {alpha}, are unbiased and that the uncertainty, {Delta}{alpha}, depends primarily on the number of observed stars and on the range of stellar masses they span, assuming that the uncertainties on individual masses and the completeness are both well characterized. Using idealized mock data, we compute the theoretical precision, i.e., lower limits, on {alpha}, and provide an analytic approximation for {Delta}{alpha} as a function of the observed number of stars and mass range. Comparison with literature studies shows that {approx}3/4 of quoted uncertainties are smaller than the theoretical lower limit. By correcting these uncertainties to the theoretical lower limits, we find that the literature studies yield ({alpha}) = 2.46, with a 1{sigma} dispersion of 0.35 dex. We verify that it is impossible for a power-law MF to obtain meaningful constraints on the upper mass limit of the initial mass function, beyond the lower bound of the most massive star actually observed. We show that avoiding substantial biases in the MF slope requires (1) including the MF as a prior when deriving individual stellar mass estimates, (2) modeling the uncertainties in the individual stellar masses, and (3) fully characterizing and then explicitly modeling the completeness for stars of a given mass. The precision on MF slope recovery in this paper are lower limits, as we do not explicitly consider all possible sources of uncertainty, including dynamical effects (e.g., mass segregation), unresolved binaries, and non-coeval populations. We briefly discuss how each of these effects can be incorporated into extensions of the present framework. Finally, we emphasize that the technique and lessons learned are applicable to more general problems involving power-law fitting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Radiative {phi}{yields}{sigma}{gamma} Vertex in QCD Sum Rules and Phenomenological Approach

We estimate the coupling constant g{sub {phi}{sigma}{gamma}} by using the three-point correlation function and we also investigate this coupling constant by studying the effective lagrangians. Here we consider the effect the {phi}-{omega} mixing effect to the calculation. Besides, a comparison of our predictions with the theoretical results existing in literature is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Measurement of beauty-strange meson production in Pb–Pb collisions at s NN = 5.02 TeV via non-prompt D s + mesons

The production yields of non-prompt ${D}^{+}_{s}$ mesons, namely ${D}^{+}_{s}$ mesons from beauty-hadron decays, were measured for the first time as a function of the transverse momentum ($p_T$) at midrapidity (|$y$| < 0.5) in central and semi-central Pb–Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s{NN}}$ = 5.02 TeV with the ALICE experiment at the LHC. The ${D}^{+}_{s}$ mesons and their charge conjugates were reconstructed from the hadronic decay channel ${D}^{+}_{s}$ → $φπ^+$, with $φ$ → K - K + , in the 4 < $p_T$ < 36 GeV/c and 2 < $p_T$ < 24 GeV/ $c$ intervals for the 0–10% and 30–50% centrality classes, respectively. The measured yields of non-prompt ${D}^{+}_{s}$ mesons are compared to those of prompt ${D}^{+}_{s}$ and non-prompt $D^0$ mesons by calculating the ratios of the production yields in Pb–Pb collisions and the nuclear modification factor $R_{AA}$. The ratio between the $R_{AA}$ of non-prompt ${D}^{+}_{s}$ prompt ${D}^{+}_{s}$ mesons, and that between the $R_{AA}$ of non-prompt ${D}^{+}_{s}$ and non-prompt $D^0$ mesons in central Pb–Pb collisions are found to be on average higher than unity in the 4 < $p_T$ < 12 GeV/ $c$ interval with a statistical significance of about 1.6σ and 1.7σ, respectively. The measured $R_{AA}$ ratios are compared with the predictions of theoretical models of heavy-quark transport in a hydrodynamically expanding QGP that incorporate hadronisation via quark recombination.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Alkaline earth-based coordination polymers derived from a cyclotriphosphazene-functionalized hexacarboxylate

Combination of hexakis(4-carboxylatephenoxy)cyclotriphosphazene with alkaline earth ions of increasing ionic radii (Mg{sup 2+}, Ca{sup 2+} and Ba{sup 2+}) under different solvothermal conditions yielded three new coordination polymers, and their crystal structures were determined by single-crystal X-ray diffraction analysis. The magnesium compound displays a three dimensional (3D) network structure constructed from the deprotonated ligand and the secondary building block Mg(COO){sub 4}, which can be rationalized as a (4,6)-connected topological net with the Schläfli symbol of (4{sup 4}·6{sup 2}){sub 3}(4{sup 9}·6{sup 6}){sub 2}. The calcium compound consists of 1D infinite “Ca-O” inorganic chains connected by the deprotonated ligands to from a 3D framework. The barium compound exhibits a 3D framework in which 1D “Ba-O” inorganic chains are connected together by the deprotonated organic linkers. Due to the semi-rigid nature, the ligand adopts distinctly different conformations in the three compounds. The metal ions’ influence exerted on the final structure of the resulting coordination polymers is also discussed. When the radii of alkaline earth ions increases descending down the group from Mg(II) to Ba(II), the coordination number becomes larger and more versatile: from 6 in the magnesium compound, to 6,7 and 10 in the calcium compound, and to 8 and 9 in the barium compound, thus substantially influencing the resulting final framework structures. Also, the photophysical properties were investigated systematically, revealing that the three compounds are photoluminscent in the solid state at room temperature. This work demonstrates that although the multiplicity of conformation in the hexacarboxylate ligand based on the inorganic scaffold cyclotriphosphazene makes it difficult to predict how this ligand will form extended network, but provides unique opportunities for the formation of diverse inorganic-organic hybrids exhibiting rich structural topologies. - Graphical abstract: Three alkaline earth-based coordination polymers were constructed from a semirigid cyclotriphosphazene-functionalized hexacarboxylate exhibiting different inorganic nodes as well as distinct ligand conformations relying on the metal ions, which presents the first example of such a ligand incorporated into alkaline earth–based coordination polymers. - Highlights: • Three alkaline earth-based coordination polymers were synthesized. • The three compounds exhibit different inorganic nodes and ligand conformations. • The three compounds are photoluminscent in the solid state.

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

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗