Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear equations”

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 343 records · Page 19

Aeroelastic Analysis of a Trimmed Generic Hypersonic Vehicle

The aeroelastic equations of motion governing a hypersonic vehicle in free flight are derived. The equations of motion for a translating and rotating flexible body using Lagrange's equations in terms of quasi-coordinates are presented. These equations are simplified for the case of a vehicle with pitch and plunge rigid body degrees of freedom and small elastic displacements. The displacements are approximated by a truncated series of the unrestrained mode shapes, which are obtained using equivalent plate theory. Subsequently, the nonlinear equations of motion are linearized about the trim state, which is obtained using a rigid body trim model and steady hypersonic aerodynamics. The appropriate flutter derivatives are calculated from piston theory. Results describing mode shapes, trim behavior, and aeroelastic stability of a generic hypersonic vehicle are presented.

Nydick, I.↗

Reentry Vehicle Flight Controls Design Guidelines: Dynamic Inversion

This report addresses issues in developing a flight control design for vehicles operating across a broad flight regime and with highly nonlinear physical descriptions of motion. Specifically it addresses the need for reentry vehicles that could operate through reentry from space to controlled touchdown on Earth. The latter part of controlled descent is achieved by parachute or paraglider - or by all automatic or a human-controlled landing similar to that of the Orbiter. Since this report addresses the specific needs of human-carrying (not necessarily piloted) reentry vehicles, it deals with highly nonlinear equations of motion, and then-generated control systems must be robust across a very wide range of physics. Thus, this report deals almost exclusively with some form of dynamic inversion (DI). Two vital aspects of control theory - noninteracting control laws and the transformation of nonlinear systems into equivalent linear systems - are embodied in DI. Though there is no doubt that the mathematical tools and underlying theory are widely available, there are open issues as to the practicality of using DI as the only or primary design approach for reentry articles. This report provides a set of guidelines that can be used to determine the practical usefulness of the technique.

Ito, Daigoro↗

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗

Pseudosteady-state analysis of nonlinear aircraft maneuvers

An analytical method was developed for studying the combined effects of rotational coupling and nonlinear aerodynamics on aircraft response for specified control inputs. The method involves the simultaneous solution of two nonlinear equations which are functions of angle attack, roll rate, and control inputs. The method was applied to a number of maneuvers for a fighter-type aircraft. Time history responses verified the usefulness of the analysis for predicting a variety of response characteristics caused by interacting nonlinear aerodynamic and inertial effects, including spin conditions.

Young, J. W.↗

Optimal feedback control of a bioreactor with a remote sensor

Sensors used to monitor bioreactor conditions directly often perform poorly in the face of adverse nonphysiological conditions. One way to circumvent this is to use a remote sensor block. However, such a configuration usually causes a significant time lag between measurements and the actual state values. Here, the problem of implementing feedback control strategies for such systems, described by nonlinear equations, is addressed. The problem is posed as an optimal control problem with a linear quadratic performance index. The linear control law so obtained is used to implement feedback. A global linearization technique as well as an expansion using Taylor series is used to linearize the nonlinear system, and the feedback is subsequently implemented.

Niranjan, S. C.↗

Dynamics of spinning spacecraft during the deployment of flexible boom

The effect of flexibility on the dynamics of a spinning symmetrical spacecraft during the extension (or retraction) of boom-type appendages is treated by modeling the appendages as compound spherical pendula with varying lengths. Both boom stiffness and structural damping are included. For constant boom length the resulting linearized equations of motion contain periodic coefficients involving the spin frequency. A bounded transformation converts this system into a kinematically equivalent one (with only constant coefficients), whose stability is analyzed using the Kelvin-Tait-Chetaev theorem. The dynamics during extension is predicted by numerical simulation of the nonlinear equations.

James, P. K.↗

Pseudosteady state analysis of nonlinear aircraft maneuvers

