Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space plasma”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Calibration of a versatile multi-energy soft x-ray diagnostic for WEST long pulse plasmas

A compact multi-energy soft x-ray diagnostic is being installed on the W Environment in Steady-state Tokamak (WEST), which was designed and built to test ITER-like tungsten plasma facing components in a long pulse (~1000 s) scenario. The diagnostic consists of a pinhole camera fielded with the PILATUS3 photon-counting Si-based detector (≲100 kpixel). The detector has sensitivity in the range 1.6–30 keV and enables energy discrimination, providing a higher energy resolution than conventional systems with metal foils and diodes with adequate space and time resolution (≲1 cm and 2 ms). The lower-absorption cut-off energy is set independently on each one of the ~100 kpixels, providing a unique opportunity to measure simultaneously the plasma emissivity in multiple energy ranges and deduce a variety of plasma parameters (e.g., T e , n Z , and ΔZ eff ). The energy dependence of each pixel is calibrated here over the range 3–22 keV. The detector is exposed to a variety of monochromatic sources—fluorescence emission from metallic targets—and for each pixel, the lower energy threshold is scanned to calibrate the energy dependence. The data are fit to a responsivity curve (“S-curve”) that determines the mapping between the possible detector settings and the energy response for each pixel. Here, the calibration is performed for three energy ranges: low (2.3–6 keV), medium (4.5–13.5 keV), and high (5.4–21 keV). We determine the achievable energy resolutions for the low, medium, and high energy ranges as 330 eV, 640 eV, and 950 eV, respectively. The main limitation for the energy resolution is found to be the finite width of the S-curve.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Multidisciplinary benchmarks of a conservative spectral solver for the nonlinear Boltzmann equation

The Boltzmann equation describes the evolution of the phase-space probability distribution of classical particles under binary collisions. Approximations to it underlie the basis for several scholarly fields, including aerodynamics and plasma physics. While these approximations are appropriate in their respective domains, they can be violated in niche but diverse applications which require direct numerical solution of the original nonlinear Boltzmann equation. An expanded implementation of the Galerkin–Petrov conservative spectral algorithm is employed to study a wide variety of physical problems. Enabled by distributed precomputation, solutions of the spatially homogeneous Boltzmann equation can be achieved in seconds on modern personal hardware, while spatially-inhomogeneous problems are solvable in minutes. Here, several benchmarks are presented focusing on accuracy compared to both analytic theoretical predictions and other Boltzmann solvers. These benchmarks span several physical domains including weakly ionized plasma, gaseous fluids, and atomic-plasma interaction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Resilient stellarator divertor characteristics in the Helically Symmetric eXperiment

Resilient divertor features connected to open chaotic edge structures in the Helically Symmetric eXperiment are investigated. For the first time, an expanded vessel wall was considered that would give space for implementation of a physical divertor target structure. The analysis was done for four different magnetic configurations with very different chaotic plasma edges. A resilient plasma wall interaction pattern was identified across all configurations. This manifests as qualitatively very similar footprint behavior across the different plasma equilibria. Overall, the resilient field lines of interest with high connection length L C lie within a helical band along the wall for all configurations. This resiliency can be used to identify the best location of a divertor. The details of the magnetic footprint’s resilient helical band is subject to specific field line structures which are linked to the penetration depth of field lines into the plasma and directly influence the heat and particle flux patterns. The differences arising from these details are characterized by introducing a new metric, the minimum radial connection min $(\delta_N)$ of a field line from the last closed flux surface. The relationship, namely the deviation from a scaling law, between min $(\delta_N)$ and L C of the field lines in the plasma edge field line behavior suggests that the field lines are associated with structures such as resonant islands, cantori, and turnstiles. This helps determine the relevant magnetic flux channels based on the radial location of these chaotic edge structures and the divertor target footprint. These details will need to be taken into account for resilient divertor design.

divertor↗
DOE OSTI2024

