A Projection Method for the Conservative Discretizations of Parabolic Partial Differential Equations

Darae Jeong, Junseok Kim

Research output: Contribution to journalArticle

Abstract

We present a projection method for the conservative discretizations of parabolic partial differential equations. When we solve a system of discrete equations arising from the finite difference discretization of the PDE, we can use iterative algorithms such as conjugate gradient, generalized minimum residual, and multigrid methods. An iterative method is a numerical approach that generates a sequence of improved approximate solutions for a system of equations. We repeat the iterative algorithm until a numerical solution is within a specified tolerance. Therefore, even though the discretization is conservative, the actual numerical solution obtained from an iterative method is not conservative. We propose a simple projection method which projects the non-conservative numerical solution into a conservative one by using the original scheme. Numerical experiments demonstrate the proposed scheme does not degrade the accuracy of the original numerical scheme and it preserves the conservative quantity within rounding errors.

Original languageEnglish
Pages (from-to)1-18
Number of pages18
JournalJournal of Scientific Computing
DOIs
Publication statusAccepted/In press - 2017 Aug 17

    Fingerprint

Keywords

  • Conservative discretization
  • Iterative methods
  • Projection method

ASJC Scopus subject areas

  • Theoretical Computer Science
  • Software
  • Engineering(all)
  • Computational Theory and Mathematics

Cite this