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 919 records · Page 51

Simple robust control laws for robot manipulators. Part 1: Non-adaptive case

A new class of exponentially stabilizing control laws for joint level control of robot arms is introduced. It has been recently recognized that the nonlinear dynamics associated with robotic manipulators have certain inherent passivity properties. More specifically, the derivation of the robotic dynamic equations from the Hamilton's principle gives rise to natural Lyapunov functions for control design based on total energy considerations. Through a slight modification of the energy Lyapunov function and the use of a convenient lemma to handle third order terms in the Lyapunov function derivatives, closed loop exponential stability for both the set point and tracking control problem is demonstrated. The exponential convergence property also leads to robustness with respect to frictions, bounded modeling errors and instrument noise. In one new design, the nonlinear terms are decoupled from real-time measurements which completely removes the requirement for on-line computation of nonlinear terms in the controller implementation. In general, the new class of control laws offers alternatives to the more conventional computed torque method, providing tradeoffs between robustness, computation and convergence properties. Furthermore, these control laws have the unique feature that they can be adapted in a very simple fashion to achieve asymptotically stable adaptive control.

Wen, J. T.↗

Nonlinear hybrid simulation of internal kink with beam ion effects in DIII-D

In DIII-D sawteething plasmas, long-lived (1,1) kink modes are often observed between sawtooth crashes. The saturated kink modes have two distinct frequencies. The mode with higher frequency transits to a fishbone-like mode with sufficient on-axis neutral beam power. In this work, hybrid simulations with the global kinetic-magnetohydrodynamic (MHD) hybrid code M3D-K have been carried out to investigate the linear stability and nonlinear dynamics of the n = 1 mode with effects of energetic beam ions for a typical DIII-D discharge where both saturated kink mode and fishbone were observed. Linear simulation results show that the n = 1 internal kink mode is unstable in MHD limit. However, with kinetic effects of beam ions, a fishbone-like mode is excited with mode frequency about a few kHz depending on beam pressure profile. The mode frequency is higher at higher beam power and/or narrower radial profile consistent with the experimental observation. Nonlinear simulations have been performed to investigate mode saturation as well as energetic particle transport. The nonlinear MHD simulations show that the unstable kink mode becomes a saturated kink mode after a sawtooth crash. With beam ion effects, the fishbone-like mode can also transit to a saturated kink mode with a small but finite mode frequency. These results are consistent with the experimental observation of saturated kink mode between sawtooth crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗

Bilinear dynamic mode decomposition for quantum control

Abstract Data-driven methods for establishing quantum optimal control (QOC) using time-dependent control pulses tailored to specific quantum dynamical systems and desired control objectives are critical for many emerging quantum technologies. We develop a data-driven regression procedure, bilinear dynamic mode decomposition (biDMD), that leverages time-series measurements to establish quantum system identification for QOC. The biDMD optimization framework is a physics-informed regression that makes use of the known underlying Hamiltonian structure. Further, the biDMD can be modified to model both fast and slow sampling of control signals, the latter by way of stroboscopic sampling strategies. The biDMD method provides a flexible, interpretable, and adaptive regression framework for real-time, online implementation in quantum systems. Further, the method has strong theoretical connections to Koopman theory, which approximates nonlinear dynamics with linear operators. In comparison with many machine learning paradigms minimal data is needed to construct a biDMD model, and the model is easily updated as new data is collected. We demonstrate the efficacy and performance of the approach on a number of representative quantum systems, showing that it also matches experimental results.

97 MATHEMATICS AND COMPUTING↗

Some Comments on the Status of Shell Theory at the End of the 20th Century: Complaints and Correctives

