Feb 28 2017

math.NA arXiv:1702.08385v1

A coercivity property of temporal convolution operators is an essential tool in the analysis of time-dependent boundary integral equations and their space and time discretisations. It is known that this coercivity property is inherited by convolution quadrature time discretisation based on A-stable multistep methods, which are of order at most two. Here we study the question as to which Runge--Kutta-based convolution quadrature methods inherit the convolution coercivity property. It is shown that this holds without any restriction for the third-order Radau IIA method, and on permitting a shift in the Laplace domain variable, this holds for all algebraically stable Runge--Kutta methods and hence for methods of arbitrary order. As an illustration, the discrete convolution coercivity is used to analyse the stability and convergence properties of the time discretisation of a non-linear boundary integral equation that originates from a non-linear scattering problem for the linear wave equation. Numerical experiments illustrate the error behaviour of the Runge--Kutta convolution quadrature time discretisation.

In a series of papers the present authors and their coworkers have developed a family of algebraic techniques to solve a number of problems in the theory of discrete or continuous dynamical systems and to analyze numerical integrators. Given a specific problem, those techniques construct an abstract, \em universal version of it which is solved algebraically; then, the results are tranferred to the original problem with the help of a suitable morphism. In earlier contributions, the abstract problem is formulated either in the dual of the shuffle Hopf algebra or in the dual of the Connes-Kreimer Hopf algebra. In the present contribution we extend these techniques to more general Hopf algebras, which in some cases lead to more efficient computations.

Feb 28 2017

math.NA arXiv:1702.08340v1

In this paper we will discuss different coupling methods suitable for use in the framework of the recently introduced CutFEM paradigm, cf. Burman, Erik; Claus, Susanne; Hansbo, Peter; Larson, Mats G.; Massing, André . CutFEM: discretizing geometry and partial differential equations. Internat. J. Numer. Methods Engrg. 104 (2015), no. 7, 472-501. In particular we will consider mortaring using Lagrange multipliers on the one hand and Nitsche's method on the other. For simplicity we will first discuss these method in the setting of uncut meshes, and end with some comments on the extension to CutFEM. We will, for comparison, discuss some different types of problems such as high contrast problems and problems with stiff coupling or adhesive contact. We will review some of the existing methods for these problems and propose some alternative methods resulting from crossovers from the Lagrange multiplier framework to Nitsche's method and vice versa.

Feb 28 2017

math.NA arXiv:1702.08241v1

In this paper two types of multgrid methods, i.e., the Rayleigh quotient iteration and the inverse iteration with fixed shift, are developed for solving the Maxwell eigenvalue problem with discontinuous relative magnetic permeability and electric permittivity. With the aid of the mixed form of source problem associated with the eigenvalue problem, we prove the uniform convergence of the discrete solution operator to the solution operator in $L^2(\Omega)$ using discrete compactness of edge element space. Then we prove the asymptotically optimal error estimates for both multigrid methods. Numerical experiments confirm our theoretical analysis.

Feb 28 2017

math.NA arXiv:1702.08224v1

We propose a novel Hybrid High-Order method for the Cahn-Hilliard problem with convection. The proposed method is valid in two and three space dimensions, and it supports arbitrary approximation orders on general meshes containing polyhedral elements and nonmatching interfaces. An extensive numerical validation is presented, which shows robustness with respect to the Péclet number.

Feb 28 2017

math.NA arXiv:1702.08179v1

The biharmonic operator plays a central role in a wide array of physical models, notably in elasticity theory and the streamfunction formulation of the Navier-Stokes equations. The need for corresponding numerical simulations has led, in recent years, to the development of a discrete biharmonic calculus. The primary object of this calculus is a high-order compact discrete biharmonic operator (DBO). The numerical results have been remarkably accurate, and have been corroborated by some rigorous proofs. However, there remained the "mystery" of the "underlying reason" for this success. This paper is a contribution in this direction, expounding the strong connection between cubic spline functions (on an interval) and the DBO. It is shown in particular that the (scaled) fourth-order distributional derivative of the cubic spline is identical to the action of the DBO on grid functions. The DBO is constructed in terms of the discrete Hermitian derivative. A remarkable fact is that the kernel of the inverse of the discrete operator is (up to scaling) equal to the grid evaluation of the kernel of $\Big[\Big(\frac{d}{dx}\Big)^4\Big]^{-1} .$ Explicit expressions are presented for both kernels. The relation between the (infinite) set of eigenvalues of the fourth-order Sturm-Liouville problem and the finite set of eigenvalues of the discrete biharmonic operator is studied, and the discrete eigenvalues are proved to converge (at an "optimal" $O(h^4)$ rate) to the continuous ones. Another remarkable consequence is the validity of a comparison principle. It is well known that there is no maximum principle for the fourth-order equation. However, a positivity result is derived, both for the continuous and the discrete biharmonic equation, showing that in both cases the kernels are order preserving.

