Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time discretization”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

A mixed pseudospectral/finite difference method for the axisymmetric flow in a heated, rotating spherical shell

For a Spacelab flight, a model experiment of the earth's atmospheric circulation has been proposed. This experiment is known as the Atmospheric General Circulation Experiment (AGCE). In the experiment concentric spheres will rotate as a solid body, while a dielectric fluid is confined in a portion of the gap between the spheres. A zero gravity environment will be required in the context of the simulation of the gravitational body force on the atmosphere. The present study is concerned with the development of pseudospectral/finite difference (PS/FD) model and its subsequent application to physical cases relevant to the AGCE. The model is based on a hybrid scheme involving a pseudospectral latitudinal formulation, and finite difference radial and time discretization. The advantages of the use of the hybrid PS/FD method compared to a pure second-order accurate finite difference (FD) method are discussed, taking into account the higher accuracy and efficiency of the PS/FD method.

Macaraeg, M. G.↗
NASA NTRS1985

Topics in spectral methods

After detailing the construction of spectral approximations to time-dependent mixed initial boundary value problems, a study is conducted of differential equations of the form 'partial derivative of u/partial derivative of t = Lu + f', where for each t, u(t) belongs to a Hilbert space such that u satisfies homogeneous boundary conditions. For the sake of simplicity, it is assumed that L is an unbounded, time-independent linear operator. Attention is given to Fourier methods of both Galerkin and pseudospectral method types, the Galerkin method, the pseudospectral Chebyshev and Legendre methods, the error equation, hyperbolic partial differentiation equations, and time discretization and iterative methods.

Gottlieb, D.↗
NASA NTRS1986

Linear-quadratic stationkeeping for the STS Orbiter

An automatic controller for Space Shuttle Orbiter stationkeeping or formationkeeping is presented. The controller uses a feed-forward loop for precise control at non-equilibrium set points, and a discrete-time linear quadratic regulator for disturbance rejection. Selection of sample time and quadratic cost matrices is done to provide good limit cycling performance in the presence of sensor noise and a control threshold nonlinearity. Large position changes are effected by feed-forward of a minimum fuel or straightline maneuver path. Results show that 20 ft precision can be achieved using current sensors, and that 2 ft precision would be possible with enhanced sensors.

Redding, D. C.↗
NASA NTRS1986

A comparison of methods for DPLL loop filter design

Four design methodologies for loop filters for a class of digital phase-locked loops (DPLLs) are presented. The first design maps an optimum analog filter into the digital domain; the second approach designs a filter that minimizes in discrete time weighted combination of the variance of the phase error due to noise and the sum square of the deterministic phase error component; the third method uses Kalman filter estimation theory to design a filter composed of a least squares fading memory estimator and a predictor. The last design relies on classical theory, including rules for the design of compensators. Linear analysis is used throughout the article to compare different designs, and includes stability, steady state performance and transient behavior of the loops. Design methodology is not critical when the loop update rate can be made high relative to loop bandwidth, as the performance approaches that of continuous time. For low update rates, however, the miminization method is significantly superior to the other methods.

Aguirre, S.↗
NASA NTRS1987

Efficient implementation of essentially non-oscillatory shock capturing schemes

In the computation of discontinuous solutions of hyperbolic conservation laws, TVD (total-variation-diminishing), TVB (total-variation-bounded) and the recently developed ENO (essentially non-oscillatory) schemes have proven to be very useful. In this paper two improvements are discussed: a simple TVD Runge-Kutta type time discretization, and an ENO construction procedure based on fluxes rather than on cell averages. These improvements simplify the schemes considerably -- especially for multi-dimensional problems or problems with forcing terms. Preliminary numerical results are also given.

Shu, Chi-Wang↗
NASA NTRS1987

Analysis of the Continuous Stellar Tracking Attitude Reference (CSTAR) attitude rate processor

