Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eigenvalue”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Pulsed-Neutron Die-Away Response of H 2 O Targets to a D-T Generator Pulse

Pulsed-neutron die-away (PNDA) experiments were completed at Lawrence Livermore National Laboratory (LLNL). The goal of these experiments was to provide a benchmark to validate thermal neutron scattering laws of H 2 O. The experiment was conducted with a deuterium-tritium (D-T) neutron generator producing pulses of 14.1 MeV neutrons that impinged on a moderating target. The neutrons scattered and thermalized within the target and were counted as a function of time using Helium-3 ( 3 He) detectors that surrounded the target. These data provide a time-decay profile of the neutron population from which the time eigenvalue of the experiment is calculated. The time eigenvalue quantity α represents the integral parameter of interest. This report describes the measurements for H 2 O targets. The measurements were conducted over several days beginning on October 23rd, 2023. All measurements were completed at the low-scatter facility at LLNL. The evaluation identifier is FUND-LLNL-DT-H2O-PNDA-001.

3He Detectors↗
DOE OSTI2024

Bilocal Field Theory for Composite Scalar Bosons

We give a bilocal field theory description of a composite scalar with an extended binding potential that reduces to the Nambu–Jona-Lasinio (NJL) model in the pointlike limit. This provides a description of the internal dynamics of the bound state and features a static internal wave function, $\phi(\vec{r})$, in the center-of-mass frame that satisfies a Schrödinger–Klein–Gordon equation with eigenvalues m 2 . We analyze the “coloron” model (single perturbative massive gluon exchange) which yields a UV completion of the NJL model. This has a BCS-like enhancement of its interaction, ∝N c the number of colors, and is classically critical with g critical remarkably close to the NJL quantum critical coupling. Negative eigenvalues for m 2 lead to spontaneous symmetry breaking, and the Yukawa coupling of the bound state to constituent fermions is emergent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Pulsed-Neutron Experiments at the Inherently Safe Subcritical Assembly

The pulsed neutron technique is a powerful, dynamic method to assay the reactivity of a multiplying system. This work presents the novel application of the pulsed neutron technique to the Inherently Safe Subcritical Assembly, an experimental configuration accepted by the International Criticality Safety Benchmark Evaluation Project Handbook. The experiments were replicated with COG11.3, a continuous-energy Monte Carlo code. The pulsed neutron data were analyzed using the Sjöstrand and Gozani area-ratio methods and by extracting the prompt neutron decay constant. Subsequent static k-eigenvalue and 𝛼-eigenvalue simulations were also performed for the same configurations. Neutron detector dead-time effects from the experiments were corrected using the Backwards Extrapolation Method and shown to have a negligible impact on the estimated reactivities. The results highlight that capturing time-dependent effects like delayed neutron precursor buildup are essential to accurately reproduce experimental results. They also highlight the importance of shielding the detectors from generator source neutrons in deeply subcritical configurations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Filtered Rayleigh-Ritz is all you need

Recent work has shown that the (block) Lanczos algorithm can be used to extract approximate energy spectra and matrix elements from (matrices of) correlation functions in quantum field theory, and identified exact coincidences between Lanczos analysis methods and others. In this work, we note another coincidence: the Lanczos algorithm is equivalent to the well-known Rayleigh-Ritz method applied to Krylov subspaces. Rayleigh-Ritz provides optimal eigenvalue approximations within subspaces; we find that spurious-state filtering allows these optimality guarantees to be retained in the presence of statistical noise. We explore the relation between Lanczos and Prony's method, their block generalizations, generalized pencil of functions (GPOF), and methods based on the generalized eigenvalue problem (GEVP), and find they all fall into a larger "Prony-Ritz equivalence class", identified as all methods which solve a finite-dimensional spectrum exactly given sufficient correlation function (matrix) data. This equivalence allows simpler and more numerically stable implementations of (block) Lanczos analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

FOS: Physical sciences↗
DOE OSTI2021

On the Stochastic Stability of Deep Markov Models

Deep Markov models (DMM) are generative models which are scalable and expressive generalization of Markov models for representation, learning, and inference problems. DMMs using deep neural networks to parametrize the transition of Markov probability distributions have recently been shown to provide more expressiveness in modeling sequential data and dynamical system responses. However, the fundamental stochastic stability guarantees of such models have not been thoroughly investigated. In this paper, we present a rigorous analytical method to prove the necessary and sufficient conditions of DMM's stochastic stability. This task is achieved by spectral analysis of the efficiently computed Jacobians of probabilistic maps modeled by deep neural networks. We make theoretical connections between the eigenvalues of neural network's weights and the different activation function types used on the stability and overall dynamic behavior of DMMs with Gaussian distributions. We empirically substantiate our theoretical results on stochastic stability and eigenvalue spectra via several numerical experiments. Formal stability guarantees of DMMs can substantially improve their robustness and trustworthiness, necessary for reliable use in safety-critical real-world applications.

Drgona, Jan↗
DOE OSTI2022

