Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear 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 451 records · Page 25

Experimental study from linear to chaotic regimes on a terahertz-frequency gyrotron oscillator

Basic wave-particle interaction dynamics from linear to chaotic regimes is experimentally studied on a frequency tunable gyrotron generating THz radiation in continuous mode (200 W) at 263 GHz which will be used for dynamic nuclear polarization nuclear magnetic resonance spectroscopy applications. In the studied system, the nonlinear dynamics associated to the wave-particle interaction is dominated by longitudinal mode competition of a given transverse TE{sub m,p} cavity-mode. This study covers a wide range of control parameter from gyro-traveling wave tube (gyro-TWT) to gyro-backward wave oscillator (gyro-BWO) like interactions for which extensive theoretical studies have been performed in the past on a simplified system. Besides the common route to chaos characterized by period doubling, other routes have been identified among which some are characterized by line-width frequency-broadening on the side-bands. The complex nonlinear dynamics is in good agreement with the theory and the experimental results are discussed on the basis of the prediction obtained with the nonlinear time-dependent self-consistent codes TWANG and EURIDICE both based on a slow-time scale formulation of the self-consistent equations governing the wave-particle dynamics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗

Nonlinear propagating modes beyond the phonons in fluorite-structured crystals

The vibrational energy of crystals is known to propagate in a fixed number of phonon branches allowed by symmetry. In the realm of nonlinear dynamics, however, additional nonlinear propagating modes are possible. Nonlinear propagating modes have unique properties that are important in many disciplines including optical communications, conducting polymers, biology, magnetism, and nuclear physics. Yet, despite the crucial importance of crystal lattice vibrations in fundamental and applied science, such additional propagating modes have not been observed in ordinary crystals. Here, we show that propagating modes exist beyond the phonons in fluorite-structured thoria, urania, and natural calcium fluoride using neutron scattering and first-principles calculations. These modes are observed at temperatures ranging from 5 K up to 1200 K, extend to frequencies 30–40% higher than the maximum phonon frequency, and travel at velocities comparable to or higher than the fastest phonon. The nonlinear origin of the modes is explained in part via perturbation theory, which approximately accounts for nonlinearity. Given that these modes are still clearly observed at 5 K, they are likely an inherent feature of the quantum ground state. The existence of these waves in three-dimensional crystals may have ramifications for material properties.

36 MATERIALS SCIENCE↗

Theoretical and numerical modelling of chaotic electrostatic ion cyclotron (EIC) oscillations by Jerk equation

In the last few years, third order explicit autonomous differential equations, known as jerk equations, have generated great interest as they show features of regular and chaotic motion. In this paper, we have modelled chaotic electrostatic ion cyclotron oscillations using a third order nonlinear ordinary differential equation (ODE) and investigated its nonlinear dynamical properties. The nonlinear ODE has been derived for a plasma system from a two fluid model in the presence of a source term, under the influence of an external magnetic field, which is perpendicular to the direction of the wave vector. It is seen that the equation does not require an external forcing term to obtain chaotic behaviour. The stability of the solutions of the equation has been investigated analytically as well as numerically, and the bifurcation diagram obtained shows a number of interesting phenomena for various regimes of parameters. The coexisting attractors as well as their corresponding basins are shown and the phase space portraits at different conditions are obtained numerically and shown here. The results obtained here are in agreement with preliminary experiments conducted for a similar configuration of a plasma system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extracting collective motions underlying nucleosome dynamics via nonlinear manifold learning

The identification of effective collective variables remains a challenge in molecular simulations of complex systems. Here, we use a nonlinear manifold learning technique known as the diffusion map to extract key dynamical motions from a complex biomolecular system known as the nucleosome: a DNA-protein complex consisting of a DNA segment wrapped around a disc-shaped group of eight histone proteins. We show that without any a priori information, diffusion maps can identify and extract meaningful collective variables that characterize the motion of the nucleosome complex. We find excellent agreement between the collective variables identified by the diffusion map and those obtained manually using a free energy-based analysis. Notably, diffusion maps are shown to also identify subtle features of nucleosome dynamics that did not appear in those manually specified collective variables. For example, diffusion maps identify the importance of looped conformations in which DNA bulges away from the histone complex that are important for the motion of DNA around the nucleosome. As a result, this work demonstrates that diffusion maps can be a promising tool for analyzing very large molecular systems and for identifying their characteristic slow modes.

