Search NASASearch

SEARCH · Search NASA

Results for “Taylor series”

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 19 records

High Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

Taylor series integration is implemented in a spacecraft trajectory analysis code-the Spacecraft N-body Analysis Program (SNAP) - and compared with the code s existing eighth-order Runge-Kutta Fehlberg time integration scheme. Nine trajectory problems, including near Earth, lunar, Mars and Europa missions, are analyzed. Head-to-head comparison at five different error tolerances shows that, on average, Taylor series is faster than Runge-Kutta Fehlberg by a factor of 15.8. Results further show that Taylor series has superior convergence properties. Taylor series integration proves that it can provide rapid, highly accurate solutions to spacecraft trajectory problems.

Scott, James R.

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.

High-Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

It has been known for some time that Taylor series (TS) integration is among the most efficient and accurate numerical methods in solving differential equations. However, the full benefit of the method has yet to be realized in calculating spacecraft trajectories, for two main reasons. First, most applications of Taylor series to trajectory propagation have focused on relatively simple problems of orbital motion or on specific problems and have not provided general applicability. Second, applications that have been more general have required use of a preprocessor, which inevitably imposes constraints on computational efficiency. The latter approach includes the work of Berryman et al., who solved the planetary n-body problem with relativistic effects. Their work specifically noted the computational inefficiencies arising from use of a preprocessor and pointed out the potential benefit of manually coding derivative routines. In this Engineering Note, we report on a systematic effort to directly implement Taylor series integration in an operational trajectory propagation code: the Spacecraft N-Body Analysis Program (SNAP). The present Taylor series implementation is unique in that it applies to spacecraft virtually anywhere in the solar system and can be used interchangeably with another integration method. SNAP is a high-fidelity trajectory propagator that includes force models for central body gravitation with N X N harmonics, other body gravitation with N X N harmonics, solar radiation pressure, atmospheric drag (for Earth orbits), and spacecraft thrusting (including shadowing). The governing equations are solved using an eighth-order Runge-Kutta Fehlberg (RKF) single-step method with variable step size control. In the present effort, TS is implemented by way of highly integrated subroutines that can be used interchangeably with RKF. This makes it possible to turn TS on or off during various phases of a mission. Current TS force models include central body gravitation with the J2 spherical harmonic, other body gravitation, thrust, constant atmospheric drag from Earth's atmosphere, and solar radiation pressure for a sphere under constant illumination. The purpose of this Engineering Note is to demonstrate the performance of TS integration in an operational trajectory analysis code and to compare it with a standard method, eighth-order RKF. Results show that TS is 16.6 times faster on average and is more accurate in 87.5% of the cases presented.

Scott, James R.

Two-point Taylor series expansions

Coefficients calculated for Taylor series expansion about two points - application of Taylor expansion to two-body problem

TWO-BODY PROBLEM

The numerical solution of ordinary differential equations by the Taylor series method

A programming implementation of the Taylor series method is presented for solving ordinary differential equations. The compiler is written in PL/1, and the target language is FORTRAN IV. The reduction of a differential system to rational form is described along with the procedures required for automatic numerical integration. The Taylor method is compared with two other methods for a number of differential equations. Algorithms using the Taylor method to find the zeroes of a given differential equation and to evaluate partial derivatives are presented. An annotated listing of the PL/1 program which performs the reduction and code generation is given. Listings of the FORTRAN routines used by the Taylor series method are included along with a compilation of all the recurrence formulas used to generate the Taylor coefficients for non-rational functions.

Silver, A. H.

Taylor Series Trajectory Calculations Including Oblateness Effects and Variable Atmospheric Density

Taylor series integration is implemented in NASA Glenn's Spacecraft N-body Analysis Program, and compared head-to-head with the code's existing 8th- order Runge-Kutta Fehlberg time integration scheme. This paper focuses on trajectory problems that include oblateness and/or variable atmospheric density. Taylor series is shown to be significantly faster and more accurate for oblateness problems up through a 4x4 field, with speedups ranging from a factor of 2 to 13. For problems with variable atmospheric density, speedups average 24 for atmospheric density alone, and average 1.6 to 8.2 when density and oblateness are combined.

Scott, James R.

Solving ODE Initial Value Problems With Implicit Taylor Series Methods

