Автор благодарит О. Чулуунбаатара, С.И. Виницкого, В.Л. Дербова, В.П. Гердта, А. Гужджа и Л.А. Севастьянова за сотрудничество.
A new computational scheme of the finite element method of a high order of accuracy for solving boundary value problems for an elliptic partial differential equation that preserves the continuity of the derivatives of the approximate solution in a bounded domain of a multidimensional Euclidean space is proposed. A piecewise continuous basis of the finite element method is generated using interpolation Hermite polynomials of several variables and ensures the continuity of not only the approximate solution but also its derivatives up to a given order on the boundaries of finite elements, depending on the smoothness of the variable coefficients of the equation and the boundary of the domain. The efficiency and accuracy order of the computational scheme, algorithm and program are demonstrated by the example of an exactly solvable boundary-value problem for a triangular membrane depending on the number of finite elements of the partition of the domain and the dimension of the eigenvector of the algebraic problem. It was shown that, in order to achieve a given accuracy of the approximate solution, for schemes of the finite element method with Hermite interpolation polynomials the dimension of the eigenvector is approximately two times smaller than for schemes with Lagrange interpolation polynomials that preserve on the boundaries of finite elements only the continuity of the approximate solution. The high-accuracy computational scheme of the finite element method is oriented to calculations of the spectral and optical characteristics of quantum-mechanical systems.
Предложена новая вычислительная схема метода конечных элементов высокого порядка точности решения краевых задач для эллиптического уравнения в частных производных, сохраняющая непрерывность производных приближенного решения в ограниченной области многомерного евклидова пространства. Кусочно-непрерывный базис метода конечных элементов генерируется с помощью интерполяционных полиномов Эрмита нескольких переменных и обеспечивает непрерывность не только приближенного решения, но и его производных до заданного порядка на границах конечных элементов в зависимости от гладкости переменных коэффициентов уравнения и границы области. Эффективность и порядок точности вычислительной схемы, алгоритма и программы демонстрируется на примере точно-решаемой краевой задачи на собственные значения для треугольной мембраны в зависимости от числа конечных элементов разбиения области и от размерности собственного вектора алгебраической задачи. Показано, что для достижения заданной точности приближённого решения схемой метода конечных элементов с интерполяционными полиномами Эрмита длина собственного вектора примерно в два раза меньше, чем для схем с интерполяционными полиномами Лагранжа, сохраняющих на границах конечных элементов только непрерывность приближённого решения. Вычислительная схема метода конечных элементов высокого порядка точности ориентирована на расчёты спектральных и оптических характеристик квантовомеханических систем.