Abstract:We formulate a new family of high order on-surface radiation conditions to approximate the outgoing solution to the Helmholtz equation in exterior domains. Motivated by the pseudo-differential expansion of the Dirichlet-to-Neumann operator developed by Antoine et al. (J. Math. Anal. Appl. 229:184-211, 1999), we design a systematic procedure to apply pseudo-differential symbols of arbitrarily high order. Numerical results are presented to illustrate the performance of the proposed method for solving both the D… Show more
“…Local boundary conditions can be derived by approximating the square root in the symbol of the DtN operator (1). In their seminal paper, Engquist and Majda [31] derived a family of local boundary conditions by using a Padé approximation of the square root.…”
Section: Rational Approximation Of the Square Rootmentioning
This paper deals with the design and validation of accurate local absorbing boundary conditions set on convex polygonal and polyhedral computational domains for the finite element solution of highfrequency acoustic scattering problems. While high-order absorbing boundary conditions (HABCs) are accurate for smooth fictitious boundaries, the precision of the solution drops in the presence of corners if no specific treatment is applied. We present and analyze two strategies to preserve the accuracy of Padé-type HABCs at corners: first by using compatibility relations (derived for right angle corners) and second by regularizing the boundary at the corner. Exhaustive numerical results for two-and three-dimensional problems are reported in the paper. They show that using the compatibility relations is optimal for domains with right angles. For the other cases, the error still remains acceptable, but depends on the choice of the corner treatment according to the angle.
“…Local boundary conditions can be derived by approximating the square root in the symbol of the DtN operator (1). In their seminal paper, Engquist and Majda [31] derived a family of local boundary conditions by using a Padé approximation of the square root.…”
Section: Rational Approximation Of the Square Rootmentioning
This paper deals with the design and validation of accurate local absorbing boundary conditions set on convex polygonal and polyhedral computational domains for the finite element solution of highfrequency acoustic scattering problems. While high-order absorbing boundary conditions (HABCs) are accurate for smooth fictitious boundaries, the precision of the solution drops in the presence of corners if no specific treatment is applied. We present and analyze two strategies to preserve the accuracy of Padé-type HABCs at corners: first by using compatibility relations (derived for right angle corners) and second by regularizing the boundary at the corner. Exhaustive numerical results for two-and three-dimensional problems are reported in the paper. They show that using the compatibility relations is optimal for domains with right angles. For the other cases, the error still remains acceptable, but depends on the choice of the corner treatment according to the angle.
“…The subject of nonreflecting or absorbing boundary conditions is beyond the scope of this paper. We refer to [9,13,6,20,1] for some relevant articles on that topic. We consider…”
In the mathematical analysis of photoacoustic imaging, it is usually assumed that the acoustic pressure (Dirichlet data) is measured on a detection surface. However, actual ultrasound detectors gather data of a different type. In this paper, we propose a more realistic mathematical model of ultrasound measurements acquired by the Fabry-Perot sensor. This modeling incorporates directional response of such sensors. We study the solvability of the resulting photoacoustic tomography problem, concluding that the problem is well-posed under certain assumptions. Numerical reconstructions are implemented using the Landweber iterations, after discretization of the governing equations using the finite element method.
“…Figure 1 displays (a) the sparsity pattern of the matrix A 1,2 for the configuration of two spheres from subsection 6.1, (b) the sparsity pattern of the discrete Laplace-Beltrami operator ∆ T 1 , and (c) CPU time for the application of the propagator (30) as a function of DOF. For the numerical results presented in the next section, we apply a fixed-point iteration to solve the system (32). Due to the convexity of the surfaces ∂Ω j for all j = 1, 2, ..., J, the norm of the propagation operator P i, j decreases as the distance between the i th and j th obstacles increases.…”
Section: Complexity O(n)mentioning
confidence: 99%
“…Possible future work includes the extension to electromagnetic and elastic waves that govern important engineering applications. It is also possible to increase the order of approximations of the Dirichlet-to-Neumann map to improve the accuracy of the OSRC [32]. Here we have considered the second order differential approximation (16).…”
Section: Final Remarksmentioning
confidence: 99%
“…Other improvements have been accomplished by Jones [11,12,13], Ammari [14,15], Calvo et al [16,17], Barucq et al [18,19,20], Chaillat et al [21,22] and Darbas et al [2,3,4,5,23]. See also [24,25,26,27,28,29,30,31,25,32].…”
The formulation of the on-surface radiation condition (OSRC) is extended to handle wave scattering problems in the presence of multiple obstacles. The new multiple-OSRC simultaneously accounts for the outgoing behavior of the wave fields, as well as, the multiple wave reflections between the obstacles. Like boundary integral equations (BIE), this method leads to a reduction in dimensionality (from volume to surface) of the discretization region. However, as opposed to BIE, the proposed technique leads to boundary integral equations with smooth kernels. Hence, these Fredholm integral equations can be handled accurately and robustly with standard numerical approaches without the need to remove singularities. Moreover, under weak scattering conditions, this approach renders a convergent iterative method which bypasses the need to solve single scattering problems at each iteration.Inherited from the original OSRC, the proposed multiple-OSRC is generally a crude approximate method. If accuracy is not satisfactory, this approach may serve as a good initial guess or as an inexpensive pre-conditioner for Krylov iterative solutions of BIE.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.