Numerical analysis of a method for solving 2D linear isotropic elastodynamics with traction free boundary condition using potentials and finite elements

Albella Martínez, Jorge; Imperiale, Sébastien; Joly, Patrick; Rodríguez, Jerónimo

doi:10.1090/mcom/3613

Numerical analysis of a method for solving 2D linear isotropic elastodynamics with traction free boundary condition using potentials and finite elements
HTML articles powered by AMS MathViewer

by Jorge Albella Martínez, Sébastien Imperiale, Patrick Joly and Jerónimo Rodríguez HTML | PDF

Math. Comp. 90 (2021), 1589-1636 Request permission

Abstract:

When solving 2D linear elastodynamic equations in homogeneous isotropic media, a Helmholtz decomposition of the displacement field decouples the equations into two scalar wave equations that only interact at the boundary. It is then natural to look for numerical schemes that independently solve the scalar equations and couple the solutions at the boundary. The case of rigid boundary condition was treated by Burel [Ph.D. thesis, Université Paris Sud-Paris XI (2014)] and Burel et al. [Numer. Anal. Appl. 5 (2012), pp. 136–143]. However the case of traction free boundary condition was proven by Martinez et al. [J. Sci. Comput. 77 (2018), pp. 1832–1873] to be unstable if a straightforward approach is used. Then an adequate functional framework as well as a time domain mixed formulation to circumvent these issues was presented. In this work we first review the formulation presented by Martinez et al. [J. Sci. Comput. 77 (2018), pp. 1832–1873] and propose a subsequent discretised formulation. We provide the complete stability analysis of the corresponding numerical scheme. Numerical results that illustrate the theory are also shown.

References

Adi Ben-Israel and Thomas N. E. Greville, Generalized inverses, 2nd ed., CMS Books in Mathematics/Ouvrages de Mathématiques de la SMC, vol. 15, Springer-Verlag, New York, 2003. Theory and applications. MR 1987382
Burel A., Contributions à la simulation numérique en élastodynamique: découplage des ondes P et S, modèles asymptotiques pour la traversée de couches minces, Ph.D. thesis, Université Paris Sud-Paris XI, 2014.
A. Burel, S. Imperiale, and P. Joly, Solving the homogeneous isotropic linear elastodynamics equations using potentials and finite elements. The case of the rigid boundary condition, Numer. Anal. Appl. 5 (2012), no. 2, 136–143.
P. Lailly, A. Bamberger, and G. Chavent, Etude de schémas numériques de l’élastodynamique linéaire, Rapport de recherche Inria 41 (1980).
Jorge Albella Martínez, Sébastien Imperiale, Patrick Joly, and Jerónimo Rodríguez, Solving 2D linear isotropic elastodynamics by means of scalar potentials: a new challenge for finite elements, J. Sci. Comput. 77 (2018), no. 3, 1832–1873. MR 3874796, DOI 10.1007/s10915-018-0768-9
Susanne C. Brenner and L. Ridgway Scott, The mathematical theory of finite element methods, 3rd ed., Texts in Applied Mathematics, vol. 15, Springer, New York, 2008. MR 2373954, DOI 10.1007/978-0-387-75934-0
Franco Brezzi and Michel Fortin, Mixed and hybrid finite element methods, Springer Series in Computational Mathematics, vol. 15, Springer-Verlag, New York, 1991. MR 1115205, DOI 10.1007/978-1-4612-3172-1
C. Amrouche, C. Bernardi, M. Dauge, and V. Girault, Vector potentials in three-dimensional non-smooth domains, Math. Methods Appl. Sci. 21 (1998), no. 9, 823–864 (English, with English and French summaries). MR 1626990, DOI 10.1002/(SICI)1099-1476(199806)21:9<823::AID-MMA976>3.0.CO;2-B
Robert Dautray and Jacques-Louis Lions, Mathematical analysis and numerical methods for science and technology. Vol. 5, Springer-Verlag, Berlin, 1992. Evolution problems. I; With the collaboration of Michel Artola, Michel Cessenat and Hélène Lanchon; Translated from the French by Alan Craig. MR 1156075, DOI 10.1007/978-3-642-58090-1
Philippe G. Ciarlet, Introduction à l’analyse numérique matricielle et à l’optimisation, Masson, Paris, 1982 (French). Collection Mathématiques Appliquées pour la Maîtrise. [Collection of Applied Mathematics for the Master’s Degree]. MR 680778
P. G. Ciarlet, The finite element method for elliptic problems, SIAM, Philadelphia, 2002.
G. Cohen, P. Joly, J. E. Roberts, and N. Tordjman, Higher order triangular finite elements with mass lumping for the wave equation, SIAM J. Numer. Anal. 38 (2001), no. 6, 2047–2078. MR 1856242, DOI 10.1137/S0036142997329554
I. Babuška, J. Osborn, and J. Pitkäranta, Analysis of mixed methods using mesh dependent norms, Math. Comp. 35 (1980), no. 152, 1039–1062. MR 583486, DOI 10.1090/S0025-5718-1980-0583486-7
R. Glowinski and O. Pironneau, Numerical methods for the first biharmonic equation and the two-dimensional Stokes problem, SIAM Rev. 21 (1979), no. 2, 167–212. MR 524511, DOI 10.1137/1021028
Vivette Girault and Pierre-Arnaud Raviart, Finite element methods for Navier-Stokes equations, Springer Series in Computational Mathematics, vol. 5, Springer-Verlag, Berlin, 1986. Theory and algorithms. MR 851383, DOI 10.1007/978-3-642-61623-5
J. Virieux, P-SV wave propagation in heterogeneous media: Velocity-stress finite-diference method, Geophysics 51 (1986), no. 4, 889–901.
Barbara I. Wohlmuth, Discretization methods and iterative solvers based on domain decomposition, Lecture Notes in Computational Science and Engineering, vol. 17, Springer-Verlag, Berlin, 2001. MR 1820470, DOI 10.1007/978-3-642-56767-4