Analysis of hydrogen and lithium isotopes using Laser-Induced Breakdown Spectroscopy

The detection and quantification of hydrogen (1,2,3H) and lithium (6,7Li) isotopes are critical to several energy and defense application areas, including hydrogen storage, nuclear forensics, and safeguards/non-proliferation. In this context, laser-induced breakdown spectroscopy (LIBS) is very promising and it comes with certain advantages such as rapid detection, and standoff capability. Although LIBS provides experimental simplicity and is capable of detecting all elements in the periodic table in any phase (solid, liquid, or gas), there exist certain challenges for isotopic analyses of H and Li due to spectral broadening, the presence of closely spaced fine and hyperfine structures, and line distortion effects (e.g., self-absorption and self-reversal) present in laser-produced plasmas. This article reports recent developments in H and Li isotopic analysis using LIBS, and existing challenges.

(020.3260) Isotope shifts, (140.3440) Laser-induce↗
DOE OSTI2021

Photoinjector Generation of High-Charge Magnetized Beams for Electron-Cooling Applications

The recently announced Electron-Ion Collider (EIC), to be built at Brookhaven National Laboratory, will collide high-energy ion-beams and electron beams to pursue nuclear physics research on quark-gluon plasma's. The facility performances (luminosity) would be improved if the ion-beam emittance degradation is mitigated via a phase-space-cooling technique. One potential cooling method uses a bright electron beam to cool the ion beams. The cooling rate for this electron-cooling method depends on the transverse emittance of the cooling electron beam and could benefit from using a beam with significant canonical angular momentum dubbed as a magnetized beam. This research focuses on simulation and experimental generation and characterization of high charge magnetized electron beams with parameters comparable to those required for electron cooling at an electron-ion collider. The experiment uses the 50-MeV photo-injector available at the Fermilab Accelerator Science and Technology (FAST) Facility. The measurements are benchmarked against simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Synchronization of dust acoustic waves in a forced Korteweg–de Vries–Burgers model

The synchronization of dust acoustic waves to an external periodic source is studied in the framework of a driven Korteweg–de Vries–Burgers equation that takes into account the appropriate nonlinear and dispersive nature of low-frequency waves in a dusty plasma medium. For a spatiotemporally varying source term, the system is shown to demonstrate harmonic (1:1) and superharmonic (1:2) synchronized states. Here, the existence domains of these states are delineated in the form of Arnold tongue diagrams in the parametric space of the forcing amplitude and forcing frequency and their resemblance to some past experimental results is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Velocity-space signatures of shock-drift acceleration at quasi-perpendicular collisionless shocks

The shock-drift acceleration of ions at quasi-perpendicular shocks is a well-known kinetic mechanism for the acceleration of a small fraction of incoming ions to high energy. Here, we use a suite of sixteen hybrid simulations of quasi-perpendicular collisionless shocks over the range of Alfvén Mach number 4.3 ≤ M A ≤ 15.8 (corresponding to a range of fast magnetosonic Mach numbers 2.6 ≤ M f ≤ 9.4) and shock-normal angle 45°≤ θ Bn ≤ 90° to identify the velocity-space signature of shock-drift acceleration using the field-particle correlation technique. We show that the features of the ion velocity distribution in the shock foot and ramp regions can be clearly interpreted by analysis of the single-particle trajectory of a reflected ion through the full 3D-3V phase space. The characteristic features of the velocity-space signature of shock-drift acceleration remain qualitatively robust over the full parameter range of our simulations, providing a potential means for its identification using single-point spacecraft measurements. At higher Alfvén Mach numbers M A ≳8 (M f ≳5), kinetic instabilities generate fluctuations of the electromagnetic fields within the shock transition region, leading to fluctuations in and smearing out of the resulting velocity-space signatures, but the signature remains generally robust and identifiable. The results on the shock-drift acceleration of ions presented here represent a novel means to determine more completely the partitioning of upstream bulk flow kinetic energy into plasma heating, particle acceleration, and electromagnetic fields in collisionless shocks.

