Научная статья на тему 'Сходимость в сильных нормах проекционно-разностного метода решения параболического уравнения с периодическим условием на решение'

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

CC BY
95
13
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГИЛЬБЕРТОВО ПРОСТРАНСТВО / ПАРАБОЛИЧЕСКОЕ УРАВНЕНИЕ / ГЛАДКАЯ РАЗРЕШИМОСТЬ / ПЕРИОДИЧЕСКОЕ УСЛОВИЕ / НЕЯВНАЯ СХЕМА ЭЙЛЕРА / HILBERT SPACE / PARABOLIC EQUATION / SMOOTH SOLVABILITY / PERIODIC CONDITION / IMPLICIT EULER METHOD

Аннотация научной статьи по математике, автор научной работы — Бондарев Андрей Сергеевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Бондарев Андрей Сергеевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

THE STRONG-NORM CONVERGENCE OF A PROJECTION-DIFFERENCE METHOD OF SOLUTION OF A PARABOLIC EQUATION WITH THE PERIODIC CONDITION ON THE SOLUTION

A smooth soluble abstract linear parabolic equation with the periodic condition on the solution is treated in a separable Hilbert space. This problem is solved approximately by a projection-difference method using the Galerkin method in space and the implicit Euler scheme in time. Effective both in time and in space strong-norm error estimates for approximate solutions, which imply convergence of approximate solutions to the exact solution and order of convergence rate depending of the smoothness of the exact solution, are obtained.

Текст научной работы на тему «Сходимость в сильных нормах проекционно-разностного метода решения параболического уравнения с периодическим условием на решение»

ISSN 1810-0198. Вестник Тамбовского университета. Серия: естественные и технические науки

Том 23, № 124

2018

DOI: 10.20310/1810-0198-2018-23-124-617-623 УДК 517.988.8

СХОДИМОСТЬ В СИЛЬНЫХ НОРМАХ ПРОЕКЦИОННО-РАЗНОСТНОГО МЕТОДА РЕШЕНИЯ ПАРАБОЛИЧЕСКОГО УРАВНЕНИЯ С ПЕРИОДИЧЕСКИМ УСЛОВИЕМ НА РЕШЕНИЕ

< А. С. Бондарев

ФГБОУ ВО «Воронежский государственный университет» 394018, Российская Федерация, г. Воронеж. Университетская пл., 1 Е-таЛ: [email protected]

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

Ключевые слова: гильбертово пространство; параболическое уравнение; гладкая разрешимость; периодическое условие; неявная схема Эйлера

Введение

В настоящей работе решение параболического уравнения с периодическим условием на решение находится полностью дискретным проекционно-разностным методом с использованием метода Галеркина по пространству и неявной схемы Эйлера по времени. Заметим, что энергетические оценки погрешности данного метода в условиях слабой разрешимости задачи были установлены в работе [1]. В настоящей работе представлены оценки в сильных нормах погрешности приближенного решения при более жестких ограничениях в предположениях гладкой разрешимости исходной задачи.

Отметим, что в случае, когда параболическое уравнение рассматривается с начальным условием (задача Коши), оценки в сильных нормах погрешности для проекционно-разностного метода с неявной схемой Эйлера по времени установлены в работе [2].

618

A.C. Бондарев

1. Основные понятия

