SEARCH · Engineering Papers

Results for “Numerical benchmark”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2026

ENDF/B-VIII.1: Updated Nuclear Reaction Data Library for Science and Applications

The ENDF/B-VIII.1 library is the newest recommended evaluated nuclear data file by the Cross Section Evaluation Working Group (CSEWG) for use in nuclear science and technology applications, and incorporates advances made in the six years since the release of ENDF/B-VIII.0. Among key advances made are that the 239 Pu file was reevaluated by a joint international effort and that updated 16,18 O, 19 F, 28–30 Si, 50–54 Cr, 55 Mn, 54,56,57 Fe, 63,65 Cu, 139 La, 233,235,238 U, and 240,241 Pu neutron nuclear data from the IAEA coordinated INDEN collaboration were adopted. Over 60 neutron dosimetry cross sections were adopted from the IAEA's IRDFF-II library. In addition, the new library includes significant changes for 3 He, 6 Li, 9 Be, 51 V, 88 Sr, 103 Rh, 140,142 Ce, Dy, 181 Ta, Pt, 206–208 Pb, and 234,236 U neutron data, and new nuclear data for the photonuclear, charged-particle and atomic sublibraries. Numerous thermal neutron scattering kernels were reevaluated or provided for the very first time. On the covariance side, work was undertaken to introduce better uncertainty quantification standards and testing for nuclear data covariances. The significant effort to reevaluate important nuclides has reduced bias in the simulations of many integral experiments with particular progress noted for fluorine, copper, and stainless steel containing benchmarks. Data issues hindered the successful deployment of the previous ENDF/B-VIII.0 for commercial nuclear power applications in high burnup situations. These issues were addressed by improving the 238 U and 239,240,241 Pu evaluated data in the resonance region. The new library performance as a function of burnup is similar to the reference ENDF/B-VII.1 library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Agentic Diagrammatica: Towards Autonomous Symbolic Computation in High Energy Physics

We present Diagrammatica, a symbolic computation extension to the HEPTAPOD agentic framework, which enables LLM agents to plan and execute multi-step theoretical calculations. Symbolic computation poses a distinctive reliability challenge for LLM agents, as correctness is governed by implicit mathematical conventions that are not encoded in a form that can be easily checked in the computational backend. We identify two complementary remedies, tool-constrained computation and targeted knowledge grounding, and pursue the first as the primary architecture. Concretely, we concentrate the agent's action distribution onto tool calls with convention-fixing semantics, in which the agent specifies a compact, human-auditable diagram specification and a trusted backend performs the symbolic or numerical manipulations exactly. The toolkit provides two complementary calculation paths consuming a shared diagram specification: Naive Dimensional Analysis (NDA) for order-of-magnitude rate estimates and Exact Diagrammatic Analysis (EDA) for tree-level symbolic calculations via automatic FeynCalc code generation, both supplemented by automatic Feynman diagram enumeration and a navigable theory knowledge base. The architecture is validated on two benchmarks: (1) an exhaustive catalog of all tree-level, single-vertex $1\to 2$ partial decay widths across scalar, fermion, and vector parents, with complete massless and threshold limits and Standard Model validation; and (2) an NDA sensitivity study of the muon decay multiplicity $μ^+ \to ν_μ\barν_e + n(e^+e^-) + e^-$, determining the maximum observable $n$ at current and planned muon experiments.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗

DOE OSTI2026

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗

DOE OSTI2011

COMPARISON OF THE ACCURACY OF VARIOUS SPATIAL DISCRETIZATION SCHEMES OF THE DISCRETE ORDINATES EQUATIONS IN 2D CARTESIAN GEOMETRY

We present a comprehensive error estimation of four spatial discretization schemes of the two-dimensional Discrete Ordinates (SN) equations on Cartesian grids utilizing a Method of Manufactured Solution (MMS) benchmark suite based on variants of Larsen’s benchmark featuring different orders of smoothness of the underlying exact solution. The considered spatial discretization schemes include the arbitrarily high order transport methods of the nodal (AHOTN) and characteristic (AHOTC) types, the discontinuous Galerkin Finite Element method (DGFEM) and the recently proposed higher order diamond difference method (HODD) of spatial expansion orders 0 through 3. While AHOTN and AHOTC rely on approximate analytical solutions of the transport equation within a mesh cell, DGFEM and HODD utilize a polynomial expansion to mimick the angular flux profile across each mesh cell. Intuitively, due to the higher degree of analyticity, we expect AHOTN and AHOTC to feature superior accuracy compared with DGFEM and HODD, but at the price of potentially longer grind times and numerical instabilities. The latter disadvantages can result from the presence of exponential terms evaluated at the cell optical thickness that arise from the semianalytical solution process. This work quantifies the order of accuracy and the magnitude of the error of all four discretization methods for different optical thicknesses, scattering ratios and degrees of smoothness of the underlying exact solutions in order to verify or contradict the aforementioned intuitive expectation.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2016

