2020/04/19 by Sofia Eriksson, Eriksson, Sofia
Computer Science · Engineering · Mathematics · #Advanced Numerical Methods in Computational Mathematics #Electromagnetic Simulation and Numerical Methods #FOS: Mathematics #Matrix Theory and Algorithms #Numerical Analysis (math.NA) #cs.NA #math.NA
paper · pdf · doi:10.48550/arxiv.2004.08938
arxiv created 2020/04/19 · openalex publication_date 2020/04/19 · arxiv updated 2020/04/21 · openalex created_date 2022/07/26 · openalex updated_date 2026/07/28
The scalar, one-dimensional advection equation and heat equation are considered. These equations are discretized in space, using a finite difference method satisfying summation-by-parts (SBP) properties. To impose the boundary conditions, we use a penalty method called simultaneous approximation term (SAT). Together, this gives rise to two semi-discrete schemes where the discretization matrices approximate the first and the second derivative operators, respectively. The discretization matrices depend on free parameters from the SAT treatment. We derive the inverses of the discretization matrices, interpreting them as discrete Green's functions. In this direct way, we also find out precisely which choices of SAT parameters that make the discretization matrices singular. In the second derivative case, it is shown that if the penalty parameters are chosen such that the semi-discrete scheme is dual consistent, the discretization matrix can become singular even when the scheme is energy stable. The inverse formulas hold for SBP-SAT operators of arbitrary order of accuracy. For second and fourth order accurate operators, the inverses are provided explicitly.