close

Вход

Забыли?

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

?

Представление уравнения теплопроводности в интегральной форме для задач измерения теплофизических величин..pdf

код для вставкиСкачать
Управление, вычислительная техника и информационные технологии
УДК 681.2:536.083
ПРЕДСТАВЛЕНИЕ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ
В ИНТЕГРАЛЬНОЙ ФОРМЕ ДЛЯ ЗАДАЧ ИЗМЕРЕНИЯ
ТЕПЛОФИЗИЧЕСКИХ ВЕЛИЧИН
Ю.И. Азима
Рассмотрена новая интегральная форма уравнения теплопроводности для
различных одно- и многомерных моделей объектов, применяемых в качестве образцов
для измерения теплофизических свойств. Представлены результаты исследований на
адекватность полученных уравнений. Показано применение данных уравнений для задач измерения тепловых величин.
Ключевые слова: интегральная форма, уравнение теплопроводности, теплофизические свойства, метод измерения, дискретная модель.
Потребность в измерении теплофизических свойств (ТФС) материалов существует во многих областях науки и промышленности: в частности, это строительство, теплоэнергетика, материаловедение, электроника.
Особенностью теплофизических измерений является большое разнообразие методов и, как правило, сравнительно невысокая точность, низкая производительность и сложная методика выполнения измерения. Это
объясняется сложностью измерительной задачи, поскольку объектом исследования является динамическая система с распределенными параметрами.
Теоретическую основу большинства современных методов измерения теплофизических свойств составляют аналитические решения краевых
задач теплопроводности во временной области или в области комплексного переменного. Данные методы обеспечивают измерение тепловых величин во многих сферах их применения. Тем не менее, они имеют ряд недостатков, во многом обусловленных общим подходом к их проектированию:
сравнительно невысокая точность, низкая производительность измерений,
сложное аппаратурное оформление и трудоемкая методика выполнения
измерения.
Многие из указанных недостатков могут быть устранены, используя
подход к их проектированию, основанный на обеспечении адекватности
интегрального закона сохранения тепла. Его математическое описание
представляет собой одну из интегральных форм уравнения теплопроводности [1]. Для одномерного объекта толщиной L это баланс тепла за некоторый интервал времени [τ1, τ2]:
τ2
L
∫ [q1 ( τ ) − q2 ( τ )]dτ = C ∫ t ( x, τ )dx τ
τ2
1
τ1
0
123
.
(1)
Известия ТулГУ. Технические науки. 2014. Вып. 2
где q1(τ), q2(τ)– тепловые потоки на противоположных гранях образца; С–
объемная теплоемкость.
Исходя из адекватности данного уравнения, из него можно получить расчетную формулу для определения объемной теплоемкости, которая используется в методе измерения данной величины [2]:
τy
[
]
C = 2 ∫ [q1 ( τ ) − q2 ( τ )]dτ L t1 (τ у ) + t2 (τ у ) ,
0
где t1(τ), t2(τ)– температура на противоположных гранях образца; τу– время, начиная с которого обеспечивается достаточная точность определения
средней на интервале [0, L] температуры по двум точкам.
Достоинством данного подхода является отсутствие необходимости
определения температурного поля модели объекта путем решения краевой
задачи, что автоматически снимает все проблемы математического и технического характера, связанные с этим. Вместе с тем, ограниченность интегральной формы (1), связанная с ее структурой и малым числом коэффициентов,
определяющих
тепловые
свойства,
ограничивает
ее
применимость.
В связи с этим возникает потребность поиска более информативного интегрального представления уравнения теплопроводности. Такое
уравнение можно получить из интегральной формы (1), записанной для
интервала с переменной правой границей, путем интегрирования ее по
данной переменной в пределах заданных границ. В результате такого преобразования для одномерного объекта в прямоугольной системе координат
получим интегральную форму следующего вида:
τ
L x
LQ (0, τ ) = λ ∫ [t (0, τ ) − t (L, τ )] dτ + C ∫ ∫ t ( x, τ )dxdx 0τ ,
τ
0
0
(2)
0 0
где Q (0, τ ) – количество тепла, поступившего в одномерное тело через
τ
0
границу с координатой х=0 за интервал времени [0, τ]; t(0, τ), t(L, τ) – температура, измеренная в точках х=0 и х=L.
Данную модель можно рассматривать как баланс средних на отрезке [0, L] количеств тепла: поступившего через границу х=0, прошедшего
через данный отрезок вследствие теплопроводности и аккумулированного
на нем за счет теплоемкости. В дальнейшем будем ее называть второй интегральной формой уравнения теплопроводности (ИФУТ-2).
Графическая иллюстрация получения данного уравнения показана
на рис.1. Повторный интеграл по координате х может быть интерпретирован как средневзвешенное значение средних температур на отрезках [0,
k⋅∆x] (k =1,…m→∞), умноженное на половину квадрата длины L:
L
x
m
m
∫ dx ∫ t( x )dx =0,5L k∑=1 k ⋅ tk / k∑=1 k ,
2
0
0
124
где t k = k
−1
k
∑ ti .
i =1
Управление, вычислительная техника и информационные технологии
Разность температур является величиной пропорциональной средней на отрезке [0, L] плотности теплового потока:
L
λ
t (L) − t (0 ) = ∫ t ′x ( x )dx = q ( L ) .
L
0
Рис. 1. Иллюстрация к ИФУТ-2:
ИФУТ q(0) – тепловой поток, действующий
на границе х=0 исследуемого участка [0, L] одномерного объекта
в направлении оси Ох;
Ох qλ(xi) – тепловой поток, проходящий через
сечения с координатой хi; qc(xi) – тепловой поток, аккумулированный
на i-ом отрезке [0, i· ∆x], где ∆х– толщина элементарного слоя
Можно показать, что если представить температурную зависимость
теплопроводности и объемной теплоемкости в виде линейной функции:
λ(t)=λ0+ λ1(t−t0), С(t)=С0+ С1(t−t0), , где t0 – начальная температура образца,
то уравнение (2) примет следующий вид:
τ
τ
0
0
[
]
LQ(0, τ ) = λ 0 ∫ [θ(0, τ ) − θ(L, τ )] dτ + 0.5λ1 ∫ θ 2 (0, τ ) − θ 2 (L, τ ) dτ +
τ
0
L x
L x
0 0
0 0
+ C0 ∫ ∫ θ( x, τ )dxdx 0τ + 0.5C1 ∫ ∫ θ 2 ( x, τ )dxdx 0τ ,
где θ = t– t0.
Для цилиндрических и сферических координат ИФУТ
ИФУТ-2 получим
аналогичным образом: запишем уравнение баланса количества тепла, приведенного к единичной площади
площади, для элементарного цилиндрического или
сферического слоя, и усредним его составляющие на заданном отрезке [r0,
R] по координате r. Получим:
Получим
Lκ Q(0, τ)
τ
0
τ
L
r
dr
= λ ∫ [t (r0 , τ) − t (R, τ)]dτ + C ∫ κ ∫ r κ t (r, τ)dr 0τ ,
0
0 r r0
125
Известия ТулГУ. Технические науки. 2014. Вып. 2
где Lκ =
r0κ
R
∫r
-κ
dr – приведенная толщина, соответственно цилиндриче-
r0
(
)
ской (κ=1) и сферической стенки (κ=2): L1 = r0 ln Rr0−1 , L2 = r02 r0−1 − R −1 ;
Q ( r0 , τ ) – количество тепла, прошедшего за время τ через сечение с координатой r0 .
Аналогичный вид имеют ИФУТ-2 для линейного (κ=1) и точечного
(κ=2) источников:
rm
r
τ
dr κ
∗
τ
Lκ Q ( τ) 0 = λ ∫ [t(r0 , τ) − t( R, τ)] dτ + C ∫ κ ∫ r t(r, τ)d r 0τ ,
r0 r 0
0
R
L
1
= κ κ = ∫ r - κ dr ; Q(τ) – количество тепла, выделяемого, соотгде
κπr0 κπ r0
ветственно в линейном, на единицу длины, и точечном источнике за время τ.
Покажем ИФУТ-2 для некоторых многомерных моделей.
L∗κ
1. Прямоугольник с высотой Lz образован прямыми: z1=0, z2=Lz, x1=0,
x2=Lx.
ИФУТ по координате х имеет вид:
( Lz )
Qx
(0, τ ) + 1
S xz
=
( Lz )
где Q x
2 Lx
x
j =1 0
0
∑
λ
Lx
τ
∫ dx ∫ Q z (x, z j , τ )dx =
∫ [t
0
( Lz )
(0, τ ) = − λ ∫ t x
′ ( Lz )
]
( 0, τ ) − t ( Lz ) ( Lx , τ ) dτ +
C
Lx
Lx
(3)
x
∫ dx ∫ t
0
( Lz )
( x , τ ) dx ,
0
(0 , τ )d τ – среднее по высоте Lz прямоугольника ко-
τ
личество
тепла,
поступившего
через
границу
х=0
за
время
τ;
τ
Q z (x, z j , τ ) = − λ ∫ t z′ (x, z j , τ )dτ – распределение по координате х количества
0
тепла, протекающего в направление координаты z через сечения: z1=0 и
z2= Lz;
t
( Lz )
( xi , τ) =
L−z1
Lz
∫ t( xi , z, τ)dz –
средняя по высоте Lz температура
0
границ прямоугольника с координатами x1 = 0, x2 = Lx .
Данное уравнение является частным случаем ИФУТ трехмерного
параллелепипеда при отсутствии потока тепла по координате у.
2. Круглый прямой цилиндр с радиусом основания R и высотой Lz.
ИФУТ по координате z:
126
Управление, вычислительная техника и информационные технологии
( Sr )
Qz
(0, τ) + 2
RLz
Lz
z
0
0
∫ dz∫ Q (R, z, τ)d z =
r
C
λ
( Sr )
( Sr )
=
[
t
(
0
,
τ
)
−
t
(
L
,
τ
)]
d
τ
+
z
Lz τ∫
Lz
Lz
z
∫ dz∫ t
0
( Sr )
(4)
( z, τ)dz,
0
− 2λ
dτ ∫ rt z′ (r ,0, τ )dr = − λ ∫ t z′(Sr ) (zi , τ )dτ – количество теп2 ∫
R τ 0
τ
ла, поступившего в цилиндр через основание z=0 цилиндра;
R
2
(Sr )
средняя
температура
основания
z=0;
t (0 ) = 2 ∫ rt (r , zi )dr –
R 0
где
Q z(Sr ) (0, τ ) =
R
Qr (R, z, τ) = −λ ∫ tr′ (R, z, τ)dτ – количество тепла поступившего через бокоτ
вую поверхность цилиндра; zi = 0, Lz ; Sr– площадь основания цилиндра.
Аналогичным образом были получены ИФУТ-2 для других моделей, в частности: полого цилиндра, трапеции, усеченного конуса.
Для проверки адекватности интегральных форм (3) и (4) использовались дискретные модели, представленные в среде Mathcad: прямоугольника размером 30h×10h, где h =0,2мм – дискретность по координате, и цилиндра, размером n×m, где n =10h – высота цилиндра; m=20h – радиус
основания цилиндра (h =0,2мм).
Температурные поля данных моделей были вычислены при идентичных краевых условиях с использованием функции Radau, использующей алгоритм RADAUS, для решения жестких систем дифференциальных
уравнений. Для прямоугольника они имеют следующий вид:
θτ
x ∈ [0, khx ] (k = 1, 5),
q( x, 0, τ ) =
t( x, z, 0) = 0 ,
,
− α t ( x, 0, τ ) x ∈ hx , Lx
(
)
(
)
[
]
q (0, z , τ ) = −α t (0, z , τ ) , q Lx , z, τ = −α t Lx , z, τ , q( x, Lz , τ) = −α t ( x, Lz , τ ) .
На рис.2 представлены результаты исследований данных моделей.
В вычислительном эксперименте были использованы следующие данные:
скорость изменения теплового потока принималась равной θ=100 Вт/(м2·с),
коэффициент
теплообмена
α=50
Вт/(м2·К),
теплопроводность
-7
2
λ=0,5 Вт/(м·К), температуропроводность а=5·10 м /с; для цилиндра:
λ=0,2 Вт/(м·К), а=1,2·10-7 м2/с, αr=20 Вт/(м2·К), αz=100 Вт/(м2·К),
Как видно из графиков, суммарная погрешность адекватности данных уравнений теплопроводности и дискретной модели на временном участке от 20 до 100 с не превышает 0,2 %. Увеличение погрешности на интервале [4, 20] с до 0,8 % связано с погрешностью, обусловленной шагом
дискретизации по координате.
127
Известия ТулГУ
ТулГУ. Технические науки. 2014. Вып. 2
Рис.2. Изменение во времени суммарной погрешности δ (
)
дискретной модели и адекватности ИФУТ для двухмерных моделей:
а – прямоугольника (координата х): б – цилиндра (координата
координата z);
составляющие ИФУТ,
ИФУТ отнесенные к среднему на интервале [0, L]
количеству тепла от источника Qис( L ) , где L=Lx – для прямоугольника
и L=R
R – для цилиндра, обусловленные:
теплопроводностью – Qλ(L ) / Qис(L ) (
),
теплоемкостью – QC( L ) / Qис(L ) ( ),
конвективным теплообменом – Qα( L ) / Q ис(L ) (
)
3. ИФУТ для полуограниченного тела с плоским источником конечных
размеров с учетом теплообмена свободной поверхности получена на базе
ИФУТ точечного источника.
источника
Для этого поверхность полуограниченного тела представим в виде
приемника тепла от источника с плотностью qis ( x, y, z = 0, τ) , и источника
тепла с плотностью qs ( x, y , z = 0, τ ) . Тогда можно считать, что данная поверхность подвергается тепловому воздействию от источника площадью
S = Sis + S s , где Sis, Ss – площадь, соответственно источника тепла и стоков.
Для получения
олучения математической модели в интегральной форме воспользуемся принципом суперпозиции температурных полей
полей, вызываемых
действием точечных источников мощностью qxy(τ), действующим на поверхности (z=0) в точке c координатой (х, у). На поверхности образца выберем две точки с координатами (х0; у0) и (хm; уm), лежащими на поверхностях полусфер с радиусами
радиусами:
rm ( x , y ) = ( x m − x ) 2 + ( y m − y ) 2 , r0 ( x , y ) = ( x0 − x ) 2 + ( y 0 − y ) 2 .
и приведенной толщиной полусферической стенки L∗2 ( x, y ) .
128
Управление, вычислительная техника и информационные технологии
Суммируя действие всех точечных источников тепла в выбранных
точках (х0, у0) и (хm, уm) с учетом независимости от координат плотность
источника тепла, получим ИФУТ в следующем виде:
τ
K is Qis ( τ ) − ∫∫ Lr ( x , y ) ⋅ Qs ( x , y , τ )d x dy = λ ∫ [t1 (τ ) − t2 (τ )] dτ + C ⋅ Ip S ( τ ) ,
Ss
(5)
0
где Qis(τ) – количество тепла, поступившего в образец на единицу площади за интервал времени [0, τ] от источника тепла; Qs ( x, y, τ) – функция распределения количества тепла, теряемого с единицы поверхности;
t1 ( τ) = t ( x0 , y0 , τ ) ; t2 ( τ) = t ( xm , ym , τ ) ; K is = ∫∫ L∗2 ( x, y )d x dy − коэффициент
S is
влияния источника тепла, площадью Sis на температуру в точках (х0, у0),
(хm, уm);
1
Ip S ( τ ) =
S
rm
r
dr 2
dS
∫∫S r∫ r 2 ∫0 r t xy ( r , τ) dr .
0
Исследование ИФУТ для плоского источника конечных размеров
проводилось для двух случаев: при отсутствии и наличии теплообмена поверхности полуограниченного тела с окружающей средой. Для определе)
ния адекватности данного уравнения при Qs ( τ ) = 0 использовалось точное
решение уравнения теплопроводности для точечного источника постоян-
[
(
)]
ной тепловой мощности q [3]: t ( r , τ ) = q (2 π λ r ) 1 - erf r / 2 a τ .
Для квадратного нагревателя погрешность адекватности уравнения
(5) на заданном временном интервале, а также его составляющие, отнесенные к количеству тепла от источника, показана на графике рис. 3,а.
С учетом конвективного теплообмена свободной поверхности полуограниченного тела исследование погрешности адекватности возможно
только для приближенной ИФУТ, в которой интегральный параметр
I~
pS ( τ) вычисляется по приближенной формуле, коэффициенты которой не
имеют аналитического выражения (в сравнении с ранее рассмотренными),
а могут быть определены только по методу наименьших квадратов для заданного временного интервала и диапазонов тепло- и температуропроводности.
В качестве дискретной модели полуограниченного тела использовался двумерный цилиндр размером 40h×40h, где h = 0,3мм. Тепловой поток источника действовал на поверхности, имеющую форму круга с радиусом: ris=2h, 5h, 8h, 10h и описывался следующими функциями времени:
qis(τ) = wτ, qis(τ)=q. Свободная от источника тепла поверхность имеет конвективный теплообмен с коэффициентом α =10 Вт/(м2·К).
Вычислительный эксперимент проводился для интервала времени
[20, 120] с, координатах (х0 = 0,5h; у0=0; z= 0) и (хm = 20h; 0; 0) (h= 0,3 мм)
−1
129
Известия ТулГУ. Технические науки. 2014. Вып. 2
точек определения разности температур и следующих значений ТФС: λ=
(0,03…0,6) Вт/(м·К), а = (0,8…5)·10-7 м2/с . Вычисление приближенного
pS ( τ) осуществлялось по температурам, определяемым в трех
значения I~
точках.
а
б
Рис. 3. Результаты исследования на адекватность ИФУТ для
полуограниченного тела с источником тепла конечных размеров:
а – аналитическая модель: источник квадратной формы со стороной
2 мм, λ=0,2 Вт/(м·К), а=10-7 м2/с; х0=у0=уm =0, хm =5 мм;
б – дискретная модель цилиндра 40h×40h (h=0,3 мм)с источником
qis= wτ в виде круга радиусом ris= 8h ( ) ; ris=2h ( ); δ –суммарная
погрешность адекватности ИФУТ и дискретной модели
для диапазона а = (0,8 ÷ 5)·10-7 м2/с
Весовые коэффициенты при данных температурах определялись из
решения системы уравнений:
−1
−1
~
(τ i ) ⋅ Ip S (τ i ) = tMT (τi ) ⋅ t M (τ i ) ⋅ Ip S (τ i ) ,
pk = t M
[
]
где tМ(τi) – матрица, составленная из значений температур в трех точках
поверхности полуограниченного тела для i–ых моментов времени; IpS(τ) –
действительное значение интегрального параметра, которое определялось
из ИФУТ:
τ
2


