Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PHASE SPACE”

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 739 records · Page 41

Study of Integrable and Quasi-Integrable Sextupole Lattice

In order to maximize beam lifetime in circular particle accelerators, the nonlinear beam optics are optimized to maximize the dynamic aperture of the beam. The dynamic aperture (DA), which is a 6-D phase space volume of stable trajectories, depends on the strength of the nonlinearities in the machine, and is calculated via particle tracking. Current DA optimization processes include multi-objective genetic algorithm optimizers, and relies on minimizing the magnitudes of resonance driving terms (RDT), which are calculated from the nonlinear contribution to the one-turn-map. The process of searching through the parameter space for an ideal combination that maximizes DA is computationally strenuous. By setting up the sextupole channel such that it is resembles a symplectic integrator of a smooth Hamiltonian, with only a few sextupoles we are able to closely reproduce phase space trajectories of a smooth Hamiltonian up to the hyperbolic point. No chaos and resonances are observed if phase advance per one sextupole magnet in the channel does not exceed ~0.12x2 pi. Therefore, an important property of the suggested approach is the intrinsic elimination of the resonances, and minimization of corresponding RDTs.

43 PARTICLE ACCELERATORS↗

Measurement of the transverse momentum spectrum of the Higgs boson produced in pp collisions at $ \sqrt{s}=8 $ TeV using $H \to WW$ decays

The cross section for Higgs boson production in pp collisions is studied using the H → W$^{+}$W$^{−}$ decay mode, followed by leptonic decays of the W bosons to an oppositely charged electron-muon pair in the final state. The measurements are performed using data collected by the CMS experiment at the LHC at a centre-of-mass energy of 8 TeV, corresponding to an integrated luminosity of 19.4 fb$^{−1}$. The Higgs boson transverse momentum (p$_{T}$) is reconstructed using the lepton pair p$_{T}$ and missing p$_{T}$. The differential cross section times branching fraction is measured as a function of the Higgs boson p$_{T}$ in a fiducial phase space defined to match the experimental acceptance in terms of the lepton kinematics and event topology. The production cross section times branching fraction in the fiducial phase space is measured to be 39 ± 8 (stat) ± 9 (syst) fb. The measurements are found to agree, within experimental uncertainties, with theoretical calculations based on the standard model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

BEAM OPTIMIZATION STUDY FOR AN X-RAY FEL OSCILLATOR AT THE LCLS-II

The 4 GeV LCLS-II superconducting linac with high repetition beam rate enables the possibility to drive an X-Ray FEL oscillator at harmonic frequencies *. Compared to the regular LCLS-II machine setup, the oscillator mode requires a much longer bunch length with a relatively lower current. Also a flat longitudinal phase space distribution is critical to maintain the FEL gain since the X-ray cavity has extremely narrow bandwidth. In this paper, we study the longitudinal phase space optimization including shaping the initial beam from the injector and optimizing the bunch compressor and dechirper parameters. We obtain a bunch with a flat energy chirp over 400 fs in the core part with current above 100 A. The optimization was based on LiTrack and Elegant simulations using LCLS-II beam parameters.

Qin, Weilun↗

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗

The Design of a High Charge Polarized Preinjector for the Electron-Ion Collider

The design of the electron pre-injector of the Electron-Ion Collider (EIC) project to generate 4 x 7 nC bunch has been advancing to meet the requirements for injection into the Rapid Cycling Synchrotron (RCS). The major challenges are high charge transport and achieving small energy spread from 3 GHz traveling-wave plate(TWP). The designed preinjector includes the polarized electron source, bunching section, TWP Linac, zigzag phase space manipulation and spin rotator. In this report, we will discuss the RF frequency selection and the way to reduce energy spread down to 0.2% by longitudinal phase space manipulate. We will also report the results of beamline simulation using space charge code and the conceptual design of spin rotator.

Wang, E.↗

Differential memory in the trilinear model magnetotail

The previously proposed concept of 'differential memory' is quantitatively demonstrated using an idealized analytical model of particle dynamics in the magnetotail geometry. In this model (the 'trilinear' tail model) the magnetotail is divided into three regions. The particle orbits are solved exactly in each region, thus reducing the orbit integration to an analytical mapping. It is shown that the trilinear model reproduces the essential phase space features of the earlier model (Chen and Palmadesso, 1986), possessing well-defined entry and exit regions, and stochastic, integrable (regular), and transient orbits, occupying disjoint phase space regions. Different regions have widely separated characteristic time scales corresponding to different types of particle motion. Using the analytical model, the evolution of single-particle distribution functions is calculated.

Chen, James↗

Beyond RMS: Understanding the Evolution of Beam Distributions in High Intensity Linacs

