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

Galerkin algorithm for multidimensional plasma simulation codes. Informal report

Technical Report ·
DOI:https://doi.org/10.2172/6159808· OSTI ID:6159808
A Galerkin finite element differencing scheme has been developed for a computer simulation of plasmas. The new difference equations identically satisfy an equation of continuity. Thus, the usual current correction procedure, involving inversion of Poisson's equation, is unnecessary. The algorithm is free of many numerical Cherenkov instabilities. This differencing scheme has been implemented in CCUBE, an already existing relativistic, electromagnetic, two-dimensional PIC code in arbitrary separable, orthogonal coordinates. The separability constraint is eliminated by the new algorithm. The new version of CCUBE exhibits good stability and accuracy with reduced computer memory and time requirements. Details of the algorithm and its implementation are presented.
Research Organization:
Los Alamos Scientific Lab., NM (USA)
DOE Contract Number:
W-7405-ENG-36
OSTI ID:
6159808
Report Number(s):
LA-7687-MS
Country of Publication:
United States
Language:
English