Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

''Hybrid'' Monte Carlo simulation of ripple transport in stellarators

Technical Report ·
OSTI ID:5223712
A rapid simulation has been developed to accurately describe ripple effects on particle motion in stellarators at low collision frequencies. The majority of ripple trapped particles are followed through the iterative conservation of their bounce action J, while those particles in the region of phase space on either side of the ripple trap/detrap boundary are followed using guiding center equations. This ''hybrid'' formulation provides the most accurate numerical description of ripple trapped particle orbits to date, short of a fully guiding center code. This is important since these orbits are often substantially different from those allowed for in analytic theory and in other transport codes. Further, the hybrid simulation is much faster than a fully guiding center treatment, making all collision frequency regimes of interest accessible at reasonable costs of computer time. The methods employed allow the examination of stellarator configurations for which the magnitude of the toroidal well, epsilon/sub t/, is larger than that of the helical well, epsilon/sub h/, as well as configurations more often treated in which epsilon/sub h/ > epsilon/sub t/. Results are obtained for the usual analytic model of the helical ripple, epsilon/sub h/ = epsilon/sub h/(r), as well as a more realistic model for which epsilon/sub h/ = epsilon/sub h/(r,theta). The results are largely explainable in terms of existing analytic theories, although some slight modifications seem to be necessary.
Research Organization:
Wisconsin Univ., Madison (USA). Torsatron/Stellarator Lab.
DOE Contract Number:
AC02-84ER53166
OSTI ID:
5223712
Report Number(s):
DOE/ER/53166-T8; ON: DE86015582
Country of Publication:
United States
Language:
English