

1.Механизм процесса обусловлен видом носителя тепла. В металлах электроны при перемещении переносят эл. заряд и тепловую энергию теплового движения-процесс теплопроводности. Чем выше температура чистого метала, тем хуже он проводит тепло. В диэлектриках механизм тепла осущ за счет колебания узловых решеток. Если в каком то месте диэлектрик нагревают, то усиливающиеся в нем колебания узлов распространяются в теле в виде затухающей волны. В полупроводниках, при низких температурах тепло распространяется как в диэлектриках. При средних вклад электронной проводимости возрастает, а при высоких температурах она становится преобладающей. В неподв.газовых слоях носителями тепла явл хаотически движ молекулы.Если газ в каком то месте нагреть то увеличивается среднеквадрт скорость движ молекул,которые сталкиваясь с более удаленными молекулами увел их кинетическую энергию и увел и температура газа. Жидкости в течении малого времени(релаксации) можно рассматривать как структуры. Затем эти структуры разрушаются и поведение жид соответствует газовому состоянии. В каждом состоянии действует механизм тепла.
2. Гипотеза Фурье. Фурье |
|
предположил, что существует прямая |
|
пропорциональность между величинами |
и gradT, т.е. |
||
q |
|q |~| gradT |
|
|
имеем также |
q ~ |
|
. (1.1)Учитывая разнонаправленность указанных в (1.1) векторов, gradT . (1.1 )Чтобы перейти в (1.1 ) от пропорции к равенству, Ж.-Б.
Фурье ввел коэффициент пропорциональности и получил зависимость
q
gradT
,
(1.2)представляющую собой математическую запись его гипотезы. Величина численно
равна
|
|
|
|
|
q |
, |
|
gradT |
|||
|
|
|
Вт |
|
|
|
м К |
(1.3)и совпадает с плотностью теплового потока при значении
|gradT|, равном 1 К/м. Фурье назвал коэффициентом теплопроводности материала тела. Зависит величина от вида материала тела, его пористости, влажности и, что очень существенно, от самой температуры Т.
7. Краевая задача нестационарной теплопроводности для определения нестационарного температурного поля решают уравнение Фурье совместно с присоединенными к нему начальным условием и граничными условиями. Краевая задача
нестационарной |
теплопроводности |
имеет |
вид. T (x, ) |
|
a |
s 1 |
T (x, ) |
||||
|
|
|
|
|
|
|
x |
|
|
, |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
xs 1 x |
|
x |
|
(1.28)
(1.29)
s 1, 2, 3; |
||
|
T (x, ) |
|
x |
||
|
0; |
0 |
T (l |
|
|
0 |
x l |
|
0 |
|
x l |
, T (x) T |
, |
0 |
0 |
|
, ) T |
|
, |
0 |
, |
f ,0 |
|
0, |
0 x l |
0 |
, |
|
|
|
(1.30) T |
0, |
0 |
. |
|
|
x |
|
||
|
x 0 |
|
|
(1.31)В
записи краевой задачи отражен факт симметричного развития температурного поля относительно плоскости (оси, центра) симметрии тела, т.е. относительно x = 0.В краевой задаче известны форма тела (величина s), его характерный размер l0, а также величины a,, T0, , Tf ,0 , т.е. известны параметры задачи, изменяющиеся от одной конкретной
ситуации к другой, и отыскивается температурное поле T(x, ), так что в итоге получаем, что температура T отыскивается из решения задачи (1.28) – (1.31) в виде зависимости от аргументов x, и от параметров s, a, , T0, , Tf ,0 , l0:
T f (x, , s, a, ,T0 , ,Tf ,0 , l0 ). (1.32)Таким образом, подлежит определению функция Т
девяти переменных, теорема существования и единственности которой для краевой