An analytical method has been developed for studying the combined effects of rotational coupling and nonlinear aerodynamics on aircraft response for specified control inputs. The method involves the simultaneous solution of two nonlinear equations which are functions of angle of attack, roll rate, and the control inputs. The method is applied to a number of maneuvers for a fighter type aircraft. Time history responses verified the usefulness of the analysis for predicting a variety of response characteristics caused by interacting nonlinear aerodynamic and inertial effects, including spin conditions.

Young, J. W.↗

Nonlinear evolution of modes in a compressible boundary layer

The nonlinear stability of an oblique mode propagating in a 2D compressible boundary layer is considered under the long wave-length approximation. The growth rate of the wave is assumed to be small so that the ideas of nonlinear critical layers can be applied. It is shown that the spatial/temporal evolution of the mode is governed by a pair of coupled unsteady nonlinear equations for the disturbance vorticity and density. Expressions for the linear growth rate shows clearly the effects of wall heating and cooling, and in particular how heating destabilizes the boundary layer for these long wavelength inviscid modes at 0(1) Mach numbers. The numerical solution of the nonlinear problem is discussed and some results are presented.

Gajjar, J. S. B.↗

Implicit Multirate GARK Methods

This article considers multirate generalized-structure additively partitioned Runge–Kutta methods for solving stiff systems of ordinary differential equations with multiple time scales. These methods treat different partitions of the system with different timesteps for a more targeted and efficient solution compared to monolithic single rate approaches. With implicit methods used across all partitions, methods must find a balance between stability and the cost of solving nonlinear equations for the stages. In order to characterize this important trade-off, we explore multirate coupling strategies, problems for assessing linear stability, and techniques to efficiently implement Newton iterations for stage equations. Unlike much of the existing multirate stability analysis which is limited in scope to particular methods, we present general statements on stability and describe fundamental limitations for certain types of multirate schemes. New implicit multirate methods up to fourth order are derived, and their accuracy and efficiency properties are verified with numerical tests.

97 MATHEMATICS AND COMPUTING↗

The dynamics and control of large flexible space structures - 12, supplement 11

The rapid 2-D slewing and vibrational control of the unsymmetrical flexible SCOLE (Spacecraft Control Laboratory Experiment) with multi-bounded controls is considered. Pontryagin's Maximum Principle is applied to the nonlinear equations of the system to derive the necessary conditions for the optimal control. The resulting two point boundary value problem is then solved by using the quasilinearization technique, and the near minimum time is obtained by sequentially shortening the slewing time until the controls are near the bang-bang type. The tradeoff between the minimum time and the minimum flexible amplitude requirements is discussed. The numerical results show that the responses of the nonlinear system are significantly different from those of the linearized system for rapid slewing. The SCOLE station-keeping closed loop dynamics are re-examined by employing a slightly different method for developing the equations of motion in which higher order terms in the expressions for the mast modal shape functions are now included. A preliminary study on the effect of actuator mass on the closed loop dynamics of large space systems is conducted. A numerical example based on a coupled two-mass two-spring system illustrates the effect of changes caused in the mass and stiffness matrices on the closed loop system eigenvalues. In certain cases the need for redesigning control laws previously synthesized, but not accounting for actuator masses, is indicated.

Bainum, Peter M.↗

Equations of Motion for the g-LIMIT Microgravity Vibration Isolation System

A desirable microgravity environment for experimental science payloads may require an active vibration isolation control system. A vibration isolation system named g-LIMIT (GLovebox Integrated Microgravity Isolation Technology) is being developed by NASA Marshall Space Flight Center to support microgravity science experiments using the microgravity science glovebox. In this technical memorandum, the full six-degree-of-freedom nonlinear equations of motion for g-LIMIT are derived. Although the motivation for this model development is control design and analysis of g-LIMIT, the equations are derived for a general configuration and may be used for other isolation systems as well.

Kim, Y. K.↗

Nonlinear dynamic analysis of deploying flexible space booms