This review is divided into complaints and correctives. Complaints are directed at: sloppy refereeing and editing; authors who fail to read or acknowledge what others have done; the mis-naming or mis-crediting of results; the misunderstanding and misuse of the Kirchhoff hypothesis; inflated claims of accuracy based on overly-simplified benchmark problems; the failure to appreciate the inherent errors in various shell models; the failure to appreciate that the physical response and mathematical structure of shell theory are fundamentally different from 3-dimensional elasticity; and the irrelevance of Cosserat-type theories. Correctives include a simple, straight-forward derivation of a general nonlinear dynamic shell theory with the following features: (1) the equations of motion and kinematics (and those of thermodynamics, if desired) are exact consequences of their 3-dimensional counterparts; (2) there are no asymptotic or series expansions through the thickness; (3) all approximations (including the Kirchhoff Hypothesis) occur in the constitutive relations; (4) in static problems, there is a mixed form of the governing equations involving a mixed-energy density and exhibiting remnants of the well-known static-geometric duality of linear theory which is numerically robust because the limiting cases of nonlinear membrane theory and inextensional bending theory fall out naturally. (These latter two special cases are known to produce numerical nightmares unless treated with great care); and (5) all equations may be expressed in coordinate-free form (although, sometimes, a hybrid form is shown to be superior).

Simmonds, James G.↗

Cylindrically confined pair-ion-electron and pair-ion plasmas having axial sheared flow and radial gradients

The linear and nonlinear dynamics of pair-ion (PI) and pair-ion-electron plasmas (PIE) have been investigated in a cylindrical geometry with a sheared plasma flow along the axial direction having radial dependence. The coupled linear dispersion relation of low frequency electrostatic waves has been presented taking into account the Guassian profile of density and linear gradient of sheared flow. It is pointed out that the quasi-neutral cold inhomogeneous pure pair ion plasma supports only the obliquely propagating convective cell mode. The linear dispersion relation of this mode has been solved using boundary conditions. The nonlinear structures in the form of vortices formed by different waves have been discussed in PI and PIE plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A simulation study on chaotic behavior of light propagation in graded index fiber under random perturbation

Chaotic behavior of optical rays in a graded index fiber is studied with the help of nonlinear dynamical method under random perturbation of refractive index. The study indicates that behavior of the rays evolves from regular to chaotic due to the random perturbation. This causes irregularity of propagation of the rays in the graded index fiber, and then the performance of the fiber is reduced.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hybrid simulation of fishbone instabilities in the EAST tokamak

Hybrid simulations with the global kinetic-magnetohydrodynamic (MHD) code M3D-K have been carried out to investigate the linear stability and nonlinear dynamics of beam-driven fishbone in the experimental advanced superconducting tokamak (EAST) experiment. Linear simulations show that a low frequency fishbone instability is excited at experimental value of beam ion pressure. The mode is mainly driven by low energy beam ions via precessional resonance. Our results are consistent with the experimental measurement with respect to mode frequency and mode structure. When the beam ion pressure is increased to exceed a critical value, the low frequency mode transits to a beta-induced Alfven eigenmode (BAE) with much higher frequency. This BAE is driven by higher energy beam ions. Nonlinear simulations show that the frequency of the low frequency fishbone chirps up and down with corresponding hole-clump structures in phase space, consistent with the Berk-Breizman theory. In addition to the low frequency mode, the high frequency BAE is excited during the nonlinear evolution. Furthermore, for the transient case of beam pressure fraction where the low and high frequency modes are simultaneously excited in the linear phase, only one dominant mode appears in the nonlinear phase with frequency jumps up and down during nonlinear evolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bistable dynamics and Hopf bifurcation in a refined model of early stage HIV infection

Recent clinical studies have shown that HIV disease pathogenesis can depend strongly on many factors at the time of transmission, including the strength of the initial viral load and the local availability of CD4+ T-cells. In this article, a new within-host model of HIV infection that incorporates the homeostatic proliferation of T-cells is formulated and analyzed. Due to the effects of this biological process, the influence of initial conditions on the proliferation of HIV infection is further elucidated. The identifiability of parameters within the model is investigated and a local stability analysis, which displays additional complexity in comparison to previous models, is conducted. The current study extends previous theoretical and computational work on the early stages of the disease and leads to interesting nonlinear dynamics, including a parameter region featuring bistability of infectious and viral clearance equilibria and the appearance of a Hopf bifurcation within biologically relevant parameter regimes.

