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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Nonlinear fishbone dynamics in spherical tokamaks

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. Our results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturates due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. Furthermore, a substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Finally, our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Nonlinear fishbone dynamics in spherical tokamaks

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. The results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturates due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. A substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Analytical nonlinear collisional dynamics of near-threshold eigenmodes

The nonlinear dynamics of isolated, near-threshold waves is studied in the presence of strong resonant particle scattering. In this limit, a simple closed-form analytical solution is found and compared with numerical simulations, which shows reasonable agreement when the wave evolves quasi-steadily to saturation, even for moderate tokamak-relevant collisionality levels. As a result, the solution can be useful in verifying codes that deal with Alfvénic instabilities and thermal plasma turbulence in fusion plasmas, as well as a rapid means for predicting and analyzing experimental outcome.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗
NASA NTRS1983

Nonlinear structural dynamics analysis using a modified modal method

The procedure for predicting the nonlinear dynamic response of structural components subjected to a step loading is presented. The procedure is a modified modal method that involves a change of dependent variables from the unknown nodal degrees of freedom of the finite element model of the structure to a smaller set of generalized coordinates. This change of dependent variables uses a combination of the nonlinear static solution and some selected vibration mode shapes. The vibration mode shapes correspond to the eigenvectors obtained by solving a standard free vibration eigenvalue problem wherein the stiffness matrix is expanded about the nonlinear static solution. A strategy is also presented for determining which and how many vibration mode shapes to include in the transformation. The effect of inaccurate representation of the spatial distribution of the applied load on the nonlinear dynamic response is discussed for two classes of structural behavior. Application of the procedure to structures which exhibit a stiffening behavior and to those with a softening behavior is presented.

Knight, N. F., Jr.↗
DOE OSTI2021

Constrained Block Nonlinear Neural Dynamical Models

Neural network modules conditioned by known priors can be effectively trained and combined to represent systems with nonlinear dynamics. This work explores a novel formulation for data-efficient learning of deep control-oriented nonlinear dynamical models by embedding local model structure and constraints. The proposed method consists of neural network blocks that represent input, state, and output dynamics with constraints placed on the network weights and system variables. For handling partially observable dynamical systems, we utilize a state observer neural network to estimate the states of the system's latent dynamics. We evaluate the performance of the proposed architecture and training methods on system identification tasks for three nonlinear systems: a continuous stirred tank reactor, a two tank interacting system, and an aerodynamics body. Models optimized with a few thousand system state observations accurately represent system dynamics in open loop simulation over thousands of time steps from a single set of initial conditions. Experimental results demonstrate an order of magnitude reduction in open-loop simulation mean squared error for our constrained, block-structured neural models when compared to traditional unstructured and unconstrained neural network models.

Skomski, Elliott↗
NASA NTRS1983

Proposed solution methodology for the dynamically coupled nonlinear geared rotor mechanics equations

The equations which describe the three-dimensional motion of an unbalanced rigid disk in a shaft system are nonlinear and contain dynamic-coupling terms. Traditionally, investigators have used an order analysis to justify ignoring the nonlinear terms in the equations of motion, producing a set of linear equations. This paper will show that, when gears are included in such a rotor system, the nonlinear dynamic-coupling terms are potentially as large as the linear terms. Because of this, one must attempt to solve the nonlinear rotor mechanics equations. A solution methodology is investigated to obtain approximate steady-state solutions to these equations. As an example of the use of the technique, a simpler set of equations is solved and the results compared to numerical simulations. These equations represent the forced, steady-state response of a spring-supported pendulum. These equations were chosen because they contain the type of nonlinear terms found in the dynamically-coupled nonlinear rotor equations. The numerical simulations indicate this method is reasonably accurate even when the nonlinearities are large.

Mitchell, L. D.↗
DOE OSTI2019

An alternative approach for MLE calculation in nonlinear continuous dynamic systems

As an important metric to tell if a nonlinear dynamic system has a singular attractor or divergent trajectory, the maximal Lyapunov exponent (MLE) can be calculated from either system models or time series of state variable measurement. However, in the real world, due to inaccurate models, measurement noise, and the fact that sometimes state variables cannot be measured directly, it is very difficult to get an accurate MLE, which limits its application in, for example, in prediction of a nonlinear physical system (e. g. power systems, etc.) behavior. To overcome these factors, this paper proposed a trajectory-estimation-based MLE calculation approach. The proposed approach addressed how to calculate the MLE when state variables cannot be accessed directly, and uncertainties in system models, as well as noise in measurements. The simulation results shows that the proposed approach is able to handle well the nonlinear measurement functions between state variables and measurements, and get better results than pure model-based approaches or measurement-based approaches in front of measurement noise and model uncertainties.

maximal Lyapunov exponents, nonlinear differential↗
NASA NTRS1989

Parallel processors and nonlinear structural dynamics algorithms and software

A nonlinear structural dynamics finite element program was developed to run on a shared memory multiprocessor with pipeline processors. The program, WHAMS, was used as a framework for this work. The program employs explicit time integration and has the capability to handle both the nonlinear material behavior and large displacement response of 3-D structures. The elasto-plastic material model uses an isotropic strain hardening law which is input as a piecewise linear function. Geometric nonlinearities are handled by a corotational formulation in which a coordinate system is embedded at the integration point of each element. Currently, the program has an element library consisting of a beam element based on Euler-Bernoulli theory and trianglar and quadrilateral plate element based on Mindlin theory.

