skip to main content

Title: Convergence of an adaptive finite element DtN method for the elastic wave scattering by periodic structures
Consider the scattering of a time-harmonic elastic plane wave by a periodic rigid surface. The elastic wave propagation is governed by the two-dimensional Navier equation. Based on a Dirichlet-to-Neumann (DtN) map, a transparent boundary condition (TBC) is introduced to reduce the scattering problem into a boundary value problem in a bounded domain. By using the finite element method, the discrete problem is considered, where the TBC is replaced by the truncated DtN map. A new duality argument is developed to derive the a posteriori error estimate, which contains both the finite element approximation error and the DtN truncation error. An a posteriori error estimate based adaptive finite element algorithm is developed to solve the elastic surface scattering problem. Numerical experiments are presented to demonstrate the effectiveness of the proposed method.
Award ID(s):
Publication Date:
Journal Name:
Computer methods in applied mechanics and engineering
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We design and analyze a $$C^0$$ C 0 interior penalty method for the approximation of classical solutions of the Dirichlet boundary value problem of the Monge–Ampère equation on convex polygonal domains. The method is based on an enhanced cubic Lagrange finite element that enables the enforcement of the convexity of the approximate solutions. Numerical results that corroborate the a priori and a posteriori error estimates are presented. It is also observed from numerical experiments that this method can capture certain weak solutions.
  2. A column-supporting system for embankments on soft soils is analyzed using the Finite Element Method. The numerical problem has boundary conditions similar to that of a trap-door, but, contrary to the classical trap-door studies, the focus in this paper is on the load transfer to the supporting columns measured by the system efficacy, and on the critical embankment height that prevents differential settlements on the embankment surface. The stress state simulated in the embankment rapidly evolves in the first stages of the soft soil settlement, measured by millimeters, but the changes to the elastic-plastic stress field developed are marginal in the subsequent stages of settlement simulated up to 10 cm. The displacement field in low embankments with a height comparable to the column spacing is dominated by the failure mechanism with shear bands reaching the embankment crown, causing differential settlements on the embankment surface. In higher embankments, the failure mechanism is confined to the lower portion of the fill, and no differential settlements were detected on the embankment surface. It was suggested that a hypothetical diffused soil arch formed above the failure mechanism, with a shape approximately following the principal stress trajectories. This conjecture was made based on the presencemore »of elevated stress along the symmetry plane of the embankment-column periodic cell. Numbers related to system efficacy and critical height are reported for a limited set of model parameters, but the qualitative outcome of the study is likely applicable to a wider variety of embankments.« less
  3. Consider the scattering of a time-harmonic acoustic plane wave by a bounded elastic obstacle which is immersed in a homogeneous acoustic medium. This paper is concerned with an inverse acoustic-elastic interaction problem, which is to determine the location and shape of the elastic obstacle by using either the phased or phaseless far-field data. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The jump relations are studied for the second derivatives of the single-layer potential in order to deduce the corresponding boundary integral equations. The well-posedness is discussed for the solution of the coupled boundary integral equations. An efficient and high order Nyström-type discretization method is proposed for the integral system. A numerical method of nonlinear integral equations is developed for the inverse problem. For the case of phaseless data, we show that the modulus of the far-field pattern is invariant under a translation of the obstacle. To break the translation invariance, an elastic reference ball technique is introduced. We prove that the inverse problem with phaseless far-field pattern has a unique solution under certain conditions. In addition, a numerical method of the reference ball technique based nonlinear integralmore »equations is proposed for the phaseless inverse problem. Numerical experiments are presented to demonstrate the effectiveness and robustness of the proposed methods.« less
  4. Consider the elastic scattering of a time-harmonic wave by multiple well-separated rigid particles with smooth boundaries in two dimensions. Instead of using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert the boundary value problem of the elastic wave equation into a coupled boundary value problem of the Helmholtz equation. Based on single, double, and combined layer potentials with the simpler Green's function of the Helmholtz equation, we present three different boundary integral equations for the coupled boundary value problem. The well-posedness of the new integral equations is established. Computationally, a scattering matrix based method is proposed to evaluate the elastic wave for arbitrarily shaped particles. The method uses the local expansion for the incident wave and the multipole expansion for the scattered wave. The linear system of algebraic equations is solved by GMRES with fast multipole method (FMM) acceleration. Numerical results show that the method is fast and highly accurate for solving elastic scattering problems with multiple particles.
  5. Abstract We derive a residual-based a posteriori error estimator for the conforming hp -Adaptive Finite Element Method ( hp -AFEM) for the steady state Stokes problem describing the slow motion of an incompressible fluid. This error estimator is obtained by extending the idea of a posteriori error estimation for the classical h -version of AFEM. We also establish the reliability and efficiency of the error estimator. The proofs are based on the well-known Clément-type interpolation operator introduced in [28] in the context of the hp -AFEM. Numerical experiments show the performance of an adaptive hp-FEM algorithm using the proposed a posteriori error estimator.