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

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

Е.А. Омельченко1, М.В. Плеханова2, П.Н. Давыдов3

Для линеаризованной квазистационарной системы уравнений фазового поля с запаздыванием представлен численный метод решения, исследована сходимость явной разностной схемы, учитывающей эффект запаздывания в исследуемой системе. На основе полученных результатов осуществлена программная реализации метода.

Ключевые слова: уравнение соболевского типа, квазистационарная система уравнении фазового поля, разностная схема.

1. Введение

В настоящей работе предложен численный метод решения линеаризованной квазистацио-

нарной системы уравнений фазового поля [1] с запаздыванием

Vt(X,t) = Av(X,t) - Aw(X,t) + ФіУ (X,) + Фі2ww(X,), (X,t)є [0,n]x[0,T], (1)

0 = v(X, t) + (fi + A)w(X, t) + Ф21V (X, •) + Ф22wf (X, •), (X, t) є [0, n] x[0,T], (2)

v(0, t) = v(n, t) = w(0, t) = w(n, t) = 0, t є [0,T], (3)

v(X, t) = ^(X, t), w(X, t) = y(X, t), (X, t) є [0,n] x[-r,0], (4)

где V (X, я) = у(X, t + 5), wt (X, 5) = w(X, t + 5) При 5 е[-г ,0] , г > 0. При этом отображения

Фг1: V(х,•) ^ га(х,*), Фг2 : wt(х,•) ^ 2(х,*) при каждом хе[0,^], *е[0,Т] линейно и непре-

рывно действуют из пространства С ([—г ,0]; Я) в Я.

Исследование разрешимости этой задачи в рамках начальной задачи

и ^) = ¿(0, t е[-г ,0], (5)

для операторно-дифференциального уравнения соболевского типа с запаздыванием

Ьи( 0 = Ми (0 + Фи + /(0, t е[0, Т ], (6)

было проведено ранее в работах В.Е. Федорова и Е.А. Омельченко [2, 3]. Здесь и, Г - банаховы пространства, ^ (я) = и^ + я) при 5 е [-г ,0], операторы Ь : и ^ Г , Ф: С ([-г, 0] ;и) ^ Г линейны

и непрерывны, кегЬ ^{0], оператор М: ^шМ^ Г линеен, замкнут и плотно определен в и, /: [0,Т] Г . Особенность линейного эволюционного уравнения с запаздыванием (6) в том, что

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

Численным аспектам исследования задач для уравнений с последействием, в том числе задач для функционально-дифференциально-алгебраических уравнений, которые относятся к классу уравнений вида (6), посвящены работы В.Г. Пименова и его учеников [4-6]. В этих работах, в частности, сконструировано семейство сеточных методов для численного решения эволюцион-

1 Омельченко Екатерина Александровна - старший преподаватель, кафедра гуманитарных и социально-экономических дисциплин, Уральский филиал «Российская академия правосудия».

E-mail: omelchenko_ea@mail.ru

2 Плеханова Марина Васильевна - доцент, кафедра дифференциальных и стохастических уравнений, Южно-Уральский государственный университет.

E-mail: mariner79@mail.ru

3 Давыдов Павел Николаевич - аспирант, кафедра математического анализа, Челябинский государственный университет.

E-mail: davydov@csu.ru

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

С помощью результатов работ [2, 3, 5, 6] авторами данной статьи исследована сходимость явной разностной схемы для задачи (1)-(4). В первом параграфе, следуя идее разделения конечномерной и бесконечномерной фазовых составляющих, показана сходимость сеточного метода для соответствующей задачи без запаздывания (Фу = 0,1, у = 1, 2). Во втором параграфе доказана

сходимость разностной схемы, учитывающей запаздывание. При этом соответствующая схеме дискретная модель определяется стартовыми значениями, формулой продвижения на шаг и оператором интерполяции и поэтому относится к классу моделей, исследованному в [5, 6] в связи с рассмотрением невырожденных эволюционных уравнений. И, наконец, третий параграф посвящен конкретной программной реализации метода.

2. Линеаризованная квазистационарная система уравнений фазового поля

Рассмотрим начально-краевую задачу для системы уравнений

где в< 0, а V, w - искомые функции. Заметим, что начальное значение w(х,0) для одной из искомых функций w не задано. Однако, как показано в [7], задача (7)-(10) однозначно разрешима. В случае же задания обеих начальных функций V(х,0), w(х,0) задача оказывается переопределенной и необходимо выполнение условий согласования данных задачи (7)-(10) для ее разрешимости.

Разобьем отрезок пространственной переменной [0,п] на части с шагом к = п/N, определив тем самым точки х п = пк, п = 0,...,N. Аналогичным образом разобьем временной отрезок [0,Т] на части с шагом т > 0, получив точки разбиения = тт, т = 0,...,М. Приближенные значения

функций V, w в узлах с координатами (хп, tm) будем обозначать через vm , wrm . Рассмотрим сеточный метод

через v(хп, tm), w(хп, tm) обозначены истинные значения решения V, w задачи (7)-(10) в соответствующих точках. Будем говорить, что невязка имеет порядок тР1 + кР2, если существует такая

V((х, t) =Дv(х, t) -Дw(х, t), (х, t) е [0,п] х[0,Т],

0 = v(х, t) + (в + Д)w(х, t), (х, t) е [0, п] х [0,Т], v(0, t) = v(п, t) = w(0, t) = w(п, t) = 0, t е [0,Т], v( х,0) = ^( х), х е[0,п],

(7)

(8) (9)

(10)

(11)

(12)

где п = 1,...,N — 1, с начальными условиями

=9(хп), п = 0,..., N,

и граничными условиями

Невязкой метода (11), (12)

£т = ^хп , ^^т+1) — V(хп , tm ) — v(хп +1, tm ) — Мхп , tm ) + ^хп—1, tm ) ^п _ ,2

+

+

^хп+1, tm ) — 2w(хп , tm ) + w(хп—1, tm ) к2

