Search NASA⌕ Search

SEARCH · Search NASA

Results for “System of 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗
NASA NTRS2024

Derivation of Equations of Motion and Model-based Control for 3D Rigid Body Approximation of LUVOIR

The Large UV/Optical/IR Surveyor (LUVOIR) is a conceptual space observatory under development by NASA Goddard Space Flight Center for the Astronomy and Astrophysics Decadal Survey. As with similar space telescopes, e.g., James Webb, this spacecraft will have ultra-stringent pointing stability requirements, i.e., on the order of miliarcseconds. Thus, the design of the LUVOIR attitude control system (ACS) may benefit from a model-based approach. In this document, we present an approximate model of the LUVOIR composed of two 3D rigid bodies connected by a 1-DOF rotary joint. From this model, we derive the nonlinear equations of motion, following the approach of Eric T. Stoneking [1], and linearize them about an arbitrary reference orientation. Using the subsequent state space representation, it is straightforward to compute a set of control gains (e.g., via the ”lqr” or ”place” commands in MATLAB). We conclude this memo with a comparison of our LQR approach with the existing PID slew controller designed by a former NASA intern.

William Bentz↗
NASA NTRS1989

Application of direct solvers to unstructured meshes for the Euler and Navier-Stokes equations using upwind schemes

The application of Newton iteration to inviscid and viscous airfoil calculations on unstructured meshes is examined. A cell-centered finite volume scheme is employed on an unstructured mesh consisting of triangles. Roe's flux difference splitting scheme is used to compute the inviscid fluxes. Higher order accuracy is achieved by an interpolation procedure that makes use of auxiliary gradients. The efficient solution of the sparse linear system of equations which arises upon linearization in time is addressed. Results are presented for inviscid and viscous test cases. The complications which arise due to the introduction of nonlinear limiters are addressed.

Venkatakrishnan, V.↗
NASA NTRS1989

Efficient loads analyses of Shuttle-payloads using dynamic models with linear or nonlinear interfaces

An efficient method for the load analysis of Shuttle-payload systems with linear or nonlinear attachment interfaces is presented which allows the kinematics of the interface degrees of freedom at a given time to be evaluated without calculating the combined system modal representation of the Space Shuttle and its payload. For the case of a nonlinear dynamic model, an iterative procedure is employed to converge the nonlinear terms of the equations of motion to reliable values. Results are presented for a Shuttle abort landing event.

Spanos, P. D.↗
NASA NTRS1990

A control problem for Burgers' equation with bounded input/output

A stabilization problem for Burgers' equation is considered. Using linearization, various controllers are constructed which minimize certain weighted energy functionals. These controllers produce the desired degree of stability for the closed-loop nonlinear system. A numerical scheme for computing the feedback gain functional is developed and several numerical experiments are performed to show the theoretical results.

Burns, John A.↗
NASA NTRS1990

NASA-UCLA Workshop on Computational Techniques in Identification and Control of Flexible Flight Structures, Lake Arrowhead, CA, Nov. 2-4, 1989, Proceedings

The present conference on the identification and control of flexible flight structures with computational techniques encompasses existing and planned testbeds, modeling techniques, adaptive control, and numerical computation and simulation. Specific issues addressed include the NASA/MFSC ground test facility, an active control test on the NASA Minimast, the Astrex testbed for large precision space structures, the dynamic analysis of truss-beam modeling, combined structure-controls-integrated optimization, distributed parameter modeling applicable to the control of flexible flight structures, and approaches to nonclassically damped dynamic systems. Also addressed are adaptive control concepts, techniques for vibration and decentralized slew maneuver control, the discretization of second-order-in-time differential equations which are applicable to nonlinear wave problems, the finite-dimensional approximations of Riccati equations, and numerical simulations of actively controlled space structures.

Balakrishnan, A. V.↗
NASA NTRS2000

Nonlinear Optical Properties and Applications of Polydiacetylene

