Electrostatic Simulation - Discontinuous Galerkin - FEM - Interior Penalty - Matlab/Python/C++
From left to right: Electric Potential, Electric Field, Electric Charge Density
Charged particles in electric field , two way coupled, elastic collision
Lagrange stuff:
Simplectic integrator, 1-form conservative
Particles with same mass and same absolute charge
Shape function: Gaussian function
Euler Stuff:
Cartesian mesh with Nx x Ny = 15 x 15
Interior Penalty discontinuous Galerkin Method (Flux + Primal form)
Matrix solver: Direct LU-Decomposition
Gauss-Legrende Distribution
Polynomial Degree=8
Order of Convergence=9