36 MATERIALS SCIENCE↗

State estimation with small nonlinearities

A variety of techniques is available for estimating the states of nonlinear dynamic systems from noisy data. The differences among several of these procedures in the presence of small dynamic and observational nonlinearities are investigated. Four discrete estimation algorithms are analyzed. The first is a strictly least square estimator, while the others are recursive algorithms similar to the Kalman filter used for estimating the states of linear systems. A group of analytic expressions is developed for the mean and covariance of the error in each of these estimators so that they may be compared without lengthy Monte Carlo simulations. The covariance expressions show that, to first order, all the estimators have the same covariance. Expressions for the means show that each estimator has a different bias. Several examples are carried out demonstrating that the relative magnitudes of the bias errors in the various estimators can be a strong function of such parameters as initial covariances and number of data points. Under some circumstances, more complicated algorithms can have larger biases than smaller ones.

Conrad, B.↗

A comparison of the effects of small nonlinearities on several estimation schemes.

A variety of techniques are available for estimating the states of nonlinear dynamic systems from noisy data. These procedures are generally equivalent when applied to linear systems. This paper investigates the difference between several of these procedures in the presence of small dynamic and observational nonlinearities. In particular, it examines one least square batch processing algorithm, and three recursive algorithms similar to the Kalman filter. To first order, all the estimators have the same covariance. Expressions for the means, however, show that each estimator has a different bias. The examples presented show that the biases can be a strong function of such parameters as initial covariances and number of data points being considered.

Conrad, B.↗

An approximation technique for predicting the transient response of a second order nonlinear equation

An approximation technique has been developed for determining the transient response of a nonlinear dynamic system. The nonlinearities in the system which has been considered appear in the system's dissipation function. This function was expressed as a second order polynomial in the system's velocity. The developed approximation is an extension of the classic Kryloff-Bogoliuboff technique. Two examples of the developed approximation are presented for comparative purposes with other approximation methods.

Laurenson, R. M.↗

Nonlinear systems approach to control system design

Consider some of the control system design methods for plants with nonlinear dynamics. If the nonlinearity is weak relative to the size of the operating region, then the linear methods apply directly. Fixed-gain design may be feasible even for significant nonlinearities. It may be possible to find a single gain which provides adequate control of the linear models at several perturbation points. If the nonlinearity is restricted to a sector, that fact may be used to obtain a fixed-gain controller. Otherwise, a gain may have to be associated with each perturbation point Pi. A gain schedule K(p(v)) is obtained by connecting the perturbation points by a function, say p(v), of the scheduling parameter v (i.e., speed). When the scheduling parameter must be multidimensional, this approach is difficult; the objective is to develop an easier procedure.

Meyer, G.↗

A methodology for designing robust multivariable nonlinear feedback control systems

This paper presents an outline of a methodology for the design of nonlinear dynamic compensators for nonlinear multivariable systems with guarantees of closed-loop stability, robustness, and performance. The method is an extension of the Linear-Quadratic-Gaussian with Loop-Transfer-Recovery (LQG/LTR) methodology for linear systems, thus hinging upon the idea of constructing an approximate nonlinear inverse operator for the plant. A major feature of the method is a unification of both the state-space and Input-Output formulations. The major results of the paper are: (1) Properties of the solution of certain optimal control problems; (2) New results on nonlinear estimation, including a new nonlinear observer with certain guaranteed properties; (3) The Loop Operator Recovery (LOR) procedure, and (4) Inclusion of command-following and disturbance rejection with a goal of a closed-loop stable, robust system.

Grunberg, D. B.↗

