Simulation: Plasma

1. Abstract This report details the design, implementation, and results of a computational simulation of a collisionless plasma. The Particle-in-Cell (PIC) method is employed to model the self-consistent motion of charged particles under the influence of internal and external electromagnetic fields. The simulation resolves fundamental plasma phenomena, including Debye shielding, plasma oscillations, and the onset of the two-stream instability. The code solves the coupled Vlasov-Maxwell system of equations on a discrete grid while tracking macro-particles. Results demonstrate the characteristic frequency of electron plasma oscillations (( \omega_pe )) and verify the dispersion relation for electrostatic waves. This framework provides a foundation for studying more complex kinetic effects in fusion, space, and laboratory plasmas. 2. Introduction Plasmas, ionized gases exhibiting collective behavior, are ubiquitous—from astrophysical jets to fusion reactors. Due to the non-linear and kinetic nature of plasma phenomena, analytical solutions are limited to highly idealized cases. Numerical simulation is therefore essential.

The PIC algorithm scales as ( O(N_p \cdot N_steps) ) where ( N_p ) is particle count. For this 1D case, runtime on a single CPU was under 10 minutes. A 2D or 3D simulation would require parallelization (MPI/OpenMP) and advanced field solvers. plasma simulation