The dynamics of the planar deployment and lock-up of two flexible boom-type appendages on a spacecraft are analyzed by establishing nonlinear equations of motion for the deployment phase and linear for the post-lock-up phase. Nondimensional parameters and threshold values are identified in terms of boom tip deflections, deployment time and root moments for appendages with tip masses attached to a central rigid body through a rotational spring; an account is taken of the nonlinear kinematic and structural terms. It was found that, in general, the solution is influenced only by two mass ratios and one nondimensional stiffness parameter.

Mcgowan, P. E.↗

A Computing Method for Sound Propagation Through a Nonuniform Jet Stream

Understanding the principles of jet noise propagation is an essential ingredient of systematic noise reduction research. High speed computer methods offer a unique potential for dealing with complex real life physical systems whereas analytical solutions are restricted to sophisticated idealized models. The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions and a more suitable approach was needed. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Padula, S. L.↗

A computing method for sound propagation through a nonuniform jet stream

The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Padula, S. L.↗

Dynamic analysis and control of lightweight manipulators with flexible parallel link mechanisms

The flexible parallel link mechanism is designed for increased rigidity to sustain the buckling when it carries a heavy payload. Compared to a one link flexible manipulator, a two link flexible manipulator, especially the flexible parallel mechanism, has more complicated characteristics in dynamics and control. The objective of this research is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model. The step response of the analytical model and the TREETOPS model match each other well. The nonlinear dynamics is studied using a sinusoidal excitation. The actuator dynamic effect on a flexible robot was investigated. The effects are explained by the root loci and the Bode plot theoretically and experimentally. For the base performance for the advanced control scheme, a simple decoupled feedback scheme is applied.

Lee, Jeh Won↗

Plasma diffusion at the magnetopause? The case of lower hybrid drift waves

The diffusion expected from the quasilinear theory of the lower hybrid drift instability at the Earth's magnetopause is recalculated. The resulting diffusion coefficient is in principle just marginally large enough to explain the thickness of the boundary layer under quiet conditions, based on observational upper limits for the wave intensities. Thus, one possible model for the boundary layer could involve equilibrium between the diffusion arising from lower hybrid waves and various low processes. However, some recent data and simulations seems to indicate that the magnetopause is not consistent with such a soft diffusive equilibrium model. Furthermore, investigation of the nonlinear equations for the lower hybrid waves for magnetopause parameters indicates that the quasilinear state may never arise because coalescence to large wavelengths, followed by collapse once a critical wavelengths is reached, occur on a time scale faster than the quasilinear diffusion. In this case, an inhomogeneous boundary layer is to be expected. More simulations are required over longer time periods to explore whether this nonlinear evolution really takes place at the magnetopause.

Treumann, R. A.↗

MATLAB Stability and Control Toolbox Trim and Static Stability Module

MATLAB Stability and Control Toolbox (MASCOT) utilizes geometric, aerodynamic, and inertial inputs to calculate air vehicle stability in a variety of critical flight conditions. The code is based on fundamental, non-linear equations of motion and is able to translate results into a qualitative, graphical scale useful to the non-expert. MASCOT was created to provide the conceptual aircraft designer accurate predictions of air vehicle stability and control characteristics. The code takes as input mass property data in the form of an inertia tensor, aerodynamic loading data, and propulsion (i.e. thrust) loading data. Using fundamental nonlinear equations of motion, MASCOT then calculates vehicle trim and static stability data for the desired flight condition(s). Available flight conditions include six horizontal and six landing rotation conditions with varying options for engine out, crosswind, and sideslip, plus three take-off rotation conditions. Results are displayed through a unique graphical interface developed to provide the non-stability and control expert conceptual design engineer a qualitative scale indicating whether the vehicle has acceptable, marginal, or unacceptable static stability characteristics. If desired, the user can also examine the detailed, quantitative results.

Kenny, Sean P.↗

Numerical computations on one-dimensional inverse scattering problems

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

Dunn, M. H.↗