VMS spectral solution of two-dimensional advection-diffusion problems with anisotropic velocity
This work provides a more accurate stabilization method for finite element simulations of advection-diffusion problems with anisotropic velocity, which is relevant for computational fluid dynamics.
The authors extend the Variational Multi-scale method with spectral approximation to 2D advection-diffusion problems with anisotropic velocity, computing stabilized coefficients for triangular and quadrilateral grids. Numerical tests show a relevant accuracy gain for moderately large grid Péclet numbers with variable advection velocity.
In this article, we extend the Variational Multi-scale method with spectral approximation of the sub-scales to two-dimensional advection-diffusion problems. The spectral VMS method is cast for low-order elements as a standard VMS method with specific stabilized coefficients associated to a component of the advection velocity. We compute the stabilized coefficients for grids of isosceles right triangles and right quadrilaterals, based upon the explicit computation of the eigen-pairs of the advection-diffusion operator with Dirichlet boundary conditions. To reduce the computing time, the stabilized coefficients are computed at the nodes of a grid in an off-line step, and then interpolated by a fast procedure in the on-line computation. Finally, we present some numerical tests, first with constant velocity and after that with anisotropic variable velocity, in order to compare our results with those provided by other stabilization coefficients. We observe a relevant accuracy gain for moderately large grid Péclet numbers for variable advection velocity.