Ортогональные финитные функции для нелинейного пространственного уравнения вязких трансзвуковых течений Текст научной статьи по специальности «Математика»

70

УДК 519.63

Вестник СамГУ — Естественнонаучная серия. 2011. № 2(83)

ОРТОГОНАЛЬНЫЕ ФИНИТНЫЕ ФУНКЦИИ ДЛЯ НЕЛИНЕЙНОГО ПРОСТРАНСТВЕННОГО УРАВНЕНИЯ ВЯЗКИХ ТРАНСЗВУКОВЫХ ТЕЧЕНИЙ

© 2011 Г.Ю. Северин1

Описан алгоритм построения системы ортогональных финитных функций для начальной задачи для нелинейного пространственного уравнения вязких трансзвуковых течений. Результат проектирования на соответствующее конечномерное подпространство — система обыкновенных дифференциальных уравнений с диагональной матрицей.

Ключевые слова: метод Галеркина, ортогональные финитные функции Леонтьева, условие Стренга — Фикса.

Введение

В этой работе предлагается непосредственное построение двух систем ортогональных финитных функций вида }. , У ^ ) } . — сдвигов

функций /(х/Ъ) и у>(у/Ъ), специально адаптированных для вязкого трансзвукового уравнения, так что после проектирования на соответствующее тестовое пространство матрица Грамма диагональная. На этом примере показано, как строить такие системы ортогональных финитных функций (ОФФ) для решения достаточно широкого класса начальных задач с частными производными. Описанная ниже конструкция ОФФ имеет два существенных отличия от ОФФ В.Л. Леонтьева [5]. Построенные в этой статье функции / и удовлетворяют не одному условию ортогональности сдвигов {/(х),/(х — 1)} = 0, а нескольким условиям ортогональности, благодаря чему при проектировании на тестовое пространство исходного нелинейного уравнения нейтрализуется нелинейность иХиХХ. Платой за несколько дополнительных условий ортогональности выступает замена точного условия Стренга — Фикса на приближенное.

1. Постановка задачи

Многие нелинейные волновые процессы достаточно точно описывают вязкое трансзвуковое уравнение [2]:

1 Северин Григорий Юрьевич (akg.77@mail.ru), кафедра математического и прикладного анализа Воронежского государственного университета, 394006, Российская Федерация, г. Воронеж, Университетская пл., 1.

д2 и д3и д2и д2п , ди д2и тт/

" — ТГл — = ЧГ* + Н(x,У,z,t), (!)

дЬдх дх3 ду2 дг2 дх дх2 где а> 0, (х,у,г,Ь) € (—те, те) х (—те, те) х (—те, те) х [0, ж); Н —кусочно непрерывна по Ь, дважды непрерывно дифференцируема по х и один раз по у и Пусть неизвестная функция и(х,у, г,Ь) удовлетворяет начальному условию:

и(х, у, г, 0) = ио(х,у,г), (2)

где ио — непрерывно дифференцируема.

2. Оистемы ортогональных финитных функций для уравнения вязких трансзвуковых течений

Пусть Н фиксированное достаточно малое число — шаг равномерной сетки с узлами {хт}^, {уи}^, . Определим функции /(х) и <р(у) с носителями

[—1,1], порождающие системы функций } , } ' на основе ко-

торых методом Галеркина будет строиться приближенное решение задачи (1), (2). Пусть

/ (х) = (ах1 + а—^х1-1 + ... + а\х + 1)ир(х). (3)

Здесь ир(х) — атомарная функция В.Л. Рвачева. Отметим, что ир(х) четная, неотрицательная, бесконечно дифференцируемая функция с носителем [—1,1]. Ее график имеет характерный колоколообразный вид. Справедлива формула ир' (х) = = 2ир(2х — 1) — 2ир(2х + 1). В [7] приведены явные формулы значений атомарной функции в двоично-рациональных точках отрезка [0,1]. Неизвестные коэффициенты «1,...,«1 определим далее из условий ортогональности и Стренга — Фикса. Пусть функция / удовлетворяет следующим условиям ортогональности в смысле

Ь2 (К):

{/(х),/ (х — 1)) =0, (4)

