CC BY
37
УДК 519.6
Об устойчивом продолжении решений уравнения теплопроводности с неточно заданной границы
Е. Б. Ланеев, М. Н. Муратов, Адель Салех Абдулхак Табет
Кафедра нелинейного анализа и оптимизации
Российский университет дружбы народов ул.Миклухо-Маклая, 6, 117198, Москва, Россия
Рассматривается задача продолжения решения уравнения теплопроводности с неточно заданной границы как некорректно поставленная задача Коши для уравнения теплопроводности с данными Коши на приближённо заданной поверхности.
Ключевые слова: некорректная задача, задача Коши, уравнение теплопроводности, метод регуляризации.
1. Введение
Продолжение поля ньютоновского потенциала давно зарекомендовало себя как эффективный инструмент исследования его плотности по косвенной информации [1]. Хорошо известна связь этой задачи с задачей Коши для уравнения Лапласа [2]. Существенной проблемой при решении этих задач является их некорректность, которая затрудняет непосредственное использование формул «точного» решения для вычислений. Эта проблема принципиально решена построением теории методов регуляризации некорректно поставленных задач [3]. Тем не менее, практические задачи требуют усложнения модели. В работах [4-6] рассматривается задача Коши для уравнения Лапласа с данными на поверхности общего вида в отличие от [1], где граница плоская. При этом решение задачи строится таким образом, что расчётные формулы сводятся к рядам Фурье, как и в [1], и вычислительные алгоритмы в основе своей сохраняют стандартную структуру, ставшую уже классической. Вместе с тем задача требовала следующего естественного обобщения постановки. В самом деле, «неплоская» поверхность в прикладных задачах сама представляет собой результат измерений, конечно, содержащих погрешность, и, таким образом, возникает задача продолжения гармонической функции с поверхности, заданной приближённо. Следует отметить, что метод построения решения задачи продолжения с произвольной поверхности предполагает вычисление нормали к поверхности, что приводит к необходимости решения частной задачи устойчивого построения нормали к поверхности, заданной приближённо. В [7-9] построено устойчивое решение как такой частной задачи, так и задачи построения продолжения решения. В данной работе эти результаты распространяются на некорректные задачи Коши для уравнения теплопроводности [10].
2. Постановка задачи
Пусть в теплопроводящем теле цилиндрической формы прямоугольного сечения имеются источники тепла плотности Будем считать, что на боковых гранях цилиндра поддерживается нулевая температура, а на поверхности 5 поддерживается конвективный теплообмен со средой нулевой температуры и. Будем
Статья поступила в редакцию 8 ноября 2010 г.
Работа выполнена при поддержке Российского фонда фундаментальных исследований (РФФИ грант №09-01-00633^).
для простоты считать, что начальная температура равна нулю. Получим смешанную краевую задачу для уравнения теплопроводности
ди
ди дп
^(М,Ь) = а2Аи(М,Ь) + р(М,Ь), М е Б(Р, то), 0 < Ь < то Н(и - и) ,
где
и\х=0,1х = о, и1У=0,1у = 0, и\г=0 = 0 и ^ г ^ то
Б(Р,Н) = {(х,у,г) : 0 <х <1х, 0 <у<1у, ^ (х,у) <г<Н}, Б = {(х, у,г) : 0 < х <1х, 0 < у < 1у, х = Р(х, у)}, Р е С2(П(0)), П(г) = {(х, у,г):0 <х<1х, 0 <у<1у, г = со^}.
(1)
(2)
Пусть теперь плотность источников р неизвестна и подлежит определению. Заданной (измеренной) будем считать функцию
и1з = /, 0 <г < то
(3)
В области И(Р, Н) ® (0, то), считая, что носитель плотности источников расположен в области г > Н, получаем задачу
(М,€) = а2Аи(М,г), М е Б(Р,Н), -то <г< то и\в = /,
ди дп
= = g,
в
(4)
и
х=0,1Х = ° и\у=0,1у и\ 4=0 = 0 и —у х —у то
Функции f и д считаем непрерывными на 5 (8> (0, то). Задача (4) некорректно поставлена [5] по условиям Коши на поверхности 5. В [5] приведён метод построения точного и регуляризованного решения аналогичной задачи Коши для уравнения Лапласа, устойчивого по отношению к погрешностям в функциям f и д. Основной элемент этой схемы — сведение задачи к интегральному уравнению.
3. Точное решение
Применением формул Грина приближенное решение задачи (4) строится в виде
д+г ж
1(М,1) = У(М,€) - Ф(М,€) = £ ^ I еР'Фпт(Ь,р)х
п,т= 1
х е~
р+х2 2
(2 + ¿2
( 'ар 81П--- 81П
ъпхм ■ ктум
1у
- Ф(м,г). (5)
в
в
0
х
где
С
${M,t) = j j[g(P,r)V(M,P,t,r) - f{P,r)^-p (M,P,t,r)]dadr, (6)
0 s
Qit-T) — р )2 — a2n
) = l^twk' ^ £ e
V v / n,m= 1
n__i m
L '£ + 'У J
(í-r)
. ■кпхм . ктум . "nnxp . жтур х sin —-— sin —--sin —-— sin —-—.
ix iy ix Ly
Здесь Фпт(Ь,р) — образ коэффициентов Фурье функции Ф(М,Щмеп(б)
С 1у
~ I 'К'Л/Х
$пт(Ь,р)= dte~р ТТ dx dy $(x,y,b,t)sin-¡— sin .
J *>x(>y J J iX iy
00
Перепишем поверхностный интеграл в (6) в виде интеграла по прямоугольнику П(0), используя нормаль к поверхности 5 вида (2):
Ф(м,г) = I I [9(р,тмм,р,г,т)-
0 П(0)
- /(Р, т)(Ур<р(М, Р, г, т), п(Р))] щ (Р)6хрАур&Г, (7)
где
П1 = %г&А(Р(х,у) - ¿) = ,-1} = + ^ - к, п = ^. (8)
Отсюда получаем
Ф(М,1)= ! у [д(Р,т)^(М,Р,1,т)т(Р)-
0 П(0)
- /(Р, Т)(Ур<р(М, Р, I, т), П1 (Р))] АхрАур&т. (9)
Таким образом, для вычисления функции Ф, на основе которой формируется решение (5) задачи (4), необходимо вычислить нормаль П1 вида (8) к поверхности 5.
4. Приближенное решение
В том случае, когда функция Р, определяющая поверхность 5, известна с некоторой погрешностью, задача вычисления градиента этой функции — некорректно поставлена. Для получения её устойчивого решения воспользуемся постановкой [7-9], то есть рассмотрим задачу вычисления нормали к поверхности (градиента приближённо заданной функции) как задачу восстановления значений неограниченного оператора [11].
Будем считать, что вместо точной функции Р задана функция Р'л такая, что
11^ - ^ 1к(п(0)) < (10)
В качестве приближения к функции УР, вычисляемого по известной функции Р^, связанной с Р условием (10), будем рассматривать градиент от экстремали функционала
т = - + Р ||^2(п(0)) . (11)
При этом будем рассматривать такие поверхности для которых Р|ш=о, 1Х = Р|у=о, 1У = 0. Это условие, в частности, имеет место в случае, когда поверхность
5 можно рассматривать как «возмущение» плоскости г = 0.
Как и в [8, 9], в качестве приближенного значения градиента функции Р1Л
будем рассматривать вектор-функцию
v * у w>;(x,y)= £ — ^
п,ш=1 1 + /3
пт
2 / \2 ■кп \ | / пт \
^х J у ly J
кпх . -кшу -кшу . кпх\ х I i— cos—,— sm —--+ cos—,— sm —— I . (12)
у ^X by ly ly lx J
Ряд (12) равномерно сходится на П(°).
В [8,9] показано, что если F £ С2 и fi(^) = ,, ^ ,, , то для приближенного
Ар
вектора нормали
n t = ' = v * у wfc ' - k (13)
имеет место оценка
К - ^^(ЩО)) < VWAF ^ (14)
Пусть теперь функции f и д заданы приближённо, а именно: пусть заданы функции f5 и д5, такие, что
max II fs — f || r m,n,, < ё, max \\qs — а II T < S. (15)
t 11 ^ J ИЬ2(П(0)) - ' t ^2(П(0)) v '
Тогда функция Ф вида (9) при точно заданной поверхности S может быть вычислена с некоторой погрешностью как функция Ф*5:
Ф5(М,1) = I I уд (Р,тММ,Р,г,т)щ(Р)-
0 П(0)
- /5(Р,т)(ур<р(м,р,г,т),П1 (Р))]ёЖрйурйт, Р е (16)
При приближенном задании поверхности 5, определяемом условием (10), функция Ф может быть вычислена приближённо по формуле (16) с использованием (13):
СЮ
Ф5,,л(М,€)= I I [д5(Р,тММ,Р,г,т)п»(Р)-
0 П(0)
- /5(Р,т)(Ур<р(М,Р,1,т), П*(Р)) }йхрйурйт, Р е Б. (17)
X
Оценим разность, считая, что t £ (0, Т)
Ф5'>л(M,t) - Ф(М,$
<
<
Ф5'1(м,г) - ф5(M,t) + ф5(M,t) - Ф(м,г) , м е п(ь), (18)
где функции Фг, Ф — функции вида (17), (16), (9) соответственно. Оценим первую разность в (18), используя неравенство Коши-Буняковского:
ф5'i\M,t) - Ф5(м,г)
f (Р,т)(п»(Р) - т(Р), v рcp(M, p, t,т))-
0 П(0)
- д5(Р, т)<р(М, Р, t, т)(n't(P) - m(P))
р es
dxpdypdr
<
<
f5 (Р,т)| ■ \nt(P) - n1(P)| ■ | Vp<p(M,P,t,r)| +
0П(0)
+ | g5(P,T) | ■ | <p(M,P,t,r) | ■ | nt(P) - щ(Р)
p es
dxPdyPdr <
< ¡Uni, I vp«'{M<p<t-T)| ■ II;
0 р e s + max | ip(M,P,t,r) | ■ ||g5
1^2(П(0))
+
м en(o)
p es
L2 (п(0)) I dr 11 n! ni|L2(n(0)) ^
(мтаь| v p ^(м,р,',т) | I|^2(n(0))+*) +
м en(o) 0 \ pes
+ MmeaaXo) | V(M,p,f,T) | {Ы Ь2(П(0)) + ) 11 < - n
p es
11 n4lb2 (n(0)).
Так как нас интересует поведение регуляризованного решения задачи (1) при S ^ 0, то можно считать, что S ^ ¿0, и, таким образом, с учётом (14)
| Ф5'1 (M,t) - Ф5 (M,t) | <
< Const Щ - ni ||МЩ0)) < Ci M е П(Ь), 0 <t<T. (19)
Для оценки второй разности получаем — так же, как в [5]:
Ф5(M,t) - Ф(М,Ь)
(f5(Р, г) - f (Р, Т))(ni (Р), Vp<р(М, Р, t, т))-
0П(0)
- (д 5 (Р,т) - д(Р,т )MM,P,t,T )m(P)
p es
dxpdypdr
<
/ mais HVv(M,t,r)H
L2(s) ■ 11 f5 - f 1 l2(n(0)) dr+
5
+ GonstiJ Mmax ИМ^т)\^) ■ \\gs — з\^2(Щ0)) dr < C2Ö, (20)
о
Из (19) и (20) для оценки (18) получаем:
тах |Ф5''(М,г) - Ф(М,г)\ < Сг^р + С26 = А(р,5) -► 0. (21)
М €П(а)
Таким образом, функция Ф известна с некоторой погрешностью А, имеющей структуру (21). В соответствии со схемой [5] устойчивое приближенное решение задачи (1) строится введением регуляризирующего множителя [11] в (5) и может быть получено в виде
и^(м,г) = - Ф5'^(м,г), м е о(Р,н), о <г<т, (22)
где Ф5'' — функция вида (17), а V"'5 имеет вид:
д+г ж
оо
/М(м t)- V — f ppt_Ф™(ь,р)
V'Z (M,t)- у — i е
4 1 + aexp(2^-f--^(b — H))
I 2 2
x e~
(—dp sin ^ sin . (23)
lX ly
Здесь Ф^(b) — коэффициенты Фурье функции Ф^'*5(M,í)Imgn(ь):
CO lx ly
фtä (b,p) - í dte-pt -¡4- fdx f dyФ/'s (x,y,b,t)sin ^ Sin , (24)
J bXl>y J J lx (>y
0 0 0
а a — параметр регуляризации. В соответствии с обозначениями, введёнными выше, величина b выбирается так, чтобы b < min F(х,у).
(х,у)€П(0)
Формулируя теорему [10] для погрешности (21) в задании функции Ф, получим следующую теорему.
Теорема 1. Пусть решение задачи (4) существует в области D(F,H) ® (0, то), а — а(А), а(А) ^ 0, А^а(А) ^ 0 при А ^ 0. Тогда функция иа(А) вида (5), где А — C1yrfi + C2Ö, равномерно сходится к точному решению задачи при 6 ^ 0, ß ^ 0 в области D(F + е,Н — е) ® (0, то), где е > 0 — некоторое фиксированное сколь угодно малое число.
х
5. Заключение
Приближенное решение вида (22) может быть использовано для построения вычислительных алгоритмов решения задачи (4), или продолжения нестационарного температурного поля с поверхности, заданной приближённо. Отметим, что продолжение осуществляется в область г > Н, где источники (неизвестной плотности) отсутствуют и уравнение теплопроводности однородно. Получив решение (температурное поле) продолжением в область, непосредственно примыкающую к источникам, можно получить более или менее достоверную информацию об их
плотности. Как правило, эта информация формируется особенностями продолженного решения. Следует иметь в виду, что в прикладных задачах измеряемая функция f = u|s представляет собой тепловизионный снимок, уже содержащий «изображение» плотности источников. Предполагается (и это подтверждается вычислительным экспериментом), что продолженное поле даёт уточнённое изображение. Продолжение нестационарного температурного поля даёт уточнение эволюции изображения плотности источников во времени.
Литература
1. О продолжении потенциала в сторону возмущающих масс на основе метода регуляризации / А. Н. Тихонов, В. Б. Гласко, О. К. Литвиненко, В. Р. Мелихов // Изв. АН СССР. Физика Земли. — 1968. — Т. 1. — С. 38-40. [O prodolzhenii potencíala v storonu vozmuthayuthikh mass na osnove metoda regulyarizacii / A. N. Tikhonov, V. B. Glasko, O. K. Litvinenko, V. R. Melikhov // Izv. AN SSSR. Fizika Zemli. — 1968. — T. 1. — S. 38-40.]
2. Лаврентьев М. М. О некоторых некорректных задачах математической физики. — Новосибирск: СО АН СССР, 1962. — 92 с. [Lavrentjev M. M. O nekotorihkh nekorrektnihkh zadachakh matematicheskoyj fiziki. — Novosibirsk: SO AN SSSR, 1962. — 92 s.]
3. Тихонов А. Н., Арсенин В. Я. Методы решения некорректных задач. — М.: Наука, 1979. — 288 с. [Tikhonov A. N., Arsenin V. Ya. Metodih resheniya nekorrektnihkh zadach. — M.: Nauka, 1979. — 288 s.]
4. Ланеев Е. Б. О задаче Коши для уравнения Лапласа в неодносвязной области // Статистическая и квантовая физика и ее приложения. Сборник научных трудов. — М.: Изд-во УДН, 1968. — С. 49-56. [Laneev E. B. O zadache Koshi dlya uravneniya Laplasa v neodnosvyaznoyj oblasti // Statisticheskaya i kvantovaya fizika i ee prilozheniya. Sbornik nauchnihkh trudov. — M.: Izd-vo UDN, 1968. — S. 49-56.]
5. Ланеев Е. Б., Бхувана Васудеван. Об устойчивом решении одной смешанной задачи для уравнения Лапласа // Вестник РУДН. Серия «Прикладная математика и информатика». — 1999. — Т. 1. — С. 128-133. [Laneev E. B., Bkhuvana Vasudevan. Ob ustoyjchivom reshenii odnoyj smeshannoyj zadachi dlya uravneniya Laplasa // Vestnik RUDN. Seriya «Prikladnaya matematika i informatika». — 1999. — T. 1. — S. 128-133.]
6. Ланеев Е. Б. О некоторых постановках задачи продолжения потенциального поля // Вестник РУДН. Серия «Физика». — 2000. — Т. 8(1). — С. 21-28. [Laneev E. B. O nekotorihkh postanovkakh zadachi prodolzheniya potencialjnogo polya // Vestnik RUDN. Seriya «Fizika». — 2000. — T. 8(1). — S. 21-28.]
7. Ланеев Е. Б. О регуляризации некоторых операций векторного анализа // Методы функционального анализа в математической физике. Сборник научных трудов. — М.: Изд-во УДН, 1987. — С. 101-106. [Laneev E. B. O regulyarizacii nekotorihkh operaciyj vektornogo analiza // Metodih funkcionaljnogo analiza v matematicheskoyj fizike. Sbornik nauchnihkh trudov. — M.: Izd-vo UDN, 1987. — S. 101-106.]
8. Ланеев Е. Б., Муратов М. Н. Об устойчивом решении одной смешанной краевой задачи для уравнения Лапласа с приближенно заданной границей // Вестник РУДН. Серия «Математика». — 2002. — Т. 9(1). — С. 102-111. [Laneev E. B., Muratov M. N. Ob ustoyjchivom reshenii odnoyj smeshannoyj kraevoyj zadachi dlya uravneniya Laplasa s priblizhenno zadannoyj graniceyj // Vestnik RUDN. Seriya «Matematika». — 2002. — T. 9(1). — S. 102-111.]
9. Ланеев Е. Б., Муратов М. Н. Об одной обратной задаче к краевой задаче для уравнения Лапласа с условием третьего рода на неточно заданной границе // Вестник РУДН. Серия «Математика». — 2003. — Т. 10(1). — С. 100-110. [Laneev E. B., Muratov M. N. Ob odnoyj obratnoyj zadache k kraevoyj zadache
dlya uravneniya Laplasa s usloviem tretjego roda na netochno zadannoyj granice // Vestnik RUDN. Seriya «Matematika». — 2003. — T. 10(1). — S. 100-110.]
10. Ланеев Е. Б., Муратов М. Н., Табет Адель Салех Абдулхак. Задача продолжения нестационарного температурного поля с произвольной поверхности // Вестник РУДН. Серия «Математика. Информатика. Физика». — 2010. — Т. 2(1). — С. 109-111. [Laneev E. B., Muratov M. N., Tobet Adelj Salekh Abdulkhak. Zadacha prodolzheniya nestacionarnogo temperaturnogo polya s proizvoljnoyj poverkhnosti // Vestnik RUDN. Seriya «Matematika. Informatika. Fizika». — 2010. — T. 2(1). — S. 109-111.]
11. Морозов В. А. Об одном устойчивом методе вычисления неограниченных операторов // ДАН СССР. — 1969. — Т. 185(2). — С. 267-270. [Morozov V. A. Ob odnom ustoyjchivom metode vihchisleniya neogranichennihkh operatorov // DAN SSSR. — 1969. — T. 185(2). — S. 267-270.]
UDC 519.6
On Stable Continuation of the Heat Conduction Equation Solution from Approximately Defined Boundary
E. B. Laneev, M. N. Mouratov, Adel Saleh Abdulhak Tabet
Department of Differential Equations and Functional Analysis Peoples' Friendship University of Russia Miklykho-Maklaya str., 6, 117198, Moscow, Russia
A problem of stable continuation of the heat conduction equation solution is considered. The continuation of the solution is performed from approximately defined boundary. The whole problem is treated as an incorrect Cauchy problem for the heat conduction equation with Cauchy data given on an arbitrary approximately defined surface.
Key words and phrases: incorrect problem, Cauchy problem, heat conduction equation, regularization.
