Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rational approximations”

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

Adaptive Power Flow Approximations With Second-Order Sensitivity Insights

The power flow equations are fundamental to power system planning, analysis, and control. However, the inherent non-linearity and non-convexity of these equations present formidable obstacles in problem-solving processes. To mitigate these challenges, recent research has proposed adaptive power flow linearizations that aim to achieve accuracy over wide operating ranges. The accuracy of these approximations inherently depends on the curvature of the power flow equations within these ranges, which necessitates considering second-order sensitivities. In this paper, we leverage second-order sensitivities to both analyze and improve power flow approximations. We evaluate the curvature across broad operational ranges and subsequently utilize this information to inform the computation of various sample-based power flow approximation techniques. Additionally, we leverage second-order sensitivities to guide the development of rational approximations that yield linear constraints in optimization problems. In conclusion, this approach is extended to enhance accuracy beyond the limitations of linear functions across varied operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-like Calculations

Here, this Letter introduces a unified emulation framework for studying continuum physics in finite quantum systems. Using a reduced basis method, we construct powerful emulators for the inhomogeneous Schrödinger equation that operate in a combined parameter space of complex energy (𝐸) and other inputs (𝜽). Within the space, the emulators simultaneously perform analytical continuation in 𝐸—extracting continuum physics from numerically simpler bound-state-like calculations—and interpolate this entire process across 𝜽. This yields a small, non-Hermitian system whose properties (e.g., resonances and scattering observables) can be rapidly predicted for any 𝜽. Crucially, the complex-𝐸 emulation provides a pathway to compute continuum observables for complex systems where advanced bound-state methods exist but direct continuum calculations are yet to be developed, while the 𝜽 emulation enables rapid parameter-space exploration and can be adapted to accelerate other existing continuum calculations. Demonstrations with two- and three-body systems highlight the method’s effectiveness and suggest its connection to (near-)optimal rational approximation. This Letter presents the key results, with further details reserved for a companion paper.

ab initio calculations↗
DOE OSTI2026

pyDecay: A CRAM-Based Isotope Decay Solver

This module (pyDecay) implements a Chebyshev Rational Approximation Method (CRAM) for solving isotope decay equations, based on the work of M. Pusa. It provides a numerically stable and efficient method for evaluating the matrix exponential involved in nuclear decay calculations. This implementation of CRAM relies on the incomplete partial factorization (IPF) algorithm published by Pusa [3], with corrections noted by Romano et al.

