A Numerical Approach to Space-Time Finite Elements for the Wave Equation

We study a space-time finite element approach for the nonhomogeneous wave equation using a continuous time Galerkin method. We present fully implicit examples in 1+1, 2+1, and 3+1 dimensions using linear quadrilateral, hexahedral, and tesseractic elements. Krylov solvers with additive Schwarz preconditioning are used for solving the linear system. We introduce a time decomposition strategy in preconditioning which significantly improves performance when compared with unpreconditioned cases.Comment: 9 pages, 5 figures, 5 table

A finite element data assimilation method for the wave equation

We design a primal-dual stabilized finite element method for the numerical approximation of a data assimilation problem subject to the acoustic wave equation. For the forward problem, piecewise affine, continuous, finite element functions are used for the approximation in space and backward differentiation is used in time. Stabilizing terms are added on the discrete level. The design of these terms is driven by numerical stability and the stability of the continuous problem, with the objective of minimizing the computational error. Error estimates are then derived that are optimal with respect to the approximation properties of the numerical scheme and the stability properties of the continuous problem. The effects of discretizing the (smooth) domain boundary and other perturbations in data are included in the analysis.Comment: 23 page

Stability of space-time isogeometric methods for wave propagation problems

This thesis aims at investigating the first steps toward an unconditionally stable space-time isogeometric method, based on splines of maximal regularity, for the linear acoustic wave equation. The unconditional stability of space-time discretizations for wave propagation problems is a topic of significant interest, by virtue of the advantages of space-time methods compared with more standard time-stepping techniques. In the case of continuous finite element methods, several stabilizations have been proposed. Inspired by one of these works, we address the stability issue by studying the isogeometric method for an ordinary differential equation closely related to the wave equation. As a result, we provide a stabilized isogeometric method whose effectiveness is supported by numerical tests. Motivated by these results, we conclude by suggesting an extension of this stabilization tool to the space-time isogeometric formulation of the acoustic wave equation.Comment: Masters thesi

Unified error analysis for non-conforming space discretizations ofwave-type equations

This paper provides a unified error analysis for non-conforming space discretizations of linear wave equations in time-domain. We propose a framework which studies wave equations as first-order evolution equations in Hilbert spaces and their space discretizations as differential equations in finite dimensional Hilbert spaces. A lift operator maps the semi-discrete solution from the approximation space to the continuous space. Our main results are a priori error bounds in terms of interpolation, data and conformity errors of the method. Such error bounds are the key to the systematic derivation of convergence rates for a large class of problems. To show that this approach significantly eases the proof of new convergence rates, we apply it to an isoparametric finite element discretization of the wave equation with acoustic boundary conditions in a smooth domain. Moreover, our results reproduce known convergence rates for already investigated conforming and non-conforming space discretizations in a concise and unified way. The examples discussed in this paper comprise discontinuous Galerkin discretizations of Maxwell’s equations and finite elements with mass lumping for the acoustic wave equation

Higher-Order Space-Time Continuous Galerkin Methods for the Wave Equation

We consider a space-time variational formulation of the second-order wave equation, where integration by parts is also applied with respect to the time variable. Conforming tensor-product finite element discretisations with piecewise polynomials of this space-time variational formulation require a CFL condition to ensure stability. To overcome this restriction in the case of piecewise multilinear, continuous ansatz and test functions, a stabilisation is well-known, which leads to an unconditionally stable space-time finite element method. In this work, we generalise this stabilisation idea from the lowest-order case to the higher-order case, i.e. to an arbitrary polynomial degree. We give numerical examples for a one-dimensional spatial domain, where the unconditional stability and optimal convergence rates in space-time norms are illustrated.Comment: 11 page

A finite element data assimilation method for the wave equation

We design a primal-dual stabilized finite element method for the numerical approximation of a data assimilation problem subject to the acoustic wave equation. For the forward problem, piecewise affine, continuous, finite element functions are used for the approximation in space and backward differentiation is used in time. Stabilizing terms are added on the discrete level. The design of these terms is driven by numerical stability and the stability of the continuous problem, with the objective of minimizing the computational error. Error estimates are then derived that are optimal with respect to the approximation properties of the numerical scheme and the stability properties of the continuous problem. The effects of discretizing the (smooth) domain boundary and other perturbations in data are included in the analysis

A unified error analysis for spatial discretizations of wave-type equations with applications to dynamic boundary conditions

This thesis provides a unified framework for the error analysis of non-conforming space discretizations of linear wave equations in time-domain, which can be cast as symmetric hyperbolic systems or second-order wave equations. Such problems can be written as first-order evolution equations in Hilbert spaces with linear monotone operators. We employ semigroup theory for the well-posedness analysis and to obtain stability estimates for the space discretizations. To compare the finite dimensional approximations with the original solution, we use the concept of a lift from the discrete to the continuous space. Time integration with the Crank–Nicolson method is analyzed. In this framework, we derive a priori error bounds for the abstract space semi-discretization in terms of interpolation and discretization errors. These error bounds yield previously unkown convergence rates for isoparametric finite element discretizations of wave equations with dynamic boundary conditions in smooth domains. Moreover, our results allow to consider already investigated space discretizations in a unified way. Here it successfully reproduces known error bounds. Among the examples which we dicuss in this thesis are discontinuous Galerkin discretizations of Maxwell’s equations and finite elements with mass lumping for the scalar wave equation

WAVEx: Stabilized Finite Elements for Spectral Wind Wave Models Using FEniCSx

Several potential FEM discretizations of the Wave Action Balance Equation are discussed. The methods, which include streamline upwind Petrov-Galerkin (SUPG), least squares, and discontinuous Galerkin, are implemented using the open source finite element library FEniCSx for simplified 2-D cases. Open source finite element libraries, such as FEniCSx, typically only support geometries up to dimension of 3. The Wave Action Balance Equation is 4 dimensions in space so this presents difficulties. A method to use a FEM library, such as FEniCSx, to solve problems in domains with dimension larger than 4 using the product basis is discussed. A new spectral wind wave model, WAVEx, is formulated and implemented using the new finite element library FEniCSx. WAVEx is designed to allow for construction of multiple FEM discretizations with relatively small modifications in the Python code base. An example implementation is then demonstrated with WAVEx using continuous finite elements and SUPG stabilization in geographic/spectral space. For propagation in time, a generalized one step implicit finite difference method is used. When source terms are active, the second order operator splitting scheme known as Strang splitting is used. In the splitting scheme, propagation is solved using the aforementioned implicit method and the nonlinear source terms are treated explicitly using second order Runge-Kutta. Several test cases which are part of the Office for Naval Research Test Bed (ONR Test Bed) are demonstrated both with and without 3rd generation source terms and results are compared to analytic solutions, observations, and SWAN output