Пусть даны вложенные сепарабельные гильбертовы пространства У —гН , где пространство V - двойственное к У, а пространство Н отождествляется со своим двойственным. Оба вложения плотны и непрерывны. Рассмотрим полуторалинейпую по и, V { V форму a(u7v). Пусть для u,v { V

v)\\> ц и у v V-, Re а(и, и) С а и у, (1)

где а > 0. Форма а(и, v) порождает линейный ограниченный оператор А : V Е V , такой, что (Au^v) = а(и, v), где выражение типа (z, v) есть значение функционала z { V на элементе и{ V. Если z{ Н7 то (z, v) - скалярное произведение в Н [3, гл. 2].

Для нормированного пространства X далее будем обозначать через С([0,Т],Х) пространство непрерывных функций, действующих на отрезке [О, Т], со значениями в пространстве X ; через Lp(0. T]X) - пространство измеримых на отрезке [О, Т] функций со значениями в X, суммируемых с р-й степенью по норме пространства X.

Рассмотрим в V на [О, Т] параболическую задачу:

и (t) + Au(t) = f{t), и(0) = и(Т). (2)

Здесь и далее производные функций понимаются в обобщенном смысле, В [4, с. 289] показано, что для заданного / { ¿2(0, Т; V ) существует (и притом единственное) решение и{ L2(Q,T,V)r\C([0,T\,H), и { L2(0:T;V) задачи (2).

Далее будем считать, что выполнены условия гладкой разрешимости, то есть справедлива следующая теорема [5].

Теорема 1. Пусть вложение V —>Н компактно, а форма а(и, v) удовлетворяет требованиям (1). Пусть функция t Е f(t) { V дифференцируема, f { ¿2(0,Т;У ), и выполняется равенство /(0) = f(T). Тогда решение задачи (2) будет таким, что и { L2(0,T;V)f)C([0, Т].Н), и { L2(0,T,V).

Пусть V/t, где h - положительный параметр, есть конечномерное подпространство пространства V. Определим пространство Vh: задав на элементах щ { Vь двойственную норму Uh v¿ = sup IK'tí/i, i>h) где точная верхняя граница берется по всем Vh { Vh, vh у = 1. Отметим, что uh > Щ v,.

Пусть Ph - ортопроектор в пространстве Н на V/, ■ Как замечено в [6], оператор Ph допускает расширение по непрерывности до : V Е Vh, причем для и { V справедливо Р^и v¿ > и у.

Для построения приближенных решений возьмем равномерное разбиение 0 = t0 < ti < Í2 < - - - < íjv = T отрезка [0. Т], где N { N. В подпространстве рассмот-

рим периодическую разностную задачу: для k = 1,N

(uhk uhk х)т 1 + PhAuhk = Phf(tk), uh0=uhN, (3)

где tN = T, tk = кт.

Таким образом, процесс нахождения приближенного решения задачи (2) сводится к нахождению решения конечной линейной алгебраической системы уравнений (3), однозначная разрешимость которой установлена в 111.

Далее будем предполагать, что форма а(и, v) является симметричной, то есть а(и, v) = a(v,u), где черта над комплексным числом означает переход к сопряженному числу

Определим гильбертово пространство

V(A) = }u,v { V\\u,v)v(A) = а(и, и)|. Из (1) следует, что нормы в пространствах V и V(A) эквивалентны, то есть

а1'2 и у > и v{A) > и v { V).

Лемма 1. Для и { V справедлива оценка

[I Qh(A)]u v>a ^ [I Qh(A)]u V(A) > a ^У'2 (I Qh)u у, где Qh(A) - ортопроектор в V(A) на Vh, a Qh - ортопроектор в V на Vh.

2. Основные результаты

Теорема 2. Пусть для задачи (2) выполнены условия теоремы 1 и форма а(и, v) симметрична. Пусть u(t) - решение задачи (2), обладающее дополнительной гладкостью и { Lp(0, Т; íf), где 1 > р > 2, a и^,..., u^¡) решение задачи (3). Тогда справедлива оценка погрешности

N

шах u{tk) v- + V О/ kl N —'

N

Е

к=1

■щ, U

к= 1 2

1 ftk Т Jtk-1

и (í) dt

и:,: и

к i

т+

н

Т > С jr3 2/р Qf и (í) pHd?j +

jf [I Qh(A)]u(t) (/ Qh)v(t) ^J.

(4)

Сходимость погрешности приближенных решений к нулю естественным образом следует из оценки (4). Предположим для этого, что задана предельно плотная в V при И Е 0 последовательность подпространств }Уь\ ■ Это означает, что (/ V € О

при /г € 0 для любого V { V. Тогда при те 0 и ¡г Е О

N

max u(tk) щ I- + У^ о г ш N К у ^

fc=i

1 fth

- «ю

т Jtk-1

dt

u¡: и

к i

Т+

н

N

Е

t=i

и (tk)

Ч Ч 1

Т е 0.

н

В случае, когда известна числовая характеристика предельной плотности последовательности подпространств , оценка (4) позволяет получить скорость сходимости и по пространственным переменным. Для этого введем в рассмотрение множество О ( А) = } и { V : Аи { Н\ . Пусть существует гильбертово пространство Е такое, что О (А) —>£' —»V и выполняется типичная для эллиптических операторов оценка

V Е > s Av н (V { D(A)),

(5)

где 5 С 0. Например, если параболическое уравнение в области £1 определено равномерно эллиптическим дифференциальным оператором второго порядка и краевым

<

условием Дирихле, то рассмотрим пространства: Н = V (О), У = \¥2 1 (О),

< '

Е = И^22(0)П Ш [4, с. 275]. Если же на границе области £1 задается условие Неймана, то пространства следующие: Н = V = И^(П), Е = 1У|(£1) [4, с. 276].

Пусть подпространства обладают следующим аппроксимационным свойством, типичным для подпространств типа конечных элементов (напр., [7, с. 143-144])

(I Qh)v v >rh V Е (v { Е), где константа г > 0 не зависит от v и h.

(б)

Следствие 1. Пусть выполнены условия теоремы 2. Пусть также выполнено условие (5) и решение и(1) задачи (2) такое, что и { С([0, Т]. Е). Пусть подпространства обладают свойством (6). Тогда справедлива оценка

N

шах u{tk) v- + V О/ ki N —'

N

Е

k=1

U (tk)

Щ

и

к 1

fc=l 2

1 rtk

r Jtk-1

и (t) dt

ui

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

и

к 1

T+

н

т > С

н

2/р

+

ti

и:

■и (¿) b dt + max u(t) % v 0/ и т

СПИСОК ЛИТЕРАТУРЫ

1. Бондарев A.C., Смагин В.В. Сходимость проекционно-разностного метода приближенного решения параболического уравнения с периодическим условием на решение // Вестник Воронежского государственного университета. Серия: Физика. Математика. 2014. № 2. С. 81-94.

2. Смагин В.В. Оценки в сильных нормах погрешности проекционно-разностного метода приближенного решения абстрактного параболического уравнения // Математические заметки. 1997. Т. 62. № 6. С. 898-909.

3. Обэн Ж.-П. Приближенное решение эллиптических краевых задач. М.: Мир, 1977. 384 с.

4. Лионе Ж.-Л., Мадженес Э. Неоднородные граничные задачи и их приложения. М.: Мир, 1971. 372 с.

5. Бондарев А.С. Разрешимость вариационного параболического уравнения с периодическим условием на решение // Вестник Воронежского государственного университета. Серия: Физика. Математика. 2015. № 4. С. 78-88.

6. Вайникко Г.М., Оя П.Э. О сходимости и быстроте сходимости метода Галеркина для абстрактных эволюционных уравнений // Дифференциальные уравнения. 1975. Т. 11. № 7. С. 1269-1277.

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

Поступила в редакцию 20 апреля 2018 г.

Прошла рецензирование 24 мая 2018 г.

Принята в печать 26 июня 2018 г.

Бондарев Андрей Сергеевич, Воронежский государственный университет, г. Воронеж, Российская Федерация, аспирант, кафедра функционального анализа и операторных уравнений, e-mail: [email protected]

Для цитирования: Бондарев А.С. Сходимость в сильных нормах проекционно-разностного метода решения параболического уравнения с периодическим условием на решение // Вестник Тамбовского университета. Серия: естественные и технические науки. Тамбов, 2018. Т. 23. № 124. С. 617-623. DOI: 10.20310/1810-0198-2018-23-124-617-623

DOI: 10.20310/1810-0198-2018-23-124-617-623

THE STRONG-NORM CONVERGENCE OF A PROJECTION-DIFFERENCE METHOD OF SOLUTION OF A PARABOLIC EQUATION WITH THE PERIODIC CONDITION ON THE SOLUTION

A. S. Bondarev

Voronezh State University 1 Universitetskaya St., Voronezh 394018, Russian Federation E-mail: [email protected]

Abstract. A smooth soluble abstract linear parabolic equation with the periodic condition on the solution is treated in a separable Hilbert space. This problem is solved approximately by a projection-difference method using the Galerkin method in space and the implicit Euler scheme in time. Effective both in time and in space strong-norm error estimates for approximate solutions, which imply convergence of approximate solutions to the exact solution and order of convergence rate depending of the smoothness of the exact solution, are obtained.

Keywords: Hilbert space; parabolic equation; smooth solvability; periodic condition; implicit Euler method

REFERENCES

1. Bondarev A.S., Smagin V.V. Skhodimost' proyektsionno-raznostnogo metoda priblizhennogo resheniya parabolicheskogo uravneniya s periodicheskim usloviyem na resheniye [The convergence of projection-difference method of approximate solution of parabolic equation with a periodic condition on the solution]. Vestnik Voronezhskogo gosudarstvennogo universiteta. Seriya: Fizika. Matematika - Proceedings of Voronezh State University. Series: Physics. Mathematics, 2014, no. 2, pp. 81-94. (In Russian).

2. Smagin V.V. Otsenki v sil'nykh normakh pogreshnosti proyektsionno-raznostnogo metoda priblizhennogo resheniya abstraktnogo parabolicheskogo uravneniya [Strong-norm error estimates for the projective-difference method for approximately solving abstract parabolic equations]. Mate-maticheskie zametki - Mathematical Notes, 1997, vol. 62, no. 6, pp. 898-909. (In Russian).

3. Aubin J.-P. Priblizhennoye resheniye ellipticheskikh krayevykh zadach [Approximate Solution of the Elliptic Boundary Problems]. Moscow, Mir Publ., 1977, 384 p. (In Russian).

4. Lions J.-L., Magenes E'. Neodnorodnyye granichnyye zadachi i ikh prilozheniya [Inhomoge-neous Boundary Problems and Its Applications]. Moscow, Mir Publ., 1971, 372 p. (In Russian).

5. Bondarev A.S. Razreshimost' variatsionnogo parabolicheskogo uravneniya s periodicheskim usloviyem na resheniye [The solvability of the variational parabolic equation with a periodic condition on the solution]. Vestnik Voronezhskogo gosudarstvennogo universiteta. Seriya: Fizika. Matematika - Proceedings of Voronezh State University. Series: Physics. Mathematics, 2015, no. 4, pp. 78-88. (In Russian).

6. Vaynikko G.M., Oya P.E. O skhodimosti i bystrote skhodimosti metoda Galerkina dlya abstraktnykh evolyutsionnykh uravneniy [About the convergence and the velocity of convergence

of the Galerkin's method for abstract evolutionary equations]. Differentsial'nye uravneniya - Differential Equations, 1975, vol. 11, no. 7, pp. 1269-1277. (In Russian).

7. Marchuk G.I., Agoshkov V.I. Vvedeniye v proyektsionno-setochnyye metody [Introduction to Projection-Difference Methods]. Moscow, Nauka Publ., 1981, 416 p. (In Russian).

Received 20 April 2018 Reviewed 24 May 2018 Accepted for press 26 June 2018

Bondarev Andrei Sergeevich, Voronezh State University, Voronezh, the Russian Federation, Post-Graduate Student, Functional Analysis and Operator Equations Department, e-mail: [email protected]

For citation: Bondarev A.S. Skhodimost' v sil'nykh normakh proekcionno-raznostnogo metoda resheniya parabolicheskogo uravneniya s periodicheskim usloviem na reshenie [The strong-norm convergence of a projection-difference method of solution of a parabolic equation with the periodic condition on the solution]. Vestnik Tambovskogo universiteta,. Seriya: estestvennye i tekhnicheskie nauki - Tambov University Reports. Series: Natural and Technical Sciences, 2018, vol. 23, no. 124, pp. 617—623. DOI: 10.20310/1810-0198-2018-23-124-617-623 (In Russian, Abstr. in Engl.).

i Надоели баннеры? Вы всегда можете отключить рекламу.