Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAcoustic scattering: high frequency boundary element methods and unified transform methods
220
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We describe some recent advances in the numerical solution of acoustic scattering problems. A major focus of the paper is the efficient solution of high frequency scattering problems via hybrid numerical-asymptotic boundary element methods. We also make connections to the unified transform method due to A.S. Fokas and co-authors, analysing particular instances of this method, proposed by J.A. DeSanto and co-authors, for problems of acoustic scattering by diffraction gratings.
rate research
Read More
We consider a standard elliptic partial differential equation and propose a geometric multigrid algorithm based on Dirichlet-to-Neumann (DtN) maps for hybridized high-order finite element methods. The proposed unified approach is applicable to any locally conservative hybridized finite element method including multinumerics with different hybridized methods in different parts of the domain. For these methods, the linear system involves only the unknowns residing on the mesh skeleton, and constructing intergrid transfer operators is therefore not trivial. The key to our geometric multigrid algorithm is the physics-based energy-preserving intergrid transfer operators which depend only on the fine scale DtN maps. Thanks to these operators, we completely avoid upscaling of parameters and no information regarding subgrid physics is explicitly required on coarse meshes. Moreover, our algorithm is agglomeration-based and can straightforwardly handle unstructured meshes. We perform extensive numerical studies with hybridized mixed methods, hybridized discontinuous Galerkin method, weak Galerkin method, and a hybridized version of interior penalty discontinuous Galerkin methods on a range of elliptic problems including subsurface flow through highly heterogeneous porous media. We compare the performance of different smoothers and analyze the effect of stabilization parameters on the scalability of the multigrid algorithm.
Collocation boundary element methods for integral equations are easier to implement than Galerkin methods because the elements of the discretization matrix are given by lower-dimensional integrals. For that same reason, the matrix assembly also requires fewer computations. However, collocation methods typically yield slower convergence rates and less robustness, compared to Galerkin methods. We explore the extent to which oversampled collocation can improve both robustness and convergence rates. We show that in some cases convergence rates can actually be higher than the corresponding Galerkin method, although this requires oversampling at a faster than linear rate. In most cases of practical interest, oversampling at least lowers the error by a constant factor. This can still be a substantial improvement: we analyze an example where linear oversampling by a constant factor $J$ (leading to a rectangular system of size $JN times N$) improves the error at a cubic rate in the constant $J$. Furthermore, the oversampled collocation method is much less affected by a poor choice of collocation points, as we show how oversampling can lead to guaranteed convergence. Numerical experiments are included for the two-dimensional Helmholtz equation.
In this paper we discuss a hybridised method for FEM-BEM coupling. The coupling from both sides use a Nitsche type approach to couple to the trace variable. This leads to a formulation that is robust and flexible with respect to approximation spaces and can easily be combined as a building block with other hybridised methods. Energy error norm estimates and the convergence of Jacobi iterations are proved and the performance of the method is illustrated on some computational examples.
We study time-harmonic scattering in $mathbb{R}^n$ ($n=2,3$) by a planar screen (a crack in the context of linear elasticity), assumed to be a non-empty bounded relatively open subset $Gamma$ of the hyperplane $mathbb{R}^{n-1}times {0}$, on which impedance (Robin) boundary conditions are imposed. In contrast to previous studies, $Gamma$ can have arbitrarily rough (possibly fractal) boundary. To obtain well-posedness for such $Gamma$ we show how the standard impedance boundary value problem and its associated system of boundary integral equations must be supplemented with additional solution regularity conditions, which hold automatically when $partialGamma$ is smooth. We show that the associated system of boundary integral operators is compactly perturbed coercive in an appropriate function space setting, strengthening previous results. This permits the use of Mosco convergence to prove convergence of boundary element approximations on smoother prefractal screens to the limiting solution on a fractal screen. We present accompanying numerical results, validating our theoretical convergence results, for three-dimensional scattering by a Koch snowflake and a square snowflake.
We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are present. The proposed coupling technique requires minimal changes in the existing schemes while maintaining strict stability, accuracy, and energy conservation. Results are demonstrated on linear and nonlinear scalar conservation laws in two spatial dimensions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
