Search NASA⌕ Search

SEARCH · Search NASA

Results for “Recurrence relations”

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

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗
NASA NTRS1971

Power series representation of partial derivatives required in orbit determination.

A program to integrate the equations of motion by series in powers of the time step can be easily modified to furnish the elements of the matrizant in power series of the time step. In particular, if the series representing the motion are obtained recursively, differentiation of the recurrence relations will provide immediately a recursive scheme for computing the coefficient in the power series for the elements of the matrizant.

Schanzle, A. F.↗
NASA NTRS1974

Operations on power series and higher derivatives of functions

Equations are derived for the general terms of the series obtained when certain operations are carried out on one or several given series. Explicit and recurrence relations are given for the literal groups which are common to the given series. A general method is shown to obtain the multipliers which depend on the operation being performed, and the method is illustrated for several examples.

Chang, F.-C.↗
NASA NTRS1976

Pines nonsingular gravitational potential derivation, description and implementation

An engineering interpretation of and some minor corrections to previously completed work on a uniform representation of the gravitational potential and its derivatives (with emphasis on avoiding the usual singularity at the pole) were presented. The physical meaning of the variables was explained, the derivation of results was separated into smaller parts for easier reading, some additional recurrence relations for the derived Legendre polynomials were included and compared, and a computer program implementing this formulation was presented. Numerical experiments have shown that the use of this representation, besides removing the singularity, substantially increases the speed of the computation.

Spencer, J. L.↗
NASA NTRS1976

Prolongation structures of nonlinear evolution equations. II

The prolongation structure of a closed ideal of exterior differential forms is further discussed, and its use illustrated by application to an ideal (in six dimensions) representing the cubically nonlinear Schroedinger equation. The prolongation structure in this case is explicitly given, and recurrence relations derived which support the conjecture that the structure is open - i.e., does not terminate as a set of structure relations of a finite-dimensional Lie group. We introduce the use of multiple pseudopotentials to generate multiple Baecklund transformation, and derive the double Baecklund transformation. This symmetric transformation concisely expresses the (usually conjectured) theorem of permutability, which must consequently apply to all solutions irrespective of asymptotic constraints.

Estabrook, F. B.↗
NASA NTRS1977

Asymptotic behavior of vector recurrences with applications

An investigation is conducted of the behavior of a vector recurrence relation under very weak assumptions. The obtained results are used in an analysis of the iterations for nonlinear equations and a study concerning the composition of such iterations. The results can be employed in studies of one-point iterations involving the use of memory elements. Iterations for solving nonlinear equations on parallel computers can be investigated.

Feldstein, A.↗
NASA NTRS1980

The method of lines in three dimensional fracture mechanics

A review of recent developments in the calculation of design parameters for fracture mechanics by the method of lines (MOL) is presented. Three dimensional elastic and elasto-plastic formulations are examined and results from previous and current research activities are reported. The application of MOL to the appropriate partial differential equations of equilibrium leads to coupled sets of simultaneous ordinary differential equations. Solutions of these equations are obtained by the Peano-Baker and by the recurrance relations methods. The advantages and limitations of both solution methods from the computational standpoint are summarized.

Gyekenyesi, J.↗
NASA NTRS1980

The method of lines in three-dimensional fracture mechanics

A review is presented of recent developments in the calculation of design parameters for fracture mechanics by the method of lines (MOL). Three-dimensional elastic and elasto-plastic formulations are examined and results are reported from previous and current research activities. The application of MOL to the appropriate partial differential equations of equilibrium leads to coupled sets of simultaneous ordinary differential equations. Solutions of these equations are obtained by the Peano-Baker and by the recurrence relations methods. The advantages and limitations of both solution methods from the computational standpoint are summarized.

Gyekenyesi, J.↗
NASA NTRS1988

An O(N squared) method for computing the eigensystem of N by N symmetric tridiagonal matrices by the divide and conquer approach

An efficient method is proposed to solve the eigenproblem of N by N Symmetric Tridiagonal (ST) matrices. Unlike the standard eigensolvers which necessitate O(N cubed) operations to compute the eigenvectors of such ST matrices, the proposed method computes both the eigenvalues and eigenvectors with only O(N squared) operations. The method is based on serial implementation of the recently introduced Divide and Conquer (DC) algorithm. It exploits the fact that by O(N squared) of DC operations, one can compute the eigenvalues of N by N ST matrix and a finite number of pairs of successive rows of its eigenvector matrix. The rest of the eigenvectors--all of them or one at a time--are computed by linear three-term recurrence relations. Numerical examples are presented which demonstrate the superiority of the proposed method by saving an order of magnitude in execution time at the expense of sacrificing a few orders of accuracy.

Gill, Doron↗
NASA NTRS1988

Block-Krylov component synthesis method for structural model reduction

