Indexed on: 19 Jan '11Published on: 19 Jan '11Published in: Journal of Scientific Computing
We present a model problem for benchmarking codes that investigate magma migration in the Earth’s interior. This system retains the essential features of more sophisticated models, yet has the advantage of possessing solitary wave solutions. The existence of such exact solutions to the nonlinear problem make it an excellent benchmark problem for combinations of solver algorithms. In this work, we explore a novel algorithm for computing high quality approximations of the solitary waves in 1-, 2- and 3 dimensions and use them to benchmark a semi-Lagrangian Crank-Nicolson scheme for a finite element discretization of the time dependent problem.