1887

Abstract

Summary

Spectral elements with mass lumping allow for explicit time stepping and are therefore attractive for modelling seismic wave propagation. Their formulation on rectangular elements is straighforward, but for tetrahedra only elements up to degree 3 are known. To preserve accuracy after mass lumping, these elements require additional nodes that make them computationally more expensive. Here, we propose a new, less restrictive accuracy condition for the construction of these continuous mass-lumped elements. This enables us to construct several new tetrahedral elements. The new degree-2 and degree-3 elements require 15 and 32 nodes, while the existing ones have 23 and 50 nodes per element, respectively. We also developed degree-4 tetrahedral elements with 60, 61, or 65 nodes per element. Numerical examples confirm that the various mass-lumped elements maintain the optimal order of accuracy and show that the new elements are significantly more efficient in terms of accuracy versus compute time than the existing ones.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201800964
2018-06-11
2024-04-26
Loading full text...

Full text loading...

References

  1. Chin-Joe-Kong, M.J.S., Mulder, W.A. and van Veldhuizen, M.
    [1999] Higher-order triangular and tetra-hedral finite elements with mass lumping for solving the wave equation. Journal of Engineering Mathematics, 35, 405–426.
    [Google Scholar]
  2. Ciarlet, P.G.
    [1978] The finite element method for elliptic problems, Studies in mathematics and its applications, 4. North-Holland.
    [Google Scholar]
  3. Cohen, G., Joly, P., Roberts, J.E. and Tordjman, N.
    [2001] Higher order triangular finite elements with mass lumping for the wave equation. SIAM Journal on Numerical Analysis, 38(6), 2047–2078.
    [Google Scholar]
  4. Cohen, G., Joly, P. and Tordjman, N.
    [1995] Higher order triangular finite elements with mass lumping for the wave equation. In: Cohen, G., Bécache, E., P.J. and Roberts, J.E. (Eds.) Proceedings of the Third International Conference on Mathematical and Numerical Aspects of Wave Propagation. SIAM, Philadelphia, 270–279.
    [Google Scholar]
  5. Cui, T., Leng, W., Lin, D., Ma, S. and Zhang, L.
    [2017] High order mass-lumping finite elements on simplexes. Numerical Mathematics: Theory, Methods and Applications, 10(2), 331–350.
    [Google Scholar]
  6. Dablain, M.A
    . [1986] The application of high-order differencing to the scalar wave equation. Geophysics, 51(1), 54–66.
    [Google Scholar]
  7. Fried, I. and Malkus, D.S.
    [1975] Finite element mass matrix lumping by numerical integration with no convergence rate loss. International Journal of Solids and Structures, 11(4), 461–466.
    [Google Scholar]
  8. Kononov, A., Minisini, S., Zhebel, E. and Mulder, W.A.
    [2012] A 3D tetrahedral mesh generator for seismic problems. In: Proceedings of the 74th EAGE Conference & Exhibition. B006.
    [Google Scholar]
  9. Liu, Y., Teng, J., Xu, T. and Badal, J.
    [2017] Higher-order triangular spectral element method with optimized cubature points for seismic wavefield modeling. Journal of Computational Physics, 336, 458–480.
    [Google Scholar]
  10. Mulder, W.A.
    [1996] A comparison between higher-order finite elements and finite differences for solving the wave equation. In: Désidéri, J.A., LeTallec, P., Oñate, E., Périaux, J. and Stein, E. (Eds.) Proceedings of the Second ECCOMAS Conference on Numerical Methods in Engineering. John Wiley & Sons, Chichester, 344–350.
    [Google Scholar]
  11. [2013] New triangular mass-lumped finite elements of degree six for wave propagation. Progress In Electromagnetics Research, 141, 671–692.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201800964
Loading
/content/papers/10.3997/2214-4609.201800964
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error