Construction of Locally Conservative Fluxes for the SUPG Method. Deng, Q. & Ginting, V. Numerical Methods for Partial Differential Equations, 31(6):1971-1994, 2015.
Construction of Locally Conservative Fluxes for the SUPG Method [link]Paper  doi  abstract   bibtex   
We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov-Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift-diffusion equations.
@article {NUM:NUM21975,
author = {Deng, Q. and Ginting, V.},
title = {Construction of {L}ocally {C}onservative {F}luxes for the {SUPG} {M}ethod},
journal = {Numerical Methods for Partial Differential Equations},
volume = {31},
number = {6},
issn = {1098-2426},
url = {http://dx.doi.org/10.1002/num.21975},
doi = {10.1002/num.21975},
pages = {1971-1994},
keywords = {advection diffusion, advection dominated, CGFEM, conservative flux, postprocessing, SUPG},
year = {2015},
abstract="
We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov-Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift-diffusion equations."
}

Downloads: 0