Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pade 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 37 records · Page 2

On High-Order Radiation Boundary Conditions

In this paper we develop the theory of high-order radiation boundary conditions for wave propagation problems. In particular, we study the convergence of sequences of time-local approximate conditions to the exact boundary condition, and subsequently estimate the error in the solutions obtained using these approximations. We show that for finite times the Pade approximants proposed by Engquist and Majda lead to exponential convergence if the solution is smooth, but that good long-time error estimates cannot hold for spatially local conditions. Applications in fluid dynamics are also discussed.

Hagstrom, Thomas↗

Rational approximations from power series of vector-valued meromorphic functions

Let F(z) be a vector-valued function, F: C yields C(sup N), which is analytic at z = 0 and meromorphic in a neighborhood of z = 0, and let its Maclaurin series be given. In this work we developed vector-valued rational approximation procedures for F(z) by applying vector extrapolation methods to the sequence of partial sums of its Maclaurin series. We analyzed some of the algebraic and analytic properties of the rational approximations thus obtained, and showed that they were akin to Pade approximations. In particular, we proved a Koenig type theorem concerning their poles and a de Montessus type theorem concerning their uniform convergence. We showed how optical approximations to multiple poles and to Laurent expansions about these poles can be constructed. Extensions of the procedures above and the accompanying theoretical results to functions defined in arbitrary linear spaces was also considered. One of the most interesting and immediate applications of the results of this work is to the matrix eigenvalue problem. In a forthcoming paper we exploited the developments of the present work to devise bona fide generalizations of the classical power method that are especially suitable for very large and sparse matrices. These generalizations can be used to approximate simultaneously several of the largest distinct eigenvalues and corresponding eigenvectors and invariant subspaces of arbitrary matrices which may or may not be diagonalizable, and are very closely related with known Krylov subspace methods.

Sidi, Avram↗

Applications of Laplace transform methods to airfoil motion and stability calculations

This paper reviews the development of generalized unsteady aerodynamic theory and presents a derivation of the generalized Possio integral equation. Numerical calculations resolve questions concerning subsonic indicial lift functions and demonstrate the generation of Kutta waves at high values of reduced frequency, subsonic Mach number, or both. The use of rational function approximations of unsteady aerodynamic loads in aeroelastic stability calculations is reviewed, and a reformulation of the matrix Pade approximation technique is given. Numerical examples of flutter boundary calculations for a wing which is to be flight tested are given. Finally, a simplified aerodynamic model of transonic flow is used to study the stability of an airfoil exposed to supersonic and subsonic flow regions.

Edwards, J. W.↗

On the convergence of local approximations to pseudodifferential operators with applications

We consider the approximation of a class pseudodifferential operators by sequences of operators which can be expressed as compositions of differential operators and their inverses. We show that the error in such approximations can be bounded in terms of L(1) error in approximating a convolution kernel, and use this fact to develop convergence results. Our main result is a finite time convergence analysis of the Engquist-Majda Pade approximants to the square root of the d'Alembertian. We also show that no spatially local approximation to this operator can be convergent uniformly in time. We propose some temporally local but spatially nonlocal operators with better long time behavior. These are based on Laguerre and exponential series.

Hagstrom, Thomas↗

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.↗

On Time Delay Margin Estimation for Adaptive Control and Optimal Control Modification

This paper presents methods for estimating time delay margin for adaptive control of input delay systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent an adaptive law by a locally bounded linear approximation within a small time window. The time delay margin of this input delay system represents a local stability measure and is computed analytically by three methods: Pade approximation, Lyapunov-Krasovskii method, and the matrix measure method. These methods are applied to the standard model-reference adaptive control, s-modification adaptive law, and optimal control modification adaptive law. The windowing analysis results in non-unique estimates of the time delay margin since it is dependent on the length of a time window and parameters which vary from one time window to the next. The optimal control modification adaptive law overcomes this limitation in that, as the adaptive gain tends to infinity and if the matched uncertainty is linear, then the closed-loop input delay system tends to a LTI system. A lower bound of the time delay margin of this system can then be estimated uniquely without the need for the windowing analysis. Simulation results demonstrates the feasibility of the bounded linear stability method for time delay margin estimation.

