Abstract:In this work we explore the fidelity of numerical approximations to continuous spectra of hyperbolic partial differential equation systems. We are particularly interested in the ability of discrete methods to accurately discover sources of physical instabilities. By focusing on the perturbed equations that arise in linearized problems, we apply high-order accurate summation-by-parts finite difference operators, with weak enforcement of boundary conditions through the simulataneous-approximation-term technique, which leads to a provably stable numerical discretization with formal order of accuracy given by $p = 2, 3, 4$ and $5$. We derive analytic solutions using Laplace transform methods, which provide important ground truth for ensuring numerical convergence at the correct theoretical rate. We find that the continuous spectrum is better captured with mesh refinement, although dissipative strict stability (where the growth rate of the discrete problem is bounded above by the continuous) is not obtained. However, we also find that sole reliance on mesh refinement can be a problematic means for determining physical growth rates as some eigenvalues emerge (and persist with mesh refinement) based on spatial order of accuracy but are non-physical. We suggest that numerical methods be used to approximate discrete spectra when numerical stability is guaranteed and convergence of the discrete spectra is evident with both mesh refinement and increasing order of accuracy.

Spectral Convergence of a Semi-discretized Numerical System for the Spatially Homogeneous Boltzmann Equation with Uncertainties

Spectral convergence of a semi-discretized numerical system for the spatially homogeneous Boltzmann equation with uncertainties

Convergence of the Fourier-Galerkin spectral method for the Boltzmann equation with uncertainties

A stochastic Galerkin method for the Boltzmann equation with uncertainty

Spectral convergence of probability densities for forward problems in uncertainty quantification

Uncertainty Quantification for the Homogeneous Landau-Fokker-Planck Equation via Deterministic Particle Galerkin methods

A stochastic Galerkin method for general system of quasilinear hyperbolic conservation laws with uncertainty

Sparse grid approximation of stochastic parabolic PDEs: The Landau--Lifshitz--Gilbert equation

On Semi-classical Limit of Spatially Homogeneous Quantum Boltzmann Equation: Weak Convergence

On Deterministic Numerical Methods for the Quantum Boltzmann-Nordheim Equation. I. Spectrally Accurate Approximations, Bose-Einstein Condensation, Fermi-Dirac Saturation

Uncertainty Quantification for the BGK Model of the Boltzmann Equation Using Multilevel Variance Reduced Monte Carlo Methods

The Discrete Stochastic Galerkin Method for Hyperbolic Equations with Non-smooth and Random Coefficients

Preserving Hyperbolicity in Stochastic Galerkin Method for Uncertainty Quantification

Sensitivity analysis and incompressible Navier-Stokes-Poisson limit of Vlasov-Poisson-Boltzmann equations with uncertainty

Absolute continuity and numerical approximation of stochastic Cahn--Hilliard equation with unbounded noise diffusion

Energy estimates and hypocoercivity analysis for a multi-phase Navier-Stokes-Vlasov-Fokker-Planck system with uncertainty

Numerical analysis of a full discretization for stochastic Cahn--Hilliard equation driven by additive noise

Uncertainty quantification for random domains using periodic random variables

New High-Order Numerical Methods for Hyperbolic Systems of Nonlinear PDEs with Uncertainties

A sample-based spectral method approach for solving high-dimensional SPDEs

Approximating Continuous Spectra of Hyperbolic Systems with Summation-by-Parts Finite Difference Operators