Recently, we have demonstrated a picosecond all-optical switch, which also functions as a partial all-optical NAND logic gate using a novel polydiacetylene that is synthesized in our laboratory. The nonlinear optical properties of the polydiacetylene material are measured using the Z-scan technique. A theoretical model based on a three level system is investigated and the rate equations of the system are solved. The theoretical calculations are proven to match nicely with the experimental results. The absorption cross-sections for both the first and higher excited states are estimated. The analyses also show that the material suffers a photochemical change beyond a certain level of the laser power and its physical properties suffer radical changes. These changes are the cause for the partial NAND gate function and the switching mechanism.

Abdeldayem, Hossin↗
NASA NTRS2005

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗
NASA NTRS2021

A Functional Interpolation Approach to Compute Period Orbits in the Circular Restricted Three-body Problem

In this paper, we develop a method to solve for periodic orbits, i.e., Lyapunov and Halo orbits, using a functional interpolation scheme called the Theory of Func- tional Connections (TFC). Using this technique, a periodic constraint is analyti- cally embedded into the TFC constrained expression. By doing this, the system of differential equations governing the three-body problem is transformed into an unconstrained optimization problem where simple numerical schemes can be used to find a solution, e.g., nonlinear least-squares is used. This allows for a simpler numerical implementation with comparable accuracy and speed to the traditional differential corrector method.

Mortari, Daniele↗
NASA NTRS2007

Numerical Technique for Analyzing Rotating Rake Mode Measurements in a Duct With Passive Treatment and Shear Flow

A technique is presented for the analysis of measured data obtained from a rotating microphone rake system. The system is designed to measure the interaction modes of ducted fans. A Fourier analysis of the data from the rotating system results in a set of circumferential mode levels at each radial location of a microphone inside the duct. Radial basis functions are then least-squares fit to this data to obtain the radial mode amplitudes. For ducts with soft walls and mean flow, the radial basis functions must be numerically computed. The linear companion matrix method is used to obtain both the eigenvalues of interest, without an initial guess, and the radial basis functions. The governing equations allow for the mean flow to have a boundary layer at the wall. In addition, a nonlinear least-squares method is used to adjust the wall impedance to best fit the data in an attempt to use the rotating system as an in-duct wall impedance measurement tool. Simulated and measured data are used to show the effects of wall impedance and mean flow on the computed results.

Dahl, Milo D.↗
NASA NTRS1981

Iterative optical vector-matrix processors (survey of selected achievable operations)

An iterative optical vector-matrix multiplier with a microprocessor-controlled feedback loop capable of performing a wealth of diverse operations was described. A survey and description of many of its operations demonstrates the versatility and flexibility of this class of optical processor and its use in diverse applications. General operations described include: linear difference and differential equations, linear algebraic equations, matrix equations, matrix inversion, nonlinear matrix equations, deconvolution and eigenvalue and eigenvector computations. Engineering applications being addressed for these different operations and for the IOP are: adaptive phased-array radar, time-dependent system modeling, deconvolution and optimal control.

Casasent, D.↗
NASA NTRS1982

Identification of aerodynamic indicial functions using flight data

It is pointed out that the use of indicial function representation provides a model superior to the aerodynamic derivative model. Specific derivatives can be approximated from the indicial models. The model can also be used to compute equivalent stability and control parameters not usually available from flight data. It is shown that derivatives regarding the angle-of-attack and the side slip angle can be derived directly from the indicial functions without any identifiability problem. Attention is given to the pitch moment coefficient, linear indicial function representation, the identification problem for the pitch moment equation, the identifiability of linear systems, parametric representations of the indicial functions, an identification technique, angle-of-attack and pitch rate dynamics in the pitch plane, multivariate linear models, nonlinear aerodynamic indicial functions, measurement system accuracy, and poststall and spin-entry data from a scaled research vehicle.

Gupta, N. K.↗
NASA NTRS1985

Elastohydrodynamic lubrication calculations used as a tool to study scuffing

