Computer simulation using particles
Computer simulation using particles
Numerical computation of internal & external flows: fundamentals of numerical discretization
Numerical computation of internal & external flows: fundamentals of numerical discretization
A vectorized particle tracer for unstructured grids
Journal of Computational Physics
Simultaneous potential and circuit solution for 1D bounded plasma particle simulation codes
Journal of Computational Physics
Dynamic load balancing for a 2D concurrent plasma PIC code
Journal of Computational Physics
A finite element formulation of the Darwin PIC model for use on unstructured grids
Journal of Computational Physics
Simultaneous potential and circuit solution for two-dimensional bounded plasma simulation codes
Journal of Computational Physics
Journal of Computational Physics
Plasma Physics Via Computer
Iterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems
High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids
Journal of Computational Physics
Algorithms for accurate collection, ejection, and loading in particle simulations
Journal of Computational Physics
Journal of Computational Physics
Journal of Computational Physics
Hi-index | 31.46 |
The mathematical formulation and computational implementation of a three-dimensional particle-in-cell methodology on unstructured Delaunay-Voronoi tetrahedral grids is presented. The method allows simulation of plasmas in complex domains and incorporates the duality of the Delaunay-Voronoi in all aspects of the particle-in-cell cycle. Charge assignment and field interpolation weighting schemes of zero- and first-order are formulated based on the theory of long-range constraints. Electric potential and fields are derived from a finite-volume formulation of Gauss' law using the Voronoi-Delaunay dual. Boundary conditions and the algorithms for injection, particle loading, particle motion, and particle tracking are implemented for unstructured Delaunay grids. Error and sensitivity analysis examines the effects of particles/cell, grid scaling, and timestep on the numerical heating, the slowing-down time, and the deflection times. The problem of current collection by cylindrical Langmuir probes in collisionless plasmas is used for validation. Numerical results compare favorably with previous numerical and analytical solutions for a wide range of probe radius to Debye length ratios, probe potentials, and electron to ion temperature ratios. The versatility of the methodology is demonstrated with the simulation of a complex plasma microsensor, a directional micro-retarding potential analyzer that includes a low transparency micro-grid.