Разработка криотерапевтического комплекса по индивидуальным параметрам объекта воздействия, страница 2

2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ОБЪЕКТА ОХЛАЖДЕНИЯ

Для описания переноса теплоты через слои охлажденного объекта используем уравнение энергии в одномерном виде [1]:

                                                ,                                                    (1)

где   h – энтальпия (теплосодержание) материала, образующего слой, qх – тепловой поток вдоль координаты х, qv –теплота, выделяемая внутренними источниками в единице объема.

При замене в уравнении (1) производных конечно–разностными приближениями получено:

                                   ,                    (2)

где  hi  –энтальпия в  ί-точке в момент времени τ+∆τ, qi, qi+1 – соответственно подвод теплоты от узловых точек.

Решая  уравнение (2) относительно hi  , получаем: 

        .                  (3)

1.3. Расчет индивидуальных характеристик объекта  криотерапевтического воздействия.

Толщина оболочки определяется по известному весу объекту криовоздействия. Вес оболочки:

 кг.

Средняя площадь наружной поверхности тела  м2.

Толщина оболочки:

0,0164 м.

Толщина эпителиального слоя: 0,00205 м.

Толщина жирового слоя: 0,01025 м.

Окончательно принимаем:

 2 мм,  10,5 мм.

1.4. Порядок выполнения исследований

Исследования выполняются методом вычислительного эксперимента. В эксперименте используются алгоритмы и программные продукты, разработанные при выполнении предыдущих работ курса, в частности методы расчета коэффициента конвективной теплоотдачи, математическая модель оболочки человеческого тела, условия однозначности.

Вычислительный эксперимент сводится к последовательному выполнению численного моделирования процесса конвективного охлаждения объекта ОКВ при 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. Изменения плотности теплового потока с поверхности объекта охлаждения.

По результатам численного эксперимента строим диаграмму изменения температуры поверхности ;

Наряду с тепловыми характеристиками следует построить и рассчитать параметры, характеризующие стимулирующее действие ОКВ: