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 937 records · Page 52

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

The global wind system in the thermosphere.

The upper atmosphere wind distribution is calculated by solving the Navier-Stokes equations with models of upper atmospheric densities and ion distribution. All non-linear terms and viscous terms are included. It is found that the velocities from the non-linear solutions are smaller than those obtained from a linear treatment and that in the equatorial region the magnitude of the average meridional velocity is strongly dependent upon the latitudinal coupling. No appreciable super-rotation is obtained by the inclusion of the nonlinear terms; the small value of the average zonal velocity is mainly due to the ion distribution.

Blum, P. W.↗

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗

The use of a movable telescoping end mass system for the time-optimal control of spinning spacecraft

The time-optimal control of a spin-stabilized spacecraft with a movable telescoping appendage (boom) is considered analytically and numerically. The motion of a control mass at the end of the boom is determined such that the terminal time will be minimized for two-axis control of a symmetric spacecraft. The equations of rotational motion are linearized about the desired state of spin about the symmetry axis. The equations for the transverse angular velocity components have the form of a coupled two dimensional harmonic oscillator with boom motion as a control force. The control function which brings the system to the desired state is known to be a series of positive and negative pulses. If the initial state is such that the system can be driven to rest in a single switch, the responses, switching and final times, and required boom motion may be determined analytically. Some typical numerical results based on these solutions are discussed.

Bainum, P. M.↗

Robustness of linear quadratic state feedback designs in the presence of system uncertainty

The paper deals with the problem of expressing the robustness (stability) property of a linear quadratic state feedback (LQSF) design quantitatively in terms of bounds on the perturbations (modeling errors or parameter variations) in the system matrices so that the closed-loop system remains stable. Nonlinear time-varying and linear time-invariant perturbations are considered. The only computation required in obtaining a measure of the robustness of an LQSF design is to determine the eigenvalues of two symmetric matrices determined when solving the algebraic Riccati equation corresponding to the LQSF design problem. Results are applied to a complex dynamic system consisting of the flare control of a STOL aircraft. The design of the flare control is formulated as an LQSF tracking problem.

Patel, R. V.↗

Stability Limits of a PD Controller for a Flywheel Supported on Rigid Rotor and Magnetic Bearings

Active magnetic bearings are used to provide a long-life, low-loss suspension of a high-speed flywheel rotor. This paper describes a modeling effort used to understand the stability boundaries of the PD controller used to control the active magnetic bearings on a high speed test rig. Limits of stability are described in terms of allowable stiffness and damping values which result in stable levitation of the nonrotating rig. Small signal stability limits for the system is defined as a nongrowth in vibration amplitude of a small disturbance. A simple mass-force model was analyzed. The force resulting from the magnetic bearing was linearized to include negative displacement stiffness and a current stiffness. The current stiffness was then used in a PD controller. The phase lag of the control loop was modeled by a simple time delay. The stability limits and the associated vibration frequencies were measured and compared to the theoretical values. The results show a region on stiffness versus damping plot that have the same qualitative tendencies as experimental measurements. The resulting stability model was then extended to a flywheel system. The rotor dynamics of the flywheel was modeled using a rigid rotor supported on magnetic bearings. The equations of motion were written for the center of mass and a small angle linearization of the rotations about the center of mass. The stability limits and the associated vibration frequencies were found as a function of nondimensional magnetic bearing stiffness and damping and nondimensional parameters of flywheel speed and time delay.

Kascak, Albert F.↗

Legendre-tau approximation for functional differential equations. Part 2: The linear quadratic optimal control problem

The numerical scheme based on the Legendre-tau approximation is proposed to approximate the feedback solution to the linear quadratic optimal control problem for hereditary differential systems. The convergence property is established using Trotter ideas. The method yields very good approximations at low orders and provides an approximation technique for computing closed-loop eigenvalues of the feedback system. A comparison with existing methods (based on averaging and spline approximations) is made.

Ito, K.↗

A study of the control problem of the shoot side environment delivery system of a closed crop growth research chamber

