High-performance finite-element simulations of seismic wave propagation in three-dimensional non linear inelastic geological media

Fabrice Dupros (1) , Florent de Martin (1) , Evelyne Foerster (1) , Dimitri Komatitsch (2, 3) , Jean Roman (4)

1 BRGM - Bureau de Recherches Géologiques et Minières (BRGM)
2 Magique 3D - Advanced 3D Numerical Modeling in Geophysics
3 IUF - Institut Universitaire de France
4 HiePACS - High-End Parallel Algorithms for Challenging Numerical Simulations

We present finite-element numerical simulations of seismic wave propagation in non linear inelastic geologicalmedia. We demonstrate the feasibility of large scale modeling based on an implicit numerical scheme and a nonlinear constitutive model. We illustrate our methodology with an application to regional scale modeling in the French Riviera, which is prone to earthquakes. The PaStiX direct solver is used to handle large matrix numerical factorizations based on hybrid parallelism to reduce memory overhead. A specific methodology is introduced for the parallel assembly in the context of soil nonlinearity. We analyse the scaling of the parallel algorithms on large-scale configurations and we discuss the physical results.

Keywords: seismic numerical simulation; finite-element method; parallel sparse direct solver; nonlinear soil behaviour

Domains: Earth Sciences

Fabrice Dupros, Florent de Martin, Evelyne Foerster, Dimitri Komatitsch, Jean Roman. High-performance finite-element simulations of seismic wave propagation in three-dimensional non linear inelastic geological media. Parallel Computing, 2010, 36, pp.308-325. ⟨10.1016/j.parco.2009.12.011⟩. ⟨hal-00537930⟩