49,280 research outputs found
Accurate and efficient spin integration for particle accelerators
Accurate spin tracking is a valuable tool for understanding spin dynamics in
particle accelerators and can help improve the performance of an accelerator.
In this paper, we present a detailed discussion of the integrators in the spin
tracking code gpuSpinTrack. We have implemented orbital integrators based on
drift-kick, bend-kick, and matrix-kick splits. On top of the orbital
integrators, we have implemented various integrators for the spin motion. These
integrators use quaternions and Romberg quadratures to accelerate both the
computation and the convergence of spin rotations. We evaluate their
performance and accuracy in quantitative detail for individual elements as well
as for the entire RHIC lattice. We exploit the inherently data-parallel nature
of spin tracking to accelerate our algorithms on graphics processing units.Comment: 43 pages, 17 figure
New numerical approaches for modeling thermochemical convection in a compositionally stratified fluid
Seismic imaging of the mantle has revealed large and small scale
heterogeneities in the lower mantle; specifically structures known as large low
shear velocity provinces (LLSVP) below Africa and the South Pacific. Most
interpretations propose that the heterogeneities are compositional in nature,
differing in composition from the overlying mantle, an interpretation that
would be consistent with chemical geodynamic models. Numerical modeling of
persistent compositional interfaces presents challenges, even to
state-of-the-art numerical methodology. For example, some numerical algorithms
for advecting the compositional interface cannot maintain a sharp compositional
boundary as the fluid migrates and distorts with time dependent fingering due
to the numerical diffusion that has been added in order to maintain the upper
and lower bounds on the composition variable and the stability of the advection
method. In this work we present two new algorithms for maintaining a sharper
computational boundary than the advection methods that are currently openly
available to the computational mantle convection community; namely, a
Discontinuous Galerkin method with a Bound Preserving limiter and a
Volume-of-Fluid interface tracking algorithm. We compare these two new methods
with two approaches commonly used for modeling the advection of two distinct,
thermally driven, compositional fields in mantle convection problems; namely,
an approach based on a high-order accurate finite element method advection
algorithm that employs an artificial viscosity technique to maintain the upper
and lower bounds on the composition variable as well as the stability of the
advection algorithm and the advection of particles that carry a scalar quantity
representing the location of each compositional field. All four of these
algorithms are implemented in the open source FEM code ASPECT
An optically actuated surface scanning probe
We demonstrate the use of an extended, optically trapped probe that is capable of imaging surface topography with nanometre precision, whilst applying ultra-low, femto-Newton sized forces. This degree of precision and sensitivity is acquired through three distinct strategies. First, the probe itself is shaped in such a way as to soften the trap along the sensing axis and stiffen it in transverse directions. Next, these characteristics are enhanced by selectively position clamping independent motions of the probe. Finally, force clamping is used to refine the surface contact response. Detailed analyses are presented for each of these mechanisms. To test our sensor, we scan it laterally over a calibration sample consisting of a series of graduated steps, and demonstrate a height resolution of ∼ 11 nm. Using equipartition theory, we estimate that an average force of only ∼ 140 fN is exerted on the sample during the scan, making this technique ideal for the investigation of delicate biological samples
Kassiopeia: A Modern, Extensible C++ Particle Tracking Package
The Kassiopeia particle tracking framework is an object-oriented software
package using modern C++ techniques, written originally to meet the needs of
the KATRIN collaboration. Kassiopeia features a new algorithmic paradigm for
particle tracking simulations which targets experiments containing complex
geometries and electromagnetic fields, with high priority put on calculation
efficiency, customizability, extensibility, and ease of use for novice
programmers. To solve Kassiopeia's target physics problem the software is
capable of simulating particle trajectories governed by arbitrarily complex
differential equations of motion, continuous physics processes that may in part
be modeled as terms perturbing that equation of motion, stochastic processes
that occur in flight such as bulk scattering and decay, and stochastic surface
processes occuring at interfaces, including transmission and reflection
effects. This entire set of computations takes place against the backdrop of a
rich geometry package which serves a variety of roles, including initialization
of electromagnetic field simulations and the support of state-dependent
algorithm-swapping and behavioral changes as a particle's state evolves. Thanks
to the very general approach taken by Kassiopeia it can be used by other
experiments facing similar challenges when calculating particle trajectories in
electromagnetic fields. It is publicly available at
https://github.com/KATRIN-Experiment/Kassiopei
Development of a GPU-based Monte Carlo dose calculation code for coupled electron-photon transport
Monte Carlo simulation is the most accurate method for absorbed dose
calculations in radiotherapy. Its efficiency still requires improvement for
routine clinical applications, especially for online adaptive radiotherapy. In
this paper, we report our recent development on a GPU-based Monte Carlo dose
calculation code for coupled electron-photon transport. We have implemented the
Dose Planning Method (DPM) Monte Carlo dose calculation package (Sempau et al,
Phys. Med. Biol., 45(2000)2263-2291) on GPU architecture under CUDA platform.
The implementation has been tested with respect to the original sequential DPM
code on CPU in phantoms with water-lung-water or water-bone-water slab
geometry. A 20 MeV mono-energetic electron point source or a 6 MV photon point
source is used in our validation. The results demonstrate adequate accuracy of
our GPU implementation for both electron and photon beams in radiotherapy
energy range. Speed up factors of about 5.0 ~ 6.6 times have been observed,
using an NVIDIA Tesla C1060 GPU card against a 2.27GHz Intel Xeon CPU
processor.Comment: 13 pages, 3 figures, and 1 table. Paper revised. Figures update
A Parallel Histogram-based Particle Filter for Object Tracking on SIMD-based Smart Cameras
We present a parallel implementation of a histogram-based particle filter for object tracking on smart cameras based on SIMD processors. We specifically focus on parallel computation of the particle weights and parallel construction of the feature histograms since these are the major bottlenecks in standard implementations of histogram-based particle filters. The proposed algorithm can be applied with any histogram-based feature sets—we show in detail how the parallel particle filter can employ simple color histograms as well as more complex histograms of oriented gradients (HOG). The algorithm was successfully implemented on an SIMD processor and performs robust object tracking at up to 30 frames per second—a performance difficult to achieve even on a modern desktop computer
- …