Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 289 records · Page 16

Second harmonic Bessel-Gauss beam shaping with elliptic axicon aberrations

The second harmonic (SH) of an axicon generated Bessel-Gauss beam is created through the nonlinear interaction of photons with a crystal, resulting in the energy doubling of the output photons. In this work, we show experimentally that in addition to frequency doubling, the SH of Bessel-Gauss beams under asymmetric aberrations from an elliptic axicon exhibit intriguing beam formation. Particularly, the central region of the SH beam profile is composed of two central spots of various geometries surrounded by nested ellipses; one of which is the configuration of two central gamma dots with similar radius knotted by nested ellipses for a zeroth-order Bessel-Gauss pump. These SH beams consistently maintain their spatial profile throughout propagation, reminiscing the behavior of screw dislocations in wave patterns. Our numerical simulations produce beam dynamics consistent with that of experiments and further implicate the remarkable interweaving of bright spots with dark vortices. This is especially noticeable when the beams dynamically oscillate along the optical axis, resulting in the genesis of spatially polarized beams with a knotted framework. The insights gained from our study establish a novel paradigm for exploring interactions of Bessel-like beams with vortex dynamics. This, in turn, has the potential to spark innovations in optical applications, fostering new methodologies to harness and manipulate complex light structures. Our experimental findings could spur new methods of generating logical states of light and new opportunities for material processing control. Published by the American Physical Society 2025

Wang, Tianhong (ORCID:0000000169104201)↗

Prediction and experimental observation of damage dependent damping in laminated composite beams

The equations of motion are developed for laminated composite beams with load-induced matrix cracking. The damage is accounted for by utilizing internal state variables. The net result of these variables on the field equations is the introduction of both enhanced damping, and degraded stiffness. Both quantities are history dependent and spatially variable, thus resulting in nonlinear equations of motion. It is explained briefly how these equations may be quasi-linearized for laminated polymeric composites under certain types of structural loading. The coupled heat conduction equation is developed, and it is shown that an enhanced Zener damping effect is produced by the introduction of microstructural damage. The resulting equations are utilized to demonstrate how damage dependent material properties may be obtained from dynamic experiments. Finaly, experimental results are compared to model predictions for several composite layups.

Allen, D. H.↗

Widening and distortion of the particle energy distribution by chromaticity in quasi-isochronous rings

This paper reports the observation of beam energy widening and distortion from Gaussian due to a nonvanishing horizontal chromaticity at a quasi-isochronous storage ring. The result originates from an average path-length dependence on the betatron oscillation amplitudes, which is intimately correlated to the transverse chromaticities. It is the first experimental validation of the impact of such a nonlinear transverse-longitudinal coupling effect on the equilibrium beam characteristics in a storage ring. The results could be important for quasi-isochronous rings, steady-state microbunching, nonscaling fixed-field alternate gradient accelerators, etc., where very small momentum compaction or large chromaticity is required.

43 PARTICLE ACCELERATORS↗

Ion heating by strong electrostatic ion cyclotron turbulence

A theory of the ion heating due to electrostatic ion cyclotron (EIC) waves in the auroral zone is presented. Due to the slowly convecting nature of the EIC mode, quasi-linear plateau formation cannot stabilize the waves, and growth occurs until the nonlinear mechanisms of ion resonance broadening and electron trapping provide saturation. The large amplitude and coherent nature of the resulting wave imply that quasi-linear theory provides only a lower limit to the ion heating. An upper bound on the heating rate is derived using a time-average model of ion dynamics in the coherent waves. The effects of ion heating in the presence of the magnetic gradient force and parallel electric fields are considered, with the result that perpendicular energies over 100 eV are easily attainable from a 1 eV source plasma. Perpendicular heating in the absence of a parallel electric field yields conical ion distributions, which in the presence of an electric field become field-aligned beams.

Lysak, R. L.↗

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗

The dynamics and control of large-flexible space structures, part 10

A mathematical model is developed to predict the dynamics of the proposed orbiting Spacecraft Control Laboratory Experiment (SCOLE) during the station keeping phase. The equations of motion are derived using a Newton-Euler formulation. The model includes the effects of gravity, flexibility, and orbital dynamics. The control is assumed to be provided to the system through the Shuttle's three torquers, and through six actuators located by pairs at two points on the mast and at the mass center of the reflector. The modal shape functions are derived using the fourth order beam equation. The generic mode equations are derived to account for the effects of the control forces on the modal shape and frequencies. The equations are linearized about a nominal equilibrium position. The linear regulator theory is used to derive control laws for both the linear model of the rigidized SCOLE as well as that of the actual SCOLE including the first four flexible modes. The control strategy previously derived for the linear model of the rigidized SCOLE is applied to the nonlinear model of the same configuration of the system and preliminary single axis slewing maneuvers conducted. The results obtained confirm the applicability of the intuitive and appealing two-stage control strategy which would slew the SCOLE system, as if rigid to its desired position and then concentrate on damping out the residual flexible motions.

