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 -