In this paper we introduce a new class of numerical methods for integrating ODE initial value problems. Specifically, we propose an extension of the Taylor series method which significantly improves its accuracy and stability while also increasing its range of applicability. To advance the solution from t (sub n) to t (sub n+1), we expand a series about the intermediate point t (sub n+mu):=t (sub n) + mu h, where h is the stepsize and mu is an arbitrary parameter called an expansion coefficient. We show that, in general, a Taylor series of degree k has exactly k expansion coefficients which raise its order of accuracy. The accuracy is raised by one order if k is odd, and by two orders if k is even. In addition, if k is three or greater, local extrapolation can be used to raise the accuracy two additional orders. We also examine stability for the problem y'= lambda y, Re (lambda) less than 0, and identify several A-stable schemes. Numerical results are presented for both fixed and variable stepsizes. It is shown that implicit Taylor series methods provide an effective integration tool for most problems, including stiff systems and ODE's with a singular point.

Scott, James R.

The ATOMFT integrator - Using Taylor series to solve ordinary differential equations

This paper discusses the application of ATOMFT, an integration package based on Taylor series solution with a sophisticated user interface. ATOMFT has the capabilities to allow the implementation of user defined functions and the solution of stiff and algebraic equations. Detailed examples, including the solutions to several astrodynamics problems, are presented. Comparisons with its predecessor ATOMCC and other modern integrators indicate that ATOMFT is a fast, accurate, and easy method to use to solve many differential equation problems.

Berryman, Kenneth W.

On Taylor-Series Approximations of Residual Stress

Although subgrid-scale models of similarity type are insufficiently dissipative for practical applications to large-eddy simulation, in recently published a priori analyses, they perform remarkably well in the sense of correlating highly against exact residual stresses. Here, Taylor-series expansions of residual stress are exploited to explain the observed behavior and "success" of similarity models. Until very recently, little attention has been given to issues related to the convergence of such expansions. Here, we re-express the convergence criterion of Vasilyev [J. Comput. Phys., 146 (1998)] in terms of the transfer function and the wavenumber cutoff of the grid filter.

Pruett, C. David

Application of AWE Along with a Combined FEM/MoM Technique to Compute RCS of a Cavity-Backed Aperture in an Infinite Ground Plane Over a Frequency Range

A hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique in conjunction with the Asymptotic Waveform Evaluation (AWE) technique is applied to obtain radar cross section (RCS) of a cavity-backed aperture in an infinite ground plane over a frequency range. The hybrid FEM/MoM technique when applied to the cavity-backed aperture results in an integro-differential equation with electric field as the unknown variable, the electric field obtained from the solution of the integro-differential equation is expanded in Taylor series. The coefficients of the Taylor series are obtained using the frequency derivatives of the integro-differential equation formed by the hybrid FEM/MoM technique. The series is then matched via the Pade approximation to a rational polynomial, which can be used to extrapolate the electric field over a frequency range. The RCS of the cavity-backed aperture is calculated using the electric field at different frequencies. Numerical results for a rectangular cavity, a circular cavity, and a material filled cavity are presented over a frequency range. Good agreement between AWE and the exact solution over the frequency range is obtained.

Reddy, C.J.

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

Application of Taylor's series to trajectory propagation

This paper describes the propagation of trajectories by the application of the preprocessor ATOMCC which uses Taylor's series to solve initial value problems in ordinary differential equations. Comparison of the results obtained with those from other methods are presented. The current studies indicate that the ATOMCC preprocessor is an easy, yet fast and accurate method for generating trajectories.

Stanford, R. H.

Approximate techniques of structural reanalysis

A study is made of two approximate techniques for structural reanalysis. These include Taylor series expansions for response variables in terms of design variables and the reduced-basis method. In addition, modifications to these techniques are proposed to overcome some of their major drawbacks. The modifications include a rational approach to the selection of the reduced-basis vectors and the use of Taylor series approximation in an iterative process. For the reduced basis a normalized set of vectors is chosen which consists of the original analyzed design and the first-order sensitivity analysis vectors. The use of the Taylor series approximation as a first (initial) estimate in an iterative process, can lead to significant improvements in accuracy, even with one iteration cycle. Therefore, the range of applicability of the reanalysis technique can be extended. Numerical examples are presented which demonstrate the gain in accuracy obtained by using the proposed modification techniques, for a wide range of variations in the design variables.

Noor, A. K.

Frequency Response Calculations of Input Characteristics of Cavity-Backed Aperture Antennas Using AWE with Hybrid FEM/MoM Technique

