On the Projections of Multivalued Maps

This paper considers analogues of the Helmholtz projections of the set of selections of a piecewise smooth multivalued map , n2. It is shown that, for mn–1 (m=1), the closure of the projection of on the subspace of gradient fields (solenoidal vector fields) is a convex set. For the general case, there are given point-wise conditions on the values of the map which ensure that the closure of the projection of contains the zero element. Possible applications to optimal control problems are discussed.     

A numerical scheme based on mean value solutions for the helmholtz equation on triangular grids

A numerical treatment for the Dirichlet boundary value problem on regular triangular grids for homogeneous Helmholtz equations is presented, which also applies to the convection-diffusion problems. The main characteristic of the method is that an accuracy estimate is provided in analytical form with a better evaluation than that obtained with the usual finite difference method. Besides, this classical method can be seen as a truncated series approximation to the proposed method. The method is developed from the analytical solutions for the Dirichlet problem on a ball together with an error evaluation of an integral on the corresponding circle, yielding accuracy. Some numerical examples are discussed and the results are compared with other methods, with a consistent advantage to the solution obtained here.

     

求解二维多区域声波散射问题的边界元方法

对于多散射区域的声波散射问题的外Neumann边值问题,用单层位势来逼近每个散射域上的散射波,再利用位势理论的跳跃关系将问题转换为第二类边界积分方程组的求解问题,然后用Nystrom方法进行了求解.对多个随机散射区域的声波散射问题,数值例子体现了该求解方法的可行性和准确性.     

Finite element method to fluid–solid interaction problems with unbounded periodic interfaces

Consider a time‐harmonic acoustic plane wave incident onto a doubly periodic (biperiodic) surface from above. The medium above the surface is supposed to be filled with a homogeneous compressible inviscid fluid of constant mass density, whereas the region below is occupied by an isotropic and linearly elastic solid body characterized by its Lamé constants. This article is concerned with a variational approach to the fluid–solid interaction problems with unbounded biperiodic Lipschitz interfaces between the domains of the acoustic and elastic waves. The existence of quasiperiodic solutions in Sobolev spaces is established at arbitrary frequency of incidence, while uniqueness is proved only for small frequencies or for all frequencies excluding a discrete set. A finite element scheme coupled with Dirichlet‐to‐Neumann mappings is proposed and the convergence analysis is performed. The Dirichlet‐to‐Neumann mappings are approximated by truncated Rayleigh series expansions. Finally, numerical tests in 2D are presented to confirm the convergence of solutions and the energy balance formula. In particular, the frequency spectrum of normally reflected signals is plotted for water–brass and water–brass–water interfaces. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 5–35, 2016     

Analysis of a finite PML approximation for the three dimensional time-harmonic Maxwell and acoustic scattering problems

We consider the approximation of the frequency domain three-dimensional Maxwell scattering problem using a truncated domain perfectly matched layer (PML). We also treat the time-harmonic PML approximation to the acoustic scattering problem. Following work of Lassas and Somersalo in 1998, a transitional layer based on spherical geometry is defined, which results in a constant coefficient problem outside the transition. A truncated (computational) domain is then defined, which covers the transition region. The truncated domain need only have a minimally smooth outer boundary (e.g., Lipschitz continuous). We consider the truncated PML problem which results when a perfectly conducting boundary condition is imposed on the outer boundary of the truncated domain. The existence and uniqueness of solutions to the truncated PML problem will be shown provided that the truncated domain is sufficiently large, e.g., contains a sphere of radius . We also show exponential (in the parameter ) convergence of the truncated PML solution to the solution of the original scattering problem inside the transition layer.

Our results are important in that they are the first to show that the truncated PML problem can be posed on a domain with nonsmooth outer boundary. This allows the use of approximation based on polygonal meshes. In addition, even though the transition coefficients depend on spherical geometry, they can be made arbitrarily smooth and hence the resulting problems are amenable to numerical quadrature. Approximation schemes based on our analysis are the focus of future research.

     

