Search NASA⌕ Search

SEARCH · Search NASA

Results for “Polynomials”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Verification of EvaluateFLux Utility Program

The EvaluateFlux program is a post-processing utility program for DIF3D, specifically DIF3D-VARIANT which handles Cartesian and hexagonal geometries. The EvaluateFlux program was developed to allow users to obtain flux and power traverses through the geometry domain, and its initial purpose was to facilitate foil analysis by evaluating the flux solution from DIF3D-VARIANT and combining it with foil cross section data. The EvaluateFlux program can calculate the neutron flux, as well as the reaction rates, at any user provided evaluation point. It does this by identifying the spatial mesh associated with the evaluation point and then evaluates the polynomial based neutron flux moments stored in the NHFLUX file at that point. The output of EvaluateFlux varies depending on the input setup. The maximum output includes the neutron flux and microscopic and macroscopic reaction rates at each evaluation point. The purpose of this work is to verify the outputs of EvaluateFlux. Simple models that have hand calculatable results are first defined and used to verify the EvaluateFlux outputs. More complex cases are then added where a duplicate program of EvaluateFlux that uses PrintTables outputs of the binary files is used to verify the EvaluateFlux outputs. In those complex cases, hand calculations of selected evaluation points were also displayed to confirm the software verification. For all the tests done, the hand calculations agreed well with those calculated by EvaluateFlux. For the larger complex problems, the duplicate program that can process hundreds of evaluation points was able to identify that zero points within some meshes have large errors. This aspect was attributed to the truncation error on the input provided to the duplicate program and is not a concern for the accuracy of the EvaluateFlux software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Software Verification of VARPOW

The VARPOW program is a post-processing utility program for DIF3D, specifically DIF3DVARIANT, and it was developed to provide interface files for the thermal analysis program DASSH. The basic methodology of VARPOW is to take the neutron and gamma flux (moments) calculated by DIF3D (or GAMSOR) and combine them with the heating (coefficient) cross sections to calculate the spatial power distributions using the DIF3DVARIANT spatial basis. VARPOW can use the output from GAMSOR (both steady state neutron and gamma flux calculations) or standard DIF3D/REBUS calculations (neutron flux only). The correct approach for defining the power distribution is to use GAMSOR as its purpose was to properly compute the gamma heating throughout the modeled domain. The power densities calculated by VARPOW are broken into fuel, cladding and coolant terms for which isotope-wise categorization is needed. VARPOW has built in options the user can select for the isotope categorization or VARPOW can import a file that details the isotope categorization. VARPOW can export the solution in the polynomial basis of DIF3D-VARIANT or the monomial basis of DIF3D-VARIANT. The purpose of this work is to verify the power distribution results calculated by VARPOW from both the GAMSOR and DIF3D input options and verify that the input and output options are consistent with the manual. Hand calculation and independent numerical calculation are used for this verification work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The linear thermal expansion of 11 polymers from approximately -100 to +100°C

The linear coefficient of thermal expansion has been measured for 11 polymers from approximately -100 to +100°C. The polymers tested were poly(tetrafluoroethylene) (PTFE), poly(chlorotrifluoroethylene) (PCTFE), three polyethylenes; high density (HDPE), ultra-high molecular weight (UHMWPE) and cross-linked (XPE). The commercial fluoropolymers Kel-F 800 (also known as FK-800®) and THV 500®were tested together with samples of poly(ether ether ketone) (PEEK), poly(methyl methyl methacrylate) (PMMA), polycarbonate (PC) and poly(vinylidene fluoride) (PVDF). Polynomial fits were made to the data over appropriate ranges allowing the expansivity (α) to be estimated.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum Monte Carlo Calculations of Chemical Binding and Reactions

The auxiliary field quantum Monte Carlo method developed by the PIs has been shown to provide the most accurate description of strongly correlated electronic systems, from molecules to solids. Unlike other explicitly many‐body approaches, the quantum Monte Carlo method scales as a low order polynomial of systems size, similar to mean‐field methods such as density functional theory. However, the auxiliary field quantum Monte Carlo algorithm is significantly more expensive than traditional density functional calculations. This creates a bottleneck for applications to extended systems, such as large molecules and solids. One principal objective of this proposal was to develop new auxiliary field quantum Monte Carlo computational strategies to achieve improved scaling with system size, using downfolding and localization schemes, without sacrificing the predictive power of the calculations. A second goal is to extend the reach of auxiliary field quantum Monte Carlo to calculate excited states. This final report summarizes what has been achieved during the course the project toward these goals.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

NEAMS Technical Area Support in MOOSE

The MOOSE framework is a foundational capability used by the NEAMS program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and other new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include enabling selective polynomial basis refinement, implementing a custom convergence system, building a scalable preconditioner for saddle-point problems, and much more.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An Analytical Tool to Evaluate Defect Thermodynamics of (La,Ba)Fe1-xMxO3-δ Perovskites for Solid-Oxide Cell Applications

