A discontinuous enrichment method for capturing evanescent waves in multiscale fluid and fluid/solid problems

Abstract An evanescent wave is produced when a propagating incident wave impinges on an interface between two media or materials at a sub-critical angle. The sub-scale nature of such a wave makes it difficult to be captured computationally. In this paper, the Discontinuous Enrichment Method (DEM) developed in [C. Farhat, I. Harari, L.P. Franca, The discontinuous enrichment method, Comput. Methods Appl. Mech. Engrg. 190 (2001) 6455–6479; C. Farhat, I. Harari, U. Hetmaniuk, A discontinuous Galerkin method with Lagrange multipliers for the solution of Helmholtz problems in the mid-frequency regime, Comput. Methods Appl. Mech. Engrg. 192 (2003) 1389–1419; C. Farhat, I. Harari, U. Hetmaniuk, The discontinuous enrichment method for multiscale analysis, Comput. Methods Appl. Mech. Engrg. 192 (2003) 3195–3210; C. Farhat, R. Tezaur, P. Weidemann-Goiran, Higher-order extensions of a discontinuous Galerkin method for mid-frequency Helmholtz problems, Int. J. Numer. Methods Engrg. 61 (2004) 1938–1956; C. Farhat, P. Wiedemann-Goiran, R. Tezaur, A discontinuous Galerkin method with plane waves and Lagrange multipliers for the solution of short wave exterior Helmholtz problems on unstructured meshes. New computational methods for wave propagation, Wave Motion 39(4) (2004) 307–317; L. Zhang, R. Tezaur, C. Farhat, The discontinuous enrichment method for elastic wave propagation in the medium-frequency regime, Int. J. Numer. Methods Engrg. 66 (2006) 2086–2114; R. Tezaur, C. Farhat, Three-dimensional discontinuous Galerkin elements with plane waves and Lagrange multipliers for the solution of mid-frequency Helmholtz problems, Int. J. Numer. Methods Engrg. 66 (2006) 796–815] is extended to a class of evanescent wave problems. New DEM elements are constructed by enriching polynomial approximations with free-space evanescent solutions in order to achieve high accuracy for problems with fluid/fluid or fluid/solid interfaces on which evanescent waves may occur. The new DEM elements are applied to the solution of various two-dimensional model problems with evanescent waves. Their performance is observed to be better than that of the basic Helmholtz type DEM elements, and superior to that of the classical higher-order polynomial finite element method.

[1]  I. Babuska,et al.  The generalized finite element method , 2001 .

[2]  Ted Belytschko,et al.  Structured extended finite element methods for solids defined by implicit surfaces , 2002 .

[3]  K. G. Budden.  Wave Propagation , 1969, Nature.

[4]  J. Tinsley Oden,et al.  An hp Adaptive Method Using Clouds C , 2006 .

[5]  Charbel Farhat,et al.  A discontinuous Galerkin method with plane waves and Lagrange multipliers for the solution of short wave exterior Helmholtz problems on unstructured meshes , 2004 .

[6]  D. S. Drumheller,et al.  Introduction to elastic wave propagation , 1994 .

[7]  Ivo Babuška,et al.  A posteriori error estimation for generalized finite element methods , 2006 .

[8]  Pierre Ladevèze,et al.  Calculation of medium-frequency vibrations over a wide frequency range , 2005 .

[9]  C. Farhat,et al.  The Discontinuous Enrichment Method , 2000 .

[10]  F. Ihlenburg.  Finite Element Analysis of Acoustic Scattering , 1998 .

[11]  M. Hong,et al.  Evanescent Wave Vibrational Microscopy , 2004 .

[12]  Charbel Farhat,et al.  Higher‐order extensions of a discontinuous Galerkin method for mid‐frequency Helmholtz problems , 2004 .

[13]  Charbel Farhat,et al.  A discontinuous Galerkin method with Lagrange multipliers for the solution of Helmholtz problems in the mid-frequency regime , 2003 .

[14]  Peter Monk,et al.  A least-squares method for the Helmholtz equation , 1999 .

[15]  C. Farhat,et al.  The discontinuous enrichment method for elastic wave propagation in the medium‐frequency regime , 2006 .

[16]  Jari P. Kaipio,et al.  The Ultra-Weak Variational Formulation for Elastic Wave Problems , 2004, SIAM J. Sci. Comput..

[17]  Harry J Simpson,et al.  Laboratory measurements of sound scattering from a buried sphere above and below the critical angle. , 2003, The Journal of the Acoustical Society of America.

[18]  James Shipman,et al.  Wave Motion , 2006 .

[19]  Charbel Farhat,et al.  The discontinuous enrichment method for multiscale analysis , 2003 .

[20]  Maurice Roseau,et al.  Asymptotic Wave Theory , 1975 .

[21]  L. Thompson,et al.  COMPLEX WAVE-NUMBER DISPERSION ANALYSIS OF STABILIZED FINITE ELEMENT METHODS FOR ACOUSTIC FLUID – STRUCTURE INTERACTION , 2007 .

[22]  O. C. Zienkiewicz,et al.  A new cloud-based hp finite element method , 1998 .

[23]  R. A. Rhodes,et al.  Propagation of waves , 1969 .

[24]  John A. DeSanto,et al.  Scalar wave theory , 1992 .

[25]  T. Belytschko,et al.  Arbitrary branched and intersecting cracks with the eXtended Finite Element Method , 2000 .

[26]  Ivo Babuška,et al.  p‐version of the generalized FEM using mesh‐based handbooks with applications to multiscale problems , 2004 .

[27]  Frédérique de Fornel,et al.  Evanescent waves : from Newtonian optics to atomic optics , 2001 .

[28]  Charbel Farhat,et al.  Three‐dimensional discontinuous Galerkin elements with plane waves and Lagrange multipliers for the solution of mid‐frequency Helmholtz problems , 2006 .