Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear beam 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 163 records · Page 9

The influence of and the identification of nonlinearity in flexible structures

Several models were built at NASA Langley and used to demonstrate the following nonlinear behavior: internal resonance in a free response, principal parametric resonance and subcritical instability in a cantilever beam-lumped mass structure, combination resonance in a parametrically excited flexible beam, autoparametric interaction in a two-degree-of-freedom system, instability of the linear solution, saturation of the excited mode, subharmonic bifurcation, and chaotic responses. A video tape documenting these phenomena was made. An attempt to identify a simple structure consisting of two light-weight beams and two lumped masses using the Eigensystem Realization Algorithm showed the inherent difficulty of using a linear based theory to identify a particular nonlinearity. Preliminary results show the technique requires novel interpretation, and hence may not be useful for structural modes that are coupled by a guadratic nonlinearity. A literature survey was also completed on recent work in parametrically excited nonlinear system. In summary, nonlinear systems may possess unique behaviors that require nonlinear identification techniques based on an understanding of how nonlinearity affects the dynamic response of structures. In this was, the unique behaviors of nonlinear systems may be properly identified. Moreover, more accutate quantifiable estimates can be made once the qualitative model has been determined.

Zavodney, Lawrence D.↗

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of Quadrupole Fringe Fields in the Interaction Region of the Hadron Storage Ring of the Electron Ion Collider

Fringe fields in quadrupole magnets are usually neglected in studies of beam dynamics at accelerators. However, the extreme optical parameters present in the final focus of a collider such as the Electron–Ion Collider (EIC) may give rise to effects that should not be overlooked. The calculation of quadrupole fringe fields presented in this study follows the procedure outlined in Ref. [1], specialized to the case of a straight reference orbit (i.e., with no dipole field component). A right-handed Cartesian coordinate system is employed, with the $z$-axis aligned with the quadrupole axis and $x$ and $y$ denoting the horizontal and vertical transverse coordinates, respectively. The magnetic quadrupole field gradient, $G = \partial B_y / \partial x$, transitions from its peak value inside the quadrupole—where it is nearly independent of the longitudinal coordinate $z$—to zero at some distance beyond the magnet edge. Consequently, $G$ is treated as a function of $z$. The region over which this variation occurs is defined as the quadrupole fringe field region. The study begins with the development of a description of the magnetic field in the fringe region using a power-series expansion in the transverse coordinates $x$and $y$, consistent with the longitudinally varying gradient. A model for the $z$-dependence of the gradient is then proposed and adjusted to reproduce magnetic field data obtained from three-dimensional field calculations. To evaluate the impact of the fringe fields on beam dynamics, the corresponding vector potential is derived and incorporated into a Hamiltonian formulation of particle motion. The significance of the fringe fields is quantified by calculating the amplitude-dependent tune shift from the Hamiltonian. Using linear beam optics parameters of the Hadron Storage Ring (HSR) of the EIC, the tune shift due to the fringe fields of all quadrupole magnets in the IR-6 interaction region is evaluated. Finally, the resulting tune shifts are compared with those arising from other nonlinear field components present in the HSR.

43 PARTICLE ACCELERATORS↗

Structural model optimization using statistical evaluation

The results of research in applying statistical methods to the problem of structural dynamic system identification are presented. The study is in three parts: a review of previous approaches by other researchers, a development of various linear estimators which might find application, and the design and development of a computer program which uses a Bayesian estimator. The method is tried on two models and is successful where the predicted stiffness matrix is a proper model, e.g., a bending beam is represented by a bending model. Difficulties are encountered when the model concept varies. There is also evidence that nonlinearity must be handled properly to speed the convergence.

Collins, J. D.↗

Modeling of composite beams and plates for static and dynamic analysis

