2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ОБЪЕКТА ОХЛАЖДЕНИЯ
Для описания переноса теплоты через слои охлажденного объекта используем уравнение энергии в одномерном виде [1]:
, (1)
где h – энтальпия (теплосодержание) материала, образующего слой, qх – тепловой поток вдоль координаты х, qv –теплота, выделяемая внутренними источниками в единице объема.
При замене в уравнении (1) производных конечно–разностными приближениями получено:
, (2)
где h′i –энтальпия в ί-точке в момент времени τ+∆τ, qi, qi+1 – соответственно подвод теплоты от узловых точек.
Решая уравнение (2) относительно h′i , получаем:
. (3)
1.3. Расчет индивидуальных характеристик объекта криотерапевтического воздействия.
Толщина оболочки определяется по известному весу объекту криовоздействия. Вес оболочки:
кг.
Средняя площадь наружной поверхности тела м2.
Толщина оболочки:
0,0164 м.
Толщина эпителиального слоя: 0,00205 м.
Толщина жирового слоя: 0,01025 м.
Окончательно принимаем:
2 мм, 10,5 мм.
Исследования выполняются методом вычислительного эксперимента. В эксперименте используются алгоритмы и программные продукты, разработанные при выполнении предыдущих работ курса, в частности методы расчета коэффициента конвективной теплоотдачи, математическая модель оболочки человеческого тела, условия однозначности.
Вычислительный эксперимент сводится к последовательному выполнению численного моделирования процесса конвективного охлаждения объекта ОКВ при 10 значениях температуры газа Т1. Схема процесса охлаждения приведена на рис.1.
Рис.1 Схема эксперимента по оптимизации температурного режима охлаждения в условиях естественной конвекции.
Таблица 1
Результаты уточняющего численного эксперимента.
Вычисляемые параметры, единицы измерений |
Температура теплоносителя, К |
|||||||||
90 |
100 |
110 |
120 |
130 |
140 |
150 |
160 |
170 |
180 |
|
Продолжительность охлаждения τmax, с |
45,0 |
53,0 |
64,0 |
76,0 |
90,0 |
110,0 |
134,0 |
164,0 |
169,0 |
174,0 |
Минимальная температура на поверхности кожи ti=1 при τ=τmax, К |
271,2 |
271,3 |
271,0 |
271,0 |
271,2 |
271,0 |
271,0 |
271,1 |
274,2 |
277,2 |
Минимальная температура на границе слоев ti=n2 при τ=τmax, К |
309,7 |
309,7 |
309,6 |
309,6 |
309,5 |
309,5 |
309,3 |
309,1 |
309,1 |
309,1 |
Максимальный отвод теплоты с поверхности qнар при τ=0, Вт/м2 |
6916,4 |
6240,9 |
5622,7 |
5081,3 |
4578,2 |
4087,4 |
3655,0 |
3260,0 |
2895,2 |
2556,3 |
Максимальный отвод теплоты через границу слоёв qn2 при τ=τmax, Вт/м2 |
48,4 |
49,4 |
50,8 |
52,8 |
56,4 |
64,4 |
78,0 |
99,1 |
98,2 |
96,8 |
Эффективное время, мин |
46,8 |
49,9 |
88,7 |
115,1 |
114,7 |
176,0 |
235,9 |
274,3 |
41,2 |
21,5 |
Общая потеря теплоты с поверхности Qнар при τ=τmax, Дж/м2 |
270,5 |
285,5 |
307,8 |
327,2 |
345,9 |
372,7 |
400,5 |
430,4 |
396,2 |
362,4 |
Общий отвод теплоты через границу слоёв Qn2 при τ=τmax, кДж/м2 |
1,8 |
2,2 |
2,7 |
3,3 |
4,1 |
5,3 |
7,0 |
9,5 |
9,8 |
10,1 |
По результатам численного эксперимента строим диаграмму изменения величины теплового потока
Рис. 2. Изменения плотности теплового потока с поверхности объекта охлаждения.
По результатам численного эксперимента строим диаграмму изменения температуры поверхности ;
Наряду с тепловыми характеристиками следует построить и рассчитать параметры, характеризующие стимулирующее действие ОКВ:
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.