On the Application of the Continuous Galerkin Finite Element Method for Conservation Problems. Bush, L. & Ginting, V. SIAM Journal on Scientific Computing, 35(6):A2953-A2975, 2013.
On the Application of the Continuous Galerkin Finite Element Method for Conservation Problems [link]Paper  doi  abstract   bibtex   
One major drawback that prevents the use of the standard continuous Galerkin finite element method in solving conservation problems is its lack of a locally conservative flux. Our present work has developed a simple postprocessing for the continuous Galerkin finite element method resulting in a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. The postprocessing requires an auxiliary fully Neumann problem to be solved on each finite element. These local problems are independent of each other and in two dimensions involve solving only a 3-by-3 system in the case of triangular elements and a 4-by-4 system for quadrilateral elements. A convergence analysis for the method is provided and its performance is demonstrated through numerical examples of multiphase flow with triangular and quadrilateral elements along with a description of its parallel implementation.
@article {MR3141756,
    AUTHOR = {Bush, L. and Ginting, V.},
     TITLE = {On the {A}pplication of the {C}ontinuous {G}alerkin {F}inite {E}lement
              {M}ethod for {C}onservation {P}roblems},
  JOURNAL = {SIAM Journal on Scientific Computing},
    VOLUME = {35},
      YEAR = {2013},
    NUMBER = {6},
     PAGES = {A2953-A2975},
      ISSN = {1064-8275},
   MRCLASS = {65N30 (65N15 65N50 65Y05)},
  MRNUMBER = {3141756},
MRREVIEWER = {Marius Ghergu},
       DOI = {10.1137/120900393},
       URL = {http://dx.doi.org/10.1137/120900393},
       ABSTRACT="
       One major drawback that prevents the use of the standard continuous Galerkin finite element method in solving conservation problems is its lack of a locally conservative flux. Our present work has developed a simple postprocessing for the continuous Galerkin finite element method resulting in a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. The postprocessing requires an auxiliary fully Neumann problem to be solved on each finite element. These local problems are independent of each other and in two dimensions involve solving only a 3-by-3 system in the case of triangular elements and a 4-by-4 system for quadrilateral elements. A convergence analysis for the method is provided and its performance is demonstrated through numerical examples of multiphase flow with triangular and quadrilateral elements along with a description of its parallel implementation."
}

Downloads: 0