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...
Gespeichert in:
Veröffentlicht in: | arXiv.org 2016-02 |
---|---|
Hauptverfasser: | , , |
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 & 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 |