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

удк 532.546+519.6

DOI: 10.14529/mmp 180113

численный метод решения коэффициентнои

обратной задачи для нелинейного уравнения диффузии-реакции

Х.М. Гамзаев, Азербайджанский государственный университет нефти и промышленности, г. Баку, Азербайджан

Рассматриваются две обратные задачи по определению коэффициентов для нелинейного уравнения диффузии-реакции типа Фишера - Колмогорова - Петровского

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

Ключевые слова: уравнение диффузии-реакции; уравнение Фишера - Колмогорова

- Петровского - Пискунова; коэффициентная обратная задача; интегральное условие; явно-неявные схемы.

Введение. Одним из нелинейных уравнений типа «диффузия-реакция> является логистическое уравнение с диффузией

du(x,t) ^d2u(x,t) . . ..„ . .. , ,

dt = D dx2 + ku(x,t)(l - u(x,t)). (1)

Это уравнение предложено H. Колмогоровым, И.Г. Петровским и Н.С. Пискуновым [1] и Фишером [2] для моделирования процесса распространения генной волны. Уравнение (1), также называемое уравнением Фишера - Колмогорова - Петровского -Пискунова (ФКПП), находит применение во многих областях: в задачах тепло- и массообмена, теории горения, биологии и экологии, в физике плазмы и задачах теории фазовых переходов и т.д. [1-3].

Необходимо отметить, что многие характеристики физических, химических, биологических, экологических и т.д. процессов, описываемых уравнением (1), во многом зависит от коэффициентов данного уравнения. В связи с этим важными считаются задачи по определению этих коэффициентов с целью обеспечения желаемого протекания процессов.

1. Постановка задачи и метод решения. Задача А. Пусть рассматривается уравнение ФКПП

du(x't) = D(t)д 2un(x>'t) + k(t)u(x,t)(l - u(x,t)), 0 <x<l, 0 <t < T, (2) dt dx2

где D(t) > 0, со следующими начальным и граничным условиями

u(x, 0) = p(x), (3)

u(0,t) = p(t), u(l,t) = f (t). (4)

Предположим, что помимо функции и(х,Ь) неизвестном является так^ке функция к(Ь). Требуется восстановление этой функции по следующему интегральному условию

I

J и(х,Ь)"х = г(Ь), (5)

0

где г(Ь) - заданная функция. Предполагается, что при этом выполняется условие согласования

I

ф(0) = р(0), ф(1) = !(0), у ф(х)"х = г(0). (6)

0

Таким образом, задача заключается в определении функций и(х,Ь) и к(Ь), удовлетворяющих уравнению (2) и условиям (3) - (5). Данная задача относится к классу коэффициентных обратных задач [4, 51. Отметим, что вопросы существования^ единственности и разрешимости, а также некоторые подходы к численному решению коэффициентных обратных задач для параболических уравнений исследованы в [6-8].

Для решения задачи (2) - (5) сначала дискретизируем уравнение (2) по переменной Ь. Уравнение (2) и условия (2) - (5) запишем при Ь*, ] = 1,ш. Производ-

ди(х, Ь) ди(х, Ь) ■

ную —д^— в уравнении (2) аппроксимируем разностью «назад> '

и (х, Ь*) — и (х, Ь*-\)

Х

используем неявную аппроксимацию по времени, а для нелинейного младшего члена - полуявную аппроксимацию. Обозначив и*(х) ~ и (х,Ь*), задачу (2) - (5) запишем в следующем виде

^ (х) — и- (х) = Б* "^их! + к*и*-1(х)(1 — и*-1(х)), 0 <х<1, (7)

/ \ Ь

и*(0) = р*, и*(I) = ¡*, (8)

I

J и* (х)"х = г*= 1, 2,...,ш, (9)

0

и0 (х) = ф(х), (10)

где к* « к(Ь*),г* = г(Ь*),р* = р(Ь*),!* = f (Ь*),Б* = Б(Ь*).

Предположим, что решение полученной дифференциально-разностной задачи (7) - (10) на каждом временном слое ] = 1, 2,... ,ш можно представить в виде

и* (х) = и* (х) + к*ф> (х), (11)

t=ti ~