59 BASIC BIOLOGICAL SCIENCES↗

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗

Linear and nonlinear heavy ion-acoustic waves in a strongly coupled plasma

A theoretical study on the propagation of linear and nonlinear heavy ion-acoustic (HIA) waves in an unmagnetized, collisionless, strongly coupled plasma system has been carried out. The plasma system is assumed to contain adiabatic positively charged inertial heavy ion fluids, nonextensive distributed electrons, and Maxwellian light ions. The normal mode analysis is used to study the linear behaviour. On the other hand, the well-known reductive perturbation technique is used to derive the nonlinear dynamical equations, namely, Burgers equation and Korteweg-de Vries (K-dV) equation. They are also numerically analyzed in order to investigate the basic features of shock and solitary waves. The adiabatic effects on the HIA shock and solitary waves propagating in such a strongly coupled plasma are taken into account. It has been observed that the roles of the adiabatic positively charged heavy ions, nonextensivity of electrons, and other plasma parameters arised in this investigation have significantly modified the basic features (viz., polarity, amplitude, width, etc.) of the HIA solitary/shock waves. The findings of our results obtained from this theoretical investigation may be useful in understanding the linear as well as nonlinear phenomena associated with the HIA waves both in space and laboratory plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The charged inflaton and its gauge fields: preheating and initial conditions for reheating

We calculate particle production during inflation and in the early stages of reheating after inflation in models with a charged scalar field coupled to Abelian and non-Abelian gauge fields. A detailed analysis of the power spectra of primordial electric fields, magnetic fields and charge fluctuations at the end of inflation and preheating is provided. We carefully account for the Gauss constraints during inflation and preheating, and clarify the role of the longitudinal components of the electric field. We calculate the timescale for the back-reaction of the produced gauge fields on the inflaton condensate, marking the onset of non-linear evolution of the fields. We provide a prescription for initial conditions for lattice simulations necessary to capture the subsequent nonlinear dynamics. On the observational side, we find that the primordial magnetic fields generated are too small to explain the origin of magnetic fields on galactic scales and the charge fluctuations are well within observational bounds for the models considered in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Prominence condensation and magnetic levitation in a coronal loop

The results of a model dynamic simulation of the formation and support of a narrow prominence at the apex of a coronal magnetic loop or arcade are described. The condensation process proceeds via an initial radiative cooling and pressure drop, and a secondary siphon flow from the dense chromospheric ends. The antibuoyancy effect as the prominence forms causes a bending of the confining magnetic field, which propagates toward the semirigid ends of the magnetic loop. Thus, a wide magnetic 'hammock' or well (of the normal-polarity Kippenhahn-Schlueter-type) is formed, which supports the prominence at or near the field apex. The simplicity of this 1.5-dimensional model, with its accompanying diagnostics, elucidates the various contributions to the nonlinear dynamics of prominence condensation and levitation.

Van Hoven, G.↗

Bifurcation theory for the L-H transition in magnetically confined fusion plasmas

The mathematical field of bifurcation theory is extended to be applicable to 1-dimensionally resolved systems of nonlinear partial differential equations, aimed at the determination of a certain specific bifurcation. This extension is needed to be able to properly analyze the bifurcations of the radial transport in magnetically confined fusion plasmas. This is of special interest when describing the transition from the low-energy-confinement state to the high-energy-confinement state of the radial transport in fusion plasmas (i.e., the L-H transition), because the nonlinear dynamical behavior during the transition corresponds to the dynamical behavior of a system containing such a specific bifurcation. This bifurcation determines how the three types (sharp, smooth, and oscillating) of observed L-H transitions are organized as function of all the parameters contained in the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic bifurcation and fractal and chaos control of a giant magnetostrictive film-shape memory alloy composite cantilever plate subjected to in-plane harmonic and stochastic excitation