A modeling tool of the defect thermodynamics of (La,Ba)Fe1-xMxO3-δ perovskites which includes energetic information about oxygen vacancy formation, hydration, hydride formation, and charge disproportionation reactions has been developed1. This tool incorporates defect energies and entropies expressed as sixth order polynomial functions to allow refinements of the defect reaction equilibrium constants in the thermodynamic analysis. Calculation of (La,Ba)Fe1-xMxO3-δ Brouwer diagrams as a function of pO2/pH2O and pH2/H2O in a range of temperatures of interest is facilitated by this modeling tool. The results obtained can provide direct guidance how the electronic and ionic defect concentrations of the triple conducting perovskite materials can be used to optimize performance of solid oxide cells for energy applications. The impact of magnetic and electronic structures of the perovskites on the defect reaction energies and entropies as obtained from density function theory modeling and the role played by hydride defect species will also be discussed.

Lee, Yueh-Lin↗
DOE OSTI2024

Symbolic diagnostics to interpret and analyze neural network models

Embedded machine-learned models (EMLMs) have the promise to improve the predictive accuracy of engineering simulators in environments of national interest. EMLMs often comprise complex input-output maps (e.g., neural networks), which make them unamenable to rigorous analysis and generally difficult to interpret. In the face of decades of theory, this lack of interpretability is a significant barrier to building confidence in these models. This work outlines an approach to interpret EMLMs using sparse polynomial regression for comparison with theoretical understanding. To do so, we build on the concept of Locally Interpretable Model-agnostic Explanations (LIME) using physics-informed clustering, prototype selection, and library construction. While general, we demonstrate our method on tensor-basis neural networks used in Reynolds-Averaged Navier-Stokes simulations of hypersonic fluid flows. Results are presented for a simulated toy model and for direct numerical simulations (DNS) of turbulent flows over a flat plate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Software Verification of VARPOW

The VARPOW program is a post-processing utility program for DIF3D, specifically DIF3DVARIANT, and it was developed to provide interface files for the thermal analysis program DASSH. The basic methodology of VARPOW is to take the neutron and gamma flux (moments) calculated by DIF3D (or GAMSOR) and combine them with the heating (coefficient) cross sections to calculate the spatial power distributions using the DIF3DVARIANT spatial basis. VARPOW can use the output from GAMSOR (both steady state neutron and gamma flux calculations) or standard DIF3D/REBUS calculations (neutron flux only). The correct approach for defining the power distribution is to use GAMSOR as its purpose was to properly compute the gamma heating throughout the modeled domain. The power desnities calculated by VARPOW are broken into fuel, cladding and coolant terms for which isotope-wise categorization is needed. VARPOW has built in options the user can select for the isotope categorization or VARPOW can import a file that details the isotope categorization. VARPOW can export the solution in the polynomial basis of DIF3D-VARIANT or the monomial basis of DIF3D-VARIANT. The purpose of this work is to verify the power distribution results calculated by VARPOW from both the GAMSOR and DIF3D input options and verify that the input and output options are consistent with the manual. Hand calculation and independent numerical calculation are used for this verification work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Verification of the EvaluateFlux Utility Program

The EvaluateFlux program is a post-processing utility program for DIF3D, specifically DIF3D-VARIANT which handles Cartesian and hexagonal geometries. The EvaluateFlux program was developed to allow users to obtain flux and power traverses through the geometry domain, and its initial purpose was to facilitate foil analysis by evaluating the flux solution from DIF3D-VARIANT and combining it with foil cross section data. The EvaluateFlux program can calculate the neutron flux, as well as the reaction rates, at any user provided evaluation point. It does this by identifying the spatial mesh associated with the evaluation point and then evaluates the polynomial based neutron flux moments stored in the NHFLUX file at that point. The output of EvaluateFlux varies depending on the input setup. The maximum output includes the neutron flux and microscopic and macroscopic reaction rates at each evaluation point. The purpose of this work is to verify the outputs of EvaluateFlux. Simple models that have hand calculatable results are first defined and used to verify the EvaluateFlux outputs. More complex cases are then added where a duplicate program of EvaluateFlux that uses PrintTables outputs of the binary files is used to verify the EvaluateFlux outputs. In those complex cases, hand calculations of selected evaluation points were also displayed to confirm the software verification. For all the tests done, the hand calculations agreed well with those calculated by EvaluateFlux. For the larger complex problems, the duplicate program that can process hundreds of evaluation points was able to identify that zero points within some meshes have large errors. This aspect was attributed to the truncation error on the input provided to the duplicate program and is not a concern for the accuracy of the EvaluateFlux software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

On the Computation of the Beam Position in the Arcs of the Hadron Storage Ring

To compensate for the velocity differences of the hadrons and electrons in the Electron-Ion Collider, the hadron beam will be operated in a large range of the radial displacement (±20 mm in the 64 mm wide vacuum chamber). Machine protection and beam optics measurement require high-accuracy calculations of the beam position from the signal induced on the pickup electrodes. For the HSR, we need at least the fifth-order polynomial fit, which provides 70 microns r.m.s. error in the 40 by 20 mm region of interest as shown in Fig. 1. The calculations are made using a MATLAB script solving the 2D electrostatic problem.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Construction of approximate invariants for non-integrable Hamiltonian systems

