Метод конечных элементов высокого порядка точности решения краевых задач для эллиптического уравнения в частных производных

Обложка

Цитировать

Полный текст

Аннотация

Предложена новая вычислительная схема метода конечных элементов высокого порядка точности решения краевых задач для эллиптического уравнения в частных производных, сохраняющая непрерывность производных приближенного решения в ограниченной области многомерного евклидова пространства. Кусочно-непрерывный базис метода конечных элементов генерируется с помощью интерполяционных полиномов Эрмита нескольких переменных и обеспечивает непрерывность не только приближенного решения, но и его производных до заданного порядка на границах конечных элементов в зависимости от гладкости переменных коэффициентов уравнения и границы области. Эффективность и порядок точности вычислительной схемы, алгоритма и программы демонстрируется на примере точно-решаемой краевой задачи на собственные значения для треугольной мембраны в зависимости от числа конечных элементов разбиения области и от размерности собственного вектора алгебраической задачи. Показано, что для достижения заданной точности приближённого решения схемой метода конечных элементов с интерполяционными полиномами Эрмита длина собственного вектора примерно в два раза меньше, чем для схем с интерполяционными полиномами Лагранжа, сохраняющих на границах конечных элементов только непрерывность приближённого решения. Вычислительная схема метода конечных элементов высокого порядка точности ориентирована на расчёты спектральных и оптических характеристик квантовомеханических систем.

Об авторах

Александр Александрович Гусев

Объединённый институт ядерных исследований

Автор, ответственный за переписку.
Email: gooseff@jinr.ru

Автор благодарит О. Чулуунбаатара, С.И. Виницкого, В.Л. Дербова, В.П. Гердта, А. Гужджа и Л.А. Севастьянова за сотрудничество.

ул. Жолио-Кюри, д. 6, г. Дубна, Московская область, Россия, 141980

Список литературы

  1. A.A. Gusev, O. Chuluunbaatar, S.I. Vinitsky, et al., Symbolic-Numerical Solution of Boundary-Value Problems with Self-Adjoint Second-Order Differential Equation using the Finite Element Method with Interpolation Hermite Polynomials, Lecture Notes in Computer Science 8660 (2014) 138–154.
  2. A.A. Gusev, L.L. Hai, O. Chuluunbaatar, S.I. Vinitsky, Program KANTBP 4M for Solving Boundary-Value Problems for Systems of Ordinary Differential Equations of the Second Order. URL http://wwwinfo.jinr.ru/programs/jinrlib.
  3. P. Ciarlet, The Finite Element Method for Elliptic Problems, North-holland Publ. Comp, Amsterdam, 1978.
  4. L. Ramdas Ram-Mohan, Finite Element and Boundary Element Applications in Quantum Mechanics, Oxford Univ. Press, New York, 2002.
  5. A.A. Gusev, O. Chuluunbaatar, S.I. Vinitsky, V.L. Derbov, A.G´o´zd´z, Algorithms for Solving the Parametric Self-Adjoint 2D Elliptic Boundary-Value Problem Using High-Accuracy Finite Element Method, RUDN Journal of Mathematics, Information Sciences and Physics 25 (1) (2017) 36–55.
  6. O.A. Ladyzhenskaya, The Boundary Value Problems of Mathematical Physics, Vol. 49, Springer, Berlin, 1985.
  7. V.V. Shaidurov, Multigrid Methods for Finite Elements, Springer, Berlin, 1995.
  8. K.J. Bathe, Finite Element Procedures in Engineering Analysis, Englewood Cliffs, Prentice Hall, New York, 1982.
  9. N.S. Bakhvalov, N.P. Zhidkov, G.M. Kobelkov, Numerical Methods, Nauka, Moscow, 1987, in Russian.
  10. E.B. Becker, G.F. Carey, J. Tinsley Oden, Finite elements. An introduction. Vol. I, Englewood Cliffs, Prentice Hall, New Jersey, 1981.
  11. G. Strang, G.J. Fix, An Analysis of the Finite Element Method, Prentice-Hall, Englewood Cliffs, New York, 1973.

Дополнительные файлы

Доп. файлы
Действие
1. JATS XML

Согласие на обработку персональных данных

 

Используя сайт https://journals.rcsi.science, я (далее – «Пользователь» или «Субъект персональных данных») даю согласие на обработку персональных данных на этом сайте (текст Согласия) и на обработку персональных данных с помощью сервиса «Яндекс.Метрика» (текст Согласия).