Highorder nodal discontinuous Galerkin particleincell method on unstructured grids
Abstract
We present a highorder particleincell (PIC) algorithm for the simulation of kinetic plasmas dynamics. The core of the algorithm utilizes an unstructured grid discontinuous Galerkin Maxwell field solver combining highorder accuracy with geometric flexibility. We introduce algorithms in the Lagrangian framework that preserve the favorable properties of the field solver in the PIC solver. Fast fullorder interpolation and effective search algorithms are used for tracking individual particles on the general grid and smooth particle shape functions are introduced to ensure low noise in the charge and current density. A precomputed levelset distance function is employed to represent the geometry and facilitates complex particleboundary interaction. To enforce charge conservation we consider two different techniques, one based on projection and one on hyperbolic cleaning. Both are found to work well, although the latter is found be too expensive when used with explicit time integration. Examples of simple plasma phenomena, e.g., plasma waves, instabilities, and Landau damping are shown to agree well with theoretical predictions and/or results found by other computational methods. We also discuss generic well known problems such as numerical Cherenkov radiation and grid heating before presenting a few twodimensional tests, showing the potential of the current method to handle fullymore »
 Authors:
 Division of Applied Mathematics, Brown University, Box F, Providence, RI 02912 (United States). Email: gjacobs2@dam.brown.edu
 Division of Applied Mathematics, Brown University, Box F, Providence, RI 02912 (United States). Email: Jan.Hesthaven@brown.edu
 Publication Date:
 OSTI Identifier:
 20767040
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 214; Journal Issue: 1; Other Information: DOI: 10.1016/j.jcp.2005.09.008; PII: S00219991(05)004250; Copyright (c) 2005 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; ALGORITHMS; CHARGE CONSERVATION; CHERENKOV RADIATION; CURRENT DENSITY; FORECASTING; GEOMETRY; INSTABILITY; INTERPOLATION; LAGRANGIAN FUNCTION; LANDAU DAMPING; MAXWELL EQUATIONS; NOISE; PLASMA WAVES; RELATIVISTIC PLASMA; SIMULATION; TWODIMENSIONAL CALCULATIONS
Citation Formats
Jacobs, G.B., and Hesthaven, J.S.. Highorder nodal discontinuous Galerkin particleincell method on unstructured grids. United States: N. p., 2006.
Web. doi:10.1016/j.jcp.2005.09.008.
Jacobs, G.B., & Hesthaven, J.S.. Highorder nodal discontinuous Galerkin particleincell method on unstructured grids. United States. doi:10.1016/j.jcp.2005.09.008.
Jacobs, G.B., and Hesthaven, J.S.. Mon .
"Highorder nodal discontinuous Galerkin particleincell method on unstructured grids". United States.
doi:10.1016/j.jcp.2005.09.008.
@article{osti_20767040,
title = {Highorder nodal discontinuous Galerkin particleincell method on unstructured grids},
author = {Jacobs, G.B. and Hesthaven, J.S.},
abstractNote = {We present a highorder particleincell (PIC) algorithm for the simulation of kinetic plasmas dynamics. The core of the algorithm utilizes an unstructured grid discontinuous Galerkin Maxwell field solver combining highorder accuracy with geometric flexibility. We introduce algorithms in the Lagrangian framework that preserve the favorable properties of the field solver in the PIC solver. Fast fullorder interpolation and effective search algorithms are used for tracking individual particles on the general grid and smooth particle shape functions are introduced to ensure low noise in the charge and current density. A precomputed levelset distance function is employed to represent the geometry and facilitates complex particleboundary interaction. To enforce charge conservation we consider two different techniques, one based on projection and one on hyperbolic cleaning. Both are found to work well, although the latter is found be too expensive when used with explicit time integration. Examples of simple plasma phenomena, e.g., plasma waves, instabilities, and Landau damping are shown to agree well with theoretical predictions and/or results found by other computational methods. We also discuss generic well known problems such as numerical Cherenkov radiation and grid heating before presenting a few twodimensional tests, showing the potential of the current method to handle fully relativistic plasma dynamics in complex geometries.},
doi = {10.1016/j.jcp.2005.09.008},
journal = {Journal of Computational Physics},
number = 1,
volume = 214,
place = {United States},
year = {Mon May 01 00:00:00 EDT 2006},
month = {Mon May 01 00:00:00 EDT 2006}
}

An entropy stable nodal discontinuous Galerkin method for the two dimensional shallow water equations on unstructured curvilinear meshes with discontinuous bathymetry
We design an arbitrary highorder accurate nodal discontinuous Galerkin spectral element approximation for the nonlinear two dimensional shallow water equations with nonconstant, possibly discontinuous, bathymetry on unstructured, possibly curved, quadrilateral meshes. The scheme is derived from an equivalent flux differencing formulation of the split form of the equations. We prove that this discretization exactly preserves the local mass and momentum. Furthermore, combined with a special numerical interface flux function, the method exactly preserves the mathematical entropy, which is the total energy for the shallow water equations. By adding a specific form of interface dissipation to the baseline entropy conserving schememore » 
Third order maximumprinciplesatisfying direct discontinuous Galerkin methods for time dependent convection diffusion equations on unstructured triangular meshes
We develop 3rd order maximumprinciplesatisfying direct discontinuous Galerkin methods [8], [9], [19] and [21] for convection diffusion equations on unstructured triangular mesh. We carefully calculate the normal derivative numerical flux across element edges and prove that, with proper choice of parameter pair (β _{0},β _{1}) in the numerical flux formula, the quadratic polynomial solution satisfies strict maximum principle. The polynomial solution is bounded within the given range and third order accuracy is maintained. There is no geometric restriction on the meshes and obtuse triangles are allowed in the partition. As a result, a sequence of numerical examples are carried outmore »Cited by 5