Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eigenvalue problem”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Convective mixing in vertically-layered porous media: The linear regime and the onset of convection

In this paper, we study the effect of permeability heterogeneity on the stability of gravitationally unstable, transient, diffusive boundary layers in porous media. Permeability is taken to vary periodically in the horizontal plane normal to the direction of gravity. In contrast to the situation for vertical permeability variation, the horizontal perturbation structures are multimodal. We therefore use a two-dimensional quasi-steady eigenvalue analysis as well as a complementary initial value problem to investigate the stability behavior in the linear regime, until the onset of convection. We find that thick permeability layers enhance instability compared with thin layers when heterogeneity is increased. On the contrary, for thin layers the instability is weakened progressively with increasing heterogeneity to the extent that the corresponding homogeneous case is more unstable. For high levels of heterogeneity, we find that a small change in the permeability field results in large variations in the onset time of convection, similar to the instability event in the linear regime. However, this trend does not persist unconditionally because of the reorientation of vorticity pairs due to the interaction of evolving perturbation structures with heterogeneity. Consequently, an earlier onset of instability does not necessarily imply an earlier onset of convection. A resonant amplification of instability is observed within the linear regime when the dominant perturbation mode is equal to half the wavenumber of permeability variation. On the other hand, a substantial damping occurs when the perturbation mode is equal to the harmonic and sub-harmonic components of the permeability wavenumber. Finally, the phenomenon of such harmonic interactions influences both the onset of instability as well as the onset of convection.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fast calculation of the tokamak vertical instability

There has been recent interest in fast calculations of the tokamak axisymmetric vertical instability for real time feedback control purposes. It is shown that the maximum eigenvalue for the basic rigid version of this stability problem can be obtained by finding the positive root to a simple scalar function. This function can be generalized to include plasma mass and has complexity linear in the number of conductive elements. The formulation is based on standard matrix decompositions of the fixed-geometry part of the eigenproblem. The calculation bottleneck is the summary of mutual inductances from the reconstructed equilibrium current density. The with-mass spectrum can be made fully real-valued by the addition of a critical amount of damping with negligible effect on the vertical growth rate. Furthermore, the calculation has been implemented in the plasma control system at the DIII-D tokamak and used in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Distorting general relativity: gravity's rainbow and f(R) theories at work

We compute the Zero Point Energy in a spherically symmetric background combining the high energy distortion of Gravity's Rainbow with the modification induced by a f(R) theory. Here f(R) is a generic analytic function of the Ricci curvature scalar R in 4D and in 3D. The explicit calculation is performed for a Schwarzschild metric. Due to the spherically symmetric property of the Schwarzschild metric we can compare the effects of the modification induced by a f(R) theory in 4D and in 3D. We find that the final effect of the combined theory is to have finite quantities that shift the Zero Point Energy. In this context we setup a Sturm-Liouville problem with the cosmological constant considered as the associated eigenvalue. The eigenvalue equation is a reformulation of the Wheeler-DeWitt equation which is analyzed by means of a variational approach based on gaussian trial functionals. With the help of a canonical decomposition, we find that the relevant contribution to one loop is given by the graviton quantum fluctuations around the given background. A final discussion on the connection of our result with the observed cosmological constant is also reported.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

(U) Radiation Transport Sensitivity Capabilities Applicable to INRAD

Intrinsic radiation (INRAD) calculations involve calculating the neutron and gamma-ray source rate densities from given compositions and densities of materials and the subsequent transport of the neutrons and gamma rays to detectors that may be stood off or internal. Related problems deal with measures of subcritical multiplication, such as the k eff eigenvalue. Realistic estimates of the uncertainties in calculated INRAD results are desired. One approach to uncertainty quantification begins with sensitivity analysis. This report describes the present capabilities that XCP-7 can apply to compute sensitivities in INRAD problems and areas where work is needed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Polynomial Preconditioned Arnoldi with Stability Control

Polynomial preconditioning can improve the convergence of the Arnoldi method for computing eigenvalues. Such preconditioning significantly reduces the cost of orthogonalization; for difficult problems, it can also reduce the number of matrix-vector products. Parallel computations can particularly benefit from the reduction of communication-intensive operations. Additoinally, the GMRES algorithm provides a simple and effective way of generating the preconditioning polynomial. For some problems high degree polynomials are especially effective, but they can lead to stability problems that must be mitigated. A two-level “double polynomial preconditioning” strategy provides an effective way to generate high-degree preconditioners.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Green’s Function of the Schrödinger Equation in the Potential Quantization Method