Bainum, Peter M.↗

Digital Miniature Cathode Ray Magnetometer

In this study, we introduce the concept and construction of an innovative Digital Miniature Cathode Ray Magnetometer designed for the precise detection of magnetic fields. This device addresses several limitations inherent to magnetic probes such as D.C. offset, nonlinearity, temperature drift, sensor aging, and the need for frequent recalibration, while capable of operating in a wide range of magnetic fields. The core principle of this device involves the utilization of a charged particle beam as the sensitivity medium. The system leverages the interaction of an electron beam with a scintillator material, which then emits visible light that is captured by an imager. The emitted scintillation light is captured by a CMOS sensor. This sensor not only records the scintillation light but also accurately determines the position of the electron beam, providing invaluable spatial information crucial for magnetic field mapping. The key innovation lies in the combination of electron beam projection, CMOS imager scintillation-based detection, and digital image signal processing. By employing this synergy, the magnetometer achieves remarkable accuracy, sensitivity and dynamic range. The precise position registration enabled by the CMOS sensor further enhances the device’s utility in capturing complex magnetic field patterns, allowing for 2D field mapping. In this work, the optimization of the probe’s performance is tailored for applications related to the characterization of insertion devices in light sources, including undulators.

Turqueti, Marcos↗

Pulsed acceleration in solar flares

We study the nonlinear dynamics of particle acceleration in solar flares by analyzing the time series of various quasi-periodic radio signatures during flares. In particular we present the radio and hard X-ray data of three flares which suppport the following tentative conclusions: (1) Particle acceleration and injection into magnetic structures occurs intrinsically in a pulsed mode (with a typical period of 1-2 s), produced by a single, spatially coherent, nonlinear system, rather than by a stochastic system with many spatially independent components ('statistical flare' produced by a fragmented primary energy release). (2) The nonlinear (quasi-periodic) mode of pulsed particle acceleration and injection into a coronal loop can be stabilized by phase locking with an MHD wave (oscillation) mode, if both periods are close to each other. (3) Pulsed injection of electron beams into a coronal loop may trigger nonlinear relaxational oscillations of wave-particle interactions. This is particularly likely when the limit cycles of both systems are similar.

Aschwanden, Markus J.↗

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

Nonlinear Optics Measurements in IOTA

Nonlinear integrable optics is a recently proposed accelerator lattice design approach which allows to generate an amplitude dependent tune shift which is needed in high brightness accelerators to mitigate fast coherent instabilities. Whereas usually octupoles are used to achieve this task, this concept allows doing so without exciting any resonances, in turn preventing any particle loss. The concept is based around a special magnet design, together with specific constraints on the optics of the accelerator. To study such a system, the Integrable Optics Test Accelerator (IOTA) was recently constructed and commissioned at Fermilab. For the assessment of the performance of this concept, good knowledge of the optics and the (non-)linear dynamics without the special magnet is of key importance. As such, measurements were conducted in the IOTA ring, using the captured turn-by-turn data by the beam position monitors after excitation to infer quantities such as amplitude detuning and resonance driving terms. In this note, first results of these measurements are presented.

43 PARTICLE ACCELERATORS↗

Vibrational sum frequency generation spectroscopy reveals the inertness of chromium oxide (001) surfaces