A flippon related singlet at the LHC II

Here, we consider the 750 GeV diphoton resonance at the 13 TeV LHC in the ℱ-SU(5) model with a Standard Model (SM) singlet field which couples to TeV-scale vector-like particles, dubbed flippons . This singlet field assumes the role of the 750 GeV resonance, with production via gluon fusion and subsequent decay to a diphoton via the vector-like particle loops. We present a numerical analysis showing that the observed 8 TeV and 13 TeV diphoton production cross-sections can be generated in the model space with realistic electric charges and Yukawa couplings for light vector-like masses. We further discuss the experimental viability of light vector-like masses in a General No-Scale ℱ-SU(5) model, offering a few benchmark scenarios in this consistent GUT that can satisfy all experimental constraints imposed by the LHC and other essential experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Probing the electronic structure of graphene near and far from the Fermi level via planar tunneling spectroscopy

Scanning tunneling spectroscopy (STS) has yielded significant insight into the electronic structure of graphene and other two-dimensional (2D) materials. STS directly measures a fundamental and directly calculable quantity: single particle density of states (SPDOS). Due to experimental setup limitations, however, STS has been unable to explore 2D materials in ultrahigh magnetic fields where electron-electron interactions can drastically change the SPDOS. Recent developments in the assembly of heterostructures composed of graphene and hexagonal boron nitride have enabled a device-based alternative to potentially overcome these roadblocks. Thus far, however, these nascent efforts are incomplete in analyzing and understanding tunneling spectra and have yet to explore graphene at high magnetic fields. Here, we report an experiment at magnetic fields up to 18 T, which uses graphene tunneling field effect transistors (TFETs) to establish a clear benchmark for the measurement and analysis of graphene planar tunneling spectroscopy. We acquire gate tunable tunneling spectra of graphene and then use these data and electrostatic arguments to develop a systematic analysis scheme. This analysis reveals that TFET devices directly probe electronic structure features near and far from the Fermi level. In particular, our study yields identification of the Dirac point and numerous Landau levels as they fill and empty with charge via the application of a gate voltage. Our work demonstrates that TFET devices are a viable platform for directly probing the electronic structure of graphene and other 2D materials in high magnetic fields, where exotic electronic states emerge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Iterative quantum optimization of spin glass problems with rapidly oscillating transverse fields

In this work, we introduce a new iterative quantum algorithm, called Iterative Symphonic Tunneling for Satisfiability problems (IST-SAT), which solves quantum spin glass optimization problems using high-frequency oscillating transverse fields. IST-SAT operates as a sequence of iterations, in which bitstrings returned from one iteration are used to set spin-dependent phases in oscillating transverse fields in the next iteration. Over several iterations, the novel mechanism of the algorithm steers the system toward the problem ground state. We benchmark IST-SAT on sets of hard MAX-3-XORSAT problem instances with exact state vector simulation, and report polynomial speedups over Trotterized adiabatic quantum computation and the best known semi-greedy classical algorithm. When IST-SAT is seeded with a sufficiently good initial approximation, the algorithm converges to exact solution(s) in a polynomial number of iterations. Our numerical results identify a critical Hamming radius, or quality of initial approximation, where the time-to-solution crosses from exponential to polynomial scaling in problem size. This work proposes IST-SAT a new quantum algorithm, which improves upon solutions obtained from initial classical or quantum optimization algorithms. The steering mechanism we introduce through IST-SAT presents a new path toward achieving quantum advantage in optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Ultrafast x-ray imaging of coherently controlled molecular dynamics in real space and time

Coherent control aims to manipulate chemical processes on the latent length and timescales of atoms and bonds, scales that are intrinsic to molecular dynamics but not directly resolved by most experimental probes. As a result, existing coherent-control experiments, which overwhelmingly rely on spectroscopic observables, leave a crucial blind spot: the direct, real-space recovery of all atomic and molecular rearrangements in response to coherently controlled excitations. Here, we overcome this limitation by integrating a Tannor-Kosloff-Rice pump-control-probe scheme with ultrafast X-ray scattering to capture snapshots of the wavepacket motion in a benchmark molecular system. We demonstrate this by photoexciting diatomic iodine vapor with a visible pump pulse, selectively steering the wavepacket toward ground-state recombination or dissociative pathways with a time-delayed near-IR control pulse, and recording the dynamics with angstrom and femtosecond precision with an ultrashort hard X-ray probe pulse. By comparing these structural observations with numerical solutions of the time-dependent Schrödinger equation, we reveal how coherent control actively reshapes the molecular charge density distribution. Our results pave the way for leveraging structural feedback as a control handle and provide a fundamental microscopic visualization of quantum decoherence and energy redistribution at the atomic level.

