Spectral solutions of PDEs on networks

by   M. Brio, et al.

To solve linear PDEs on metric graphs with standard coupling conditions (continuity and Kirchhoff's law), we develop and compare a spectral, a second-order finite difference, and a discontinuous Galerkin method. The spectral method yields eigenvalues and eigenvectors of arbitary order with machine precision and converges exponentially. These eigenvectors provide a Fourier-like basis on which to expand the solution; however, more complex coupling conditions require additional research. The discontinuous Galerkin method provides approximations of arbitrary polynomial order; however computing high-order eigenvalues accurately requires the respective eigenvector to be well-resolved. The method allows arbitrary non-Kirchhoff flux conditions and requires special penalty terms at the vertices to enforce continuity of the solutions. For the finite difference method, the standard one-sided second-order finite difference stencil reduces the accuracy of the vertex solution to O(h^3/2). To preserve overall second-order accuracy, we used ghost cells for each edge. For all three methods we provide the implementation details, their validation, and examples illustrating their performance for the eigenproblem, Poisson equation, and the wave equation.



There are no comments yet.


page 1

page 2

page 3

page 4


Arbitrary Order Solutions for the Eikonal Equation using a Discontinuous Galerkin Method

We provide a method to compute the entropy-satisfying weak solution to t...

Local Fourier Analysis of P-Multigrid for High-Order Finite Element Operators

Multigrid methods are popular for solving linear systems derived from di...

Energy-based discontinuous Galerkin difference methods for second-order wave equations

We combine the newly-constructed Galerkin difference basis with the ener...

A phase-space discontinuous Galerkin approximation for the radiative transfer equation in slab geometry

We derive and analyze a symmetric interior penalty discontinuous Galerki...

A conservative and energy stable discontinuous spectral element method for the shifted wave equation in second order form

In this paper, we develop a provably energy stable and conservative disc...

Hybrid Solver for the Radiative Transport Equation Using Finite Volume and Discontinuous Galerkin

We propose a hybrid spatial discretization for the radiative transport e...

Discontinuous Galerkin methods for the Laplace-Beltrami operator on point cloud

A generic geometric error analysis framework for numerical solution of P...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.