The Discontinuous Enrichment Method (DEM) for Multi-scale Transport Problems

Download or read book The Discontinuous Enrichment Method (DEM) for Multi-scale Transport Problems written by Irina Kalashnikova and published by Stanford University. This book was released on 2011 with total page 178 pages. Available in PDF, EPUB and Kindle. Book excerpt: A discontinuous enrichment method (DEM) for the efficient finite element solution of advection-dominated transport problems in fluid mechanics whose solutions are known to possess multi-scale features is developed. Attention is focused specifically on the two-dimensional (2D) advection-diffusion equation, the usual scalar model for the Navier-Stokes equations. Following the basic DEM methodology [1], the usual Galerkin polynomial approximation is locally enriched by the free-space solutions to the governing homogeneous partial differential equation (PDE). For the constant-coefficient advection-diffusion equation, several families of free-space solutions are derived. These include a family of exponential functions that exhibit a steep gradient in some flow direction, and a family of discontinuous polynomials. A parametrization of the former class of functions with respect to an angle parameter is developed, so as to enable the systematic design and implementation of DEM elements of arbitrary orders. It is shown that the original constant-coefficient methodology has a natural extension to variable-coefficient advection-diffusion problems. For variable-coefficient transport problems, the approximation properties of DEM can be improved by augmenting locally the enrichment space with a "higher-order" enrichment function that solves the governing PDE with the advection field a(x) linearized to second order. A space of Lagrange multipliers, introduced at the element interfaces to enforce a weak continuity of the solution and related to the normal derivatives of the enrichment functions, is developed. The construction of several low and higher-order DEM elements fitting this paradigm is discussed in detail. Numerical results for several constant as well as variable-coefficient advection-diffusion benchmark problems reveal that these DEM elements outperform their standard Galerkin and stabilized Galerkin counterparts of comparable computational complexity by a large margin, especially when the flow is advection-dominated.