~ i_/ У КЛ1Х-Т AJL ^ AAJLAJLA 1 ¿rnd i AJL V./ A W JLiiVXJEJ. У lA'J AJL V./ V _JL X^AV/ lA'J I t_t

dt dt 1 J

где At = T/m шаг по времени. Для диффузионного члена

где wj(х),фj (x), j = l,m неизвестные функции. Подставив соотношение (11) в (7), (8) и учитывая произвольности функций wj(х),ф (х), получим

w (Х) ~ ^ (Х) - В^^^М = 0, 0 <x<l, (12)

At dx2 v '

wj (0)= pj ,wj (l) = fj, (13)

_ D3 f^l _ u3-i(x)(1 _ ui-i(x)) = 0, о <x<l, (14)

At ax2

ф>(0) = 0, ф>(l) = 0,j = 1, 2,... ,m. (15)

Таким образом, решение обратной задачи (7) - (10) uj (x), kj j = 1,m определяется по следующей схеме: для каждого фиксированного значения j = 1, 2,... ,m определяются решения прямых задач (12), (13) и (14), (15); полученные решения подставляются в дополнительное соотношение (9) и из полученного уравнения определяется приближенное значение искомой функции k(t) при t = tj, т.е.

■ (■ ! ■ \ г ■

kj = I rj _ wj(x)dx I / фj (x)dx

V о ) {

по формуле (11) определяется приближение к искомой функции u(x,t) при t = tj.

Для численного решения задач (12), (13) и (14), (15) можно использовать метод конечных разностей [6].

Задача В. Предположим, что снова рассматривается уравнение ФКПП

9u(x,i) = D(t) d2u(xo,tl + k(t)u(x, t)(1 _ u(x, t)), 0 <x<l, 0 <t < T, (16) dt dx2

со следующими начальным и граничным условиями

u(x, 0) = ф), (17)

u(0,t)= p(t),u(l,t) = f (t). (18)

u(x, t) неизвестном является также D(t)

восстановление этой функции по следующему интегральному условию

i

Iu(x-t)dx =r(t)■ (19)

о

Предполагается, что выполняются условия согласования (6).

Предположим, что коэффициент диффузии представляется в виде D(t) = Do + d(t) > 0, где D0 = const > 0 заданное число, a d(t) неизвестная функция.

du(x, t)

Снова дискретизируем производную —^— разностью «назад> в уравнении (16)

при tj, j = 1, m и используем явно-неявную аппроксимацию по времени для диффузионного члена, а полуявную аппроксимацию для нелинейного младшего члена. Тогда задача (16) - (19) с учетом представления коэффициента диффузии запишется в следующем виде

u (x) _ ■ (x) = Do dj) + dj duj-1(x^ + kj uj-1(x)(1 _ uj—l (x)), 0 <x <l, (20) At dx2 dx2

uj(0) = pj, uj(l) = fj, (21)

I

! и (х)$х = г3^ =1, 2,... ,т, (22)

о

и0 (х) = ф(х), (23)

где й3 ~ ¿(Ь^).

Предположим, что решение дифференциально-разностной задачи (20)^(23) на каждом временном слое ] = 1, 2,... ,т, можно представить в виде

и3 (х)= т3 (х) + $ф3 (х) , (24)

где т3 (х), ф3 (х), ] = 1,т- неизвестные функции. Подставив соотношение (24) в (20), (21) получим следующие краевые задачи относительно функций V3 (х), ф3 (х)

^ (х) -и3-1 (х) - Бо - к3 и3-1(х)(1 - и3-1 (х)) = 0, 0 <х<1, (25)

I \ ь

V3 (0)= р3, V3 (]) = ¡3. (26)

ф3 (х) _ й2ф3 (х) ¿2и3-1(х) < < ]

"дТ - Во^ъ2---1ыг~ = 0 0 <х<], (27)

ф3 (0) = 0, ф3 (]) = 0,з = 1, 2,..., т. (28)

Подставив (24) в интегральное условие (22), получим

I \ I

