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

УДК 519.6+517.9

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ РЕШЕНИЙ ОБРАТНОЙ ГРАНИЧНОЙ ЗАДАЧИ ТЕПЛОПРОВОДНОСТИ

Н.М. Япарова

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

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

Введение

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

В статье предложены методы решения обратной задачи теплопроводности с постоянной границей, основанные на применении преобразований Лапласа и Фурье.

Преобразования Фурье по пространственной переменной нашло широкое применение в области решения обратных задач теплопроводности. Так, в работе [1] для получения регуля-ризованных решений используют преобразование Фурье и стабилизирующий функционал. В работе [2] с целью регуляризации и оценки сходимости построенных решений были применены преобразование Фурье и метод сопряженных градиентов.

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

Применения преобразования Лапласа для решения обратных задач теплопроводности также вызывает большой интерес [3, 4]. В существующих подходах, как правило, после выполнения преобразования Лапласа для полученных уравнений применяются методы регуляризации, а затем для регуляризованного решения выполняется обратное преобразование.

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

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

1. Постановка задачи. Построение операторного уравнения с использованием преобразования Лапласа

Рассмотрим обратную задачу теплопроводности:

«[ = «хх, х € (0,1), Ь > 0, (1)

«(х, 0) = 0, щ(0, Ь) = 0, х € (0,1), Ь > 0,

«х(0, Ь) = д(Ь), Ь > 0.

Известно, что при д(Ь) = до(Ь) существует точное решение «о(Ь). Но вместо до даны д$ и 5 > 0 такие, что ||д^ — до У — 5. В этой задаче требуется найти граничное значение функции

щ(1, Ь) = ^(Ь), Ь > 0. (2)

Сведем исходную задачу (1), (2) к интегральному уравнению. С этой целью решаем прямую задачу, предположив, что функция щ(1,Ь) = <^(Ь) известна:

Щ = «хх, х € (0,1), Ь > 0, (3)

щ(х, 0) = 0, х € (0,1),

«(0,Ь) = 0, «(1,Ь) = ^(Ь), Ь > 0.

Решив эту задачу, получим уравнение, связывающее функции «(1,Ь)= ^(Ь) и «х(0,Ь) = д(Ь).

Если «(х,Ь) € ^([0,1] х [0, то)), ЗС > 0, во > 0 такие, что УЬ € [0, то) выполнено

|^(Ь)| > Сево[, и, кроме того, при любом Т = то и УЬ € [0,Т] функция ^(Ь) удовлетворяет

условиям Дирихле, то для решения прямой задачи можно использовать преобразование Лапласа [5]. Из результатов, изложенных в работе [6] следует, что при |х| = 1 решение задачи (3) определяется формулой:

^ [

«(х, Ь) = 2п ^2(—1)п+1п 8т(ппх)е_п2п2[ / ^(т)еп2п2т ^т. (4)

п= о

Если при любом Т = то ЗМ > 0 такое, что УЬ € [0,Т] выполнено ||^(Ь)||я2([о,Т]) — М,

и ^>;(0) = 0, то для того, чтобы получить соотношение, связывающее «(1,Ь) = <^(Ь) и

«х(0, Ь) = д(Ь) используем следующий подход. Рассмотрим частичную сумму ряда (4):

N [

«М (х, Ь)=2п^(—1)п+1 п 81п(ппх)е-п2га2^ ^(т)еп2п2т ^т. п= о

Интегрируя по частям, получим, что «м(х,Ь) имеет вид:

« 2 ^ ( -Пп+1^ 8Ш(ппх) 8Ш(ппх) _/(ЬЛ , э1п(ппх) е_Л2* [ ,, (_ЛеЖ2п2т ^

«м(х,Ь) = 2}^( —1) + пп 1(Ь) — П3П3 ^ (Ь) + п3п3 е V (т)е йт •

п=1 \ о )

Тогда

«м(х,Ь) = 2 V(—1)п+181п(ппх) ^(Ь)+

пп

п=1

+2 V (—1)п+^ ££(=> ^ (0 — ^

“ \ п3п3 п3 п

п=1 \

7 1/^7 (3Пп3х) е_п2п2[ / V'(т)е^пт ^т ).

п3п3 п3 п3 J

о

Имеет место соотношение:

«(х, Ь) = (х, Ь)+

+2 V ( —1)п+^ 51п(ппх) ^(г) — 81п(П’3х) V' («) + 81п(!Г’3х) е-п2"2. / V''(т )е*1п>т(

^ \ пп п3п3 п3 п3 У

п=М+1 \ о

Следуя идее, предложенной в работе [7], и учитывая сходимость ряда

[

^(—1)4 8111(3пта3х) V'(Ь) — 8111(3пта3х) е_п2'п2[ I V'(т)е^пт ^т

п= 1

получаем при N ^ то

N ' [

»(х, () = х^) + 2 £(—1)п (V' М — е_п2п2[[ / V' '(т )еп2п2т йт ] + « (!).

Тогда в качестве приближенного решения задачи (3) рассмотрим функцию, заданную формулой:

N ' [

и(х, Ь) = xv(^) + 2 ^(—1)п ( 8Шп(зп;3х) V'(Ь) — 81пп(зп;3х) е_п2п2^ V''(т)еп2п2т ^т ) . (5)

п=1 \

Продифференцируем (5) по переменной х, тогда при фиксированном N получим, что функция «х(0,Ь) = д(Ь) определяется формулой:

N / '(,) _п2п2. [

д(Ь) = «х(0,Ь) = V(г) + 2 ^(—1)п+М ^-2 — V''(т)е

п=1 V п

Выполнив преобразования и принимая во внимание то, что для обеспечения сходимости регуляризованных решений необходима нечетность числа N, мы получаем следующее уравнение:

^ = —V(t) +2 ;£(—:!)^п2п2 I v(т)е_п2п2([_т) ^ = д(Ь). (6)

При этом вместо точных значений функции д(Ь) нам известны приближения д<$ (Ь) и 5 > 0 такие, что при любом Ь > 0 выполнено ||д(Ь) — до(Ь) | — 5. Требуется по д^(Ь) и 5 > 0 найти V(^) = «(1, Ь).

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

2. Метод проекционной регуляризации

Пусть Н - гильбертово пространство. А : Н ^ Н - линейный, инъективный, ограниченный оператор. А1 = А*А, спектр £р(А1) = [0, || А||2] . Рассмотрим операторное уравнение

Аи = д, и,д € Н. (7)

Из результатов, доказанных в [8], следует существование унитарного оператора Q тако-

1 /2

го, что имеет место полярное разложение оператора А и А = QA^ . Пусть Q* - оператор, сопряженный к Q и д, = Q*g,. Метод проекционной регуляризации, следуя [9], основан на том, что регуляризующее семейство € (0,5о], а € (0, ||А||]} определяется формулой:

1|А||

Т = J 1 д,, (8)

а

1/2

где {Ест| а € [0, ||А||]} - спектральное разложение Е, порожденное оператором А/ .

В качестве параметра регуляризации а выберем одно из решений уравнения:

!№% - д,||2 = 9||А||252. (9)

При использовании предложенного подхода к построению метода проекционной регуляризации приближенное решение уравнения (7) определяется формулой:

_ [ Т“д,, при ||д,|| > 3||А||5,

иа = Тад, = ^ (10)

I 0 , при ||д,|| < 3||А||5,

где оператор Т^д, задан формулой (8), а а - одно из решений уравнения (9).

В работе [9] доказано, что семейство операторов {Т,*| 5 € (0, 5о], а € (0, ||А||]} является регуляризующим семейством для уравнения (7), и погрешность регуляризованного решения иа, определяемого соотношением (10), является неулучшаемой по порядку.

Это свойство регуляризованного решения послужило основой для проведения сравнительного анализа решений, получаемых с использованием преобразований Лапласа и Фурье.

3. Численное моделирование решения обратной задачи для преобразования Лапласа

Если д(0) =0 и д(Ь) € ^([0; то)) П С 1([0, Т]) для любого Т = то, то, введя следующие обозначения:

N

Кп(* - Т) = ^(-1)га+1п2п2в-п2га2(*-т), (11)

П=1

N

Кп(0) ^(-1)П+1П2П2,

П=1

получаем, что для любого Т = то при Ь € [0; Т] уравнение (6) можно рассматривать как уравнение Вольтерра второго рода. При этом необходимо учитывать [10] положительную определенность ядра (11). Это возможно, если N - нечетное, тогда уравнение (6) примет вид

* N

А^> = -^(Ь) + 2 ]Т(-1)>(т)КП(Ь - т) йт = д,(Ь). (12)

0 га=1

Таким образом, задача (1),(2) сведена к решению уравнения (12). Из этого уравнения требуется найти функцию <^(Ь) при условии, что ||д(Ь) - д,(Ь)|| < 5.

Это решение будет зависеть от числа N, которое можно рассматривать как параметр регуляризации. Для отыскания этого параметра используем метод проекционной регуляризации, тогда количество слагаемых в частичной сумме ряда (12) определяется из уравнения:

||А^) - д,(Ь)||2 = 952. (13)

Предложенный подход послужил основой для создания следующей вычислительной схемы решения задачи (1), (2):

• Используя формулу правых прямоугольников, заменяем интегралы на интегральные суммы, введя равномерную сетку узлов на [0,Т].

• Находим значение N - количество слагаемых в ряде (12). С этой целью применяем итерационный процесс. Для обеспечения положительной определенности ядра в (12), необходимо рассматривать только нечетные значения N.

• Решаем полученную систему линейных уравнений.

• Итерационный процесс повторяем до тех пор, пока не будет выполнено условие (13).

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

4. Численное моделирование решения обратной задачи для преобразования Фурье

Если известно, что решение прямой задачи удовлетворяет условиям «(ж, Ь), м^(ж, Ь) € С([0,1] х [0, то)) П £2([0,1] х [0, то)),

и*(ж, Ь), ижж(ж, Ь) € ^2([0,1] х [0, то)),

и для функции и(Ь) € ^^([0, то)) выполнено ||«(Ь)|^1([0,^)) < Г то для решения задачи (1), (2) используем преобразование Фурье по переменной Ь, предварительно осуществив преобразование Фурье-образов функции д, с помощью метода проекционной регуляризации.

Пусть ||д,||ь2 > 35. Применяя преобразование Фурье Е* по переменной Ь для исходной задачи (1), (2), получаем, что операторное изображение этой задачи имеет вид:

Й2

гЛй(ж, Л) = ^“2 «(ж, А), ж € (0,1), Л € Л, (14)

«(0, Л) = 0, йж(0,Л)= д,(Л),

где и(ж, Л) = и(ж,Ь)), а д,(Л) = Е*(д,(Ь)) - Фурье образы функций и(ж,Ь) и д,(Ь) соответ-

ственно.

В силу того, что преобразование Фурье необходимо применить только для Ь > 0, то спроектируем Фурье-образ функции д, на Л(Е*). Тогда функции д, при Ь > 0 будет соответствовать функция д, (Л,а:(5)) такая, что:

( д,(Л), |Л| >а(5) д, (Л,а(5)) = < (15)

[ 0 , |Л| < а(5).

В качестве параметра регуляризации а(5), заданного формулой (9), выберем решение следующего уравнения:

—а(,)

У |д,(Л)|2^л |д,(Л)|2^л = 952. (16)

—^ а(,)

Так как функция д,(Л) была преобразована с помощью соотношений (15) и (16), то, решая задачу (14), получаем, что искомая функция и(1, Л) имеет вид:

«(1,Л) = д,(17)

^оУ Л

где ^о = (1 + *), параметр регуляризации а(5) - одно из решений уравнения (17), а

у2

функция д,(Л,а(5)) задана формулой (16).

Применим к функции и(1, Л) обратное преобразование Фурье Е-1. В результате получим следующую функцию:

«, (М) = Е— ^и, (1,Л)). (18)

Используя для функции (18) метод проекционной регуляризации, имеем, что регуляризо-ванное решение задачи (1), (2) определяется формулой:

( Ие й,(1,Ь) при Ь > 0, и,(Ь) = Р(И,е й,(1,Ь)) = < (19)

[ 0 , при Ь < 0.

Предложенная схема решения служит основой для разработки численного метода решения задачи (1), (2). Вычислительная схема метода имеет следующий вид:

• Применяем преобразования Фурье для элементов д,, таких, что ||д,||ь2 > 35, находим д,(Л).

• Приближенными методами решаем уравнение (16), находим параметр а(5).

• Согласно формуле (15) получаем д,(Л,а(5)).

• Используя формулу (17), находим Фурье-образ й(1,Л) регуляризованного решения.

• Применяя формулы (18) и (19), получаем регуляризованное решение.

Из результатов работы [9] следует, что погрешность полученных решений неулучшае-ма по порядку. Это свойство послужило основой для проведения сравнительного анализа численных методов, разработанных на основе преобразований Лапласа и Фурье.

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

5. Результаты вычислительного эксперимента

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

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

Вычислительный эксперимент проводился для Ь € [0, Т]. Была выбрана равномерная сетка из т + 1 узлов, Т = Лт, и ^ = (г — 1)Л,, г = 1, (т + 1). Основные этапы вычислительного эксперимента состоят в следующем:

• Для построения функции $(Ь), такой, что д(0) = 0 методом конечных разностей решаем прямую задачу (3). Используя метод численного дифференцирования, определяем значения функции д(Ь) в узлах сетки.

• Для функции д(Ь) моделируем значения (Ь^) в каждом узле сетки как значения нормально распределенной случайной величины с математическим ожиданием а = $(Ьг) и среднеквадратичным отклонением а = 5.

• Для решения обратной задачи (1), (2) по исходным данным (Ь^) и 5 > 0 реализуем вычислительные схемы, основанные на преобразовании Лапласа и Фурье.

• Вычисляем величины ||и<5 — ио||, где ио - модельная функция, и - регуляризованное решение, полученное при использовании соответствующего подхода.

Результаты проведенного эксперимента проиллюстрированы на нижеприведенных рисунках. На всех рисунках приняты общие обозначения. На нижнем графике обозначение $о(£) соответствует точному значению их(0,Ь) = $(Ь), а соответствует его возмущенному значению. На верхних графиках линия, соответствующая точному значению модельной функции и(1,Ь), обозначена как ио. Обозначение соответствует регуляризованному решению, полученному с использованием преобразования Фурье, а - регуляризованному решению, полученному с использованием преобразования Лапласа, 5 - величина погрешности, для которой проводились расчеты, а - величина параметра регуляризации для подхода, использующего преобразование Фурье.

Пример 1. В проведенной серии экспериментов рассматривались непрерывные монотонные функции. Эксперимент проводился при различных уровнях погрешности. На рис. 1 изображены результаты численного решения задачи (1), (2) для модельной функции и(1, Ь) = Ье-* при 5 = 0, 03.

их(0Д)=д(Ч

и(1,Ч=«Ч

Рис. 1. Результаты численного решения задачи (1),(2)

Пример 2. В этой серии экспериментов рассматривались непрерывные гладкие функции с одним экстремумом. Эксперимент проводился при различных уровнях погрешности. На рис. 2 изображены результаты численного решения задачи (1), (2) для модельной функции и(1, Ь) = Ь(е-* — е-1) при 5 = 0, 01.

их(0Д)=д(Ц

и(1Д)=фЩ

Рис. 2. Результаты численного решения задачи (1),(2)

Пример 3. В предложенной серии экспериментов исследовались непрерывные гладкие функции с несколькими экстремумами. Эксперимент проводился при различных уровнях погрешности. На рис. 3 изображены результаты численного решения задачи (1), (2) для модельной функции и(1,Ь) = 8ш(3пЬ)е-* при 5 = 0.1, а на рис. 4 изображены для модельной функции и(1,Ь) = 8т(10пЬ)е-* при 5 = 0,1.

их(0Д)=д(Ц

и(1Д)=фЩ

Рис. 3. Результаты численного решения задачи (1),(2)

их(0Д)=дЩ

Рис. 4. Результаты численного решения задачи (1),(2)

Пример 4. В проведенной серии экспериментов рассматривались кусочно-непрерывные функции. Эксперимент проводился при различных уровнях погрешности. На рис. 5 и рис. 6 изображены результаты численного решения задачи (1), (2) для модельных функций при 5 = 0,1.

их(0Д)=дЩ

0 0.2 0.4 0.6 0.8 1

и(1Д)=фЩ

Рис. 5. Результаты численного решения задачи (1),(2)

Для проведения сравнительного анализа численных методов необходимо было в процессе проведения вычислительного эксперимента получить величину уклонения регуляризован-ных решений от ио - заданных модельных функций. С этой целью для каждой функции осуществлялось несколько повторных расчетов, при этом каждый раз вычислялись величины ||и^ — ио| и ||ир — ио||, где - решение, построенное с использованием преобразования Фурье, а - с использованием преобразования Лапласа. Средние значения этих величин, полученных в каждой серии эксперимента, представлены в таблице.

их(0Д)=д(Ц

Рис. 6. Результаты численного решения задачи (1),(2)

Таблица

Модельная функция Значение Среднее Среднее

5 значение ||м£ - Мо|| значение ||и,Т - ио||

и = іе- 0,03 0,0435 0,0439

0,01 0,0149 0,0152

I—1 ) і 1 «-к 1 е 1 0,03 0,0409 0,0441

0,01 0,0144 0,0155

и(1,і) = 8Іп(3пі)в - 0,01 0,0520 0,0183

0,1 0,1478 0,1393

и(1,і) = 8Іп(10пі)в-* 0,01 0,1065 0,0180

0,1 0,1985 0,1875

( —40і, і Є [0; 0, 025], и (1 і) 1 -1, і є (°,025;05], и (1,і) ] 80і - 39, і Є (0, 5; 0, 525], 0,01 0,0194 0,0141

[ 1, і Є (0,525; 1].

0,1 0,2304 0,0164

(1 і) Г віп(5пі) — 2.5і2, і Є [0;0.5), и( ’ ) | еов(5пі), і Є [0.5; 1]. 0,01 0,1 0,0174 0,2004 0,0140 0,0151

На основании полученных результатов эксперимента можно сделать следующие выводы. Оба подхода позволяют находить регуляризованные решения с достаточной точностью. Подход, основанный на преобразовании Фурье и методе проекционной регуляризации, позволяет получать результаты с гарантированной точностью, что подтверждается теоретическими оценками метода проекционной регуляризации, приведенными в работе [9]. Однако, при увеличении порядка погрешности исходных данных изменяется и порядок погрешности полученных решений и численный метод решения задачи, основанный на преобразовании Фурье и методе проекционной регуляризации, более сложен, чем метод, основанный на преобразовании Лапласа.

Подход, основанный на преобразовании Лапласа и методе проекционной регуляризации, позволяет построить приближенные решения с достаточной степенью точности. С целью уточнения оценки погрешности решений, полученных при использовании этого подхода, и улучшения сходимости построенных решений в дальнейшем планируется рассмотреть численное решение уравнения такого типа с помощью квадратур правых (средних) прямоугольников по аналогии с подходом, предложенным в [10].

Заключение

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

Подход, основанный на преобразовании Лапласа, позволил свести решение исходной задачи к решению корректных интегральных операторных уравнений и исключить элементы операционного исчисления из процесса регуляризации. Численный метод, используемый при таком подходе, более прост с точки зрения компьютерной реализации, чем метод, использующий преобразования Фурье. Экспериментальные оценки погрешности полученных решений свидетельствуют о достаточной точности полученных решений. Следует заметить, что необходимо провести дополнительные исследования, связанные с уточнением оценок погрешности полученных решений.

Литература

1. Jonas, P. Approximate Inverse for a One-Dimensional Inverse Heat Conduction Problem /

P. Jonas, A.K. Louis // Inverse Problems. - 2000. - V. 16. - С. 175-185.

2. Prud’homme, M. Fourier Analysis of Conjugate Gradient Method applied to Inverse Heat Conduction Problems / M. Prud’homme, T.H. Hguyen // International J. of Heat and Mass Transfer. - 1999. - V. 42. - С. 4447-4460.

3. Cialkowski, M. Sequential and Global Method of Solving an Inverse Problem of Heat Conduction Equation / M. Cialkowski, K. Grysa // J. of Theoretical and applied Mechanics.

- 2010. - V. 48. - С. 111-134.

4. An Analytical Solution for Two-Dimensional Inverse Heat Conduction Problems Using Laplace Transform. International / M. Monde, H. Arima, W. Liu, Y. Mitutake, J.A. Hammad // J. of Heat and Mass Transfer. - 2003. - V. 46. - P. 2135-2148.

5. Деч, Г. Руководство к практическому применению преобразования Лапласа и Z-

преобразования / Г. Деч. - М.: Наука, 1971.

6. Япарова, Н.М. О различных подходах к решению обратных граничных задач тепловой диагностики / Н.М. Япарова // Вестник ЮУрГУ. Серия «Математика. Механика. Физика>. - 2012. - № 34. - С. 60-67.

7. Лаврентьев, М.М. Некорректные задачи математической физики и анализа / М.М. Лаврентьев, В.Г. Романов, С.П. Шишатский. - М.: Наука, 1980.

8. Менихес, Л.Д. Конечномерная аппроксимация в методе М.М. Лаврентьева / Л.Д. Ме-нихес, В.П. Танана // Сибирский журнал вычислительной математики. - 1998. - Т. 1, № 1. - С. 416-423.

9. Танана, В.П. Об оптимальном по порядку методе решения условно-корректных задач / В.П. Танана, Н.М. Япарова // Сибирский журнал вычислительной математики. - 2006.

- Т. 9, № 4. - С. 353-368.

10. Солодуша, С.В. Моделирование систем автоматического управления на основе полиномов Вольтерра / С.В. Солодуша // Моделирование и анализ информационных систем.

- 2012. - № 1. - С. 60-68.

Наталья Михайловна Япарова, кандидат физико-математических наук, кафедра «Прикладная математика>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), nataly7@mail.ru.

Bulletin of the South Ural State University. Series «Mathematical Modelling, Programming & Computer Software>,

2013, vol. 6, no. 3, pp. 112-124.

MSC 65N21 58J35; 35R30; 80M30

Numerical Simulation for Solving an Inverse Boundary Heat Conduction Problem

N.M. Yaparova, South Ural State University, Chelyabinsk, Russian Federation, nataly7@mail.ru

This paper proposes different approaches that help to find numerical solution to the boundary problem for heat equation. The Laplace and Fourier transforms are the basis for these approaches. The application of the Laplace transform allowed us to obtain an operator equation which connected the unknown function at one boundary with the initial data on the other boundary. The approach based on the Fourier transform for a time variable enables us to get a stable solution for the inverse problem of heat diagnostics. The obtained results are used for devising numerical methods. Comparative computational analysis of these approaches shows the limits of applications and effectiveness of each numerical method.

Keywords: boundary value problems for heat equation, regularization methods, the Laplace transform, Fourier transform, method of projective regularization.

References

1. Jonas P., Louis A.K. Approximate Inverse for a One-Dimensional Inverse Heat Conduction Problem. Inverse Problems, 2000, vol. 16, pp. 175-185.

2. Prud’homme M., Hguyen T.H. Fourier Analysis of Conjugate Gradient Method Applied to Inverse Heat Conduction Problems. International J. of Heat and Mass Transfer, 1999, vol. 42, pp. 4447-4460.

3. Kolodziej J., Mierzwiczak M., Cialkowski M. Application of the Method of Fundamental Solutions and Radial Basis Functions for Inverse Heat Source Problem in Case of Steady-State. International Communications in Heat and Mass Transfer, 2010, vol. 37, no. 2, pp. 121-124.

4. Monde M., Arima H., Liu W., Mitutake Y., Hammad J.A. An Analytical Solution for TwoDimensional Inverse Heat Conduction Problems Using Laplace Transform. International J. of Heat and Mass Transfer, 2003, vol. 46, pp. 2135-2148.

5. Doetsch G. Anleitung zum praktischen gebrauch der Laplace-transformation und der Z-transformation. R. Oldenbourg, 1961. 256 p.

6. Yaparova N.M. Different Approaches to Solve Inverse Bonduary Value Problems of Thermal Diagnostics. Bulletin of the South Ural State University. Series «.Mathematics. Mechanics. Physics>, 2012, vol. 34, pp. 60-67. (in Rusian)

7. Lavrentiev M.M., Romanov V.G., Shishatskii S.P. Nekorrektnye zadachi matematicheskoy fiziki i analiza [Ill -Posed Problems of Mathematical Physics and Analysis]. Moscow, Nauka, 1980. 287 p.

8. Menikhes L.D., Tanana V.P. The Finite Dimensional Approximation for the Lavrentiev Method. Sibirskiy zhurnal vychislitel'noy matematiki [Numerical Analysis and Applications], Novosibirsk, 1998, vol. 1, no.1, pp. 416-423. (in Russian)

9. Tanana V.P., Yaparova N.M. The Optimum in Order Method of Solving Conditionally-Correct Problems. Sibirskiy zhurnal vychislitel'noy matematiki [Numerical Analysis and Applications], Novosibirsk, 2006, vol. 9, issue 4, pp. 301-316. (in Russian)

10. Solodusha S.V. Automatic Control Systems Modeling by Volterra Polynomials. Modeling and Analysis of Information System, 2012, vol. 19, issue 1, pp. 60-68. (in Russian)

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