NDLI: Grid-Free Plasma Simulations

Content Provider	IEEE Xplore Digital Library
Author	Christlieb, A.J. Krasny, R. Boyd, I.D. Emhoff, J. Verboncoeur, J.
Copyright Year	2005
Description	Author affiliation: Michigan Univ., Ann Arbor, MI (Christlieb, A.J.; Krasny, R.; Boyd, I.D.; Emhoff, J.)
Abstract	Summary form only given. Many problems in plasma physics are modeled using a Lagrangian approach. One of the most common models is particle-in-cell (PIC), where the system is modeled as a collection of macro particles which interact through long range forces. To compute long range forces in fewer than $(N^{2})$ operations, particles are interpolated to a mesh where the field equations are solved. Afterwards, the fields are interpolated back to the macro particles. The major problems with this approach are that it can not resolve steep plasma gradients, has cumulative interpolation errors, and has difficulty describing non-conformal domains. The current work seeks to overcome these limitations in particle plasma models by developing a grid-free approach to plasma simulations. The boundary integral/treecode (BIT) method merges boundary integral formulations with treecode algorithms. Using Green's theorem, and setting the functions in Green's theorem to u=G and v=Phi, we can express Poisson's equation as an integral equation Phi(y) = $intint_{Omega}$ $rho(x)/epsi_{o}$ G(x/y) dOmega - $conint_{deltaOmega}$ (Phi(x) $grad_{x}$ G(x/y) - G(x/y) grad Phi(x)) nds where G(x/y) is the free space Green's function for the Laplace operator. If we are modeling the system as a collection of point charges, rho is given by rho = Sigmai-1 $^{N}$ $q_{i}delta(x-Zi).$ The contribution to the field from the charge density reduces to $intint_{Omega}$ rho(x)G(x/y) dOmega = $Sigma_{i-1}$ $^{N}$ $q_{i}G(Z_{i}/y).$ Choosing an appropriate collocation method for the boundary integral formulation effectively reformulates the boundary as a set of point charges, making the field evaluation mesh-free and capable of handling complicated domains. The operation count in the field evaluation is reduced from $O(N_{2})$ to O(N log N) by using a treecode algorithm, which gains its efficiency in essence by treating clusters of particles at a distance as point charges at the center of the cluster. BIT has been applied to several bounded plasmas, including 1D sheath formation in DC discharges, 1D formation of a virtual cathode, 2D planar and cylindrical ion optics and the Penning-Malmberg trap. The method has been compared with PIC for the first three examples. For static systems, such as the ion optics and DC sheath, the results of the approaches are in good agreement. In the dynamic example of a 1D virtual cathode, obvious differences between the PIC and BIT are observed. The differences are attributed to the fact that PIC does not resolve interparticle forces within a mesh cell. Currently, the method is being used to investigate a potentially new instability in Penning-Malmberg traps and is being extended to grid-free hybrid fluid-kinetic plasma simulations
Sponsorship	Plasma Sci. Appl. Comm. IEEE Nucl. Plasma Sci. Soc.
Starting Page	246
Ending Page	246
File Size	1358777
Page Count	1
File Format	PDF
ISBN	0780393007
ISSN	07309244
DOI	10.1109/PLASMA.2005.359318
Language	English
Publisher	Institute of Electrical and Electronics Engineers, Inc. (IEEE)
Publisher Date	2005-06-20
Publisher Place	USA
Access Restriction	Subscribed
Rights Holder	Institute of Electrical and Electronics Engineers, Inc. (IEEE)
Subject Keyword	Plasma simulation Integral equations Cathodes Particle beam optics Physics Lagrangian functions Interpolation Poisson equations Green's function methods Laplace equations
Content Type	Text
Resource Type	Article
Subject	Atomic and Molecular Physics, and Optics Condensed Matter Physics Electrical and Electronic Engineering

Sl.	Authority	Responsibilities	Communication Details
1	Ministry of Education (GoI), Department of Higher Education	Sanctioning Authority	https://www.education.gov.in/ict-initiatives
2	Indian Institute of Technology Kharagpur	Host Institute of the Project: The host institute of the project is responsible for providing infrastructure support and hosting the project	https://www.iitkgp.ac.in
3	National Digital Library of India Office, Indian Institute of Technology Kharagpur	The administrative and infrastructural headquarters of the project	Dr. B. Sutradhar bsutra@ndl.gov.in
4	Project PI / Joint PI	Principal Investigator and Joint Principal Investigators of the project	Dr. B. Sutradhar bsutra@ndl.gov.in Prof. Saswat Chakrabarti will be added soon
5	Website/Portal (Helpdesk)	Queries regarding NDLI and its services	support@ndl.gov.in
6	Contents and Copyright Issues	Queries related to content curation and copyright issues	content@ndl.gov.in
7	National Digital Library of India Club (NDLI Club)	Queries related to NDLI Club formation, support, user awareness program, seminar/symposium, collaboration, social media, promotion, and outreach	clubsupport@ndl.gov.in
8	Digital Preservation Centre (DPC)	Assistance with digitizing and archiving copyright-free printed books	dpc@ndl.gov.in
9	IDR Setup or Support	Queries related to establishment and support of Institutional Digital Repository (IDR) and IDR workshops	idr@ndl.gov.in