Howes, Gregory G. [Univ. of Iowa, Iowa City, IA (U↗
DOE OSTI2023

Autonomous sputter synthesis of thin film nitrides with composition controlled by Bayesian optimization of optical plasma emission

Autonomous experimentation has emerged as an efficient approach to accelerate the pace of material discovery. Although instruments for autonomous synthesis have become popular in molecular and polymer science, solution processing of hybrid materials, and nanoparticles, examples of autonomous tools for physical vapor deposition are scarce yet important for the semiconductor industry. Here, we report the design and implementation of an autonomous workflow for sputter deposition of thin films with controlled composition, leveraging a highly automated sputtering reactor custom-controlled by Python, optical emission spectroscopy (OES), and a Bayesian optimization algorithm. We modeled film composition, measured by x-ray fluorescence, as a linear function of plasma emission lines monitored during co-sputtering from elemental Zn and Ti targets in an N 2 and Ar atmosphere. A Bayesian control algorithm, informed by OES, navigates the space of sputtering power to fabricate films with user-defined compositions by minimizing the absolute error between desired and measured optical emission signals. We validated our approach by autonomously fabricating Zn x Ti 1-x N y films that deviate from the targeted cation composition by a relative ±3.5%, even for 15 nm thin films, demonstrating that the proposed approach can reliably synthesize thin films with a specific composition and minimal human interference. Moreover, the proposed method can be extended to more difficult synthesis experiments where plasma intensity lines depend non-linearly on pressure, or the elemental sticking coefficients strongly depend on the substrate temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anomalous Retrograde Drifts in Obstructed Magnetron Microdischarges: a Consequence of a Field Reversal in the Anode Sheath?

The goal of this project was to carry out experiments at the Plasma Collaborative Research Facility (Princeton Plasma Physics Laboratory) aimed at confirming or refuting the hypothesis that an electric field reversal occurs in the near-anode region of a strongly obstructed planar magnetron microdischarges. Magnetron discharges are magnetized plasmas with E x B field topologies, and variants of such discharges are widely used in materials processing. The topologies are also similar to that of Hall thrusters used in space propulsion. Understanding the structure of these discharges will aid in developing reliable simulations that are intended to facilitate the further development of such technologies, as the field reversals, presumed to be due to strong axial-driven electron transport, provides further insight into transport properties. The possible existence of such a field reversal is hinted at because of the retrograde-propagating (i.e., negative E x B) spokes or gradient-driven drift instabilities seen in the high speed videos and segmented anode probing of the plasma. The presence of this field reversal should be reflected in the ion dynamics, particularly, in an ion drift towards the anode. The proposed experiments to be conducted at the PCRF are to characterize this anode-directed ion drift using state-of-the-art laser-induced fluorescence velocimetry. In particular, we seek to see if the near-anode ion velocity reverses direction when the external circuit is varied to alter the electron current demand imposed at the anode. Such an ion velocity reversal would confirm the mechanism, which is that the highly obstructed (small gapped) discharge produces strong axial gradients in plasma density, driving more electrons towards the anode then that demanded by the external circuit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Edge fluctuation measurements in EDA H-mode and QCE plasmas in ASDEX Upgrade using the correlation electron cyclotron emission diagnostic

The requirement of a plasma solution without large edge-localized modes (ELMs) for future fusion devices motivates the study of plasmas without Type I ELMs. The enhanced D–α (EDA) and quasi-continuous exhaust (QCE) scenarios are two such Type I ELM free plasmas, but their underlying physics principles are not fully understood. A database of EDA and QCE plasmas from ASDEX Upgrade was compiled and analyzed using the correlation electron cyclotron emission diagnostic in order to study differences between properties of the quasi-coherent mode (QCM) between these two regimes. Investigations into measured spectral properties of the QCM, radiated temperature fluctuations δT rad /T rad , and radial correlation lengths L c showed no differences between the two regimes. Spectral properties using the cross-power spectral density, G xy were similar between regimes, with the values measured to be a central frequency, $f$ $QCM$ ≈ 20 - 50 kHz; spectral width, Δ $f$ ≈ 5 - 40 kHz; δT rad /T rad ≈ 4-6%; and L c ≈ 1-2 cm. Investigations into a quantity representing the relative bandwidth Δ$f$/$f$ $QCM$ revealed differences between EDA and QCE plasmas when plotted radially, suggesting some type of interaction between the QCM in QCE plasmas with the surrounding broadband turbulence that is not present in EDA phases. Additionally, investigations into various drive terms for hypothesized instability drivers for the QCM are performed showing no difference in the parameter spaces of collisionality and the ion temperature, electron temperature, or electron pressure scale lengths between each regime indicating that there is no observable difference between the drive behind the QCM in each regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus

A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve gauge symmetries and their associated conservation laws via the momentum map. We numerically demonstrate this time invariance of the momentum map and its usefulness in establishing precise initial conditions with a desired initial electric field and fixed background charge. The restriction of this canonical, finite element Poisson structure to the 1X2P $1\frac {1}{2}$ -dimensional phase space is also considered and simulated numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Characterizing core and edge turbulence regimes with fluctuation imaging diagnostics in Wendelstein 7-X

Two density fluctuation imaging systems, phase contrast imaging (PCI) and gas puff imaging (GPI) measure spatially resolved density fluctuations with high time resolution throughout the core plasma (PCI) and in the scrape-off layer (GPI) of the Wendelstein 7-X (W7-X) stellarator. Both systems combined give a comprehensive overview of overall fluctuation levels, spectral properties such as their distribution in frequency and wavenumber space as well as their spatial distribution. These tools are used to assess changes in density turbulence in three representative discharges that transition into stable divertor detachment by different strategies (impurity seeding, density ramping and power starvation). Several general trends are identified when the radiated power fraction is systematically increased: In the plasma edge, the line emission observed by GPI shifts radially inward with a drop in electron temperature, and normalized intensity fluctuation profiles follow this inward shift. Skewness and kurtosis of these edge fluctuations are reduced, indicating a reduction of large intermittent transport events, and poloidal phase velocities decrease in magnitude. These observations are consistent with a reduced power input into the plasma edge and a general reduction of turbulent activity. Core density fluctuation levels remain nearly constant in the impurity seeding scenario, indicating that detachment does not significantly impact turbulence there. However, a strong reduction in the dominant outboard fluctuation phase velocity is observed that deviates from the previous interpretation of neoclassical radial electric field changes, showing that the core plasma is not completely unaffected. In the density ramp and power starvation scenarios, undesirable and irregular large-scale events arise clearly in both diagnostic systems as the radiative fraction is increased. Impurity seeding therefore seems to be a promising strategy on W7-X to achieve detachment without significantly altering core turbulence, especially when targeting a specific operating point in core density and heating power.

Wendelstein 7-X↗
DOE OSTI2026

Supramolecular Assembly of Multielement Ribbon-like Structures Derived from Halide Perovskites

Halide perovskites are crucial materials with broad applications owing to their exceptional optoelectronic properties. Vacancy-ordered double perovskites, featuring highly tunable transition metal sites, enable controllable optoelectronic properties through multielement compositional design. Here, in this study, we introduced 18-crown-6 into the vacancy-ordered double perovskites system and developed two-dimensional ribbon-like single crystals (18C6@K) 2 {PtSnTeIrRe} 1 Cl 6 via an antisolvent supramolecular assembly method, demonstrating morphology modulation through multielement composition design. The crystals crystallize in the centrosymmetric space group $P\bar{1}$. The dumbbell-shaped structural units (crown ether@A)2MX6 pack along a and b axes to form a 2-dimensional (2D) monolayer, and these monolayers further stack along the c axis to generate the ribbon-like single crystals. Energy-dispersive X-ray spectroscopy (EDX) qualitatively confirmed the uniform distribution of the five transition metals throughout the crystal, while inductively coupled plasma atomic emission spectroscopy (ICP–AES) quantitatively verified their atomic ratios. We further investigated the origin of the morphology, distinct from the previously reported cube-like single crystals with the $R\bar{3}$ space group. When acetonitrile was used as the solvent, three-dimensional crystals with $R\bar{3}$ symmetry were obtained, whereas dimethylformamide (DMF) was essential for forming two-dimensional ribbon-like single crystals. The essential role of DMF could be ascribed to its capability to maintain a higher concentration of the building blocks. Moreover, Ir 4+ and Pt 4+ cations also played critical roles in inducing the two-dimensional ribbon-like morphology. The three-element (18C6@K) 2 {PtSnTe}1Cl 6 single crystals exhibited bright yellow emission under 375 nm laser excitation, demonstrating the tunability of optoelectronic properties of this class of material.

Zhu, Heqing [University of California, Berkeley, C↗
DOE OSTI2021

HIBP Applications to Advance Understanding of Plasma Transport Physics (2011-2017) (and subsequently titled) Development of Beam Measurements to Advance Understanding of Transients and Improve Validation (2017-2020) (Final Technical Report)

This is the Final Technical Report for the DoE Measurement Innovation (previously Diagnostic Development) grant DE-SC0006077 titled HIBP Applications to Advance Understanding of Plasma Transport Physics (2011-2017), and subsequently titled Development of Beam Measurements to Advance Understanding of Transients and Improve Validation (2017-2020). The grant was funded by the Office of Fusion Energy Sciences (FES) for the period 1 May 2011–17 July 2020. Work performed through the grant has resulted in diagnostic innovations and new measurement capabilities that will improve understanding of transport, magnetic equilibrium, and electric fields in magnetically confined plasmas. We have developed novel hardware and advanced techniques to expand beam-based diagnostic capabilities, and extend measurements to new operating regimes. Some have been realized using a Heavy Ion Beam Probe (HIBP) having traditional features, while other applications have joined our innovations with subsets of advantageous HIBP features. This work has established smaller, more economical detection systems and, in doing so, may enable substantial extension of scenarios in which beam-based diagnostics are deployed. Key Measurement Innovations include: Computer Models and Simulations of HIBP Applications-Computer modeling and simulating of HIBP applications is central to predicting interaction of beam particles with plasmas, designing diagnostic systems, and anticipating measurement characteristics. We have extended the capability and improved measurement fidelity of the HIBP diagnostic on the Madison Symmetric Torus (MST) reversed field pinch (RFP); simulated feasibility of HIBP operation in the Helically Symmetric eXperiment (HSX) optimized stellarator; and investigated the plausibility of HIBP measurements in the ASDEX Upgrade tokamak. We have also developed a new technique that uses ion optics to model a finite phase-space beam. It yields more realistic (than traditional method) estimates of sample volume characteristics, which influence measurement resolution and sensitivity. Extension of Diagnostic Capabilities and Measurements on the MST RFP-Capabilities have been extended using the first and only HIBP installation on an RFP. Challenges associated with operation of the diagnostic on MST include the three-dimensional nature of particle trajectories, temporal topology and amplitude changes in the plasma equilibrium, and strong particle and radiative emission from the plasma. We have addressed these challenges through development of detailed and higher precision diagnostic simulations, calibrated subtraction of noise to resolve secondary ion signals, and implementation of hardware that enables higher fidelity measurement of plasma fluctuations. Development of New Hardware and Measurement Techniques - We have advanced beam-based measurement capabilities via development of new detectors and techniques to determine the poloidal magnetic flux ψ in the plasma, and enable placement of detectors close to the plasma and thus allow a smaller and less expensive diagnostic. We have demonstrated measurement of ψ, an ability that was made possible by development of highly effective noise reduction techniques. We also modeled the effect of non-ideal and instrumentation effects by developing a virtual (simulated) diagnostic. Application of beam-based diagnostics on various devices and multiple magnetic configurations (e.g. stellarator, tokamak, and reversed field pinch) enables investigation of critical physics issues and measurements resulting in the broad parameter space data needed to test and validate theory & modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Compositionally-driven surface nanostructuring on refractory compositionally complex alloys under low energy helium bombardment

Additive manufacturing enables user-defined control of the compositional complexity, opening new design spaces for complex concentrated alloys (CCAs). Refractory CCAs may offer enhanced performance in the divertor region of a fusion reactor environment where plasma-facing materials will be subject to high temperatures, low energy He and D particles, and 14 MeV neutrons. In this work, specimens with the nominal composition of NbTaMoTi, NbTaMo, NbTaTi, and NbTa were fabricated via directed energy deposition (DED) then bombarded with 40 eV He ions to a fluence of 2x10 26 m −2 at ∼ 1000 K. Post irradiation, the surface morphology and composition were examined with electron microscopy and x-ray photoelectron spectroscopy to offer information on the spatial distribution of surface nano-structuring. Sub-surface He bubbles driving the nano-structuring were examined with electron microscopy techniques. Analysis indicates the local composition directly influences the He bubble and tendril size, while the region with the highest complexity showed the shortest nano-structuring. Molecular Dynamics simulations complement the experimental results, showing comparable helium bubble growth and migration as a function of compositional complexity. This work demonstrates the compositional dependence of surface nanostructure formation of compositionally complex alloys, important for future design of complex plasma facing materials in fusion reactors.

Complex concentrated alloys↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2023

An implicit-in-time DPG formulation of the 1D1V Vlasov-Poisson equations

Efficient solution of the Vlasov equation, which can be up to six-dimensional, is key to the simulation of many difficult problems in plasma physics. The discontinuous Petrov-Galerkin (DPG) finite element methodology provides a framework for the development of stable (in the sense of Ladyzhenskaya–Babuška–Brezzi conditions) finite element formulations, with built-in mechanisms for adaptivity. While DPG has been studied extensively in the context of steady-state problems and to a lesser extent with space-time discretizations of transient problems, relatively little attention has been paid to time-marching approaches. In the present work, we study a first application of time-marching DPG to the Vlasov equation, using backward Euler for a Vlasov-Poisson discretization. We demonstrate adaptive mesh refinement for two problems: the two-stream instability problem, and a cold diode problem. Furthermore, we believe the present work is novel both in its application of unstructured adaptive mesh refinement (as opposed to block-structured adaptivity, which has been studied previously) in the context of Vlasov-Poisson, as well as in its application of DPG to the Vlasov-Poisson system. We also discuss extensive additions to the Camellia library in support of both the present formulation as well as extensions to higher dimensions, Maxwell equations, and space-time formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Extracting dynamical frequencies from invariants of motion in finite-dimensional nonlinear integrable systems

Integrable dynamical systems play an important role in many areas of science, including accelerator and plasma physics. An integrable dynamical system with n degrees of freedom possesses n nontrivial integrals of motion, and can be solved, in principle, by covering the phase space with one or more charts in which the dynamics can be described using action-angle coordinates. To obtain the frequencies of motion, both the transformation to action-angle coordinates and its inverse must be known in explicit form. However, no general algorithm exists for constructing this transformation explicitly from a set of n known (and generally coupled) integrals of motion. In this study we describe how one can determine the dynamical frequencies of the motion as functions of these n integrals in the absence of explicitly known action-angle variables, and we provide several examples.

43 PARTICLE ACCELERATORS↗