Multiscale and Stabilized Methods

Abstract

This article presents an introduction to multiscale and stabilized methods, which represent unified approaches to modeling and numerical solution of fluid dynamic phenomena. Finite element applications are emphasized but the ideas are general and apply to other numerical methods as well. (They have been used in the development of finite difference, finite volume, and spectral methods, in addition to finite‐element methods.) The analytical ideas are first illustrated for time‐harmonic wave‐propagation problems in unbounded fluid domains governed by the Helmholtz equation. This leads to the well‐known Dirichlet‐to‐Neumann formulation. A general treatment of the variational multiscale method in the context of an abstract Dirichlet problem is then presented, which is applicable to advective–diffusive processes and other processes of physical interest. It is shown how the exact theory represents a paradigm for subgrid‐scale models and a posteriori error estimation. Hierarchical p ‐methods and bubble function methods are examined in order to understand and ultimately approximate the ‘fine‐scale Green's function’ that appears in the theory. Relationships among so‐called residual‐free bubbles, element Green's functions, and stabilized methods are exhibited. These ideas are then generalized to a class of nonsymmetric, linear evolution operators formulated in space–time. The variational multiscale method also provides guidelines and inspiration for the development of stabilized methods (e.g. SUPG, GLS, etc.) that have attracted considerable interest and have been extensively utilized in engineering and the physical sciences. An overview of stabilized methods for advective–diffusive equations is presented. A variational multiscale treatment of incompressible viscous flows, including turbulence is also described. This represents an alternative formulation of large eddy simulation (LES), which provides a simplified theoretical framework of LES with potential for improved modeling.

Stabilized bi-cubic Hermite Bézier finite element method with application to gas-plasma interactions occurring during massive material injection in tokamaks

“…The evaluation of the exact stabilization matrix (L −1 ) is the main ingredient of the VMS stabilized FEM [24,25]. For linear operators, it is possible to express the stabilization matrix (L −1 ) without any approximation [24]. However, in general, evaluating the exact form of L −1 poses a more difficult problem than the original one.…”

“…Indeed, Galerkin FEMs do not provide mechanisms to control the sub-scales effects on the resolved scales: stabilization. The vast literature exist for stabilized Galerkin FEM such as Taylor-Galerkin methods [19,20], bubble functions [21,22], Streamline Upwind Petrov-Galerkin (SUPG) methods [23], Variational Multiscale (VMS) [24,25] methods etc. The VMS approach provides attractive guidelines for developing stabilized FEMs where numerical stabilization is achieved by an additional contribution to Galerkin's weak formulation.…”