Nanoengineered metal oxides such as Cr(III)-oxide (chromia) films have diverse potential applications in corrosion inhibition, remediation, energy generation, catalysis, data storage, and biological and environmental systems. Concerns about material degradation or oxidation to toxic chromate necessitate an understanding of chromia/aqueous interfaces, beginning with their hydroxylation and hydration behavior. Vibrational sum-frequency generation spectroscopy (vSFG) provides specific molecular-level information about water at the oxide/aqueous junction with high surface selectivity. To overcome the strong absorber problem typical of certain metal oxides in the UV–visible range during nonlinear optical studies, we employed molecular beam epitaxy to deposit transparent, epitaxial nanofilms of Cr 2 O 3 with (001) crystalline orientation on sapphire (Al 2 O 3 (001)) substrates, as confirmed by atomic force microscopy and X-ray diffraction. vSFG spectra of the air and water interfaces of the Cr 2 O 3 (001) films reveal hydroxyl features corresponding to both dissociated and molecular water on the surface. In contrast to the dangling hydroxyls found on the bare Al 2 O 3 (001) substrate, the hydroxyl groups on the deposited Cr 2 O 3 (001) nanofilm do not readily undergo isotopic H/D exchange when exposed to varying forms of D 2 O under ambient conditions. When considering the chemistries of the corresponding trivalent cations in aqueous solution, the finding is at variance with their similar acidities and proton exchange dynamics but consistent with markedly slower inner-sphere water ligand exchange of hexaquo Cr 3+ . This finding challenges the idea that proton and water exchange at the water/Cr₂O₃(001) interface are solely correlated, driven by the strength of metal-bridging oxygen bonds and surface hydroxyl distribution, without direct causation.

