TY - JOUR
AU - MOTAILO, А.P.
PY - 2020/10/06
Y2 - 2021/08/01
TI - CUBATURE FORMULA FOR AN OCTAHEDRON OF THE SEVENTH ALGEBRAIC ORDER OF ACCURACY
JF - APPLIED QUESTIONS OF MATHEMATICAL MODELLING
JA - ППММ
VL - 3
IS - 2.2
SE -
DO - 10.32782/KNTU2618-0340/2020.3.2-2.18
UR - //ojs.kntu.net.ua/index.php/aqmm/article/view/48
SP - 184-193
AB - <p><em>When solving the problems of mathematical physics by the finite element method for volume regions using lattices of a tetrahedral-octahedral structure, there is the problem of choosing a specific basis for the octahedron and the formula for numerical integration over this polyhedron. The numerical solution of the problem is the solution of a system of linear algebraic equations with coefficients that are elements of the stiffness and mass matrices. The accuracy of the solution of the boundary problem depends on the accuracy of the cubature formulas for the octahedron.</em></p><p><em>When</em> <em>the computational domain is discretized by the linear octahedron and tetrahedron</em><em>, </em><em>the problem of numerical integration over the octahedron region is partially solved</em><em>. </em><em>Cubature</em> <em>formulas are constructed for calculating the local stiffness matrix for an octahedron with piecewise linear</em><em>, </em><em>trigonometric and</em><em> second-order </em><em>polynomial bases</em><em>. </em><em>The cubature formula for calculating the elements of the local mass matrix is constructed for an octahedron with</em><em> a </em><em>trigonometric basis</em><em>. </em><em>Cubature formulas for an octahedron with trigonometric and</em><em> second-order </em><em>polynomial bases are exact for</em><em> a </em><em>trigonometric partial form and</em><em> third-</em><em>order algebraic polynomials</em><em>, </em><em>respectively</em><em>, </em><em>and contain</em><em> a </em><em>minimal number of interpolation nodes</em><em>.</em></p><p><em>In</em> <em>this</em> <em>paper</em><em>, a </em><em>cubature</em> <em>formula for</em><em> a </em><em>quadratic octahedron with</em><em> a fourth-order </em><em>polynomial basis is constructed</em><em>. </em><em>This formula is exact for</em><em> seventh-order </em><em>algebraic polynomials and has two different sets of node coordinates and weight coefficients</em><em>. </em><em>An estimate of the remainder term of the cubature formula for integrand functions of the class </em> <em> is obtained</em><em>. </em><em>Theoretical results were verified by calculating the elements of the local stiffness matrix for</em><em> a fourth-order </em><em>polynomial system of basis functions of</em><em> a </em><em>quadratic octahedron</em><em>. </em><em>Based on the calculation results</em><em>, </em><em>the cubature formula optimal in accuracy is determined</em><em>. </em><em>The weighting coefficients of the formula are positive</em><em>; </em><em>one of the four groups of interpolation nodes does not belong to the region of the octahedron</em><em>.</em></p><p><em>This cubature formula can be used to solve the boundary problems of mathematical physics for volume regions that are discretized by the lattice of the</em><em> tetrahedral-octahedral </em><em>structure</em><em>.</em></p>
ER -