A new analytical method is presented for generating component shape vectors, or Ritz vectors, for use in component synthesis. Based on the concept of a block-Krylov subspace, easily derived recurrence relations generate blocks of Ritz vectors for each component. The subspace spanned by the Ritz vectors is called a block-Krylov subspace. The synthesis uses the new Ritz vectors rather than component normal modes to reduce the order of large, finite-element component models. An advantage of the Ritz vectors is that they involve significantly less computation than component normal modes. Both 'free-interface' and 'fixed-interface' component models are derived. They yield block-Krylov formulations paralleling the concepts of free-interface and fixed-interface component modal synthesis. Additionally, block-Krylov reduced-order component models are shown to have special disturbability/observability properties. Consequently, the method is attractive in active structural control applications, such as large space structures. The new fixed-interface methodology is demonstrated by a numerical example. The accuracy is found to be comparable to that of fixed-interface component modal synthesis.

Craig, Roy R., Jr.↗
NASA NTRS2015

Verification of a Viscous Computational Aeroacoustics Code Using External Verification Analysis

The External Verification Analysis approach to code verification is extended to solve the three-dimensional Navier-Stokes equations with constant properties, and is used to verify a high-order computational aeroacoustics (CAA) code. After a brief review of the relevant literature, the details of the EVA approach are presented and compared to the similar Method of Manufactured Solutions (MMS). Pseudocode representations of EVA's algorithms are included, along with the recurrence relations needed to construct the EVA solution. The code verification results show that EVA was able to convincingly verify a high-order, viscous CAA code without the addition of MMS-style source terms, or any other modifications to the code.

Computational↗
NASA NTRS2015

Verification of a Viscous Computational Aeroacoustics Code using External Verification Analysis

The External Verification Analysis approach to code verification is extended to solve the three-dimensional Navier-Stokes equations with constant properties, and is used to verify a high-order computational aeroacoustics (CAA) code. After a brief review of the relevant literature, the details of the EVA approach are presented and compared to the similar Method of Manufactured Solutions (MMS). Pseudocode representations of EVA's algorithms are included, along with the recurrence relations needed to construct the EVA solution. The code verification results show that EVA was able to convincingly verify a high-order, viscous CAA code without the addition of MMS-style source terms, or any other modifications to the code.

Computational↗
NASA NTRS2018

Analytic Expressions for Derivatives from Series Solutions to the Three Body Problem

This paper presents a notation system to facilitate to solution of differential equations via Taylor series expansions and applies it to solve the circular restricted three body problem. Unlike previous Taylor series methods in the astrodynamics literature, computer algebra solvers are not used. Instead the notation system allows one to solve a system of differential equations analytically “by hand” without resorting to computer algebra software. This method produces recurrence relations explicitly in terms of a sequence of derivatives of the state with respect to time for the coefficients of Taylor Series solutions that can be evaluated numerically or manipulated further to investigate properties of the solution. For example, additional derivatives with respect to other parameters may also be found, including those that describe the dependence of the solution on initial conditions.

Strange, Nathan↗
NASA NTRS1972

Recurrent active regions related to metric radio continuum emissions and the interplanetary magnetic sector structure

Active heliographic longitudes at the sun are investigated by using the observational data for long-lived metric continuum noise sources. It is shown that, for the period from 1963 to 1969, the number of such longitudes was four in general and these longitudes were very stable for this radio activity since 1963. A discussion is given on the relationship between those longitudes and the sector structure of the interplanetary magnetic field.

Sakurai, K.↗
NASA NTRS1986

Temperature-plant pigment-optical relations in a recurrent offshore mesoscale eddy near Point Conception, California

The temperature-plant pigment-optical structure of a mesoscale anticyclonic eddy consistently found in shipboard surveys and satellite-sensed data several hundred kilometers southwest of Point Conception, CA, is described on three different time scales (100-day mesoscale, annual, and several-year). The satellite coastal zone color scanner (CZCS) ocean color imagery detected the near-surface chlorophyll structure of the eddy, but in situ optical and plant pigment data suggest that such imagery does not provide a good estimate of the integrated chlorophyll field of the eddy. The temperature and plant pigment boundaries of the eddy, as determined from two-dimensional gradients of advanced very high resolution radiometer (AVHRR) and CZCS imagery, do not coincide spatially. This and in situ temperature, plant pigment, and optical structure provide additional evidence that some eddy systems in the California Current are not isolated vortex systems but rather continuously entrain waters of nonlocal origin laterally into their upper layers. Within the California Current a ratio of AVHRR/CZCS data is useful for separating inshore from oceanic water masses and following their surface entrainment by offshore vortices. The historical 28-year California Cooperative Oceanic Fisheries Investigations data for the Point Conception region of the California Current and remotely sensed data over this region show that the Point Conception eddy is a recurrent feature in the offshore California Current. Moreover, the available data provide evidence that a large number of warm-core mesoscale eddies occur simultaneously in a transition zone between coastal and oceanic regimes, that these features recur at preferred locations within the transition zone, and that this family of eddies should impose a significant offshore boundary condition on the flow of the California Current.

Simpson, James J.↗
DOE OSTI2021

Einstein–Rosen waves and the Geroch group

Under the action of the Geroch group, the Minkowski metric can be transformed into any vacuum metric with two commuting Killing vectors. In principle, this reduces the problem of deriving vacuum metrics with two commuting Killing vectors to pure algebra. In this short note, we use these facts to give a purely algebraic derivation of the Einstein-Rosen metric, which describes a cylindrical gravitational wave. Furthermore, our derivation has a straightforward extension to gravitational pulse waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