Understanding the evolution of beams with space charge is crucial to design and operation of high intensity linacs. While the community holds a broad understanding of the mechanisms leading to emittance growth and halo formation, there is outstanding discrepancy between measurements and beam evolution models that precludes prediction of halo losses. This may be due in part to insufficient information of the initial beam distribution. This talk will describe work at the SNS Beam Test Facility to directly measure the 6D beam distribution. Full-and-direct 6D measurement has revealed hidden but physically significant dependence between the longitudinal distribution and transverse coordinates. This nonlinear correlation is driven by space charge and reproduced by self-consistent simulation of the RFQ. Omission of this interplane correlation, common when bunches are reconstructed from lower-dimensional measurements, degrades downstream predictions. This talk will also describe the novel diagnostics supporting this work. This includes ongoing improvements to efficiency of the 6D phase space measurement as well as recent achievement of six orders of dynamic range in 2D phase space.

Ruisard, Kiersten↗

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗

Shearless effective barriers to chaotic transport induced by even twin islands in nontwist systems

For several decades now it has been known that systems with shearless invariant tori, nontwist Hamiltonian systems, possess barriers to chaotic transport. These barriers are resilient to breakage under perturbation and therefore regions where they occur are natural places to look for barriers to transport. Here we describe a kind of effective barrier that persists after the shearless torus is broken. Because phenomena are generic, for convenience we study the standard nontwist map (SNM), an area-preserving map that violates the twist condition locally in the phase space. The barrier occurs in nontwist systems when twin even period islands are present, which happens for a broad range of parameter values in the SNM. With a phase space composed of regular and irregular orbits, the movement of chaotic trajectories is hampered by the existence of shearless curves, total barriers, and a network of partial barriers formed by the stable and unstable manifolds of the hyperbolic points. Being a degenerate system, the SNM has twin islands and, consequently, twin hyperbolic points. We show that the structures formed by the manifolds intrinsically depend on period parity of the twin islands. For this even scenario the structure that we call a torus free barrier occurs because the manifolds of different hyperbolic points form an intricate chain atop a dipole configuration and the transport of chaotic trajectories through the chain becomes a rare event. This structure impacts the emergence of transport, the escape basin for chaotic trajectories, the transport mechanism, and the chaotic saddle. The case of odd periodic orbits is different: we find for this case the emergence of transport immediately after the breakup of the last invariant curve, and this leads to a scenario of higher transport, with intricate escape basin boundary and a chaotic saddle with nonuniformly distributed points.

classical mechanics↗

Exact Maps to Second Order for Selected Elements in MAD-X Variables

This report presents calculations for exact maps to second order about an arbitrary orbit for certain elements: drifts, solenoids, and rotations about a transverse axis. These expressions were used in recent updates to MAD-X, and thus the phase space variables used are those of MAD-X. Formulas are given here for the final orbit and its first and second derivatives with respect to the incoming phase space coordinates for a drift, solenoid, and coordinate system rotations about transverse axes. I do not claim these results to be new: these expressions are well-known and are presented here only for reference. This document does not collect all maps that could be expressed exactly, in particular dipoles are omitted.

43 PARTICLE ACCELERATORS↗

Extraction of Drell-Yan Angular Parameters in $pp$ Collisions with a 120 GeV Beam Energy Using a Deep-Learning Unfolding Algorithm

Dilepton production in pp collisions through the Drell-Yan process provides a crucial tool for studying the internal quark-gluon structure of the nucleon. By precisely measuring the $\cos2\phi$ asymmetry, where $\phi$ represents the azimuthal angle of the $l^{+}l^{-}$ pair in the Collins-Soper frame, we can gain valuable insights into the proton’s structure and the transverse momentum ($q_{T}$) dependence of the $\cos2\phi$ asymmetry. SeaQuest, a fixed-target Drell-Yan experiment at Fermilab, involved an unpolarized proton beam colliding with unpolarized LH$_{2}$ and LD$_{2}$ targets. Measurements obtained from experiments typically require corrections for detector inefficiencies, smearing, and acceptance. Traditionally, these corrections involve “unfolding” the detector-level measurements through matrix operations. However, in higher-dimensional phase space, these conventional methods fail to scale effectively. To overcome these limitations, we employ an unbinned unfolding method that utilizes deep neural networks for unfolding higher-dimensional phase space. In this presentation, we will explain the design of the neural network architecture, our training strategies, and outline our plans to achieve conclusive results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Non-ergodic extended phase of the Quantum Random Energy model

The concept of non-ergodicity in quantum many body systems can be discussed in the context of the wave functions of the many body system or as a property of the dynamical observables, such as time-dependent spin correlators. In the former approach the non-ergodic delocalized state is defined as the one in which the wave functions occupy a volume that scales as a non-trivial power of the full phase space. In this work we study the simplest spin glass model and find that in the delocalized non-ergodic regime the spin–spin correlators decay with the characteristic time that scales as non-trivial power of the full Hilbert space volume. The long time limit of this correlator also scales as a power of the full Hilbert space volume. We identify this phase with the glass phase whilst the many body localized phase corresponds to a ’hyperglass’ in which dynamics is practically absent. We discuss the implications of these findings to quantum information problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rapid Loss of Radiation Belt Relativistic Electrons by EMIC Waves