Boamah-Agyemang, Mavis D. [Pacific Northwest Natio↗

Explicit modeling of composite plates and beams in the dynamics of multibody systems

The state of the art dynamic response analysis of flexible multibody systems is currently restricted to elastic bodies with homogeneous materials. The requirements for high speed operation has made it necessary to use lightweight multi layered composite bodies in robotic systems and space structure applications. Dynamic modeling and analysis of such systems are particularly important since the effects of body flexibility to the performance are likely to be more pronounced. The eight-noded isoperimetric quadrilateral element with independent rotational and displacement degrees of freedom is extended to laminated composite elements. The element includes an arbitrary number of bonded layers, each of which may have a different thickness. The transverse shear deformation which is a predominant factor in the analysis of laminated composite structures is taken into account in developing the stiffness and mass matrices. The corresponding 3-D mode shapes are then incorporated to the multibody system dynamical equations. Floating body reference frames allow the selection of different boundary conditions, and the dynamical equations contain all the nonlinear interactions between the rigid and elastic motion. Example simulations are presented to illustrate the methods proposed.

Amirouche, F. M. L.↗

Shifting and Splitting of Resonance Lines due to Dynamical Friction in Plasmas

We report a quasilinear plasma transport theory that incorporates Fokker-Planck dynamical friction (drag) and pitch angle scattering is self-consistently derived from first principles for an isolated, marginally unstable mode resonating with an energetic minority species. It is found that drag fundamentally changes the structure of the wave-particle resonance, breaking its symmetry and leading to the shifting and splitting of resonance lines. In contrast, scattering broadens the resonance in a symmetric fashion. Comparison with fully nonlinear simulations shows that the proposed quasilinear system preserves the exact instability saturation amplitude and the corresponding particle redistribution of the fully nonlinear theory. Even in situations in which drag leads to a relatively small resonance shift, it still underpins major changes in the redistribution of resonant particles. This novel influence of drag is equally important in plasmas and gravitational systems. In fusion plasmas, the effects are especially pronounced for fast-ion-driven instabilities in tokamaks with low aspect ratio or negative triangularity, as evidenced by past observations. The same theory directly maps to the resonant dynamics of the rotating galactic bar and massive bodies in its orbit, providing new techniques for analyzing galactic dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bilinear modeling and nonlinear estimation

New methods are illustrated for online nonlinear estimation applied to the lateral deflection of an elastic beam on board measurements of angular rates and angular accelerations. The development of the filter equations, together with practical issues of their numerical solution as developed from global linearization by nonlinear output injection are contrasted with the usual method of the extended Kalman filter (EKF). It is shown how nonlinear estimation due to gyroscopic coupling can be implemented as an adaptive covariance filter using off-the-shelf Kalman filter algorithms. The effect of the global linearization by nonlinear output injection is to introduce a change of coordinates in which only the process noise covariance is to be updated in online implementation. This is in contrast to the computational approach which arises in EKF methods arising by local linearization with respect to the current conditional mean. Processing refinements for nonlinear estimation based on optimal, nonlinear interpolation between observations are also highlighted. In these methods the extrapolation of the process dynamics between measurement updates is obtained by replacing a transition matrix with an operator spline that is optimized off-line from responses to selected test inputs.

Dwyer, Thomas A. W., III↗

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

The Stability of Radiatively Cooling Jets. 2: Nonlinear Evolution

We use two-dimensional time-dependent hydrodynamical simulations to follow the growth of the Kelvin-Helmholtz (K-H) instability in cooling jets into the nonlinear regime. We focus primarily on asymmetric modes that give rise to transverse displacements of the jet beam. A variety of Mach numbers and two different cooling curves are studied. The growth rates of waves in the linear regime measured from the numerical simulations are in excellent agreement with the predictions of the linear stability analysis presented in the first paper in this series. In the nonlinear regime, the simulations show that asymmetric modes of the K-H instability can affect the structure and evolution of cooling jets in a number of ways. We find that jets in which the growth rate of the sinusoidal surface wave has a maximum at a so-called resonant frequency can be dominated by large-amplitude sinusoidal oscillations near this frequency. Eventually, growth of this wave can disrupt the jet. On the other hand, nonlinear body waves tend to produce low-amplitude wiggles in the shape of the jet but can result in strong shocks in the jet beam. In cooling jets, these shocks can produce dense knots and filaments of cooling gas within the jet. Ripples in the surface of the jet beam caused by both surface and body waves generate oblique shock "spurs" driven into the ambient gas. Our simulations show these shock "spurs" can accelerate ambient gas at large distances from the jet beam to low velocities, which represents a new mechanism by which low-velocity bipolar outflows may be driven by high-velocity jets. Rapid entrainment and acceleration of ambient gas may also occur if the jet is disrupted. For parameters typical of protostellar jets, the frequency at which K-H growth is a maximum (or highest frequency to which the entire jet can respond dynamically) will be associated with perturbations with a period of - 200 yr. Higher frequency (shorter period) perturbations excite waves associated with body modes that produce internal shocks and only small-amplitude wiggles within the jet. The fact that most observed systems show no evidence for large-amplitude sinusoidal oscillation leading to disruption is indicative that the perturbation frequencies are generally large, consistent with the suggestion that pro- tostellar jets arise from the inner regions (r less than 1 AU) of accretion disks.

Stone, James M.↗

Feasibility of a multigroup Boltzmann–Fokker–Planck solution for electron beam dose calculations

Legacy nuclear-reactor Boltzmann solvers start clinical deployment as an alternative to Monte Carlo (MC) codes and Fermi–Eyges semiemprical models in radiation oncology treatment planning. Today’s certified clinical solvers are limited to photon beams. In this paper, ELECTR, a state-of-the-art multigroup electron cross sections generation module in NJOY is presented and validated against Lockwood’s calorimetric measurements, EGS-nrc and GEANT-4 for 1–20 MeV unidirectional electron beams. The nuclear-reactor DRAGON-5 solver is upgraded to access the library and solve the Boltzmann–Fokker–Planck (BFP) equation. A variety of heterogeneous radiotherapy and radiosurgery phantom configurations were used for validation purpose. Case studies include a thorax benchmark, that of a typical breast Intra-Operative Radiotherapy and a high-heterogeneity patient-like benchmark. For all beams, 100% of the water voxels satisfied the American Association of Physicists in Medicine accuracy criterion for a BFP-MC dose error below 2%. At least, 97.0% of adipose, muscle, bone, lung, tumor and breast voxels satisfied the 2% criterion. The average BFP-MC relative error was about 0.56% for all voxels, beams and materials combined. By irradiating homogeneous slabs from Z = 1 (hydrogen) to Z = 99 (einsteinium), we reported performance and defects of the CEPXS mode [US. Sandia National Lab., SAND-89-1685] in ELECTR for the entire periodic table. For all Lockwood’s benchmarks, NJOY-DRAGON dose predictions are within the experimental data precision for 98% of voxels.

42 ENGINEERING↗

Coherent homodyne detection for amplified cross-beam electric-field induced second harmonic

The electric-field induced second harmonic (E-FISH) signal sensing is commonly used as a nonlinear optical technique to probe electric fields inside a plasma region. Cross-beam E-FISH is being investigated to improve spatial resolution by defining the interaction volume via a controlled geometry of two overlapping noncolinear optical beams. This drastic reduction in interaction length of the electric field and laser region results in a significant signal reduction. Here, to overcome this signal reduction, we introduce coherent amplification of the cross-beam E-FISH signal by mixing the low E-FISH signal with a phase-locked bright local oscillator. We demonstrate enhancement of the signal. By introducing a local oscillator we can now derive the polarity of the measured electric field through the phase of the homodyne signal. To illustrate the technique, we, for the first time, to the best of our knowledge, measure the magnitude and the direction of the electric field in a cold atmospheric pressure plasma jet, which dynamically follows the profile of the applied bias current.

42 ENGINEERING↗