The problem for determining Green’s function G(r, r') for the time-independent Schrödinger equation is considered using the potential quantization method. The expression for G(r, r') is obtained in terms of the eigenfunctions and eigenvalues of this method. The equivalence of the used and traditional methods for solving the given problem has been proven.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Spectral properties of a confined nonlinear quantum oscillator in one and three dimensions

We analyze the spectral behaviour of a nonlinear quantum oscillator model under confinement. The underlying potential is given by a harmonic oscillator interaction plus a nonlinear term that can be weakened or strengthened through a parameter. Numerical eigenvalues of the model in one and three dimensions are presented. The asymptotic behaviour of the eigenvalues for confinement relaxation and for vanishing nonlinear term in the potential is investigated. Our findings are compared with existing results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

On a theory of stability for nonlinear stochastic chemical reaction networks

We present elements of a stability theory for small, stochastic, nonlinear chemical reaction networks. Steady state probability distributions are computed with zero-information (ZI) closure, a closure algorithm that solves chemical master equations of small arbitrary nonlinear reactions. Stochastic models can be linearized around the steady state with ZI-closure, and the eigenvalues of the Jacobian matrix can be readily computed. Eigenvalues govern the relaxation of fluctuation autocorrelation functions at steady state. Autocorrelation functions reveal the time scales of phenomena underlying the dynamics of nonlinear reaction networks. In accord with the fluctuation-dissipation theorem, these functions are found to be congruent to response functions to small perturbations. Significant differences are observed in the stability of nonlinear reacting systems between deterministic and stochastic modeling formalisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Eigenmode Analysis of Pulsed Neutron Transport Simulations

We discuss the time dependent behavior of some simple pulsed neutron simulations of subcritical problems in slab geometry. Our intent is to investigate the eigenvalue structure of the discretized neutron transport equation and to show under some reasonable assumptions that a dominant time eigenvalue exists that has a nonnegative eigenvector that determines the long time dependent behavior of the solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Resonance Extraction from the Finite Volume

The spectrum of excited hadrons becomes accessible in simulations of Quantum Chromodynamics on the lattice. Extensions of Lüscher's method allow to address multi-channel scattering problems using moving frames or modified boundary conditions to obtain more eigenvalues in finite volume. As these are at different energies, interpolations are needed to relate different eigenvalues and to help determine the amplitude. Expanding the T- or the K-matrix locally provides a controlled scheme by removing the known non-analyticities of thresholds. This can be stabilized by using Chiral Perturbation Theory. Different examples to determine resonance pole parameters and to disentangle resonances from thresholds are dis- cussed, like the scalar meson f0(980) and the excited baryons N(1535)1/2^- and Lambda(1405)1/2^-.

Doring, Michael↗
DOE OSTI2018

Localized Modes Due to Defects in High Contrast Periodic Media Via Two-Scale Homogenization

The spectral problem for an infinite periodic medium perturbed by a compact defect is considered. For a high contrast small ε-size periodicity and a finite size defect we consider the critical ε{sup 2}-scaling for the contrast. We employ two-scale homogenization for deriving asymptotically explicit limit equations for the localized modes and associated eigenvalues. Those are expressed in terms of the eigenvalues and eigenfunctions of a perturbed version of a two-scale limit operator introduced by V. V. Zhikov with an emergent explicit nonlinear dependence on the spectral parameter for the spectral problem at the macroscale. Using the method of asymptotic expansions supplemented by a high contrast boundary layer analysis, we establish the existence of the actual eigenvalues near the eigenvalues of the limit operator, with “ε square root” error bounds. An example for circular or spherical defects in a periodic medium with isotropic homogenized properties is given.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A Galerkin method with two-dimensional Haar basis functions for the computation of the Karhunen–Loève expansion

We study the numerical approximation of a homogeneous Fredholm integral equation of second kind associated with the Karhunen–Loève expansion of Gaussian random fields. We develop, validate, and discuss an algorithm based on the Galerkin method with two-dimensional Haar wavelets as basis functions. The shape functions are constructed from the orthogonal decomposition of tensor-product spaces of one-dimensional Haar functions, and a recursive algorithm is employed to compute the matrix of the discrete eigenvalue system without the explicit calculation of integrals, allowing the implementation of a fast and efficient algorithm that provides considerable reduction in CPU time, when compared with classical Galerkin methods. Numerical experiments confirm the convergence rate of the method and assess the approximation error and the sparsity of the eigenvalue system when the wavelet expansion is truncated. We illustrate the numerical solution of a diffusion problem with random input data with the present method. In this problem, accuracy was retained after dropping the coefficients below a threshold value that was numerically determined. A similar method with scaling functions rather than wavelet functions does not need a discrete wavelet transform and leads to eigenvalue systems with better conditioning but lower sparsity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implementation of a Pressure Correction Term in the CTF Subchannel Code

