Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn extended Filon--Clenshaw--Curtis method for high-frequency wave scattering problems in two dimensions
455
0
0.0
(
0
)
Added by
Georg Maierhofer
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We study the efficient approximation of integrals involving Hankel functions of the first kind which arise in wave scattering problems on straight or convex polygonal boundaries. Filon methods have proved to be an effective way to approximate many types of highly oscillatory integrals, however finding such methods for integrals that involve non-linear oscillators and frequency-dependent singularities is subject to a significant amount of ongoing research. In this work, we demonstrate how Filon methods can be constructed for a class of integrals involving a Hankel function of the first kind. These methods allow the numerical approximation of the integral at uniform cost even when the frequency $omega$ is large. In constructing these Filon methods we also provide a stable algorithm for computing the Chebyshev moments of the integral based on duality to spectral methods applied to a version of Bessels equation. Our design for this algorithm has significant potential for further generalisations that would allow Filon methods to be constructed for a wide range of integrals involving special functions. These new extended Filon methods combine many favourable properties, including robustness in regard to the regularity of the integrand and fast approximation for large frequencies. As a consequence, they are of specific relevance to applications in wave scattering, and we show how they may be used in practice to assemble collocation matrices for wavelet-based collocation methods and for hybrid oscillatory approximation spaces in high-frequency wave scattering problems on convex polygonal shapes.
rate research
Read More
We introduce a new method for the numerical approximation of time-harmonic acoustic scattering problems stemming from material inhomogeneities. The method works for any frequency $omega$, but is especially efficient for high-frequency problems. It is based on a time-domain approach and consists of three steps: emph{i)} computation of a suitable incoming plane wavelet with compact support in the propagation direction; emph{ii)} solving a scattering problem in the time domain for the incoming plane wavelet; emph{iii)} reconstruction of the time-harmonic solution from the time-domain solution via a Fourier transform in time. An essential ingredient of the new method is a front-tracking mesh adaptation algorithm for solving the problem in emph{ii)}. By exploiting the limited support of the wave front, this allows us to make the number of the required degrees of freedom to reach a given accuracy significantly less dependent on the frequency $omega$, as shown in the numerical experiments. We also present a new algorithm for computing the Fourier transform in emph{iii)} that exploits the reduced number of degrees of freedom corresponding to the adapted meshes.
This paper proposes an interface/boundary-unfitted eXtended hybridizable discontinuous Galerkin (X-HDG) method for Darcy-Stokes-Brinkman interface problems in two and three dimensions. The method uses piecewise linear polynomials for the velocity approximation and piecewise constants for both the velocity gradient and pressure approximations in the interior of elements inside the subdomains separated by the interface, uses piecewise constants for the numerical traces of velocity on the inter-element boundaries inside the subdomains, and uses piecewise constants or linear polynomials for the numerical traces of velocity on the interface. Optimal error estimates are derived for the interface-unfitted X-HDG scheme. Numerical experiments are provided to verify the theoretical results and the robustness of the proposed method.
An interface/boundary-unfitted eXtended hybridizable discontinuous Galerkin (X-HDG) method of arbitrary order is proposed for linear elasticity interface problems on unfitted meshes with respect to the interface and domain boundary. The method uses piecewise polynomials of degrees $k (>= 1)$ and $k-1$ respectively for the displacement and stress approximations in the interior of elements inside the subdomains separated by the interface, and piecewise polynomials of degree $k$ for the numerical traces of the displacement on the inter-element boundaries inside the subdomains and on the interface/boundary of the domain. Optimal error estimates in $L^2$-norm for the stress and displacement are derived. Finally, numerical experiments confirm the theoretical results and show that the method also applies to the case of crack-tip domain.
Consider the electromagnetic scattering of a time-harmonic plane wave by an open cavity which is embedded in a perfectly electrically conducting infinite ground plane. This paper is concerned with the numerical solutions of the transverse electric and magnetic polarizations of the open cavity scattering problems. In each polarization, the scattering problem is reduced equivalently into a boundary value problem of the two-dimensional Helmholtz equation in a bounded domain by using the transparent boundary condition (TBC). An a posteriori estimate based adaptive finite element method with the perfectly matched layer (PML) technique is developed to solve the reduced problem. The estimate takes account both of the finite element approximation error and the PML truncation error, where the latter is shown to decay exponentially with respect to the PML medium parameter and the thickness of the PML layer. Numerical experiments are presented and compared with the adaptive finite element TBC method for both polarizations to illustrate the competitive behavior of the proposed method.
Homogenization in terms of multiscale limits transforms a multiscale problem with $n+1$ asymptotically separated microscales posed on a physical domain $D subset mathbb{R}^d$ into a one-scale problem posed on a product domain of dimension $(n+1)d$ by introducing $n$ so-called fast variables. This procedure allows to convert $n+1$ scales in $d$ physical dimensions into a single-scale structure in $(n+1)d$ dimensions. We prove here that both the original, physical multiscale problem and the corresponding high-dimensional, one-scale limiting problem can be efficiently treated numerically with the recently developed quantized tensor-train finite-element method (QTT-FEM). The method is based on restricting computation to sequences of nested subspaces of low dimensions (which are called tensor ranks) within a vast but generic virtual (background) discretization space. In the course of computation, these subspaces are computed iteratively and data-adaptively at runtime, bypassing any offline precomputation. For the purpose of theoretical analysis, such low-dimensional subspaces are constructed analytically to bound the tensor ranks vs. error $tau>0$. We consider a model linear elliptic multiscale problem in several physical dimensions and show, theoretically and experimentally, that both (i) the solution of the associated high-dimensional one-scale problem and (ii) the corresponding approximation to the solution of the multiscale problem admit efficient approximation by the QTT-FEM. These problems can therefore be numerically solved in a scale-robust fashion by standard (low-order) PDE discretizations combined with state-of-the-art general-purpose solvers for tensor-structured linear systems. We prove scale-robust exponential convergence, i.e., that QTT-FEM achieves accuracy $tau$ with the number of effective degrees of freedom scaling polynomially in $log tau$.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
