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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

A pepper-pot emittance meter for low-energy heavy-ion beams

A novel emittance meter has been developed to measure the four-dimensional, transverse phase-space distribution of a low-energy ion beam using the pepper-pot technique. A characteristic feature of this instrument is that the pepper-pot plate, which has a linear array of holes in the vertical direction, is scanned horizontally through the ion beam. This has the advantage that the emittance can also be measured at locations along the beam line where the beam has a large horizontal divergence. A set of multi-channel plates, scintillation screen, and ccd camera is used as a position-sensitive ion detector allowing a large range of beam intensities that can be handled. This paper describes the design, construction, and operation of the instrument as well as the data analysis used to reconstruct the four-dimensional phase-space distribution of an ion beam. Measurements on a 15 keV He{sup +} beam are used as an example.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Slip-stacking Dynamics for High-Power Proton Beams at Fermilab

Slip-stacking is a particle accelerator configuration used to store two particle beams with different momenta in the same ring. The two beams are longitudinally focused by two radiofrequency (RF) cavities with a small frequency difference between them. Each beam is synchronized to one RF cavity and perturbed by the other RF cavity. Fermilab uses slip-stacking in the Recycler so as to double the power of the 120 GeV proton beam in the Main Injector. This dissertation investigates the dynamics of slip-stacking beams analytically, numerically and experimentally. In the analytic analysis, I find the general trajectory of stable slip-stacking particles and identify the slip-stacking parametric resonances. In the numerical analysis, I characterize the stable phase-space area and model the particle losses. In particular, I evaluate the impact of upgrading the Fermilab Booster cycle-rate from 15 Hz to 20 Hz as part of the Proton Improvement Plan II (PIP-II). The experimental analysis is used to verify my approach to simulating slip-stacking loss. I design a study for measuring losses from the longitudinal single-particle dynamics of slip-stacking as a function of RF cavity voltage and RF frequency separation. I further propose the installation of a harmonic RF cavity and study the dynamics of this novel slip-stacking configuration. I show the harmonic RF cavity cancels out parametric resonances in slip-stacking, reduces emittance growth during slip-stacking, and dramatically enhances the stable phase-space area. The harmonic cavity is expected to reduce slip-stacking losses to far exceed PIP-II requirements. These results raise the possibility of extending slip-stacking beyond the PIP-II era.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

A time-implicit numerical method and benchmarks for the relativistic Vlasov–Ampere equations

We present a time-implicit numerical method to solve the relativistic Vlasov–Ampere system of equations on a two dimensional phase space grid. The time-splitting algorithm we use allows the generalization of the work presented here to higher dimensions keeping the linear aspect of the resulting discrete set of equations. The implicit method is benchmarked against linear theory results for the relativistic Landau damping for which analytical expressions using the Maxwell-Jüttner distribution function are derived. We note that, independently from the shape of the distribution function, the relativistic treatment features collective behaviours that do not exist in the nonrelativistic case. The numerical study of the relativistic two-stream instability completes the set of benchmarking tests.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

SU-E-T-242: Monte Carlo Simulations Used to Test the Perturbation of a Reference Ion Chamber Prototype Used for Small Fields

Purpose: To study the perturbation due to the use of a novel Reference Ion Chamber designed to measure small field dosimetry (KermaX Plus C by IBA). Methods: Using the Phase-space files for TrueBeam photon beams available by Varian in IAEA-compliant format for 6 and 15 MV. Monte Carlo simulations were performed using BEAMnrc and DOSXYZnrc to investigate the perturbation introduced by a reference chamber into the PDDs and profiles measured in water tank. Field sizes ranging from 1×1, 2×2,3×3, 5×5 cm2 were simulated for both energies with and without a 0.5 mm foil of Aluminum which is equivalent to the attenuation equivalent of the reference chamber specifications in a water phantom of 30×30×30 cm3 and a pixel resolution of 2 mm. The PDDs, profiles, and gamma analysis of the simulations were performed as well as a energy spectrum analysis of the phase-space files generated during the simulation. Results: Examination of the energy spectrum analysis performed shown a very small increment of the energy spectrum at the build-up region but no difference is appreciated after dmax. The PDD, profiles and gamma analysis had shown a very good agreement among the simulations with and without the Al foil, with a gamma analysis with a criterion of 2% and 2mm resulting in 99.9% of the points passing this criterion. Conclusion: This work indicates the potential benefits of using the KermaX Plus C as reference chamber in the measurement of PDD and Profiles for small fields since the perturbation due to in the presence of the chamber the perturbation is minimal and the chamber can be considered transparent to the photon beam.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Lagrangian and Hamiltonian constraints for guiding-center Hamiltonian theories

