A convolution quadrature Galerkin boundary element method for the exterior Neumann problem of the wave equation

The numerical solution of the Neumann problem of the wave equation on unbounded three‐dimensional domains is calculated using the convolution quadrature method for the time discretization and a Galerkin boundary element method for the spatial discretization. The mathematical analysis that has been built up for the Dirichlet problem is extended and developed for the Neumann problem, which is important for many modelling applications. Numerical examples are then presented for one of these applications, modelling transient acoustic radiation. Copyright © 2009 John Wiley & Sons, Ltd.     

Quasi‐local multitrace boundary integral formulations

In the context of time harmonic wave scattering by piecewise homogeneous penetrable objects, we present a new variant of the multitrace boundary integral formulation, that differs from the local multitrace approach of (Jerez‐Hanckes and Hiptmair, 2012) [13] by the presence of regularization terms involving boundary integral operators and localized at junctions i.e. points where at least three subdomains abut. We prove well‐posedness and quasioptimal convergence of conforming Galerkin discretizations for this new formulation, and present numerical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2043–2062, 2015     

Convolution quadrature time-domain boundary element method for 2-D fluid-saturated porous media

For large-scale wave analyses of fluid-saturated porous media, a conventional time-domain boundary element method (BEM) cannot be applied because of the following reasons: (1) no time-domain fundamental solutions are known for some problems, (2) the method sometimes suffers from instability, and (3) the analyses require large amounts of computational time and memory. In this study, an innovative time-domain BEM is developed for a fluid-saturated porous medium. The formulation presented herein overcomes the above disadvantages using a convolution quadrature method (CQM), first proposed by Lubich, and hybrid-parallelization with both MPI and OpenMP. Problems involving the scattering of an incident plane wave by cavities in a 2-D poroelastic medium are solved as a means of validating the proposed method.     

Estimating quadrature errors for an efficient method for quasi-singular boundary integral

The numerical resolution of the boundary integral equations applied to the differential equations of Laplace, Helmholtz and Maxwell requires the handling of quasi-singular integrals with different order of singularity. The numerical approximation of the integral equations of different kinds is made by boundary finite elements. In this paper, we present a complete survey for estimating quadrature errors for the numerical techniques proposed by Huang and Cruse [Q. Huang, T.A. Cruse, Some notes on singular integral techniques in boundary element analysis, Int. J. Numer. Methods Eng. 36 (15) (1993) 2643-2659], to calculate the quasi-singular integrals. To validate the accuracy and efficiency of these techniques and approve our study some numerical examples are presented and discussed.     

Efficient long-time computations of time-domain boundary integrals for 2D and dissipative wave equation

Linear hyperbolic partial differential equations in a homogeneous medium, e.g., the wave equation describing the propagation and scattering of acoustic waves, can be reformulated as time-domain boundary integral equations. We propose an efficient implementation of a numerical discretization of such equations when the strong Huygens’ principle does not hold.For the numerical discretization, we make use of convolution quadrature in time and standard Galerkin boundary element method in space. The quadrature in time results in a discrete convolution of weights W_j with the boundary density evaluated at equally spaced time points. If the strong Huygens’ principle holds, W_j converge to 0 exponentially quickly for large enough j. If the strong Huygens’ principle does not hold, e.g., in even space dimensions or when some damping is present, the weights are never zero, thereby presenting a difficulty for efficient numerical computation.In this paper we prove that the kernels of the convolution weights approximate in a certain sense the time domain fundamental solution and that the same holds if both are differentiated in space. The tails of the fundamental solution being very smooth, this implies that the tails of the weights are smooth and can efficiently be interpolated. Further, we hint on the possibility to apply the fast and oblivious convolution quadrature algorithm of Schädle et al. to further reduce memory requirements for long-time computation. We discuss the efficient implementation of the whole numerical scheme and present numerical experiments.     

Numerical treatment of retarded boundary integral equations by sparse panel clustering

Stefan Sauter Institute for Mathematics, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland ^{We consider the wave equation in a boundary integral formulation.The discretization in time is done by using convolution quadraturetechniques and a Galerkin boundary element method for the spatialdiscretization. In a previous paper, we have introduced a sparseapproximation of the system matrix by cut-off, in order to reducethe storage costs. In this paper, we extend this approach byintroducing a panel clustering method to further reduce thesecosts.}     