{/'''(х — 1),/ '(х)) =0, (5)

{/'(х),/ (х — 1)) =0, (6)

{/'(х — 1)/''(х — 1),/' (х)) =0, (7)

{/'(х)/''(х — 1),/'(х)) =0. (8)

Необходимо, чтобы система сдвигов х—Х}. обеспечивала достаточную точность аппроксимации. Известно [6; 8], что если функция /(х) с носителем [—1,1] принадлежит , то эквивалентны следующие условия:

1) /(0) = 0, но / имеет нули порядка р в других точках, кратных 2п, т. е.

Ва!(2пз)=0, 0= € Z, М < р;

2) если |а| ^ р, то

Т.за/(Ь — з)

зеъ

является полиномом от ¿1,...,^ с главным членом еЬа, е = 0;

3) для каждой функции и(х) € ^Р+1(М) существуют такие коэффициенты Wj, что при Н —> 0+ выполнены неравенства

1|и — £ II™? < еН+^Ии^г+г, £ ^^ |2 < е\\и\\12, 0 < * < р,

где постоянные с3 и с не зависят от и(х) и от к. Показатель р +1 — в является наилучшим для рассматриваемых классов и .

Для точности порядка к в £2 (К) необходимо и достаточно условие разбиения единицы:

Ух е К ]Г /(х — з) = 1.

В конструкции (3) условие Ух е [0,1] : /(х) + /(х — 1) = 1, вообще говоря, невыполнимо.

Потребуем, чтобы финитная функция /(х) удовлетворяла неточному условию типа Стренга — Фикса на отрезке [0,1]:

/(х) + /(х — 1) « 1 Ух е [0,1].

Пусть

/ (х<) + / (х — 1) = 1, г = 1,...,д, (9)

где х1,...,хч — сетка (необязательно равномерная) на отрезке [0,1].

Итак, объединяя условия ортогональности (4)—(8) и приближенное условие Стренга — Фикса (9), получим алгебраическую систему уравнений для чисел ах,...,а1. Степень I нужно положить равной 5 + д. Численно решая полученную систему, в которой только 5 уравнений нелинейные (три квадратичных и два кубических), находим неизвестные коэффициенты а,1,...,щ. После чего нормируем производную:

11/'(х)||ь2 = 1 (10)

и введем обозначения:

1 1 I/"(х)Ш2 = а, I/(х)3,х = 5, ¡х/'(х)3,х = (11)

-1 -1 Функция /(х) построена (с точностью до знака). Пусть

ф(у) = (Ь у" + Ь—1У8-1 + ... + Ьх у + 1)ир(у).

Неизвестные коэффициенты Ь1,...,Ь8 определим аналогично из одного условия ортогональности сдвигов и неточных условий Стренга — Фикса:

1

У ^(у)^(у — 1)3,у = 0, (12)

о

^ — 1) = 1, г = 1,...,д, (13)

где у1,...,уч — сетка (необязательно равномерная) на отрезке [0,1].

Объединяя условия (12), (13), получим алгебраическую систему уравнений для чисел Ь1,...,Ь8. То есть степень в нужно положить равной 1 + д. Численно решая полученную алгебраическую систему, в которой только 1 уравнение нелинейное, находим неизвестные коэффициенты Ь1,...,Ь8. После этого нормируем найденную функцию:

Му)1к =1 (14)

и введем обозначение:

1

У <р(у)3у = Л. (15)

-1

Функция у>(у) построена (с точностью до знака).

3. Метод Галеркина на построенных системах ортогональных финитных функций

Приближенное решение начальной задачи (1), (2) будем искать в виде:

Атпр = 53 Атпр/тфпРр- (16)

тпр тпр

Подставим (16) в (2) и умножим скалярно в смысле ¿2 обе части полученного равенства на /^ х■ Затем полученное равенство умножим скалярно на

У н-^) ■ Результат умножим скалярно на ^ (^,3,к = -ж,..., ж). Исполь-

зуя (4), (10)-(12), (14), (15), получим: 1 1 1

А^к(0) = J ! J uo(xi + Нх,у^ + Ну,хи + Нг)/(х)ф(у)ф(г)йхйуйг =

-1-1-1

= ио(хн,Уз ,гк)5\2 + О(Н). (17)

Подставим (16) в (1) и умножим скалярно в смысле ¿2 обе части полученного равенства на /х(^- Из (5)-(8) и (10), (11) следует:

Н ^.Аъпрфпфр + А™р¥п¥р - Н1Н'х (х^ у, г, £) + О(Н2) = 0.

пр пр

Из-за условий (7), (8) на этом шаге обнулилось нелинейное слагаемое ихихх исходного уравнения (1), а благодаря (6) — слагаемое иуу + игг.

Теперь умножим скалярно обе части полученного равенства на ^У-уУ:

+ А-чр^р - Н27\Н'х X,уч ,г,t) + О(Н3) = 0.

рр

Умножим скалярно обе части полученного равенства на Н ^ :

■^А^к + а-НаА^к - Н2^Х2Н'х,гк=0. (18)

Из (17), (18) получим:

г

А^к(¿) = 5\2ио^,уч,гк)е-аа - Н2^Х2 е-аа(г-3')Н'хх,уч,гк,з)йз.

о

В зависимости от конкретной физической задачи сетки хт, уп, гр могут иметь разные шаги Н1, Н2, Нз. Построенные ОФФ / и будут теми же.

Литература

[1] Блатов И.А., Стрыгин В.В. Элементы теории сплайнов и метод конечных элементов для задач с погранслоем. Воронеж: Изд-во ВГУ, 1997.

[2] Солитоны и нелинейные волновые уравнения / Р. Додд [и др.]. М.: Мир, 1988.

[3] Лаевский Ю.М. Метод конечных элементов. Новосибирск: Изд-во НГУ, 1999.

[4] Ильин В.П. Методы и технологии конечных элементов. Новосибирск, 2007.

[5] Леонтьев В.Л. Ортогональные финитные функции и численные методы. Ульяновск: УлГУ, 2003.

[6] Агошков В.И., Марчук Г.И. Введение в проекционно-сеточные методы. М.: Наука, 1981.

[7] Рвачев В.А. Теория приближений и атомарные функции. М.: Знание, 1978.

[8] Стренг Г., Фикс Дж. Теория метода конечных элементов. М: Мир, 1977.

Поступила в редакцию 4///2011; в окончательном варианте — 4///2011.

ORTOGONAL COMPACT FUNCTIONS FOR NONLINEAR SPATIAL EQUATION OF VISCOUS TRANSONIC FLOW

© 2011 G.Yu. Sewerin2

In the article the unique algorithm of construction of the system of orthogonal compact V.L. Leontev's type functions for the nonlinear spatial equation of viscous transonic flow is described by the method of Galerkin. The result of designing on corresponding subspaces is system of ordinary differential equations with a diagonal matrix.

Key words: method of Galerkin, orthogonal Leontiev compact functions, Streng-Fix condition.

Paper received 4///2011. Paper accepted 4///2011.

2Severin Grigoriy Yurievich (akg.77Smail.ru), the Dept. of Mathematics and Applied Analysis, Voronezh State University, Voronezh, 394006, Russian Federation.