A methodology for designing robust multivariable nonlinear control systems

A new methodology is described for the design of nonlinear dynamic controllers for nonlinear multivariable systems providing guarantees of closed-loop stability, performance, and robustness. The methodology is an extension of the Linear-Quadratic-Gaussian with Loop-Transfer-Recovery (LQG/LTR) methodology for linear systems, thus hinging upon the idea of constructing an approximate inverse operator for the plant. A major feature of the methodology is a unification of both the state-space and input-output formulations. In addition, new results on stability theory, nonlinear state estimation, and optimal nonlinear regulator theory are presented, including the guaranteed global properties of the extended Kalman filter and optimal nonlinear regulators.

Grunberg, D. B.↗

Phase-Space Detection of Cyber Events

Energy Delivery Systems (EDS) are a network of processes that produce, transfer and distribute energy. EDS are increasingly dependent on networked computing assets, as are many Industrial Control Systems. Consequently, cyber-attacks pose a real and pertinent threat, as evidenced by Stuxnet, Shamoon and Dragonfly. Hence, there is a critical need for novel methods to detect, prevent, and mitigate effects of such attacks. To detect cyber-attacks in EDS, we developed a framework for gathering and analyzing timing data that involves establishing a baseline execution profile and then capturing the effect of perturbations in the state from injecting various malware. The data analysis was based on nonlinear dynamics and graph theory to improve detection of anomalous events in cyber applications. The goal was the extraction of changing dynamics or anomalous activity in the underlying computer system. Takens' theorem in nonlinear dynamics allows reconstruction of topologically invariant, time-delay-embedding states from the computer data in a sufficiently high-dimensional space. The resultant dynamical states were nodes, and the state-to-state transitions were links in a mathematical graph. Alternatively, sequential tabulation of executing instructions provides the nodes with corresponding instruction-to-instruction links. Graph theorems guarantee graph-invariant measures to quantify the dynamical changes in the running applications. Results showed a successful detection of cyber events.

cyber anomaly detection↗

Formulation of the linear model from the nonlinear simulation for the F18 HARV

The F-18 HARV is a modified F-18 Aircraft which is capable of flying in the post-stall regime in order to achieve superagility. The onset of aerodynamic stall, and continued into the post-stall region, is characterized by nonlinearities in the aerodynamic coefficients. These aerodynamic coefficients are not expressed as analytic functions, but rather in the form of tabular data. The nonlinearities in the aerodynamic coefficients yield a nonlinear model of the aircraft's dynamics. Nonlinear system theory has made many advances, but this area is not sufficiently developed to allow its application to this problem, since many of the theorems are existance theorems and that the systems are composed of analytic functions. Thus, the feedback matrices and the state estimators are obtained from linear system theory techniques. It is important, in order to obtain the correct feedback matrices and state estimators, that the linear description of the nonlinear flight dynamics be as accurate as possible. A nonlinear simulation is run under the Advanced Continuous Simulation Language (ACSL). The ACSL simulation uses FORTRAN subroutines to interface to the look-up tables for the aerodynamic data. ACSL has commands to form the linear representation for the system. Other aspects of this investigation are discussed.

Hall, Charles E., Jr.↗

Probabilistic learning on manifolds constrained by nonlinear partial differential equations for small datasets

A novel extension of the Probabilistic Learning on Manifolds (PLoM) is presented. It makes it possible to synthesize solutions to a wide range of nonlinear stochastic boundary value problems described by partial differential equations (PDEs) for which a stochastic computational model (SCM) is available and which depend on a vector-valued random control parameter. The cost of a single numerical evaluation of this SCM is assumed to be such that only a limited number of points can be computed for constructing the training dataset (small data). Each point of the training dataset is made up of realizations from a vector-valued stochastic process (the stochastic solution) and the associated random control parameter on which it depends. The presented PLoM constrained by PDE allows for generating a large number of learned realizations of the stochastic process and its corresponding random control parameter. These learned realizations are generated so as to minimize the vector-valued random residual of the PDE in the mean-square sense. Appropriate novel methods are developed to solve this challenging problem. Three applications are presented. The first one is a simple uncertain nonlinear dynamical system with a nonstationary stochastic excitation. The second one concerns the 2D nonlinear unsteady Navier–Stokes equations for incompressible flows in which the Reynolds number is the random control parameter. Here, the last one deals with the nonlinear dynamics of a 3D elastic structure with uncertainties. The results obtained make it possible to validate the PLoM constrained by stochastic PDE but also provide further validation of the PLoM without constraint.

