Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Non linear dynamics”

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 235 records · Page 13

A seal test facility for the measurement of isotropic and anisotropic linear rotordynamic characteristics

A new seal test facility for measuring high-pressure seal rotor-dynamic characteristics has recently been made operational at Case Western Reserve University (CWRU). This work is being sponsored by the Electric Power Research Institute (EPRI). The fundamental concept embodied in this test apparatus is a double-spool-shaft spindle which permits independent control over the spin speed and the frequency of an adjustable circular vibration orbit for both forward and backward whirl. Also, the static eccentricity between the rotating and non-rotating test seal parts is easily adjustable to desired values. By accurately measuring both dynamic radial displacement and dynamic radial force signals, over a wide range of circular orbit frequency, one is able to solve for the full linear-anisotropic model's 12 coefficients rather than the 6 coefficients of the more restrictive isotropic linear model. Of course, one may also impose the isotropic assumption in reducing test data, thereby providing a valid qualification of which seal configurations are well represented by the isotropic model and which are not. In fact, as argued in reference (1), the requirement for maintaining a symmetric total system mass matrix means that the resulting isotropic model needs 5 coefficients and the anisotropic model needs 11 coefficients.

Adams, M. L.↗

Expected Deployment Dynamics of Proseds

The control law for The Propulsive Small Expendable Deployment System (ProSEDS) deployment is a modification of the control routine that was successfully used in the flight of SEDS-II. Unlike SEDS, the tether of ProSEDS consists of different sections with different mechanical characteristics. A non-linear control trajectory in phase-space (i.e., the reference profile) is fed forward to the controller to guide the satellite, at the tether tip, to the desired final state under nominal conditions and no external perturbations. A linear feedback control is applied by the brake to keep the actual trajectory as close as possible to the reference. The paper also shows the results of simulations of deployment dynamics with and without noise. The control law has thus far been developed and tested on the ground for the original ProSEDS tether configuration of 15 km. A new reference will have to be designed and tested for other tether configurations.

Lorenzini, E. C.↗

On Kinetics Modeling of Vibrational Energy Transfer

Two models of vibrational energy exchange are compared at equilibrium to the elementary vibrational exchange reaction for a binary mixture. The first model, non-linear in the species vibrational energies, was derived by Schwartz, Slawsky, and Herzfeld (SSH) by considering the detailed kinetics of vibrational energy levels. This model recovers the result demanded at equilibrium by the elementary reaction. The second model is more recent, and is gaining use in certain areas of computational fluid dynamics. This model, linear in the species vibrational energies, is shown not to recover the required equilibrium result. Further, this more recent model is inconsistent with its suggested rate constants in that those rate constants were inferred from measurements by using the SSH model to reduce the data. The non-linear versus linear nature of these two models can lead to significant differences in vibrational energy coupling. Use of the contemporary model may lead to significant misconceptions, especially when integrated in computer codes considering multiple energy coupling mechanisms.

Gilmore, John O.↗

Mechanics of fracture - Fundamentals and some recent developments

An overview is presented of the fundamental aspects of and recent developments in fracture mechanics. Reference is made to linear elastic fracture mechanics including the state of stresses and displacements in the vicinity of cracks, effects of crack geometry and orientation on stress intensity factors, energy balance of Griffith, Irwin's stress intensity concept, and linear elastic fracture mechanics testing for fracture toughness. Other aspects of this paper include the non-linear behavior of materials and their influence on fracture mechanics parameters, consideration of viscoelasticity and plasticity, non-linear fracture toughness parameters as C.O.D., R-curve and J-integral, and a non-linear energy method, proposed by Liebowitz. Finite element methods applied to fracture mechanics problems are indicated. Also, consideration has been given to slow crack growth, dynamic effects on K(IC), Sih's criterion for fracture, Lee and Liebowitz's criterion relating crack growth with plastic energy, and applications of fracture mechanics to aircraft design. Suggestions are offered for future research efforts to be undertaken in fracture mechanics.

Liebowitz, H.↗

Dynamics of Permanent-Magnet Biased Active Magnetic Bearings

Active magnetic radial bearings are constructed with a combination of permanent magnets to provide bias forces and electromagnets to generate control forces for the reduction of cost and the operating energy consumption. Ring-shaped permanent magnets with axial magnetization are attached to a shaft and share their magnet stators with the electromagnets. The magnet cores are made of solid iron for simplicity. A simplified magnetic circuit of the combined magnet system is analyzed with linear circuit theory by approximating the characteristics of permanent magnets with a linear relation. A linearized dynamical model of the control force is presented with the first-order approximation of the effects of eddy currents. Frequency responses of the rotor motion to disturbance inputs and the motion for impulsive forces are tested in the non-rotating state. The frequency responses are compared with numerical results. The decay of rotor speed due to magnetic braking is examined. The experimental results and the presented linearized model are similar to those of the all-electromagnetic design.

