Research output: Contribution to journal › Article › peer-review
hp-Вариант метода коллокации и наименьших квадратов с интегральными коллокациями решения бигармонического уравнения. / Shapeev, Vasiliy Pavlovich; Брындин, Лука Сергеевич; Беляев, Василий Алексеевич.
In: Vestnik Samarskogo Gosudarstvennogo Tekhnicheskogo Universiteta, Seriya Fiziko-Matematicheskie Nauki, Vol. 26, No. 3, 8, 2022, p. 556-572.Research output: Contribution to journal › Article › peer-review
}
TY - JOUR
T1 - hp-Вариант метода коллокации и наименьших квадратов с интегральными коллокациями решения бигармонического уравнения
AU - Shapeev, Vasiliy Pavlovich
AU - Брындин, Лука Сергеевич
AU - Беляев, Василий Алексеевич
N1 - Шапеев В.П., Брындин Л.С., Беляев В.А. hp-Вариант метода коллокации и наименьших квадратов с интегральными коллокациями решения бигармонического уравнения // Вестник Самарского государственного технического университета. Серия: Физико-математические науки. – 2022. – Т. 26. - № 3. – С. 556-572. Институт теоретической и прикладной механики им. С. А. Христиановича СО РАН. Работа выполнена в рамках государственного задания (№ госрегистрации 121030500137-5).
PY - 2022
Y1 - 2022
N2 - Разработан новый алгоритм численного решения бигармонического уравнения. Он основан на впервые реализованном hp-варианте метода коллокации и наименьших квадратов (hp-МКНК) с интегральными коллокациями для эллиптического уравнения четвертого порядка в комбинации с современными способами ускорения итерационных процессов решения систем линейных алгебраических уравнений (СЛАУ). В hp-МКНК использовались его возможности измельчать шаги расчетной сетки (h-подход) и увеличивать степень базисных аппроксимирующих полиномов до произвольного порядка (p-подход). На примере численного моделирования изгиба шарнирно закрепленной изотропной пластины проведен анализ сходимости приближенных решений, полученных реализованным вариантом метода. Показано достижение высокой точности и повышенного порядка сходимости решений при применении полиномов высоких вплоть до десятой степеней в hp-МКНК. Исследована эффективность комбинированного применения сочетающихся с МКНК алгоритмов ускорения итерационных процессов решения СЛАУ. Применены предобуславливание матриц СЛАУ; алгоритм ускорения итераций, основанный на подпространствах Крылова; операция продолжения на многосеточном комплексе; распараллеливание вычислительной программы с помощью OpenMP; модифицированный алгоритм решения локальных СЛАУ, определяющих решение задачи в каждой ячейке сетки. Последний, применимый в случае решения линейного дифференциального уравнения, позволяет более эффективно решать переопределенные СЛАУ в МКНК, реализуемом итерациями по подобластям, в которых вид матриц локальных СЛАУ не изменяется на каждой итерации. Комбинированное применение всех перечисленных способов ускорения уменьшило время расчетов на персональном компьютере более, чем в 350 раз по сравнению со случаем, когда использовалось только предобуславливание.
AB - Разработан новый алгоритм численного решения бигармонического уравнения. Он основан на впервые реализованном hp-варианте метода коллокации и наименьших квадратов (hp-МКНК) с интегральными коллокациями для эллиптического уравнения четвертого порядка в комбинации с современными способами ускорения итерационных процессов решения систем линейных алгебраических уравнений (СЛАУ). В hp-МКНК использовались его возможности измельчать шаги расчетной сетки (h-подход) и увеличивать степень базисных аппроксимирующих полиномов до произвольного порядка (p-подход). На примере численного моделирования изгиба шарнирно закрепленной изотропной пластины проведен анализ сходимости приближенных решений, полученных реализованным вариантом метода. Показано достижение высокой точности и повышенного порядка сходимости решений при применении полиномов высоких вплоть до десятой степеней в hp-МКНК. Исследована эффективность комбинированного применения сочетающихся с МКНК алгоритмов ускорения итерационных процессов решения СЛАУ. Применены предобуславливание матриц СЛАУ; алгоритм ускорения итераций, основанный на подпространствах Крылова; операция продолжения на многосеточном комплексе; распараллеливание вычислительной программы с помощью OpenMP; модифицированный алгоритм решения локальных СЛАУ, определяющих решение задачи в каждой ячейке сетки. Последний, применимый в случае решения линейного дифференциального уравнения, позволяет более эффективно решать переопределенные СЛАУ в МКНК, реализуемом итерациями по подобластям, в которых вид матриц локальных СЛАУ не изменяется на каждой итерации. Комбинированное применение всех перечисленных способов ускорения уменьшило время расчетов на персональном компьютере более, чем в 350 раз по сравнению со случаем, когда использовалось только предобуславливание.
UR - https://www.scopus.com/inward/record.url?eid=2-s2.0-85147779776&partnerID=40&md5=ef1b1c93685f1424c64806ba2b82d4ef
UR - https://elibrary.ru/item.asp?id=49802895
UR - https://www.mendeley.com/catalogue/9ce70ec5-2230-339c-a3d4-b4ce4370928a/
U2 - 10.14498/vsgtu1936
DO - 10.14498/vsgtu1936
M3 - статья
VL - 26
SP - 556
EP - 572
JO - Vestnik Samarskogo Gosudarstvennogo Tekhnicheskogo Universiteta, Seriya Fiziko-Matematicheskie Nauki
JF - Vestnik Samarskogo Gosudarstvennogo Tekhnicheskogo Universiteta, Seriya Fiziko-Matematicheskie Nauki
SN - 1991-8615
IS - 3
M1 - 8
ER -
ID: 45655876