The Continuous Stellar Tracking Attitude Reference (CSTAR) system is an in-house project for Space Station to provide high accuracy, drift free attitude and angular rate information for the GN&C system. The outputs of the solid state star trackers are processed to provide attitude information; rate data is then derived from the attitude. Rate derivation is based on discrete time polynomial approximation techniques. This gives simple algorithms which allow for interpolation by other users. Attitude rate is modeled as a constant with low amplitude, low frequency sinusoids superimposed. The rate processor is parameterized to account for the effects of random errors, sample rate, data processing rate and perturbation frequency. The baseline system may be characterized as follows: the three sigma attitude accuracy is 0.01 degrees, the three sigma rate accuracy is 0.0001 degrees per second, the sample rate is 100 Hertz, the sampled signal is bandlimited to 0.5 Hertz, and the data processing rate is 10 Hertz. The above system requires a differentiator of length 127. This will track rate perturbations of frequencies less than 0.01 Hertz with low systematic errors.

Uhde-Lacovara, J.↗
NASA NTRS1987

The dynamics and control of large flexible space structures X, part 1

The effect of delay in the control system input on the stability of a continuously acting controller which is designed without considering the delay is studied. The stability analysis of a second order plant is studied analytically and verified numerically. For this example it is found that the system becomes unstable for a delay which is equivalent to only 16 percent of its natural period of motion. It is also observed that even a small amount of natural damping in the system can increase the amount of delay that can be tolerated before the onset of instability. The delay problem is formulated in the discrete time domain and an analysis procedure suggested. The maximum principle from optimal control theory is applied to minimize the time required for the slewing of a general rigid spacecraft. The slewing motion need not be restricted to a single axis maneuver. The minimum slewing time is calculated based on a quasi-linearization algorithm for the resulting two point boundary value problem. Numerical examples based on the rigidized in-orbit model of the SCOLE also include the more general reflector line-of-sight slewing maneuvers.

Bainum, Peter M.↗
NASA NTRS1986

Single step optimization of feedback-decoupled collision avoidance manipulator maneuvers

Simultaneous robot path planning and path following is shown to be achievable in the presence of motor saturation and obstacle avoidance requirements. The discrete time algorithm derived performs one step ahead mean square optimization of commanded joint accelerations, subject to present actuator force or torque constraints and N step ahead prediction of configuration constraints.

Chen, N.↗
NASA NTRS1986

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid points are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, Hillel↗
NASA NTRS1986

Discrete Dirac equation on a finite half-integer lattice

The formulation of the Dirac equation on a discrete lattice with half-integer spacing and periodic boundary conditions is investigated analytically. The importance of lattice formulations for problems in field theory and quantum mechanics is explained; the concept of half-integer Fourier representation is introduced; the discrete Dirac equation for the two-dimensional case is derived; dispersion relations for the four-dimensional case are developed; and the spinor formulation for the Dirac fields on the half-integer lattice and the discrete time variable for the four-dimensional time-dependent Dirac equation are obtained. It is argued that the half-integer lattice, because it takes the Dirac Lagrangian into account, is more than a mere relabeling of the integer lattice and may have fundamental physical meaning (e.g., for the statistics of fermions). It is noted that the present formulation does not lead to species doubling, except in the continuum limit.

Smalley, L. L.↗
NASA NTRS1987

A modern control theory based algorithm for control of the NASA/JPL 70-meter antenna axis servos

A digital computer-based state variable controller was designed and applied to the 70-m antenna axis servos. The general equations and structure of the algorithm and provisions for alternate position error feedback modes to accommodate intertarget slew, encoder referenced tracking, and precision tracking modes are descibed. Development of the discrete time domain control model and computation of estimator and control gain parameters based on closed loop pole placement criteria are discussed. The new algorithm was successfully implemented and tested in the 70-m antenna at Deep Space Network station 63 in Spain.

Hill, R. E.↗
NASA NTRS1988

Efficient implementation of essentially non-oscillatory shock capturing schemes, 2

Earlier work on the efficient implementation of ENO (essentially non-oscillatory) shock capturing schemes is continued. A new simplified expression is provided for the ENO construction procedure based again on numerical fluxes rather than cell averages. Also considered are two improvements which are labeled ENO-LLF (local Lax-Friedrichs) and ENO-Roe, which yield sharper shock transitions, improved overall efficiency, and lower computational cost than previous implementation of the ENO schemes. Two methods of sharpening contact discontinuities, i.e., the subcell resolution idea of Harten and the artificial compression idea of Yang, which those authors used originally in the cell-average framework, are supplied to the current ENO schemes using numerical fluxes and TVD Runge-Kutta time discretizations. The implementation for nonlinear systems and multi-dimensions is given. Finally, many numerical examples, including a compressible shock turbulence interaction flow calculation, are given.

