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 37 records · Page 2

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.

Mixing control in a plane shear layer

An investigation of mixing processes in a plane shear layer by a direct numerical simulation with a high-order numerical technique, the isoparametric spectral element method, is presented. The computational domain includes a region with the splitter plate that allows to start the flow as a laminar (Blasius) boundary layers with different velocities and concentrations, merging in the shear layer and undergoing transition to highly-unsteady structures with intensive mixing. The augmentation of mixing as a function of a downstream location, the intensity and type of vortical structures in the layer is investigated by an introduction of disturbance at the high-velocity inflow. The character of the dynamic mixing is displayed along with time averaged profiles and statistical characteristics.

Korczak, K. Z.

Control of large flexible spacecraft by the independent modal-space control method

The problem of control of a large-order flexible structure in the form of a plate-like lattice by the Independent Modal-Space Control (IMSC) method is presented. The equations of motion are first transformed to the modal space, thus obtaining internal (plant) decoupling of the system. Then, the control laws are designed in the modal space for each mode separately, so that the modal equations of motion are rendered externally (controller) decoupled. This complete decoupling applies both to rigid-body modes and elastic modes. The application of linear optimal control, in conjunction with a quadratic performance index, is first reviewed. A solution for high-order systems is proposed here by the IMSC method, whereby the problem is reduced to a number of modal minimum-fuel problems for the controlled modes.

Meirovitch, L.

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

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Harten, A.

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.

Implicit flux-corrected transport scheme for low-speed flow

A new high-order nonlinear Eulerian convection algorithm, the barely implicit correction (BIC)-FCT method, is presented. The BIC-FCT makes it possible to perform accurate, multidimensional, slow-flow calculations in which the fluid expansion is important. Owing to its low cost even in two dimensions, fairly detailed chemistry models and other physical processes can be included. Areas of BIC-FCT applicability include premixed flames, diffusion flames, and turbulent jet flames.

Patnaik, G.

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.

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.

Harten, A.

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.

Application of a sensitivity analysis technique to high-order digital flight control systems

A sensitivity analysis technique for multiloop flight control systems is studied. This technique uses the scaled singular values of the return difference matrix as a measure of the relative stability of a control system. It then uses the gradients of these singular values with respect to system and controller parameters to judge sensitivity. The sensitivity analysis technique is first reviewed; then it is extended to include digital systems, through the derivation of singular-value gradient equations. Gradients with respect to parameters which do not appear explicitly as control-system matrix elements are also derived, so that high-order systems can be studied. A complete review of the integrated technique is given by way of a simple example: the inverted pendulum problem. The technique is then demonstrated on the X-29 control laws. Results show linear models of real systems can be analyzed by this sensitivity technique, if it is applied with care. A computer program called SVA was written to accomplish the singular-value sensitivity analysis techniques. Thus computational methods and considerations form an integral part of many of the discussions. A user's guide to the program is included. The SVA is a fully public domain program, running on the NASA/Dryden Elxsi computer.

Paduano, James 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.

The unique determination of density from high-order potentials

The existence of higher-order potentials is postulated, and that the density function can be determined from this set of potentials is demonstrated. Moreover, it is shown that there is a one-to-one correspondence between such potentials and the density function, provided it is continuous. Finally, it is noted that this demonstration includes a method of constructing the density function from the set of potentials as well as indicating how to approximate the density from incomplete knowledge of the potentials.

Johnson, G. G.

A spectral multi-domain technique with application to generalized curvilinear coordinates

Spectral collocation methods have proven to be efficient discretization schemes for many aerodynamic and fluid mechanic problems. The high order accuracy and resolution shown by these methods allows one to obtain engineering accuracy solutions on coarse meshes, or alternatively, to obtain solutions with very small error. One drawback to these techniques was the requirement that a complicated physical domain must map into a simple computational domain for discretization. This mapping must be smooth if the high order accuracy and expontential convergence rates associated with spectral methods are to be preserved. Additionally even smooth stretching transformations can decrease the accuracy of a spectral method, if the stretching is severe. A further difficulty with spectral methods was in their implementation on parallel processing computers, where efficient spectral algorithms were lacking. The above restrictions are overcome by splitting the domain into regions, each of which preserve the advantages of spectral collocation, and allow the ratio of the mesh spacing between regions to be several orders of magnitude higher than allowable in a single domain. Such stretchings would be required to resolve the thin viscous region in an external aerodynamic problem. Adjoining regions are interfaced by enforcing a global flux balance which preserves high-order continuity of the solution, regardless of the type of the equations being solved.

Macaraeg, M. G.

Minimum-fuel control of high-order systems

The minimum-fuel control problem is of special interest in various space systems. To date, solutions of minimum-fuel control problems have been carried out for relatively low-order systems. Space structures, however, are generally characterized by a large number of degrees of freedom, so that minimum-fuel control of such systems requires a new approach. In the independent modal-space control (IMSC) method, the control laws are designed in the modal space for each mode independently. The minimum-fuel problem reduces to that of a set of independent second-order systems, so that minimum-fuel control is possible. This paper shows how the IMSC method can be used to control a space structure with a minimum amount of fuel. A numerical example is presented.

Shenhar, J.