CTF is a modernized and improved version of the legacy subchannel code COBRA-TF, being jointly developed and maintained by Oak Ridge National Laboratory (ORNL) and North Carolina State University. The code was adopted by ORNL for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) in 2012 to help address CASL challenge problems. Since that time, activities related to CTF have included implementing software quality assurance measures, implementing new features and models, performing validation and verification testing, establishing and supporting a CTF user group, and developing the code for use in coupled applications. CTF is applicable to single- and two-phase flows in light water reactor geometries at normal and accident operating conditions. The two-phase flow model consists of the six-equation two-phase flow model augmented by equations for the droplet field, and it also includes appropriate source terms. This work focuses on the numerical stability of the six-equation two-phase flow model and is part of a recent effort to improve CTF capabilities for simulating flows in boiling water reactors (BWRs). The six-equation two-phase flow model is not strongly hyperbolic, meaning that some of its eigenvalues are either not real or not distinct. In a Cauchy sense, and for an initial-value problem, strong hyperbolicity was found to be equivalent to well-posedness, which ensures the uniqueness of the solution in the space-time used. This concept is also of importance when designing high-order numerical methods and also when implementing boundary conditions that rely on real eigenvalues such as Riemann solvers. A tremendous amount of research is available in the literature on how to hyperbolize the six-equation two-phase flow model by adding ad hoc terms to the system of equations. Such terms are often referred to as virtual mass or interfacial pressure correction terms. From a mathematical perspective, this approach adds terms to the phasic momentum equation containing first-order partial derivatives of the unknowns to modify the Jacobian matrix, thus the imaginary part of the eigenvalues. For this work, an interfacial pressure correction (IPC) term on the same model as what was used for the CATHARE system code, is implemented in CTF. This allows for recovery of a two-pressure formulation. In the following, the six-equation two-phase flow model is recalled, along with the IPC term and its influence on the eigenvalues. Then, numerical results and convergence studies are presented for a water-faucet problem at a pressure of 10 bar. Finally, conclusions and future work are addressed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Topology and many-body localization

We discuss the problem of localization in two dimensional electron systems in the quantum Hall (single Landau level) regime. After briefly summarizing the well-studied problem of Anderson localization in the non-interacting case, we concentrate on the problem of disorder induced many-body localization (MBL) in the presence of electron–electron interactions using numerical exact diagonalization and eigenvalue spacing statistics as a function of system size. We provide evidence showing that MBL is not attainable in a single Landau level with short range (white noise) disorder in the thermodynamic limit. We then study the interplay of topology and localization, by contrasting the behavior of topological and nontopological subbands arising from a single Landau level in two models — (i) a pair of extremely flat Hofstadter bands with an optimally chosen periodic potential, and (ii) a Landau level with a split-off nontopological impurity band. Both models provide convincing evidence for the strong effect of topology on the feasibility of many-body localization as well as slow dynamics starting from a nonequilibrium state with charge imbalance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Analytical solutions for quantum walks on 1D chain with different shift operators

In this paper, we study the discrete-time quantum walks on 1D Chain with the moving and swapping shift operators. We derive analytical solutions for the eigenvalues and eigenstates of the evolution operator U{sup -hat} using the Chebyshev polynomial technique, and calculate the long-time averaged probabilities for the two different shift operators respectively. It is found that the probability distributions for the moving and swapping shift operators display completely different characteristics. For the moving shift operator, the probability distribution exhibits high symmetry where the probabilities at mirror positions are equal. The probabilities are inversely proportional to the system size N and approach to zero as N→∞. On the contrary, for the swapping shift operator, the probability distribution is not symmetric, the probability distribution approaches to a power-law stationary distribution as N→∞ under certain coin parameter condition. We show that such power-law stationary distribution is determined by the eigenstates of the eigenvalues ±1 and calculate the intrinsic probability for different starting positions. Our findings suggest that the eigenstates corresponding to eigenvalues ±1 play an important role for the swapping shift operator. - Highlights: • QWs on 1D chain with the moving and swapping operators are studied for the first time. • We derive analytical results for the probability distribution for the two operators. •We compare the dynamics of QWs with two different shift operators. • We find the particular eigenvalues ±1 play an important role for the dynamics. • We use the Chebyshev technique to treat the problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