Machine learning↗

Neural Networks for Rapid Design and Analysis

Artificial neural networks have been employed for rapid and efficient dynamics and control analysis of flexible systems. Specifically, feedforward neural networks are designed to approximate nonlinear dynamic components over prescribed input ranges, and are used in simulations as a means to speed up the overall time response analysis process. To capture the recursive nature of dynamic components with artificial neural networks, recurrent networks, which use state feedback with the appropriate number of time delays, as inputs to the networks, are employed. Once properly trained, neural networks can give very good approximations to nonlinear dynamic components, and by their judicious use in simulations, allow the analyst the potential to speed up the analysis process considerably. To illustrate this potential speed up, an existing simulation model of a spacecraft reaction wheel system is executed, first conventionally, and then with an artificial neural network in place.

Sparks, Dean W., Jr.↗

Application of Contraction Mappings to the Control of Nonlinear Systems

The theoretical and applied aspects of successive approximation techniques are considered for the determination of controls for nonlinear dynamical systems. Particular emphasis is placed upon the methods of contraction mappings and modified contraction mappings. It is shown that application of the Pontryagin principle to the optimal nonlinear regulator problem results in necessary conditions for optimality in the form of a two point boundary value problem (TPBVP). The TPBVP is represented by an operator equation and functional analytic results on the iterative solution of operator equations are applied. The general convergence theorems are translated and applied to those operators arising from the optimal regulation of nonlinear systems. It is shown that simply structured matrices and similarity transformations may be used to facilitate the calculation of the matrix Green functions and the evaluation of the convergence criteria. A controllability theory based on the integral representation of TPBVP's, the implicit function theorem, and contraction mappings is developed for nonlinear dynamical systems. Contraction mappings are theoretically and practically applied to a nonlinear control problem with bounded input control and the Lipschitz norm is used to prove convergence for the nondifferentiable operator. A dynamic model representing community drug usage is developed and the contraction mappings method is used to study the optimal regulation of the nonlinear system.

Killingsworth, W. R., Jr.↗

Physics-Informed Active Learning With Simultaneous Weak-Form Latent Space Dynamics Identification

The parametric greedy latent space dynamics identification (gLaSDI) framework has demonstrated promising potential for accurate and efficient modeling of high-dimensional nonlinear physical systems. However, it remains challenging to handle noisy data. Here, to enhance robustness against noise, we incorporate the weak-form estimation of nonlinear dynamics (WENDy) into gLaSDI. In the proposed weak-form gLaSDI (WgLaSDI) framework, an autoencoder and WENDy are trained simultaneously to discover intrinsic nonlinear latent-space dynamics of high-dimensional data. Compared with the standard sparse identification of nonlinear dynamics (SINDy) employed in gLaSDI, WENDy enables variance reduction and robust latent space discovery, therefore leading to more accurate and efficient reduced-order modeling. Furthermore, the greedy physics-informed active learning in WgLaSDI enables adaptive sampling of optimal training data on the fly for enhanced modeling accuracy. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including viscous and inviscid Burgers' equations, time-dependent radial advection, and the Vlasov equation for plasma physics. With data that contains 5%–10% Gaussian white noise, WgLaSDI outperforms gLaSDI by orders of magnitude, achieving 1%–7% relative errors. Compared with the high-fidelity models, WgLaSDI achieves 121 to 1779x speed-up.

97 MATHEMATICS AND COMPUTING↗