IDEAS home Printed from https://ideas.repec.org/a/eee/matcom/v67y2004i1p55-66.html
   My bibliography  Save this article

Analysis and optimization of inner products for mimetic finite difference methods on a triangular grid

Author

Listed:
  • Liska, Richard
  • Shashkov, Mikhail
  • Ganzha, Victor

Abstract

The support operator method designs mimetic finite difference schemes by first constructing a discrete divergence operator based on the divergence theorem, and then defining the discrete gradient operator as the adjoint operator of the divergence based on the Gauss theorem connecting the divergence and gradient operators, which remains valid also in the discrete case. When evaluating the discrete gradient operator, one needs to define discrete inner products of two discrete vector fields. The local discrete inner product on a given triangle is defined by a 3×3 symmetric positive definite matrix M defined by its six independent elements–parameters. Using the Gauss theorem over our triangle, we evaluate the discrete gradient in the triangle. We require the discrete gradient to be exact for linear functions, which gives us a system of linear equations for elements of the matrix M. This system, together with inequalities which guarantee positive definiteness of the matrix M, results in a one parameter family of inner products which give exact gradients for linear functions. The traditional inner product is a member of this family. The positive free parameter can be used to improve another property of the discrete method. We show that accuracy of the method for quadratic functions improves with decreasing this parameter, however, at the same time, the condition number of the matrix M, which is the local matrix of the linear system for computing the discrete gradient, increases to infinity when the parameter goes to zero, so one needs to choose a compromise between accuracy and solvability of the local system. Our analysis has been performed by computer algebra tools which proved to be essential.

Suggested Citation

  • Liska, Richard & Shashkov, Mikhail & Ganzha, Victor, 2004. "Analysis and optimization of inner products for mimetic finite difference methods on a triangular grid," Mathematics and Computers in Simulation (MATCOM), Elsevier, vol. 67(1), pages 55-66.
  • Handle: RePEc:eee:matcom:v:67:y:2004:i:1:p:55-66
    DOI: 10.1016/j.matcom.2004.05.008
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0378475404001545
    Download Restriction: Full text for ScienceDirect subscribers only

    File URL: https://libkey.io/10.1016/j.matcom.2004.05.008?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    As the access to this document is restricted, you may want to search for a different version of it.

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:matcom:v:67:y:2004:i:1:p:55-66. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    We have no bibliographic references for this item. You can help adding them by using this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Catherine Liu (email available below). General contact details of provider: http://www.journals.elsevier.com/mathematics-and-computers-in-simulation/ .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.