Skutnik, SteveEugene [Oak Ridge National Laborator↗
NASA NTRS1972

System evaluation from impulse response fragments

By use of rational approximation, a procedure was found which allows calculation of the z-transfer function from a restricted number of time sequence values. The order of the system need not be known. Transformation to the s-domain is easily provided.

Baumeister, J. K.↗
NASA NTRS1975

On the analytic structure of the Taylor-Maccoll conical-flow solution

The conical flow solution for axisymmetric supersonic flow past cones has been found to be virtually independent of the ratio of specific heats when normalized in a certain way. A simple rational approximation to this flow is derived. The important singularities and the limiting behavior of the solution are also discussed.

Schwartz, L. W.↗
NASA NTRS1981

Design for active and passive flutter suppression and gust alleviation

Analytical design techniques for active and passive control of aeroelastic systems are based on a rational approximation of the unsteady aerodynamic loads in the entire Laplace domain, which yields matrix equations of motion with constant coefficients. Some existing schemes are reviewed, the matrix Pade approximant is modified, and a technique which yields a minimal number of augmented states for a desired accuracy is presented. The state-space aeroelastic model is used to design an active control system for simultaneous flutter suppression and gust alleviation. The design target is for a continuous controller which transfers some measurements taken on the vehicle to a control command applied to a control surface. Structural modifications are formulated in a way which enables the treatment of passive flutter suppression system with the same procedures by which active control systems are designed.

Karpel, M.↗
NASA NTRS1982

Real time prediction of marine vessel motions using Kalman filtering techniques

The present investigation is concerned with the prediction of the future behavior of a vessel within some confidence bounds at a specific instant of time, taking into account an interval of a few seconds. The ability to predict accurately the motions of a vessel can reduce significantly the probability of failure of operations in rough seas. The investigation was started as part of an effort to ensure safe landing of aircraft on relatively small vessels. However, the basic principles involved in the study are the same for any offshore operation, such as carbo transfer in the open sea, structure installation, and floating crane operation. The Kalman filter is a powerful tool for achieving the goals of the prediction procedure. Attention is given to a linear optimal predictor, the equations of motion of the vessel, the wave spectrum, rational approximation, the use of Kalman filter and predictor in an application for a ship, and the motions of a semisubmersible.

Triantafyllou, M. S.↗
NASA NTRS1985

Adding computationally efficient realism to Monte Carlo turbulence simulation

Frequently in aerospace vehicle flight simulation, random turbulence is generated using the assumption that the craft is small compared to the length scales of turbulence. The turbulence is presumed to vary only along the flight path of the vehicle but not across the vehicle span. The addition of the realism of three-dimensionality is a worthy goal, but any such attempt will not gain acceptance in the simulator community unless it is computationally efficient. A concept for adding three-dimensional realism with a minimum of computational complexity is presented. The concept involves the use of close rational approximations to irrational spectra and cross-spectra so that systems of stable, explicit difference equations can be used to generate the turbulence.

Campbell, C. W.↗
NASA NTRS1986

A simple method for computing the relativistic Compton scattering kernel for radiative transfer

Correct computation of the Compton scattering kernel (CSK), defined to be the Klein-Nishina differential cross section averaged over a relativistic Maxwellian electron distribution, is reported. The CSK is analytically reduced to a single integral, which can then be rapidly evaluated using a power series expansion, asymptotic series, and rational approximation for sigma(s). The CSK calculation has application to production codes that aim at understanding certain astrophysical, laser fusion, and nuclear weapons effects phenomena.

Prasad, M. K.↗
NASA NTRS1987

On implicit Runge-Kutta methods for parallel computations

Implicit Runge-Kutta methods which are well-suited for parallel computations are characterized. It is claimed that such methods are first of all, those for which the associated rational approximation to the exponential has distinct poles, and these are called multiply explicit (MIRK) methods. Also, because of the so-called order reduction phenomenon, there is reason to require that these poles be real. Then, it is proved that a necessary condition for a q-stage, real MIRK to be A sub 0-stable with maximal order q + 1 is that q = 1, 2, 3, or 5. Nevertheless, it is shown that for every positive integer q, there exists a q-stage, real MIRK which is I-stable with order q. Finally, some useful examples of algebraically stable MIRKs are given.

Keeling, Stephen L.↗
NASA NTRS1987

An algorithm for a generalization of the Richardson extrapolation process

The paper presents a recursive method, designated the W exp (m)-algorithm, for implementing a generalization of the Richardson extrapolation process. Compared to the direct solution of the linear sytems of equations defining the extrapolation procedure, this method requires a small number of arithmetic operations and very little storage. The technique is also applied to solve recursively the coefficient problem associated with the rational approximations obtained by applying a d-transformation to power series. In the course of development a new recursive algorithm for implementing a very general extrapolation procedure is introduced, for solving the same problem. A FORTRAN program for the W exp (m)-algorithm is also appended.

Ford, William F.↗
NASA NTRS1991

Response of a confined gas to volumetric heating in the absence of gravity. I - Slow transients

A one-dimensional model for bulk motion induced by a transient volumetric heat source in a confined gas at zero gravity is considered. Rational approximation methods are used to derive a quantitative theory for the gas response to a spatiallly distributed, time-dependent internal power deposition. The resulting low Mach number compressible flow equations are solved by using perturbation methods. Solutions are given for a conduction-free core and thin conductive boundary layers adjacent to the end walls. It is found that may any spatially nonuniform power deposition will cause fluid motion. Net mass transport in the closed container will occur for certain spatially distributed heating. The model mimics the thermal effects of an exothermic gas phase reaction in vapor transport experiments conducted in space. The solutions demonstrate that thermally induced mass transport can be as large as diffusive mass transport in a typical experiment.

Herczynski, A.↗
NASA NTRS1990

Efficient parallel solution of parabolic equations - Implicit methods on the Cedar multicluster

A class of implicit methods for the parallel solution of linear parabolic differential equations based on Pade and Chebyshev rational approximations to the matrix exponential are presented. It is pointed out that this approach incorporates both natural hierarchical parallelism, improved intrinsic efficiency, and fewer timesteps. These advantages lead to an extremely fast family of methods for the solution of certain time-dependent problems. These techniques are illustrated with numerical experiments on the University of Illinois Cedar multicluster architecture. The experiments indicate that implicit methods of very high degree offer great promise for the solution of certain parabolic problems when in computational environment with parallel resources. Hierarchically organized parallel computers, such as the Cedar multicluster, are found to be especially attractive for these schemes.

Gallopoulos, E.↗
NASA NTRS1993

A highly accurate Voigt function algorithm

A complex Voigt lineshape algorithm is presented whose maximum relative error over the complex plane is less than 1 x 10(exp -8). The algorithm consists of series, rational approximations and Gauss-Hermite integrations which makes it suitable as a general purpose software module for a wide variety of uses, including a Voigt function standard.

Shippony, Z.↗
NASA NTRS2001

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗
NASA NTRS1993

A Highly Accurate Voight Function Algorithm

A complex Voight lineshape algorithm is presented whose maximum relative error over the complex plane is less than 1 x 10 super minus 8. The algorithm consists of series, rational approximations and Gauss-Hermite integrations which makes it suitable as a general purpose software module for a wide variety of uses, including a Voight function standard.

Voight lineshape Doppler velocity distribution col↗
NASA NTRS2008

Role of Hf on Phase Formation in Ti45Zr(38-x)Hf(x)Ni17 Liquids and Solids

Hafnium and zirconium are very similar, with almost identical sizes and chemical bonding characteristics. However, they behave differently when alloyed with Ti and Ni. A sharp phase formation boundary near 18-21 at.% Hf is observed in rapidly-quenched and as-cast Ti45Zr38-xHfxNi17 alloys. Rapidly-quenched samples that contain less than 18 at.% Hf form the icosahedral quasicrystal phase, whiles samples containing more than 21 at.% form the 3/2 rational approximant phase. In cast alloys, a C14 structure is observed for alloys with Hf lower than the boundary concentration, while a large-cell (11.93 ) FCC Ti2Ni-type structure is found in alloys with Hf concentrations above the boundary. To better understand the role of Hf on phase formation, the structural evolution with supercooling and the solidification behavior of liquid Ti45Zr38-xHfxNi17 alloys (x=0, 12, 18, 21, 38) were studied using the Beamline Electrostatic Levitation (BESL) technique using 125keV x-rays on the 6ID-D beamline at the Advanced Photon Source, Argonne National Laboratory. For all liquids primary crystallization was to a BCC solid solution phase; interestly, an increase in Hf concentration leads to a decrease in the BCC lattice parameter in spite of the chemical similarity between Zr and Hf. A Reitveld analysis confirmed that as in the cast alloys, the secondary phase that formed was the C14 below the phase formation boundary and a Ti2Ni-type structure at higher Hf concentrations. Both the liquidus temperature and the reduced undercooling change sharply on traversing the phase formation boundary concentration, suggesting a change in the liquid structure. Structural information from a Honeycutt-Anderson index analysis of reverse Monte Carlo fits to the S(q) liquid data will be presented to address this issue.

Wessels, V.↗