Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equation systems”

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 883 records · Page 49

Proportional and Integral Thermal Control System for Large Scale Heating Tests

The National Aeronautics and Space Administration Armstrong Flight Research Center (Edwards, California) Flight Loads Laboratory is a unique national laboratory that supports thermal, mechanical, thermal/mechanical, and structural dynamics research and testing. A Proportional Integral thermal control system was designed and implemented to support thermal tests. A thermal control algorithm supporting a quartz lamp heater was developed based on the Proportional Integral control concept and a linearized heating process. The thermal control equations were derived and expressed in terms of power levels, integral gain, proportional gain, and differences between thermal setpoints and skin temperatures. Besides the derived equations, user's predefined thermal test information generated in the form of thermal maps was used to implement the thermal control system capabilities. Graphite heater closed-loop thermal control and graphite heater open-loop power level were added later to fulfill the demand for higher temperature tests. Verification and validation tests were performed to ensure that the thermal control system requirements were achieved. This thermal control system has successfully supported many milestone thermal and thermal/mechanical tests for almost a decade with temperatures ranging from 50 F to 3000 F and temperature rise rates from -10 F/s to 70 F/s for a variety of test articles having unique thermal profiles and test setups.

Thermal Control System↗

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗

A three-dimensional dynamic analysis and libration study of a tethered satellites system

By a three-dimensional Lagrangian formulation, an exact set of ten coupled nonlinear second-order differential equations has been derived for a system with an extensible tether connecting two end satellites of distributed mass. The effects of tether mass, small orbital eccentricity, central body oblateness, aerodynamic drag force, and solar radiation pressure are also included in the formulation. By linearizing the exact differential equations, the in-plane (orbital plane) differential equations are found to be decoupled from the out-of-plane ones. The characteristic equation of the in-plane differential equations is derived and some associated stability constraints are shown.

Tang, C. C. H.↗

The stability of numerical methods for second order ordinary differential equations

An important characterization of a numerical method for first order ODE's is the region of absolute stability. If all eigenvalues of the linear problem dy/dt = Ay are inside this region, the numerical method is stable. If the second order system d/dt(dy/dt) = 2Ady/dt - By is solved as a first order system, the same result applies to the eigenvalues of the generalized eigenvalue problem (lambda-squared)I 2(lambda)A + B. No such region exists for general methods for second order equations, but in some cases a region of absolute stability can be defined for methods for the single second order equation d/dt(dy/dt) = 2ady/dt - by. The absence of a region of absolute stability can occur when different members of a system of first order equations are solved by different methods.

Gear, C. W.↗

Complex generalized minimal residual algorithm for iterative solution of quantum-mechanical reactive scattering equations

Complex dense matrices corresponding to the D + H2 and O + HD reactions were solved using a complex generalized minimal residual (GMRes) algorithm described by Saad and Schultz (1986) and Saad (1990). To provide a test case with a different structure, the H + H2 system was also considered. It is shown that the computational effort for solutions with the GMRes algorithm depends on the dimension of the linear system, the total energy of the scattering problem, and the accuracy criterion. In several cases with dimensions in the range 1110-5632, the GMRes algorithm outperformed the LAPACK direct solver, with speedups for the linear equation solution as large as a factor of 23.

Chatfield, David C.↗

Linear and nonlinear aspects of snow albedo feedbacks in atmospheric models

Namias' hypothesis, that anomalous snowcover on the eastern side of the North American continent can generate an anomalous east coast low pressure system and an anomalous inland high pressure system, is consistent with the time-averaged anomalous response from a nonlinear, primitive equation channel model with an idealized, flat land-sea arrangement. An attempt to understand and describe this anomalous response in the nonlinear model as a linear response to anomalous diabatic heating was largely unsuccessful, primarily because the anomalous eddy fluxes were also important. This unsuccessful attempt to describe the nonlinear model's time averages by linear theory then motivated several comparisons between linear and nonlinear severely truncated quasi-geostrophic models. It was also found in these models that the eddy fluxes were extremely important for forcing or dissipating the stationary eddies.

Roads, J. O.↗

Linearized theory of inhomogeneous multiple 'water-bag' plasmas

Equations are derived for describing the inhomogeneous equilibrium and small deviations from the equilibrium, giving particular attention to systems with trapped particles. An investigation is conducted of periodic systems with a single trapped-particle water bag, taking into account the behavior of the perturbation equations at the turning points. An outline is provided concerning a procedure for obtaining the eigenvalues. The results of stability calculations connected with the sideband effects are considered along with questions regarding the general applicability of the multiple water-bag approach in stability calculations.

Bloomberg, H. W.↗

Eigenmode analysis of unsteady one-dimensional Euler equations

The initial boundary value problem describing the evolution of unsteady linearized perturbations of a steady, uniform subsonic flow is analyzed. The eigenmodes and eigenfrequencies of the system are derived and several examples are presented to illustrate the effect of different boundary conditions on the exponential decay rate of the eigenmodes. The resultant implications for the stability and convergence rates of finite difference computations are discussed.

Giles, M.↗

Propagating Qualitative Values Through Quantitative Equations

In most practical problems where traditional numeric simulation is not adequate, one need to reason about a system with both qualitative and quantitative equations. In this paper, we address the problem of propagating qualitative values represented as interval values through quantitative equations. Previous research has produced exponential-time algorithms for approximate solution of the problem. These may not meet the stringent requirements of many real time applications. This paper advances the state of art by producing a linear-time algorithm that can propagate a qualitative value through a class of complex quantitative equations exactly and through arbitrary algebraic expressions approximately. The algorithm was found applicable to Space Shuttle Reaction Control System model.

