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

Mayank Drolia, M. Shadi Mohamed, Omar Laghrouche, Mohammed Seaid, Abdellah El-Kacimi

Research output: Contribution to journalArticle

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.
Original languageEnglish
Pages (from-to)1273-1293
Number of pages21
JournalApplied Mathematical Modelling
Volume77
Early online date12 Aug 2019
DOIs
Publication statusPublished - 1 Jan 2020

Fingerprint

Program processors
Wave propagation
Polynomials
Finite element method
Processing

Keywords

  • Field enrichment
  • Finite elements
  • Lumped mass
  • Partition of unity
  • Time domain
  • Wave equation

ASJC Scopus subject areas

  • Modelling and Simulation
  • Applied Mathematics

Cite this

@article{a594d73ce13c4fd6971f1b8a5cfa1e26,
title = "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.",
keywords = "Field enrichment, Finite elements, Lumped mass, Partition of unity, Time domain, Wave equation",
author = "Mayank Drolia and Mohamed, {M. Shadi} and Omar Laghrouche and Mohammed Seaid and Abdellah El-Kacimi",
year = "2020",
month = "1",
day = "1",
doi = "10.1016/j.apm.2019.07.054",
language = "English",
volume = "77",
pages = "1273--1293",
journal = "Applied Mathematical Modelling",
issn = "0307-904X",
publisher = "Elsevier",

}

Explicit time integration with lumped mass matrix for enriched finite elements solution of time domain wave problems. / Drolia, Mayank; Mohamed, M. Shadi; Laghrouche, Omar; Seaid, Mohammed; El-Kacimi, Abdellah.

In: Applied Mathematical Modelling, Vol. 77, 01.01.2020, p. 1273-1293.

Research output: Contribution to journalArticle

TY - JOUR

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

AU - Drolia, Mayank

AU - Mohamed, M. Shadi

AU - Laghrouche, Omar

AU - Seaid, Mohammed

AU - El-Kacimi, Abdellah

PY - 2020/1/1

Y1 - 2020/1/1

N2 - 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.

AB - 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.

KW - Field enrichment

KW - Finite elements

KW - Lumped mass

KW - Partition of unity

KW - Time domain

KW - Wave equation

UR - http://www.scopus.com/inward/record.url?scp=85074439956&partnerID=8YFLogxK

U2 - 10.1016/j.apm.2019.07.054

DO - 10.1016/j.apm.2019.07.054

M3 - Article

VL - 77

SP - 1273

EP - 1293

JO - Applied Mathematical Modelling

JF - Applied Mathematical Modelling

SN - 0307-904X

ER -