(i )
−1
IpS ( τ ) = C  Lis Qis ( τ ) − α ∫ dτ ∑ Ls ∆ti (τ ) − λ ∫ [t1 (τ ) − t2 (τ )] dτ  ,
i =1

τ
0

где
130
Управление, вычислительная техника и информационные технологии
L(si ) =
ri2 − x 2
L∗2
− ri −1
∫
− ri
ri −1
∫ ( x, y )d y + ∫
dx
− ri2 − x 2
− ri −1
− ri2−1 − x 2
dx
L∗2
ri −1
∫ ( x, y )d y + ∫
− ri2 − x 2
− ri −1
dx
ri2 − x 2
L∗2
∫ ( x, y )dy +
ri2−1 − x 2
(xm − x )2 + y 2 − (x0 − x )2 + y 2
( x,y ) =
+ ∫ dx ∫ ( x, y )dy;
;
2
2
2
2
2 π ⋅ ( x m − x ) + y ⋅ ( x0 − x ) + y
r
− r −x
t1 ( τ) = t ( x0 ,0, τ) ; t2 ( τ) = t ( xm ,0, τ ) ; r1=ris, r2=12h, r3=20h.
ri
i −1
ri2 − x 2
L∗2
2
i
L∗2
2
Суммарная погрешность адекватности ИФУТ и дискретной модели
определялась из выражения:
где


