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 1,279 records · Page 71

A variational assimilation method for satellite and conventional data: Development of basic model for diagnosis of cyclone systems

A three-dimensional diagnostic model for the assimilation of satellite and conventional meteorological data is developed with the variational method of undetermined multipliers. Gridded fields of data from different type, quality, location, and measurement source are weighted according to measurement accuracy and merged using least squares criteria so that the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation are satisfied. The model is used to compare multivariate variational objective analyses with and without satellite data with initial analyses and the observations through criteria that were determined by the dynamical constraints, the observations, and pattern recognition. It is also shown that the diagnoses of local tendencies of the horizontal velocity components are in good comparison with the observed patterns and tendencies calculated with unadjusted data. In addition, it is found that the day-night difference in TOVS biases are statistically different (95% confidence) at most levels. Also developed is a hybrid nonlinear sigma vertical coordinate that eliminates hydrostatic truncation error in the middle and upper troposphere and reduces truncation error in the lower troposphere. Finally, it is found that the technique used to grid the initial data causes boundary effects to intrude into the interior of the analysis a distance equal to the average separation between observations.

Achtemeier, G. L.↗

Modeling, Simulation and Control of a Propellant Mixer

This paper addresses some modeling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer must combine high pressure liquid hydrogen (LH2) and gaseous hydrogen (GH2) to produce and output flow that meets certain thermodynamic properties before it is fed into a test article. More precisely, this paper considers that the quantities to be tracked and/or regulated are mixer internal pressure, exit mass flow, and exit temperature. The available control inputs are given by three value positions, namely those of the GH2, LH2 and exit valves. The mixer is modelled by a system of two nonlinear ordinary differential equations having density and internal energy as states. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as controllability and uniqueness of the equilibrium point are shown to hold. Validation of the model against real data is also provided. As a first control approach, a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions. A feedback linearization controller is designed and used to achieve tracking and regulation of the outputs over an extended range of the variables of interest.

Richter, Hanz↗

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↗

A nonlinear analysis of pulsatile flow in arteries.

An approximate numerical method for calculating flow profiles in arteries is developed. The theory takes into account the nonlinear terms of the Navier-Stokes equations as well as the nonlinear behaviour and large deformations of the arterial wall. Through the locally measured values of the pressure, pressure gradient, and pressure-radius function, the velocity distribution and wall shear at a given location along the artery can be determined. The computed results agree well with the corresponding experimental data.

Ling, S. C.↗

Fully Nonlinear Modeling and Analysis of Precision Membranes

High precision membranes are used in many current space applications. This paper presents a fully nonlinear membrane theory with forward and inverse analyses of high precision membrane structures. The fully nonlinear membrane theory is derived from Jaumann strains and stresses, exact coordinate transformations, the concept of local relative displacements, and orthogonal virtual rotations. In this theory, energy and Newtonian formulations are fully correlated, and every structural term can be interpreted in terms of vectors. Fully nonlinear ordinary differential equations (ODES) governing the large static deformations of known axisymmetric membranes under known axisymmetric loading (i.e., forward problems) are presented as first-order ODES, and a method for obtaining numerically exact solutions using the multiple shooting procedure is shown. A method for obtaining the undeformed geometry of any axisymmetric membrane with a known inflated geometry and a known internal pressure (i.e., inverse problems) is also derived. Numerical results from forward analysis are verified using results in the literature, and results from inverse analysis are verified using known exact solutions and solutions from the forward analysis. Results show that the membrane theory and the proposed numerical methods for solving nonlinear forward and inverse membrane problems are accurate.

Pai, P. Frank↗

Nonlinear analysis and performance evaluation of the Annular Suspension and Pointing System (ASPS)

The Annular Suspension and Pointing System (ASPS) can provide high accurate fine pointing for a variety of solar-, stellar-, and Earth-viewing scientific instruments during space shuttle orbital missions. In this report, a detailed nonlinear mathematical model is developed for the ASPS/Space Shuttle system. The equations are augmented with nonlinear models of components such as magnetic actuators and gimbal torquers. Control systems and payload attitude state estimators are designed in order to obtain satisfactory pointing performance, and statistical pointing performance is predicted in the presence of measurement noise and disturbances.