Nguyen, Nhan T.↗

Well-posedness of one-way wave equations and absorbing boundary conditions

A one-way wave equation is a partial differential which, in some approximate sense, behaves like the wave equation in one direction but permits no propagation in the opposite one. The construction of such equations can be reduced to the approximation of the square root of (1-s sup 2) on -1, 1 by a rational function r(s) = p sub m (s)/q sub n(s). Those rational functions r for which the corresponding one-way wave equation is well-posed are characterized both as a partial differential equation and as an absorbing boundary condition for the wave equation. We find that if r(s) interpolates the square root of (1-s sup 2) at sufficiently many points in (-1,1), then well-posedness is assured. It follows that absorbing boundary conditions based on Pade approximation are well-posed if and only if (m, n) lies in one of two distinct diagonals in the Pade table, the two proposed by Engquist and Majda. Analogous results also hold for one-way wave equations derived from Chebyshev or least-squares approximation.

Trefethen, L. N.↗

Diagrammatic perturbation theory - The ground state of the carbon monosulfide molecule

Diagrammatic many-body perturbation theory is employed in a study of the ground state of the carbon monosulfide molecule for bond lengths close to the equilibrium value. The calculations are complete through third order in the energy within the algebraic approximation. Two different zero-order Hamiltonians are considered, and all two-, three-, and four-body terms are determined for the corresponding perturbation expansions. Many-body effects are found to be very important. Pade approximants to the energy expansion are constructed, and upper bounds evaluated. Almost 53 percent of the estimated correlation energy is recovered. The variation of components of the correlation energy with nuclear separation is investigated. Spectroscopic constants are also calculated.

Wilson, S.↗

The Investigation of Optimal Discrete Approximations for Real Time Flight Simulations

The results are presented of an investigation of discrete approximations for real time flight simulation. Major topics discussed include: (1) consideration of the particular problem of approximation of continuous autopilots by digital autopilots; (2) use of Bode plots and synthesis of transfer functions by asymptotic fits in a warped frequency domain; (3) an investigation of the various substitution formulas, including the effects of nonlinearities; (4) use of pade approximation to the solution of the matrix exponential arising from the discrete state equations; and (5) an analytical integration of the state equation using interpolated input.

Parrish, E. A.↗

Tortuosity measurement and the effects of finite pulse widths on xenon gas diffusion NMR studies of porous media

We have extended the utility of NMR as a technique to probe porous media structure over length scales of approximately 100-2000 microm by using the spin 1/2 noble gas 129Xe imbibed into the system's pore space. Such length scales are much greater than can be probed with NMR diffusion studies of water-saturated porous media. We utilized Pulsed Gradient Spin Echo NMR measurements of the time-dependent diffusion coefficient, D(t), of the xenon gas filling the pore space to study further the measurements of both the pore surface-area-to-volume ratio, S/V(p), and the tortuosity (pore connectivity) of the medium. In uniform-size glass bead packs, we observed D(t) decreasing with increasing t, reaching an observed asymptote of approximately 0.62-0.65D(0), that could be measured over diffusion distances extending over multiple bead diameters. Measurements of D(t)/D(0) at differing gas pressures showed this tortuosity limit was not affected by changing the characteristic diffusion length of the spins during the diffusion encoding gradient pulse. This was not the case at the short time limit, where D(t)/D(0) was noticeably affected by the gas pressure in the sample. Increasing the gas pressure, and hence reducing D(0) and the diffusion during the gradient pulse served to reduce the previously observed deviation of D(t)/D(0) from the S/V(p) relation. The Pade approximation is used to interpolate between the long and short time limits in D(t). While the short time D(t) points lay above the interpolation line in the case of small beads, due to diffusion during the gradient pulse on the order of the pore size, it was also noted that the experimental D(t) data fell below the Pade line in the case of large beads, most likely due to finite size effects.

Non-NASA Center↗

Estimation of nonlinear pilot model parameters including time delay.

