.
После завершения вычислений на очередном участке переменная приобретает новое значение:
.
Начальный участок объекта замораживания () находится в тепловом контакте с теплоотводящим элементом криодеструктора . В соответствии с принятой моделью объекта замораживания теплообменник представляет собой объект с полусферической поверхностью. В простейшем случае температура теплообменника постоянна .Тепловой поток между объектом и теплообменником описывается выражением:
,
где – эффективный коэффициент теплоотдачи.
На практике часто используют инструменты с охлаждаемым теплообменником, т.е. температура меняется во времени. Может быть задан закон изменения температуры теплообменника или закон изменения мощности системы охлаждения .
РАЗРАБОТКА ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ ДЛЯ ПРОВЕДЕНИЯ ЧИСЛЕННОГО ЭКСПЕРИМЕНТА.
Для проведения численного эксперимента используем редактор Excel [3]. Для сокращения затрат времени по созданию программного обеспечения можно воспользоваться уже разработанными программными продуктами. Наиболее близким по структуре программным модулем является Excel книга “Криоконсервация” разработанная при подготовке и выполнении лабораторной работы №4 [4]. Для редактирования открываем файл “Криоконсервация” и сохраняем его под новым именем “Криодеструкция”.
Производим редактирование листа “Эксперимент”. В таблице характеристика материалов приводим характеристики теплообменника и объекта криохирургического воздействия.
Свойства объекта вносим из табл. 1. Свойства материала теплообменника определяем по справочной литературе [5]. Следует выбирать металлы с высокой теплопроводностью: медь, латунь, серебро.
Таблица 2
Физические свойства металлов
Метал |
Свойства |
||
, кг/м3 |
, Вт/м·К |
, Дж/кг·К |
|
Медь |
8940 |
400 |
390 |
Латунь |
8500 |
110 |
377 |
Серебро |
10500 |
430 |
235 |
Аллюминий |
2700 |
240 |
800 |
Полученная таблица используется для описания элементов криохирургической системы.
B |
C |
D |
|
1 |
ХАРАКТЕРИСТИКА МАТЕРИАЛОВ |
||
2 |
ТЕПЛООБМЕННИК |
ОБЪЕКТ |
|
3 |
Теплоемкость, Дж/кг·К |
390 |
4000 |
4 |
Теплопроводность, Вт/м·К |
400 |
0,5 |
5 |
Плотность, кг/м3 |
8940 |
1010 |
6 |
Влажность |
0 |
0,9 |
7 |
Температура дефростации, К |
270 |
|
8 |
Температура замораживания, К |
240 |
|
9 |
Удельное тепловыделение, Вт/кг |
10,6 |
|
10 |
Теплота дефростации |
0 |
300600 |
Рис. 3 Характеристика элементов криохирургической системы
На листе “Эксперимент” необходимо сформировать блок условий однозначности, которые начинаются с геометрического описания элементов. Задается (см. рис. 4) радиус поверхности теплообменника, пассивной границы объекта криовоздействия и удаляемого фрагмента.
А |
B |
C |
D |
||
11 |
|||||
12 |
ХАРАКТЕРИСТИКА КРИОХИРУРГИЧЕСКОЙ СИСТЕМЫ |
||||
13 |
|||||
14 |
Радиус теплообменника, мм |
10 |
0,01 |
||
15 |
Радиус пассивной границы, мм |
50 |
0,05 |
||
16 |
Радиус удаляемого фрагмента,мм |
20 |
0,02 |
||
17 |
|||||
18 |
ПАРАМЕТРЫ СЕТКИ РАЗБИЕНИЯ |
||||
19 |
Шаг по времени, с |
0,05 |
|||
20 |
Шаг по координате, мм |
2,5 |
|||
21 |
Число участков |
20 |
|||
Рис.4 Геометрическое описание криохирургической системы и сетки разбиения.
Затем определяются параметры сетки разбиения: шаг по времени и число участков в объекте.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.