Under the strongly convex assumption, several recent works studied the global linear convergence rate of the proximal incremental aggregated gradient (PIAG) method for minimizing the sum of a large number of smooth component functions and a non-smooth convex function. In this paper, under \textslthe quadratic growth condition--a strictly weaker condition than the strongly convex assumption, we derive a new global linear convergence rate result, which implies that the PIAG method attains global linear convergence rates in both the function value and iterate point errors. The main idea behind is to construct a certain Lyapunov function.

Feb 28 2017

math.NA arXiv:1702.08111v1

We present a new version of the Method of Local Corrections (MLC) \citemlc, a multilevel, low communications, non-iterative, domain decomposition algorithm for the numerical solution of the free space Poisson's equation in 3D on locally-structured grids. In this method, the field is computed as a linear superposition of local fields induced by charges on rectangular patches of size $O(1)$ mesh points, with the global coupling represented by a coarse grid solution using a right-hand side computed from the local solutions. In the present method, the local convolutions are further decomposed into a short-range contribution computed by convolution with the discrete Green's function for an $Q^{th}$-order accurate finite difference approximation to the Laplacian with the full right-hand side on the patch, combined with a longer-range component that is the field induced by the terms up to order $P-1$ of the Legendre expansion of the charge over the patch. This leads to a method with a solution error that has an asymptotic bound of $O(h^P) + O(h^Q) + O(\epsilon h^2) + O(\epsilon)$, where $h$ is the mesh spacing, and $\epsilon$ is the max norm of the charge times a rapidly-decaying function of the radius of the support of the local solutions scaled by $h$. Thus we have eliminated the low-order accuracy of the original method (which corresponds to $P=1$ in the present method) for smooth solutions, while keeping the computational cost per patch nearly the same with that of the original method. Specifically, in addition to the local solves of the original method we only have to compute and communicate the expansion coefficients of local expansions (that is, for instance, 20 scalars per patch for $P=4$). Several numerical examples are presented to illustrate the new method and demonstrate its convergence properties.

Feb 28 2017

math.NA arXiv:1702.07962v1

We propose boundary conditions for the diffusion equation that maintain the initial mean and the total mass of a discrete data sample in the density estimation process. A complete study of this framework with numerical experiments using the finite element method is presented for the one dimensional diffusion equation, some possible applications of this results are presented as well. We also comment on a similar methodology for the two-dimensional diffusion equation for future applications in two-dimensional domains.

Feb 28 2017

math.NA arXiv:1702.07894v1

We present an analysis of the ensemble Kalman filter for inverse problems based on the continuous time limit of the algorithm. The analysis of the dynamical behaviour of the ensemble allows to establish well-posedness and convergence results for a fixed ensemble size. We will build on the results presented in [Schillings, Stuart 2017] and generalise them to the case of noisy observational data, in particular the influence of the noise on the convergence will be investigated, both theoretically and numerically.

Feb 28 2017

math.NA arXiv:1702.07757v1

In this work, several multilevel decoupled algorithms are proposed for a mixed Navier-Stokes/Darcy model. These algorithms are based on either successively or parallelly solving two linear subdomain problems after solving a coupled nonlinear coarse grid problem. Error estimates are given to demonstrate the approximation accuracy of the algorithms. Experiments based on both the first order and the second order discretizations are presented to show the effectiveness of the decoupled algorithms.

We formulate a general criterion for the exact preservation of the "lake at rest" solution in general mesh-based and meshless numerical schemes for the strong form of the shallow-water equations with bottom topography. The main idea is a careful mimetic design for the spatial derivative operators in the momentum flux equation that is paired with a compatible averaging rule for the water column height arising in the bottom topography source term. We prove consistency of the mimetic difference operators analytically and demonstrate the well-balanced property numerically using finite difference and RBF-FD schemes in the one- and two-dimensional cases.