A P-Adaptive Hybridizable Discontinuous Galerkin Spectral Element Method for Electrostatic Particle-in-Cell Simulations

Marcel Pfeiffer; Stephen Copplestone; Tobias Ott

arxiv: 2604.00782 · v2 · submitted 2026-04-01 · ⚛️ physics.comp-ph

A P-Adaptive Hybridizable Discontinuous Galerkin Spectral Element Method for Electrostatic Particle-in-Cell Simulations

Tobias Ott , Stephen Copplestone , Marcel Pfeiffer This is my paper

Pith reviewed 2026-05-13 22:08 UTC · model grok-4.3

classification ⚛️ physics.comp-ph

keywords p-adaptivehybridizable discontinuous Galerkinspectral element methodPoisson equationparticle-in-cellelectrostatic plasmadegree of freedom reductionion optic

0 comments

The pith

P-adaptive HDG spectral elements reduce global degrees of freedom for electrostatic Poisson solves in particle-in-cell plasma models.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a p-adaptive hybridizable discontinuous Galerkin spectral element method for the Poisson equation inside electrostatic particle-in-cell simulations. Polynomial degree is raised only inside elements that contain strong gradients, so computational effort concentrates where the solution varies most. This produces fewer total degrees of freedom than a uniform high-order mesh while the solution accuracy stays comparable on standard test problems. The approach is coded inside the open-source PICLas framework and demonstrated on a dielectric sphere, a one-dimensional plasma sheath, and a two-dimensional axisymmetric ion-optic geometry.

Core claim

Element-local p-refinement inside the HDG-SEM discretization of the Poisson equation allows the global number of degrees of freedom to be reduced relative to fixed high-order schemes while solution accuracy on benchmark electrostatic PIC problems is preserved.

What carries the argument

Element-wise p-adaptation driven by local gradient indicators inside the hybridizable discontinuous Galerkin spectral element Poisson solver.

If this is right

The total number of degrees of freedom drops compared with uniform high-order meshes while accuracy on the dielectric-sphere and one-dimensional sheath benchmarks remains comparable.
The method can be applied directly to two-dimensional axisymmetric ion-optic geometries inside existing PIC codes.
Computational savings scale with the fraction of the domain that requires high polynomial order.
The implementation is already available inside the open-source PICLas framework.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same local adaptation logic could be tested on three-dimensional multi-scale discharges where gradients occupy only a small fraction of the volume.
Coupling the p-adaptive field solver with dynamic particle weighting or mesh motion might further reduce cost in time-dependent problems.
Interface flux continuity across polynomial-degree jumps will need explicit verification in regimes with strong particle-induced source terms.

Load-bearing premise

Gradient-based local p-refinement keeps solution error controlled and prevents artifacts at interfaces between elements of different polynomial degree.

What would settle it

A simulation of a plasma sheath or ion optic in which the error or charge conservation error grows measurably at p-refinement interfaces compared with a uniform high-order run would falsify the claim.

Figures

Figures reproduced from arXiv: 2604.00782 by Marcel Pfeiffer, Stephen Copplestone, Tobias Ott.