Joshi, S. M.↗

Discrete time learning control in nonlinear systems

In this paper digital learning control methods are developed primarily for use in single-input, single-output nonlinear dynamic systems. Conditions for convergence of the basic form of learning control based on integral control concepts are given, and shown to be satisfied by a large class of nonlinear problems. It is shown that it is not the gross nonlinearities of the differential equations that matter in the convergence, but rather the much smaller nonlinearities that can manifest themselves during the short time interval of one sample time. New algorithms are developed that eliminate restrictions on the size of the learning gain, and on knowledge of the appropriate sign of the learning gain, for convergence to zero error in tracking a feasible desired output trajectory. It is shown that one of the new algorithms can give guaranteed convergence in the presence of actuator saturation constraints, and indicate when the requested trajectory is beyond the actuator capabilities.

Longman, Richard W.↗

Singular perturbation equations for flexible satellites

The dynamic model of a flexible satellite with n generalized structural position coordinates requires the solution of a set of N coupled nonlinear ordinary or partial differential equations. For single composite body satellites, N is equal to (n + 3). For dual-spin systems, N is equal to (n + 9). These equations usually involve time derivatives up to the second order. For large values of n, linearization of the system has so far been the only practicable way of solution. The present study shows a method of avoiding this linearization by reducing the N equations to three or nine 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 then be solved by known methods.

Huang, T. C.↗

The nonlinear interaction of Tollmien-Schlichting waves and Taylor-Goertler vortices in curved channel flows

It is known that a viscous fluid flow with curved streamlines can support both Tollmien-Schlichting and Taylor-Goertler instabilities. In a situation where both modes are possible on the basis of linear theory a nonlinear theory must be used to determine the effect of the interaction of the instabilities. The details of this interaction are of practical importance because of its possible catastrophic effects on mechanisms used for laminar flow control. This interaction is studied in the context of fully developed flows in curved channels. A part form technical differences associated with boundary layer growth the structures of the instabilities in this flow are very similar to those in the practically more important external boundary layer situation. The interaction is shown to have two distinct phases depending on the size of the disturbances. At very low amplitudes two oblique Tollmein-Schlichting waves interact with a Goertler vortex in such a manner that the amplitudes become infinite at a finite time. This type of interaction is described by ordinary differential amplitude equations with quadratic nonlinearities.

Hall, P.↗

A feedback linearization approach to spacecraft control using momentum exchange devices

Recent developments in the area of nonlinear control theory have shown how coordiante changes in the state and input spaces can be used with nonlinear feedback to transform certain nonlinear ordinary differential equations into equivalent linear equations. These feedback linearization techniques are applied to resolve two problems arising in the control of spacecraft equipped with control moment gyroscopes (CMGs). The first application involves the computation of rate commands for the gimbals that rotate the individual gyroscopes to produce commanded torques on the spacecraft. The second application is to the long-term management of stored momentum in the system of control moment gyroscopes using environmental torques acting on the vehicle. An approach to distributing control effort among a group of redundant actuators is described that uses feedback linearization techniques to parameterize sets of controls which influence a specified subsystem in a desired way. The approach is adapted for use in spacecraft control with double-gimballed gyroscopes to produce an algorithm that avoids problematic gimbal configurations by approximating sets of gimbal rates that drive CMG rotors into desirable configurations. The momentum management problem is stated as a trajectory optimization problem with a nonlinear dynamical constraint. Feedback linearization and collocation are used to transform this problem into an unconstrainted nonlinear program. The approach to trajectory optimization is fast and robust. A number of examples are presented showing applications to the proposed NASA space station.

Dzielski, John Edward↗

Nonlinear analysis of rotor-bearing systems using component mode synthesis

