CC BY
150
= Вычислительная математика =
УДК 517.96, 517.956, 519.63 DOI: 10.14529/cmse160204
ЧИСЛЕННЫЙ МЕТОД РЕШЕНИЯ ОБРАТНОЙ ЗАДАЧИ С НЕИЗВЕСТНЫМИ НАЧАЛЬНЫМИ УСЛОВИЯМИ ДЛЯ НЕЛИНЕЙНОГО ПАРАБОЛИЧЕСКОГО УРАВНЕНИЯ
Н.М. Япарова
В статье рассмотрена обратная задача для нелинейного параболического уравнения с неизвестными начальными условиями. Для решения обратной задачи предложен метод дискретной регуляризации, основанный на использовании конечно-разностных уравнений и применении регуляризирующих функционалов. Построенная вычислительная схема позволяет одновременно найти численное решение внутри рассматриваемой области и неизвестную граничную функцию. В статье проведено исследование устойчивости вычислительной схемы. Выявлено влияние величин шагов дискретизации и погрешности исходных данных на устойчивость численных решений. Предложенная схема послужила основой для разработки численного метода и проведения вычислительного эксперимента. Результаты эксперимента для серии тестовых функций также представлены в данной работе и свидетельствуют о достаточной эффективности предложенного метода дискретной регуляризации.
Ключевые слова: обратные задачи, численный метод, метод регуляризации, оценка погрешности, вычислительная схема
ОБРАЗЕЦ ЦИТИРОВАНИЯ
Япарова Н.М. Численный метод решения обратной задачи с неизвестными начальными условиями для нелинейного параболического уравнения // Вестник ЮУрГУ. Серия: Вычислительная математика и информатика. 2016. Т. 5, № 2. С. 43-58. DOI: 10.14529/^е160204.
Введение
Исследования многих процессов тепло- и массопереноса приводят к формированию математических моделей, представленных нелинейными параболическими уравнениями с граничными условиями, известными на части границы. В этих задачах требуется найти неизвестную функцию, удовлетворяющую имеющимся уравнениям внутри рассматриваемой области, а также граничным условиям. Примеры подобного рода задач можно найти в работах [1-6]. Один из подходов, позволяющий найти требуемые функции, заключается в решении обратных граничных задач. При этом сначала на оставшейся части границы восстанавливают граничную функцию, а уже затем, используя полученные результаты, находят решение поставленной задачи во внутренних точках. Исследованиям в этой области посвящены, например, работы [7-13]. Другое динамично развивающееся направление решения задач тепло- и массопереноса связано с разработкой и исследованием численных методов, основанных на использовании конечно-разностных и дифференциально-разностных схем. Различные методы, относящиеся к этому направлению, представлены, например, в работах [15-19]. Общей чертой существующих подходов является то, все эти методы разрабатываются с учетом известных начальных условий. С другой стороны, протекание многих процессов тепло- и массооб-мена связано с ситуацией, когда осуществить измерение во всем теле требуемой температурной или концентрационной функции не представляется возможным. Это приводит
к возникновению математических моделей с неизвестными начальными условиями. Один из путей решения данной проблемы заключается, например, во введении фиктивных начальных условий [1-2, 14], что неизбежно отражается на точности получаемых решений.
В данном исследовании полагают, что начальные условия неизвестны исследователям и для построения численного решения обратной граничной задачи с неизвестными начальными условиями предлагают метод дискретной регуляризации, основанный на использовании конечно-разностных схем. Основная идея метода заключается в том, что, в отличие от существующих подходов, где вычисление требуемой функции во внутренних точках рассматриваемой области осуществляется по переменной, характеризующей время, мы находим неизвестную функцию относительно пространственной переменной. Такой подход позволяет значительно уменьшить влияние неизвестных начальных условий на результаты численного решения задач с неизвестными начальными условиями. Подобный подход был рассмотрен в работах [20-21] для линейных параболических уравнений.
Статья имеет следующую структуру. В первом разделе приведена постановка задачи, а также дополнительные свойства и ограничения на рассматриваемый класс функций. Во втором разделе предложена вычислительная схема метода дискретной регуляризации. Третий и четвертый разделы посвящены исследованию устойчивости основной вычислительной схемы метода. В третьем разделе установлена зависимость погрешности общей вычислительной схемы от величин шагов дискретизации. В четвертом разделе исследована локальная устойчивость метода при соответствующем выборе соотношения шагов дискретизации и параметра регуляризации, и получены некоторые оценки погрешности численного решения рассматриваемой задачи. В пятой части статьи приведены результаты вычислительного эксперимента, основными целями которого являлись проверка принципиальной возможности построения численного решения задачи, устойчивого в некоторой подобласти и получение экспериментальных оценок погрешности метода. В заключении сформулированы итоговые выводы, полученные в результате проведенных исследований, а также указаны основные направления дальнейших исследований.
1. Постановка задачи
Рассмотрим цилиндрическую область R3 и 0 = Sí х R, где Sí — кольцо, ограниченное окружностями д1Sl и д2Sl с радиусами 10 и I соответственно. Множество О ограничено поверхностями д10 = д1£1 х R и д20 = д2Se х R , а Q = Se х (0, Т) . Рассмотрим уравнение
с(и)р(и)иг = div (А(и)Уи), (х, у, z, t )еПТ, (1)
где t е (0,Т) , z е R , От =0х(0,Т) , а функция и(х,у,z,t) - осесимметрична и принимает постоянные значения вдоль любой прямой, параллельной оси Oz. Дополнительно известны граничные условия:
и кох[0,Т ]= Р (t)' | кох[0,Т ]= Ц (t)' (2)
где вектор п является внешней нормалью к д10, при 10 ^ 0. Коэффициентные функции с(и), р(и),Л(и) в (1) положительны, дважды непрерывно дифференцируемы в со-
ответствующей области. Известно, что при некоторых р(1:) = р0(:) и g (t) = gQ(t) существует функция и0, удовлетворяющая (1), (2) но вместо точных значений р0(:) и qQ(t') известны некоторые приближения р8 , qs и некоторый уровень погрешности 8 такие, что ||ро(0-р81с([0г])-8 и -qs||C([QГ])-8. в данной задаче необходимо найти
и(х, у, z, t) , удовлетворяющую (1), (2) и граничное значение функции
и 1э2ПХ(о,г ) = ф) (3)
Из результатов, представленных в работе [16, 23] следует существование функции иОа е С2Д (DT )|П С (Dт ) , удовлетворяющей в некоторой области DT ^О.Т соотношениям
(1), (2), а также ф8 е Н2'1+в(0, Т) , где в е (0,1) , удовлетворяющей (3).
Из результатов измерений величин с(и), р(и),Л(и) в используемом материале известно, что зависимость этих коэффициентов от значений функции и( х, у, z, t) допускает представление в виде полиномов. На практике наиболее часто используют линейную зависимость. При этом точность такого представления такова, что погрешность полиномиальных аппроксимаций не оказывает значительного влияния на величину погрешности определения требуемой функции и( х, у, z, t) в исследуемом теле. Учитывая свойства с(и), р(и),Л(и) , для вычисления этих величин используем разложение по формуле Тейлора в окрестности точки и* = и(0, 0,0,0) :
с(и) = с(и*) + с'(и*)(и(х, у, г, {) - и*), р(и) = р(и*) + р'(и*)(и(х, у, г, {) - и*),
Л(и) = Л(и*) + Л'(и*)(и(х,у,г,{) - и*),
где ^х2 + у2 - 10, а значения с(и*),р(и*),Л(и*) , с'(и*), р'(и*), Л'(и*) получены из результатов измерений соответствующих величин в материале, используемом в исследуемом теле. Тогда полагаем
kc = с'(и*) , kp = р'(и*) , kЛ = Л'(и*) , с0 = с(и*) - kcu*, р0 = р(и*) - kcu*, Л0 = Л(и*) - kcu*.
Получаем, что с(и), р(и),Л(и) определяются формулами
с(и) = с0 + kcu, р(и) = р0 + kрu, Л(и) = Л0 + kЛu .
Перейдем к цилиндрической системе координат, выполняя сдвиг по радиусу. Тогда область Q преобразуется в Ql = {(г, t ):0 < г < R,Q < t < Т} , где R = I -10, а уравнение
(1) примет вид:
ди 1 д ( ди 1
с(и)р(и)— =--1 Л(и)г— I, (г, t)е Qf.
дt г дг { дг ) У } 1
д
Раскрывая оператор —, получаем следующее представление задачи (1), (2):
дг
с(и)р(и) ^ =дЛии1 +1 Л(и) д-и + Л(и) ^, (г, t) е Qt. (4)
дt ди \дг) г дг дг2 к ' " '
и (0, ^ = р (t), иг (0, t) = q , t е[0, Т], (5)
где коэффициенты с(и), р(и), Л(и) определим формулами
с(и) = с0 + kc •и, р(и) = р0 + k •и, Л(и) = Л0 + kЛ и. (6)
Требуется по ps , qs и 8 данным найти функцию uas(r, t) , являющуюся решением задачи (4)—(6) и определить граничную функцию (fr^it) . Если из теплофизических характеристик процесса известно, что функция u(r, t) достаточно гладкая, и мы можем полагать, что существуют константы Ф, в, R, C такие, что max |u(r, t)| < Фев(r+t) в каждой точке (r, t) e Qf , а
max { max |u.. I, max \urr I, max \urrrr |} < C , (7)
l(r,t)egj ttl (r,t)egj rrl (r,t)egj 1) то для решения задачи (4)—(6) используем метод дискретной регуляризации. Вычислительная схема метода представлена в следующем разделе.
2. Метод дискретной регуляризации
Рассмотрим конечно-разностную сетку G в области Qt,
(r, ): r = (i - 1)h, = (j — 1)т,
G =
h = R / N1 т = Т / М; I = 1, N +1; j = 1,М +1
где значение М не превосходит количества измерений, выполненных для р(Х) и q(t). Обозначим щ - = и(г, tj) . Составим конечно-разностный аналог задачи (4)-(6). С
этой целью используем конечно-разностное представление частных производных по г и по t в каждой точке G , предложенное в [22] и рассмотрим множество дискретных функций Ук = {у(г,^-) = V..} , заданных на G . Конечно-разностный аналог уравнения (4) имеет вид:
V. , . — v..
1+1,1 1,1
h2 (i — 1) h h 1,1 h2
В отличие от существующих методов, использующих данный шаблон, мы выразим из этого уравнения +1 ■, что позволит значительно уменьшит влияние неизвестных
начальных условий на результаты численного решения. Получаем при г = 2, N и
- = 1, М
V. 1 . = V. . + --У- + 4 (г -1) (у - у, 1 _) (8)
Граничные условия в конечно-разностной форме представлены следующим образом:
= ии = Р- У2, - = + ^ = Р- + ^ (9)
Непосредственное использование уравнения (8) не позволяет получать устойчивые решения задачи (4)—(6), более того, использование уравнения (8) невозможно при
- = М +1. Поэтому для построения численного решения рассматриваемой задачи используем метод дискретной регуляризации. Суть метода заключается в том, что коэффициенты в уравнении (8) вычисляются по формулам
= С0 + kcVгj , Рг, - = Ро + kp\- , 4 - = 40 + КУ, - , (10)
4+и- = 4 + (а)У,- = -,
т
где а — некоторый параметр регуляризации и шаги дискретизации выбираются согласовано с уровнем погрешности, а для вычисления значений viM +1 используется следующая схема:
(i - Л (с п I
с,,М+1Иi М+1
V
Vi+1,М+1 Vi,M+1 +
К+1,М+1 (* - 1) + Лг
,,М+1
Т(1 - а) ^ м+1 - VМ) + Л, М+1 (V М+1 - у,-1м+1) I +
+ х лм (* - 0 + лм)(^+1М - VМ) - Лм 0'- ])(\м - У-1М) (11)
(1 -а) (I -1) И2* Л+1М+1 (I - \) + Л1М+1 ' 1 ]
где ае(0,1), а коэффициенты в уравнении (11) также вычисляются согласно формулам (10) и шаги дискретизации также согласуются.
3. Алгоритмические особенности метода дискретной регуляризации
Исследуем устойчивость вычислительной схемы (8)-(10). С этой целью, мы используем подход, принятый в теории регуляризации, см., например, [24]. Введем в рассмотрение функции wi при всех / = 1,N , определяемые формулой:
4+1 = ПМ, Ь+1,1 - ,
1 =1,М+1 1 ^ ^ 1
где ■ удовлетворяет (8)-(10) и определим условия, при соблюдении которых будет иметь место неравенство wN+1 — К8 , где К — некоторая константа.
Так как значение М не превосходит количества измерений, выполненных исследователями для функций р(Х) и q(t), то мы полагаем, что М и, следовательно, величина шага дискретизации т по переменной, характеризующей время, нам известна. Выберем значения N , характеризующей шаг дискретизации по пространственной переменной и параметр регуляризации а так, чтобы при известном значении М имело место требуемое неравенство.
К у 0 -1)
Сначала оценим а и коэффициенты в (8). Обозначим ц = тах
1 К? -а(1 -1) 1 V,у Г
П = тах
i, 1
с- П у ^ -1) И2
(К+11 О - О + К)т
Полагаем а таким, что при всех i = 1, N +1 и
1 = 1, М +1 выполнено
Имеем:
Л
а < тт
1-
1 (! - 1)1Ч.у \
К (! - 1)
N +1 ' Ли1 -а(I -1) V,
<1,
(12)
(13)
тогда ц < 1. С учетом (10) получаем следующую оценку:
4+1 —птах1 +1 - 1+М4.
1 =1М
Так как нам не известны начальные условия, то мы можем учитывать только ограниченность функции и(г, ^) . Имеем:
mM vi,j+1 - Ч j <
j =i m j j
Полагая S = max \q,\, учитывая (9) и выбирая N > R +1, получаем
1=1,M+11 1
w2 = max | v2 i - v1 1 |< max | qh |< h(S + 5) < S + 5,
1=1 M+1 1=1 M+1 1
wi+1 < 20eR+T (n + ßn + ß2n +... + П-2) + (S + 5)ßi-1 < 20eR+T+ (S + 5)ßi-1, v ' 1 - ß
i+1
ß
при i = 3, N +1. Тогда
Wn+1 <n20eR+T--- + (S + 5)ßN-2. (14)
1 -ß
Из (10) следует, что имеет место оценка n< — max__ыРй—. Отсюда и того,
т i=2,N+1,1=1M avij (i -1)
что a(i - 1)vi 1 <a(i -1)V 1 < 41 следует, что для получения устойчивого решения
h2 a(i -1) v , . 4,
необходимо выполнение условия —< max__— < min__:— . Учитывая
т i=2,N+1,1=1M С Р i=2,N+1,1=1 M c. P .
i, j i, j i, j i, j
эти особенности, получаем, что шаги дискретизации должны удовлетворять соотношению
h2 1 - ( r+t )„ . Л
— < ^Фе(R+T5 2 min м 1 —^-. (15)
Тогда
т 2 i=2,N+1,1=1 M+1 c p 1
с
^^+1 < т-5 + (^ + 5)Ц-2 < К5. 1 -ц
Принимая во внимание, что при предложенном подходе ц < 1, получаем, что для выполнения последнего неравенства необходимо, чтобы величина N, характеризующая количество частей разбиения отрезка [0, Л], удовлетворяла соотношению:
(( 1 ^ 5 ^ N > 2 + ^ К--- — , (16)
И 1 8+5)
( 1 ^ 5 . „ _ 1
K-
< 1. Последнее достигается, если положить K = 1 + -
где ____________ _______ _
^ 1 - ц )8 + 5 1 -ц
Таким образом, соотношение wN+1 < К5 имеет место при известном значении М, если выполнено условие (12), (15) и (16).
4. Локальная устойчивость метода дискретной регуляризации
Выполним оценку погрешностей конечномерных аппроксимаций при известном М, выбранном из условия (16) значении N и а, удовлетворяющем полагая, что коэффициентные погрешности не оказывают существенного влияния на суммарную погрешность метода.
Используя конечномерные аппроксимации частных производных, получаем, что при
1 = 2, N +1 и j = 1 ,М имеет место следующее представление:
*+1 > = и > + (и- У -,) +
г+1,1 ^ 1 Л+1, у (- - О + Л, т к+,у (- - + 1 1-1, 1
+ с- Р • (- -1)И п(т)+2(к+и)(- -1)И + -Т- л—Г~п(т) + 2~1-Т- л—\-П(И) +
Л+1И (1 - Ч + ку Л+1И (1 - Ч + ку
1 \ / ',]
(к, -к,) (1 -1) и2 ч 2 ЛИ к, (1 -1) И2 2
(и+11- )(к+1,1(1 -1)+К, У ' к+1,1(1 -1)+Л к+1,1(1 -1)+Л
Введем в рассмотрение функции si, определяемые формулой:
£ = тах V., - и А, 1 = 1,N +1.
' 1=1 М+1 '' 1 ''1
где 1 удовлетворяет (8)—(10), а ui л является решением задачи (4)-(6). Для исследования устойчивости необходимо оценить величину sN+1 при выбранных значениях М, N и а. Имеем:
£ < £ + Счрч( -1) И 2 2Фе(*+т) + 2 Л 1(1 -1) £ + 1 (к+, 1(1 - 1) + к)т к+1, 1(1 - 1) + к1 -
С, Р 1 (1 -1)И2 (Л,+, 1 - к Л (1 -1)И (Л+ г - к Л (1 -1)И2 Ч2
+ к-1 (д 1) +л п(т)+2К л1,1 г \ Л п(и) + ' ^ \(л])V ;) , ,(п(И))2 +
Л+и (1-1) + к, 1 к+ц (1-1) + к, 1 (и1+1,1 - ии 1 )(ЛМ1 (1-1) + к, 1-У
к ии к Л1 -1) И 2 2
П(И) + Л '-1}. П(И2) (17)
к+1,1(1 - 1) + Л„ к+1.1(1 - ч + к, 1
Оценим слагаемые в (17). С учетом (14) и (15) имеем:
С
Р 1(1 -1) И
1Н1,1 V_Л!_пъЛ К+т)
2Фе(Л+т) (18)
к+1, А1 - 1) + к1т
Из соотношения (11), определения /и и того, что а1 V ^ | (1 - 1)<к1-, следует
■(1 -1) И К 1 (19)
к+1.1 (1 -1)+Л к+1, 1 (1 -1)+Л к+и (1 -1)+к. 1
Л- ^л .
Из (6) следует, что ' + ^-г— = кЛ. С другой стороны, из (10) следует, что
и1+1, - и1,
Л1+1 , = к 1 - аУ1 ,, тогда, полагая ЛЛ = _min_Л1 : и учитывая (12) и (19), имеем:
22
к+1,1 - к1 (1 - 1) И2
и1+1,1 иг,1 к+11 (1 - 1) + к,1
И1 к,И1
"Л '
к,1
< кл^-<^ , (20)
1 - 1) <„<! Л1 < /
</< 1, -т71-< — </< 1. (21)
к+11(1 - 1) + к1 ' л+и(1 - 1) + Ли 1 -1
Из соотношений (7), (17)-(21) следует, что
£1+1
<(1 + 2/) £г +5 + {пСг2 + 3/СИ2 +Л С2 +/С И4 , 1 = 2, N. (22)
V V1 у У
2
Полагаем *¥ = цСт2 + ЪцСЪ1 + | C2 +^C| h4 . Из того, что s2 = max |v2 , - u2 J <
I Л ) j=l-M+i1 1
u1 . + hq, - (u0,. + hqj) + Ch <S + hS + Ch , где U0., q0 — точные значения соот-
< max
j=l M+1
ветствующих функций, а u1 .„q. — приближенные значения, следует, что соотношение (22) имеет вид
S+j < (1 + 2М) st + S + ^ < (1 + 2^)г 1 S2 + (1 + 2М)к j (S + .
Тогда при i = 2, Ж выполнено
<(1 + 2^)N-1 S + hS + Ch + (S + Y)|. (23)
Г
SN+1
1
^ 4 ^^
Из соотношений (15), (16), (23) следует, что устойчивость метода зависит от величины шагов дискретизации и для получения приближенного решения необходимо выбирать достаточно небольшие значения N .
5. Вычислительный эксперимент
Основная цель вычислительного эксперимента заключалась в проверке принципиальной возможности построения численного решения во всей области рассматриваемой обратной нелинейной задачи методом дискретной регуляризации и получение соответствующей граничной функции (4). Кроме того, для получения экспериментальных оценок погрешностей были вычислены величины отклонения найденных граничных функций от тестовых значений. Основные этапы эксперимента заключаются в следующем. На начальном этапе эксперимента находим решение прямой задачи:
с(и)р(и) ди = ^(М +1 Л(и) ди + Л(и) % (г,/)е б,. (24)
дг ди ^Г J Г дг дг2 К ' " '
и(0,г) = р(г), и(Я, г) = ф(г),г е[0, Г], (25)
и(г,0) = / (г), г е[0, R] (26)
где коэффициенты определены формулами:
с(и) = с0 + kc •и, р(и) = р0 + k и, Л(и) = Л0 + ^-и, (27)
и 2 1 и •
Далее, определяем функцию q(t) по формуле д(г:) = —:-— , а затем находим
h
р3(г}.), qs(tj) как значения случайных величин, равномерно распределенных на соответствующих интервалах р ()-£, р (^ ) + ^ и q (q (^ ) + £ .
На следующем этапе с помощью предложенной вычислительной схемы получаем численное решение задачи (4)—(6) для р5, qs ,как внутри области, так и на границе.
Результаты численного решения задачи (4)—(6) для некоторых тестовых функций представлены на нижеприведенных рисунках и табл.1. Так как исходное уравнение (1) не содержит дополнительной функции источника, то в качестве тестовых функций использовались функции, представимые в экспоненциальном виде. Все рисунки имеют одинаковую структуру и используют одинаковые обозначения. На левых рисунках изображены графики граничной функции. Обозначение и0 соответствует графику тестовой
функции, а uas — численному решению задачи идентификации граничной функции.
Многомерные поверхности, представленные на правых рисунках, иллюстрируют решения прямой задачи (24)—(27) и обратной задачи (4)—(6), полученные во всей области Qt.
Поверхности, имеющие название "Exact solution" соответствуют численному решению прямой задачи (24)-(27), а поверхности "Regularized solution" — численному решению задачи (4)-(6).
Рис. 1. Результаты численного решения нелинейной обратной граничной задачи для тестовых функций и (0, г) = и (R, г) = ге г
Рис. 2. Результаты численного решения нелинейной обратной граничной задачи для тестовых функций и ( 0, г) = г(е~* - е-1), и (R, г) =
Рис. 3. Результаты численного решения нелинейной обратной граничной задачи для тестовых функций и (0, t) = te~t, и (Я, t) = 1 - ^ -1)2
Рис. 4. Результаты численного решения нелинейной обратной граничной задачи
1 - е~Ч7 -
для тестовых функций и ( 0, t) =-—^ -t, и (Я, t) = te г
1 - е
Рис. 5. Результаты численного решения нелинейной обратной граничной задачи
- 1 - е~*1г
для тестовых функций и (0, t) = te г, и (Я, t) =-—— -1 , где у> 0
1 - е
С целью получения экспериментальных оценок погрешностей были найдены величи-
ны А5( N) =
и0 и8
С(0,Т)
, где и0 — тестовое значение функции и (R, г) , а и£ — чис-
ленное решение задачи (4)—(6). Для краткости изложения в табл. 1 приведены некоторые средние значения АД N) , полученные при 5 = 0,05 и соответствующей реализации
условий (5) в тестовых функциях.
Таблица 1
Экспериментальные оценки погрешности численных решений
Тестовые функции N А5 (N)
и (0,г ) = геи (R, г ) = ге- 5 0,0377
6 0,0346
7 0,0282
8 0,0252
9 0,5 -1081
10 5 •10208
/п ч 1 - е~ф и (0, г )= ,, г, V ; 1 - е~1г у = 0.05 и (R, г ) = гв-, 5 0,1084
6 0,1026
7 0,0805
8 0,0765
9 0,5 -105
10 2 • 10133
и (0, г ) = ге-, 1 в-1 г 7 = °.°5 и (R'г )= 1 е-У г- ^ 5 0,2795
6 0,2688
7 0,2556
8 0,2555
9 2 -1056
10 0,5 •Ю123
Результаты эксперимента свидетельствуют о том, что предложенный метод дискретной регуляризации позволяет с достаточной точностью определять неизвестную граничную функцию, а также находить требуемую функцию и(г, г) во всей рассматриваемой области при условии согласования величин шагов дискретизации и параметра а.
Заключние
В статье рассмотрена обратная граничная задача для нелинейного параболического уравнения с неизвестными начальными условиями. Для численного решения задачи предложен метод дискретной регуляризации и проведено исследование устойчивости основной вычислительной схемы метода. С целью проверки эффективности предложенного метода и получения экспериментальных оценок погрешностей был проведен вычислительный эксперимент. Результаты эксперимента согласуются с полученными теоретическими оценками и свидетельствуют о том, что метод обладает свойством саморегуляризации. Предложенная вычислительная схема позволяет значительно уменьшить влия-
ние неизвестных начальных условий и, соответственно, получать численные решения в соответствующей подобласти с достаточной точностью. Результаты данного исследования могут быть использованы при решении прикладных задач, возникающих, например, в металлургии, в метрологии, а также при исследовании наноструктур.
Для уточнения оценок погрешностей предложенного метода в дальнейшем планируется провести исследование устойчивости полной вычислительной схемы и оценить влияния саморегуляризующих свойств алгоритма на величину погрешностей численных решений.
Литература
1. Алифанов О.М. Обратные задачи теплообмена. М.: Машиностроение, 1988. 280 с.
2. Самарский А.А., Вабищевич П.Н. Вычислительная теплопередача. М.: Едиториал УРСС, 2009. - 784 с.
3. Прокудина Л.А., Вяткин Г.П. Самоорганизация возмущений в жидких пленках // Доклады Академии наук. 2011. Т. 439, № 4. С. 481-484.
4. Булгакова Г.Т., Кондратьева Н.Р. Аналитическая модель вертикального вытеснения нефти водой с учетом вязкостных, гравитационных и капиллярных // Вестник Самарского государственного технического университета. Серия: Физико-математические науки. 2012. № 1. С. 208-213.
5. Шестаков А.Л. Методы теории автоматического управления в динамических измерениях. Министерство образования и науки Российской Федерации, Челябинск: Изд-во ЮУрГУ, 2013. 257 с.
6. Мартинсон Л.К., Чигирева О.Ю. Температурное поле цилиндрического тела в режиме периодического разогрева // Вестник Московского государственного технического университета им. Н.Э. Баумана. Серия: Естественные науки. 2015. № 3 (60). С. 88-98.
7. Dorofeev K.Y., Nikolaeva N.N., Titarenko V.N., Yagola A.G. New approaches to error estimation to Ill-posed problems with application to inverse problems of heat conductivity // Journal of Inverse and Ill-posed problems. 2002. Vol. 10. No 2. P.155-169.
8. Танана В.П. Об оценке погрешности метода решения одной обратной задачи для параболического уравнения // Сибирский журнал вычислительной математики. 2010. Т. 13, № 4. С. 451-465.
9. Табаринцева Е.В. О решении граничной задачи для параболического уравнения методом вспомогательных граничных условий // Вестник Южно-Уральского государственного университета. Серия: Математика. Механика. Физика. 2011. №32 (249). С. 68-76.
10. Танана В.П., Гайнова И.А., Сидикова А.И. Об оценке погрешности приближенного решения одной переопределенной обратной задачи тепловой диагностики // Сибирский журнал индустриальной математики. 2012. Т. XV, № 1. С. 145-154.
11. Кабанихин С.И., Шишленин М.А. Прямые и итерационные методы решения обратных и некорректных задач // Сибирские электронные математические известия. 2008. Т. 5. С. 595-608.
12. Zhang Y., Lukyanenko D.V., Yagola A.G. Using Lagrange principle for solving two-dimensional integral equation with a positive kernel // Inverse Problems in Science and Engineering. 2015. DOI: 10.1080/17415977.2015.1077445.
13. Солодуша С.В., Япарова Н.М. Численное решение обратной граничной задачи теплопроводности с помощью уравнений Вольтерра I рода // Сибирский журнал вычислительной математики. 2015. Т. 18, № 3. С. 327-335.
14. Дрозин А.Д., Дудоров М.В., Рощин В.Е., Гамов П.А., Менихес Л.Д. Математическая модель образования кристаллическиз зардышей в переохлажденном расплаве эвтектического сплава // Вестник Южно-Уральского государственного университета. Серия: Математика. Механика. Физика. 2012. № 11 (270). С. 66-77.
15. Марчук Г.И. Методы вычислительной математики. СПб.: Лань. 2009. 608 с.
16. Вабищевич П.Н. Монотонные разностные схемы для задач конвекции / диффузии // Дифференциальные уравнения. 1994 Т. 30, №3. С. 503-515.
17. Тихонов А.Н., Самарский А.А. Уравнения математической физики. Изд-во МГУ, 1999. 799 с.
18. Камонт З., Кропельницка К. Неявные разностные методы для эволюционных функционально-дифференциальных уравнений // Сибирский журнал вычислительной математики. 2011. Т. 14, № 4. С. 361-379.
19. Глазырина О.В., Павлова М.Ф. Исследование сходимости метода конечных элементов для решения параболических уравнений с нелинейным нелокальным пространственным оператором // Дифференциальные уравнения. 2015. Т. 51, №7. С. 876.
20. Япарова Н.М. Численный метод решения некоторых обратных задач теплопроводности с неизвестными начальными условиями // Вестник ЮжноУральского государственного университета. Серия: Компьютерные технологии, управление, радиоэлектроника. 2015. Т. 15, № 2. С. 55-65.
21. Япарова Н.М. Метод решения некоторых многомерных обратных граничных задач параболического типа без начальных условий // Вестник ЮжноУральского государственного университета. Серия: Компьютерные технологии, управление, радиоэлектроника. 2015. Т. 15, № 2. С. 97-108.
22. Самарский А.А. Введение в теорию разностных схем. М.: Наука, 1971. 552 с.
23. Ладыженская О.А., Солонников В.А., Уральцева Н.Н. Линейные и квазилинейные уравнения параболического типа. M.: Наука. 1967. 736 с.
24. Тихонов А.Н., Гончарский А.В., Степанов В.В., Ягола А.Г. Численные методы решения некорректных задач. Изд-во МГУ, 1990. 115 с.
Япарова Наталья Михайловна, к.ф.-м.н., доцент, кафедра прикладной математики, Южно-Уральский государственный университет (Челябинск, Российская Федерация), [email protected]
Поступила в редакцию 30 декабря 2015 г.
Bulletin of the South Ural State University Series "Computational Mathematics and Software Engineering"
2016, vol. 5, no. 2, pp. 43-58
DOI: 10.14529/cmse160204
NUMERICAL METHOD FOR SOLVING AN INVERSE PROBLEM FOR NONLINEAR PARABOLIC EQUATION WITH UNKNOWN INITIAL CONDITIONS
N.M. Yaparova, South Ural State University, Chelyabinsk, Russian Federation
The paper is devoted to the inverse problem for a nonlinear parabolic equation with unknown initial conditions. A computational scheme for solving this problem is proposed. This approach allows obtain the numerical solution in internal points of domain and the unknown boundary function. The proposed scheme is based on the using of finite-difference equations and regularization technique. We investigate the stability of computational method. We obtained the dependence of stability on the discretization steps and level error of the initial data The proposed scheme proved the basis for development of numerical method and for the computational experiment. The experimental results are also presented in this paper, and confirm the effectiveness of the method.
Keywords inverse problem, numerical method, regularization method, error estimate, computational scheme.
FOR CITATION
Yaparova N.M. Numerical Method for Solving an Inverse Problem for Nonlinear Parabolic Equation with Unknown Initial Conditions. Bulletin of the South Ural State University. Series: Computational Mathematics and Software Engineering. 2016. vol. 5, no. 2. pp. 43-58. (in Russian) DOI: 10.14529/cmse160204.
References
1. Alifanov O.M. Obratnye zadachi teploobmena [Inverse Heat Transfer Problems]. Moscow, Mashinostroenie, 1988. 280 p. (in Russian)
2. Samarskij A.A., Vabishchevich P.N. Vychislitel'naya teploperedacha [Computational Heat Transfer]. Moscow, Editorial URSS, 2003. 784 p. (in Russian)
3. Prokudina L.A. Vyatkin G.P Self-Organization of Perturbations in Fluid Films. Dokla-dy Akademii nauk [Doklady Physics]. 2011. vol. 56, no. 8. pp. 444-447.
4. Bulgakova G.T., Kondrat'eva N.R. Analiticheskaya model' vertikal'nogo vytesneniya nefti vodoj s uchetom vyazkostnyh, gravitacionnyh i kapillyarnyh sil [Analytical Model of Vertical Oil-Water Displacement with the Account of Viscous, Capillary and Gravity Forces]. Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta. seriya: Fizi-ko-matematicheskie nauki [Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences]. 2012. no. 1. pp. 208-213. (in Russian)
5. Shestakov A.L. Metody teorii avtomaticheskogo upravleniya v dinamicheskih izmereni-yah [The Methods of Control Theory in Dynamic Measurements]. Chelyabinsk, Publishing of the South Ural State University, 2013. 257 p. (in Russian)
6. Martinson L.K., Chigiryova O.Y. Temperaturnoe pole cilindricheskogo tela v rezhime pe-riodicheskogo razogreva [Thermal Fields of a Cylindrical Body During Cyclic Heating].
H.M. 5InapoBa
Vestnik Moskovskogo gosudarstvennogo tekhnicheskogo universiteta im. N.EH. Baumana. Seriya: Estestvennye nauki [Herald of the Bauman Moscow State Technical University, Series: Natural Sciences]. 2015. no. 3(60). pp. 88-98. (in Russian)
7. Dorofeev, K.Y., Nikolaeva N.N., Titarenko V.N., Yagola A.G. New Approaches to Error Estimation to Ill-Posed Problems with Application to Inverse Problems of Heat Conductivity. Journal of Inverse and Ill-Posed Problems. 2002. vol. 10, no. 2. pp. 155-169.
8. Tanana V.P. An order-optimal method for solving an inverse problem for a parabolic equation. Sibirskij zhurnal vychislitelnoj matematiki [Numerical Analysis and Applications]. 2010. vol. 3, no. 4. pp. 367-371.
9. Tabarinceva E.V. O reshenii granichnoj zadachi dlya parabolicheskogo uravneniya metodom vspomogatel'nyh granichnyh uslovij [About Solution of the Boundary Inverse Problem for Parabolic Equation by Means of Subsidiary Boundary Conditions Method]. Vestnik Yuzho-Uralskogo gosudarstvennogo universiteta. Seriya "Matematika. Mekhani-ka. Fizika" [Bulletin of the South Ural State University, Series: Mathematics. Mechanics. Physics]. 2011. no. 32(249). pp. 68-76. (in Russian)
10. Tanana V.P., Gajnova I.A., Sidikova A.I. Ob otsenke pogreshnosti priblizhennogo resh-eniya odnoy pereopredelennoy obratnoy zadachi teplovoy diagnostiki [The Estimation of Approximate Solutions an Overdetermined Inverse Problem of Thermal Diagnostics]. Sibirskij zhurnal industrialnoj matematiki [Journal of Applied and Industrial Mathematics]. 2012. vol. XV, no. 1. pp. 145-154. (in Russian)
11. Kabanihin S.I., Shishlenin M.A. Pryamye i iteratsionnye metody resheniya obratnykh i nekorrektnykh zadach [Direct and Iterati on Methods for Solving Inverse and Ill-Posed Problems]. Sibirskie ehlektronnye matematicheskie izvestiya [Siberian Electronic Mathematical Reports]. 2008. vol. 5, pp. 595-608. (in Russian)
12. Zhang Y., V. Lukyanenko D.V., Yagola A.G. Using Lagrange Principle for Solving Two-Dimensional Integral Equation with a Positive Kernel. Inverse Problems in Science and Engineering. 2015. DOI: 10.1080/17415977.2015.1077445.
13. Solodusha S.V., Yaparova N.M. Numerical Solving an Inverse Boundary Value Problem of Heat Conduction Using Volterra Equations of the First Kind. Sibirskij zhurnal vychislitelnoj matematiki [Numerical Analysis and Applications]. 2015. vol. 8, no. 3. pp. 267274.
14. Drozin A.D., Dudorov M.V., Roshchin V.E., Gamov P.A., Menihes L.D. Matematich-eskaya model' obrazovaniya kristallicheskiz zardyshey v pereokhlazhdennom rasplave ev-tekticheskogo splava [Mathematical Description of the Nucleation in Supercooled Eutec-tic Melt]. Vestnik Yuzho-Uralskogo gosudarstvennogo universiteta. Seriya "Matematika. Mekhanika. Fizikd' [Bulletin of the South Ural State University, Series: Mathematics. Mechanics. Physics]. 2011. no. 11(270). pp. 66-77. (in Russian)
15. Marchuk G.I. Metody vychislitel'noj matematiki [Methods of Computational Mathemat-ics].St. Petersburg, Publishing of the Lan', 2009. 608 p. (in Russian)
16. Vabishchevich P.N. Monotone Finite-Difference Schemes on Triangular Grids for Convection-Diffusion Problems. Differencialnye Uravneniya [Differential Equations]. 1994. vol. 42, no. 9. pp. 1317-1330.
17. Tihonov A.N. Samarskij A.A. Uravneniya matematicheskoj fiziki [Equations of Mathematical Physics]. Moscow, Publishing of the Moscow State University, 1999. 799 p. (in Russian)
18. Kamont Z., Kropelnicka K. Neyavnye raznostnye metody dlya evolyutsionnykh funktsional'no-differentsial'nykh uravneniy [Implicit Difference Methods for Evolution Functional Differential Equations]. Sibirskij zhurnal vychislitelnoj matematiki [Numerical Analysis and Applications]. 2011. vol. 14, no. 4. pp. 361-379. (in Russian)
19. Glazyrina O.V. Pavlova M.F. Study of the Convergence of the Finite-Element Method for Parabolic Equations with a Nonlinear Nonlocal Spatial Operator. Differencialnye Uravneniya [Differential Equations]. 2015. vol. 51, no. 7. pp. 872-885.
20. Yaparova N.M. Chislennyy metod resheniya nekotorykh obratnykh zadach teploprovod-nosti s neizvestnymi nachal'nymi usloviyami [Numerical Method for Solving Some Inverse Heat Conduction Problems with Unknown Initial Conditions]. Vestnik Yuzho-Uralskogo gosudarstvennogo universiteta. Seriya "Komp'yuternye tekhnologii, upravlenie, radioehlektronika" [Bulletin of the South Ural State University, Series: Computer Technologies, Automatic Control and Radioelectronics]. 2015. no. 2. pp. 55-65. (in Russian)
21. Yaparova N.M. Metod resheniya nekotorykh mnogomernykh obratnykh granichnykh zadach parabolicheskogo tipa bez nachal'nykh usloviy [Method for Solving Some Multidimensional Inverse Boundary Value Problems for Parabolic PDEs Without Conditions]. Vestnik Yuzho-Uralskogo gosudarstvennogo universiteta. Seriya "Komp'yuternye tekhnologii, upravlenie, radioehlektronika" [Bulletin of the South Ural State University, Series: Computer Technologies, Automatic Control and Radioelectronics]. 2015. no. 3. pp. 97-108. (in Russian)
22. Samarskij A.A. Vvedenie v teoriyu raznostnyh skhem [The Theory of Difference Shemes]. Moscow, Nauka, 1971. 552 p. (in Russian)
23. Ladyzhenskaya O.A., Solonnikov V.A., Uralceva N.N. Linejnye i kvazilinejnye uravneniya parabolicheskogo tipa [Linear and Quasi-Linear Parabolic Equations]. Moscow, Nauka, 1967. 736 p. (in Russian)
24. Tihonov A.N., Goncharskij A.V., Stepanov V.V., YAgola A.G. Chislennye metody resh-eniya nekorrektnyh zadach [Numerical Methods for the Solution of Ill-Posed Problems]. Moscow, Publishing of the Moscow State University, 1990. 115 p. (in Russian)
Received December 30, 2015.
