An a priori error analysis for the coupling of local discontinuous Galerkin and boundary element methods

In this paper we analyze the coupling of local discontinuous Galerkin (LDG) and boundary element methods as applied to linear exterior boundary value problems in the plane. As a model problem we consider a Poisson equation in an annular polygonal domain coupled with a Laplace equation in the surrounding unbounded exterior region. The technique resembles the usual coupling of finite elements and boundary elements, but the corresponding analysis becomes quite different. In particular, in order to deal with the weak continuity of the traces at the interface boundary, we need to define a mortar-type auxiliary unknown representing an interior approximation of the normal derivative. We prove the stability of the resulting discrete scheme with respect to a mesh-dependent norm and derive a Strang-type estimate for the associated error. Finally, we apply local and global approximation properties of the subspaces involved to obtain the a priori error estimate in the energy norm.

     

An acoustic scattering problem for periodic,inhomogeneous media

We study the problem of the scattering by a periodic, inhomogeneous, penetrable medium. Using the Dirichlet-to-Neumann operator from the classical formulation of the problem we derive a variational equation and give regularity result to show the equivalence of both formulations. We present certain uniqueness results, which by the Fredholm alternative yield existence of the solution and its continuous dependence on the incoming wave. We prove existence of a solution for special incident waves even if there is no uniqueness. A result about analytical dependence of the solution on the wave number and the incident angle is given. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

An improved On-Surface Radiation Condition for acoustic scattering problems in the high-frequency spectrum

This Note addresses the derivation of an improved On-Surface Radiation Condition for the numerical solution of the exterior Helmholtz equation at high-frequencies. This condition is built as an approximation of the Neumann-to-Dirichlet map by using a local regularization of its principal classical symbol in the gliding zone for modelling the creeping waves. The numerical simulation of this pseudodifferential operator is efficiently realized with a linear cost according to the dimension of the boundary element approximation space using suitable complex Padé approximants. A numerical example is provided. To cite this article: X. Antoine et al., C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

On inf-sup conditions and projections in the theory of mixed finite-element methods

We indicate constraints on the space of finite elements providing the validity of discrete inf-sup conditions and the existence of projections specific for mixed finite-element methods. We consider both conformal and nonconformal approximations. We suggest a definition of special projections onto the vector space of finite elements which provides their existence under quite general conditions without determining the degrees of freedom of the elements.     

Stable methods for an inverse problem in acoustic scattering by an obstacle and an inhomogeneous medium

We consider (in two-dimensional Euclidean space) the scattering of a plane, time-harmonic acoustic wave by an inhomogeneous medium Ω with compact support and a bounded obstacle D lying completely outside of the inhomogeneous medium. We show that one may determine the shape of D and the local speed of sound in Ω from a knowledge of the asymptotic behavior of the scattered wave (i.e. the far field). This is done by considering a constrained optimization problem and employing integral equation and conformal mapping techniques. By assuming a priori that the functions which determine the shape of D and the local speed of sound in Ω lie in given compact sets, we show that the problem is stable, in the sense that the solution of the inverse scattering problem depends continuously on the far field data.     

Low frequency acoustic and electromagnetic scattering

This paper deals with two classes of problems arising from acoustics and electromagnetic scattering in the low frequency situations. The first class of problem involves solving the Helmholtz equation with Dirichlet boundary conditions on an arbitrary two-dimensional body whereas the second one is an interior-exterior interface problem with the Helmholtz equation in the exterior. Low frequency analysis shows that there are two intermediate problems shich solve the above problems to accuracy O(k²(logk)²) where k is the wave number. These solutions are more accurate than the simpler zero frequency approximations but require very little more work to compute. For the Dirichlet problem numerical examples are shown to verify our theoretical estimates.     

Stability of penalty finite-element methods for nonconforming problems

Penalty methods have been proposed as a viable method for enforcing interelement continuity constraints on nonconforming elements. Particularly for fourth-order problems in which C¹-continuity leads to elements of high degree or complex composite elements, the use of penalty methods to enforce the C¹-continuity constraint appears promising. In this study we demonstrate equivalence of the finite-element penalty method to a hybrid method and provide a stability analysis which implies that the penalty method is stable only if reduced integration of a certain order is used. Numerical experiments confirm that the penalty method fails if this condition is not met.     

On regularized full- and partial-cloaks in acoustic scattering

The aim of this work is to derive sharp quantitative estimates of the qualitative convergence results developed by Li et al. (2015) for regularized full- and partial-cloaks through transformation optics approach. Let Γ₀ be a compact set in ?³ and Γ_δ be a δ-neighborhood of Γ₀ for δ∈?₊. Γ_δ represents the virtual domain used for the blow-up construction. By incorporating suitably designed lossy layers, it is shown that if the generating set Γ₀ is a generic curve, then one would have an approximate full-cloak within δ² to the perfect full-cloak, whereas if Γ₀ is the closure of an open subset on a flat surface, then one would have an approximate partial-cloak within δ to its perfect counterpart. The estimates derived are independent of the contents being cloaked; that is, the cloaking devices are capable of nearly cloaking an arbitrary content. Furthermore, as a significant by-product, our argument allows the relaxation of the convexity requirement on Γ₀ by Li et al. (2015), which is critical for the Mosco convergence argument therein.     