In this paper, a consistent guiding-center Hamiltonian theory is derived by Lie-transform perturbation method, with terms up to second order in magnetic-field nonuniformity. Consistency is demonstrated by showing that the guiding-center transformation presented here satisfies separate Jacobian and Lagrangian constraints that have not been explored before. A new first-order term appearing in the guiding-center phase-space Lagrangian is identified through a calculation of the guiding-center polarization. It is shown that this new polarization term also yields a simpler expression of the guiding-center toroidal canonical momentum, which satisfies an exact conservation law in axisymmetric magnetic geometries. Finally, an application of the guiding-center Lagrangian constraint on the guiding-center Hamiltonian yields a natural interpretation for its higher-order corrections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

SU-E-T-512: Electromagnetic Simulations of the Dielectric Wall Accelerator

Purpose: To characterize and parametrically study the key components of a dielectric wall accelerator through electromagnetic modeling and particle tracking. Methods: Electromagnetic and particle tracking simulations were performed using a commercial code (CST Microwave Studio, CST Inc.) utilizing the finite integration technique. A dielectric wall accelerator consists of a series of stacked transmission lines sequentially fired in synchrony with an ion pulse. Numerous properties of the stacked transmission lines, including geometric, material, and electronic properties, were analyzed and varied in order to assess their impact on the transverse and axial electric fields. Additionally, stacks of transmission lines were simulated in order to quantify the parasitic effect observed in closely packed lines. Particle tracking simulations using the particle-in-cell method were performed on the various stacks to determine the impact of the above properties on the resultant phase space of the ions. Results: Examination of the simulation results show that novel geometries can shape the accelerating pulse in order to reduce the energy spread and increase the average energy of accelerated ions. Parasitic effects were quantified for various geometries and found to vary with distance from the end of the transmission line and along the beam axis. An optimal arrival time of an ion pulse relative to the triggering of the transmission lines for a given geometry was determined through parametric study. Benchmark simulations of single transmission lines agree well with published experimental results. Conclusion: This work characterized the behavior of the transmission lines used in a dielectric wall accelerator and used this information to improve them in novel ways. Utilizing novel geometries, we were able to improve the accelerating gradient and phase space of the accelerated particle bunch. Through simulation, we were able to discover and optimize design issues with the device at low cost. Funding: Morgridge Institute for Research, Madison WI; Conflict of Interest: Dr. Mackie is an investor and board member at CPAC, a company developing compact accelerator designs similar to those discussed in this work, but designs discussed are not directed by CPAC. Funding: Morgridge Institute for Research, Madison WI; Conflict of Interest: Dr. Mackie is an investor and board member at CPAC, a company developing compact accelerator designs similar to those discussed in this work, but designs discussed are not directed by CPAC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Clathrates and beyond: Low-density allotropy in crystalline silicon

In its common, thermodynamically stable state, silicon adopts the same crystal structure as diamond. Although only a few alternative allotropic structures have been discovered and studied over the past six decades, advanced methods for structure prediction have recently suggested a remarkably rich low-density phase space that has only begun to be explored. The electronic properties of these low-density allotropes of silicon, predicted by first-principles calculations, indicate that these materials could offer a pathway to improving performance and reducing cost in a variety of electronic and energy-related applications. In this focus review, we provide an introduction and overview of recent theoretical and experimental results related to low-density allotropes of silicon, highlighting the significant potential these materials may have for technological applications, provided substantial challenges to their experimental preparation can be overcome.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Equilibrium properties of blackbody radiation with an ultraviolet energy cut-off

We study various equilibrium thermodynamic properties of blackbody radiation (i.e. a photon gas) with an ultraviolet energy cut-off. We find that the energy density, specific heat etc. follow usual acoustic phonon dynamics as have been well studied by Debye. Other thermodynamic quantities like pressure, entropy etc. have also been calculated. The usual Stefan–Boltzmann law gets modified. We observe that the values of the thermodynamic quantities with the energy cut-off is lower than the corresponding values in the theory without any such scale. The phase-space measure is also expected to get modified for an exotic spacetime appearing at Planck scale, which in turn leads to the modification of Planck energy density distribution and the Wien’s displacement law. We found that the non-perturbative nature of the thermodynamic quantities in the SR limit (for both unmodified and modified cases), due to nonanalyticity of the leading term, is a general feature of the theory accompanied with an ultraviolet energy cut-off. We have also discussed the possible modification in the case of Big Bang and the Stellar objects and have suggested a table top experiment for verification in effective low energy case. - Highlights: • We study the equilibrium thermodynamic properties of a photon gas with an ultraviolet energy cut-off. The values of the thermodynamic quantities are lower than the corresponding values in the theory without any such scale. • The usual Stefan–Boltzmann law gets modified. Also, there is a one to one correspondence with photons in Doubly Special Relativity and the acoustic phonons. • The phase space measure modification at Planck scale has been explored. This leads to modification of Planck energy density distribution and the Wien’s displacement law. • We found that non-perturbative nature of the thermodynamic quantities in the SR limit is a general feature of the theory accompanied with an ultraviolet energy cut-off. • At very high energies, the Hubble’s parameter gets a modification in the radiation dominated epoch of big bang, while at low energies, a table top experiment for verification of our theory has been suggested.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Two-nucleon emission in neutrino and electron scattering from nuclei: The modified convolution approximation

