Реферат: Расчет радиаторов

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

Рассмотрим некоторые конкретные примеры составления конечно-разностных схем для узлов двумерной задачи теплопроводности. В этом случае уравнение (2) принимает вид

dT/dx + dT/dy = 0 . (3)

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

Рис.1. Расположение узла внутри двумерного тела толщиной б.

Каждый элементарный прямоугольник (ячейка сетки) имеет длину ‑х и высоту ‑у в направлениях осей х и у. Внутренний узел, обозначенный символом 0, окружен четырьмя соседними узлами: 1,2,3,4. Кондуктивный перенос теплоты, который в действительности происходит в твердом теле через поверхности y*б и x*б (б -толщина тела) будем считать как перенос теплоты от соответствующих узлов к центральному. В установившихся условиях уравнение баланса тепловых потоков для узла 0 при отсутствии внутреннего тепловыделения будет иметь вид

Q(1-0) + Q(2-0) + Q(3-0) + Q(4-0) = 0 , (4)

где Q(I-0) - тепловой поток; индекс (I-0) указывает направление переноса в узлах.

Для определения кондуктивного теплового потока может быть применен закон Фурье

Q = - lamda * F * dT/dn, (5)

где Т - температура, n - направление переноса теплового потока, F - поверхность, через которую переносится тепловой поток.

Для построения расчетной схемы градиент температуры в выражении (5) заменим разностью температур в соседних узлах. В этом случае первый член выражения (4) примет вид

Q(1-0) = y*б*(T[1] - T[0])/x. (6)

Здесь градиент температуры определяется на границе двух узлов 1 и 0, имеющих температуры соответственно Т[1] и Т[0].

Аналогичные уравнения могут быть получены и для остальных трех членов уравнения (1):

Q(2-0) = x*б*(T[2] - T[0])/y, (7)

Q(3-0) = y*б*(T[3] - T[0])/x, (8)

Q(4-0) = x*б*(T[4] - T[0])/y . (9)

Точность аппроксимации градиента зависит от размера ячейки. Если ячейка имеет квадратную форму, то уравнение теплового потока становится независимым от формы тела.

Подставляя зависимости (6)...(9) в выражение (4), можно увидеть, что при постоянном коэффициенте теплопроводности для квадратной сетки (x = y) оно сводится к соотношению между температурами в рассматриваемом узле и близлежащих:

T[1]+ T[2] + T[3] + T[4] - 4*T[0] = 0. (10)

Выражение (10) применимо ко всем внутренним узлам.

Рассмотрим узел, расположенный на поверхности твердого тела, толщиной б в двухмерной задаче (рис.2).

Рис.2.Расположение узлов на поверхности

двумерного тела, омываемого жидкостью

Пусть узел 0, расположенный на границе твердого тела, контактирует с окружающей средой, имеющей температуру Тc. Интенсивность теплообмена с окружающей средой характеризуется коэффициентом теплоотдачи alfa. Узел 0 может также обмениваться кондуктивным потоком теплоты с тремя соседними узлами: 1,2,3. В этом случае тепловой баланс для узла 0 запишется следующим образом:

Q(1-0) + Q(2-0) + Q(3-0) + Q(c-0) = 0, (11)

где Q(c 0)-тепловой поток, передаваемый от среды узлу 0 конвекцией.

По закону Ньютона - Рихмана