The main purpose of this research was to develop a rigorous theory and corresponding computational algorithms for through-the-thickness analysis of composite plates. This type of analysis is needed in order to find the elastic stiffness constants for a plate and to post-process the resulting plate solution in order to find approximate three-dimensional displacement, strain, and stress distributions throughout the plate. This also requires the development of finite deformation plate equations which are compatible with the through-the-thickness analyses. After about one year's work, we settled on the variational-asymptotical method (VAM) as a suitable framework in which to solve these types of problems. VAM was applied to laminated plates with constant thickness in the work of Atilgan and Hodges. The corresponding geometrically nonlinear global deformation analysis of plates was developed by Hodges, Atilgan, and Danielson. A different application of VAM, along with numerical results, was obtained by Hodges, Lee, and Atilgan. An expanded version of this last paper was submitted for publication in the AIAA Journal.

Hodges, Dewey H.↗

Configuration design sensitivity analysis of built-up structures

In the design of complex built-up structures that are made of truss, beam, membrane, shell, and solid, there are five different kinds of design variables: material property, sizing, shape, configuration, and topological variables. Previous research has shown that the improvement in performances obtained by altering the configuration of structural components can be much more significant than those obtained when the geometry is assumed to be fixed. Using the variational approach, a unified design sensitivity was developed for the first three kinds of design variables, and was further extended recently in many structural analysis problems such as nonlinear, structural dynamics, and frequency response analysis. A continuum design sensitivity analysis method is developed for the configuration design variable of built-up structures.

Choi, Kyung K.↗

Dynamic responses of railroad car models to vertical and lateral rail inputs

Simplified dynamic models were applied in a study of vibration in a high-speed railroad car. The mathematical models used were a four-degree-of-freedom model for vertical responses to vertical rail inputs and a ten-degree-of-freedom model for lateral response to lateral or rolling (cross-level) inputs from the rails. Elastic properties of the passenger car body were represented by bending and torsion of a uniform beam. Rail-to-car (truck) suspensions were modeled as spring-mass-dashpot oscillators. Lateral spring nonlinearities approximating certain complicated truck mechanisms were introduced. The models were excited by displacement and, in some cases, velocity inputs from the rails by both deterministic (including sinusoidal) and random input functions. Results were obtained both in the frequency and time domains. Solutions in the time domain for the lateral model were obtained for a wide variety of transient and random inputs generated on-line by an analog computer. Variations in one of the damping properties of the lateral car suspension gave large fluctuations in response over a range of car speeds for a given input. This damping coefficient was significant in reducing lateral car responses that were higher for nonlinear springs for three different inputs.

Sewall, J. L.↗

Load-shortening behavior of an initially curved eccentrically loaded column

To explore the feasibility of using buckled columns to provide a soft support system for simulating a free-free boundary condition in dynamic testing, the nonlinear load-shortening behavior of initially imperfect, eccentrically loaded slender columns is analyzed. Load-shortening curves are obtained for various combinations of load eccentricity and uniform initial curvature and are compared, for reference purposes, with the limiting case of the classical elastica. Results for numerous combinations of initial curvature and load eccentricity show that, over a wide range of shortening, an axially loaded slender column exhibits load-deflection compliance which is of the same order as that of a straight but otherwise identical cantilever beam under lateral tip loading.

Fichter, W. B.↗

Dynamics and efficiency of type III solar radio emission

Existing calculations of nonlinear coupling coefficients for fundamental and harmonic emission via three-wave interactions are refined and used to calculate the conversion efficiency of Langmuir energy into electromagnetic waves when the Langmuir waves have the observed bursty form. Resulting field strengths for harmonic emission are found to be consistent with typical International Sun Earth Explorer 3 (ISSE 3) observations at 1 AU. Fundamental emission at 1 AU is found to proceed only when stimulated by the presence of a source of ion sound waves. However, it is argued that electrostatic decay of Langmuir waves can supply the necessary waves provided the driving electron beam is sufficiently fast. Under these conditions, the predicted fundamental field strengths can account for both the highest and typical fields observed; they also dominate the predicted harmonic fields, consistent with observations. This mechanism is also consistent with previous observations that fundamental emission generally occurs early in type III events, when the beam is fastest. For typical parameters it is shown that neither fundamental nor harmonic emission saturates its respective source instability, contrary to previous assumptions. However, saturation cannot be ruled out under particularly favorable conditions.

Robinson, P. A.↗

The Dynamics of Coupled Oscillator Phase Control

