Search NASASearch

SEARCH · Search NASA

Results for “High-order methods”

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-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster

On the practical use of high-order methods for hyperbolic systems

The paper tests a number of high order methods on a variety of dynamic problems in one, two, and three space dimensions. The problems covered include wave propagation phenomena as well as an asymptotic approach to a steady state. Consideration is given to both smooth and shocked flows. It is shown that the methods compared require only minor modifications of many existing second-order schemes. Further, the results show that significant gains can be expected from the use of fourth-order methods. Finally, spectral methods are also considered for some of the problems presented.

Turkel, E.

Higher-Order Panel Method for Aerodynamic Flow Analysis

PANAIR uses high-order panel method to predict inviscid subsonic or supersonic flows about arbitrary configuration. Panel method solves linear partial differential equation numerically by approximating configuration surface with panels on which unknown "singularity strengths" are defined. PANAIR includes advanced software technology as well as advanced aerodynamic technology.

Erickson, L.

Comparison of numerical integration techniques for orbital applications

The present work gives a brief comparison of the performance of programs for integrating differential equations for orbital applications. The evaluation criteria and the method of testing are described, and the results of the test problem set are included. Integration methods that were chosen for comparison include high-order Runge-Kutta methods; a rational extrapolation method (Bulirsch and Stoer, 1966); a variable step, variable order, multistep method (Krogh, 1969); classical multistep methods of Adams and Cowell; and modified multistep methods. The high-order Runge-Kutta methods used in the comparison include RKF 7(8) and RKF 8(9) (Fehlberg, 1968), and RKS 8-10 (Shanks, 1966).

Moore, H.

A method for obtaining reduced-order control laws for high-order systems using optimization techniques

A method of synthesizing reduced-order optimal feedback control laws for a high-order system is developed. A nonlinear programming algorithm is employed to search for the control law design variables that minimize a performance index defined by a weighted sum of mean-square steady-state responses and control inputs. An analogy with the linear quadractic Gaussian solution is utilized to select a set of design variables and their initial values. To improve the stability margins of the system, an input-noise adjustment procedure is used in the design algorithm. The method is applied to the synthesis of an active flutter-suppression control law for a wind tunnel model of an aeroelastic wing. The reduced-order controller is compared with the corresponding full-order controller and found to provide nearly optimal performance. The performance of the present method appeared to be superior to that of two other control law order-reduction methods. It is concluded that by using the present algorithm, nearly optimal low-order control laws with good stability margins can be synthesized.

Mukhopadhyay, V.

Multistep methods of numerical integration using back-corrections

A class of linear multistep methods is proposed for the solution of the equations of motion of certain dynamical systems encountered in celestial mechanics and astrodynamics. These methods are distinguished from the classical predictor-corrector methods in that they permit 'back-corrections' of the solution to be made. As the integration advances in time, the numerical solution is corrected or improved at certain points in the past. The enhanced numerical stability of these methods allows the meaningful application of high-order algorithms. Consequently, step sizes larger than those attainable with the classical methods may be adopted, and greater overall efficiency may be realized. These methods are applied to the problem of determining the orbit of an artificial satellite, and the results are compared with those obtained using classical methods.

Feagin, T.

Evaluation of the 'mean frequency' technique

The 'mean frequency' technique, a simple procedure introduced by Bebb and Gold for the approximate evaluation of sums occurring in high-order perturbation theory, represents a useful approximation method. Its predictions compare favorably to exact results obtained by Gontier and Trahin for multiphoton bound-bound transitions in hydrogen. However, the technique can be in error if the 'mean frequency' lies near certain integers.

Gaddy, E. M.

On the improvement of deconvolution with digitized data using a Runge-Kutta integration scheme

A relatively simple change in the treatment of the input function in numerical integration of high-order differential equations by Runge-Kutta methods provides substantial improvements in accuracy, particularly when the forcing function is in digitized form. The Runge-Kutta-Gill coefficients are modified to incorporate the changes; with pulse-type excitations, improvements on the order of 2 to 50 times greater accuracy are demonstrated.

Houghton, J. R.

Modeling with Liapunov functions.

Behavior of high-order linear control systems analyzed using liapunov second method, by finding low-order model with closely approximate response behavior of high order linear control systems analyzed, using Liapunov second method, by finding low order model with closely approximate response

LIAPUNOV FUNCTION

Determination of pole sensitivities by Danilevskii's method

In control theory, a synonymous term for pole sensitivity is eigenvalue sensitivity. Existing methods of calculating eigenvalues are cumbersome, and cannot be trusted for systems roughly greater than tenth order. The method proposed in the present paper is applicable to high-order system. (It has been routinely used to generate eigenvalue sensitivities for systems up to 26th order, using a UNIVAC 1106.) Danilevskii's method is shown to be suitable for performing the necessary evaluations. The result is a rational function that can be used to evaluate the sensitivities for all distinct poles.