We present a method to construct high-order polynomial approximate invariants (AI) for non integrable Hamiltonian dynamical systems, and apply it to a modern ring-based particle accelerator. Taking advantage of a special property of one-turn transformation maps in the form of a square matrix, AIs can be constructed order-by-order iteratively. Evaluating AI with simulation data, we observe that AI’s fluctuation is actually a measure of chaos. Through minimizing the fluctuations, the stable region of long-term motions, i.e., the dynamic aperture of the accelerator, could be enlarged.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Reassessment of Boiling Water Reactor Rod Cusping Capability of VERA-MPACT

This technical report assesses the current modeling capabilities of the Virtual Environment for Reactor Applications (VERA) for minimizing control rod cusping in boiling water reactors. Assessment is performed on a 10 × 10 General Electric (GE)-14 assembly as well as a simplified version of the same assembly. Runs were performed using VERA for three available cusping treatments in the Michigan PArallel Characteristics Transport (MPACT) code: no treatment, subplane treatment, and fitted polynomial correction treatment. Subplane treatment is discussed, but specific results are not provided because of failures in the method of implementation for boiling water reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Optimising HEP Simulations on 3D Architecture: Primitive Gates for BT

In this report, the first attempt for an alternative approach for the simulation of the non-abelian group BT has been done. First, we analysed why SNAP and Displacement decomposition is not feasible for this problem with the current technology, and then we explored a new approach for the optimisation of the pulses. In particular, we studied how the expansion of the pulses in Chebyshev Polynomials could affect the optimisation process. Even if the study is not conclusive, we believe that it might be worth continuing to study this approach in the future, in order to be able to efficiently implement quantum circuits on computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Reassessment of Boiling Water Reactor Rod Cusping Capability of VERA-MPACT

This technical report assesses the current modeling capabilities of the Virtual Environment for Reactor Applications (VERA) for minimizing control rod cusping in boiling water reactors. Assessment is performed on a 10 × 10 General Electric (GE)-14 assembly as well as a simplified version of the same assembly. Runs were performed using VERA for three available cusping treatments in the Michigan PArallel Characteristics Transport (MPACT) code: no treatment, subplane treatment, and fitted polynomial correction treatment. Subplane treatment is discussed, but specific results are not provided because of failures in the method of implementation for boiling water reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Assessment of BPM options for the EIC Beam Accumulator Ring

The electron injection system for the Electron-Ion Collider (EIC) at BNL is designed to provide a beam of polarized electrons, which is crucial for studying the structure of protons and atomic nuclei. The Beam Accumulator Ring (BAR) is a part of the injection chain between the 750 MeV linear accelerator and the Rapid Cycling Synchrotron (RCS), which accelerates the beam up to full energy (5–18 GeV). The functional role of the BAR is to accumulate the charge injected from the linear accelerator in order to achieve the high intensity of the polarized electron beam required by the specifications for injection into the RCS. This objective will be realized through the sequential injection of bunches with a charge of 1.1 nC at a repetition rate of 30 Hz. Once the charge of a single bunch reaches 28 nC, the beam will be extracted from the BAR and injected into the RCS. The beam instrumentation needs to provide reliable measurements with the required accuracy over the dynamic range from 0.1 nC (one tenth of a typical injected bunch charge) to 32 nC – the maximum expected accumulated current. There are 10 Beam Position Monitors (BPMs) in the ring, 1 extra button-electrode assembly for the RF system, and 8 BPMs in the beam transport lines. The BPM locations are shown in Fig. 1, marked by blue rectangles. The ring BPMs will be capable of both average orbit and turn-by-turn measurement modes. Accurate measurement of the beam position with a large horizontal offset requires polynomial correction of the BPM nonlinearity.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Charge And Dynamic Current On Tubular Antennas For Various Drive Conditions

The mixed boundary value problem of a tubular conductor is solved using an approximate representation of its Fourier coefficients. A two term solution is derived, which represents the solution over an extremely broad range of aspect ratios. This representation is used to find the electrostatic solution and capacitance of a charged tube as well as the solution of a tube in a uniform field and its dipole moment. This second case is directly useful as a model for a monopole electric field probe. This approximation is a special case of a representation using a combination of Chebyshev and Legendre polynomials. Combining the charged tube and tube in a uniform field allows the solution of voltage driven tubular antennas. Comparisons are made with numerical solutions using piecewise sinusoidal representations of the current. The results are also generalized to the dynamic case and up to and beyond the first resonance. Simple corrections for finite gap and delta gap drives to magnetic frill drives are examined using infinite tube integral transform representations. Corrections between magnetic frill drives and coaxial drives are also given. Approximate drive corrections using conformal mapping and an effective radius are also discussed. Finally, this efficient current representation is applied to the magnetic problem involving simple tubular solenoids.

97 MATHEMATICS AND COMPUTING↗