Accession Number:

ADA357700

Title:

A Study of Stability and Energy Conservation of a 3-D Electromagnetic PIC Code for Non-Orthogonal Meshes

Descriptive Note:

Final rept. 1 Jan 94-31 Dec 96

Corporate Author:

CALIFORNIA INST OF TECH PASADENA

Report Date:

1997-01-01

Pagination or Media Count:

15.0

Abstract:

A 3-D electromagnetic particle-in cell code EMPIC was developed for large-scale plasma simulation on parallel computers. To simulate plasma problems with complex geometries such as high-power microwave generation devices, curvilinear coordinates were used. A logically connected Cartesian grid consisted of hexahedral cells that could be deformed to body-fit complex shapes. A finite-volume method for a non-orthogonal grid was used to calculate the electromagnetic field. This method is based on Gedney-Lansing 1 and Madsen 2 algorithms, and is reduced to a standard FDTD algorithm for an orthogonal grid. Particle updates use current deposit formulation of Villasenor and Buneman and 3 generalized to non-orthogonal meshes, preserving charge and current. The numerical stability of the electromagnetic algorithm was analyzed for a planar EM wave propagation in a distorted periodic box. The influence of face-to edge transformation properties on the EM algorithm stability is analyzed. Energy conservation for a full PIC simulation of two-stream instability on non-orthogonal meshes is studied.

Subject Categories:

  • Electricity and Magnetism
  • Plasma Physics and Magnetohydrodynamics

Distribution Statement:

APPROVED FOR PUBLIC RELEASE