CC BY
65
УДК 519.63
Сеттиев Ш.Р., Ражабов Ж.Ш. Некоторые применения спектральных методов Галеркина
Сеттиев Шамсуддин Ражабович,
и.о. доцент кафедры, кандидат физико-математических наук. Филиал ФГБОУ ВПО «Российский экономический университет им.
Г.В.Плеханова» в г. Ташкенте. 100003. г. Ташкент, Узбекистанский проспект 49, корпус.
E- mail: [email protected] Ражабов Жалолиддин Шамсуддин угли, студент ФГАОУ ВПО «Северо-Кавказский федеральный университет».
355009, Ставропольский край, г. Ставрополь, ул. Пушкина, 1.
E- mail: [email protected] Аннотация. В данной работе рассматриваются способы применение спектральных методов к решениям классических задач механики. Излагается алгоритм решение уравнений Бюргерса с применением метода Галеркина. Результаты численных расчетов сопоставляются с точными решениями уравнений Бюргерса при различных значениях чисел Рейнольдса Re.
Ключевые слова: спектральный метод, метод Галеркина, уравнение Навье-Стокса, многочлены Чебышева, численное решение, быстрое преобразование Фурье.
Settiev Sh.R., Rajabov J.Sh. Some applications of Galerkin's spectral methods
Settiev Sh.R.,
acting associate Professor, candidate of physical and mathematical Sciences. Branch FGBOU VPO "Russian economic University named after. G. V. Plekhanov"
in Tashkent.
100003. Tashkent, Uzbekistan Avenue, 49, building.
E - mail: [email protected] Rajabov J.Sh.
student fsaei HPE "North-Caucasian Federal University". 355009, Stavropol territory, Stavropol, Pushkin street, 1.
E - mail: Jaloliddin.raj abov@mail .ru
Abstract. In this work there considered methods of application of spectral methods to the solution of classical tasks of mechanics. There recited algorithm of solution of Byurgers' equations with application of Galerkin method. The results of numerical calculations are compared with exact solution of Byurgers' equations at various values of numbers of Reynolds Re.
Keywords: spectral method, Galerkin method, equation of Navye-Stoks, multinomial of Chebyshev, numerical solution, fast transformation of Furye.
В данной работе приведем некоторые алгоритмы решения и численные реализации применения спектральных методов Галеркина к некоторым задачам, которые описываются уравнениями типа Навье-Стокса.
Уравнения Бюргерса. Рассмотрим одномерный вариант уравнения Бюргерса
ди , ди 1 д2и _
— + и------ = 0 (1)
dt дх Reдх2 v '
Это уравнение представляет собой одну из моделей уравнений Навье-Стокса и показывает баланс между нелинейным конвективным членом и диссипативным процессом. Будем рассматривать задачи о распространение волны. Волна представляет собой разрыв, расположенный в плоскости. В дальнейшем эта волна распространяется вправо, однако, будучи сначала разрывной, сглаживается за счет воздействия диссипативного вязкого процесса.
Будем искать решение уравнения Бюргерса при в области —1 <
х < 1.
Начальные условия:
и0 = и(х, 0) = 1 при — 1 < х < 0 ,
и0 = и(х, 0) = 0 при 0 < х < 1. (2)
Граничные условия:
u(—1,t) = 1 , и(1, t) = 0. (3)
Для численного решения задачи (1)-(3) применим спектральный метод Галеркина. Пробное решение представим в таком виде:
ua(x,t)=ZIJ=0aj(t)Tj(x) (4)
где Tj(x) - многочлены Чебышева, a.j- искомые коэффициенты разложения. Подставив (4) в (1), получим невязку R
R = Ij äj Tj + Zj aj & щ Z aJ 4? (5)
Система уравнений относительно коэффициентов разложения a.j получается из условия (внутреннее произведение равно нулю)
(R,Tk) = 0, k = 0,l,...,N -2 (6)
т.е.
f_±R• Tk(x)dx = 0.
Два коэффициента, входящие в пробное решение, выражаются через другие алгебраические коэффициенты за счет требования об удовлетворении граничных условий. Выполнение условия (6) дает:
MÀ + (B + C)A = 0,
где
А = Œj , M = mkj = (Tj ,Тк) ,
в=ьы=Ит1 di^ ,
C = ckj=1 ( Ц1 . (7)
KJ Re V дх ' dx ) v '
К уравнению (6) применялась теорема Грина. Начальные значения a.j определяются путем применения метода Галеркина к начальным данным
(иа-щ,Тк) = 0 , к = 0,1, ...,N — 2 .
Получим
MA = D , D = (щ, Tk) = f\ Tkdx (8)
Мы сформулировали общий подход применения спектрального метода Галеркина. Надо особо отметить, что если пробные функции составлены из
ортогональных базисных функций, решение поставленной задачи значительно
упрощается. Известно, что многочлены Чебышева ортогональные в отрезке
1 1 [-1 , 1], с весом ^ 2, мы берем и 2 вместо , где
VI-X2
Г0 — 1 , 7\ — х , 7^+1 — 2х7\ — 7^-
к — '/с-1
Далее опишем алгоритм решения уравнения Бюргерса, применяя псевдоспектральный подход метода Галеркина. Предположим, что на п -м слое времени нам известно решение и™. 1) Вычислим коэффициенты а,™, применяя быстрое преобразование Фурье к
<(**) — ^ , / — 1,2, ... ,2Я.
3) Вычислим
дх
2) Вычислим Ь™ по а]2 из рекуррентных соотношений (для —).
— ^¿^(хО , г — 1,2, .,2^.
4) Вычислим
5) Вычислим
6) Вычислим
^(хг) — <(хг)
Зх
^^ — /_11 ^ ^^
р!?-
I — 1,2,...,2М.
дх2 / ;
из рекуррентных соотношений (для ). 7) Вычислим
,п+1
^г1 - _„п , — ^ + .
8) Вычислим
а^1 — + /
"а?
Производные по времени аппроксимируем по следующим разностным схемам (9) и (10) [2], чтобы можно было сравнивать между собой:
[/п+1 — + ДtF"+I ( )
и
^П — ^п +
(10)
^п+1 — ^п + (1
1
При / — - схема устойчива. В нашем случае известно решение в нулевом слое
времени, т.е. выполнение начального условия. Далее выполняются требования
всех пунктов. Для последнего пункта используется двухшаговая схема (9). Во
избежание побочных погрешностей для вычисления и™(хг) использовались
точек. Интегрирование проводилось по формуле Симпсона высокого порядка
точности [1]. В пробном решении мы брали 10 членов разложения (4). Запишем
рекуррентные соотношения, которые используются в пунктах 2 и 6:
^ — ^+2 + , 2<£<М
— 0.5Ьз + «2 , Ь^+1 — Ь^+2 — 0, ( )
^ — ^+2 + 2^+1 , 2<£<М-1
51 —0.5ьз + ь2 , Ь^ —Ь^+1 —0 . (12)
Приведем полученные результаты в табличной форме, чтобы легко было их
сопоставлять с точными решениями [3].
Сравнение приближенного и точного решения при t — 0,23 , Де — 1,0
X Приближенное решение Точное решение
-1.00000 1.00000 0.9593
-0.89999 0.97498 0.9447
-0.79999 0.94698 0.9253
-0.69999 0.91565 0.9008
-0.59999 0.88062 0.8709
-0.49999 0.84155 0.8350
-0.39999 0.79824 0.7931
-0.29999 0.75063 0.7452
-0.19999 0.69887 0.6919
-0.99999 0.64337 0.6340
0.000000 0.58476 0.5728
0.100000 0.52388 0.5097
0.200000 0.46169 0.4464
0.300000 0.39918 0.3844
0.400000 0.33729 0.3254
0.500000 0.27675 0.2706
0.600000 0.21806 0.2210
0.700000 0.16137 0.1772
0.80000 0.10650 0.1395
0.90000 0.52955 0.1077
1.00000 0.00000 0.0817
Не очень хорошая сходимость в окрестностях граничных точек объясняется тем, что при численном расчете на каждом шаге по времени контролировалось выполнение граничных условий, и трудностями в нулевом слое времени при одновременном удовлетворении начальных и граничных условий. Численные расчеты проводились при различных значениях числа Рейнолдьса Я е = 1.0,10.0, 100.0. В восьмом пункте применялись разностные схемы (9) и (10). Надо отметить, что схема (9) дает хорошее приближение, но требует в два раза больше машинного времени, чем схема (10). Остальные результаты приведем в виде графических зависимостей и от х при различных значениях Яе , Ь.
Рис. 1
u
Рис. 2
u
1.20
0.80 —
0.40
0.00 —
. /......\ Re= 100.0
' "" ' \ t=0.92
-0.40
i i i i i i i i i | i i i i 1.00 -0.50
0.00
0.50
.....I X
1.00
Рис. 3
точное решение уравнения Бюргерса приближенное решение по схеме (9)
• • • *.....* • • • приближенное решение по схеме (10)
Список литературы
1. Демидович Б.П., Марон И.А. Основы вычислительной математики. -Москва.: Изд-во Лань, 2006. 672с.
2. Роуч П. Вычислительная гидродинамика. - М.: Мир, 1980. -612c.
3. Флетчер К. Численные методы на основе метода Галеркина. -М.: Мир, 1988. -352с.
4. Settiev Sh.R. Volnoobrazovaniya nad peschanim dnom. https://www.morebooks.de/store/gb/book/Волнообразование-над-песчаным-дном^Ьп/978-3-8443-5920-6
References
1. Demidovich B. P., Maron I. A. Foundations of computational mathematics. -Moscow.: Publishing house LAN, 2006. S.
2. Roache P. Computational fluid dynamics. - M.: Mir, 1980. -612c.
3. Fletcher K., Numerical methods based on Galerkin method. -M.: Mir, 1988. -352 C.
4. Settiev Sh.R. Volnoobrazovaniya peschanim nad dnom.
https://www.morebooks.de/store/gb/book/Волнообразование-над-песчаным-
дном/isbn/978-3-8443-5920-6