Arrays of coupled oscillators have been proposed as means of realizing high power rf sources via coherent spatial power combining. In such applications, a uniform phase distribution over the aperture is usually desired. However, it has been shown that by detuning some of the oscillators away from the oscillation frequency of the ensemble of oscillators, one may achieve other useful aperture phase distributions. Of particular interest among those achievable are linear phase distributions because these result in steering of the output rf beam away from the broadside direction. The theory describing the behavior of such arrays of coupled oscillators is quite complicated since the phenomena involved are inherently nonlinear. However, a simplified theory has been developed which facilitates intuitive understanding. This simplified theory is based on a "continuum model" in which the aperture phase is represented by a continuous function of the aperture coordinates. A challenging aspect of the development of this theory is the derivation of appropriate boundary conditions at the edges or ends of the array.

Pogorzelski, R. J.↗

Assessment of BPM options for the EIC Beam Accumulator Ring

The electron injection system for the Electron-Ion Collider (EIC) at BNL is designed to provide a beam of polarized electrons, which is crucial for studying the structure of protons and atomic nuclei. The Beam Accumulator Ring (BAR) is a part of the injection chain between the 750 MeV linear accelerator and the Rapid Cycling Synchrotron (RCS), which accelerates the beam up to full energy (5–18 GeV). The functional role of the BAR is to accumulate the charge injected from the linear accelerator in order to achieve the high intensity of the polarized electron beam required by the specifications for injection into the RCS. This objective will be realized through the sequential injection of bunches with a charge of 1.1 nC at a repetition rate of 30 Hz. Once the charge of a single bunch reaches 28 nC, the beam will be extracted from the BAR and injected into the RCS. The beam instrumentation needs to provide reliable measurements with the required accuracy over the dynamic range from 0.1 nC (one tenth of a typical injected bunch charge) to 32 nC – the maximum expected accumulated current. There are 10 Beam Position Monitors (BPMs) in the ring, 1 extra button-electrode assembly for the RF system, and 8 BPMs in the beam transport lines. The BPM locations are shown in Fig. 1, marked by blue rectangles. The ring BPMs will be capable of both average orbit and turn-by-turn measurement modes. Accurate measurement of the beam position with a large horizontal offset requires polynomial correction of the BPM nonlinearity.

43 PARTICLE ACCELERATORS↗

ICIASF '85 - International Congress on Instrumentation in Aerospace Simulation Facilities, 11th, Stanford University, CA, August 26-28, 1985, Record

Developments related to laser Doppler velocimetry are discussed, taking into account a three-component dual beam laser-Doppler-anemometer to be operated in large wind tunnels, a new optical system for three-dimensional laser-Doppler-anemometry using an argon-ion and a dye laser, and a two-component laser Doppler velocimeter by switching fringe orientation. Other topics studied are concerned with facilities, instrumentation, control, hot wire/thin film measurements, optical diagnostic techniques, signal and data processing, facilities and adaptive wall test sections, data acquisition and processing, ballistic instrument systems, dynamic testing and material deformation measurements, optical flow measurements, test techniques, force measurement systems, and holography. Attention is given to nonlinear calibration of integral wind tunnel balances, a microcomputer system for real time digitized image compression, and two phase flow diagnostics in propulsion systems.

Source record↗

Nonlinear strain-displacement relations and flexible multibody dynamics

Dynamics of chains of flexible bodies undergoing large rigid body motions, but small elastic deflections are considered. The role of nonlinear strain-displacement relations in the development of the motion equations correct to first order in elastic deflections is investigated. The general form of these equations linearized only in the small elastic deflections is presented, and the relative significance of various nonlinear terms is studied both analytically and through the use of the numerical simulations. Numerical simulations are performed for a two link chain constrained to move in the plane, subject to hinge torques. Each link is modeled as a thin beam. Slew maneuver simulation results are compared for models with and without properly modeled kinematics of deformation. The goal of this case study is to quantify the importance of the terms in the equations of motion which arise from the inclusion of nonlinear strain-displacement relations. It is concluded that unless the consistently linearized equations in elastic deflections and speeds are available and necessary, the inconsistently (prematurely) linearized equations should be replaced in all cases by ruthlessly linearized equations: equations in which all nonlinear terms involving the elastic deflections and speeds are ignored.