The method of component mode synthesis is developed to determine the forced response of nonlinear, multishaft, rotor-bearing systems. The formulation allows for simulation of system response due to blade loss, distributed unbalance, base shock, maneuver loads, and specified fixed frame forces. The motion of each rotating component of the system is described by superposing constraint modes associated with boundary coordinates and constrained precessional modes associated with internal coordinates. The precessional modes are truncated for each component and the reduced component equations are assembled with the nonlinear supports and interconnections to form a set of nonlinear system equations of reduced order. These equations are then numerically integrated to obtain the system response. A computer program, which is presently restricted to single shaft systems, has been written and results are presented for transient system response associated with blade loss dynamics with squeeze film dampers, and with interference rubs.

Nelson, H. D.↗

Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

Beam, R. M.↗

An exact solution of an augmented Burgers equation and amplitude-dependent acoustic propagation speed

Nonlinear sound propagation in the atmosphere is usually modeled using an augmented Burgers equation accounting for a weak nonlinearity and atmospheric absorption. Because the absorption includes the molecular vibrational relaxation, such a Burgers equation is more complex than the regular Burgers equation that only accounts for the thermoviscous dissipation in the absorption. Although an exact solution of the regular Burgers equation has long been derived using the Cole- Hopf transform, an exact solution of the augmented Burgers equation has not been derived previously. Thus, this paper presents an exact solution of the augmented Burgers equation. This novel solution is shown to be equivalent to the solution using the Cole-Hopf transform when the absorption only involves thermoviscous dissipation. It can also be reduced to the known solution of an N-wave when the absorption is ignored. The augmented Burgers equation is an approximation valid for weak nonlinearity. However, this assumption may not be accurate for acoustic signals propagating from the lower atmosphere and which are subsequently refracted downward from the upper atmosphere (e.g., stratosphere and thermosphere) due to the decreasing air density with increasing altitude [Lonzaga, et al., Geophysical Journal International, 200(3), pp.1347-1361]. Consequently, the current paper also discusses the effects of a strong nonlinearity that lead to an amplitude-dependent increase in signal propagation speed. For an impulsive signal such as a sonic boom, these effects cause a dispersion of the signal similar to the observed dispersion of acoustic signals from supersonic Concorde as well as from large explosions.

Joel B Lonzaga↗

Nonlinear cellular motions in Poiseuille channel flow

In a number of nonlinear solutions to the equations of channel flow the velocity is decomposed into a mean part plus a nonlinear disturbance. The idea that nonlinear effects place a limitation on the amplitudes of the disturbance flow is considered. In the reported investigation the disturbance flow is represented by drastically limited Fourier expansions in the downstream coordinate. The resulting equations are solved numerically with high accuracy to obtain a good representation of the cross-stream structure of the solution. The results of the investigation show that indeed the nonlinear terms always limit the amplitude of the disturbance flow in this approximation.

Zahn, J.-P.↗

Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies

The particular problems investigated in the present study arise from nonlinear transient heat conduction. One of two types of nonlinearities considered is related to a material temperature dependence which is frequently needed to accurately model behavior over the range of temperature of engineering interest. The second nonlinearity is introduced by radiation boundary conditions. The finite element equations arising from the solution of nonlinear transient heat conduction problems are formulated. The finite element matrix equations are temporally discretized, and a nonlinear iterative solution algorithm is proposed. Algorithms for solving the linear problem are discussed, taking into account the form of the matrix equations, Gaussian elimination, cost, and iterative techniques. Attention is also given to approximate factorization, implementational aspects, and numerical results.

Winget, J. M.↗

Nonlinear (time domain) and linearized (time and frequency domain) solutions to the compressible Euler equations in conservation law form

Two linearized solvers (time and frequency domain) based on a high resolution numerical scheme are presented. The basic approach is to linearize the flux vector by expressing it as a sum of a mean and a perturbation. This allows the governing equations to be maintained in conservation law form. A key difference between the time and frequency domain computations is that the frequency domain computations require only one grid block irrespective of the interblade phase angle for which the flow is being computed. As a result of this and due to the fact that the governing equations for this case are steady, frequency domain computations are substantially faster than the corresponding time domain computations. The linearized equations are used to compute flows in turbomachinery blade rows (cascades) arising due to blade vibrations. Numerical solutions are compared to linear theory (where available) and to numerical solutions of the nonlinear Euler equations.

Sreenivas, Kidambi↗