On Spurious Mesh-Scale Oscillations Associated with Resonance-Cone Regimes in Wave-Plasma Simulations
Abstract
Spurious mesh-scale oscillations can compromise frequency-domain finite-element simulations of electromagnetic wave propagation in magnetized cold plasmas. We investigate their origin and mitigation through a verification framework combining a mixed Nédélec-Lagrange solver, a waveport accommodating non-orthogonal modes in anisotropic media, and an analytical generalized Fresnel benchmark at a vacuum-plasma interface. Comparisons with this benchmark show that the observed oscillations occur in resonance-cone regimes and are absent from the analytical fields, confirming their numerical origin. Systematic tests assess the effects of approximation order, variational formulation, mesh geometry, and numerical collisionality. Plane-wave residual and Fourier spectral analyses suggest that local mesh variations induce spectral coupling to additional wavevector components, predominantly populating the dispersion branches in collisionless resonance-cone regimes where the reduced scalar equation loses ellipticity. Uniform rectangular grids recover the analytical wave structure, further supporting the role of mesh geometry and symmetry. However, strict symmetry requirements limit the applicability of this approach, motivating numerical regularization. In resonance-cone regimes, introducing sufficient numerical collisionality suppresses the mesh-scale oscillations, yielding smooth electric-field solutions and reducing the error relative to the analytical benchmark. For one resonance-cone benchmark configuration, the error-minimizing damping decreases with mesh refinement according to an empirical Zn∝ h1.48 scaling, sustaining error reduction consistent with second-order convergence and avoiding the accuracy limit imposed by fixed collisionality. ...
Related
- Attention Is All You Need (to Avoid Spurious Oscillations)
- Experimental observation of Higgs oscillations and their power-law decay across pairing regimes in a cavity QED quantum simulator
- Probing hairy Kerr black holes through quasi-periodic oscillations II: a study based on the resonance models
- Chirpons: elementary local phase oscillations associated with complex zeros
- Spurious spin nutation modes arising from truncated memory effects
BibTeX
@misc{bergerioux2026spurious,
title = {{On Spurious Mesh-Scale Oscillations Associated with Resonance-Cone Regimes in Wave-Plasma Simulations}},
author = {Hugo Bergerioux and Mario Merino},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.31021},
eprint = {2610.31021},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.31021},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - On Spurious Mesh-Scale Oscillations Associated with Resonance-Cone Regimes in Wave-Plasma Simulations
AU - Hugo Bergerioux
AU - Mario Merino
PY - 2026
DA - 2026/10/22
PB - arXiv
DO - 10.48550/arxiv.2610.31021
UR - https://arxiv.org/abs/2610.31021
AB - Spurious mesh-scale oscillations can compromise frequency-domain finite-element simulations of electromagnetic wave propagation in magnetized cold plasmas. We investigate their origin and mitigation through a verification framework combining a mixed Nédélec-Lagrange solver, a waveport accommodating non-orthogonal modes in anisotropic media, and an analytical generalized Fresnel benchmark at a vacuum-plasma interface. Comparisons with this benchmark show that the observed oscillations occur in resonance-cone regimes and are absent from the analytical fields, confirming their numerical origin. Systematic tests assess the effects of approximation order, variational formulation, mesh geometry, and numerical collisionality. Plane-wave residual and Fourier spectral analyses suggest that local mesh variations induce spectral coupling to additional wavevector components, predominantly populating the dispersion branches in collisionless resonance-cone regimes where the reduced scalar equation loses ellipticity. Uniform rectangular grids recover the analytical wave structure, further supporting the role of mesh geometry and symmetry. However, strict symmetry requirements limit the applicability of this approach, motivating numerical regularization. In resonance-cone regimes, introducing sufficient numerical collisionality suppresses the mesh-scale oscillations, yielding smooth electric-field solutions and reducing the error relative to the analytical benchmark. For one resonance-cone benchmark configuration, the error-minimizing damping decreases with mesh refinement according to an empirical Z_n∝ h^1.48 scaling, sustaining error reduction consistent with second-order convergence and avoiding the accuracy limit imposed by fixed collisionality. ...
ER -