Théophile Chaumont: Monday 18 September at 1 pm, A415 Inria Paris.
Time-harmonic wave propagation problems are costly to solve numerically
since the corresponding PDE operators are not strongly elliptic, and as a result,
discretization methods might become unstable. Specifically, the finite element solution is
quasi-optimal (almost as good as the best approximation the finite element space
can provide) only under restrictive assumptions on the mesh size. If the mesh
size is too large, stability is lost, and the finite element solution can become
completely inaccurate, even when the best approximation is. This phenomenon is called the
« pollution effect » and becomes more important for larger frequencies.
For the case of wave propagation problems in homogeneous media, it is known that
high order finite element methods are less sensitive to the pollution effect. For this
reason, they are employed in a wide range of applications, as the corresponding linear
systems are smaller and easier to solve.
In this talk, we investigate the use of high order finite element methods to solve
wave propagation problems in highly heterogeneous media. Since the heterogeneities
of the medium can exhibit small scale features, we consider « non-fitting » meshes,
that are not aligned with the physical interfaces of the medium. Instead, the parameters
defining the medium of propagation can be discontinuous inside each element. We propose
a convergence analysis and draw two main conclusions:
– the asymptotic convergence rate of the proposed finite element method is suboptimal
due to the lack of regularity of the solution inside each cell
– the pollution effect is greatly reduced by increasing the order of discretization.
We illustrate our main conclusions with geophysical application benchmarks.
These examples confirm that higher order methods are more efficient than
linear finite elements.