Exponential convergence through linear finite element discretization of stratified subdomains

DOI10.1016/J.JCP.2016.06.045zbMATH Open1352.65565arXiv1507.05038OpenAlexW2281052490WikidataQ57430343 ScholiaQ57430343MaRDI QIDQ729601FDOQ729601

Ali Vaziri Astaneh, V. Druskin, Murthy N. Guddati

Publication date: 20 December 2016

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: Motivated by problems where the response is needed at select localized regions in a large computational domain, we devise a novel finite element discretization that results in exponential convergence at pre-selected points. The two key features of the discretization are (a) use of midpoint integration to evaluate the contribution matrices, and (b) an unconventional bending of the mesh into complex space. Named complex-length finite element method (CFEM), the technique is linked to Pade approximants that provide exponential convergence of the Dirichlet-to-Neumann maps and thus the solution at specified points in the domain. Exponential convergence facilitates drastic reduction in the number of elements. This, combined with sparse computation associated with linear finite elements, results in significant reduction in the computational cost. The paper presents the basic ideas of the method as well as illustration of its effectiveness for a variety of problems involving Laplace, Helmholtz and elastodynamic equations.

Full work available at URL: https://arxiv.org/abs/1507.05038

zbMATH Keywords

rational approximation[https://portal.mardi4nfdi.de/w/index.php?title=+Special%3ASearch&search=Pad%EF%BF%BD%EF%BF%BD+approximants&go=Go Pad�� approximants]spectral element methodsoptimal gridsDirichlet-to-Neumann maps

Mathematics Subject Classification ID

Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Spectral, collocation and related methods for boundary value problems involving PDEs (65N35)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Spectral Methods in MATLAB
• A Practical Guide to Pseudospectral Methods
• Gaussian spectral rules for second order finite-difference schemes
• Finite element analysis of acoustic scattering
• EQUILIBRIUM DISTRIBUTIONS AND DEGREE OF RATIONAL APPROXIMATION OF ANALYTIC FUNCTIONS
• On the zeros and poles of Padè approximants to \(e^z\). III
• Application of the difference Gaussian rules to solution of hyperbolic problems
• Gaussian Spectral Rules for the Three-Point Second Differences: I. A Two-Point Positive Definite Problem in a Semi-Infinite Domain
• Three-point finite-difference schemes, Padé and the spectral Galerkin method. I. One-sided impedance approximation
• Optimal finite difference grids and rational approximations of the square root. I: Elliptic problems
• Continued fraction absorbing boundary conditions for convex polygonal domains
• Efficient computation of dispersion curves for multilayered waveguides and half-spaces
• Padded continued fraction absorbing boundary conditions for dispersive waves
• Arbitrarily wide-angle wave equations for complex media
• Relative rational approximation of the function \(e^x\)
• Application of the difference Gaussian rules to solution of hyperbolic problems. II: Global expansion
• Optimal grids for anisotropic problems

Cited In (1)

• Efficient computation of dispersion curves for multilayered waveguides and half-spaces

Uses Software

• Matlab

This page was built for publication: Exponential convergence through linear finite element discretization of stratified subdomains

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q729601)