задачи (1.28)–(1.31).
12. Два типа инженерных задач С помощью графиков решают два типа инженерных задач:1) определение температуры в указанной точке тела по истечении
которое будет достигнута заданная температура T(x, |
|
|
|
|
При этом |
|||||||
считаются |
известными величины |
a, , , T0 , Tf ,0 , l0 , s .Для решения первой |
|
и второй |
||||||||
задачи сначала вычисляют величину критерия Био Bi |
l0 / .Далее, при |
решении |
||||||||||
|
|
|
|
|
|
|
|
Fo a / l 2 |
и по |
|||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
графику для указанной точки тела находят температуру |
|
|
|
|
|
, |
||||||
|
|
T (x, ) T |
, откуда следует T (x, ) (Tf ,0 |
T0 ) T0 .При решении второй |
||||||||
равную |
|
|
0 |
|||||||||
T |
|
T |
||||||||||
|
|
f ,0 |
|
|
|
|
|
|
|
|
||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
обезразмеривают по правилу |
T (x, ) T |
|
||||
задачи заданную температуру T(x, |
|
|
|
0 |
и по |
|||||||
T |
|
T |
||||||||||
|
|
|
|
|
|
|
|
f ,0 |
|
|||
|
|
|
|
|
|
|
|
|
|
0 |
|
графику находят безразмерное время Fо ее достижения (путь б на рис. 1.9):
Fo
a /
l 2 0
, откуда
Fo l 2
0
/
a
.
3. 5. Вывод уравнения Фурье для одномерной задачи теплопроводности.
Пусть за время d температура в выделенном объеме изменится на величину
так что изменение внутренней энергии за единицу времени составляет |
U |
cm |
|
|
|||
|
|
dT,
T |
, |
|
|
||
|
(1.9)Единственной причиной изменения внутренней энергии во времени является
разность «втекающего» q(x) и «вытекающего» q(x + dx) количества тепла, т.е. верно
также равенство U q(x) q(x dx) . |
(1.10) U q (x)dx . (1.10 )Приравнивая правые |
|||||
|
|
|
|
|
|
|
части (1.9 ) и (1.10 ), получаем c |
T |
q (x ) . (1.11)Таким образом, температура в том |
||||
|
||||||
|
|
|
|
|
||
или ином месте пластины изменяется во времени |
T |
|
||||
|
0 лишь в том случае, когда |
|||||
|
|
|
|
|
|
изменяется от места к месту плотность теплового потока (q (x ) 0).На основании (1.3)
имеем q(x ) |
T |
, |
так что уравнение (1.11) при = const принимает |
c |
T |
|
2T . |
|
x |
|
|
|
|
|
x 2 |
(1.11 )Это и есть |
|
уравнение Фурье, описывающее нестационарное |
одномерное |
(изменяющееся лишь по 0x) температурное поле в пластине.
4. 6.Вывод уравнения Фурье для двумерного температурного поля. двухмерное температурное поле, формируется в том случае, когда вектор плотности теплового потока «втекающего» и «вытекающего» в элементарный участок сечения призмы имеет ненулевые компоненты в направлениях 0x и 0y
U c T dx dy L,

c T dx dy L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q(x ) q( y ) dx dy L. |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
2 |
|
|
Наконец, учитывая, что при |
= const q(x ) |
, |
|
q(x ) |
T |
, |
|
|||||||
x |
|
x |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
q( y) T , |
q ( y) 2T |
, имеем в окончательном виде уравнение Фурье для |
||||||||||||
y |
|
|
y2 |
|
|
|
|
|
|
|
|
|
|
|
описания двухмерного температурного поля |
|
|
c |
T |
|
2T |
|
2T |
|
|||||
|
|
|
|
x 2 |
y2 |
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
8.ГУ I-го и II-го рода
Если из физических соображений или в результате проведенных измерений известна температура TW на поверхности Г тела, то мы располагаем граничными условиями первого рода в форме
T (M , ) T |
(M , ) f |
(M , ), M |
W |
1 |
|
В простейшем случае в течение всего поверхностях тела температура одинакова, и тогда
, 0. |
(1.20) |
процесса во всех точках на всех вместо (1.20) имеем ГУ-I в виде
T (M , ) TW , M , 0. |
(1.20 ) |
Если известна плотность теплового потока q на поверхности тела, то к уравнению Фурье присоединяют ГУ-II в форме
q(M , ) f2 (M , ), M , 0. |
(1.21) |
С учетом (1.4) формула (1.21) принимает вид
gradT |
|
0 |
f |
|
(M , ), M , 0 |
||||
n |
|
2 |
|||||||
|
|
|
|
|
|
|
|
||
или |
|
|
|
|
|
|
|
|
|
|
0 |
|
f |
|
(M , ) / , M , 0. |
||||
gradT n |
|
|
2 |
||||||
|
|
|
|
|
|
|
|
(1.21 )
(1.21 )

9.ГУ III-го и IV-го рода
Граничные условия третьего рода (ГУ-III) присоединяют к уравнению Фурье в том частном случае, когда тело омывается потоком жидкости (газа), температура которого
Tf ,0 |
на удалении от тела известна При этом плотность теплового потока, передаваемого |
от движущейся среды к поверхности тела, полагают пропорциональной разности
температур TW –Tf ,0 |
(температура |
поверхности TW неизвестна и сама подлежит |
|
определению) |
|
|
|
q qW ~ TW |
Tf ,0 . |
(1.22) |
Чтобы перейти в (1.22) от пропорции к равенству, вводится коэффициент пропорциональности , называемый коэффициентом теплоотдачи, так что имеем
q |
|
qW (TW Tf ,0 ), 0. |
(1.22 ) |
|
|||
В формуле (1.22 ) считаются известными лишь величины Tf ,0 |
и . Величина |
численно равна плотности теплового потока,
TW Tf ,0 |
= 1K: |
|
q |
, |
|
W |
|||
|
|
||
T |
T |
f ,0 |
|
W |
|
передаваемого от поверхности тела при
|
Вт |
|
(1.23) |
||
|
|
2 |
|
. |
|
м |
К |
|
|||
|
|
Граничные условия четвертого рода относятся к специфическому случаю теплового контакта между двумя твердыми телами (При этом возможен случай идеального теплового контакта и неидеального теплового контакта когда поверхности Г тел № 1 и № 2 разделены газовой прослойкой, слоем окислов, слоем масла и т.п.Ясно, что в обоих случаях плотности теплового потока, пересекающего поверхность Г слева (Г–0) направо (Г+0), совпадают, так что с привлечением (1.4) имеем
|
gradT |
(1) |
|
0 |
|
|
|
gradT |
(2) |
|
0 |
|
, 0. |
|
n |
|
0 |
2 |
|
n |
|
0 |
|||||
1 |
|
|
|
|
|
|
|
|
|
|
(1.25)
В случае идеального теплового контакта на поверхностях Г–0 и Г+0 в течение всего процесса совпадают и температуры контактирующих тел:
T |
(1) |
|
T |
(2) |
|
, |
0, |
|
0 |
|
0 |
||||
|
|
|
|
|
|
(1.26)
а в случае неидеального теплового контакта имеет место скачок температуры T, формирующийся на термическом сопротивлении, разделяющем оба тела, т.е. выполняется равенство
T |
(1) |
T |
(2) |
|
|
||
|
0 |
|
0 |
T , 0.
(1.27)

10.нестационарная теплопроводность при Если известна плотность теплового потока
Фурье присоединяют ГУ-II в форме
q(M , ) f |
(M , ), M , |
2 |
|
ГУ 2
q на поверхности тела, то к уравнению
0. |
(1.21) |
С учетом (1.4) формула (1.21) принимает вид
|
0 |
f |
(M , ), M , 0 |
gradT n |
|
||
|
2 |
|
(1.21 )
или
|
0 |
f |
|
(M , ) / , M , 0. |
(1.21 ) |
gradT n |
|
2 |
|||
|
|
|
|
11. Безразмерная форма краевой задачи теплопроводности при ГУ-III-го рода
Вместо «размерной» температуры T(x, ) [T0; Tf ,0 ] вводится безразмерная
относительная температура (x, ) по правилу
|
T T |
||
|
|
0 |
|
|
|
|
|
|
T |
f ,0 |
T |
|
|
0 |
[0;1],
так что в задаче (1.28)–
(1.31) надо везде заменить T на , подставив T (Tf ,0 T0 ) T0 .
Далее, вместо размерной протяженности x [0; l0] вводится безразмерная протяженность = x/l0 [0;1], так что в исходной задаче надо везде заменить x на x = l0.
Задача (1.28) – (1.31) принимает в результате таких подстановок вид
|
|
|
|
a |
|
|
s 1 |
|
|
(1.28 ) |
s 1, 2, 3; |
|
|
l |
2 |
|
s 1 |
|
|
, |
|||||
|
|
|
|
|
|
|
|
|||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
( ) 0, |
0, 0 |
1, |
(1.29 ) |
0,
0
1
l0
,
|
|
|
1 |
|
|
1 |
, |
0, |
|
|
|||||||
|
|
|||||||
|
|
|
||||||
|
|
1 |
|
|
|
|
||
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
(1.30 )
|
0, |
0. |
|
|
|||
0 |
|
||
|
|
(1.31 ) |
|
|
|
|
|
1 1 , |
0, |
(1.30 ) |
|
|
|
|
|
||||||
|
l0 |
|
|||||||
|
|
|
|
1 |
|
|
|||
|
|
0, |
0. |
|
|
(1.31 ) |
|
||
|
|
|
|
|
|||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сформируем безразмерные комплексы |
a |
в (1.28 ) и |
l0 |
в (1.30 ). Безразмерный |
|
l 2 |
|
||||
|
|
|
|||
|
0 |
|
|
|
комплекс |
a / |
(Fo a / l02 ) ,
l 2 0
а
представляет собой безразмерное время и называется числом Фурье безразмерный комплекс l0 / представляет собой известную
безразмерную интенсивность внешнего теплообмена потока с поверхностью тела и называется критерием Био (Bi l0 / ) . Число Фурье Fо содержит в себе аргумент
задачи и поэтому является ее безразмерным аргументом, а критерий Био Bi составлен

из известных при постановке задачи параметров.
В конечном виде имеем следующую задачу нестационарной теплопроводности относительно искомой температуры ( , Fо)
Fo
( )
|
|
1 |
|
|
|
s 1 |
|
|
s 1 |
|
|
, |
|||
|
|
|
|
|
|
||
|
|
|
|
|
|||
0, |
Fo 0, 0 |
|
(1.28 )
1, |
(1.29 ) |
s 1, 2, 3; |
Fo 0, 0 1, |
|||||
|
Bi 1 |
|
, |
Fo 0, |
||
|
1 |
|||||
|
|
|
|
|||
1 |
|
|
|
|
||
|
|
|
|
|
(1.30 )
0, Fo 0.0
(1.31 )
Решение задачи (1.28 )–(1.31 ) переменных (вместо девяти в (1.32)) как
отыскивается в виде функции |
от четырех |
f ( , Fo, s, Bi). |
(1.33) |
14. Сеточный метод решения одно- и двумерных задач нестационарной теплопроводности
Метод сеток практически совпадает с методом элементарного теплового баланса. Отличие между ними состоит в том, что, во-первых, метод сеток обосновывается формальной дискретизацией уравнений исходной краевой задачи нестационарной теплопроводности и, во-вторых, полученный таким образом разностный аналог уравнения нестационарной теплопроводности относят ко всем элементарным слоям, на которые мысленно разбивается исходная геометрическая область протекания процесса, т.е. не рассматривают отдельно, как в методе элементарного теплового баланса, пристенные слои Исходное нелинейное уравнение нестационарной теплопроводности, имеющее для неограниченной пластины в
одномерной постановке вид |
c T T |
T |
|
|
|
T |
T |
, |
|
|
|
|
|||||
|
|
|
|
x |
|
x |
|
использовании |
неявной безытерационной |
схемы |
|
дискретизируется |
|||||||||
середины |
|
|
i-го |
|
|
слоя |
(рис. |
1.11) |
|||||
образом: c T T |
T |
c T |
T |
T |
(n 1) T (n) |
c(n) (n) |
T (n 1) |
T |
(n) |
||||
|
|
i |
i |
|
i |
i |
|
, |
|||||
|
|
|
|
|
|
|
|
||||||
|
i |
i |
i |
|
|
i |
i |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
(1.68) при
относительно
следующим
(1.69)
Приравнивая правые части формул (1.69) , (1.70), получаем разностный аналог уравнения (1.68) в виде
|
|
T (n 1) T (n) |
|
|
|
|
|
T (n 1) T |
(n 1) |
|
|
|
T (n 1) |
T (n 1) |
|
|
|||||||||
(n) |
(n) |
i |
i |
|
|
|
(n) |
|
i 1 |
i |
|
|
|
|
(n) |
|
i |
|
i 1 |
|
|
||||
ci |
i |
|
|
|
i |
1 |
|
|
|
|
|
i |
1 |
|
|
|
|
|
|
||||||
|
|
|
|
x 2 |
|
|
|
x 2 |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
T (n 1) T (n 1) |
|
|
|
|
T (n 1) T (n 1) |
|
|
|
|
||||||||||
|
|
|
(n) |
|
|
i 1 |
|
i |
|
(n) |
|
i |
|
i 1 |
|
|
|
|
|||||||
|
|
|
i |
1 |
|
|
|
|
|
i |
1 |
|
|
|
|
|
, |
|
(1.71)где |
||||||
|
|
|
|
|
|
|
x 2 |
|
|
|
x 2 |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
||
i 2,3,...,m 1.Метод |
сеток |
позволяет |
|
решать и |
|
многомерные |
нелинейные |
задачи |
|||||||||||||||||
нестационарной теплопроводности. В этом случае на тело наносится сетка, |
т.е. его, |
например при рассмотрении двухмерного температурного поля, мысленно делят на элементарные прямоугольники со сторонами x и y , которые представляют собой
шаги по пространственным переменным x и y, при выборе в качестве шага по

времени.Можно
нестационарной
показать, что в |
этом случае исходное |
|
нелинейное |
уравнение |
||||||||
теплопроводности |
вида c T T |
T |
|
|
|
T |
T |
|
|
|
T |
T |
|
|
|
|
|
|
|
||||||
|
|
|
|
x |
|
x |
|
|
|
|
|
|
|
|
|
|
|
y |
|
y |
|
|
|
(1.78)при использовании безытерационной неявной схемы имеет следующий |
|||||||||||||||||||||
конечно-разностный |
|
|
|
T |
|
T |
|
|
|
T |
|
|
T |
|
аналог: |
|||||||||
|
|
|
|
T |
(n 1) |
T |
(n) |
|
|
|
(n 1) |
(n 1) |
|
|
(n 1) |
(n 1) |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
c |
(n) |
|
(n) |
i, j |
i, j |
(n) |
|
|
i 1, j |
|
i, j |
(n) |
|
i, j |
|
|
i 1, j |
|
||||||
|
i, j |
|
|
|
|
|
1 |
|
|
|
|
2 |
|
|
1 |
|
|
|
|
2 |
|
|||
i, j |
|
|
|
|
i |
, j |
|
x |
|
i |
, j |
|
x |
|
|
|||||||||
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i
|
T |
(n 1) |
T |
(n 1) |
|
|||
|
|
|
|
|
|
|||
(n) |
i, j 1 |
|
|
|
i, j |
(n |
||
|
1 |
|
y |
2 |
|
|
|
|
i, j |
|
|
|
i, |
||||
2 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
1,2,...,M ; |
|
|
j |
1,2,..., N |
|
T |
(n 1) |
T |
(n 1) |
||
|
|
|
|
|||
) |
i, j |
|
|
i, j 1 |
||
j |
1 |
|
y |
2 |
|
|
2 |
|
|
|
|||
|
|
|
|
|
|
.
,
(1.79)
15. Явная и неявная схемы численного определения температурных полей
В зависимости от того, по какому распределению температуры в теле вычисляются
величины qi 1, i и qi, i 1 |
, различают явную и неявную схемы численного решения задачи |
|||||||||||||||
теплопроводности. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При явной схеме |
qi 1, i |
и |
qi, |
i 1 |
определяют по предшествующему распределению |
|||||||||||
температуры (по отношению к искомому), т.е. следующим образом: |
|
|||||||||||||||
q |
|
|
|
T |
(n) |
T |
(n) |
q |
|
|
|
T |
(n) |
T |
(n) |
, |
|
|
|
|
|
, |
|
|
|
|
|
|
|||||
|
(n) |
(n) |
i 1 |
i |
|
|
(n) |
(n) |
i |
|
i 1 |
|
||||
|
i 1, i |
|
i 1/ 2 |
|
x |
|
|
i, i 1 |
|
i 1/ 2 |
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
так что имеем
(n) |
|
(n) |
T |
(n |
1) |
T |
(n) |
|
|
(n) |
T |
(n) |
T |
(n) |
|
(n) |
T |
(n) |
T |
(n) |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
i |
|
|
|
i |
|
|
|
i 1 |
|
i |
|
|
|
i |
|
|
|
i 1 |
, |
||||||||
ci |
i |
|
|
|
|
|
x i 1 / 2 |
|
x |
|
|
i 1 / 2 |
|
x |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где коэффициенты теплопроводности вычисляются как |
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
T |
(n) |
T |
(n) |
|
|
|
|
|
|
|
|
T |
(n) |
T |
(n) |
|
|
|||
|
|
|
(n) |
|
|
f |
|
i 1 |
i |
|
|
|
(n) |
|
|
|
f |
i |
|
|
i 1 |
|
|
|||||
|
|
|
|
i 1 / 2 |
|
|
|
|
|
, |
|
|
i 1 / 2 |
|
|
|
|
|
|
. |
|
|||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
При |
неявной |
схеме |
величины |
qi 1, |
i |
и |
|
qi, |
i 1 |
|
определяют |
распределению температуры, так что вместо (1.63 ) получают формулу
(1.63 )
по искомому
c |
|
|
|
T |
(n 1) |
T |
(n) |
|
|
|
T |
(n 1) |
T |
(n 1) |
|
|
T |
(n 1) |
T |
(n 1) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
||||
|
(n) |
|
(n) |
i |
|
i |
|
(n) |
|
i 1 |
|
i |
|
(n) |
i |
|
|
|
i 1 |
|
||||
i |
|
i |
|
|
|
|
i 1 / 2 |
|
x |
2 |
|
|
i 1 / 2 |
|
x |
2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
i [2, m 1]. |
|
|
|
|
|
|
|
|
(1.63 ) |
Явная схема счета является устойчивой при выборе соотношения между шагами x и по правилу
(in) |
|
1 |
. |
c(n) (n) x2 |
2 |
||
i i |
|
|
|
Неявная же схема счета является абсолютно устойчивой, так что, вообще говоря, не имеется ограничений на выбор величин x и .

13. Метод элементарного теплового баланса при численном решении задач
теплопроводности. |
|
|
|
|
Определение температурных полей в телах сложной формы |
при зависящих от |
|||
температуры характеристиках материала тела (с = с(T |
= |
T |
= |
T |
связано с решением проблемы многомерности и нелинейности, преодолеваемой, в общем случае, с привлечением возможностей ЭВМ.
При численном решении температуру определяют в дискретных точках x, и в дискретные моменты
При этом полагают известным распределение температуры в n-й момент времени, т.е. поле температуры
T |
(n) |
, T |
(n) |
, T |
(n) |
,...,T |
(n) |
, T |
(n) |
, T |
(n) |
,...,T |
(n) |
, T |
(n) |
, T |
(n) |
, |
|
|
|
|
|
|
|
|
|
||||||||||
W 1 |
1 |
2 |
i 1 |
i |
|
i 1 |
m 1 |
m |
W 2 |
|
(n+1)-й
момент времени.
Изменение внутренней энергии в i-м пространственном слое за 1 с будет равно
c(n) (n) |
T |
(n 1) |
T |
(n) |
|
|
|
|
|
|
|
|
|
||||
i |
|
i |
1 x q |
q |
i, i 1 |
, |
||
i |
i |
|
|
i 1, i |
|
|
||
|
|
|
|
|
|
|
(1.63)
где Ti |
(n) |
и Ti |
(n 1) |
- предыдущее и последующее значения температуры в середине i- |
|
|
го слоя пластины; |
qi 1, i |
и |
qi, i 1 |
- плотности теплового потока, ―втекающего‖ из (i–1)-го |
слоя в i-й слой и «вытекающего» из него в сторону (i+1)-го слоя; |
(n) |
и |
(n) |
– значения |
|
ci |
i |
|
|||
удельной теплоемкости и плотности вещества, выбранные из таблицы |
(массива) их |
зависимости от температуры
T |
(n) |
|
|
i |
|
в середине i-го слоя:
c |
(n) |
f (T |
(n) |
) |
|
|
|||
i |
i |
|
|
и
|
(n) |
|
|
|
|
|
i |
|
f
(T |
(n) |
) |
|
||
i |
|
|
.
В зависимости от того, по какому распределению температуры в теле вычисляются
величины qi 1, i и qi, i 1 |
, различают явную и неявную схемы численного решения задачи |
теплопроводности. |
|

17. Неявная схема численного решения двухмерной задачи теплопроводности
Для решения системы уравнений нестационарной теплопроводности, ее разностный аналог ее нужно замкнуть путем присоединения двух дополнительных уравнений, представляющих собой конечно-разностный аналог граничных условий. При их формулировке полагают, что между серединами первого фиктивного слоя и примыкающего к нему первого (по оси Ох) слоя пластины и между серединами последнего слоя пластины и примыкающего к нему второго фиктивного слоя температуры распределены в пространстве линейно, т.е. выполняются равенства
T |
(n 1) |
T |
(n 1) |
(n 1) |
|
|
|
|
|||
1 |
2 |
, |
|||
|
|
2 |
|
TW1 |
|
|
|
|
|
|
T |
(n 1) |
T |
(n 1) |
(n 1) |
|
|
|
|
|||
m 1 |
m 2 |
. |
|||
|
|
2 |
|
TW 2 |
|
|
|
|
|
|
(1.72)
(1.73)
С учетом (1.72), (1.73) граничные условия записываются как
|
|
T |
(n 1) |
T |
(n 1) |
|
(n 1) |
|
|
(n 1) |
|
T |
(n 1) |
T |
(n 1) |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
, |
|
|
|||||
(n) |
1 |
|
|
|
2 |
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|||||||||||
|
11 |
|
|
|
|
x |
|
|
|
|
1 |
|
|
|
f 1,0 |
|
|
|
|
2 |
|
|
|
|
|
|
|
||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
T |
|
|
T |
|
|
|
|
T |
|
|
T |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
(n 1) |
(n 1) |
|
|
|
(n 1) |
(n 1) |
T |
|
|
|
. |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
(n) |
|
|
m 1 |
|
|
m 2 |
|
|
(n 1) m 1 |
|
m 2 |
|
|
(n 1) |
|
|
|||||||||||||
|
m 11 |
|
|
|
|
x |
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
f 2,0 |
|
|
|
|||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
третьего рода (1.66), (1.67)
(1.74)
(1.75)
Если же заданы граничные условия второго рода (плотности тепловых потоков
(n 1) |
(n 1) |
qW1 |
и qW 2 |
на ограничивающих поверхностях пластины), то с учетом (1.72), (1.73) они
будут представлены в виде
|
|
T |
(n 1) |
T |
(n 1) |
|
|
|
|
|
|
(n) |
1 |
2 |
|||
1 |
1 |
|
|
x |
|
2 |
|
|
|
||
|
|
|
|
|
q(n 1) W ,1
,
(1.76)
|
T |
(n 1) |
T |
(n 1) |
|
|
|||
(n) |
m 1 |
m 2 |
||
m 11 |
2 |
|
x |
|
|
|
|
q(n 1) W ,2
.
(1.77)
Если направления векторов |
|
(n 1) |
и |
|
(n 1) |
совпадают с направлением оси Ox, то в |
q |
q |
|||||
|
W1 |
|
W 2 |
|
(1.76), (1.77) их численные значения положительны.

18.Физический принцип стационарности температурных полей.
По истечении достаточно длительного времени от начала процесса теплопроводности (теоретически при ) температурные изменения в теле во времени прекращаются и наступает режим стационарной теплопроводности, когда
T |
0 |
. В этом режиме при одинаковой температуре омывающей тело среды в нем |
|
|
|||
|
|
отсутствуют градиенты температуры (устанавливается однородное температурное поле) и отсутствует теплоперенос.
При рассмотрении стационарной теплопроводности обычно решаются два вопроса: 1) определение температуры в любом месте тела, 2) нахождение величины стационарного теплового потока через конструкцию.
Физический принцип: в стационарном тепловом режиме одинаков тепловой
поток, пересекающий любую изотермическую поверхность в теле и любую его часть, ограниченную изотермическими поверхностями.
19. Стационарное температурное поле в неограниченной пластине
Стационарный тепловой поток Q(x) через отстоящий на расстоянии х участок изотермической поверхности площадью F(x) (рис. 1.4) за единицу времени равен
Q(x)
имеем
q(x)F (x) |
|||
dT |
Q(x) |
|
|
F (x) |
|||
|
|
dT |
F (x) |
. Разделяя переменные, |
||
dx |
||||
|
|
|
||
dx (1.81). |
Принимаем, что |
коэффициент теплопроводности одинаков, т.е. = const. Кроме того, в стационарном тепловом режиме всегда Q(x) = const, а в пластине и F(x) = const в случае одномерного температурного поля. Интегрирование (1.81) дает
T |
Q |
|
F |
||
|
x
C1
, где C1 - произвольная постоянная.
Таким образом, в указанных выше предположениях распределение температуры Т по координате х в неограниченной пластине подчиняется линейному закону (рис. 1.13).