Пт = ^ хп , tm ) + ^4 хп , tm ) +

^хп+1, tm ) — 2w(хп, tm ) + w(хп—1, tm ) к 2

константа C, не зависящая от т и h, что

< C (тр + hP2) для всех,

п = 1, ...,N — 1, т = 0,...,М — 1.

Лемма 1. Пусть точное решение V, w задачи (7)-(10) таково, что функция V дважды непрерывно дифференцируема по t, функции V, w четырежды непрерывно дифференцируемы по х. Тогда невязка метода (11), (12) имеет порядок т + к2.

Доказательство. С помощью тейлоровского разложения функций V(х,t), w(х,t) получим выражения для невязки

1

1

1

(*т 12/ \ т 12

Ьп =—2Т +12к {ухххх +Wхххх ^ Пп =—12к wхххх. С

Исследуем устойчивость этой схемы методом разделения переменных.

Теорема 1. Пусть в < 0. Тогда схема (11), (12) устойчива, если выполнено условие т < к2. Доказательство. Обозначив через р™, (Г^ коэффициенты д-х гармоник на т-м слое, имеем

v(xn,tm) = pqeiqxn, w(xn ,tm) = aqe'qxn, ,

Подстановка этих функций в (11), (12) приведет к равенствам

- = 0, ± 1, ± 2,..

РР+ Рр =^(2pm (cosqh -1) -2&m (cosqh -1)):

(13)

2^m

0 =pqm +Pamq + -f-(cos qh -1).

h2

Из последнего равенства получим выражение

pm

q

aq =Т

h2

(1 - cos qh) -в

и подставим его в (13), тогда

Pqm+1 =pq

^ 2Т f 1--------(1 - cos qh) 1 --

Из того, что в < 0, следует неравенство

2Т

rq =1 - —2(1 - cos qh)

1 -

2(1 - cos qh) - h в

< 1

2(1 - cos qh) - h в

(14)

(15)

сразу для всех д е N . Покажем, что гд > — 1 при любом д е N . Для этого надо показать, что

—~(1 - cos qh)

h

Это неравенство выполняется при т < h2 . В силу равенств (14) и (15)

2(1 - cos qh) - h2в

< 1.

^_m+1

°q =T

P

m+1

Pq

f 2T f

1 -^2(1 - cos qh) 1 -

2(1 - cos qh) - h2в

J J

(

= oq,

(1 - cos qh) -в

(

1 - 2^T(1 - cos qh)

(1 - cos qh)-в

1

1

(1 - cos qh) -в

Поэтому коэффициент перехода для второй неизвестной функции также равен rq и поэтому рассматриваемая разностная схема является устойчивой прит < h2 . □

Замечание 1. Анализ доказательства теоремы 1 приводит к выводу, что при в — 0 не представляется возможным доказать устойчивость предложенной разностной схемы и необходима ее модификация или выбор другого метода вычисления.

В качестве итога полученных результатов сформулируем теорему о сходимости, которая сразу следует из критерия Куранта [8], леммы 1 и теоремы 1.

Теорема 2. Пусть в < 0 и для разностной схемы (11), (12) т, h таковы, что т < h2 . Тогда разностное решение сходится к точному решению задачи (7) - (10) с порядком не ниже т + h2.

3. Система уравнений с запаздыванием

Вернемся к начально-краевой задаче с запаздыванием (1)-(4). Для отрезка пространственной переменной [0,п] имеем прежнее разбиение с шагом h = п/N и точками x n = nh, n = 0,...,N .

Временной отрезок теперь имеет вид [-r,Т]. Для определенности считаем, что Т/r - рациональное число, шаг разбиения т = Т /M = r / L , где L, M е N, точки разбиения t m = тт, m = -L,...,0,...,M. Помимо приближений vm, wm искомых функций v, w понадобится также дискретная предыстория в точке xn к моменту tk :

{(Щ, wm )={(, wm ) - L < m < k}, n = 0,..., N, к = 0,..., M.

Предполагается также, что задан оператор интерполяции дискретной предыстории

I: {(m, wm , h^ )е Q [-r,0]x Q [-r ,0]. Здесь n = 0,..., N, Q [-r ,0] - множество кусочно-

непрерывных функций на [-r,0] с конечным числом точек разрыва первого рода в точках разрыва непрерывных справа. Зададим норму: ||g||Q[- 0] = sup |g(s)| для gе Q[-r,0]. Будем считать,

"Q[ r,0] ,e[-r,0]'

что линейные операторы Ф^ продолжимы на пространство Q [-r,0] ограниченным образом.

Следуя работам [4, 5], будем говорить, что оператор интерполяции имеет порядок погрешности тр на точном решении v, w, если существуют такие положительные константы Сь С2,что для всех n = 0, ..., N, к = 0,...,M и t е [tk - r, tk ] выполняются неравенства

gn (t) - v(xn, t) < Q max v'm - v(xn, tm) + C2^p,

k-L<m<k

hkn (t) - w(xn, t) < Cx max wm - w(xn, tm)

k - L<m<k

+ C2tp .

Рассмотрим разностную схему

т +1 т т ^ т , т т 7лл,тМлл>т

Уп - Уп _ Уп+1 - 2Уп + Уп-1 ™п+1 - 2^п + ™п-1 + Ф „т + Ф ит (16)

_ ,2 ,2 "'"ча16п + ^1Гп ’ Vі0,)

Т И2 И2

^т - 2^т + ^т ,

о _ < +^№; + п+1 >2п—^+ф 21§Пі +ф пит, (17)

И2

где п _ 1,...,N-1, т _ 0,...,М -1 в (16) и т _ 0,...,М в (17), с начальными условиями

уп _^(Хп ,0), ^ _у(Хп ,0), п _ 0,..., Ы, (18)

&0(0 _Ф(хп,0, И°(0 _^(Хп,і), п _0,...,Ы, іє[-г,0], (19)

и граничными условиями

У0т _у; _< _ w; _ 0, т _ 0,...,М. (20)

Согласно результатам работ [2, 3] необходимым условием разрешимости задачи (1) - (4) является выполнение условия согласования начальных данных ри у

0 _ р(х, 0) + (в + А)у(х, 0) + Ф21р(х, •) + Ф22у(х, •), х є (0, п).

Омельченко Е.А., Плеханова М.В., Численное решение линеаризованной квазистационарной

Давыдов П.Н. системы уравнений фазового поля с запаздыванием

По умолчанию считаем, что оно выполняется.

Невязкой метода (16), (17) назовем сеточную функцию /m = (¿m П), где

¿m = v(Xn , W) - V(Xn , tm ) - V(Xn+1, tm ) - 2v(Xn , tm ) + v(Xn-1, tm ) +

->n 2

t h¿

+ W( xn+1, ,m ) - 2w( V >m ) + W( xn-1. tm ) -фпГ,т (Xn , .)-Ф12 Wm (^,.),

h

П = v(Xn , tm ) + fiw{xn, tm ) + w(Xn +Ltm ) - 2w(^ ^ ) + ^(^tm ) + Ф2/m (Xn, •) + Ф22Wm (Xn, •).

h

Так же, как теорема 4 в [5], с помощью теоремы 1 настоящей работы об устойчивости схемы (11) , (12), теорем 1 и 3 из [5] доказывается следующий результат.

2 2 L+2 г ~\

Теорема 3. Пусть в< 0, т< h , оператор I: R ^ Q [-r,0] Липшицев и имеет порядок

погрешности тР0 на точном решении, невязка / имеет порядок tP1 + hP2. Тогда разностное

решение сходится к точному решению задачи (1)-(4) с порядком Tmin{Р0,Р1} + hP2.

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

4. Численный эксперимент

Пусть Фyg = üjjg(-1) для gе Q[-1,0], i,j = 1, 2. Другими словами, система (1)-(4) имеет вид

vt(X, t) =Av(X, t) - Aw(X, t) + anv(X, t -1) + a12w(X, t -1), (X, t) е [0,^] x[0,T], (21)

0 = v(X, t) + (в +A)w(X, t) + a21v(X, t -1) + a22w(X, t -1), (X, t) е [0,^] x[0,T], (22)

и снабжена краевыми и начальными условиями (3), (4). Разностная схема (16), (17) для нее при выборе, например, кусочно-постоянной или кусочно-линейной интерполяции будет иметь вид

vm+1 vm vm 2vm + vm wm 2wm + wm

vn - vn = vn+1 - 2 n + vn-1 Wn+1 - 2wn + Wn-1 , a vm-L , a w,m-L (23)

----------=--------Г2-----------------------------T"2-+ a11vn + a12 Wn , (23)

t h2 h2

0 = vm I ewm . Wn+1 - 2wn + Wn-1 . a vm-L . a wm-L (24)

0 = vn + PWn +---------—2-+ a21vn + a22 Wn , (24)

h2

где n = 1,..., N -1, m = 0, ...,M -1 в (23) и m = 0,...,M в (24), с начальными условиями (18), (19) и граничными условиями (20). При этом порядок погрешности оператора кусочно-постоянной интерполяции на точном решении равен т, для кусочно-линейной интерполяции - т [9]. По аналогии с леммой 1, нетрудно доказать следующее утверждение.

Лемма 2. Пусть точное решение v, w задачи (21), (22), (3), (4) таково, что функция v дважды непрерывно дифференцируема по t, функции v, w четырежды непрерывно дифференцируемы по

X. Тогда невязка метода (23), (24) имеет порядок т + h2 .

Теорема 3 влечет

Следствие 1. Пусть в < 0, т < h2 , интерполяция кусочно-постоянна или кусочно-линейна.

Тогда разностное решение сходится к точному решению задачи (1)-(4) с порядком т + h2.

На рис. 1 представлено решение (v,w) для параметров в = -0,75, T = 10, M = 1000, N = 16 с начальным условием v(X,0) = sinX, Xе[0,^] для задачи (7)-(10).

Рис. 1

Для задачи с запаздыванием (1)-(4) при значениях ап = 1, а12 = а21 = а22 = 0 , в = -0,75, Т = 10, М = 1000, N = 16, у(х,I) = (I +1)з1пх, (х,I)е[0,^]х[—1,0] решение {у,м>) показано на рис. 2.

Рис. 2

Литература

1. Плотников, П.И. Уравнения фазового поля и градиентные потоки маргинальных функций / П.И. Плотников, A.B. Клепачева // Сиб. мат. журн. - 2001. - Т. 42, № 3. - С. 651-669.

2. Fedorov, V.E. On solvability of some classes of Sobolev type equations with delay / V.E. Fedorov, E.A. Omelchenko // Functional Differential Equations. - 2011. - Vol. 18, № 3-4. -P.187-199.

3. Федоров, B.E. Неоднородные линейные уравнения соболевского типа с запаздыванием / В.Е. Федоров, Е.А. Омельченко // Сиб. мат. журн. - 2012. - Т. 53, № 2. - С. 418-429.

4. Лекомцев, A.B. Полуявный метод для численного решения функционально-дифференциально-алгебраических уравнений / A.B. Лекомцев, В.Г. Пименов // Изв. вузов. Математика. - 2009. - № 5. - С. 62-67.

5. Пименов, В.Г. Разностные схемы в моделировании эволюционных управляемых систем с последействием / В.Г. Пименов // Тр. Ин-та математики и механики УрО РАН. - 2010. - Т. 16, № 5.- C. 151-158.

6. Пименов, В.Г. Разностные схемы численного решения уравнения теплопроводности с последействием / В.Г. Пименов, А.Б. Ложников // Тр. Ин-та математики и механики УрО РАН. -2011. - Т. 17, № 1. - C. 178-189.

7. Федоров, В.Е. Обратная задача для одного класса сингулярных линейных операторнодифференциальных уравнений / В.Е. Федоров, A.B. Уразаева // Тр. Воронежск. зимн. мат. шк. Воронеж: ВГУ. - 2004. - С. 161-172.

8. Рихтмайер, Р. Разностные методы решения краевых задач / Р. Рихтмайер, K.M. Мортон. -M.: Мир, 1972. - 420 с.

9. Ким, A.B. i-Гладкий анализ и численные методы решения функциональнодифференциальных уравнений / A.B. Ким, В.Г. Пименов. - Ижевск: РХД, 2004.- 256 с.

NUMERICAL SOLUTION OF DELAYED LINEARIZED QUASISTATIONARY PHASE-FIELD SYSTEM OF EQUATIONS

E.A. Omelchenko1, M.V. Plekhanova2, P.N. Davydov*

For delayed linearized quasistationary phase-field system of equations the numerical method of solution was proposed. The convergence of explicit difference scheme that takes account of delay in the system under investigation was thoroughly studied. On the basis of the results obtained the implementation of the method was realized.

Keywords: Sobolev type equation, quasistationary phase-field system of equations, difference scheme.

References

1. Plotnikov P.I., Klepacheva A.V. Siberian Mathematical Journal. 2001. Vol. 42, no. 3. pp. 551—

567.

2. Fedorov V.E., Omelchenko E.A. On solvability of some classes of Sobolev type equations with delay. Functional Differential Equations. 2011. Vol. 18, no. 3-4. pp. 187-199.

3. Fedorov V.E., Omel’chenko E.A. Siberian Mathematical Journal. 2012. Vol. 53, no. 2. pp. 335344.

4. Lekomtsev A.V., Pimenov V.G. Russian Mathematics (Izvestiya VUZ. Matematika). 2009. Vol. 53, no. 5. pp. 54-58.

5. Pimenov V.G. Trudy Inst. Mat. i Mekh. UrO RAN. 2010. Vol. 16, no. 5. pp. 151-158. (in Russ.).

6. Pimenov V.G., Lozhnikov A.B. Proceedings of the Steklov Institute of Mathematics (Supplementary issues). 2011. Vol. 275. Suppl. 1. pp. 137-148.

7. Fedorov V.E., Urazaeva A.V. Trudy Voronezhskoy zimney matematicheskoy shkoly (Proc. of the Voronezh Winter Mathematical School). Voronezh: VGU. 2004. pp. 161-172.

8. Richtmyer R.D., Morton K.W. Difference Methods for Initial-Value Problems. Interscience, New York, 1967. 405 p.

9. Kim A.V., Pimenov V.G. i-Gladkiy analiz i chislennye metody resheniya funktsional'no-differentsial'nykh uravneniy (i-differential analysis and numerical methods of solutions of functional and differential equations). Izhevsk: RKhD, 2004. 256 p. (in Russ.).

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

1 Omelchenko Ekaterina Aleksandrovna is Senior Lecturer, Department of Humanities and socio-economic disciplines, Ural Branch of Russian Academy of Justice.

E-mail: omelchenko_ea@mail.ru

2 Plekhanova Marina Vasilyevna is Associate Professor, Differential and Stochastic Equations Department, South Ural State University.

E-mail: mariner79@mail.ru

3 Davydov Pavel Nikolaevich is Post-graduate student, Mathematical Analysis Department, Chelyabinsk State University.

E-mail: davydov@csu.ru