Investigation of the feasibility of using a Kalman filter estimator for the identification of unknown parameters in nonlinear dynamic systems with a time delay. The problem considered is the application of estimation theory to determine the parameters of a family of pilot models containing delayed states. In particular, the pilot-plant dynamics are described by differential-difference equations of the retarded type. The pilot delay, included as one of the unknown parameters to be determined, is kept in pure form as opposed to the Pade approximations generally used for these systems. Problem areas associated with processing real pilot response data are included in the discussion.

Schiess, J. R.↗

Recent results of nonlinear estimators applied to hereditary systems.

An application of the extended Kalman filter to delayed systems to estimate the state and time delay is presented. Two nonlinear estimators are discussed and the results compared with those of the Kalman filter. For all the filters considered, the hereditary system was treated with the delay in the pure form and by using Pade approximations of the delay. A summary of the convergence properties of the filters studied is given. The results indicate that the linear filter applied to the delayed system performs inadequately while the nonlinear filters provide reasonable estimates of both the state and the parameters.

Schiess, J. R.↗

A comparison of computational methods and algorithms for the complex gamma function

A survey and comparison of some computational methods and algorithms for gamma and log-gamma functions of complex arguments are presented. Methods and algorithms reported include Chebyshev approximations, Pade expansion and Stirling's asymptotic series. The comparison leads to the conclusion that Algorithm 421 published in the Communications of ACM by H. Kuki is the best program either for individual application or for the inclusion in subroutine libraries.

Ng, E. W.↗

Computer extension and analytic continuation of Stokes' expansion for gravity waves

Stokes' infinitesimal-wave expansion for steady progressive free-surface waves has been extended to high order using a computer to perform the coefficient arithmetic. Stokes' expansion has been found to be incapable of yielding the highest wave for any value of the water depth since convergence is limited by a square-root branch-point some distance short of the maximum. By reformulating the problem using a different independent parameter, the highest waves are obtained correctly. Series summation and analytic continuation are facilitated by the use of Pade approximants. The method is valid in principle for any finite value of the wavelength and solutions of high accuracy can be obtained for most values of the wave height and water depth.

Schwartz, L. W.↗

Hypersonic flows generated by parabolic and paraboloidal shock waves

A computer algorithm has been developed to determine the blunt-body flowfields supporting symmetric parabolic and paraboloidal shock waves at infinite free-stream Mach number. Solutions are expressed in an analytic form as high-order power series, in the coordinate normal to the shock, whose coefficients can be determined exactly. Analytic continuation is provided by the use of Pade approximations. Test cases provide solutions of very high accuracy. In the axisymmetric case for gamma equals 715 the solution has been found far downstream, where it agrees with the modified blast-wave results. For plane flow, on the other hand, a limit line appears within the shock layer, a short distance past the sonic line, suggesting the presence of an imbedded shock. Local solutions in the downstream limit are discussed.

Schwartz, L. W.↗

A technique for accelerating iterative convergence in numerical integration, with application in transonic aerodynamics

A technique is described for the efficient numerical solution of nonlinear partial differential equations by rapid iteration. In particular, a special approach is described for applying the Aitken acceleration formula (a simple Pade approximant) for accelerating the iterative convergence. The method finds the most appropriate successive approximations, which are in a most nearly geometric sequence, for use in the Aitken formula. Simple examples are given to illustrate the use of the method. The method is then applied to the mixed elliptic-hyperbolic problem of steady, inviscid, transonic flow over an airfoil in a subsonic free stream.

Martin, E. D.↗

A technique for accelerating iterative convergence in numerical integration, with application in transonic aerodynamics

A technique is described for the efficient numerical solution of nonlinear partial differential equations by rapid iteration. In particular, a special approach is described for applying the Aitken acceleration formula (a simple Pade approximant) for accelerating the iterative convergence. The method finds the most appropriate successive approximations, which are in a most nearly geometric sequence, for use in the Aitken formula. Simple examples are given to illustrate the use of the method. The method is then applied to the mixed elliptic-hyperbolic problem of steady, inviscid, transonic flow over an airfoil in a subsonic free stream.

Martin, E. D.↗