Explicit time integration with lumped mass matrix for enriched finite elements solution of time domain wave problems

Abstract We present a partition of unity finite element method for wave propagation problems in the time domain using an explicit time integration scheme. Plane wave enrichment functions are introduced at the finite elements nodes which allows for a coarse mesh at low order polynomial shape functions even at high wavenumbers. The initial condition is formulated as a Galerkin approximation in the enriched function space. We also show the possibility of lumping the mass matrix which is approximated as a block diagonal system. The proposed method, with and without lumping, is validated using three test cases and compared to an implicit time integration approach. The stability of the proposed approach against different factors such as the choice of wavenumber for the enrichment functions, the spatial discretization, the distortions in mesh elements or the timestep size, is tested in the numerical studies. The method performance is measured for the solution accuracy and the CPU processing times. The results show significant advantages for the proposed lumping approach which outperforms other considered approaches in terms of stability. Furthermore, the resulting block diagonal system only requires a fraction of the CPU time needed to solve the full system associated with the non-lumped approaches.

[1]  Gyu-Hong Kang,et al.  Thermal and electromagnetic characteristics for cross-sectional design optimization of the integrated production umbilical , 2015, 2015 IEEE Energy Conversion Congress and Exposition (ECCE).

[2]  Jesse Chan,et al.  Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: Explicit time-stepping and efficient mass matrix inversion , 2017, 1708.02972.

[3]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[4]  Marc Alexander Schweitzer,et al.  Variational Mass Lumping in the Partition of Unity Method , 2013, SIAM J. Sci. Comput..

[5]  Jian-Ming Jin,et al.  The Finite Element Method in Electromagnetics , 1993 .

[6]  Ralf Hiptmair,et al.  A Survey of Trefftz Methods for the Helmholtz Equation , 2015, 1506.04521.

[7]  Mohammed Seaïd,et al.  An enriched finite element model with q-refinement for radiative boundary layers in glass cooling , 2014, J. Comput. Phys..

[8]  G.F. Ross,et al.  Time-domain electromagnetics and its applications , 1978, Proceedings of the IEEE.

[9]  O. Laghrouche,et al.  Implementation and computational aspects of a 3D elastic wave modelling by PUFEM , 2017 .

[10]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[11]  Thomas de Quincey.  [C] , 2000, The Works of Thomas De Quincey, Vol. 1: Writings, 1799–1820.

[12]  Mohammed Seaïd,et al.  Time-independent hybrid enrichment for finite element solution of transient conduction-radiation in diffusive grey media , 2013, J. Comput. Phys..

[13]  M. S. Mohamed,et al.  Identifying the wavenumber for the inverse Helmholtz problem using an enriched finite element formulation , 2018, Computer Methods in Applied Mechanics and Engineering.

[14]  C. Schwab,et al.  Boundary element methods for Maxwell's equations on non-smooth domains , 2002, Numerische Mathematik.

[15]  Klaus-Jürgen Bathe,et al.  A finite element method enriched for wave propagation problems , 2012 .

[16]  Shaofan Li,et al.  Reproducing kernel hierarchical partition of unity, Part II—applications , 1999 .

[17]  Ilaria Perugia,et al.  A priori error analysis of space–time Trefftz discontinuous Galerkin methods for wave problems , 2015, 1501.05253.

[18]  Michael I. Mishchenko,et al.  Numerically exact computer simulations of light scattering by densely packed, random particulate media , 2011 .

[19]  Jean-Fu Kiang,et al.  Novel technologies for microwave and millimeter-wave applications , 2004 .

[20]  David Jones,et al.  Coupled electromagnetic and heat transfer model for microwave heating in domestic ovens , 2012 .

[21]  Jon Trevelyan,et al.  High-order finite elements for the solution of Helmholtz problems , 2017 .

[22]  Omar Laghrouche,et al.  Numerical modelling of elastic wave scattering in frequency domain by the partition of unity finite element method , 2009 .

