Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective

We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter-gat...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:arXiv.org 2016-02
Hauptverfasser: Moon, Haksu, Teixeira, Fernando L, Omelchenko, Yuri A
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
container_end_page
container_issue
container_start_page
container_title arXiv.org
container_volume
creator Moon, Haksu
Teixeira, Fernando L
Omelchenko, Yuri A
description We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter-gather algorithms presented recently by Campos-Pinto et al. [1] (to first-order) and by Squire et al. [2], but it is derived here in a streamlined fashion from a geometric viewpoint. Some ingredients reflecting this viewpoint are (1) the use of (discrete) differential forms of various degrees to represent fields, currents, and charged particles and provide localization rules for the degrees of freedom thereof on the various grid elements (nodes, edges, facets), (2) use of Whitney forms as basic interpolants from discrete differential forms to continuum space, and (3) use of a Galerkin formula for the discrete Hodge star operators (i.e., "mass matrices" incorporating the metric datum of the grid) applicable to generally irregular, unstructured grids. The expressions obtained for the scatter charges and scatter currents are very concise and do not involve numerical quadrature rules. Appropriate fractional areas within each grid element are identified that represent scatter charges and scatter currents within the element, and a simple geometric representation for the (exact) charge conservation mechanism is obtained by such identification. The field update is based on the coupled first-order Maxwell's curl equations to avoid spurious modes with secular growth (otherwise present in formulations that discretize the second-order wave equation). Examples are provided to verify preservation of discrete Gauss' law for all times.
doi_str_mv 10.48550/arxiv.1409.0854
format Article
fullrecord <record><control><sourceid>proquest_arxiv</sourceid><recordid>TN_cdi_arxiv_primary_1409_0854</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><sourcerecordid>2080304123</sourcerecordid><originalsourceid>FETCH-LOGICAL-a513-43578f2ba24652f398c84fcf7f93ac535ba1e97c45067abd23c0a25d0709191a3</originalsourceid><addsrcrecordid>eNotkEtrwkAUhYdCoWLdd1UGuh47T5N0J8GmBaGCdh2uk4mO5GFvJmJ3_emNtauz-c7h8BHyIPhUx8bwZ8CzP02F5smUx0bfkJFUSrBYS3lHJl134JzLWSSNUSPysziDDTTdA-4cS9umc3jyzY6uLYTgkGUQ9g7pvNq16MO-pmWLdAUYvK0c8w1LXVXRta_7CoIf-rRt6GfTBext6NEVNENfdC90TjPX1i6gt3TlsDs6G_zJ3ZPbEqrOTf5zTDavi036xpYf2Xs6XzIwQjGtTBSXcgtSz4wsVRLbWJe2jMpEgTXKbEG4JLLa8FkE20Iqy0Gagkc8EYkANSaP19k_O_kRfQ34nV8s5RdLA_B0BY7YfvWuC_mh7bEZLuWSx1xxLQaNvyH0a3Q</addsrcrecordid><sourcetype>Open Access Repository</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype><pqid>2080304123</pqid></control><display><type>article</type><title>Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective</title><source>arXiv.org</source><source>Free E- Journals</source><creator>Moon, Haksu ; Teixeira, Fernando L ; Omelchenko, Yuri A</creator><creatorcontrib>Moon, Haksu ; Teixeira, Fernando L ; Omelchenko, Yuri A</creatorcontrib><description>We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter-gather algorithms presented recently by Campos-Pinto et al. [1] (to first-order) and by Squire et al. [2], but it is derived here in a streamlined fashion from a geometric viewpoint. Some ingredients reflecting this viewpoint are (1) the use of (discrete) differential forms of various degrees to represent fields, currents, and charged particles and provide localization rules for the degrees of freedom thereof on the various grid elements (nodes, edges, facets), (2) use of Whitney forms as basic interpolants from discrete differential forms to continuum space, and (3) use of a Galerkin formula for the discrete Hodge star operators (i.e., "mass matrices" incorporating the metric datum of the grid) applicable to generally irregular, unstructured grids. The expressions obtained for the scatter charges and scatter currents are very concise and do not involve numerical quadrature rules. Appropriate fractional areas within each grid element are identified that represent scatter charges and scatter currents within the element, and a simple geometric representation for the (exact) charge conservation mechanism is obtained by such identification. The field update is based on the coupled first-order Maxwell's curl equations to avoid spurious modes with secular growth (otherwise present in formulations that discretize the second-order wave equation). Examples are provided to verify preservation of discrete Gauss' law for all times.</description><identifier>EISSN: 2331-8422</identifier><identifier>DOI: 10.48550/arxiv.1409.0854</identifier><language>eng</language><publisher>Ithaca: Cornell University Library, arXiv.org</publisher><subject>Algorithms ; Charged particles ; Computer simulation ; Conservation ; Datum (elevation) ; First principles ; Formulations ; Galerkin method ; Mathematical analysis ; Mathematics - Numerical Analysis ; Operators (mathematics) ; Particle in cell technique ; Post-processing ; Scattering ; Unstructured grids (mathematics) ; Wave equations</subject><ispartof>arXiv.org, 2016-02</ispartof><rights>2016. This work is published under http://arxiv.org/licenses/nonexclusive-distrib/1.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.</rights><rights>http://arxiv.org/licenses/nonexclusive-distrib/1.0</rights><oa>free_for_read</oa><woscitedreferencessubscribed>false</woscitedreferencessubscribed></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><link.rule.ids>228,230,776,780,881,27904</link.rule.ids><backlink>$$Uhttps://doi.org/10.1016/j.cpc.2015.04.014$$DView published paper (Access to full text may be restricted)$$Hfree_for_read</backlink><backlink>$$Uhttps://doi.org/10.48550/arXiv.1409.0854$$DView paper in arXiv$$Hfree_for_read</backlink></links><search><creatorcontrib>Moon, Haksu</creatorcontrib><creatorcontrib>Teixeira, Fernando L</creatorcontrib><creatorcontrib>Omelchenko, Yuri A</creatorcontrib><title>Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective</title><title>arXiv.org</title><description>We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter-gather algorithms presented recently by Campos-Pinto et al. [1] (to first-order) and by Squire et al. [2], but it is derived here in a streamlined fashion from a geometric viewpoint. Some ingredients reflecting this viewpoint are (1) the use of (discrete) differential forms of various degrees to represent fields, currents, and charged particles and provide localization rules for the degrees of freedom thereof on the various grid elements (nodes, edges, facets), (2) use of Whitney forms as basic interpolants from discrete differential forms to continuum space, and (3) use of a Galerkin formula for the discrete Hodge star operators (i.e., "mass matrices" incorporating the metric datum of the grid) applicable to generally irregular, unstructured grids. The expressions obtained for the scatter charges and scatter currents are very concise and do not involve numerical quadrature rules. Appropriate fractional areas within each grid element are identified that represent scatter charges and scatter currents within the element, and a simple geometric representation for the (exact) charge conservation mechanism is obtained by such identification. The field update is based on the coupled first-order Maxwell's curl equations to avoid spurious modes with secular growth (otherwise present in formulations that discretize the second-order wave equation). Examples are provided to verify preservation of discrete Gauss' law for all times.</description><subject>Algorithms</subject><subject>Charged particles</subject><subject>Computer simulation</subject><subject>Conservation</subject><subject>Datum (elevation)</subject><subject>First principles</subject><subject>Formulations</subject><subject>Galerkin method</subject><subject>Mathematical analysis</subject><subject>Mathematics - Numerical Analysis</subject><subject>Operators (mathematics)</subject><subject>Particle in cell technique</subject><subject>Post-processing</subject><subject>Scattering</subject><subject>Unstructured grids (mathematics)</subject><subject>Wave equations</subject><issn>2331-8422</issn><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>2016</creationdate><recordtype>article</recordtype><sourceid>BENPR</sourceid><sourceid>GOX</sourceid><recordid>eNotkEtrwkAUhYdCoWLdd1UGuh47T5N0J8GmBaGCdh2uk4mO5GFvJmJ3_emNtauz-c7h8BHyIPhUx8bwZ8CzP02F5smUx0bfkJFUSrBYS3lHJl134JzLWSSNUSPysziDDTTdA-4cS9umc3jyzY6uLYTgkGUQ9g7pvNq16MO-pmWLdAUYvK0c8w1LXVXRta_7CoIf-rRt6GfTBext6NEVNENfdC90TjPX1i6gt3TlsDs6G_zJ3ZPbEqrOTf5zTDavi036xpYf2Xs6XzIwQjGtTBSXcgtSz4wsVRLbWJe2jMpEgTXKbEG4JLLa8FkE20Iqy0Gagkc8EYkANSaP19k_O_kRfQ34nV8s5RdLA_B0BY7YfvWuC_mh7bEZLuWSx1xxLQaNvyH0a3Q</recordid><startdate>20160226</startdate><enddate>20160226</enddate><creator>Moon, Haksu</creator><creator>Teixeira, Fernando L</creator><creator>Omelchenko, Yuri A</creator><general>Cornell University Library, arXiv.org</general><scope>8FE</scope><scope>8FG</scope><scope>ABJCF</scope><scope>ABUWG</scope><scope>AFKRA</scope><scope>AZQEC</scope><scope>BENPR</scope><scope>BGLVJ</scope><scope>CCPQU</scope><scope>DWQXO</scope><scope>HCIFZ</scope><scope>L6V</scope><scope>M7S</scope><scope>PIMPY</scope><scope>PQEST</scope><scope>PQQKQ</scope><scope>PQUKI</scope><scope>PRINS</scope><scope>PTHSS</scope><scope>AKZ</scope><scope>GOX</scope></search><sort><creationdate>20160226</creationdate><title>Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective</title><author>Moon, Haksu ; Teixeira, Fernando L ; Omelchenko, Yuri A</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-a513-43578f2ba24652f398c84fcf7f93ac535ba1e97c45067abd23c0a25d0709191a3</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>2016</creationdate><topic>Algorithms</topic><topic>Charged particles</topic><topic>Computer simulation</topic><topic>Conservation</topic><topic>Datum (elevation)</topic><topic>First principles</topic><topic>Formulations</topic><topic>Galerkin method</topic><topic>Mathematical analysis</topic><topic>Mathematics - Numerical Analysis</topic><topic>Operators (mathematics)</topic><topic>Particle in cell technique</topic><topic>Post-processing</topic><topic>Scattering</topic><topic>Unstructured grids (mathematics)</topic><topic>Wave equations</topic><toplevel>online_resources</toplevel><creatorcontrib>Moon, Haksu</creatorcontrib><creatorcontrib>Teixeira, Fernando L</creatorcontrib><creatorcontrib>Omelchenko, Yuri A</creatorcontrib><collection>ProQuest SciTech Collection</collection><collection>ProQuest Technology Collection</collection><collection>Materials Science &amp; Engineering Collection</collection><collection>ProQuest Central (Alumni Edition)</collection><collection>ProQuest Central UK/Ireland</collection><collection>ProQuest Central Essentials</collection><collection>ProQuest Central</collection><collection>Technology Collection</collection><collection>ProQuest One Community College</collection><collection>ProQuest Central Korea</collection><collection>SciTech Premium Collection</collection><collection>ProQuest Engineering Collection</collection><collection>Engineering Database</collection><collection>Publicly Available Content Database</collection><collection>ProQuest One Academic Eastern Edition (DO NOT USE)</collection><collection>ProQuest One Academic</collection><collection>ProQuest One Academic UKI Edition</collection><collection>ProQuest Central China</collection><collection>Engineering Collection</collection><collection>arXiv Mathematics</collection><collection>arXiv.org</collection><jtitle>arXiv.org</jtitle></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext</fulltext></delivery><addata><au>Moon, Haksu</au><au>Teixeira, Fernando L</au><au>Omelchenko, Yuri A</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective</atitle><jtitle>arXiv.org</jtitle><date>2016-02-26</date><risdate>2016</risdate><eissn>2331-8422</eissn><abstract>We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter-gather algorithms presented recently by Campos-Pinto et al. [1] (to first-order) and by Squire et al. [2], but it is derived here in a streamlined fashion from a geometric viewpoint. Some ingredients reflecting this viewpoint are (1) the use of (discrete) differential forms of various degrees to represent fields, currents, and charged particles and provide localization rules for the degrees of freedom thereof on the various grid elements (nodes, edges, facets), (2) use of Whitney forms as basic interpolants from discrete differential forms to continuum space, and (3) use of a Galerkin formula for the discrete Hodge star operators (i.e., "mass matrices" incorporating the metric datum of the grid) applicable to generally irregular, unstructured grids. The expressions obtained for the scatter charges and scatter currents are very concise and do not involve numerical quadrature rules. Appropriate fractional areas within each grid element are identified that represent scatter charges and scatter currents within the element, and a simple geometric representation for the (exact) charge conservation mechanism is obtained by such identification. The field update is based on the coupled first-order Maxwell's curl equations to avoid spurious modes with secular growth (otherwise present in formulations that discretize the second-order wave equation). Examples are provided to verify preservation of discrete Gauss' law for all times.</abstract><cop>Ithaca</cop><pub>Cornell University Library, arXiv.org</pub><doi>10.48550/arxiv.1409.0854</doi><oa>free_for_read</oa></addata></record>
fulltext fulltext
identifier EISSN: 2331-8422
ispartof arXiv.org, 2016-02
issn 2331-8422
language eng
recordid cdi_arxiv_primary_1409_0854
source arXiv.org; Free E- Journals
subjects Algorithms
Charged particles
Computer simulation
Conservation
Datum (elevation)
First principles
Formulations
Galerkin method
Mathematical analysis
Mathematics - Numerical Analysis
Operators (mathematics)
Particle in cell technique
Post-processing
Scattering
Unstructured grids (mathematics)
Wave equations
title Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective
url https://sfx.bib-bvb.de/sfx_tum?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2025-01-27T06%3A51%3A50IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-proquest_arxiv&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=Exact%20Charge-Conserving%20Scatter-Gather%20Algorithm%20for%20Particle-in-Cell%20Simulations%20on%20Unstructured%20Grids:%20A%20Geometric%20Perspective&rft.jtitle=arXiv.org&rft.au=Moon,%20Haksu&rft.date=2016-02-26&rft.eissn=2331-8422&rft_id=info:doi/10.48550/arxiv.1409.0854&rft_dat=%3Cproquest_arxiv%3E2080304123%3C/proquest_arxiv%3E%3Curl%3E%3C/url%3E&disable_directlink=true&sfx.directlink=off&sfx.report_link=0&rft_id=info:oai/&rft_pqid=2080304123&rft_id=info:pmid/&rfr_iscdi=true