Hopper, Thomas R. [SLAC National Accelerator Labor↗

DOE OSTI2015

Excitonic effects in two-dimensional semiconductors: Path integral Monte Carlo approach

The most striking features of novel two-dimensional semiconductors (e.g., transition metal dichalcogenide monolayers or phosphorene) is a strong Coulomb interaction between charge carriers resulting in large excitonic effects. In particular, this leads to the formation of multicarrier bound states upon photoexcitation (e.g., excitons, trions, and biexcitons), which could remain stable at near-room temperatures and contribute significantly to the optical properties of such materials. In our work we have used the path integral Monte Carlo methodology to numerically study properties of multicarrier bound states in two-dimensional semiconductors. Specifically, we have accurately investigated and tabulated the dependence of single-exciton, trion, and biexciton binding energies on the strength of dielectric screening, including the limiting cases of very strong and very weak screening. Our results of this work are potentially useful in the analysis of experimental data and benchmarking of theoretical and computational models.

2D semiconductor, exciton, trion, biexciton, quant↗

DOE OSTI2012

Winchester/Camberley Homes New Construction Test House Design, Construction, and Short-Term Testing in a Mixed-Humid Climate

The NAHB Research Center partnered with production builder Winchester/Camberley Homes to build a DOE Building America New Construction Test House (NCTH). This single family, detached house, located in the mixed-humid climate zone of Silver Spring, MD, was completed in June 2011. The primary goal for this house was to improve energy efficiency by 30% over the Building America B10 benchmark bydeveloping and implementing an optimized energy solutions package design that could be cost effectively and reliably constructed on a production basis using quality management practices. The intent of this report is to outline the features of this house, discuss the implementation of the energy efficient design, and report on short-term testing results. During the interactive design process ofthis project, numerous iterations of the framing, air sealing, insulation, and space conditioning systems were evaluated for energy performance, cost, and practical implementation. The final design featured numerous advanced framing techniques, high levels of insulation, and the HVAC system entirely within conditioned space. Short-term testing confirmed a very tight thermal envelope andefficient and effective heating and cooling. In addition, relevant heating, cooling, humidity, energy, and wall cavity moisture data will be collected and presented in a future long-term report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2015

Casting Process Modeling and Validation

Computational simulations of gravity casting for metallic U-Pu-Zr nuclear fuel rods were performed with the free surface flow solver Flow3d. Two different casting apparatus designs were explored for studying the effects on the final cast rod profiles. Experimental data for the casting problem obtained from Idaho National Laboratory (INL) was validated with numerical simulations on one of the designs considered. The volume of the rods cast, length of the rods, and temperature profiles were predicted through the transient casting process and were found to match closely with experimental data. A sensitivity analysis with vent locations, pressure, mesh refinement, and effects of surface tension at the interface were explored. This work serves as a validation benchmark for the casting problem and provides guidelines when performing simulation exercises for modeling casting. This paper also provides the mathematical framework to enabling modeling of the casting process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1977

Guide star targeting success for the HEAO-B observatory

The statistics associated with the successful selection and acquisition of guide stars as attitude benchmarks for use in reorientation maneuvers of the HEAO-B observatory are considered as a function of the maneuver angle, initial attitude uncertainties, and the pertinent celestial region. Success likelihoods in excess of 0.99 are predicted assuming anticipated gyro and star tracker error sources. The maneuver technique and guide star selection constraints are described in detail. The results presented are specialized numerically to the HEAO-B observatory. However, the analytical techniques developed are considered applicable to broader classes of spacecraft requiring celestial targeting.

Farrenkopf, R. L.↗

NASA NTRS2003

Computational Hemodynamic Simulation of Human Circulatory System under Altered Gravity

A computational hemodynamics approach is presented to simulate the blood flow through the human circulatory system under altered gravity conditions. Numerical techniques relevant to hemodynamics issues are introduced to non-Newtonian modeling for flow characteristics governed by red blood cells, distensible wall motion due to the heart pulse, and capillary bed modeling for outflow boundary conditions. Gravitational body force terms are added to the Navier-Stokes equations to study the effects of gravity on internal flows. Six-type gravity benchmark problems are originally presented to provide the fundamental understanding of gravitational effects on the human circulatory system. For code validation, computed results are compared with steady and unsteady experimental data for non-Newtonian flows in a carotid bifurcation model and a curved circular tube, respectively. This computational approach is then applied to the blood circulation in the human brain as a target problem. A three-dimensional, idealized Circle of Willis configuration is developed with minor arteries truncated based on anatomical data. Demonstrated is not only the mechanism of the collateral circulation but also the effects of gravity on the distensible wall motion and resultant flow patterns.

Kim. Chang Sung↗

NASA NTRS2010

TPSAS-NF1676L-10819-DND

Additional aerodynamic data for the benchmark configuration of partially-dressed cavity-closed nose landing gear (PDCC-NLG), obtained at a flow Mach number of 0.166, are provided in this presentation. A detailed account of the unsteady surface pressure fluctuations acquired with a Kulite transducer designed to be mobile is also given. The roving nature of the transducer allowed the sensor to be placed at numerous locations where high-amplitude pressure fluctuations are expected. The measured surface pressure spectra are mainly broadband in nature, lacking any local peaks associated with coherent vortex shedding. This finding is in agreement with off-surface flow measurements based on particle image velocimetry, which revealed the flow field to be a collection of separated shear layers without any dominant vortex shedding processes.

Mehdi R Khorrami↗

NASA NTRS2010

TPSAS-NF1676L-10756-DND

Mehdi R Khorrami↗

DOE OSTI2013

Critical node treatment in the analytic function expansion method for Pin Power Reconstruction

Pin Power Reconstruction (PPR) was implemented in PARCS using the eight term analytic function expansion method (AFEN). This method has been demonstrated to be both accurate and efficient. However, similar to all the methods involving analytic functions, such as the analytic node method (ANM) and AFEN for nodal solution, the use of AFEN for PPR also has potential numerical issue with critical nodes. The conventional analytic functions are trigonometric or hyperbolic sine or cosine functions with an angular frequency proportional to buckling. For a critic al node the buckling is zero and the sine functions becomes zero, and the cosine function become unity. In this case, the eight terms of the analytic functions are no longer distinguishable from ea ch other which makes their corresponding coefficients can no longer be determined uniquely. The mode flux distribution of critical node can be linear while the conventional analytic functions can only express a uniform distribution. If there is critical or near critical node in a plane, the reconstructed pin power distribution is often be shown negative or very large values using the conventional method. In this paper, we propose a new method to avoid the numerical problem wit h critical nodes which uses modified trigonometric or hyperbolic sine functions which are the ratio of trigonometric or hyperbolic sine and its angular frequency. If there are no critical or near critical nodes present, the new pin power reconstruction method with modified analytic functions are equivalent to the conventional analytic functions. The new method is demonstrated using the L336C5 benchmark problem. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

Locality of correlation in density functional theory

The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi (TF) approximation in the non-relativistic semiclassical (or large-Z) limit for all matter, i.e., the kinetic energy becomes local. Exchange also becomes local in this limit. Numerical data on the correlation energy of atoms support the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around a large particle number are equivalent to local density approximations and their strong relevance to density functional approximations. Analyzing highly accurate atomic correlation energies, we show that E{sub C} → −A{sub C} ZlnZ + B{sub C}Z as Z → ∞, where Z is the atomic number, A{sub C} is known, and we estimate B{sub C} to be about 37 mhartree. The local density approximation yields A{sub C} exactly, but a very incorrect value for B{sub C}, showing that the local approximation is less relevant for the correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-Z limit generally takes this form, but with B{sub C} a functional of the TF density for the system. The implications for the construction of approximate density functionals are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Scaled boundary cubature scheme for numerical integration over planar regions with affine and curved boundaries

This paper introduces the scaled boundary cubature (SBC) scheme for accurate and efficient integration of functions over polygons and two-dimensional regions bounded by parametric curves. Over two-dimensional domains, the SBC method reduces integration over a region bounded by m curves to integration over m regions (referred to as curved triangular regions), where each region is bounded by two line segments and a curve. With proper (counterclockwise) orientation of the boundary curves, the scheme is applicable to convex and nonconvex domains. Additionally, for star-convex domains, a tensor-product cubature rule with positive weights and integration points in the interior of the domain is obtained. If the integrand is homogeneous, we show that this new method reduces to the homogeneous numerical integration scheme; however, the SBC scheme is more versatile since it is equally applicable to both homogeneous and non-homogeneous functions. This paper also introduces several methods for smoothing integrands with point singularities and near-singularities. When these methods are used, highly efficient integration of weakly singular functions is realized. Finally, the SBC method is applied to a number of benchmark problems, which reveal its broad applicability and superior performance (in terms of time to generate a rule and accuracy per cubature point) when compared to existing methods for integration.

42 ENGINEERING↗