The theoretical formalism of inclusive lepton–nucleus scattering in the two-nucleon emission channel is discussed in the context of a simplified approach, the modified convolution approximation. This allows one to write the 2p2h responses of the relativistic Fermi gas as a folding integral of two 1p1h responses with the energies and momenta transferred to each nucleon. The idea behind this method is to introduce different average momenta for the two initial nucleons in the matrix elements of the two-body current, with the innovation that they depend on the transferred energies and momenta. This method treats exactly the two-body phase space kinematics, and reduces the formulae of the response functions from seven-dimensional integrals over momenta to much simpler three-dimensional ones. The applicability of the method is checked by comparing with the full results within a model of electroweak meson-exchange currents. The predictions are accurate enough, especially in the low-energy threshold region where the average momentum approximation works the best.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Synchronous detection of longitudinal and transverse bunch signals at a storage ring

To understand and control dynamics in the longitudinal phase space, time-resolved measurements of different bunch parameters are required. For a reconstruction of this phase space, the detector systems have to be synchronized. This reconstruction can be used for example for studies of the microbunching instability which occurs if the interaction of the bunch with its own radiation leads to the formation of substructures on the longitudinal bunch profile. These substructures can grow rapidly—leading to a sawtooth-like behavior of the bunch. At KARA, we use a fast-gated intensified camera for energy spread studies, Schottky diodes for coherent synchrotron radiation studies as well as electro-optical spectral decoding for longitudinal bunch profile measurements. For a synchronization, a synchronization scheme is used which compensates for hardware delays. Here in this paper, the different experimental setups and their synchronization are discussed and first results of synchronous measurements presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

A Monte Carlo calculation model of electronic portal imaging device for transit dosimetry through heterogeneous media

Purpose: To develop and evaluate a fast Monte Carlo (MC) dose calculation model of electronic portal imaging device (EPID) based on its effective atomic number modeling in the XVMC code. Methods: A previously developed EPID model, based on the XVMC code by density scaling of EPID structures, was modified by additionally considering effective atomic number (Z{sub eff}) of each structure and adopting a phase space file from the EGSnrc code. The model was tested under various homogeneous and heterogeneous phantoms and field sizes by comparing the calculations in the model with measurements in EPID. In order to better evaluate the model, the performance of the XVMC code was separately tested by comparing calculated dose to water with ion chamber (IC) array measurement in the plane of EPID. Results: In the EPID plane, calculated dose to water by the code showed agreement with IC measurements within 1.8%. The difference was averaged across the in-field regions of the acquired profiles for all field sizes and phantoms. The maximum point difference was 2.8%, affected by proximity of the maximum points to penumbra and MC noise. The EPID model showed agreement with measured EPID images within 1.3%. The maximum point difference was 1.9%. The difference dropped from the higher value of the code by employing the calibration that is dependent on field sizes and thicknesses for the conversion of calculated images to measured images. Thanks to the Z{sub eff} correction, the EPID model showed a linear trend of the calibration factors unlike those of the density-only-scaled model. The phase space file from the EGSnrc code sharpened penumbra profiles significantly, improving agreement of calculated profiles with measured profiles. Conclusions: Demonstrating high accuracy, the EPID model with the associated calibration system may be used for in vivo dosimetry of radiation therapy. Through this study, a MC model of EPID has been developed, and their performance has been rigorously investigated for transit dosimetry.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

RECONSTRUCTING THE ACCRETION HISTORY OF THE GALACTIC STELLAR HALO FROM CHEMICAL ABUNDANCE RATIO DISTRIBUTIONS