Application of Asymptotic Waveform Evaluation (AWE) is presented in conjunction with a hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique to calculate the input characteristics of cavity-backed aperture antennas over a frequency range. The hybrid FEM/MoM technique is used to form an integro-partial-differential equation to compute the electric field distribution of the cavity-backed aperture antenna. The electric field, thus obtained, is expanded in a Taylor series around the frequency of interest. The coefficients of 'Taylor series (called 'moments') are obtained using the frequency derivatives of the integro-partial-differential Equation formed by the hybrid FEM/MoM technique. Using the moments, the electric field in the cavity is obtained over a frequency range. Using the electric field at different frequencies, the input characteristics of the antenna are obtained over a wide frequency band. Numerical results for an open coaxial line, probe fed cavity, and cavity-backed microstrip patch antennas are presented. Good agreement between AWE and the exact solution over the frequency range is observed.

Reddy, C. J.

Application of AWE for RCS Frequency Response Calculations Using Method of Moments

An implementation of the Asymptotic Waveform Evaluation (AWE) technique is presented for obtaining the frequency response of the Radar Cross Section (RCS) of arbitrarily shaped, three-dimensional perfect electric conductor (PEC) bodies. An Electric Field Integral Equation (EFIE) is solved using the Method of Moments (MoM) to compute the RCS. The electric current, thus obtained, is expanded in a Taylor series around the frequency of interest. The coefficients of the Taylor series (called 'moments') are obtained using the frequency derivatives of the EFIE. Using the moments, the electric current on the PEC body is obtained over a frequency band. Using the electric current at different frequencies, RCS of the PEC body is obtained over a wide frequency band. Numerical results for a square plate, a cube, and a sphere are presented over a bandwidth. A good agreement between AWE and the exact solution over the bandwidth is observed.

Reddy, C. J.

TSEUQLib

TSEUQLib (Taylor Series Expansion Uncertainty Quantification Library): Conduct UQ and sensitivity analysis with the Taylor series expansion surrogate model.

Balcer, Matthew

Structural Reliability Analysis and Optimization: Use of Approximations

This report is intended for the demonstration of function approximation concepts and their applicability in reliability analysis and design. Particularly, approximations in the calculation of the safety index, failure probability and structural optimization (modification of design variables) are developed. With this scope in mind, extensive details on probability theory are avoided. Definitions relevant to the stated objectives have been taken from standard text books. The idea of function approximations is to minimize the repetitive use of computationally intensive calculations by replacing them with simpler closed-form equations, which could be nonlinear. Typically, the approximations provide good accuracy around the points where they are constructed, and they need to be periodically updated to extend their utility. There are approximations in calculating the failure probability of a limit state function. The first one, which is most commonly discussed, is how the limit state is approximated at the design point. Most of the time this could be a first-order Taylor series expansion, also known as the First Order Reliability Method (FORM), or a second-order Taylor series expansion (paraboloid), also known as the Second Order Reliability Method (SORM). From the computational procedure point of view, this step comes after the design point identification; however, the order of approximation for the probability of failure calculation is discussed first, and it is denoted by either FORM or SORM. The other approximation of interest is how the design point, or the most probable failure point (MPP), is identified. For iteratively finding this point, again the limit state is approximated. The accuracy and efficiency of the approximations make the search process quite practical for analysis intensive approaches such as the finite element methods; therefore, the crux of this research is to develop excellent approximations for MPP identification and also different approximations including the higher-order reliability methods (HORM) for representing the failure surface. This report is divided into several parts to emphasize different segments of the structural reliability analysis and design. Broadly, it consists of mathematical foundations, methods and applications. Chapter I discusses the fundamental definitions of the probability theory, which are mostly available in standard text books. Probability density function descriptions relevant to this work are addressed. In Chapter 2, the concept and utility of function approximation are discussed for a general application in engineering analysis. Various forms of function representations and the latest developments in nonlinear adaptive approximations are presented with comparison studies. Research work accomplished in reliability analysis is presented in Chapter 3. First, the definition of safety index and most probable point of failure are introduced. Efficient ways of computing the safety index with a fewer number of iterations is emphasized. In chapter 4, the probability of failure prediction is presented using first-order, second-order and higher-order methods. System reliability methods are discussed in chapter 5. Chapter 6 presents optimization techniques for the modification and redistribution of structural sizes for improving the structural reliability. The report also contains several appendices on probability parameters.

Grandhi, Ramana V.