1887

Abstract

Summary

The representation of a force or of a moment point source in a spectral finite-element code for modelling elastic wave propagation becomes fundamentally different in degenerate cases where the source is located on the boundary of an element. This difference is related to the fact that the finite-element basis functions are continuous across element boundaries, but their derivatives are not. A method is presented that effectively deals with this problem. Tests on 1-D elements show that the numerical errors for a force source follow the expected convergence rate in terms of the element size, apart from isolated cases where superconvergence occurs. For a moment source, the method also converges but one order of accuracy is lost, probably because of the reduced regularity of the problem. Numerical tests in 3D on continuous mass-lumped tetrahedral elements show a similar error behaviour as in the 1-D case, although in 3D, the loss of accuracy for the moment source is not a severe as a full order.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.202010044
2020-12-08
2024-04-20
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 tetrahedral finite elements with mass lumping for solving the wave equation.Journal of Engineering Mathematics, 35, 405–426.
    [Google Scholar]
  2. Fichtner, A.
    [2011] Full Seismic Waveform Modelling and Inversion. Springer-Verlag, Berlin, Heidelberg.
    [Google Scholar]
  3. Geevers, S., Mulder, W.A. and van der Vegt, J.J.W.
    [2018] New higher-order mass-lumped tetrahedral elements for wave propagation modelling.SIAM Journal on Scientific Computing, 40(5), A2830–A2857.
    [Google Scholar]
  4. [2019] Efficient quadrature rules for computing the stiffness matrices of mass-lumped tetrahedral elements for linear wave problems.SIAM Journal on Scientific Computing, 41(2), A1041–A1065.
    [Google Scholar]
  5. Komatitsch, D. and Vilotte, J.P.
    [1998] The spectral-element method: an efficient tool to simulate the seismic response of 2-D and 3-D geological structures.Bulletin of the Seismological Society of America, 88(2), 368–392.
    [Google Scholar]
  6. Lesage, A.C., Aubry, R., Houzeaux, G., Polo, M.A. and Cela, J.
    [2010] 3D spectral element method combined with h-refinement.72nd EAGE Conference & Exhibition, Barcelona, Spain, Extended Abstracts, C047.
    [Google Scholar]
  7. 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]
  8. Mulder, W.A. and Shamasundar, R.
    [2016] Performance of continuous mass-lumped tetrahedral elements for elastic wave propagation with and without global assembly.Geophysical Journal International, 207(1), 414–421.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.202010044
Loading
/content/papers/10.3997/2214-4609.202010044
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