TY - JOUR

T1 - New 25-point stencils with optimal accuracy for 2-D heat transfer problems. Comparison with the quadratic isogeometric elements

AU - Idesman, A.

AU - Dey, B.

N1 - Publisher Copyright:
© 2020 Elsevier Inc.

PY - 2020/10/1

Y1 - 2020/10/1

N2 - A new approach for the increase in the order of accuracy of high order elements used for the time dependent heat equation and for the time independent Poisson equation has been suggested on uniform square and rectangular meshes. It is based on the optimization of the coefficients of the corresponding discrete stencil equation with respect to the local truncation error. By a simple modification of the coefficients of 25-point stencils, the new approach exceeds the accuracy of the quadratic isogeometric elements by four orders for the heat equation and by twelve orders for the Poisson equation. Despite the significant increase in accuracy, the computational costs of the new technique are the same as those for the conventional quadratic isogeometric elements on a given mesh. The numerical examples are in a good agreement with the theoretical results for the new approach and also show that the new approach is much more accurate than the conventional isogeometric elements at the same number of degrees of freedom. Hybrid methods that combine the new stencils with the conventional isogeometric and finite elements and can be applied to irregular domains are also presented.

AB - A new approach for the increase in the order of accuracy of high order elements used for the time dependent heat equation and for the time independent Poisson equation has been suggested on uniform square and rectangular meshes. It is based on the optimization of the coefficients of the corresponding discrete stencil equation with respect to the local truncation error. By a simple modification of the coefficients of 25-point stencils, the new approach exceeds the accuracy of the quadratic isogeometric elements by four orders for the heat equation and by twelve orders for the Poisson equation. Despite the significant increase in accuracy, the computational costs of the new technique are the same as those for the conventional quadratic isogeometric elements on a given mesh. The numerical examples are in a good agreement with the theoretical results for the new approach and also show that the new approach is much more accurate than the conventional isogeometric elements at the same number of degrees of freedom. Hybrid methods that combine the new stencils with the conventional isogeometric and finite elements and can be applied to irregular domains are also presented.

KW - Heat equation

KW - High-order isogeometric elements

KW - Increase in the order of accuracy

KW - Local truncation error

KW - Poisson equation

UR - http://www.scopus.com/inward/record.url?scp=85086727687&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2020.109640

DO - 10.1016/j.jcp.2020.109640

M3 - Article

AN - SCOPUS:85086727687

SN - 0021-9991

VL - 418

JO - Journal of Computational Physics

JF - Journal of Computational Physics

M1 - 109640

ER -