close

Вход

Забыли?

вход по аккаунту

?

Аналитические методы решения обратных задач теплопроводности для одномерных объектов.

код для вставкиСкачать
ISSN 1810-0198 Вестник ТГУ, т.17, вып.3, 2012
УДК 66.01.011
АНАЛИТИЧЕСКИЕ МЕТОДЫ РЕШЕНИЯ ОБРАТНЫХ ЗАДАЧ
ТЕПЛОПРОВОДНОСТИ ДЛЯ ОДНОМЕРНЫХ ОБЪЕКТОВ
© Е.Н. Туголуков, В.А. Карпук, А.В. Рухов
Ключевые слова: обратная задача теплопроводности; аналитическое решение; математическое моделирование.
Рассмотрена методика аналитического решения обратных задач теплопроводности для одномерных областей в
декартовых и цилиндрических координатах.
При математическом моделировании температурных полей, как правило, основным источником погрешностей служат значения конвективных коэффициентов теплоотдачи, входящие в граничные условия
третьего рода задачи теплопроводности.
Коэффициент теплоотдачи является комплексной
характеристикой интенсивности теплообмена теплоотдающей (тепловоспринимающей) поверхности и омывающего ее потока жидкости (газа). Он зависит от
большого количества физических, геометрических и
режимных параметров теплообменного процесса. Поэтому вывод прямых аналитических зависимостей для
расчета коэффициентов конвективной теплоотдачи на
основе фундаментальных знаний о природе процессов
теплопереноса в пространстве не представляется возможным.
Существуют различные возможности для определения численных значений коэффициентов конвективной теплоотдачи в конкретных условиях протекания
теплообменного процесса.
Классическая инженерная методика расчета коэффициентов конвективной теплоотдачи, базирующаяся
на теории подобия, основана на использовании критериальных уравнений алгебраического типа, обобщающих экспериментальные данные по различным веществам, выступающим в роли теплоносителей, для каждого набора условий протекания теплообменного процесса. Поэтому использование критериальных уравнений,
являющихся, по сути, результатом многомерной аппроксимации, приводит в каждом конкретном расчете к
погрешностям, не поддающимся оценке [1]. Как правило, погрешность расчета коэффициентов конвективной теплоотдачи по критериальным уравнениям составляет 30–50 %.
Другой путь связан с непосредственным измерением температурных полей в лабораторных или промышленных условиях на действующем оборудовании для
исследуемых условий протекания теплообменных процессов и видов теплоносителей. По результатам измерений температурных полей могут быть вычислены
локальные значения коэффициентов теплоотдачи.
Существует ряд методик расчета коэффициентов
теплоотдачи по экспериментальным данным:
– расчет непосредственно по определению коэффициента теплоотдачи как удельного количества тепла,
986
приходящегося на единицу площади поверхности теплообмена в единицу времени и отнесенного к единичной разности температур поверхности и определяющей
температуры потока;
– итеративный алгоритм нахождения коэффициента теплоотдачи, в котором при каждой итерации
решается прямая задача теплопроводности, т. е. рассчитывается температурное поле в моделируемых условиях и корректируется значение коэффициента теплоотдачи, входящего в граничные условия задачи теплопроводности; итерации выполняются до приемлемого совпадения расчетного и измеренного температурных полей;
– прямой расчет коэффициента теплоотдачи по
результатам решения обратной задачи теплопроводности для исследуемого процесса.
Первый способ позволяет определить усредненные
значения коэффициентов теплоотдачи по результатам
экспериментов, выполненных на экспериментальных
установках, с помощью которых можно определять
значения тепловых потоков. Этот способ широко освещен в литературе и является своего рода классическим, хотя применим только для стационарных процессов теплообмена.
Второй способ является «кустарным», т. е. не имеет
строгого аналитического обоснования, хотя он относительно прост и нагляден. Этот способ целесообразно
использовать в оценочных расчетах.
Третий способ математически строг, но сложен и
специфичен. Решению обратных задач теплопроводности посвящено много работ, но их результаты часто
оказываются не адаптированными для решения прикладных инженерных задач.
Рассмотрим возможные относительно простые пути
определения коэффициента теплоотдачи по результатам экспериментальных исследований с использованием решения обратных задач теплопроводности.
Для решения обратных задач теплопроводности
преимущественно используются численные методы [1–
2], аналитическим подходам уделяется меньшее внимание. Это связано со спецификой аналитических методов и возможностью их использования лишь для
ряда частных случаев процессов теплообмена.
В [3] приводится методика, основанная на использовании интегральных преобразований Лапласа для
ISSN 1810-0198 Вестник ТГУ, т.17, вып.3, 2012
решения обратной задачи теплопроводности в цилиндрических координатах.
В качестве иллюстрации такого подхода (для декартовых координат) рассмотрим процесс конвективного теплообмена плоской неограниченной вертикальной пластины с окружающей средой, имеющей постоянную температуру.
Пусть измеряется температура во времени в какойто точке х0 по толщине пластины. Температурное поле
такой пластины описывается решением следующей
задачи теплопроводности:
∂ t ( x, τ )
∂ 2 t ( x, τ )
=a
,
∂τ
∂ x2
0 ≤ x ≤ R,
τ ≥ 0;
(1)
с начальным условием
t (x, 0) =t0 = сonst ,
(2)
a
λ
d 2 S ( x, p )
− p S ( x, p ) = 0 ;
(10)
⎛
t +t ⎞
∂ S (R , p )
+ α ⎜⎜ S (R, p ) − c 0 ⎟⎟ = 0 ;
p ⎠
∂x
⎝
(11)
d x2
∂ S (0, p )
=0 .
∂x
Решение этой задачи в изображениях имеет вид:
⎛ p ⎞
⎛ p ⎞
S (x, p ) = A sh ⎜
x ⎟ + B ch ⎜
x⎟ ,
⎜ a ⎟
⎜ a ⎟
⎠
⎝
⎠
⎝
∂ t ( R, τ )
+ α (t (R, τ) − tc )= 0 ;
∂x
∂ t (0, τ)
=0 ,
∂x
p
a
λ
(3)
(4)
где х – координата, направленная по нормали к поверхности пластины; τ – время; t(x, τ) – температурное
λ
– коэффициент температуроcρ
проводности материала пластины; λ,с, ρ – теплопро-
p
a
λ
водность, теплоемкость и плотность материала пластины соответственно; t0 – начальная температура пластины; tс – температура окружающей среды; R – полутолщина пластины; α – коэффициент теплоотдачи.
Путем замены переменной T (x, τ ) = t (x, τ ) − t0
можно получить задачу с нулевым начальным распределением:
откуда
λ
∂ T (R, τ)
+ α (T (R, τ) − tc + t0 )= 0 ;
∂x
∂ T (0, τ)
=0 .
∂x
(5)
⎛
⎛
⎞
⎛
⎞⎞
⎜ A ch ⎜ p R ⎟ + B sh ⎜ p R ⎟ ⎟ +
⎜ a ⎟
⎜ a ⎟⎟
⎜
⎝
⎠
⎝
⎠⎠
⎝
⎛
⎛ p ⎞
⎛ p ⎞⎞
R ⎟ + B ch ⎜
R⎟⎟ −
+ α ⎜ A sh ⎜
⎜ a ⎟
⎜ a ⎟⎟
⎜
⎝
⎠
⎝
⎠⎠
⎝
t c − t0
−α
= 0;
p
поле пластины; a =
∂ T ( x, τ )
∂ 2 T ( x, τ )
=a
;
∂τ
∂ x2
T (x, 0)= 0 ;
(8)
(14)
⎛
⎛
⎞
⎛
⎞⎞
⎜ A ch ⎜ p 0 ⎟ + B sh ⎜ p 0 ⎟ ⎟ = 0 ,
⎜ a ⎟
⎜ a ⎟⎟
⎜
⎝
⎠⎠
⎝
⎠
⎝
(15)
А = 0,
B=
(16)
α (tc − t0 )
⎛
p⎜λ
⎜
⎝
⎛ p ⎞⎞
⎛ p ⎞
p
sh ⎜
R⎟ ⎟
R ⎟ + α ch ⎜
⎟
⎜ a ⎟⎟
⎜
a
a
⎠⎠
⎠
⎝
⎝
.
(17)
Тогда
(6)
(7)
(13)
где sh ( z ), ch (z ) – гиперболические функции.
A и B определяются из граничных условий (11),
(12):
и граничными условиями
λ
(12)
S ( x, p ) =
⎛
p⎜λ
⎜
⎝
⎛ p ⎞
α (tc − t0 ) ch ⎜
x⎟
⎜ a ⎟
⎠
⎝
. (18)
⎛ p ⎞⎞
p ⎛⎜ p ⎞⎟
⎟
sh
R + α ch ⎜
R⎟
⎜ a ⎟⎟
⎜ a ⎟
a
⎠⎠
⎠
⎝
⎝
Если известно изображение температурного поля
Выполним преобразование Лапласа данной задачи
по формуле:
S (x , p ) =
∞
∫
T (x , τ ) еxp (− p τ ) d τ .
0
Получаем задачу в изображениях:
(9)
S (x0 , p ) , значение коэффициента теплоотдачи может
быть найдено аналитически:
⎛ p ⎞
p
R⎟
sh ⎜
⎜ a ⎟
a
⎝
⎠
. (19)
α=
⎛ p ⎞
⎛ p ⎞
(tc − t0 )ch ⎜⎜ x0 ⎟⎟ − p S (x0 , p)ch ⎜⎜ R ⎟⎟
⎝ a
⎠
⎝ a ⎠
S (x0 , p ) pλ
987
ISSN 1810-0198 Вестник ТГУ, т.17, вып.3, 2012
Изображение температурного поля S (x0 , p ) определяется по формуле:
S (x0 , p ) =
∞
∫
где функция S(х), являющаяся ядром интегрального
преобразования, определяется как решение вспомогательной задачи Штурма–Лиувилля:
d 2 S (x )
(ti (x0 , τi ) − t0 ) exp(− p τi ) dτi =
d x2
0
τk
=
∫
(ti (x0 , τi ) − t0 ) exp(− p τi ) dτi +
(20)
(27)
d S (R )
+ α S (R ) = 0 ;
dx
(28)
d S (0 )
=0 .
dx
0
t ( x , τ ) − t0
+ k 0 k
exp(− p τ k ).
p
Здесь ti (x0 , τ ) , i = 1, 2, ..., k – экспериментальные значения температур, измеренных в фиксированной точке
х0 в различные моменты времени.
Можно показать, что расчетный результат не зависит от численного значения параметра интегрального
преобразования Лапласа «р», которое может быть как
действительным, так и комплексным.
Однако практически метод может быть использован только для регулярного температурного режима изза его ярко выраженной неустойчивости по отношению
к погрешности измерений температур при нерегулярных температурных режимах.
Другой метод решения обратной задачи теплопроводности основан на использовании конечных интегральных преобразований [4]. Он устойчив, но требует
экспериментального определения температурного профиля по толщине образца хотя бы в единственный момент времени. Определение пространственного температурного профиля может быть сложно осуществимо
практически, особенно для действующего промышленного оборудования.
Проиллюстрируем этот метод на том же примере,
для пластины, температурное поле которой описывается решением задачи (1)–(4).
Введем другую замену переменной:
T ( x, τ) = t (x, τ) − tс ,
λ
+ μ 2 S (x ) = 0 ;
(21)
(29)
Решение этой задачи с точностью до постоянного
множителя имеет вид:
S ( x ) = sin (μ x + ϕ ) .
(30)
Тогда
d S (x )
= μ cos (μ x + ϕ) .
dx
(31)
Из граничного условия (29) μ cos ( ϕ ) = 0 ,
откуда
ϕ=
π
(минимальный положительный корень),
2
(32)
(33)
следовательно,
S (x ) = cos (μ x ) .
(34)
Из граничного условия (28) имеем:
− λ μ sin (μ R ) + α cos (μ R ) = 0 ,
(35)
откуда в конечном итоге и определяется искомый коэффициент теплоотдачи:
позволяющую перейти к задаче с однородными граничными условиями:
α = λ μ tg (μ R ) .
∂ T ( x, τ ) ∂ 2 T ( x, τ )
=a
;
∂τ
∂ x2
(22)
Весовая функция ρ определяется как решение уравнения:
(23)
a ρ' = 0,
(37)
откуда ρ = сonst, в частности ρ = 1.
Выполняем переход к изображениям:
(38)
T ( x, 0 ) = Т 0 = t0 − tc ;
∂ T ( R, τ )
λ
+ α T ( R, τ ) = 0 ;
∂x
∂ T (0, τ)
=0 .
∂x
(24)
(25)
Применим к этой задаче конечное интегральное
преобразование:
W (τ) =
0
988
d W (τ)
+ a μ2 W (τ) = 0 ;
dτ
(26)
(39)
изображение начального условия (23):
R
∫ T (x, τ)ρ S (x )dx ,
(36)
W (0) =
R
T0
∫ T cos (μ x) dx = μ
0
0
sin (μ R ) .
(40)
ISSN 1810-0198 Вестник ТГУ, т.17, вып.3, 2012
Решение задачи (22)–(23) в изображениях:
(
)
W (τ) = W (0) exp − μ 2 a τ =
T
= 0 sin (μ R ) exp − μ 2 a τ .
μ
(
)
(41)
Изображение температурного профиля находим по
формуле:
W (τ) =
R
∫ T cos (μ x )dx ,
(42)
i
0
где Ti = T (xi , τ ) – массив экспериментальных значений температур; i = 1, 2, …, N – номер точки измерения
температуры.
От точности вычисления этого интеграла, естественно, зависит погрешность конечного результата. При
использовании численной схемы интегрирования не
ниже третьего порядка точности погрешность расчета
коэффициента теплоотдачи практически определяется
погрешностью измерения температур.
Значение μ определяется как любой положительный корень уравнения:
I=
(
)
T0
sin (μ R ) exp − μ 2 a τ .
μ
(43)
⎛ ∂ 2 t ( r , τ) 1 ∂t (r , τ) ⎞
∂t (r , τ)
⎟;
= а ⋅ ⎜⎜
+ ⋅
2
∂τ
r
∂ r ⎟⎠
⎝ ∂r
(44)
t ( r , 0 ) = t 0 ( r );
(45)
∂t ( R , τ)
+ α ⋅ (t ( R , τ ) − t c ) = 0 ;
∂r
∂t (0, τ)
= 0.
∂r
(46)
(47)
Тогда:
T ( r ,0 ) = t 0 ( r ) − t c ;
(50)
λ⋅
∂T ( R , τ)
+ α ⋅ T ( R , τ ) = 0;
∂r
∂T (0, τ)
= 0.
∂r
(48)
(51)
(52)
Для исключения координаты r используется формула перехода к изображениям:
W (μ , τ) =
R
∫ ρ ( r ) ⋅ T ( r , τ ) ⋅ S ( r , μ ) dr ,
0
(53)
где ρ(r) = r – весовая функция, являющаяся решением
уравнения:
dρ(r ) 1
− ⋅ ρ(r ) = 0.
dr
r
(54)
Ядро интегрального преобразования S(r, μ) – решение уравнения Штурма–Лиувилля с однородными граничными условиями (здесь μ – параметр):
dr
2
+
1 dS ( r , μ )
⋅
+ μ 2 ⋅ S ( r , μ ) = 0;
r
dr
dS ( R , μ )
+ α ⋅ S ( R , μ ) = 0;
dr
dS ( 0 , μ )
= 0.
dr
λ⋅
(55)
(56)
(57)
Решение уравнения (55), определяемое с точностью
до постоянного множителя, имеет вид:
S ( r , μ ) = J 0 (μ ⋅ r ) + C ⋅ Y0 (μ ⋅ r ),
(58)
где J0(z) и Y0(z) – функции Бесселя нулевого порядка
первого и второго рода соответственно.
Из (58) получаем:
dS ( 0 , μ )
= − μ ( J 1 ( 0 ) + C ⋅ Y1 ( 0 ) ) = 0 ,
dr
(59)
откуда С2 = 0.
Таким образом, решение задачи (55)–(57) имеет
вид:
S ( r , μ ) = J 0 (μ ⋅ r ).
Переход к задаче с однородными граничными условиями достигается заменой:
T ( r , τ ) = t ( r, τ ) − t c .
(49)
d 2 S (r , μ)
Лучше выбирать первый корень во избежание накопления погрешности компьютерного счета.
Теперь для момента времени τ из уравнения (36)
можно найти искомое значение коэффициента теплоотдачи α.
Устойчивость метода обусловлена сглаживанием
значений экспериментальных температур при интегрировании по толщине пластины.
Аналогично может быть получено решение обратной задачи теплопроводности для одномерного объекта в цилиндрических координатах.
Постановка задачи нестационарной теплопроводности для сплошного неограниченного цилиндра имеет
вид:
λ⋅
⎛ ∂ 2T ( r , τ ) 1 ∂ T ( r , τ ) ⎞
∂T (r , τ)
⎟;
= a ⋅⎜
+ ⋅
⎟
⎜
∂r
∂τ
r
∂r 2
⎠
⎝
(60)
Выполним переход к изображениям задачи (49)–
(52):
dW ( μ , τ )
+ μ 2 a ⋅ W ( μ , τ ) = 0;
dτ
(61)
989
ISSN 1810-0198 Вестник ТГУ, т.17, вып.3, 2012
R
W (μ ,0 ) =
∫
r ⋅ ( t 0 ( r ) − t c ) ⋅ J 0 (μ ⋅ r )dr .
(62)
α = μλ⋅
0
Решение задачи (61)–(62) в изображениях имеет
вид:
(
)
W ( μ , τ ) = W ( μ ,0 ) ⋅ exp − μ 2 a τ .
(63)
1.
2.
∫
0
(
R
(64)
)∫
= exp − μ 2 a τ ⋅ r ⋅ (t 0 ( r ) − t c ) ⋅ J 0 (μ ⋅ r ) dr .
0
Искомый коэффициент теплоотдачи α определяется
из граничного условия (56):
990
4.
Кабанихин С.И. Обратные и некорректные задачи. Новосибирск:
Сибирское научное издательство, 2009. 458 с.
Алифанов О.В. Обратные задачи теплообмена. М.: Машиностроение, 1988. 280 с.
Пономарев С.В., Мищенко С.В., Дивин А.Г. Метод, устройство и
автоматизированная система для исследования зависимости теплофизических свойств жидкости от скорости сдвига // Вестник
ТГТУ. Тамбов, 1995. Т. 1. № 1–2. С. 39-52.
Туголуков Е.Н. Математическое моделирование технологического
оборудования многоассортиментных химических производств:
монография. М.: Изд-во «Машиностроение», 2004. 100 с.
Поступила в редакцию 3 апреля 2012 г.
R
r ⋅ (t ( ri , τ ) − t c ) ⋅ J 0 (μ ⋅ r ) dr =
(65)
ЛИТЕРАТУРА
3.
Значение μ определяется как любой положительный корень уравнения, записанного для изображения
экспериментально определенного температурного поля
t (ri , τ) :
J 1 (μ ⋅ r )
.
J 0 (μ ⋅ r )
Tugolukov E.N., Karpuk V.A., Rukhov A.V. ANALYTICAL
SOLUTIONS OF REVERSE TASKS OF HEAT CONDUCTIVITY FOR ONE-DIMENSIONAL OBJECTS
Methodology of analytical solutions of reverse tasks of heat
conductivity is considered for one-dimensional areas in Cartesian
and cylindrical coordinates.
Key words: inverse problem of heat conductivity; analytical
solution; mathematical modeling.
Документ
Категория
Без категории
Просмотров
7
Размер файла
573 Кб
Теги
обратный, решение, аналитическая, метод, объектов, теплопроводность, одномерных, задачи
1/--страниц
Пожаловаться на содержимое документа