Multigroup cross section generation capability in GRIFFIN

GRIFFIN is an advanced reactor multiphysics application built on the object-oriented simulation environment (MOOSE) and is jointly developed by Idaho National Laboratory and Argonne National Laboratory. The cross section application programming interface, originally developed for the PROTEUS code, has been integrated into GRIFFIN to prepare cross sections for thermal reactor applications with heterogeneous geometries. Additional improvements have been made by implementing an on-the-fly slowing down method, a double heterogeneity treatment capability, and updating the procedure to generate the fine multigroup library. The cross section preparation capability in GRIFFIN was verified for graphite-moderated TRISO fuel-based reactor benchmark problems: unit-cell problems of VHTR and EMPIRE micro reactor and HTTR assembly problems. Eigenvalues and multigroup cross sections of GRIFFIN agreed very well with those of the continuous-energy Monte Carlo code Serpent2 within 200 pcm in eigenvalue and 2% in cross sections. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

3D full core BEAVRS using OpenMOC with transport equivalence

Using an optimized implementation of the three-dimension method of characteristics for neutron transport, as well as a novel equivalence method for transport calculations de- signed to recover self-shielding errors from neglecting the angular dependence of reso- nant group absorption, a 3D full core light water reactor hybrid stochastic-deterministic eigenvalue calculation was achieved. This paper presents the optimizations developed and compares the transport solutions obtained. It realizes a runtime near 6,000 CPU hours for the statepoint, an improvement of an order of magnitude on previous works, with less than 1% error on pin fission/U-238 capture ratios and a few dozen pcms on the eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Development of a Representative Molten Chloride Fast Reactor Model to Assess the Impact of Nuclear Data

The SCALE code system was employed to conduct a preliminary investigation of nuclear data impacts for a fast spectrum molten chloride salt reactor. A computationally effective depletion model that is representative of the reactor system was successfully developed and used to conduct fuel depletion simulations. Development of this model draws from the SLICE method that was developed at Oak Ridge National Laboratory to enable generation of fuel compositions for an advanced reactor core at equilibrium operation. Eigenvalue uncertainty calculations using the ENDF/B-VII.1 nuclear data library were performed for the reactor in the fresh fuel state and an irradiated fuel state. It was determined that the primary driver of eigenvalue uncertainty was the uncertainty in the 235U (n, 𝛾) cross section. Uncertainty calculation results from this study were compared to results available for a different fast system, a sodium-cooled fast reactor, to confirm similarities and identify differences with respect to nuclear data impacts between the two fast advanced reactor systems.

Hirji, Rakim [Georgia Institute of Technology]↗
DOE OSTI2021

Relativistic levels of mesic atoms

We revisit the derivation of the covariant two-body scalar-fermion equation with a Coulomb interaction, presented in a previous paper. We show that it can be given the formal aspect of a Dirac equation, but for the fact that the eigenvalue is also contained in one of the coefficients and thus it is not linearly included. The discussion of the boundary value problem is therefore different, although some properties of the Dirac equation can be recovered in an approximation bringing back to the concept of reduced mass. We discuss a mixed analytic-numerical method of solution which allows to obtain very accurate results and we calculate the lowest levels, states and QED corrections for pionic and kaonic atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the missing magnetic flux and topological effects of a screw dislocation on a charged particle in an inhomogeneous magnetic field

Highlights: • Aharonov–Bohm-type effect associated with the missing magnetic flux. • Persistent currents associated with the missing magnetic flux. • Aharonov–Bohm-type effect associated with the screw dislocation topology. We study the interaction of an electron/hole with inhomogeneous magnetic field in the presence of a screw dislocation. We consider the internal structure of the defect, i.e., we consider the core that gives rise to a finite size to the defect. In addition, we assume that this core determines a forbidden region for the electron/hole. Then, we solve the Schrödinger equation and show that an Aharonov–Bohm-type effect arises from the influence of the topological defect and the missing magnetic flux on the eigenvalues of energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hybrid eigensolvers for nuclear configuration interaction calculations

We examine and compare several iterative methods for solving large-scale eigenvalue problems arising from nuclear structure calculations. In particular, we discuss the possibility of using block Lanczos method, a Chebyshev filtering based subspace iterations and the residual minimization method accelerated by direct inversion of iterative subspace (RMM-DIIS) and describe how these algorithms compare with the standard Lanczos algorithm and the locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm. Although the RMM-DIIS method does not exhibit rapid convergence when the initial approximations to the desired eigenvectors are not sufficiently accurate, it can be effectively combined with either the block Lanczos or the LOBPCG method to yield a hybrid eigensolver that has several desirable properties. We will describe a few practical issues that need to be addressed to make the hybrid solver efficient and robust.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗
DOE OSTI2022

Sampling electronic structure quadratic unconstrained binary optimization problems (QUBOs) with Ocean and Mukai solvers

