A {delta}f Monte Carlo method to calculate plasma parameters

Journal Article · Wed Oct 01 04:00:00 EDT 1997 · Physics of Plasmas

DOI:https://doi.org/10.1063/1.872549· OSTI ID:542592

Sasinowski, M ^[1]; Boozer, A H ^[2]

---

1. Department of Physics, College of William and Mary, Williamsburg, Virginia 23185 (United States)
2. Department of Applied Physics, Columbia University, New York, New York 10027 (United States)

A {delta}f Monte Carlo code has been developed which efficiently calculates plasma currents by finding the deviation, {cflx f}, of the exact distribution function, f, from the Maxwellian, {bold F}, with {psi} the toroidal magnetic flux enclosed by a pressure surface and H the Hamiltonian. The particles in the simulation are followed with a traditional Monte Carlo scheme consisting of an orbit step in which new values for the positions and momenta are obtained and a collision step in which a Monte Carlo equivalent of the Lorentz operator is applied to change the pitch of each particle. The dual-species code includes momentum conservation and was used to study the bootstrap and two Fourier components of the Pfirsch{endash}Schl{umlt u}ter current in nonaxisymmetric magnetic field configuration. The effects of the toroidal ripple on the orbits of the particles as well as on the bootstrap current were examined. {copyright} {ital 1997 American Institute of Physics.}

• OSTI ID: 542592
• Journal Information: Physics of Plasmas, Journal Name: Physics of Plasmas Journal Issue: 10 Vol. 4; ISSN 1070-664X; ISSN PHPAEN
• Country of Publication: United States
• Language: English

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
ELECTRIC CURRENTS
MONTE CARLO METHOD
PLASMA
PLASMA SIMULATION