Contributions to the elliptic restricted three- body problem.
Elliptic restricted three-body problem solved for eigenvalues and orbital series solutions of Lagrangian triangular point
ORBITAL CALCULATION↗
SEARCH · Search NASA
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.
The application of finite element techniques to acoustic transmission in lined ducts with flow
The finite element method (FEM) is used to analyze the propagation of sound in two-dimensional nonuniform ducts carrying a compressible subsonic mean flow. Galerkin and residual least squares (RLS) methods with natural and forced boundary conditions are considered. The accuracy of FEM results for the eigenvalue and transmission problems is assessed by comparison with alternative numerical schemes for nonuniform ducts. The results presented and those from associated investigations indicate that modal coupling is a significant feature of the acoustic field, especially at high Mach numbers. A multimodal model therefore appears to be essential if any reliable conclusions are to be drawn in the context of turbofan inlet regions. Improvements to the eigenvalue scheme following the implementation of higher-order Hermitian elements indicate a similar modification for the transmission problem.
Astley, R. J.↗
An Analytic Benchmark for the Solution to the Isotopic Fission Spectrum Mixture Problem
In neutron transport calculations, it is common to specify material compositions in terms of constituent isotopes. Material compositions may be described by isotope number densities and associated microscopic cross sections. For general reaction cross sections, the macroscopic cross sections of a composition are simply the summation of the sum of the products of isotope number densities and microscopic cross sections. A notable exception is the mixture of the neutron fission spectrum in fissile material. To demonstrate proper and improper mixture of the neutron fission spectrum, a zero-dimensional test problem is developed. Using the test problem, it is demonstrated that the improper mixing of the fission spectrum results in an eigenvalue error of approximately 65 pcm. Though this may seem small, any error in such a simple calculation is unacceptable. In similar infinite-homogeneous test problems, eigenvalue errors of more than 300 pcm have been observed. It is also shown that the error due to fission spectrum mixing can be exaggerated for computer program verification purposes. It is concluded that the proper mixture of the neutron fission spectrum is essential for accurate neutron transport simulations. The effect of improper neutron fission spectrum mixture is demonstrated and quantified and this test problem may be used for verification purposes in the future.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Bounds on lamda for positive solutions of perturbation theory boundary value problem.
Green function-eigenvalue relationship on closed smooth boundary of open dimensional region
PARTIAL DIFFERENTIAL EQUATION↗
Stochastic gradient descent for optimization for nuclear systems
The use of gradient descent methods for optimizing k-eigenvalue nuclear systems has been shown to be useful in the past, but the use of k-eigenvalue gradients have proved computationally challenging due to their stochastic nature. ADAM is a gradient descent method that accounts for gradients with a stochastic nature. This analysis uses challenge problems constructed to verify if ADAM is a suitable tool to optimize k-eigenvalue nuclear systems. ADAM is able to successfully optimize nuclear systems using the gradients of k-eigenvalue problems despite their stochastic nature and uncertainty. Furthermore, it is clearly demonstrated that low-compute time, high-variance estimates of the gradient lead to better performance in the optimization challenge problems tested here.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Comments on numerical solution of boundary value problems of the Laplace equation and calculation of eigenvalues by the grid method
The mathematics involved in numerically solving for the plane boundary value of the Laplace equation by the grid method is developed. The approximate solution of a boundary value problem for the domain of the Laplace equation by the grid method consists of finding u at the grid corner which satisfies the equation at the internal corners (u=Du) and certain boundary value conditions at the boundary corners.
Lyusternik, L. A.↗
Observations on the Computation of Eigenvalue and Eigenvector Jacobians
Many scientific and engineering problems benefit from analytic expressions for eigenvalue and eigenvector derivatives with respect to the elements of the parent matrix. While there exists extensive literature on the calculation of these derivatives, which take the form of Jacobian matrices, there are a variety of deficiencies that have yet to be addressed — including the need for both left and right eigenvectors, limitations on the matrix structure, and issues with complex eigenvalues and eigenvectors. This work addresses these deficiencies by proposing a new analytic solution for the eigenvalue and eigenvector derivatives. The resulting analytic Jacobian matrices are numerically efficient to compute and are valid for the general complex case. It is further shown that this new general result collapses to previously known relations for the special cases of real symmetric matrices and real diagonal matrices. Finally, the new Jacobian expressions are validated using forward finite differencing and performance is compared with another technique.
Jacobian↗
Controlled precision QUBO-based algorithm to compute eigenvectors of symmetric matrices
We describe an algorithm to compute the extremal eigenvalues and corresponding eigenvectors of a symmetric matrix which is based on solving a sequence of Quadratic Binary Optimization problems. This algorithm is robust across many different classes of symmetric matrices; It can compute the eigenvector/eigenvalue pair to essentially any arbitrary precision, and with minor modifications, can also solve the generalized eigenvalue problem. Performance is analyzed on small random matrices and selected larger matrices from practical applications.
97 MATHEMATICS AND COMPUTING↗
OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport
OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Stability of expanding accretion shocks for an arbitrary equation of state
We present a theoretical stability analysis for an expanding accretion shock that does not involve a rarefaction wave behind it. The dispersion equation that determines the eigenvalues of the problem and the explicit formulae for the corresponding eigenfunction profiles are presented for an arbitrary equation of state and finite-strength shocks. For spherically and cylindrically expanding steady shock waves, we demonstrate the possibility of instability in a literal sense, a power-law growth of shock-front perturbations with time, in the range of $h_c< h<1+2 {\mathcal {M}}_2$ , where $h$ is the D'yakov-Kontorovich parameter, $h_c$ is its critical value corresponding to the onset of the instability and ${\mathcal {M}}_2$ is the downstream Mach number. Shock divergence is a stabilizing factor and, therefore, instability is found for high angular mode numbers. As the parameter $h$ increases from $h_c$ to $1+2 {\mathcal {M}}_2$ , the instability power index grows from zero to infinity. This result contrasts with the classic theory applicable to planar isolated shocks, which predicts spontaneous acoustic emission associated with constant-amplitude oscillations of the perturbed shock in the range $h_c< h<1+2 {\mathcal {M}}_2$ . Examples are given for three different equations of state: ideal gas, van der Waals gas and three-terms constitutive equation for simple metals.
Huete, César (ORCID:0000000232278520)↗
Multilevel-in-Space-and-Energy CMFD in VERA
For full-core modeling in the Virtual Environment for Reactor Analysis (VERA), the three-dimensional multigroup eigenvalue neutron transport problem is solved by MPACT. To improve the efficiency of MPACT, advancements have been made in the transport accelerator. Multilevel-in-energy and multilevel-in-space coarse mesh finite difference (CMFD) solvers were developed to improve the efficiency of the CMFD accelerator. In this paper a new multilevel-in-space-and-energy CMFD solver is developed with coarsening in both space and energy on every level. Several different strategies are investigated for coarsening groups in energy. Modified V-cycle and multiple-cycle algorithms are evaluated for solving the multilevel equations. The performance of these solvers is compared for typical full-core reactor physics problems.
42 ENGINEERING↗
Aircraft model prototypes which have specified handling-quality time histories
Several techniques for obtaining linear constant-coefficient airplane models from specified handling-quality time histories are discussed. One technique, the pseudodata method, solves the basic problem, yields specified eigenvalues, and accommodates state-variable transfer-function zero suppression. The method is fully illustrated for a fourth-order stability-axis small-motion model with three lateral handling-quality time histories specified. The FORTRAN program which obtains and verifies the model is included and fully documented.
Johnson, S. H.↗
Aircraft model prototypes which have specified handling-quality time histories
Several techniques for obtaining linear constant-coefficient airplane models from specified handling-quality time histories are discussed. The pseudodata method solves the basic problem, yields specified eigenvalues, and accommodates state-variable transfer-function zero suppression. The algebraic equations to be solved are bilinear, at worst. The disadvantages are reduced generality and no assurance that the resulting model will be airplane like in detail. The method is fully illustrated for a fourth-order stability-axis small motion model with three lateral handling quality time histories specified. The FORTRAN program which obtains and verifies the model is included and fully documented.
Johnson, S. H.↗
Effect of reactive intermediates on the propagation of turbulent premixed flames
The most simple version of the Bray and Libby model of turbulent flame propagation has been extended to include a two-step reaction mechanism in which reactants decompose into reactive intermediates which then form products. A transport equation for the joint pdf of intermediate and product concentrations was derived, the pdf parameterized, and transport equations for the parameters obtained. The one-dimensional concentration equations were then solved for the parameters of the pdf and the turbulent kinetic energy. Flame speed and angle are recovered as eigenvalues of the problem. The results show the essential features observed in experiments. The product pdf is composed of peaks at zero and unity product concentration, with a continuous structure representing partially reacted gas. The intermediates pdf has a peak at zero concentration and a broad bell-shaped distribution at nonzero concentration, representing partially reacted gas. The results also show the effect of partitioning of energy release on the flame thickness and speed.
Azzazy, M.↗
Recent developments in structural sensitivity analysis
Recent developments are reviewed in two major areas of structural sensitivity analysis: sensitivity of static and transient response; and sensitivity of vibration and buckling eigenproblems. Recent developments from the standpoint of computational cost, accuracy, and ease of implementation are presented. In the area of static response, current interest is focused on sensitivity to shape variation and sensitivity of nonlinear response. Two general approaches are used for computing sensitivities: differentiation of the continuum equations followed by discretization, and the reverse approach of discretization followed by differentiation. It is shown that the choice of methods has important accuracy and implementation implications. In the area of eigenproblem sensitivity, there is a great deal of interest and significant progress in sensitivity of problems with repeated eigenvalues. In addition to reviewing recent contributions in this area, the paper raises the issue of differentiability and continuity associated with the occurrence of repeated eigenvalues.
Haftka, Raphael T.↗
Block Lanczos Algorithm For Gyroscopic Systems
Report describes details of implementation of procedure for accurate and economical computation of natural frequencies and associated vibrational modes of elastic structure rotating along arbitrary axis. Block version of Lanczos algorithm derived for solution of eigenvalue and eigenvector problems. Fully exploits sparsity of associated matrices and employs only real numbers in all relevant computations. Capable of determining multiple roots and proves to be most efficient when compared to other similar existing techniques.
Gupta, Kajal K.↗
Quantized discrete space oscillators
A quasi-canonical sequence of finite dimensional quantizations was found which has canonical quantization as its limit. In order to demonstrate its practical utility and its numerical convergence, this formalism is applied to the eigenvalue and 'eigenfunction' problem of several harmonic and anharmonic oscillators.
Uzes, C. A.↗
Quantum statistics of Raman scattering model with Stokes mode generation
The model describing three coupled quantum oscillators with decay of Rayleigh mode into the Stokes and vibration (phonon) modes is examined. Due to the Manley-Rowe relations the problem of exact eigenvalues and eigenstates is reduced to the calculation of new orthogonal polynomials defined both by the difference and differential equations. The quantum statistical properties are examined in the case when initially: the Stokes mode is in the vacuum state; the Rayleigh mode is in the number state; and the vibration mode is in the number of or squeezed states. The collapses and revivals are obtained for different initial conditions as well as the change in time the sub-Poisson distribution by the super-Poisson distribution and vice versa.
Tanatar, Bilal↗