Belytschko, Ted↗
NASA NTRS1986

Parallel processors and nonlinear structural dynamics algorithms and software

A nonlinear structural dynamics program with an element library that exploits parallel processing is under development. The aim is to exploit scheduling-allocation so that parallel processing and vectorization can effectively be treated in a general purpose program. As a byproduct an automatic scheme for assigning time steps was devised. A rudimentary form of the program is complete and has been tested; it shows substantial advantage can be taken of parallelism. In addition, a stability proof for the subcycling algorithm has been developed.

Belytschko, T.↗
NASA NTRS2002

Dynamical Approach Study of Spurious Numerics in Nonlinear Computations

The last two decades have been an era when computation is ahead of analysis and when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems and non-traditional fields. Ensuring a higher level of confidence in the predictability and reliability (PAR) of these numerical simulations could play a major role in furthering the design, understanding, affordability and safety of our next generation air and space transportation systems, and systems for planetary and atmospheric sciences, and in understanding the evolution and origin of life. The need to guarantee PAR becomes acute when computations offer the ONLY way of solving these types of data limited problems. Employing theory from nonlinear dynamical systems, some building blocks to ensure a higher level of confidence in PAR of numerical simulations have been revealed by the author and world expert collaborators in relevant fields. Five building blocks with supporting numerical examples were discussed. The next step is to utilize knowledge gained by including nonlinear dynamics, bifurcation and chaos theories as an integral part of the numerical process. The third step is to design integrated criteria for reliable and accurate algorithms that cater to the different multiscale nonlinear physics. This includes but is not limited to the construction of appropriate adaptive spatial and temporal discretizations that are suitable for the underlying governing equations. In addition, a multiresolution wavelets approach for adaptive numerical dissipation/filter controls for high speed turbulence, acoustics and combustion simulations will be sought. These steps are corner stones for guarding against spurious numerical solutions that are solutions of the discretized counterparts but are not solutions of the underlying governing equations.

Yee, H. C.↗
DOE OSTI2019

Nonlinear transport, dynamic ordering, and clustering for driven skyrmions on random pinning

Using numerical simulations, we examine the nonlinear dynamics of skyrmions driven over random pinning. For weak pinning, the skyrmions depin elastically, retaining sixfold ordering; however, at the onset of motion there is a dip in the magnitude of the structure factor peaks due to a decrease in positional ordering, indicating that the depinning transition can be detected using the structure factor even within the elastic depinning regime. At higher drives, the moving skyrmion lattice regains full ordering. For increasing pinning strength, we find a transition from elastic to plastic depinning that is accompanied by a sharp increase in the depinning threshold due to the proliferation of topological defects, similar to the peak effect found at the elastic to plastic depinning transition in superconducting vortex systems. For strong pinning and strong Magnus force, the skyrmions in the moving phase can form a strongly clustered or phase-separated state with highly modulated skyrmion density, similar to that recently observed in continuum-based simulations for strong disorder. As the Magnus force is decreased, the density phase-separated state crosses over to a dynamically phase-separated state with uniform density but with flow localized in bands of motion, while in the strongly damped limit, both types of phase-separated states are lost. In the strong pinning limit, we find highly nonlinear velocity-force curves in the transverse and longitudinal directions, along with distinct regions of negative differential conductivity in the plastic flow regime. The negative differential conductivity is absent in the overdamped limit. In conclusion, the Magnus force is responsible for both the negative differential conductivity and the clustering effects since it causes faster moving skyrmions to partially rotate around slower moving or pinned skyrmions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Multi-level optimization with the koopman operator for data-driven, domain-aware, and dynamic system security

Cyber-Physical Systems (CPSs) like the power grid are critically important but also increasingly vulnerable; ensuring reliable system operation in the face of disruptions is becoming more and more challenging. Multi-Level Optimization (MLO) is a powerful way to model adversarial interactions, which naturally makes it applicable to studying CPS security. However, MLO typically does not address underlying system dynamics, and incorporating nonlinear dynamics is generally infeasible. In this paper, we show how to combine MLO with the Koopman Operator (KO) to remedy this. The KO maps nonlinear dynamics to a lifted space in which those dynamics are linear, thus making it ideal for use with MLO. Moreover, the structure of the KO also provides convenient ways to incorporate domain knowledge into the data-driven process of learning the KO representation of a given system. Here we then demonstrate the use of MLO-KO on a small example problem taken from the power grid domain, discuss the scalability and computational cost of MLO-KO, and identify future research directions for this work.

42 ENGINEERING↗
NASA NTRS1989

Algorithms and software for nonlinear structural dynamics

The objective of this research is to develop efficient methods for explicit time integration in nonlinear structural dynamics for computers which utilize both concurrency and vectorization. As a framework for these studies, the program WHAMS, which is described in Explicit Algorithms for the Nonlinear Dynamics of Shells (T. Belytschko, J. I. Lin, and C.-S. Tsay, Computer Methods in Applied Mechanics and Engineering, Vol. 42, 1984, pp 225 to 251), is used. There are two factors which make the development of efficient concurrent explicit time integration programs a challenge in a structural dynamics program: (1) the need for a variety of element types, which complicates the scheduling-allocation problem; and (2) the need for different time steps in different parts of the mesh, which is here called mixed delta t integration, so that a few stiff elements do not reduce the time steps throughout the mesh.

Belytschko, Ted↗