Differential quadrature domain decomposition method for problems on a triangular domain

The differential quadrature method (DQM) has been studied for years and it has been shown by many researchers that the DQM is an attractive numerical method with high efficiency and accuracy. The conventional DQM is mostly effective for one‐dimensional and multidimensional problems with geometrically regular domains. But to deal with problems on a triangular domain, we will meet difficulties. In this article we will study how to solve problems on a triangular domain by using DQM combined with the domain decomposition method (DDM). © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Mechanical quadrature methods and their splitting extrapolations for boundary integral equations of first kind on open arcs

This paper presents the mechanical quadrature methods (MQMs) for solving boundary integral equations (BIEs) of the first kind on open arcs. The spectral condition number of MQMs is only O(h⁻¹), where h is the maximal mesh width. The errors of MQMs have multivariate asymptotic expansions, accompanied with for all mesh widths h_i. Hence, once discrete equations with coarse meshes are solved in parallel, the accuracy order of numerical approximations can be greatly improved by splitting extrapolation algorithms (SEAs). Moreover, a posteriori asymptotic error estimates are derived, which can be used to formulate self-adaptive algorithms. Numerical examples are also provided to support our algorithms and analysis. Furthermore, compared with the existing algorithms, such as Galerkin and collocation methods, the accuracy order of the MQMs is higher, and the discrete matrix entries are explicit, to prove that the MQMs in this paper are more promising and beneficial to practical applications.     

Modified boundary integral formulations for the Helmholtz equation

In this paper we describe some modified regularized boundary integral equations to solve the exterior boundary value problem for the Helmholtz equation with either Dirichlet or Neumann boundary conditions. We formulate combined boundary integral equations which are uniquely solvable for all wave numbers even for Lipschitz boundaries Γ=∂Ω. This approach extends and unifies existing regularized combined boundary integral formulations.     

Numerical solution of evolutionary integral equations with completely monotonic kernel by Runge–Kutta convolution quadrature

We study the numerical solutions of the initial boundary value problems for the Volterra‐type evolutionary integal equations, in which the integral operator is a convolution product of a completely monotonic kernel and a positive definite operator, such as an elliptic partial‐differential operator. The equation is discretized in time by the Runge–Kutta convolution quadrature. Error estimates are derived and numerical experiments reported. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq31: 105–142, 2015     

A Petrov-Galerkin method with quadrature for elliptic boundary value problems

We propose and analyse a fully discrete Petrov–Galerkinmethod with quadrature, for solving second-order, variable coefficient,elliptic boundary value problems on rectangular domains. Inour scheme, the trial space consists of C² splines of degreer 3, the test space consists of C⁰ splines of degree r –2, and we use composite (r – 1)-point Gauss quadrature.We show existence and uniqueness of the approximate solutionand establish optimal order error bounds in H², H¹ and L² norms.     

An improved time domain linear sampling method for Robin and Neumann obstacles

We consider inverse obstacle scattering problems for the wave equation with Robin or Neumann boundary conditions. The problem of reconstructing the geometry of such obstacles from measurements of scattered waves in the time domain is tackled using a time domain linear sampling method. This imaging technique yields a picture of the scatterer by solving a linear operator equation involving the measured data for many right-hand sides given by singular solutions to the wave equation. We analyse this algorithm for causal and smooth impulse shapes, we discuss the effect of different choices of the singular solutions used in the algorithm, and finally we propose a fast FFT-based implementation.     

A fast numerical method for a natural boundary integral equation for the Helmholtz equation

A Neumann boundary value problem of the Helmholtz equation in the exterior circular domain is reduced into an equivalent natural boundary integral equation. Using our trigonometric wavelets and the Galerkin method, the obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. Especially, our method is also efficient when the wave number k in the Helmholtz equation is very large.     

A random integral quadrature method for numerical analysis of the second kind of Volterra integral equations