[23]  M. S. Mohamed,et al.  Mixed enrichment for the finite element method in heterogeneous media , 2015 .

[24]  Peter Monk,et al.  Discretization of the Wave Equation Using Continuous Elements in Time and a Hybridizable Discontinuous Galerkin Method in Space , 2014, J. Sci. Comput..

[25]  Omar Laghrouche,et al.  Some numerical aspects of the PUFEM for efficient solution of 2D Helmholtz problems , 2010 .

[26]  Gary Cohen.  Higher-Order Numerical Methods for Transient Wave Equations , 2001 .

[27]  P. Alam.  ‘O’ , 2021, Composites Engineering: An A–Z Guide.

[28]  D. Pozar.  Microwave Engineering , 1990 .

[29]  Marc Alexander Schweitzer,et al.  Stable enrichment and local preconditioning in the particle-partition of unity method , 2011, Numerische Mathematik.

[30]  Shaofan Li,et al.  Reproducing Kernel Hierarchical Partition of Unity , 1998 .

[31]  Peter Monk,et al.  The Ultra Weak Variational Formulation Using Bessel Basis Functions , 2012 .

[32]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[33]  David B. Davidson,et al.  Computational Electromagnetics for RF and Microwave Engineering , 2005 .

[34]  Charbel Farhat,et al.  A space–time discontinuous Galerkin method for the solution of the wave equation in the time domain , 2009 .

[35]  M. S. Mohamed,et al.  Locally enriched finite elements for the Helmholtz equation in two dimensions , 2010 .

[36]  Zhao-Liang Li,et al.  Scale Issues in Remote Sensing: A Review on Analysis, Processing and Modeling , 2009, Sensors.

[37]  R. Hiptmair,et al.  Galerkin Boundary Element Methods for Electromagnetic Scattering , 2003 .

[38]  Omar Laghrouche,et al.  Improvement of PUFEM for the numerical solution of high‐frequency elastic wave scattering on unstructured triangular mesh grids , 2010 .

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

[40]  H. Zheng,et al.  A rigorous and unified mass lumping scheme for higher-order elements , 2017 .

[41]  Ralf Hiptmair,et al.  Boundary Element Methods for Maxwell Transmission Problems in Lipschitz Domains , 2003, Numerische Mathematik.

[42]  O. Laghrouche,et al.  A partition of unity FEM for time‐dependent diffusion problems using multiple enrichment functions , 2013 .

[43]  Charbel Farhat,et al.  A hybrid discontinuous in space and time Galerkin method for wave propagation problems , 2014 .

[44]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[45]  A. Taflove,et al.  Numerical Solution of Steady-State Electromagnetic Scattering Problems Using the Time-Dependent Maxwell's Equations , 1975 .

[46]  Hong Zheng,et al.  Explicit Discontinuous Deformation Analysis Method with Lumped Mass Matrix for Highly Discrete Block System , 2018, International Journal of Geomechanics.

[47]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[48]  Ivo Babuška,et al.  The generalized finite element method for Helmholtz equation: Theory, computation, and open problems , 2006 .

[49]  M. S. Mohamed,et al.  Enriched finite elements for initial-value problem of transverse electromagnetic waves in time domain , 2017 .

[50]  R. J. Astley,et al.  Modelling of short wave diffraction problems using approximating systems of plane waves , 2002 .

[51]  Jin-Fa Lee,et al.  Finite Element Methods for Microwave Engineering , 2004 .

[52]  Jon Trevelyan,et al.  Wave interpolation finite elements for Helmholtz problems with jumps in the wave speed , 2005 .

[53]  Seyyed Mehdi Mirimani,et al.  Electromagnetic and thermal analysis of underground power solid‐conductor cables under harmonic and unbalancing currents based on FEM , 2018 .

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

[55]  Ivo Babuška,et al.  The generalized finite element method for Helmholtz equation. Part II: Effect of choice of handbook functions, error due to absorbing boundary conditions and its assessment , 2008 .

[56]  R. Alferness,et al.  Guided-wave optoelectronics , 1988 .