Fukata, Satoru↗

Fast Optimization for Aircraft Descent and Approach Trajectory

We address problem of on-line scheduling of the aircraft descent and approach trajectory. We formulate a general multiphase optimal control problem for optimization of the descent trajectory and review available methods of its solution. We develop a fast algorithm for solution of this problem using two key components: (i) fast inference of the dynamical and control variables of the descending trajectory from the low dimensional flight profile data and (ii) efficient local search for the resulting reduced dimensionality non-linear optimization problem. We compare the performance of the proposed algorithm with numerical solution obtained using optimal control toolbox General Pseudospectral Optimal Control Software. We present results of the solution of the scheduling problem for aircraft descent using novel fast algorithm and discuss its future applications.

Luchinsky, Dmitry G.↗

Comparing Consider-Covariance Analysis with Sigma-Point Consider Filter and Linear-Theory Consider Filter Formulations

Recent literature in applied estimation theory reflects growing interest in the sigma-point (also called unscented ) formulation for optimal sequential state estimation, often describing performance comparisons with extended Kalman filters as applied to specific dynamical problems [c.f. 1, 2, 3]. Favorable attributes of sigma-point filters are described as including a lower expected error for nonlinear even non-differentiable dynamical systems, and a straightforward formulation not requiring derivation or implementation of any partial derivative Jacobian matrices. These attributes are particularly attractive, e.g. in terms of enabling simplified code architecture and streamlined testing, in the formulation of estimators for nonlinear spaceflight mechanics systems, such as filter software onboard deep-space robotic spacecraft. As presented in [4], the Sigma-Point Consider Filter (SPCF) algorithm extends the sigma-point filter algorithm to the problem of consider covariance analysis. Considering parameters in a dynamical system, while estimating its state, provides an upper bound on the estimated state covariance, which is viewed as a conservative approach to designing estimators for problems of general guidance, navigation and control. This is because, whether a parameter in the system model is observable or not, error in the knowledge of the value of a non-estimated parameter will increase the actual uncertainty of the estimated state of the system beyond the level formally indicated by the covariance of an estimator that neglects errors or uncertainty in that parameter. The equations for SPCF covariance evolution are obtained in a fashion similar to the derivation approach taken with standard (i.e. linearized or extended) consider parameterized Kalman filters (c.f. [5]). While in [4] the SPCF and linear-theory consider filter (LTCF) were applied to an illustrative linear dynamics/linear measurement problem, in the present work examines the SPCF as applied to nonlinear sequential consider covariance analysis, i.e. in the presence of nonlinear dynamics and nonlinear measurements. A simple SPCF for orbit determination, exemplifying an algorithm hosted in the guidance, navigation and control (GN&C) computer processor of a hypothetical robotic spacecraft, was implemented, and compared with an identically-parameterized (standard) extended, consider-parameterized Kalman filter. The onboard filtering scenario examined is a hypothetical spacecraft orbit about a small natural body with imperfectly-known mass. The formulations, relative complexities, and performances of the filters are compared and discussed.

Lisano, Michael E.↗

Wind profile estimation from point to point laser distortion data

The author's results on the problem of using laser distortion data to estimate the wind profile along the path of the beam are presented. A new model for the dynamics of the index of refraction in a non-constant wind is developed. The model agrees qualitatively with theoretical predictions for the index of refraction statistics in linear wind shear, and is approximated by the predictions of Taylor's hypothesis in constant wind. A framework for a potential in-flight experiment is presented, and the estimation problem is discussed in a maximum likelihood context.

Leland, Robert↗

FREPS - A forced response prediction system for turbomachinery blade rows

FREPS (Forced REsponse Predicition System) is a software system that integrates structural dynamic, steady and unsteady aerodynamic analyses to efficiently predict the forced dynamic stresses of turbomachinery blades to aerodynamic and mechanical excitations. The program performs flutter analysis also. The FREPS system uses a modal approach for aeroelastic analysis. The structural dynamic analysis is based on MSC/NASTRAN, the steady aerodynamic analysis is based on potential theory and the unsteady aerodynamic analysis is based on a linearization of the non-uniform potential mean flow. The capabilities of the program are described and illustrated by application to the High Pressure Oxygen Turbopump turbine of the Space Shuttle Main Engine.

Murthy, Durbha V.↗

Non-Linear System Identification for Aeroelastic Systems with Application to Experimental Data