Time harmonic acoustic scattering in anisotropic media

The scattering problem of a plane or a point source generated wave is considered for the case where both the medium of propagation and the interior of the scatterer exhibit their own anisotropies. A particular redirected gradient operator is introduced, which carries all directional characteristics of the anisotropic medium. Once the fundamental solution is obtained, integral representations for the scattered as well as for the interior and the total fields are generated. For such media even the handling of the singularities, in generating integral representations, depends on the characteristics of the particular medium. A modified, also medium dependent, radiation condition is introduced. Detailed asymptotic analysis leads to an integral representation for the scattering amplitude. The associated energy functionals are presented and the relative cross sections are also defined. Copyright © 2005 John Wiley & Sons, Ltd.     

Approximate solution methods in mechanics analysis of a class of nonlinear seepage problems by a finite-element method

We analyze the nonlinear boundary-value problem of seepage under a subsurface hydrotechnical construction over an inclined rectilinear aquifer. The method of inverse boundary-value problems is applied, using the velocity hodograph plane in which the original problem is reduced to a linear problem. The linear problem is solved in the general case using the finite-element method. A computer program realizing the proposed algorithms has been developed. We have used this program to run a series of numerical experiments, reaching certain conclusions about the behavior of the main seepage characteristics.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 55, pp. 75–80, 1985.     

Large deformation shape uncertainty quantification in acoustic scattering

We address shape uncertainty quantification for the two-dimensional Helmholtz transmission problem, where the shape of the scatterer is the only source of uncertainty. In the framework of the so-called deterministic approach, we provide a high-dimensional parametrization for the interface. Each domain configuration is mapped to a nominal configuration, obtaining a problem on a fixed domain with stochastic coefficients. To compute surrogate models and statistics of quantities of interest, we apply an adaptive, anisotropic Smolyak algorithm, which allows to attain high convergence rates that are independent of the number of dimensions activated in the parameter space. We also develop a regularity theory with respect to the spatial variable, with norm bounds that are independent of the parametric dimension. The techniques and theory presented in this paper can be easily generalized to any elliptic problem on a stochastic domain.     

Mathematical simulation of acoustic wave scattering in fractured media

The simulation of acoustic waves in fractured media is considered. A self-consistent field model is proposed that describes the formation of a scattered field and the attenuation of the incident field. For the total field, a wave equation with a complex velocity is derived and the corresponding dispersion equation is studied. A frequency-dependent field damping law and an energy variation law are established. An initial and a boundary value problem for waves in a fractured medium is addressed. A finite-difference scheme for the initial value problem is constructed, and a condition for its stability is established. Numerical results are presented.     

Error analysis for hybrid Trefftz methods coupling traction conditions in linear elastostatics

For linear elastostatics, the Lagrange multiplier to couple the displacement (i.e., Dirichlet) condition is well known in mathematics community, but the Lagrange multiplier to couple the traction (i.e., Neumann) condition is popular for elasticity problems by the Trefftz method in engineering community, which is called the Hybrid Trefftz method (HTM). However, there has not been any analysis for these Lagrange multipliers to couple the traction condition so far. New error analysis of the HTM for elasticity problems is explored in this paper, to derive error bounds with the optimal convergence rates. Numerical experiments are reported to support this analysis. The error analysis of the HTM for linear elastostatics is the main aim of this paper. In this paper, the collocation Trefftz method (CTM) without a multiplier is also introduced, accompanied with error analysis. Numerical comparisons are made for HTM and CTM using fundamental solutions (FS) and particular solutions (PS). The error analysis and numerical computations show that the accuracy of the HTM is equivalent to that of the CTM, but the stability of the CTM is good. For elasticity and other complicated problems, the simplicity of algorithms and programming grants the CTM a remarkable advantage. More numerical comparisons show that using PS is more efficient than using FS in both HTM and CTM. However, since the optimal convergence rates are the most important criterion in evaluation of numerical methods, the global performance of the HTM is as good as that of the CTM. The comparisons of HTM and CTM using FS and PS are the next aim of this article. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

用正则化方法求解声波散射反问题

研究了从声波散射场的远场模式的信息来再现散射物边界形状的反问题.首先构造表达散射物特征的指示函数,然后利用该函数之特性,建立求解该类反问题的基本方程,从而确定散射物的边界形状.在这个算法中,不需预先知道散射物的边界类型和形状等知识,从T ikhonov正则化方法进行的数值计算结果表明了该方法是有效的和实用的.     

A conforming mixed finite-element method for the coupling of fluid flow with porous media flow

Salim Meddahi ^{We consider a porous medium entirely enclosed within a fluidregion and present a well-posed conforming mixed finite-elementmethod for the corresponding coupled problem. The interfaceconditions refer to mass conservation, balance of normal forcesand the Beavers–Joseph–Saffman law, which yieldsthe introduction of the trace of the porous medium pressureas a suitable Lagrange multiplier. The finite-element subspacesdefining the discrete formulation employ Bernardi–Raugeland Raviart–Thomas elements for the velocities, piecewiseconstants for the pressures and continuous piecewise-linearelements for the Lagrange multiplier. We show stability, convergenceand a priori error estimates for the associated Galerkin scheme.Finally, we provide several numerical results illustrating thegood performance of the method and confirming the theoreticalrates of convergence.}