Observational studies of halo stars during the past two decades have placed some limits on the quantity and nature of accreted dwarf galaxy contributions to the Milky Way (MW) stellar halo by typically utilizing stellar phase-space information to identify the most recent halo accretion events. In this study we tested the prospects of using 2D chemical abundance ratio distributions (CARDs) found in stars of the stellar halo to determine its formation history. First, we used simulated data from 11 “MW-like” halos to generate satellite template sets (STSs) of 2D CARDs of accreted dwarf satellites, which are composed of accreted dwarfs from various mass regimes and epochs of accretion. Next, we randomly drew samples of ∼10{sup 3–4} mock observations of stellar chemical abundance ratios ([α/Fe], [Fe/H]) from those 11 halos to generate samples of the underlying densities for our CARDs to be compared to our templates in our analysis. Finally, we used the expectation-maximization algorithm to derive accretion histories in relation to the STS used and the sample size. For certain STSs used we typically can identify the relative mass contributions of all accreted satellites to within a factor of two. We also find that this method is particularly sensitive to older accretion events involving low-luminosity dwarfs, e.g., ultra-faint dwarfs—precisely those events that are too ancient to be seen by phase-space studies of stars and too faint to be seen by high-z studies of the early universe. Since our results only exploit two chemical dimensions and near-future surveys promise to provide ∼6–9 dimensions, we conclude that these new high-resolution spectroscopic surveys of the stellar halo will allow us to recover its accretion history—and the luminosity function of infalling dwarf galaxies—across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Activation of Water by Thorium Cation: A Guided Ion Beam and Quantum Chemical Study

The reaction of atomic thorium cations with deuterated water as a function of kinetic energy from thermal to 10 eV was studied using guided ion beam tandem mass spectrometry. At thermal energies, both ThO{sup +} + D{sub 2} and DThO{sup +} + D are formed in barrierless exothermic processes and reproduce results in the literature obtained using ion cyclotron resonance mass spectrometry. As the energy is increased, the branching ratio between these two channels changes such that the dominant product changes from ThO{sup +} to DThO{sup +} and back to ThO{sup +}, until ThD{sup +} + OD is energetically available and is the dominant product channel. To help understand these experimental results, a variety of theoretical approaches were tried and used to establish a potential energy surface, which compares well with previous theoretical studies. Utilizing the theoretical results, the kinetic energy dependent branching ratio between the ThO{sup +} + D{sub 2} and DThO{sup +} + D channels was calculated using both RRKM and phase space theory (PST). The results indicate that consideration of angular momentum conservation (as in PST) and spin-orbit corrected energies are needed to reproduce experimental results quantitatively. The PST modeling also provides relative energies for the two competing transition states that lead to the primary products, for which theory provides reasonable agreement. .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Kinetic electron and ion instability of the lunar wake simulated at physical mass ratio

The solar wind wake behind the moon is studied with 1D electrostatic particle-in-cell (PIC) simulations using a physical ion to electron mass ratio (unlike prior investigations); the simulations also apply more generally to supersonic flow of dense magnetized plasma past non-magnetic objects. In this work, a hybrid electrostatic Boltzmann electron treatment is first used to investigate the ion stability in the absence of kinetic electron effects, showing that the ions are two-stream unstable for downstream wake distances (in lunar radii) greater than about three times the solar wind Mach number. Simulations with PIC electrons are then used to show that kinetic electron effects can lead to disruption of the ion beams at least three times closer to the moon than in the hybrid simulations. This disruption occurs as the result of a novel wake phenomenon: the non-linear growth of electron holes spawned from a narrow dimple in the electron velocity distribution. Most of the holes arising from the dimple are small and quickly leave the wake, approximately following the unperturbed electron phase-space trajectories, but some holes originating near the center of the wake remain and grow large enough to trigger disruption of the ion beams. Lastly, non-linear kinetic-electron effects are therefore essential to a comprehensive understanding of the 1D electrostatic stability of such wakes, and possible observational signatures in ARTEMIS data from the lunar wake are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Communication: Wigner functions in action-angle variables, Bohr-Sommerfeld quantization, the Heisenberg correspondence principle, and a symmetrical quasi-classical approach to the full electronic density matrix