The details of our initial study of the control problem of the crop shoot environment of a hypothetical closed crop growth research chamber (CGRC) are presented in this report. The configuration of the CGRC is hypothetical because neither a physical subject nor a design existed at the time the study began, a circumstance which is typical of large scale systems control studies. The basis of the control study is a mathematical model which was judged to adequately mimic the relevant dynamics of the system components considered necessary to provide acceptable realism in the representation. Control of pressure, temperature, and flow rate of the crop shoot environment, along with its oxygen, carbon dioxide, and water concentration is addressed. To account for mass exchange, the group of plants is represented in the model by a source of oxygen, a source of water vapor, and a sink for carbon dioxide. In terms of the thermal energy exchange, the group of plants is represented by a surface with an appropriate temperature. Most of the primitive equations about an experimental operating condition and a state variable representation which was extracted from the linearized equations are presented. Next, we present the results of a real Jordan decomposition and the repositioning of an undesirable eigenvalue via full state feedback. The state variable representation of the modeling system is of the nineteenth order and reflects the eleven control variables and eight system disturbances. Five real eigenvalues are very near zero, with one at zero, three having small magnitude positive values, and one having a small magnitude negative value. A Singular Value Decomposition analysis indicates that these non-zero eigenvalues are not results of numerical error.

Blackwell, C. C.↗

The identification of a distributed parameter model for a flexible structure

A computational method is developed for the estimation of parameters in a distributed model for a flexible structure. The structure we consider (part of the RPL experiment) consists of a cantilevered beam with a thruster and linear accelerometer at the free end. The thruster is fed by a pressurized hose whose horizontal motion effects the transverse vibration of the beam. The Euler-Bernoulli theory is used to model the vibration of the beam and treat the hose-thruster assembly as a lumped or point mass-dashpot-spring system at the tip. Using measurements of linear acceleration at the tip, it is estimated that the parameters (mass, stiffness, damping) and a Voight-Kelvin viscoelastic structural damping parameter for the beam using a least squares fit to the data. Spline based approximations to the hybrid (coupled ordinary and partial differential equations) system are considered; theoretical convergence results and numerical studies with both simulation and actual experimental data obtained from the structure are presented and discussed.

Banks, H. T.↗

Frequency analysis via the method of moment functionals

Several variants are presented of a linear-in-parameters least squares formulation for determining the transfer function of a stable linear system at specified frequencies given a finite set of Fourier series coefficients calculated from transient nonstationary input-output data. The basis of the technique is Shinbrot's classical method of moment functionals using complex Fourier based modulating functions to convert a differential equation model on a finite time interval into an algebraic equation which depends linearly on frequency-related parameters.

Pearson, A. E.↗

Frequency analysis via the method of moment functionals

Several variants are presented of a linear-in-parameters least squares formulation for determining the transfer function of a stable linear system at specified frequencies given a finite set of Fourier series coefficients calculated from transient nonstationary input-output data. The basis of the technique is Shinbrot's classical method of moment functionals using complex Fourier based modulating functions to convert a differential equation model on a finite time interval into an algebraic equation which depends linearly on frequency-related parameters.

Pearson, A. E.↗

Effect of a crystal-melt interface on Taylor-vortex flow

The linear stability of circular Couette flow between concentric infinite cylinders is considered for the case that the stationary outer cylinder is a crystal-melt interface rather than a rigid surface. A radial temperature difference is maintained across the liquid gap, and equations for heat transport in the crystal and melt phases are included to extend the ordinary formulation of this problem. The stability of this two-phase system depends on the Prandtl number. For small Prandtl number the linear stability of the two-phase system is given by the classical results for a rigid-walled system. For increasing values of the Prandtl number, convective heat transport becomes significant and the system becomes increasingly less stable. Previous results in a narrow-gap approximation are extended to the case of a finite gap, and both axisymmetric and nonaxisymmetric disturbance modes are considered. The two-phase system becomes less stable as the finite gap tends to the narrow-gap limit. The two-phase system is more stable to nonaxisymmetric modes with azimuthal wavenumber n = 1; the stability of these n = 1 modes is sensitive to the latent heat of fusion.