δа ( τ) = 100Qis ( τ) − ∫ qis ( τ)dτ ∫ qis ( τ)dτ ,

 τ
τ
τ
Ls
C 2 ~
λ
Qis ( τ) = Qs ( τ) − ∫ [t1 (τ ) − t2 (τ )] d τ −
∑ pk t(xk , τ).
Lis
Lis 0
Lis k = 0
Результаты исследования данной погрешности для источников тепла с радиусом ris = 0,6 мм и ris = 2,4 мм приведены на графике рис.3, б.
Из результатов данных исследований вытекают следующие выводы: погрешность определения IpS ( τ) по приближенной формуле с использованием информации о температуре в трех точках находится в пределах
δIp= (1…20) %, что не позволит в дальнейшем обеспечить достаточную
точность определения объемной теплоемкости для объекта исследования в
виде полуограниченного тела и данного уравнения баланса тепла; погрешность определения теплопроводности и адекватности ИФУТ изменяется от
δ = (0,1…2) % при ris= 0,6 мм до δ= (1…11) % при ris = 3мм, что дает основание считать точность ИФУТ приемлемой и возможно создание метода измерения теплопроводности на ее основе при правильном выборе размера источника тепла.
При разработке методов измерения ТФС, количества тепла и теплового потока на базе ИФУТ-2 используется подход, аналогичный, применяемому при разработке методов на базе ИФУТ-1, основанному на адекватности данных математических описаний процесса теплопередачи.
Входящий в интегральные формы уравнения теплопроводности интегральный параметр: повторный и однократный интеграл, является единственным источником погрешности их адекватности, которая может быть
уменьшена путем наложения определенных условий на выполнение теплового воздействия и проведения измерения.
Методы измерения, разработанные на базе ИФУТ-2, представлены
на рис. 4, а подробное их описание приведены в литературе [4-10].
131
Известия ТулГУ
ТулГУ. Технические науки. 2014. Вып. 2
Рис. 4. Методы измерения тепловых величин, разработанные
на базе ИФУТ-2
Основным достоинством данных методов является удачное сочетание: степени сложности технической реализации – точности – времени измерения (включая
включая подготовку образца к проведению измерения),
измерения что является основным при разработке средств измерения ТФС.
ТФС Можно
предположить, что приведенный перечень разработок не исчерпывает всех
возможностей и перспективных применений ИФУТ-2 в задачах измерения
тепловых величин.
Список литературы
1. Тихонов А
А.Н.,
Н., Самарский А.А. Уравнения математической физики. М.: Наука, Изд-во.. МГУ,
МГУ 1977. 799 с.
2. Теплофизические измерения и приборы / под ред. Е.С.
Е
Платунова. Л.: Машиностроение,
Машиностроение 1986. 256 с.
3. Филиппов П.И.
П.И Приложение теории теплопроводности к теплофизическим измерениям
измерениям. Новосибирск:
Новосибирск Наука, 1973. 62 с.
4. Азима Ю.И. Применение явного метода идентификации объектов
к решению задач нестационарной теплопроводности // Измерительная техника. 2008. № 6. С. 32––38.
5. Азима Ю.И. Исследование балансного метода измерения неста132
Управление, вычислительная техника и информационные технологии
ционарного теплового потока // Измерительная техника. 2006. № 12.
С.37-42.
6. Азима Ю.И., Бобкова О.А. Измерение теплопроводности низкотеплопроводных материалов методом квазистационарной точки для полуограниченных тел // Известия ТулГУ. Технические науки. Тула: Изд-во
ТулГУ. 2013. Вып. 5. С.92-102.
7. Азима Ю.И. Влияние тепловых сопротивлений при измерении
теплопроводности в диапазоне λ=(10÷400) Вт/(м·К) методом квазистационарной точки// Известия ТулГУ. Технические науки. Тула: Изд-во ТулГУ.
2013. Вып. 5. С.108-117.
8. Азима Ю.И. Метод комплексного измерения теплопроводности и
теплоемкости на базе интегральной формы уравнения Фурье // Известия
ТулГУ. Технические науки. Тула: Изд-во ТулГУ. 2013. Вып. 10. С.261-271.
9. Азима Ю.И. Измерение теплопроводности двухслойных образцов
малой толщины методом квазистационарной точки // Известия ТулГУ.
Технические науки. Тула: Изд-во Тул.ГУ. 2013. Вып. 10. С.244-254.
10. Азима Ю.И. Методы определения температурной зависимости
теплопроводности на основе интегральной формы уравнения Фурье. // Заводская лаборатория. Диагностика материалов. 2002. Т.68. № 5. С.25-30.
Азима Юрий Иванович, канд. техн. наук, доц, juia@bmail.ru, Россия,
Новомосковск, НИ (ф) РХТУ им. Д.И. Менделеева,
REPRESENTATION OF THE HEAT EQUATION IN THE INTEGRAL FORM FOR THE
PROBLEMS OF MEASURING THERMOPHYSICAL VALUES
Yu.I. Azima
Considered new integral form of the heat equation for different single - and multidimensional models of the objects used as samples for measurement of thermophysical properties.. Presents the results of research on the adequacy of the obtained equations. Shows the
use of these equations for the problems of measuring thermal units.
Key words: integral form, the equation of heat conductivity, thermal properties, method of measurement, discrete model.
Azima Yuriy Ivanovich, candidate of technical sciences, Associate Professor,
juia@bmail.ru, Russia, Novomoskovsk, The Novomoskovsk’s Institute (subdivision) of the
Mendeleyev Russian Chemical-Technological University
133
1/--страниц
Пожаловаться на содержимое документа