Nail, J. B.

Transfer function verification and block diagram simplification of a very high-order distributed pole closed-loop servo by means of non-linear time-response simulation

Linear frequency domain methods are inadequate in analyzing the 1975 Viking Orbiter (VO75) digital tape recorder servo due to dominant nonlinear effects such as servo signal limiting, unidirectional servo control, and static/dynamic Coulomb friction. The frequency loop (speed control) servo of the VO75 tape recorder is used to illustrate the analytical tools and methodology of system redundancy elimination and high order transfer function verification. The paper compares time-domain performance parameters derived from a series of nonlinear time responses with the available experimental data in order to select the best possible analytical transfer function representation of the tape transport (mechanical segment of the tape recorder) from several possible candidates. The study also shows how an analytical time-response simulation taking into account most system nonlinearities can pinpoint system redundancy and overdesign stemming from a strictly empirical design approach. System order reduction is achieved through truncation of individual transfer functions and elimination of redundant blocks.

Mukhopadhyay, A. K.

Stabilization of large-scale systems - A spinning flexible spacecraft

Stabilization of high-order constant systems by multilevel control is proposed in the framework of the decomposition-aggregation method for stability analysis of large-scale systems. An 11-order linear model of a spinning flexible spacecraft is stabilized by the method using both the passive and the active stabilization devices. The method can take advantage of the special structural features of the model to provide the best estimate of the stability region for an important parameter which couples the wobble and the spin motion of the spacecraft.

Siljak, D. D.

The gravity field in the central Pacific from satellite-to-satellite tracking

Satellite-to-satellite Doppler tracking between the ATS 6 and the GEOS 3 spacecraft was used to measure the high-degree and high-order gravity field over an 80-deg region in the central Pacific Ocean. Forty passes of GEOS 3/ATS 6 Doppler data have been analyzed. The precision of these range rate data is about 0.3 mm/s, and the line-of-sight gravity anomalies recovered from these data have a precision of about 0.2 mGal at the GEOS 3 altitude of about 840 km. In general, the agreement between the SST-derived map and the conventional GEM method and an altimeter-derived geoid is good. Eight significant positive gravity anomalies were exposed in the central Pacific. Generally speaking, the anomalies form a roughly east-west pattern of alternating sign in the central region, and near the East Pacific they strike about north and south.

Marsh, J. G.

Discrete control of linear distributed systems with application to the deformable primary mirror of a large orbiting telescope

One of the more significant technological problems associated with the orbital operation of large astronomical telescope's is the fabrication and maintenance of the primary mirror surface to the tolerance required for diffraction-limited performance. An interesting approach to the solution of this problem involves continuously measuring and automatically correcting the optical surface of a thin deformable mirror by means of discrete actuators located on its rear surface: The realization of diffraction-limited performance from a telescope in space by this method rests on the ability of the designer to achieve extremely accurate control of a highly complex, interacting, multivariable system. This paper presents the results of a detailed study of the discrete control of linear distributed systems with specific application to the design of a practical controller for a plant representative of a telescope primary mirror for an orbiting astronomical observatory. The problem of controlling the distributed plant is treated by employing modal techniques to represent variations in the optical figure. Distortion of the mirror surface, which arises primarily from thermal gradients, is countered by actuators working against a backing structure to apply a corrective force distribution to the controlled surface. Each displacement actuator is in series with a spring attached to the mirror by means of a pad intentionally introduced to restrict the excitation of high-order modes. Control is then exerted over a finite number (equal.to the number of actuators) of the most significant modes. Through the application of the modal expansion technique the mirror equation of motion is transformed tb a set of uncoupled, linear, time-invariant, ordinary differential equations. The desired dynamic response and static accuracy may then be achieved by the application of classical single-variable design techniques. The formulation of a quadratic performance index which incorporates a measure of image quality permits determination of the trade-off between the-number of actuators and optical purity. A criterion for defining actuator placement and pad size is presented which minimizes the tendency of the controller to excite the unmonitored modes.

Jeremiah F Creedon

High-order numerical solutions using cubic splines

The cubic spline collocation procedure for the numerical solution of partial differential equations was reformulated so that the accuracy of the second-derivative approximation is improved and parallels that previously obtained for lower derivative terms. The final result is a numerical procedure having overall third-order accuracy for a nonuniform mesh and overall fourth-order accuracy for a uniform mesh. Application of the technique was made to the Burger's equation, to the flow around a linear corner, to the potential flow over a circular cylinder, and to boundary layer problems. The results confirmed the higher-order accuracy of the spline method and suggest that accurate solutions for more practical flow problems can be obtained with relatively coarse nonuniform meshes.

Rubin, S. G.