Representation and identification of a non-linear aeroelastic pitch-plunge system as a model of the NARMAX class is considered. A non-linear difference equation describing this aircraft model is derived theoretically and shown to be of the NARMAX form. Identification methods for NARMAX models are applied to aeroelastic dynamics and its properties demonstrated via continuous-time simulations of experimental conditions. Simulation results show that (i) the outputs of the NARMAX model match closely those generated using continuous-time methods and (ii) NARMAX identification methods applied to aeroelastic dynamics provide accurate discrete-time parameter estimates. Application of NARMAX identification to experimental pitch-plunge dynamics data gives a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗

Elastic continua in high frequency excitation field

The response of elastic continua to high-frequency excitations is decomposed in two parts: 'slow' motion which practically remains unchanged during a vanishingly small period of time, and 'fast' motions whose mean value during this period is negligible but whose energy contribution is essential. After such a decomposition the 'slow' and 'fast' motions become non-linearly coupled by the corresponding governing equations. This coupling leads to an additional 'effective' potential energy which changes the 'mean' stiffness characteristics. The results can be used for dynamical stiffening of flexible structural elements, for a temporary increase of their stiffness in the course of occasional loads to prevent buckling or wrinkling.

Zak, M.↗

Real-Time Simulation of the X-33 Aerospace Engine

This paper discusses the development and performance of the X-33 Aerospike Engine RealTime Model. This model was developed for the purposes of control law development, six degree-of-freedom trajectory analysis, vehicle system integration testing, and hardware-in-the loop controller verification. The Real-Time Model uses time-step marching solution of non-linear differential equations representing the physical processes involved in the operation of a liquid propellant rocket engine, albeit in a simplified form. These processes include heat transfer, fluid dynamics, combustion, and turbomachine performance. Two engine models are typically employed in order to accurately model maneuvering and the powerpack-out condition where the power section of one engine is used to supply propellants to both engines if one engine malfunctions. The X-33 Real-Time Model is compared to actual hot fire test data and is been found to be in good agreement.

Aguilar, Robert↗

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

Shih, Tsan-Hsing↗

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗

On the Importance of the Dynamics of Discretizations

It has been realized recently that the discrete maps resulting from numerical discretizations of differential equations can possess asymptotic dynamical behavior quite different from that of the original systems. This is the case not only for systems of Ordinary Differential Equations (ODEs) but in a more complicated manner for Partial Differential Equations (PDEs) used to model complex physics. The impact of the modified dynamics may be mild and even not observed for some numerical methods. For other classes of discretizations the impact may be pronounced, but not always obvious depending on the nonlinear model equations, the time steps, the grid spacings and the initial conditions. Non-convergence or convergence to periodic solutions might be easily recognizable but convergence to incorrect but plausible solutions may not be so obvious - even for discretized parameters within the linearized stability constraint. Based on our past four years of research, we will illustrate some of the pathology of the dynamics of discretizations, its possible impact and the usage of these schemes for model nonlinear ODEs, convection-diffusion equations and grid adaptations.

Sweby, Peter K.↗

Optimal Fitting of Non-linear Detector Pulses with Nonstationary Noise

Optimal extraction of pulses of constant known shape from a time series with stationary noise is well understood and widely used in detection applications. Applications where high resolution is required over a wide range of input signal amplitudes use much of the dynamic range of the sensor. The noise will in general vary over this signal range, and the response may be a nonlinear function of the energy input. We present an optimal least squares procedure for inferring input energy in such a detector with nonstationary noise and nonlinear energy response.

Fixsen, D. J.↗

Impact of Non-Gaussian Error Volumes on Conjunction Assessment Risk Analysis

An understanding of how an initially Gaussian error volume becomes non-Gaussian over time is an important consideration for space-vehicle conjunction assessment. Traditional assumptions applied to the error volume artificially suppress the true non-Gaussian nature of the space-vehicle position uncertainties. For typical conjunction assessment objects, representation of the error volume by a state error covariance matrix in a Cartesian reference frame is a more significant limitation than is the assumption of linearized dynamics for propagating the error volume. In this study, the impact of each assumption is examined and isolated for each point in the volume. Limitations arising from representing the error volume in a Cartesian reference frame is corrected by employing a Monte Carlo approach to probability of collision (Pc), using equinoctial samples from the Cartesian position covariance at the time of closest approach (TCA) between the pair of space objects. A set of actual, higher risk (Pc >= 10 (exp -4)+) conjunction events in various low-Earth orbits using Monte Carlo methods are analyzed. The impact of non-Gaussian error volumes on Pc for these cases is minimal, even when the deviation from a Gaussian distribution is significant.

Ghrist, Richard W.↗