Padilla, Carlos E.↗

The dynamics and control of large flexible space structures, 2. Part A: Shape and orientation control using point actuators

The equations of planar motion for a flexible beam in orbit which includes the effects of gravity gradient torques and control torques from point actuators located along the beam was developed. Two classes of theorems are applied to the linearized form of these equations to establish necessary conditions for controlability for preselected actuator configurations. The feedback gains are selected: (1) based on the decoupling of the original coordinates and to obtain proper damping, and (2) by applying the linear regulator problem to the individual model coordinates separately. The linear control laws obtained using both techniques were evaluated by numerical integration of the nonlinear system equations. Numerical examples considering pitch and various number of modes with different combination of actuator numbers and locations are presented. The independent model control concept used earlier with a discretized model of the thin beam in orbit was reviewed for the case where the number of actuators is less than the number of modes. Results indicate that although the system is controllable it is not stable about the nominal (local vertical) orientation when the control is based on modal decoupling. An alternate control law not based on modal decoupling ensures stability of all the modes.

Bainum, P. M.↗

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↗

Simulation study of Type 2 counterstreaming electrons along auroral field lines

The production of counterstreaming electrons associated with parallel fields along auroral field lines is examined through the use of computer simulation. A 2 1/2-dimensional (two spatial and three velocity dimensions) electrostatic particle algorithm and auroral boundary conditions are used to set up a self-consistent V potential structure. The simulation produces signatures of counterstreaming electrons resembling those observed by the Dynamics Explorer 1 satellite. The main signatures are as follows: (1) the phase space contours of the electron distribution function are elongated along the V-parallel axis, and (2) the energy of electrons streaming in the upward direction is comparable to the energy of the accelerated electron beam. The simulation indicates that a portion of the accelerated electron beam is trapped by large amplitude electrostatic waves produced through the two-stream instability. Strong wave-particle interactions then thermalize the trapped electrons to produce suprathermal electrons streaming in the direction opposite to that of the accelerated electron beam. These results suggest a possible mechanism of producing counterstreaming electron fluxes through nonlinear processes of the two-stream instability.

Wagner, J. S.↗

Transmission Telescope Optical Dynamic Alignment

The Integrated Radio Optical Communication System (iROC) is designed to transmit data between Mars and Earth by means of radio waves at 32.67 GHz (Ka Band) and laser beam at 1550 nm, both transmitted via a combined teletenna. The iROC terminal will provide ‘beaconless’ operations to allow full function from the outer planets. In order to point without the aid of an uplink beacon, the proof of concept presented here is addressing the need for an accurate determination and control of the relative position of the laser beam with respect to a reference star. The experiment presented simulates a surrogate transmission telescope system in laboratory setting and presents the model used in the correction of the outgoing beam. The results of the model show a nonlinear dependence between the outgoing and the reference beam, indicating the necessity of a minimum of two metrology instruments placed along the optical system for increased pointing precision.

optical communication↗

Hybrid model predictive control techniques for safety factor profile and stored energy regulation while incorporating NBI constraints

Abstract A novel hybrid Model Predictive Control (MPC) algorithm has been designed for simultaneous safety factor ( q ) profile and stored energy ( w ) control while incorporating the pulse-width-modulation constraints associated with the neutral beam injection (NBI) system. Regulation of the q -profile has been extensively shown to be a key factor for improved confinement as well as non-inductive sustainment of the plasma current. Simultaneous control of w is necessary to prevent the triggering of pressure-driven magnetohydrodynamic instabilities as the controller shapes the q profile. Conventional MPC schemes proposed for q -profile control have considered the NBI powers as continuous-time signals, ignoring the discrete-time nature of these actuators and leading in some cases to performance loss. The hybrid MPC scheme in this work has the capability of incorporating the discrete-time actuator dynamics as additional constraints. In nonlinear simulations, the proposed hybrid MPC scheme demonstrates improved q -profile+ w control performance for NSTX-U operating scenarios.

Physics↗