A new Reynolds equation was developed that takes into account the nonlinear viscous behavior of the fluid. The new Reynolds equation considers the nonlinear viscous fluid model of Eyring, the equilibrium equation, the constant mass flow, and the kinematic boundary condition. The new Reynolds equation and the elasticity equation are solved simultaneously by using a system approach and a Newton-Raphson technique. Comparisons are made with results obtained from the classical Reynolds equation. The effects of sliding speed and introducing a bump or a groove within the conjunction are studied. Results are shown for both moderate and heavy loads.

Houpert, L. G.↗
NASA NTRS1987

The numerical simulation of subsonic flutter

The present paper describes a numerical simulation of unsteady, subsonic aeroelastic responses. The technique accounts for aerodynamic nonlinearities associated with angles of attack, vortex-dominated flow, static deformations, and unsteady behavior. The fluid and the wing together are treated as a single dynamic system, and the equations of motion for the structure and flowfield are integrated simultaneously and interactively in the time domain. The method employs an iterative scheme based on a predictor-corrector technique. The aerodynamic loads are computed by the general unsteady vortex-lattice method and are determined simultaneously with the motion of the wing. Two models are used to demonstrate the technique: a rigid wing on an elastic support experiencing plunge and pitch about the elastic axis, and a continuous wing rigidly supported at the root chord experiencing spanwise bending and twisting. The time domain solution coupled with the unsteady vortex-lattice method provides the capability of graphically depicting wing and wake motion. Several graphs that illustrate the time domain behavior of the wing and wake are presented.

Strganac, Thomas W.↗
NASA NTRS1993

Minimum-fuel trajectories for hypersonic vehicles with aeropropulsive interactions

The class of hypersonic vehicle configurations with single stage-to-orbit (SSTO) capability reflect highly integrated airframe and propulsion systems. These designs are also known to exhibit a large degree of interaction between the airframe and engine dynamics. Consequently, even simplified hypersonic models are characterized by tightly coupled nonlinear equations of motion. It is to be determined how these features, along with the stringent mission requirements placed on the vehicle, affect its trajectory performance. In this paper, a trajectory optimization problem for a generic hypersonic lifting body will be formulated, a solution method outlined, and results presented. The solution sought is a minimum-fuel trajectorv for orbit injection. It is also desired to study the effect of heating constraints on the vehicle.

Lovell, T. A.↗
NASA NTRS1997

Application of Least-Squares Adjustment Technique to Geometric Camera Calibration and Photogrammetric Flow Visualization

Flow visualization produces data in the form of two-dimensional images. If the optical components of a camera system are perfect, the transformation equations between the two-dimensional image and the three-dimensional object space are linear and easy to solve. However, real camera lenses introduce nonlinear distortions that affect the accuracy of transformation unless proper corrections are applied. An iterative least-squares adjustment algorithm is developed to solve the nonlinear transformation equations incorporated with distortion corrections. Experimental applications demonstrate that a relative precision on the order of 40,000 is achievable without tedious laboratory calibrations of the camera.

Chen, Fang-Jenq↗
NASA NTRS1991

Predicting chaos for infinite dimensional dynamical systems: The Kuramoto-Sivashinsky equation, a case study

The results of extensive computations are presented in order to accurately characterize transitions to chaos for the Kuramoto-Sivashinsky equation. In particular, the oscillatory dynamics in a window that supports a complete sequence of period doubling bifurcations preceding chaos is followed. As many as thirteen period doublings are followed and used to compute the Feigenbaum number for the cascade and so enable, for the first time, an accurate numerical evaluation of the theory of universal behavior of nonlinear systems, for an infinite dimensional dynamical system. Furthermore, the dynamics at the threshold of chaos exhibit a fractal behavior which is demonstrated and used to compute a universal scaling factor that enables the self-similar continuation of the solution into a chaotic regime.

Smyrlis, Yiorgos S.↗