Shu, Chi-Wang↗
NASA NTRS1987

Fast simulation techniques for switching converters

Techniques for simulating a switching converter are examined. The state equations for the equivalent circuits, which represent the switching converter, are presented and explained. The uses of the Newton-Raphson iteration, low ripple approximation, half-cycle symmetry, and discrete time equations to compute the interval durations are described. An example is presented in which these methods are illustrated by applying them to a parallel-loaded resonant inverter with three equivalent circuits for its continuous mode of operation.

King, Roger J.↗
NASA NTRS1987

Linear quadratic stationkeeping on travelling ellipses

An automatic controller for Space Shuttle stationkeeping or formationkeeping is presented. The controller, which is a candidate for future on-orbit autopilot enhancement, uses a fuel 'optimized' limit cycle trajectory in a feedforward loop for precise control at nonequilibrium set points. Disturbance rejection is accomplished by a discrete time, linear quadratic regulator, which is employed for tracking of the feedforward model. Feedback gain selection is done to provide good limit cycling performance and low fuel consumption. Velocity correlations are carried out by use of a pseudo 6 degree-of-freedom jet selection scheme. Results indicate that 20 ft precision can be achieved using current sensors.

Adams, Neil J.↗
NASA NTRS1988

Generating parity relations for detecting and identifying control system component failures

The monitoring of control system sensors and actuators for failures is presently undertaken by means of an exceptionally simple form of generalized parity relations based on a discrete-time model of the dynamics of linear, time-invariant systems. These generalized parity relations are constructed by recourse to a transfer matrix-description of the system that is additionally useful in the interpretation of their properties. Attention is given to a novel method for constructing the parity relation of minimum length that depends on the output of only a single sensor.

Vander Velder, Wallace E.↗
NASA NTRS1987

Adaptive identification and control for manipulators without using joint accelerations

The authors present a new scheme for the adaptive control of mechanical manipulators along with proof of convergence. Difficulties encountered in the input-error type adaptive controller are remedied by combining an explicit identifier with discrete time updates of the parameter values along with a computed torque type control law. The new scheme requires no matrix inversion with its consequent drain on the computational power of the controller. A discussion is presented of the dynamic model of a manipulator, a computed torque adaptive scheme, an input error identifier, and a hybrid identifier.

Hsu, Ping↗
NASA NTRS1988

Vectorization of a particle simulation method for hypersonic rarefied flow

An efficient particle simulation technique for hypersonic rarefied flows is presented at an algorithmic and implementation level. The implementation is for a vector computer architecture, specifically the Cray-2. The method models an ideal diatomic Maxwell molecule with three translational and two rotational degrees of freedom. Algorithms are designed specifically for compatibility with fine grain parallelism by reducing the number of data dependencies in the computation. By insisting on this compatibility, the method is capable of performing simulation on a much larger scale than previously possible. A two-dimensional simulation of supersonic flow over a wedge is carried out for the near-continuum limit where the gas is in equilibrium and the ideal solution can be used as a check on the accuracy of the gas model employed in the method. Also, a three-dimensional, Mach 8, rarefied flow about a finite-span flat plate at a 45 degree angle of attack was simulated. It utilized over 10 to the 7th particles carried through 400 discrete time steps in less than one hour of Cray-2 CPU time. This problem was chosen to exhibit the capability of the method in handling a large number of particles and a true three-dimensional geometry.

Mcdonald, Jeffrey D.↗
NASA NTRS1988

Efficient implementation of essentially non-oscillatory shock-capturing schemes

In the computation of discontinuous solutions of hyperbolic conservation laws, TVD (total-variation-diminishing), TVB (total-variation-bounded) and the recently developed ENO (essentially non-oscillatory) schemes have proven to be very useful. In this paper two improvements are discussed: a simple TVD Runge-Kutta type time discretization, and an ENO construction procedure based on fluxes rather than on cell averages. These improvements simplify the schemes considerably - especially for multi-dimensional problems or problems with forcing terms. Preliminary numerical results are also given.

Shu, Chi-Wang↗