How relativistic electrons are lost is an important question surrounding the complex dynamics of the Earth's outer radiation belt. Radial loss to the magnetopause and local loss to the atmosphere are two main competing paradigms. Here on the basis of the analysis of a radiation belt storm event on 27 February 2014, we present new evidence for the electromagnetic ion cyclotron (EMIC) wave-driven local precipitation loss of relativistic electrons in the heart of the outer radiation belt. During the main phase of this storm, the radial profile of relativistic electron phase space density was quasi-monotonic, qualitatively inconsistent with the prediction of radial loss theory. The local loss at low L shells was required to prevent the development of phase space density peak resulting from the radial loss process at high L shells. The rapid loss of relativistic electrons in the heart of outer radiation belt was observed as a dip structure of the electron flux temporal profile closely related to intense EMIC waves. Our simulations further confirm that the observed EMIC waves within a quite limited longitudinal region were able to reduce the off-equatorially mirroring relativistic electron fluxes by up to 2 orders of magnitude within about 1.5 h.

79 ASTRONOMY AND ASTROPHYSICS↗

GalaxyFlow: upsampling hydrodynamical simulations for realistic mock stellar catalogues

ABSTRACT Cosmological N-body simulations of galaxies operate at the level of ‘star particles’ with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogues requires ‘upsampling’ the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase-space density and sample from it. Secondly, we improve on existing upsamplers based on adaptive kernel density estimation (KDE), using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighbourhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalogue. Furthermore, we introduce a novel multimodel classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow more accurately estimates the density of the underlying star particles than methods based on KDE, at the cost of being more computationally intensive.

Lim, Sung Hak (ORCID:0000000330981092)↗

Wave kinetics of drift-wave turbulence and zonal flows beyond the ray approximation

Inhomogeneous drift-wave turbulence can be modeled as an effective plasma where drift waves act as quantumlike particles and the zonal-flow velocity serves as a collective field through which they interact. This effective plasma can be described by a Wigner-Moyal equation (WME), which generalizes the quasilinear wave-kinetic equation (WKE) to the full-wave regime, i.e., resolves the wavelength scale. Unlike waves governed by manifestly quantumlike equations, whose WMEs can be borrowed from quantum mechanics and are commonly known, drift waves have Hamiltonians very different from those of conventional quantum particles. This causes unusual phase-space dynamics that is typically not captured by the WKE. We demonstrate how to correctly model this dynamics with the WME instead. Specifically, we report full-wave phase-space simulations of the zonal-flow formation (zonostrophic instability), deterioration (tertiary instability), and the so-called predator-prey oscillations. We also show how the WME facilitates analysis of these phenomena, namely, (i) we show that full-wave effects critically affect the zonostrophic instability, particularly its nonlinear stage and saturation; (ii) we derive the tertiary-instability growth rate; and (iii) we demonstrate that, with full-wave effects retained, the predator-prey oscillations do not require zonal-flow collisional damping, contrary to previous studies. In conclusion, we also show how the famous Rayleigh-Kuo criterion, which has been missing in wave-kinetic theories of drift-wave turbulence, emerges from the WME.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shortened Booster Bunch Lengths

Booster typically produces beam bunches that have a phase space profile with a low momentum spread but high time spread, correlating to longer bunch lengths. The goal of this study is to produce beam bunches that have a phase space profile with a high momentum spread but low time spread or narrow bunch lengths. This can be achieved by slowly increasing the RF cavities amplitudes (Adiabatic Excitation) or by oscillating the amplitudes (Quadratic Bunch Rotation, QBR) until the desired bunch lengths are achieved. Results show that when shortening the bunches using QBR, it showed the bunches had a length of 0.6ns (within 1 sigma) but beam losses were being incurred in the 8 GeV line which were a result of beam dispersion.

Abdelhamid, Maan↗

Near-specular reflection of ions at quasi-parallel shocks

One-dimensional hybrid simulations and a semianalytical model of the shock front are employed to investigate the source regions in the incident ion phase space of reflected and transmitted ions, the relative importance of the electric and magnetic forces in the reflection of incident ions, and how these characteristics change in time. The phase space origin of reflected particles and the fraction of incident ions that reflect are found to depend on the electromagnetic field structure of the shock front at the time the ions encounter it. The reflection fraction is maximized when the electric field along the shock normal (Ex) and the noncoplanar magnetic field (By) are at their maximum values. When Ex and By are large, the reflection process produces a beam that is cooler, more dense, and closer to specular than when they are small.

Mckean, M. E.↗