In this paper, a novel meshless technique termed the random integral quadrature (RIQ) method is developed for the numerical solution of the second kind of the Volterra integral equations. The RIQ method is based on the generalized integral quadrature (GIQ) technique, and associated with the Kriging interpolation function, such that it is regarded as an extension of the GIQ technique. In the GIQ method, the regular computational domain is required, in which the field nodes are scattered along straight lines. In the RIQ method however, the field nodes can be distributed either uniformly or randomly. This is achieved by discretizing the governing integral equation with the GIQ method over a set of virtual nodes that lies along straight lines, and then interpolating the function values at the virtual nodes over all the field nodes which are scattered either randomly or uniformly. In such a way, the governing integral equation is converted approximately into a system of linear algebraic equations, which can be easily solved.     

On a nonlinear wave equation with boundary damping

This paper is concerned with the existence and decay of solutions of the mixed problem for the nonlinear wave equation with boundary conditions Here, Ω is an open bounded set of with boundary Γ of class C²; Γ is constituted of two disjoint closed parts Γ₀ and Γ₁ both with positive measure; the functions μ(t), f(s), g(s) satisfy the conditions μ(t) ≥ μ₀ > 0, f(s) ≥ 0, g(s) ≥ 0 for t ≥ 0, s ≥ 0 and h(x,s) is a real function where x ∈ Γ₁, ν(x) is the unit outward normal vector at x ∈ Γ₁ and α, β are non‐negative real constants. Assuming that h(x,s) is strongly monotone in s for each x ∈ Γ₁, it is proved the global existence of solutions for the previous mixed problem. For that, it is used in the Galerkin method with a special basis, the compactness approach, the Strauss approximation for real functions and the trace theorem for nonsmooth functions. The exponential decay of the energy is derived by two methods: by using a Lyapunov functional and by Nakao's method. Copyright © 2013 John Wiley & Sons, Ltd.     

Convolution quadrature time discretization of fractional diffusion-wave equations

We propose and study a numerical method for time discretization of linear and semilinear integro-partial differential equations that are intermediate between diffusion and wave equations, or are subdiffusive. The method uses convolution quadrature based on the second-order backward differentiation formula. Second-order error bounds of the time discretization and regularity estimates for the solution are shown in a unified way under weak assumptions on the data in a Banach space framework. Numerical experiments illustrate the theoretical results.

     

Boundary element method and wave equation

In this paper the boundary integral expression for a one-dimensional wave equation with homogeneous boundary conditions is developed. This is done using the time dependent fundamental solution of the corresponding hyperbolic partial differential equation. The integral expression developed is a generalized function with the same form as the well-known D'Alembert formula. The derivatives of the solution and some useful invariants on the characteristics of the partial differential equation are also calculated.The boundary element method is applied to find the numerical solution. The results show excellent agreement with analytical solutions.A multi-step procedure for large time steps which can be used in the boundary element method is also described.In addition, the way in which boundary conditions are introduced during the time dependent process is explained in detail. In the Appendix the main properties of Dirac's delta function and the Heaviside unit step function are described.     

Generalized convolution quadrature with variable time stepping. Part II: Algorithm and numerical results

In this paper we address the implementation of the Generalized Convolution Quadrature (gCQ) presented and analyzed by the authors in a previous paper for solving linear parabolic and hyperbolic convolution equations. Our main goal is to overcome the current restriction to uniform time steps of Lubich's Convolution Quadrature (CQ). A major challenge for the efficient realization of the new method is the evaluation of high-order divided differences for the transfer function in a fast and stable way. Our algorithm is based on contour integral representation of the numerical solution and quadrature in the complex plane. As the main application we consider the wave equation in exterior domains, which is formulated as a retarded boundary integral equation. We provide numerical experiments to illustrate the theoretical results.     

Finite difference method combined with differential quadrature method for numerical computation of the modified equal width wave equation

The aim of this study is to improve the numerical solution of the modified equal width wave equation. For this purpose, finite difference method combined with differential quadrature method with Rubin and Graves linearizing technique has been used. Modified cubic B‐spline base functions are used as base function. By the combination of two numerical methods and effective linearizing technique high accurate numerical algorithm is obtained. Three main test problems are solved for various values of the coefficients. To observe the performance of the present method, the error norms of the single soliton problem which has analytical solution are calculated. Besides these error norms, three invariants are reported. Comparison of the results displays that our algorithm produces superior results than those given in the literature.