Acoustic absorbers at low frequency based on split-tube metamaterials

The remarkable properties of acoustic metamaterials have attracted massive researches and applications, especially on low-frequency sound absorptions. Currently, most of the acoustic metamaterial absorbers employ resonances in plastic cavities, and their structural strengths are important in many circumstances, especially in harsh environment. However, studies of metamaterials including this point are very scarce. Here, we propose an acoustic metamaterial for low-frequency (<500 Hz) absorptions, composed of three nested square split tubes with inverted opening directions. The efficiency of the absorber is investigated both numerically and experimentally, and absorptions at the peeks are found to exceed 90% and the frequency can be effectively adjusted by tuning its geometric parameters. We further test its yield strength under compression and confirm its buckling behavior happens from the outmost layer. This tunable acoustic metamaterial with a fairly good mechanical strength may lead to broad applications in noise reduction.     

A hybrid numerical–asymptotic boundary element method for high frequency scattering by penetrable convex polygons

We present a novel hybrid numerical–asymptotic boundary element method for high frequency acoustic and electromagnetic scattering by penetrable (dielectric) convex polygons. Our method is based on a standard reformulation of the associated transmission boundary value problem as a direct boundary integral equation for the unknown Cauchy data, but with a nonstandard numerical discretization which efficiently captures the high frequency oscillatory behaviour. The Cauchy data is represented as a sum of the classical geometrical optics approximation, computed by a beam tracing algorithm, plus a contribution due to diffraction, computed by a Galerkin boundary element method using oscillatory basis functions chosen according to the principles of the Geometrical Theory of Diffraction. We demonstrate with a range of numerical experiments that our boundary element method can achieve a fixed accuracy of approximation using only a relatively small, frequency-independent number of degrees of freedom. Moreover, for the scattering scenarios we consider, the inclusion of the diffraction term provides an order of magnitude improvement in accuracy over the geometrical optics approximation alone.     

Symmetric boundary integral formulations for Helmholtz transmission problems

In this work we propose and analyze numerical methods for the approximation of the solution of Helmholtz transmission problems in two or three dimensions. This kind of problems arises in many applications related to scattering of acoustic, thermal and electromagnetic waves. Formulations based on boundary integral methods are powerful tools to deal with transmission problems in unbounded media. Different formulations using boundary integral equations can be found in the literature. We propose here new symmetric formulations based on a paper by Martin Costabel and Ernst P. Stephan (1985), that uses the Calderón projector for the interior and exterior problems to develop closed expressions for the interior and exterior Dirichlet-to-Neumann operators. These operators are then matched to obtain an integral system that is equivalent to the Helmholtz transmission problem and uses Cauchy data on the transmission boundary as unknowns. We show how to simplify the aspect and analysis of the method by employing an additional mortar unknown with respect to the ones used in the original paper, writing it in an appropriate way to devise Krylov type iterations based on the separate Dirichlet-to-Neumann operators.     

An alternating method for Cauchy problems for Helmholtz-type operators in non-homogeneous medium

Kozlov & Maz'ya (1989, Algebra Anal., 1, 144–170)proposed an alternating iterative method for solving Cauchyproblems for general strongly elliptic and formally self-adjointsystems. However, in many applied problems, operators appearthat do not satisfy these requirements, e.g. Helmholtz-typeoperators. Therefore, in this study, an alternating procedurefor solving Cauchy problems for self-adjoint non-coercive ellipticoperators of second order is presented. A convergence proofof this procedure is given.     

Modified regularization method for the Cauchy problem of the Helmholtz equation

In this paper, the Cauchy problem for the Helmholtz equation is investigated. It is known that such problem is severely ill-posed. We propose a modified regularization method to solve it based on the solution given by the method of separation of variables. Convergence estimates are presented under two different a-priori bounded assumptions for the exact solution. Finally, numerical examples are given to show the effectiveness of the proposed numerical method.