CC BY
23
УДК 517.957+519.6
DOI: 10.14529/mmp170209
ЛОКАЛЬНАЯ РАЗРЕШИМОСТЬ И РАЗРУШЕНИЕ РЕШЕНИЯ ОДНОГО УРАВНЕНИЯ С КВАДРАТИЧНОЙ НЕКОЭРЦИТИВНОЙ НЕЛИНЕЙНОСТЬЮ
М.О. Корпусов, Д. В. Лукьяненко, Е.А. Овсянников, A.A. Панин
Московский государственный университет имени М.В. Ломоносова, г. Москва
Рассмотрена начально-краевая задача для уравнения ионно-звуковых волн в плазме. При этом распределение Больцмана плотности электронов приближено квадратичной функцией. Для рассмотренной задачи доказана локальная (по времени) разрешимость и проведено аналитико-численное исследование разрушения решения. Методом пробных функций получены достаточные условия разрушения решения за конечное время и оценка сверху на время разрушения. В конкретных численных примерах эти оценки уточнены численно методом сгущения сеток по Ричардсону. Промежуток времени для численного счета выбирается согласно аналитически полученной оценке сверху на время разрушения решения. В свою очередь, численное моделирование уточняет момент и характер этого разрушения. В частности, показано распространение разрушения в пространстве. Таким образом, аналитическая и численная части исследования взаимно дополняют друг друга.
Ключевые слова: 'разрушение решения; нелинейная начально-краевая задача; уравнение соболевского типа; экспоненциальная нелинейность; экстраполяция по Ричардсону.
Введение
В этой работе мы продолжаем исследования нелинейных уравнений, описывающих уравнения ионно-звуковых волн в различных модельных предположениях и Н&Чс1ТЫ6 в работах [1-4]. В данной работе мы рассмотрим следующую модельную «1+1>-мерную задачу:
д2
( £2 I Uxx — £U--—и 1
uxx — £u--u + uxx = 0 при (x,t) E (0,l) x (0,T), (1)
дг2 V 2
пх{0,г) = ил(/,г) = 0 при г е [0,Т], (2)
и(х, 0) = и0(х), щ(х, 0) = и1(х) при х е [0,/]. (3)
Решение мы будем искать в классе и(х,г) е С(2)([0,Т]; С(2)[0,/]). Поэтому нам с необходимостью нужно потребовать выполнимость условий согласования начального и граничного условий
иох(0) = иох(/) = и1х(0) = Щх(1) = 0. (4)
Уравнение (1) описывает нелинейные ионно-звуковые волны в плазме с очень большим радиусом Дебая
_( кТе )1/2 _ 1 г° ={4ПёПо) , £ = г2
То/ ГЪ
Ниже приведен вывод уравнения (1) в многомерном случае. Затем мы докажем локальную во времени классическую разрешимость и получим достаточное условие разрушения решения задачи (1) - (3) за конечное время и получим оценку сверху на это время.
1. Вывод уравнения
Рассмотрим ионно-электронную плазму в области D С RN, N ^ 1. В приближении квазистационарного электрического поля имеют место следующие уравнения:
div E = —4nen, rot E = 0, (5)
где E - это вектор напряженности электрического поля, n = ne + ni это суммарная
D
ностно односвязной, можно ввести потенциал электрического поля ф, связанный с E
E = -Уф. (6)
В этом случае плотность электронов хорошо описывается распределением Больцмана
ne = no exp ^kT^J , (7)
где n0 - равновесная плотность электронов, Te - температура электронов. А для концентрации ионов ni справедливо уравнение неразрывности
ж = divJ С)
где J - это плотность тока ионной составляющей плазмы и учитывает временную дисперсию
t
J = j a(t, r)E(t) dr, a(t,r) E C(2)([0,T] 0 [0,T]). (9)
o
Из уравнений (5) - (7) мы получим следующее дифференциальное следствие:
Дф — 4nen0 exp ^= ^neni. (10)
Из системы уравнений (8) и (9) мы получим следующее дифференциальное следствие:
t
ni
= — *(t,r)Дф(т) dr. (11)
Из уравнений (10) и (11) вытекает дифференциальное следствие
t
dt (ДФ — 4neno exp ^^ + 4ne j a(t, т)Дф(т) dr = 0.
В одном частном, но важном случае a(t,r) = а0 > 0 мы приходим к следующему дифференциальному СЛ6ДСТВИЮ !
(Дф — 4neno exp ^k^ ) ) + ^^Дф = 0. (13)
д2 {Аф 4______^ вф
Сделаем в этом уравнении замену« = Пеп0, тогДа получим следующее уравнение:
д2 U ( u\\ л
— ^Дм — exp ^ rD J J + 4пeaoДм = 0,
u
11 ' o Дм = 0, (14)
ГДе ( Т ) 1/2
то
\4пе2по /
- это так называемый радиус экранировки Дебая. В случае плазмы радиус Дебая достаточно большой настолько, что можно ограничиться рассмотрением первых трех членов разложения функции ехр(еп) то малому параметру е = 1 /т0:
£-2 =
ехр(еп) = 1 + ей + — п2 + о(е2п2). (15)
Из уравнений (14) и (15) вытекает искомое уравнение:
д 2 / е2 \
— (Ап - еп - — п2] + Aп = 0, (16)
в котором мы коэффициент 4пва0 при слагаемом Ап заменили единицей.
2. Локальная разрешимость
Будем доказывать локальную во времени разрешимость задачи (1) - (3) в классе С(2)([0,Т]; С(2)[0,/]). Введем функцию Грина дифференциального оператора
-£ + е1, е >
^ . совЫл/ех) со8Ь(л/е(у — /)), при х < у;
Се(х,у) '
л/ё81пЬ^л/ё/) 1 совЬ(^у) совЬ(^(х - /)), при х > у. Перепишем рассматриваемое уравнение для удобства в эквивалентном виде
д2
е2 2
I -пТх + еп + — п I - пхх + еп = еп.
дг2
В классе гладкости С(2)([0,Т]; С(2)[0,/]) с учетом введенной функции Грина получим интегродифференциальное равенство
д 2 ( ¿2 л \
дЙ2 \ Gеп2 ) + п = еСеп, (17)
Д I
Се! = J Се(х,У)!(У) ¿у.
Проинтегрируем равенство (17) один раз по времени, тогда получим равенство
2 г
дИ^п + ~2 Сп2) = п.(х) + е2С(поп1)(х) + J ¿в п(х, в) + еСеп(в¿в, (18) ^ о
С(поп1 )(х) = J Се(х,у)по(у)п1(у) ¿у. (19)
о
Заметим, что
|(п + е2 Се^ = (I + ке.) (20)
где I - это единичный оператор, а
i
Ks(u) * w — £2G£(x,y)u(y,t)w(y,t) dy. (21)
' £ V 0
Тогда уравнение (18) можно переписать в следующем виде:
г
+ К£(и)*^ = и1(х) + £2С(иощ)(х) + ds (-и(х,в)+ £(]£и(в)^ ds. (22)
0
Рассмотрим банахово пространство В = Ь^((0,Т) х (0,/)), оде Т > 0 - произвольное фиксированное. Точно так же, как и в [1], можно доказать, что
НКДфНв^в ^ С1(0ф||, (23)
где ||-|в^в _ это операторная норма линейного пространства всех линейных непре-
ВВ
и выпуклое подмножество банахова пространства В :
В^п = {^(х,г) € В : ^ 1^п} , (24)
||ио|и~(о,0 ^ ||и1 ^ (25)
где п = п(Т) € N некоторое фиксированное, которое мы фиксируем позже, а ц = ц(Т,е) > 0 - достаточно малое число. При фиксированных Т, I, £ и п выберем ц > 0 настолько малым, чтобы
ГЛ(1)£^П ^ 2 •
В этом случае оператор + К£(и)*^ существует и представим в виде ряда Неймана
_1
(l + Kе(п)*У =Y,[-Ke(u)*]m.
m=0
Следовательно. уравнение (22) можно переписать в виде
-U = ^I + K£(u)*j I щ(х) + £2G(u0Ui) + ds (—u(x, s) + £&eu(s) j ds I . (26)
ddu — (j- + K£(u)*^ I "lV^y c ^V^O^i; I LOO \ — uyo,, O) ~T ou£
0
Далее это уравнение можно проинтегрировать один раз по времени и получить еле-дующее равенство:
t
u(x,t) — u0(x) + У A(u)(x,s) ds, (27)
0
A(u)(x,t)= (i + K£(u)*) |ui(x) + e2(G(u0u1) + J ds(- u(x,s)+ eG£u(s)) dsj . (28)
t
Точно так же, как в [1,2], можно доказать, что при выборе N Э п > 4Т + 1 - фиксированном и при достаточно малом ц > 0 имеют место следующие свойства
А(и) : Вп ^ Вп \\А(^) - Л(у2)\\ ^ 2К - ъ\\,
т.е. оператор А(-у) является сжимающим на В^п. Используя вариант принципа сжимающих отображений работы [5], можно доказать существование единственного решения и(х,Ь) Е В^п С В интегрального уравнения (27).
Далее так же, как и в [1], нужно воспользоваться «бутстэп>-методом для повышения гладкости решений интегральных уравнений и получить, что на самом деле решение интегрального уравнения (27) принадлежит классу и(х, Ь) Е С(2) ([0, Т]; С(2)([0, /])). Стало быть, нами доказана следующая теорема:
Теорема 1. Пусть фиксированы, исходные параметры задачи: £ > 0,1 > 0. Тогда для произвольного Т > 0 найдется такое ц> 0, что при достаточно малых начальных функций ио(х) Е С(2)([0,1]) и и1(х) Е С(2)([0,/]) :
\L™(0,1) ^ \\Ui\\l™(0,1) ^ удовлетворяющих условиям
uox(0) = U0x(l) = Uix (0) = Щх(1) = 0,
существует единственное решение рассматриваемой задачи в классе u(x,t) Е C(2)([0,T]; C(2)([0,l])).
3. Разрушение и оценка времени разрушения решения
Пусть u(x,t) Е C(2)([0,T]; C(2)[0,l]) - это локальное единственное решение задачи (1)-(3). В этом разделе мы докажем, что
Г 1_1
T ^ T^ = J ^u0(x)--^ dx —2J ^u0(x)--^ u1(x) dx (29)
П L П
при достаточном условии, что
У ^u0(x)--^ u1(x) dx < 0. (30)
п
Действительно, интегрируя уравнение (1) по области (0, l), с учетом граничных условий (2) мы получим следующее равенство:
i i 0 d2 [ ( £2 \ d2 f f l \2
dt2 ( £u(x,t) —2u2(x,t)J dx = 0 ^ dtt2 \u(x,t)--j dx = 0. (31)
00 Интегрируя равенство (31) по времени, с учетом начальных условий (3) получим
С Л 6ДУ ЮЩ6 G p^BGHCTBO!
i 2 i 2
0 ^ ^u(x,t)--^ dx = J ^u0(x)--^ dx + 2tJ ^u0(x)--^ u1(x) dx. (32)
0 0 п
Потребуем выполнения условия (30). Тогда сразу же приходим к искомой оценке (29), поскольку при t > Tравенство (32) становится противоречивым.
4. Численная диагностика разрушения решения
В этом разделе мы детально обсудим методы, которые помогут нам численно диагностировать факт разрушения решения и уточнить его локализацию как во времени, так и в пространстве. Напомним, что априорная информация, полученная аналитически в разделе «Разрушение и оценка времени разрушения решения>, дает нам оценку времени разрушения решения, однако не дает детального описания процесса разрушения. Численный же подход, использующий аналитически полученную априорную информацию, может помочь детализировать процесс разрушения и уточнить момент разрушения.
Будем использовать обозначения Щпи0(х),Щпил (х) вместо «о(х),«1(х) в связи с тем, что мы используем подобные индексы для определения сеточных значений функций. Момент времени разрушения обозначим через Ты.
Для начала сведем исходное уравнение (1) - (3) к системе первого порядка по времени. Это необходимо для применения эффективных численных методов, которые будут описаны ниже. Таким образом, задача
— ( - - £2 2\ + dt2 y^xx £u 2 U J + uxx
u(x, 0) = uinito(x), ut(x, 0) = uinit1 (x), ux(0, t) = ux(l, t) = 0
0, x E (0,l), t E (0,Tbi]
(33)
примет вид
x E (0,l), t E (0, Tbl]
fд ( £2 2\
~r~ uxx — £u--u = 1
0t\ 2 )
д
—v + uxx = 0, dt
u(x, 0) = uinito (x),
v(x, 0) = uinit1 (x)
xx £uiniti (x) £ uinito (x)uiniti (x),
^ ux(0, t) = ux(l,t) = 0.
(34)
5. Жесткий метод прямых и схема Розенброка с комплексным коэффициентом
С целью численного решения системы (34) мы применяем жесткий метод прямых (ЭМОЬ) [6,7] для того, чтобы свести исходную систему уравнений в частных производных к неявной системе обыкновенных дифференциальных уравнений (к сожалению, не разрешенной относительно производных^ которая может быть эффективно решена с помощью одностадийной схемы Розенброка с комплексным коэффициентом СЕСШ [8].
Сначала мы введем равномерную сетку Х^ только по пространственной переменной х с шагом к = (/ — 0)/Ж, содержащую N + 1 узлов (что соответствует N интервалам): XN = {хп, 0 ^ п ^ N : хп = 0 + пк}. Таким образом, после конечноразностной
Замечание. Мы ставим задалу нахождения решения до момента времени Ты включительно, хотя знаем, что решение в этот момент времени не существует и, более того, даже может разрушиться ранее. Это допущение связано с тем, что мы хотим диагностировать разрушение решения численно, а значит, нам необходимо найти численное решение вплоть до этого момента времени включительно.
аппроксимации пространственных производных со вторым порядком точности в (34) мы получим следующую систему обыкновенных дифференциальных уравнений, из которой требуется определить N — 1 неизвестных функций ип = ип(Ь) = и(хп,Ь) (п = — 1, и0 и им определяются из соотношений, заданных граничными усло-
ио
4и1 — 1 и2, им = 3им-1 — 3им-2 ) и N — 1 вспомогательных функций
уп = Уп(Ь) = у(хп, Ь) (п = 1, N — 1, у0 и ум не входят в систему):
(ип—1(п , 2 , -21-2. \(ип , (ип+1
(И (ип
— (2 + £к2 + £2Н2ип) +
ип+1 — 2ип + ип— 1
к Уп
к2
ип (0) игп%1о (хп),
Уп (0)
(хп+1) — 2и (хп ) + игпИ1 (хп— 1)
к2
£игпИ\ (хп) £ игпгЬо (хп)игпИ1 (хп) ,
1
и0
—и1--и2,
3 1 3 2'
им = 3 им — 1 — 3 им—2.
Эта система может быть переписана в виде
М(ч) (Ь =/(«),
и(0) = ЩпЦ,
где и = и и2 из ... им—1 У1 У2 Уз ... Ум — 1)Т, f= (/1 /2 /э ... ¡2М—2) ИгпИ = (и1(0) и2(0) из(0) . . . им—1(0) У1(0) У2(0) Уз(0) ... Ум— 1(0)2Т Вектор-функция / имеет следующую структуру:
к2Уп,
(35)
(36)
Т
/п
—(и2 — 2и1 + (3и1 — 3щ)^ ,
— (ип—м+2 — 2ип—м+1 + ип—м^ ,
— /г? ( (зим—1 — 1 им—2) — 2им—1 + им—2
если п = 1, N — 1, если п = N.
если п = N + 1, 2N — 3, если п = 2N — 2,
а матричная функция М имеет следующие элементы:
Мп
Мп,п— 1
М
п,п+1
'з — (2 + £к2 + £2к2 и 1)
— (2 + £к2 + £2к2ип2, (2
п п
1
2,^- 2,
4
з 1,
1,
— (2 + £к2 + £2к2им—^, если п = N — 1,
если п = N 2N — 2,
п
2, N — 2,
— з + 1, если п = N — 1,
— 1 + 1, если п = 1,
1, если п = 2, N — 2,
остальные элементы матрицы М равны нулю.
4
4
1
Для численного решения системы (36) мы будем использовать схему Розенбро-ка с комплексным коэффициентом (СШЭ81), которая является наилучшим выбором для решения задач такого рода по причине порядка точности этой схемы (0(тее монотонности и устойчивости (Ь2) [9].
Важно отметить, что в случае постоянной матрицы М схема Розенброка с комплексным коэффициетом имеет порядок точности 0(т2) [10,11]. Неявная система обыкновенных дифференциальных уравнений с такой матрицей может быть получена, если в системе (34) ввести следующую вспомогательную переменную: w =
„2
ихх — еи — у и . В этом случае мы получим неявную систему обыкновенных дифференциальных уравнений (подобную (35)) с постоянной (но вырожденной) матрицей, размерность которой будет в 3/2 раза больше исходной, что приведет к увеличению времени вычислений. Какой из этих подходов более эффективен с вычислительной точки зрения - открытый вопрос. Нашей же основной целью является демонстрация того, как можно численно диагностировать факт разрушения решения по времени п п. ш пространству, - таким образом, в данной работе мы не исследуем поднятый в этом замечании вопрос.
Для того чтобы применить эту схему, мы введем равномерную сетку (также можно использовать и квазиравномерную сетку без каких-либо изменений в алгоритме) Тм по которая имеет М + 1 узлов (то есть М интервалов): Тм = {Ьт, 0 ^ т ^ М : 0 = Ь0 < <Ь2 < ... < Ьм-1 <Ьм = Т^}, где напомним, что Тж - это верхняя оценка времени разрушения (см. (29)).
После этого мы можем применить схему СНОЭ! для решения системы (36):
m(Wi) = U(tm) + (tm+1 - tm) Re W .
где w является решением системы
M
(u(tm)^ - (tm+1 - tm)fu (u(tm)^ W = f (u(tm)^
(37)
Здесь /и - якобиан, который для рассматриваемой системы имеет следующую структуру:
(fu)n,n+N -1 (fu)n,n-N
(fu)n,n-N +1 (fu)n,n-N+2
h , если n = 1, N — 1,
!
h2, __L + Л.
h2 + 3h2,
3h2'
если n если n
N + 1, 2N — 3, 2N — 2,
W — 3h2 n = N,
h2,
h2 — 3h2 ^лп = 2N — 2,
если n = N + 1, 2N — 3,
. h2
+ —
h2 + 3h2' h2,
n
N
если n = N + 1, 2N — 3.
Другие компоненты матрицы ^ в случае рассматриваемого уравнения равны нулю.
Таким образом, матрица системы (37) состоит из четырех блоков размерности (М — 1) х (М — 1) (структура матрицы представлена на рис. 1). Этот факт дает возможность применить алгоритм решения СЛАУ, который найдет решение системы (37) за 0(М) операций.
Рис. 1. Структура матрицы СЛАУ (37) 6. Вычисления на сгущающихся сетках
При численных расчетах важно не только получить приближенный численным результат, но также и выполнить некоторую оценку его точности. Метод вычисления апостериорной асимптотически точной оценки погрешности [10] позволяет это сделать. Но этот метод также может помочь и диагностировать факт разрушения точного решения [11]. Основные формулы и утверждения этого параграфа впервые были представлены в работах [10 12].
Мы аппроксимировали все пространственные производные в (34) с точностью 0(к2), а при численном интегрировании системы (36) используем схему СКОЭ1, которая имеет точность 0(тПоэтому построенный метод решения системы (34) имеет точность 0(т1 + к2).
Для начала введем базовую сетку Х^ х Тм-. {хп, 0 ^ п ^ М, 0 ^ т ^ М. После этого произведем последовательное сгущение сетки, начиная с базовой, и вычислим решения и(х,Ь) на полученных сетках. В связи с тем, что теоретический порядок точности по времени равен 1, а по пространству 2, мы выполним последовательное сгущение сетки по времени в целое число раз тг и сгущение пространственной сетки в целое число раз тх так, чтобы выполнялось условие т1 = тХ (подробности см. в [12]). Наиболее удобно /игл я с ч е тсх выбрать г = 4 и Тх = 2. В этом случае каждая последующая сетка ХГа-1И х ТГэ-1М (^ - номер сетки) имеет узлы, совпадающие с узлами базовой сетки (хп,1т). В этих узлах (х,£) мы можем выполнить апостериорную асимптотически точную оценку погрешности [11,12]
п(г%М,т!М)(х уЛ _ п(гех-1М,гег-1М)(х уЛ
ДГ )(х,1) = ^-' (Х^ 1 П-!-(М + 0(т1 + к2)
Тг - 1
и оценить эффективный порядок точности [11,12]
р6/1 (х,г) = 1о&
и(гХ'^г1 м) (х,г) - и(гх"м'гг2м)(х,1)
Гг _ (Га N ) / .. и(Г^-1^ Г'3-11
п(г|М,Г?М)(х,г) - и(г**,гГ1 м) (х, ¿)
В точках (х,£), в которых решение исходной задачи имеет непрерывные первые производные по времени и вторые по пространству, имеет место сходимость
р1 (х,1) -► рШеог = 1, (38)
и соответствующая оценка погрешности является асимптотически точной при s ^ ж (или, что то же самое, N,M ^ ж). Нарушение этой сходимости (38) говорит о потере гладкости точного решения. В частности, в случае степенной <сингулярности> u(x,t) ~ (t* — t)-ß для любого t > t* эффективный порядок точности peff (x,t) -> — ß. Это позволяет нам найти соответствующую степень ß. Если peff (x,t) -> —ж для любого t > t*, мы можем утверждать, что решение
экспоненциально возрастает, т.е. u(x,t) ^ ж; если Peff (x,t) -> 0 для любого
t > t*
u(x,t) ~ ln(t* — t). Момент разрушения решения t* может быть найден с точностью до величины интервала базовой сетки по времени.
Если нарушение гладкости решения возникает во всей области по пространственной переменной одновременно, то отклонение сходимости peff (x,t) от 1 возникает
во всех точках сетки {xn} с первого временного слоя t ^ t* (см. Пример 1а). Если
x*
позволяет проследить во времени процесс разрушения решения в остальных точках (см. Пример 16 и Пример 1в). Такая диагностика процесса разрушения решения возможна в связи с тем, что схема CROS не приводит к переполнению даже в том случае, если решение задачи устремляется к бесконечности [10,11].
Пример 1а. Для начала рассмотрим пример со следующим набором
вход н ых
данных: Uinito (x) = 0, Uiniti (x) ^ 1, l — n, s — 1.
В данном случае решение уравнения (33) может быть найдено аналитически:
u(x, t) — Vi — 2t — 1. (39)
Очевидно, что в данном случае время разрушения равно Ты — Tж — 0, 5. Кроме того, это подтверждается непосредственно
формулой (29).
Применим численный алгоритм диагностики момента и места разрушения решения, а заодно и сверим получившиеся результаты с аналитическими расчетами. Для численного решения задачи (33) мы возьмем следующий набор параметров: Тх — 0, 5, N — 50, M — 200, rx — 2, rt — 4 S — 6 (число последовательно используемых для вычислений сеток, включая начальную).
Получив приближенное численное решение на разных сетках, мы можем проверить сходимость эффективного порядка точности к теоретическому для каждого временного слоя по формуле
efn , 1 y/TNo (u(rX-lN'rt-1M)(xn,tm) — ui^r^^tm))2 Р7} (tm) — logrt " ■--, (40)
Venu (u(rx N'rsM )(xn,tm) - u(rx-iN'r't-iM )(xn,tm))2
где s - номер соответствующей сетки. Также еще раз хотим обратить внимание, что под обозначением (xn,tm) подразумеваются узлы, совпадающие с соответствующими узлами базовой сетки. После вычислений на S вложенных сетках эффективный порядок точности peff для каждого временно го слоя tm сходится к ptheor = 1 (см. рис. 2) за исключением временного слоя, соответствующего моменту времени t2oo' peff (t2oo) 0, 5, что озтачает, что Ты £ (tigg,t2oo] = (0, 498, 0, 5] является временем разрушения решения, и что в точке Ты решение имеет степенной характер роста u(x,t) ~ (Ты — t)0'5. Мы также можем оценить эффективный порядок точности для
каждой пространственной точки выбранного временного слоя, соответствующего моменту времени tm, по формуле
PХ t ) - W uUr-NTLM'(Xn-tm) - nUr-NT-M'(Xn.tm)
Можно использовать эту формулу, например, для временного слоя, соответствующего моменту времени t200> в котором было диагностировано разрушение решения, с целью определить, возникло ли разрушение на всем временном слое или только в отдельных точках по пространственной переменной. На рис. 3 мы можем видеть, что решение разрушилось во всех пространственных точках соответствующего временного слоя. Тем не менее, как мы видим на рис. 2, «точной> сходимости peff к 0,5 нет. В данном случае этого можно добиться, увеличив число сеток S, на которых будут выполняться вычисления.
i
0.5 -
f 0 -
-0.5 -
-i -
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 t
Рис. 2. Пример 1а. Эффективный порядок точности для каждого временного слоя. Разрушение решения диагностировано в момент времени Ты Е (¿199,^200] = (0, 498, 0, 5]. Здесь и далее нулевой временной слой не показан
Рис. 3. Пример 1а. Эффективный порядок точности для каждой пространственной точки временного слоя, соответствующего моменту времени t2oo = 0, 5, на котором был диагностирован факт разрушения решения. Разрушение решения диагностировано во всех пространственных точках этого временного слоя
Пример 16. Рассмотрим пример со следующим набором данных. Uinit0 (x) = 0, uinitl (x) = - (x(n — x))2 sin x, l = П, £ = 1010.
В данном случае решение уравнения (33) уже не может быть найдено аналитически. Для численной диагностики разрушения решения во времени и в пространстве опять применим описанный алгоритм. Верхняя оценка Т^ для времени разрушения Ты может быть вычислена по формуле (29): Тж ~ 3.14 • 10_11. Для численного решения задачи (33) мы будем использвать следующий набор параметров: Т^ = 3,14 • 10_11, N = 50, М = 50, гх = 2, VI = 4, Б =5 (число сеток, использовавшихся для вычислений; включая исходную). Так же, как и в Примере 1а. мы можем оценить эффективный порядок точности по формуле (40). После вычислений на Б вложенных сетках эффективный порядок точности р е^ для каждого временного слоя ¿т сходится к рЬНеог = 1 (см. рис. 4) до временного слоя, соответствующего моменту времени ¿25. Это означает, что Ты Е (¿24,£25] = (1, 44,1, 51] • 10_11 является временем разрушения решения. Мы также можем использовать формулу (41) для временного слоя, который соответствует моменту времени ¿25 и на котором было диагностировано разрушение решения, с целью определить, возникло ли разрушение на всем временном слое или только в отдельных точках по пространственной переменной. На рис. 5 мы можем видеть, что решение разрушается в одной пространственной точке соответствующего временного слоя, а затем взрыв распространяется постепенно на остальные точки отрезка.
Рис. 4. Пример 16. Эффективный порядок точности для каждого временного слоя. Разрушение решения диагностировано в точке Tbi Е (W25] = (1, 44,1, 51] • 10-11
Пример 1в. Рассмотрим теперь аналогичную Примеру 16 задачу, но с более большим е uinit0(ж) = 0, Uinitl (ж) = — (х(п — x))2 sin |, l = п, £ = 1015.
В этом случае решение уравнения (33) также не может быть найдено аналитически. Применим вновь описанный вьтттте алгоритм для разрушения решения во времени и в пространстве. Верхнюю оценку T^ для времени разрушения Tb¡ вычисляем по формуле (29): Tх ~ 3,14 • 10-16. Для численного решения задачи (33) мы будем использвать следующий набор параметров: T= 3,14 • 10-16, N = 50, M = 200, rx = 2, rt = 4 S = 5 (число сеток, использовавшихся для вычислений, включая исходную).
Так же, как и в предыдущих двух примерах, мы можем оценить эффективный порядок точности по формуле (40). После вычислений на S вложенных сетках эффективный порядок точности peff для. каждого временно го слоя tm сходится к ptheor = 1 (см. рис. 6) до временного слоя, соответствующего моменту времени t98. Это означает, что Tbi Е (t97,tg8] = (1, 51,1, 52] • 10-16 является временем разрушения решения. Используем формулу (41) для временного слоя, который соответствует моменту
Рис. 5. Пример 16. Эффективный порядок точности для каждой пространственной точки для временных слоев (¿24, ¿25, ¿26, ¿27, ¿28, ¿39, ¿50). Здесь мы видим, что сначала разрушение происходит в одной пространственной точке, а потом оно распространяется на остальные точки отрезка
Рис. 6. Пример 1в. Эффективный порядок точности для каждого временного слоя. Разрушение решения диагностировано в точке Ты Е (¿97, ¿98] = (1, 51,1, 52] • 10"16
времени ¿98 и на котором было диагностировано разрушение решения, с целью определить, возникло ли разрушение на всем временном слое или только в отдельных точках по пространственной переменной. На рис. 7 мы можем видеть, что решение разрушается в одной пространственной точке соответствующего временного слоя, а затем взрыв распространяется постепенно на остальные точки отрезка.
Рис. 7. Пример 1в. Эффективный порядок точности для каждой пространственной точки для временных слоев (¿97,¿98,¿110,¿150,¿180,¿200)- Здесь мы видим, что сначала разрушение происходит в одной пространственной точке, а потом оно распространяется на остальные точки отрезка
Литература
1. Корпусов, М.О. Критические показатели мгновенного разрушения или локальной разрешимости нелинейных уравнений соболевского типа ./' М.О. Корпусов /'/ Известия РАН. Серия математическая. - 2015. -Т. 79, № 5. - С. 103-162.
2. Korpusov, М.О. Blow-up for One Sobolev Problem: Theoretical Approach and Numerical Analysis / M.O. Korpusov, D.V. Lukyanenko, A.A. Panin, E.V. Yushkov /'/ Journal of Mathematical Analysis and Applications. - 2016. - V. 442, № 2. - P. 451-468.
3. Korpusov, M.О. Blow-up Phenomena in the Model of a Space Charge Stratification in Semiconductors: Analytical and Numerical Analysis / M.O. Korpusov, D.V. Lukyanenko, A.A. Panin, E.V. Yushkov // Mathematical Methods in the Applied Sciences. - 2017. -V. 40, № 7. - P. 2336-2346.
4. Корпусов, M.O. О разрушении за конечное время решения начально-краевой задачи для нелинейного уравнения ионно-звуковых волн / М.О. Корпусов // Теоретическая и математическая физика. - 2016. - Т. 187, № 3. - С. 447-454.
5. Панин, А.А. О локальной разрешимости и разрушении решения абстрактного нелинейного интегрального уравнения Вольтерра / А.А. Панин // Математические заметки. -2015. - Т. 97, № 6. - С. 884-903.
6. Hairer, Е. Solving of Ordinary Differential Equations / E. Hairer, G. Wanner. - Berlin; Heidelberg: Spinger-Verlag, 2002.
7. Калиткин, H.H. Численные методы решения жестких систем / Н.Н. Калиткин // Математическое моделирование. - 1995. - Т. 7, № 5. - С. 8-11.
8. Rosenbrock, Н.Н. Some General Implicit Processes for the Numerical Solution of Differential Equations / H.H. Rosenbrock // The Computer Journal. - 1963. -V. 5, № 4. - P. 329-330.
9. Алышш, А.Б. Схемы Розенброка с комплексными коэффициентами для жестких и дифференциально-алгебраических систем / А.Б. Алынин, Е.А. Алышша, Н.Н. Калиткин, А.Б. Корягина // Журнал вычислительной математики и математической физики. - 2006. - Т. 46, № 8. - С. 1392-1414.
10. Алышша, Е.А. Диагностика особенностей точного решения при расчетах с контролем точности / Е.А. Алышша, Н.Н. Калиткин, П.В. Корякин // Журнал вычислительной математики и математической физики. - 2005. - Т. 45, № 10. - С. 1837-1847.
11. APshin А.В. Numerical Diagnosis of Blow-up of Solutions of Pseudoparabolic Equations / A.B. Al'shin, E.A. Al'shina // Journal of Mathematical Sciences. - 2008. - V. 148, № 1. -C. 143-162.
12. Вычисления на квазиравномерных сетках / Н.Н. Калиткин, А.Б. Алышш, Е.А. Алыни-на, Б.В. Рогов. - М.: Физматлит, 2005.
Максим Олегович Корпусов, доктор физико-математических наук, профессор, кафедра математики, физический факультет, Московский государственный университет имени М.В. Ломоносова (г. Москва, Российская Федерация), korpusov@gmail.com.
Дмитрий Витальевич Лукьяненко, кандидат физико-математических наук, доцент, кафедра математики, физический факультет, Московский государственный университет имени М.В. Ломоносова (г. Москва, Российская Федерация), lukyanenko@physics.msu.ru.
гьвгении Алексеевич Овсянников, студент, кафедра математики, физический факультет, Московский государственный университет имени М.В. Ломоносова (г. Москва, Российская Федерация), evg.bud@yandex.ru.
Александр Анатольевич Панин, кандидат физико-математических наук, доцент, кафедра математики, физический факультет, Московский государственный университет имени М.В. Ломоносова (г. Москва, Российская Федерация), a-panin@yandex.ru.
Поступила в редакцию 1 марта 2017 г.
MSC 35Q60, 35G31, 65L04, 65L12 DOI: 10.14529/mmp170209
LOCAL SOLVABILITY AND DECAY OF THE SOLUTION OF AN EQUATION WITH QUADRATIC NONCOERCIVE NONLINEATITY
M.O. Korpusov, D.V. Lukyanenko, E.A. Ovsyannikov, A.A. Panin
Lomonosov Moscow State University, Moscow, Russian Federation, E-mail: korpusov@gmail.com, lukyanenko@physics.msu.ru, evg.bud@yandex.ru, a-panin@yandex.ru
An initial-boundary value problem for plasma ion-sound wave equation is considered. Boltzmann distribution is approximated by a quadratic function. The local (in time) solvability is proved and the analitycal-numerical investigation of the solution's decay is performed for the considered problem. The sufficient conditions for solution's decay and an upper bound of the decay moment are obtained by the test function method. In some numerical examples, the estimation is specified by Richardson's mesh refinement method. The time interval for numerical modelling is chosen according to the decay moment's analytical upper bound. In return, numerical calculations refine the moment and the spacetime pattern of the decay. Thus, analytical and numerical parts of the investigation amplify each other.
Keywords: blow-up; nonlinear initial-boundary value problem; Sobolev type equation; exponential nonlinearity; Richardson extrapolation.
References
1. Korpusov M.O. Critical Exponents of Instantaneous Blow-up or Local Solvability of NonLinear Equations of Sobolev Type. Izvestiya: Mathematics, 2015, vol. 79, no. 5, pp. 955-1012. DOI: 10.1070/IM2015v079n05ABEH002768
2. Korpusov M.O., Lukyanenko D.V., Panin A.A., Yushkov E.V. Blow-up for One Sobolev Problem: Theoretical Approach and Numerical Analysis. Journal of Mathematical Analysis and Applications, 2016, vol. 442, no. 2, pp. 451-468. DOI: 10.1016/j.jmaa.2016.04.069
3. Korpusov M.O., Lukyanenko D.V., Panin A.A., Yushkov E.V. Blow-up Phenomena in the Model of a Space Charge Stratification in Semiconductors: Analytical and Numerical Analysis. Mathematical Methods in the Applied Sciences, 2017, vol. 40, no. 7, pp. 2336-2346. DOI: 10.1002/ mma.4142
4. Korpusov M.O. The Finite-Time Blowup of the Solution of an Initial Boundary-Value Problem for the Nonlinear Equation of Ion Sound Waves. Theoretical and Mathematical Physics, 2016, vol. 187, no. 3, pp. 835-841. DOI: 10.1134/S0040577916060040
5. Panin A.A. On Local Solvability and Blow-up of the Solution of an Abstract Nonlinear Volterra Integral Equation. Mathematical Notes, 2015, vol. 97, no. 6, pp. 892-908. DOI: 10.1134/S0001434615050247
6. Hairer E. Solving of Ordinary Differential Equations. Berlin, Heidelberg, Spinger-Verlag, 2002.
7. Kalitkin N.N. [Numerical Methods for Solving Stiff Systems]. Matematicheskoe modelirovanie, 1995, vol. 7, no. 5, pp. 8-11. (in Russian)
8. Rosenbrock H.H. Some General Implicit Processes for the Numerical Solution of Differential Equations. The Computer Journal, 1963, vol. 5, issue 4, pp. 329-330. DOI: 10.1093/comjnl/5.4.329
9. Al'shin A.B., Al'shina E.A., Kalitkin N.N., Koryagina A.B. Rosenbrock Schemes with Complex Coefficients for Stiff and Differential-Algebraic Systems. Computational Mathematics and Mathematical Physics, 2006, vol. 46, no. 8, pp. 1320-1340. DOI: 10.1134/S0965542506080057
10. Al'shina E.A., Kalitkin N.N., Koryakin P.V. Diagnosis of Singularities of an Exact Solution in Computations with Accuracy Control. Computational Mathematics and Mathematical Physics, 2005, vol. 45, no. 10, pp. 1769-1779.
11. Al'shin A.B., Al'shina E.A. Numerical Diagnosis of Blow-up of Solutions of Pseudoparabolic Equations. Journal of Mathematical Sciences, 2008, vol. 148, no. 1, pp. 143-162. DOI: 10.1007/sl0958-007-0542-2
12. Kalitkin N.N., Al'shin A.B., Al'shina E.A., Rogov B.V. Vychisleniya na kvaziravnomernykh setkakh [Calculations on Quasi-Uniform Grids]. Moscow, Fizmatlit, 2005.
Received March 7, 2017
