اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدA Galerkin BEM for high-frequency scattering problems based on frequency dependent changes of variables
134
0
0.0
(
0
)
تأليف
Fatih Ecevit
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
In this paper we develop a class of efficient Galerkin boundary element methods for the solution of two-dimensional exterior single-scattering problems. Our approach is based upon construction of Galerkin approximation spaces confined to the asymptotic behaviour of the solution through a certain direct sum of appropriate function spaces weighted by the oscillations in the incident field of radiation. Specifically, the function spaces in the illuminated/shadow regions and the shadow boundaries are simply algebraic polynomials whereas those in the transition regions are generated utilizing novel, yet simple, emph{frequency dependent changes of variables perfectly matched with the boundary layers of the amplitude} in these regions. While, on the one hand, we rigorously verify for smooth convex obstacles that these methods require only an $mathcal{O}left( k^{epsilon} right)$ increase in the number of degrees of freedom to maintain any given accuracy independent of frequency, and on the other hand, remaining in the realm of smooth obstacles they are applicable in more general single-scattering configurations. The most distinctive property of our algorithms is their emph{remarkable success} in approximating the solution in the shadow region when compared with the algorithms available in the literature.
قيم البحث
اقرأ أيضاً
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.
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 ty
pes 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.
We consider quadrature formulas of high order in time based on Radau-type, L-stable implicit Runge-Kutta schemes to solve time dependent stiff PDEs. Instead of solving a large nonlinear system of equations, we develop a method that performs iterative
deferred corrections to compute the solution at the collocation nodes of the quadrature formulas. The numerical stability is guaranteed by a dedicated operator splitting technique that efficiently handles the stiffness of the PDEs and provides initial and intermediate solutions to the iterative scheme. In this way the low order approximations computed by a tailored splitting solver of low algorithmic complexity are iteratively corrected to obtain a high order solution based on a quadrature formula. The mathematical analysis of the numerical errors and local order of the method is carried out in a finite dimensional framework for a general semi-discrete problem, and a time-stepping strategy is conceived to control numerical errors related to the time integration. Numerical evidence confirms the theoretical findings and assesses the performance of the method in the case of a stiff reaction-diffusion equation.
When using spectral methods, a question arises as how to determine the expansion order, especially for time-dependent problems in which emerging oscillations may require adjusting the expansion order. In this paper, we propose a frequency-dependent $
p$-adaptive technique that adaptively adjusts the expansion order based on a frequency indicator. Using this $p$-adaptive technique, combined with recently proposed scaling and moving techniques, we are able to devise an adaptive spectral method in unbounded domains that can capture and handle diffusion, advection, and oscillations. As an application, we use this adaptive spectral method to numerically solve the Schr{o}dinger equation in the whole domain and successfully capture the solutions oscillatory behavior at infinity.
This paper presents a new fast multipole boundary element method (FM-BEM) for solving the acoustic transmission problems in 2D periodic media. We divide the periodic media into many fundamental blocks, and then construct the boundary integral equatio
ns in the fundamental block. The fast multipole algorithm is proposed for the square and hexagon periodic systems, the convergence of the algorithm is analyzed. We then apply the proposed method to the acoustic transmission problems for liquid phononic crystals and derive the acoustic band gaps of the phononic crystals. By comparing the results with those from plane wave expansion method, we conclude that our method is efficient and accurate.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