Kulkarni, Deepak↗

A macroscopic plasma Lagrangian and its application to wave interactions and resonances

The derivation of a macroscopic plasma Lagrangian is considered, along with its application to the description of nonlinear three-wave interaction in a homogeneous plasma and linear resonance oscillations in a inhomogeneous plasma. One approach to obtain the Lagrangian is via the inverse problem of the calculus of variations for arbitrary first and second order quasilinear partial differential systems. Necessary and sufficient conditions for the given equations to be Euler-Lagrange equations of a Lagrangian are obtained. These conditions are then used to determine the transformations that convert some classes of non-Euler-Lagrange equations to Euler-Lagrange equation form. The Lagrangians for a linear resistive transmission line and a linear warm collisional plasma are derived as examples. Using energy considerations, the correct macroscopic plasma Lagrangian is shown to differ from the velocity-integrated low Lagrangian by a macroscopic potential energy that equals twice the particle thermal kinetic energy plus the energy lost by heat conduction.

Peng, Y. K. M.↗

Chemoviscosity modeling for thermosetting resins - I

A new analytical model for chemoviscosity variation during cure of thermosetting resins was developed. This model is derived by modifying the widely used WLF (Williams-Landel-Ferry) Theory in polymer rheology. Major assumptions involved are that the rate of reaction is diffusion controlled and is linearly inversely proportional to the viscosity of the medium over the entire cure cycle. The resultant first order nonlinear differential equation is solved numerically, and the model predictions compare favorably with experimental data of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. The model describes chemoviscosity up to a range of six orders of magnitude under isothermal curing conditions. The extremely non-linear chemoviscosity profile for a dynamic heating cure cycle is predicted as well. The model is also shown to predict changes of glass transition temperature for the thermosetting resin during cure. The physical significance of this prediction is unclear at the present time, however, and further research is required. From the chemoviscosity simulation point of view, the technique of establishing an analytical model as described here is easily applied to any thermosetting resin. The model thus obtained is used in real-time process controls for fabricating composite materials.

Hou, T. H.↗

Parameter identification for nonlinear aerodynamic systems

Parameter identification for nonlinear aerodynamic systems is examined. It is presumed that the underlying model can be arranged into an input/output (I/O) differential operator equation of a generic form. The algorithm estimation is especially efficient since the equation error can be integrated exactly given any I/O pair to obtain an algebraic function of the parameters. The algorithm for parameter identification was extended to the order determination problem for linear differential system. The degeneracy in a least squares estimate caused by feedback was addressed. A method of frequency analysis for determining the transfer function G(j omega) from transient I/O data was formulated using complex valued Fourier based modulating functions in contrast with the trigonometric modulating functions for the parameter estimation problem. A simulation result of applying the algorithm is given under noise-free conditions for a system with a low pass transfer function.

Pearson, Allan E.↗

Satellite tracking by combined optimal estimation and control techniques.

Combined optimal estimation and control techniques are applied for the first time to satellite tracking systems. Both radio antenna and optical tracking systems of NASA are considered. The optimal estimation is accomplished using an extended Kalman filter resulting in an estimated state of the satellite and of the tracking system. This estimated state constitutes an input to the optimal controller. The optimal controller treats a linearized system with a quadratic performance index. The maximum principle is applied and a steady-state approximation to the resulting Riccati equation is obtained. A computer program, RATS, implementing this algorithm is described. A feasibility study of real-time implementation, tracking simulations, and parameter sensitivity studies are also reported.

Dressler, R. M.↗

Analysis of the Rim Inertial Measuring System (RIMS)

Nonlinear equations of motion are derived for the rim inertial measuring system, which is mounted in a strapped down configuration on a carrier vehicle. The measurement and equations derived from angular rates and linear accelerations were interfaced with the dynamic model of a user defined carried vehicle. Two methods are presented for rate measurement, and the results of the nonlinear computer simulation are presented.

Joshi, S. M.↗

Evaluation of automated decisionmaking methodologies and development of an integrated robotic system simulation. Volume 1: Study results

A variety of artificial intelligence techniques which could be used with regard to NASA space applications and robotics were evaluated. The techniques studied were decision tree manipulators, problem solvers, rule based systems, logic programming languages, representation language languages, and expert systems. The overall structure of a robotic simulation tool was defined and a framework for that tool developed. Nonlinear and linearized dynamics equations were formulated for n link manipulator configurations. A framework for the robotic simulation was established which uses validated manipulator component models connected according to a user defined configuration.

Lowrie, J. W.↗

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

True covariance simulation of the EUVE update filter

A covariance analysis of the performance and sensitivity of the attitude determination Extended Kalman Filter (EKF) used by the On Board Computer (OBC) of the Extreme Ultra Violet Explorer (EUVE) spacecraft is presented. The linearized dynamics and measurement equations of the error states are derived which constitute the truth model describing the real behavior of the systems involved. The design model used by the OBC EKF is then obtained by reducing the order of the truth model. The covariance matrix of the EKF which uses the reduced order model is not the correct covariance of the EKF estimation error. A true covariance analysis has to be carried out in order to evaluate the correct accuracy of the OBC generated estimates. The results of such analysis are presented which indicate both the performance and the sensitivity of the OBC EKF.

Bar-Itzhack, Itzhack Y.↗