**Figure 2.** Figure 2: Refinement levels and P-adaptation area used for the dielectric sphere convergence tests. The [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗

**Figure 3.** Figure 3: Convergence of the normalized L 2 error of the electric potential for the dielectric sphere test case with different polynomial degrees. The results with (solid) and without (dashed) p-adaptation are compared at different polynomial degrees. In the P-adaptive case, the polynomial degree for the elements inside the dielectric sphere and not adjacent to its boundary is set to 2 [PITH_FULL_IMAGE:figures/full… view at source ↗

**Figure 4.** Figure 4: Electric potential for two plasma sheath simulations on a mesh with four elements. In the first [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗

**Figure 5.** Figure 5: Discretized domain and boundaries used in the simulation. The values of the boundary conditions [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗

**Figure 6.** Figure 6: Polynomial degree distribution across the mesh after each successive simulation run. The poly [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗

**Figure 7.** Figure 7: Electric potential (upper plot) and Xe+ charge density (lower plot) as simulated with the padapted HDG-SEM. upstream region where the ion velocities are low. 6. Conclusion This paper presented a p-adaptive high-order hybridizable discontinuous Galerkin spectral element method (HDG-SEM) for solving the Poisson equation within electrostatic Particle-in-Cell plasma simulations. The method was successfully i… view at source ↗

read the original abstract

This paper presents a p-adaptive high-order hybridizable discontinuous Galerkin spectral element method (HDG-SEM) for solving the Poisson equation in electrostatic plasma simulations using particle-in-cell (PIC) schemes. This approach enables element-local refinement of the polynomial degree, concentrating computational effort specifically in regions with strong gradients. Thus, the method significantly reduces the global number of degrees of freedom compared to uniform high-order methods. The proposed method is implemented in the open-source framework PICLas and validated through a series of benchmark test cases, including a dielectric sphere and a one-dimensional plasma sheath. Finally, a two-dimensional axisymmetric simulation of an ion optic demonstrates the method's capability to efficiently model complex plasma phenomena but also highlights current limitations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The p-adaptive HDG-SEM for electrostatic PIC Poisson solves reduces degrees of freedom on the reported benchmarks but leaves the variable-order face coupling rules underspecified.

read the letter

The main thing here is a working p-adaptive version of hybridizable DG spectral elements applied to the Poisson solve inside an electrostatic PIC code. They implement it in the open-source PICLas framework and show that local polynomial degree increases driven by gradient indicators cut the total number of degrees of freedom relative to uniform high-order runs while keeping accuracy on a dielectric sphere and a 1-D sheath. The 2-D axisymmetric ion-optic example is presented as a realistic test, and the authors are upfront that it also reveals current limitations. That combination of adaptivity with HDG-SEM for this specific application is new enough to be worth noting. The code is public, which helps reproducibility. The benchmarks are standard and the efficiency claim is quantified in terms of global DOFs rather than just wall time, so the central result is at least falsifiable. The soft spot is exactly the one flagged in the stress test. HDG relies on a single-valued hybrid trace on element faces; when adjacent elements carry different polynomial degrees, the face space must be chosen consistently or projected. The abstract describes the adaptation as purely element-local and gradient-based, with no mention of a face-degree rule or an a-priori estimate that accounts for the resulting nonconformity. The sphere and 1-D sheath cases have relatively smooth or aligned interfaces, so they may not expose artifacts that would appear in a multi-scale sheath or optic problem with strong local refinement. If the full manuscript supplies a clear mortar or min/max rule plus numerical checks on interface error, this concern shrinks; on the supplied description it remains open. This paper is aimed at computational plasma physicists who already run high-order PIC and want to reduce cost in regions of strong gradients. A reader working on adaptive high-order methods for Poisson or Maxwell problems will get concrete implementation details and timing numbers. It is coherent on its own terms and shows honest engagement with the practical limits of the method, so it deserves a serious referee rather than a desk reject. I would bring it to a reading group for the implementation choices and the ion-optic results, but I would not cite it until the interface consistency question is closed.

Referee Report

2 major / 2 minor

Summary. The paper presents a p-adaptive high-order hybridizable discontinuous Galerkin spectral element method (HDG-SEM) for the Poisson equation in electrostatic PIC simulations. Element-local polynomial-degree refinement is driven by gradient indicators to concentrate degrees of freedom in high-gradient regions, thereby reducing the global DOF count relative to uniform high-order discretizations. The method is implemented in PICLas and validated on a dielectric sphere, a one-dimensional plasma sheath, and a two-dimensional axisymmetric ion-optic problem.

Significance. If the interface consistency of the hybrid variable is rigorously established, the approach would offer a practical route to efficient high-order electrostatic PIC simulations in multi-scale plasma configurations, directly addressing the computational cost of uniform high-order methods while preserving accuracy on standard benchmarks.

major comments (2)

[§3 and §2.2] §3 (p-adaptation strategy) and §2.2 (HDG formulation): the element-local gradient-based criterion for increasing polynomial degree does not specify how the polynomial space for the hybrid trace variable is chosen on faces shared by elements of differing degree. In HDG the numerical flux couples adjacent elements through this trace; without an explicit rule (minimum-p, maximum-p, or mortar projection) or an a priori error estimate that accounts for the resulting nonconformity, the claim that local adaptation preserves global accuracy without interface artifacts remains unproven.
[Table 2] Table 2 (error tables for dielectric sphere and 1-D sheath): quantitative comparison of total DOFs and L2 errors versus uniform p-refinement is presented, but the tables do not report the number or location of p-jumps at interfaces nor any measure of trace-variable discontinuity; this information is required to assess whether the observed error reduction is achieved without hidden interface penalties.

minor comments (2)

[Figure 4] Figure 4 (ion-optic mesh): the color scale for local polynomial degree is not labeled with the exact range of p values used; adding this would improve reproducibility.
[Eq. (12)] Equation (12) (adaptation indicator): the threshold value for gradient-based refinement is introduced without a sensitivity study; a brief remark on its selection would clarify robustness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript describing the p-adaptive HDG-SEM for electrostatic PIC simulations. We address each major comment below and will revise the manuscript accordingly to improve clarity and completeness.

read point-by-point responses

Referee: [§3 and §2.2] §3 (p-adaptation strategy) and §2.2 (HDG formulation): the element-local gradient-based criterion for increasing polynomial degree does not specify how the polynomial space for the hybrid trace variable is chosen on faces shared by elements of differing degree. In HDG the numerical flux couples adjacent elements through this trace; without an explicit rule (minimum-p, maximum-p, or mortar projection) or an a priori error estimate that accounts for the resulting nonconformity, the claim that local adaptation preserves global accuracy without interface artifacts remains unproven.

Authors: We agree that an explicit rule for the hybrid trace variable on p-nonconforming faces is essential. In the implementation, the polynomial degree for the trace variable on each shared face is set to the minimum of the degrees of the two adjacent elements. This choice maintains flux consistency and stability without requiring mortar projections. We will add a precise description of this rule to Section 2.2, together with a short justification referencing existing HDG theory for variable-degree spaces. While a full a priori error estimate for the p-adaptive case lies beyond the scope of the present work, the numerical benchmarks exhibit no visible interface artifacts, and we will expand the discussion in Section 3 to note this empirical evidence and cite supporting literature on nonconforming HDG methods. revision: yes
Referee: [Table 2] Table 2 (error tables for dielectric sphere and 1-D sheath): quantitative comparison of total DOFs and L2 errors versus uniform p-refinement is presented, but the tables do not report the number or location of p-jumps at interfaces nor any measure of trace-variable discontinuity; this information is required to assess whether the observed error reduction is achieved without hidden interface penalties.

Authors: We acknowledge that additional quantitative information on p-jumps would allow readers to better evaluate interface effects. In the revised manuscript we will augment Table 2 with the number of p-jumps present in each adapted configuration. We will also add a supplementary note (or short table) reporting the maximum L2 discontinuity of the hybrid trace variable across all interfaces for the p-adaptive runs. These data are directly available from the simulations and will demonstrate that the observed error reductions occur with only small, controlled discontinuities. revision: yes

Circularity Check

0 steps flagged

P-adaptive HDG-SEM derivation is self-contained; no load-bearing reductions to fitted inputs or self-citations

full rationale

The paper introduces an element-local gradient-based p-adaptation criterion within an HDG-SEM discretization of the Poisson equation for PIC electrostatics. The central claim (reduced global DOFs while preserving accuracy) is demonstrated via direct numerical validation on independent benchmarks (dielectric sphere, 1-D sheath, 2-D ion optic). No equations equate a derived quantity to its own fitting procedure, no uniqueness theorem is imported from prior self-work to force the formulation, and the adaptation rule is presented as a practical heuristic rather than a first-principles derivation that collapses by construction. Interface consistency for the hybrid variable when neighboring elements differ in p is addressed as an implementation detail, not as a tautological assumption. This yields a minor self-citation score only; the method remains externally falsifiable through the reported benchmark comparisons.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Numerical method paper; no new physical axioms or invented entities. Free parameters such as gradient thresholds for p-refinement are expected but not quantified in the abstract.

pith-pipeline@v0.9.0 · 5423 in / 1003 out tokens · 42958 ms · 2026-05-13T22:08:34.392192+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

iterative modal analysis strategy... squared normalized modal contribution falls below a predefined threshold ε... dynamic threshold that accounts for the noise level... ˆρ²_k < εσ²_k

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages

[1]

Adams, Steven Benson, Jed Brown, Peter Brune, Kris Buschelman, Emil M

Satish Balay, Shrirang Abhyankar, Mark F. Adams, Steven Benson, Jed Brown, Peter Brune, Kris Buschelman, Emil M. Constantinescu, Lisandro Dalcin, Alp Dener, Victor 21 Eijkhout, Jacob Faibussowitsch, William D. Gropp, Václav Hapla, Tobin Isaac, Pierre Jolivet, Dmitry Karpeev, Dinesh Kaushik, Matthew G. Knepley, Fande Kong, Scott Kruger, Dave A. May, Lois C...

work page 2025
[2]

Raumfahrt

Tilman Binder.Development and Application of PICLas for Combined Optic-/Plume- Simulation of Ion-Propulsion Systems. Raumfahrt. Verlag Dr. Hut, München, 1. auflage edition, 2019

work page 2019
[3]

Birdsall and A

Charles K. Birdsall and A. Bruce Langdon.Plasma Physics via Computer Simulation. McGraw-Hill, New York, 1985

work page 1985
[4]

Courier Corporation, 2001

John P Boyd.Chebyshev and Fourier spectral methods. Courier Corporation, 2001

work page 2001
[5]

A superconvergent LDG- hybridizable Galerkin method for second-order elliptic problems.Mathematics of Com- putation, 77(264):1887–1916, May 2008

Bernardo Cockburn, Bo Dong, and Johnny Guzmán. A superconvergent LDG- hybridizable Galerkin method for second-order elliptic problems.Mathematics of Com- putation, 77(264):1887–1916, May 2008

work page 1916
[6]

A Hybridizable Discontinuous Galerkin Method for Steady-State Convection-Diffusion- Reaction Problems.SIAM Journal on Scientific Computing, 31(5):3827–3846, January 2009

Bernardo Cockburn, Bo Dong, Johnny Guzmán, Marco Restelli, and Riccardo Sacco. A Hybridizable Discontinuous Galerkin Method for Steady-State Convection-Diffusion- Reaction Problems.SIAM Journal on Scientific Computing, 31(5):3827–3846, January 2009

work page 2009
[7]

The Derivation of Hybridizable Dis- continuous Galerkin Methods for Stokes Flow.SIAM Journal on Numerical Analysis, 47(2):1092–1125, January 2009

Bernardo Cockburn and Jayadeep Gopalakrishnan. The Derivation of Hybridizable Dis- continuous Galerkin Methods for Stokes Flow.SIAM Journal on Numerical Analysis, 47(2):1092–1125, January 2009

work page 2009
[8]

Bernardo Cockburn, Jayadeep Gopalakrishnan, and Raytcho Lazarov. Unified Hy- bridization of Discontinuous Galerkin, Mixed, and Continuous Galerkin Methods for Second Order Elliptic Problems.SIAM Journal on Numerical Analysis, 47(2):1319– 1365, January 2009

work page 2009
[9]

Fasoulas, C.-D

S. Fasoulas, C.-D. Munz, M. Pfeiffer, J. Beyer, T. Binder, S. Copplestone, A. Mirza, P. Nizenkov, P. Ortwein, and W. Reschke. Combining particle-in-cell and direct sim- ulation Monte Carlo for the simulation of reactive plasma flows.Physics of Fluids, 31(7):072006, July 2019

work page 2019
[10]

Hairer, Christian Lubich, and Gerhard Wanner.Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations

E. Hairer, Christian Lubich, and Gerhard Wanner.Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations. Number 31 in Springer Series in Computational Mathematics. Springer, Berlin ; New York, 2nd ed. edition, 2006

work page 2006
[11]

CRC Press, 1st ed

R.W Hockney and J.W Eastwood.Computer Simulation Using Particles. CRC Press, 1st ed. edition, 1988. 22

work page 1988
[12]

HYPRE: High performance preconditioners.https://llnl.gov/casc/hypre,https: //github.com/hypre-space/hypre

work page
[13]

Wiley, New York, NY, 3rd ed

John David Jackson.Classical electrodynamics. Wiley, New York, NY, 3rd ed. edition, 1999

work page 1999
[14]

Jacobs and J.S

G.B. Jacobs and J.S. Hesthaven. High-order nodal discontinuous Galerkin particle-in- cell method on unstructured grids.Journal of Computational Physics, 214(1):96–121, May 2006

work page 2006
[15]

The Bohm criterion and sheath formation.Journal of Physics D: Applied Physics, 24(4):493–518, April 1991

K -U Riemann. The Bohm criterion and sheath formation.Journal of Physics D: Applied Physics, 24(4):493–518, April 1991

work page 1991
[16]

Blind, Anna Schwarz, Marius Kurz, Felix Rodach, Stephen M

Patrick Kopper, Marcel P. Blind, Anna Schwarz, Marius Kurz, Felix Rodach, Stephen M. Copplestone, and Andrea D. Beck. Pyhope: A python toolkit for three-dimensional unstructured high-order meshes.Journal of Open Source Software, 10(115):8769, 2025

work page 2025
[17]

Adaptive mesh strategies for the spectral element method.Com- puter Methods in Applied Mechanics and Engineering, 116(1-4):77–86, January 1994

Catherine Mavriplis. Adaptive mesh strategies for the spectral element method.Com- puter Methods in Applied Mechanics and Engineering, 116(1-4):77–86, January 1994

work page 1994
[18]

Coupled Particle-In-Cell and Direct Sim- ulation Monte Carlo method for simulating reactive plasma flows.Comptes Rendus

Claus-Dieter Munz, Monika Auweter-Kurtz, Stefanos Fasoulas, Asim Mirza, Philip Or- twein, Marcel Pfeiffer, and Torsten Stindl. Coupled Particle-In-Cell and Direct Sim- ulation Monte Carlo method for simulating reactive plasma flows.Comptes Rendus. Mécanique, 342(10-11):662–670, August 2014

work page 2014
[19]

Pfeiffer, F

M. Pfeiffer, F. Hindenlang, T. Binder, S.M. Copplestone, C.-D. Munz, and S. Fasoulas. A Particle-in-Cell solver based on a high-order hybridizable discontinuous Galerkin spectral element method on unstructured curved meshes.Computer Methods in Applied Mechanics and Engineering, 349:149–166, June 2019

work page 2019

[1] [1]

Adams, Steven Benson, Jed Brown, Peter Brune, Kris Buschelman, Emil M

Satish Balay, Shrirang Abhyankar, Mark F. Adams, Steven Benson, Jed Brown, Peter Brune, Kris Buschelman, Emil M. Constantinescu, Lisandro Dalcin, Alp Dener, Victor 21 Eijkhout, Jacob Faibussowitsch, William D. Gropp, Václav Hapla, Tobin Isaac, Pierre Jolivet, Dmitry Karpeev, Dinesh Kaushik, Matthew G. Knepley, Fande Kong, Scott Kruger, Dave A. May, Lois C...

work page 2025

[2] [2]

Raumfahrt

Tilman Binder.Development and Application of PICLas for Combined Optic-/Plume- Simulation of Ion-Propulsion Systems. Raumfahrt. Verlag Dr. Hut, München, 1. auflage edition, 2019

work page 2019

[3] [3]

Birdsall and A

Charles K. Birdsall and A. Bruce Langdon.Plasma Physics via Computer Simulation. McGraw-Hill, New York, 1985

work page 1985

[4] [4]

Courier Corporation, 2001

John P Boyd.Chebyshev and Fourier spectral methods. Courier Corporation, 2001

work page 2001

[5] [5]

A superconvergent LDG- hybridizable Galerkin method for second-order elliptic problems.Mathematics of Com- putation, 77(264):1887–1916, May 2008

Bernardo Cockburn, Bo Dong, and Johnny Guzmán. A superconvergent LDG- hybridizable Galerkin method for second-order elliptic problems.Mathematics of Com- putation, 77(264):1887–1916, May 2008

work page 1916

[6] [6]

A Hybridizable Discontinuous Galerkin Method for Steady-State Convection-Diffusion- Reaction Problems.SIAM Journal on Scientific Computing, 31(5):3827–3846, January 2009

Bernardo Cockburn, Bo Dong, Johnny Guzmán, Marco Restelli, and Riccardo Sacco. A Hybridizable Discontinuous Galerkin Method for Steady-State Convection-Diffusion- Reaction Problems.SIAM Journal on Scientific Computing, 31(5):3827–3846, January 2009

work page 2009

[7] [7]

The Derivation of Hybridizable Dis- continuous Galerkin Methods for Stokes Flow.SIAM Journal on Numerical Analysis, 47(2):1092–1125, January 2009

Bernardo Cockburn and Jayadeep Gopalakrishnan. The Derivation of Hybridizable Dis- continuous Galerkin Methods for Stokes Flow.SIAM Journal on Numerical Analysis, 47(2):1092–1125, January 2009

work page 2009

[8] [8]

Bernardo Cockburn, Jayadeep Gopalakrishnan, and Raytcho Lazarov. Unified Hy- bridization of Discontinuous Galerkin, Mixed, and Continuous Galerkin Methods for Second Order Elliptic Problems.SIAM Journal on Numerical Analysis, 47(2):1319– 1365, January 2009

work page 2009

[9] [9]

Fasoulas, C.-D

S. Fasoulas, C.-D. Munz, M. Pfeiffer, J. Beyer, T. Binder, S. Copplestone, A. Mirza, P. Nizenkov, P. Ortwein, and W. Reschke. Combining particle-in-cell and direct sim- ulation Monte Carlo for the simulation of reactive plasma flows.Physics of Fluids, 31(7):072006, July 2019

work page 2019

[10] [10]

Hairer, Christian Lubich, and Gerhard Wanner.Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations

E. Hairer, Christian Lubich, and Gerhard Wanner.Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations. Number 31 in Springer Series in Computational Mathematics. Springer, Berlin ; New York, 2nd ed. edition, 2006

work page 2006

[11] [11]

CRC Press, 1st ed

R.W Hockney and J.W Eastwood.Computer Simulation Using Particles. CRC Press, 1st ed. edition, 1988. 22

work page 1988

[12] [12]

HYPRE: High performance preconditioners.https://llnl.gov/casc/hypre,https: //github.com/hypre-space/hypre

work page

[13] [13]

Wiley, New York, NY, 3rd ed

John David Jackson.Classical electrodynamics. Wiley, New York, NY, 3rd ed. edition, 1999

work page 1999

[14] [14]

Jacobs and J.S

G.B. Jacobs and J.S. Hesthaven. High-order nodal discontinuous Galerkin particle-in- cell method on unstructured grids.Journal of Computational Physics, 214(1):96–121, May 2006

work page 2006

[15] [15]

The Bohm criterion and sheath formation.Journal of Physics D: Applied Physics, 24(4):493–518, April 1991

K -U Riemann. The Bohm criterion and sheath formation.Journal of Physics D: Applied Physics, 24(4):493–518, April 1991

work page 1991

[16] [16]

Blind, Anna Schwarz, Marius Kurz, Felix Rodach, Stephen M

Patrick Kopper, Marcel P. Blind, Anna Schwarz, Marius Kurz, Felix Rodach, Stephen M. Copplestone, and Andrea D. Beck. Pyhope: A python toolkit for three-dimensional unstructured high-order meshes.Journal of Open Source Software, 10(115):8769, 2025

work page 2025

[17] [17]

Adaptive mesh strategies for the spectral element method.Com- puter Methods in Applied Mechanics and Engineering, 116(1-4):77–86, January 1994

Catherine Mavriplis. Adaptive mesh strategies for the spectral element method.Com- puter Methods in Applied Mechanics and Engineering, 116(1-4):77–86, January 1994

work page 1994

[18] [18]

Coupled Particle-In-Cell and Direct Sim- ulation Monte Carlo method for simulating reactive plasma flows.Comptes Rendus

Claus-Dieter Munz, Monika Auweter-Kurtz, Stefanos Fasoulas, Asim Mirza, Philip Or- twein, Marcel Pfeiffer, and Torsten Stindl. Coupled Particle-In-Cell and Direct Sim- ulation Monte Carlo method for simulating reactive plasma flows.Comptes Rendus. Mécanique, 342(10-11):662–670, August 2014

work page 2014

[19] [19]

Pfeiffer, F

M. Pfeiffer, F. Hindenlang, T. Binder, S.M. Copplestone, C.-D. Munz, and S. Fasoulas. A Particle-in-Cell solver based on a high-order hybridizable discontinuous Galerkin spectral element method on unstructured curved meshes.Computer Methods in Applied Mechanics and Engineering, 349:149–166, June 2019

work page 2019