Very Weak Space-Time Variational Formulation for the Wave Equation: Analysis and Efficient Numerical Solution

We introduce a very weak space-time variational formulation for the wave equation, prove its well-posedness (even in the case of minimal regularity) and optimal inf-sup stability. Then, we introduce a tensor product-style space-time Petrov-Galerkin discretization with optimal discrete inf-sup stability, obtained by a non-standard definition of the trial space. As a consequence, the numerical approximation error is equal to the residual, which is particularly useful for a posteriori error estimation. For the arising discrete linear systems in space and time, we introduce efficient numerical solvers that appropriately exploit the equation structure, either at the preconditioning level or in the approximation phase by using a tailored Galerkin projection. This Galerkin method shows competitive behavior concerning wall-clock time, accuracy and memory as compared with a standard time-stepping method in particular in low regularity cases. Numerical experiments with a 3D (in space) wave equation illustrate our findings.

[1]  Daniel Kressner,et al.  Low-Rank Tensor Krylov Subspace Methods for Parametrized Linear Systems , 2011, SIAM J. Matrix Anal. Appl..

[2]  Karsten Urban,et al.  (Parametrized) First Order Transport Equations: Realization of Optimally Stable Petrov-Galerkin Methods , 2018, SIAM J. Sci. Comput..

[3]  R. Nochetto,et al.  Theory of adaptive finite element methods: An introduction , 2009 .

[4]  Ludmil T. Zikatanov,et al.  Some observations on Babu\vs}ka and Brezzi theories , 2003, Numerische Mathematik.

[5]  I. Babuska,et al.  Analysis of mixed methods using mesh dependent norms , 1980 .

[6]  Leszek Demkowicz,et al.  A class of discontinuous Petrov-Galerkin methods. Part III , 2012 .

[7]  Valeria Simoncini,et al.  Adaptive rational Krylov subspaces for large-scale dynamical systems , 2011, Syst. Control. Lett..

[8]  Leszek Demkowicz,et al.  A class of discontinuous Petrov–Galerkin methods. II. Optimal test functions , 2011 .

[9]  Valeria Simoncini,et al.  Computational Methods for Linear Matrix Equations , 2016, SIAM Rev..

[10]  C. Mollet,et al.  Parabolic PDEs in Space-Time Formulations: Stability for Petrov-Galerkin Discretizations with B-Splines and Existence of Moments for Problems with Random Coefficients , 2016 .

[11]  Omar Ghattas,et al.  Constructively well-posed approximation methods with unity inf-sup and continuity constants for partial differential equations , 2013, Math. Comput..

[12]  Wolfgang Dahmen,et al.  Adaptive Petrov-Galerkin Methods for First Order Transport Equations , 2011, SIAM J. Numer. Anal..

[13]  Sandra May,et al.  Error Analysis for a Finite Element Approximation of Elliptic Dirichlet Boundary Control Problems , 2013, SIAM J. Control. Optim..

[14]  David J. Silvester,et al.  A Black-Box Multigrid Preconditioner for the Biharmonic Equation , 2004 .

[15]  Karsten Urban,et al.  Matrix oriented reduction of space-time Petrov-Galerkin variational problems , 2019, ArXiv.

[16]  A. Quarteroni,et al.  Reduced Basis Methods for Partial Differential Equations: An Introduction , 2015 .

[17]  Karsten Urban,et al.  A new error bound for reduced basis approximation of parabolic partial differential equations , 2012 .

[18]  Davide Palitta.  Matrix equation techniques for certain evolutionary partial differential equations , 2019, ArXiv.

[19]  Martin Berggren,et al.  Approximations of Very Weak Solutions to Boundary-Value Problems , 2003, SIAM J. Numer. Anal..

[20]  Karsten Urban,et al.  An improved error bound for reduced basis approximation of linear parabolic problems , 2013, Math. Comput..

[21]  J. Hesthaven,et al.  Certified Reduced Basis Methods for Parametrized Partial Differential Equations , 2015 .

[22]  J. Lions,et al.  Non-homogeneous boundary value problems and applications , 1972 .

[23]  Roman Andreev,et al.  Stability of sparse space–time finite element discretizations of linear parabolic evolution equations , 2013 .

[24]  Valeria Simoncini,et al.  An Efficient Reduced Basis Solver for Stochastic Galerkin Matrix , 2017 .

[25]  S. Nicaise,et al.  Discretization of the Poisson equation with non‐smooth data and emphasis on non‐convex domains , 2015, 1505.01229.