The most advanced D-Wave Advantage quantum annealer has 5000+ qubits, however, every qubit is connected to a small number of neighbors. As such, implementation of a fully-connected graph results in an order of magnitude reduction in qubit count. To compensate for the reduced number of qubits, one has to rely on special heuristic software such as qbsolv, the purpose of which is to decompose a large quadratic unconstrained binary optimization (QUBO) problem into smaller pieces that fit onto a quantum annealer. In this work, we compare the performance of the open-source qbsolv which is a part of the D-Wave Ocean tools and a new Mukai QUBO solver from Quantum Computing Inc. (QCI). The comparison is done for solving the electronic structure problem and is implemented in a classical mode (Tabu search techniques). The Quantum Annealer Eigensolver is used to map the electronic structure eigenvalue-eigenvector equation to a QUBO problem, solvable on a D-Wave annealer. We find that the Mukai QUBO solver outperforms the Ocean qbsolv with one to two orders of magnitude more accurate energies for all calculations done in the present work, both the ground and excited state calculations. This work stimulates the further development of software to assist in the utilization of modern quantum annealers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

New Cross Section Measurements and Evaluations of Zirconium Isotopes [Slides]

This presentation gives a look to new cross section measurements and evaluations of Zirconium (Zr) isotopes. This lecture includes highlights of accomplishments made in 2022. Additionally, it explains the motivation for investigation, new measurements, and re-evaluation of the Zr isotopes. This presentation concludes with a preliminary evaluation result from RPI 90 Zr evaluation and preliminary results from Eigenvalue sensitivity study on Zr sensitive critical benchmark experiments.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2026

Control Parameter Sensitivity Study for Inverter-Based-Resource Dominated Grids: A Small Signal Stability Approach and Framework

The growing adoption of renewable energy is driving the prevalence of inverter-based resources (IBRs) within power grids. Future power grids will integrate both grid-following IBRs (GFM-IBRs) and grid-forming IBRs (GFL-IBRs) alongside synchronous generators. Therefore, it is crucial to perform stability studies that account for all components and especially control interactions related to IBRs. Extensive research has performed to study the IBR-related stability, however, the sensitivity study of IBRs' control parameters on system stability has not been adequately studied yet, especially from a systematic way. Therefore, this paper conducts a small signal stability analysis for a generic grid with multiple types of resources and develops an analytical framework for assessing the sensitivity of control parameters affecting stability margins. To achieve that, the non-autonomous reduced-order non-linear dynamic model is developed for a generic power system with multiple synchronous generator-based resources (SGBRs), GFM-IBRs, and GFL-IBRs. Based on the analytic model, a systematic framework for parametric sensitivity on systems' asymptotic stability is developed. A parameter sensitivity analysis based on eigenvalue methods is proposed. The impact of the droop controllers of GFM-IBRs, PQ-dispatch and the PLL controller of GFL-IBR on the system asymptotic stability is discussed. This sensitivity study is aiming to provide deep insights on control parameters' impact on system stability, and gives direction for parameter tuning in case of instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Bounds on spectral gaps of Hyperbolic spin surfaces

We describe a method for constraining Laplacian and Dirac spectra of two dimensional compact orientable hyperbolic spin manifolds and orbifolds. The key ingredient is an infinite family of identities satisfied by the spectra. These spectral identities follow from the consistency between 1) the spectral decomposition of functions on the spin bundle into irreducible representations of SL(2,R) and 2) associativity of pointwise multiplication of functions. Applying semidefinite programming methods to our identities produces rigorous upper bounds on the Laplacian spectral gap as well as on the Dirac spectral gap conditioned on the former. In several examples, our bounds are nearly sharp; a numerical algorithm based on the Selberg trace formula shows that the [0;3,3,5] orbifold, a particular surface with signature [1;3], and the Bolza surface nearly saturate the bounds at genus 0, 1 and 2 respectively. Under additional assumptions on the number of harmonic spinors carried by the spin-surface, we obtain more restrictive bounds on the Laplacian spectral gap. In particular, these bounds apply to hyperelliptic surfaces. We also determine the set of Laplacian spectral gaps attained by all compact orientable two-dimensional hyperbolic spin orbifolds. We show that this set is upper bounded by 12.13798; this bound is nearly saturated by the [0;3,3,5] orbifold, whose first non-zero Laplacian eigenvalue is λ^(0)_1 ≈ 12.13623.

Spectral theory↗
DOE OSTI2022

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Abstract]

The SCALE code package offers several nuclear data libraries to support Monte Carlo (MC) transport, as well as MC-based derivation of $\kappa$ eff sensitivity and uncertainty (S/U) data. The CSAS sequence using the KENO MC code can utilize continuous-energy (CE) cross sections, or pre-generated multigroup (MG) cross section libraries. The use of MG libraries introduces bias into calculations in exchange for faster transport solutions. The TSUNAMI-3D sequence also utilizes KENO MC calculations. TSUNAMI-3D has two CE calculational methods: the Iterated Fission Probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications, e.g., fissionable and polyethylene reflectors

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