It is pointed out that the classical phase space distribution in action-angle (a-a) variables obtained from a Wigner function depends on how the calculation is carried out: if one computes the standard Wigner function in Cartesian variables (p, x), and then replaces p and x by their expressions in terms of a-a variables, one obtains a different result than if the Wigner function is computed directly in terms of the a-a variables. Furthermore, the latter procedure gives a result more consistent with classical and semiclassical theory—e.g., by incorporating the Bohr-Sommerfeld quantization condition (quantum states defined by integer values of the action variable) as well as the Heisenberg correspondence principle for matrix elements of an operator between such states—and has also been shown to be more accurate when applied to electronically non-adiabatic applications as implemented within the recently developed symmetrical quasi-classical (SQC) Meyer-Miller (MM) approach. Moreover, use of the Wigner function (obtained directly) in a-a variables shows how our standard SQC/MM approach can be used to obtain off-diagonal elements of the electronic density matrix by processing in a different way the same set of trajectories already used (in the SQC/MM methodology) to obtain the diagonal elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Suppression of Space Charge Induced Beam Halo in Nonlinear Focusing Channel

An intense non-uniform particle beam exhibits strong emittance growth and halo formation in focusing channels due to nonlinear space charge forces of the beam. This phenomenon limits beam brightness and results in particle losses. The problem is connected with irreversible distortion of phase space volume of the beam in conventional focusing structures due to filamentation in phase space. Emittance growth is accompanied by halo formation in real space, which results in inevitable particle losses. We discuss a new approach for solving a self-consistent problem for a matched non-uniform beam in two-dimensional geometry. The resulting solution is applied to the problem of beam transport, while avoiding emittance growth and halo formation by the use of nonlinear focusing field. Conservation of a beam distribution function is demonstrated analytically and by particle-in-cell simulation for a beam with a realistic beam distribution.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

THE CONVERGENCE OF PARTICLE-IN-CELL SCHEMES FOR COSMOLOGICAL DARK MATTER SIMULATIONS

Particle methods are a ubiquitous tool for solving the Vlasov–Poisson equation in comoving coordinates, which is used to model the gravitational evolution of dark matter (DM) in an expanding universe. However, these methods are known to produce poor results on idealized test problems, particularly at late times, after the particle trajectories have crossed. To investigate this, we have performed a series of one- and two-dimensional “Zel’dovich pancake” calculations using the popular particle-in-cell (PIC) method. We find that PIC can indeed converge on these problems provided that the following modifications are made. The first modification is to regularize the singular initial distribution function by introducing a small but finite artificial velocity dispersion. This process is analogous to artificial viscosity in compressible gas dynamics, and, as with artificial viscosity, the amount of regularization can be tailored so that its effect outside of a well-defined region—in this case, the high-density caustics—is small. The second modification is the introduction of a particle remapping procedure that periodically reexpresses the DM distribution function using a new set of particles. We describe a remapping algorithm that is third-order accurate and adaptive in phase space. This procedure prevents the accumulation of numerical errors in integrating the particle trajectories from growing large enough to significantly degrade the solution. Once both of these changes are made, PIC converges at second order on the Zel’dovich pancake problem, even at late times, after many caustics have formed. Furthermore, the resulting scheme does not suffer from the unphysical, small-scale “clumping” phenomenon known to occur on the pancake problem when the perturbation wavevector is not aligned with one of the Cartesian coordinate axes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Issues in measure-preserving three dimensional flow integrators: Self-adjointness, reversibility, and non-uniform time stepping

Properties of integration schemes for solenoidal fields in three dimensions are studied, with a focus on integrating magnetic field lines in a plasma using adaptive time stepping. It is shown that implicit midpoint (IM) and a scheme we call three-dimensional leapfrog (LF) can do a good job (in the sense of preserving KAM tori) of integrating fields that are reversible, or (for LF) have a “special divergence-free” (SDF) property. We review the notion of a self-adjoint scheme, showing that such schemes are at least second order accurate and can always be formed by composing an arbitrary scheme with its adjoint. We also review the concept of reversibility, showing that a reversible but not exactly volume-preserving scheme can lead to a fractal invariant measure in a chaotic region, although this property may not often be observable. We also show numerical results indicating that the IM and LF schemes can fail to preserve KAM tori when the reversibility property (and the SDF property for LF) of the field is broken. We discuss extensions to measure preserving flows, the integration of magnetic field lines in a plasma and the integration of rays for several plasma waves. The main new result of this paper relates to non-uniform time stepping for volume-preserving flows. We investigate two potential schemes, both based on the general method of Feng and Shang [Numer. Math. 71, 451 (1995)], in which the flow is integrated in split time steps, each Hamiltonian in two dimensions. The first scheme is an extension of the method of extended phase space, a well-proven method of symplectic integration with non-uniform time steps. This method is found not to work, and an explanation is given. The second method investigated is a method based on transformation to canonical variables for the two split-step Hamiltonian systems. This method, which is related to the method of non-canonical generating functions of Richardson and Finn [Plasma Phys. Controlled Fusion 54, 014004 (2012)], appears to work very well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