No abstract

In this paper, we present a numerical method based on the coupling between a Curved Virtual Element Method (CVEM) and a Boundary Element Method (BEM) for the simulation of wave fields scattered by obstacles immersed in homogeneous infinite media. In particular, we consider the 2D time-domain damped wave equation, endowed with a Dirichlet condition on the boundary (sound-soft scattering). To reduce the infinite domain to a finite computational one, we introduce an artificial boundary on which we impose a Boundary Integral Non-Reflecting Boundary Condition (BI-NRBC). We apply a CVEM combined with the Crank–Nicolson time integrator in the interior domain, and we discretize the BI-NRBC by a convolution quadrature formula in time and a collocation method in space. We present some numerical results to test the performance of the proposed approach and to highlight its effectiveness, especially when obstacles with complex geometries are considered.

For the solution of 2D exterior Dirichlet Poisson problems, we propose the coupling of a Curved Virtual Element Method (CVEM) with a Boundary Element Method (BEM), by using decoupled approximation orders. We provide optimal convergence error estimates, in the energy and in the weaker $$\textit{L}^\text {2}$$ L 2 -norm, in which the CVEM and BEM contributions to the error are separated. This allows for taking advantage of the high order flexibility of the CVEM to retrieve an accurate discrete solution by using a low order BEM. The numerical results confirm the a priori estimates and show the effectiveness of the proposed approach.

No abstract

We consider the non-symmetric coupling of finite and boundary elements to solve second-order nonlinear partial differential equations defined in unbounded domains. We present a novel condition that ensures that the associated semi-linear form induces a strongly monotone operator, keeping track of the dependence on the linear combination of the interior domain equation with the boundary integral one. We show that an optimal ellipticity condition, relating the nonlinear operator to the contraction constant of the shifted double-layer integral operator, is guaranteed by choosing a particular linear combination. These results generalize those obtained by Of and Steinbach [Is the one-equation coupling of finite and boundary element methods always stable?, ZAMM Z. Angew. Math. Mech. 93 (2013), 6–7, 476–484] and [On the ellipticity of coupled finite element and one-equation boundary element methods for boundary value problems, Numer. Math. 127 (2014), 3, 567–593], and by Steinbach [A note on the stable one-equation coupling of finite and boundary elements, SIAM J. Numer. Anal. 49 (2011), 4, 1521–1531], where the simple sum of the two coupling equations has been considered. Numerical examples confirm the theoretical results on the sharpness of the presented estimates.