Stochastic bifurcation and fractal and chaos control of a giant magnetostrictive film–shape memory alloy (GMF–SMA) composite cantilever plate subjected to in-plane harmonic and stochastic excitation were studied. Van der Pol items were improved to interpret the hysteretic phenomena of both GMF and SMA, and the nonlinear dynamic model of a GMF–SMA composite cantilever plate subjected to in-plane harmonic and stochastic excitation was developed. The probability density function of the dynamic response of the system was obtained, and the conditions of stochastic Hopf bifurcation were analyzed. The conditions of noise-induced chaotic response were obtained in the stochastic Melnikov integral method, and the fractal boundary of the safe basin of the system was provided. Finally, the chaos control strategy was proposed in the stochastic dynamic programming method. Numerical simulation shows that stochastic Hopf bifurcation and chaos appear in the parameter variation process. The boundary of the safe basin of the system has fractal characteristics, and its area decreases when the noise intensifies. The system reliability was improved through stochastic optimal control, and the safe basin area of the system increased.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modelling Heterogeneity in Cellulose Properties Predicts the Slowdown Phenomenon during Enzymatic Hydrolysis

Fine control of the enzymatic hydrolysis of cellulose is challenging due to complex process dynamics. Conversion rate from cellulose to valuable monomeric products is often significantly slowed down after an initial rapid but short-lived phase. Underlying mechanisms for this process have yet to be fully understood and as a result poorly represented in existing models. Here, we propose a new modelling platform termed Multi-Layered Population Balance Model (ML-PBM), which enables various key aspects of cellulose enzymatic hydrolysis occurring over the entire breakdown process to be captured. As a core component for predicting the slowdown phenomenon, the ML-PBM properly accounts for heterogeneity in cellulose crystallinity and chain lengths across the structural layers of cellulose particles with different morphologies. Beyond a decent quantitative fit to highly nonlinear dynamic experimental data collected across different conditions, the ML-PBM reveals that the rate slowdown phenomenon is potentially due to heterogeneity in cellulose properties coupled with cellulose morphology. Equipped with the unification of various process fundamentals, the ML-PBM is a rational framework with the potential to facilitate sound cellulose engineering.

Ahamed, Firnaaz↗

Stellarator microinstabilities and turbulence at low magnetic shear

Gyrokinetic simulations of drift waves in low-magnetic-shear stellarators reveal that simulation domains comprised of multiple turns can be required to properly resolve critical mode structures important in saturation dynamics. Marginally stable eigenmodes important in saturation of ion temperature gradient modes and trapped electron modes in the Helically Symmetric Experiment (HSX) stellarator are observed to have two scales, with the envelope scale determined by the properties of the local magnetic shear and an inner scale determined by the interplay between the local shear and magnetic field-line curvature. Properly resolving these modes removes spurious growth rates that arise for extended modes in zero-magnetic-shear approximations, enabling use of a zero-magnetic-shear technique with smaller simulation domains and attendant cost savings. Analysis of subdominant modes in trapped electron mode (TEM)-driven turbulence reveals that the extended marginally stable modes play an important role in the nonlinear dynamics, and suggests that the properties induced by low magnetic shear may be exploited to provide another route for turbulence saturation.

Faber, B. J. (ORCID:000000034934400X)↗

Analysis of stochastically forced quasi-periodic attractors

A problem of the analysis of stochastically forced quasi-periodic auto-oscillations of nonlinear dynamic systems is considered. A stationary distribution of random trajectories in the neighborhood of the corresponding deterministic attractor (torus) is studied. A parametric description of quadratic approximation of the quasipotential based on the stochastic sensitivity functions (SSF) technique is given. Using this technique, we analyse a dispersion of stochastic flows near the torus. For the case of two-torus in three-dimensional space, the stochastic sensitivity function is constructed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