Mcfadden, G. B.↗

Design of minimax output feedback controller for system with parameter uncertainty

The problem of controlling a time-invariant system with parameter uncertainty is considered with incomplete state feedback. The controller is designed by minimaximizing a quadratic performance criterion and a sensitivity (or loss) criterion, involving the state of the system, the control, and the uncertainty vector. The resulting optimal controller is linear and optimal feedback gain matrix must satisfy a set of nonlinear algebraic equations. some algorithms for algebraic minimax problems are presented.

Basuthakur, S.↗

A new approach to the design of time varying control systems with application to the space shuttle booster

A general approach toward the analysis and design of closed loop control for slowly time-varying linear systems was investigated. The approach is by the method of generalized multiple scales, in which slow and fast dynamics are systematically separated by employing different clocks which measure time at varying rates. The clocks, which are necessarily nonlinear functions of time, so that the system dynamics are asymptotically invariant with respect to the new time scale. A transfer function relating the output to the input for general linear slowly time-varying systems is developed and represents the actual system under certain conditions. The clock functions are shown to satisfy an algebraic characteristic equation and can be determined in general in terms of the coefficients.

Ramnath, R. V.↗

Direct numerical solution of the transonic perturbation integral equation for lifting and nonlifting airfoils

The linear transonic perturbation integral equation previously derived for nonlifting airfoils is formulated for lifting cases. In order to treat shock wave motions, a strained coordinate system is used in which the shock location is invariant. The tangency boundary conditions are either formulated using the thin airfoil approximation or by using the analytic continuation concept. A direct numerical solution to this equation is derived in contrast to the iterative scheme initially used, and results of both lifting and nonlifting examples indicate that the method is satisfactory.

Nixon, D.↗

Control of a large flexible platform in orbit

The dynamics and attitude and shape control of a large thin flexible platform in orbit are studied. Attitude and shape control is assumed to result from actuators placed perpendicular to the main surface and one edge and their effect on the rigid body and elastic modes is modelled to first order. The equations of motion are linearized about nominal orientations where the undeformed plate follows either the local vertical or local horizontal. The stability of the uncontrolled system is investigated analytically. Once controllability is established for a set of actuator locations, control law development is based on pole placement, decoupling, and linear optimal control theory.

Reddy, A. S. S. R.↗

Singular perturbation equations for flexible satellites

Force equations of motion of the individual flexible elements of a satellite were obtained in a previous paper. Moment equations of motion of the composite bodies of a flexible satellite are to be developed using two sets of equations which form the basic system for any dynamic model of flexible satellites. This basic system consists of a set of N-coupled, nonlinear, ordinary, or partial differential equations, for a flexible satellite with n generalized, structural position coordinates. For single composite body satellites, N is equal to (n + 3); for dual-spin systems, N is equal to (n + 9). These equations involve time derivatives up to the second order. The study shows a method of avoiding this linearization by reducing the N equations to 3 or 9 nonlinear, coupled, first order, ordinary, differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can be solved by known methods.

Huang, T. C.↗

Midlatitude tropical interactions during winter

Pre-FGGE and FGGE/MONEX data are used to identify short term midlatitude tropical and longitudinal interactions during the winter monsoon. These interactions occur as cold surges, which develop over the East Asian continent and penetrate deep into the tropics with fast gravity wave speed. The observed interactions that occur after a surge include cyclogenesis and enhanced convection in the equatorial region, feedback from equatorial convection to midlatitude circulation systems, tropical east-west (Walker) circulations, and cross-equatorial influence. These interactions are also studied theoretically by analytical solutions of linearized shallow water equations. Response to transient forcing (monsoon surges) are mainly in Rossby and Kelvin modes. When the forcing time scale is short, significant gravity modes are also excited. The responses closely resemble observed winter monsoon flow. Responses to stationary forcing show that deep (barotropic) motions propagate energy away into high latitudes and that shallow (baroclinic) motions are trapped around the equator. It is shown that the barotropic teleconnection-type response to tropical sources found in previous numerical studies was due to the specified vertical wind shear and surface friction.

Chang, C. P.↗