A Wavelet-Based Adaptive Finite Element Method for the Stokes Problems

In this work, we present the mathematical formulation of the new adaptive multiresolution method for the Stokes problems of highly viscous materials arising in computational geodynamics. The method is based on particle-in-cell approach-the Stokes system is solved on a static Eulerian finite element grid and material properties are carried in space by Lagrangian material points. The Eulerian grid is adapted using the wavelet-based adaptation algorithm. Both bilinear (Q(1)P(0), Q(1)Q(1)) and biquadratic (Q(2)P(-1)) mixed approximations for the Stokes system are supported. The proposed method is illustrated for a number of linear and nonlinear two-dimensional benchmark problems of geophysical relevance. The results of the adaptive numerical simulations using the proposed method are in an excellent agreement with those obtained on non-adaptive grids and with analytical solutions, while computational requirements are few orders of magnitude less compared to the non-adaptive simulations in terms of both time and memory usage.