$ = (И - J т3 (х)$х I / ! ф3 (х)$х. (29)

оо

Таким образом, решение обратной задачи (20) - (23) uj (x), dj, j = 1,т определяется по следующей схеме: для каждого фиксированного значения j = 1, 2,... ,т, определяются решения прямых задач (25), (26) и (27), (28); по формуле (29) определяется приближенное значение искомой функции d(t) при t = tj\ и, по формуле (24) определяется приближение к искомой функции u(x,t) при t = tj. Для численного решения задач (25), (26) и (27), (28) также можно использовать метод конечных разностей.

2. Результаты численных расчетов. На основе предложенных вычислительных алгоритмов были проведены численные эксперименты для модельных задач. Ниже приводятся результаты численного эксперимента для двух модельных задач. Схема численного эксперимента заключается в следующем: задаются функции k(t) (в задаче А) и D(t) (в задаче В) и определяются решения прямых задач (2) - (5) и (16)

i

- (19). Далее по формуле r(t) = f u(x,t)dx определяется функция r(t), и наиденная

0

зависимость принимается за точные данные для восстановления функции k(t) и D(t), соответственно.

Результаты численного эксперимента, проведенного для случая l = 1, T = 1, k(t) = 3 + 2cos1000t, Do = 0, 5, d(t) = 0, 3 + 0, 2sin100t,p(t) = 1, f (t) = 0, <p(x) = 1 - x представлены в таблицах 1, 2.

Таблица 1

Численные результаты по задаче А_

tj Значение функции k(t)

Точное Вычисленное при At =10 4, Ax = 4 • 10 2

0,1 4,725 4,725

0,2 3,974 3,974

0,3 2,956 2,956

0,4 1,949 1,949

0,5 1,232 1,232

0,6 1,002 1,002

0,7 1,322 1,322

0,8 2,104 2,104

0,9 3,132 3,132

1,0 4,125 4,125

Таблице Численные результаты по задаче В

tj Значение функции d(t)

Точное Вычисленное при At =10 4, Ax = 4 • 10 2

0,1 0,191 0,192

0,2 0,483 0,483

0,3 0,102 0,102

0,4 0,449 0,449

0,5 0,248 0,247

0,6 0,239 0,239

0,7 0,445 0,445

0,8 0,101 0,101

0,9 0,479 0,479

1,0 0,199 0,199

Результаты численных экспериментов показывают, что искомые функции k(t) и D(t) восстанавливаются с высокой точностью при всех расчетных сетках по времени.

Анализ результатов численного экспериментов свидетельствует, что для повышения точности решений необходимо использовать мелкие шаги разностной сетки.

Заключение. В предложенном вычислительном алгоритме эффект регуляризации обеспечивается за счет выбора разностной сетки по времени.

Литература

1. Fisher, R.A. The Wave of Advance of Advantageous Genes / R.A. Fisher // Annals of Eugenics. - 1937. - № 7. - P. 355-369.

2. Колмогоров, A.H. Исследование уравнения диффузии, соединенной с возрастанием количества вещества, и его применение к одной биологической проблеме / А.Н. Колмогоров, ИТ. Петровский, И.С. Пискунов // Бюллетень МГУ. Секция А. - 1937. - Т. 1, № 6. - Р. 1-25.

3. Маслов, В.П. Исследование уравнения диффузии, соединенной с возрастанием количества вещества, и его применение к одной биологической проблеме / В.П. Маслов, В.Г. Данилов, К.А. Волосов. - М.: Наука, 1987.

4. Самарский, A.A. Численные методы решения обратных задач математической физики / A.A. Самарский, П.Н. Вабищевич. - М.: Изд-во ЛКИ, 2009.

5. Кабанихин, С.И. Обратные и некорректные задачи / С.И. Кабанихин. - Новосибирск: Сиб. науч. изд-во, 2009.

6. Камынин, В.Л. Обратная задача определения младшего коэффициента в параболическом уравнении при условии интегрального наблюдения / В.Л. Камынин // Математические заметки. - 2013. - Т. 94, № 2. - Р. 207-217.

7. Кожанов, А.И. Параболические уравнения с неизвестными коэффициентами, зависящими от времени / А.И. Кожанов // Журнал вычислительной математики и математической физики. - 2017. - Т. 57, № 6. - С. 961-972.

8. Kerimov, N.B. An Inverse Coefficient Problem for the Heat Equation in the Case of Nonlocal Boundary Conditions / N.B. Kerimov, M.I. Ismailov // Journal of Mathematical Analysis and Applications. - 2012. - V. 396, № 2. - P. 546-554.

Хан. lap Мехвалн оглу Гамзаев, доктор технических наук, профессор, кафедра «Общая и прикладная математика:», Азербайджанский государственный университет нефти и промышленности (г. Баку, Азербайджан), xan.h@rambler.ru.

Поступила в редакцию 12 сентября 2017 г.

MSC 65N21, 80A23 DOI: 10.14529/mmpl80113

A NUMERICAL METHOD OF SOLVING THE COEFFICIENT INVERSE PROBLEM FOR THE NONLINEAR EQUATION OF DIFFUSION-REACTION

Kh.M. Gamzaev, Azerbaijan State Oil and Industry University, Baku, Azerbaijan, xan.h@rambler.ru

We consider two inverse problems for determining the coefficients for a one-dimensional nonlinear diffusion-reaction equation of the Fisher-Kolmogorov-Petrovsky-Piskunov type. The first problem consists in determining the kinetic coefficient for a nonlinear lower term, depending only on the time variable, according to a given integral condition. And the second problem consists in determining the time-dependent diffusion coefficient, again according to a given integral condition.

To solve both problems, the time derivative of the derivative is first sampled. In the first problem, the diffusion term is approximated in time according to the implicit scheme, and the nonlinear minor term in the semi-explicit scheme. And in the second problem, the diffusion term is approximated in time in an explicitly implicit scheme, and the nonlinear minor term is again in a semi-explicit scheme. As a result, both problems reduce to differential-difference problems with respect to functions depending on the spatial variable. For numerical solution of the problems obtained, a non-iterative computational algorithm is proposed, based on reduction of the differential-difference problem to two direct boundary-value problems and a linear equation with respect to the unknown coefficient. On the basis of the proposed numerical method, numerical experiments were performed for model problems.

Keywords: diffusion-reaction equation; Fisher-Kolmogorov-Petrovsky-Piskunov equation; coefficient inverse problem; integral condition; differential-difference problem; explicitly implicit schemes.

References

1. Fisher R.A. The Wave of Advance of Advantageous Genes. Annals of Eugenics, 1937, no 7, pp. 355-369. DOI: 10.1111/j.l469-1809.1937.tb02153.x

2. Kolmogorov A.N., Petrovsky I.G., Piskunov I.S. [A Study of the Diffusion Equation with Increase in the Amount of Substance, and Its Application to a Biological Problem]. Bulletin of the Moscow State University, Section A, 1937, vol. 1, no. 6, pp. 1-25.

3. Danilov V.G., Maslov V.P., Volosov K.A. Mathematical Modelling of Heat and Mass Transfer Processes, Dordrecht, Kluwer, 1995. DOI: 10.1007/978-94-011-0409-8

4. Samarskii A.A., Vabishchevich P.N. Numerical Methods for Solving Inverse Problems of Mathematical Physics. Walter de Gruyter, 2007. DOI: https:10.1515/9783110205794

5. Kabanikhin S.I. Inverse and Ill-Posed Problems. Theory and Applications. De Gruyter, Germany, 2011. DOI: 10.1515/9783110224016

6. Kamynin V.L. The Inverse Problem of Determining the Lower-Order Coefficient in Parabolic Equations with Integral Observation. Mathematical Notes, 2013, vol. 94, no. 2, pp. 205-213. DOI: 10.1134/S0001434613070201

7. Kozhanov A.I. Parabolic Equations with Unknown Time-Dependent Coefficients. Computational Mathematics and Mathematical Physics, 2017, vol. 57, no. 6, pp. 956-966. DOI: 10.1134/S0965542517060082

8. Kerimov N.B., Ismailov M.I. An Inverse Coefficient Problem for the Heat Equation in the Case of Nonlocal Boundary Conditions. Journal of Mathematical Analysis and Applications, 2012, vol. 396, no. 2, pp. 546-554. DOI: 10.1016/j.jmaa.2012.06.046

Received September 12, 2017