Исследование процесса консервации биопрепаратов: Методические указания к лабораторной работе № 4 по курсу «Криогенные технологии в медицине и биологии», страница 4

Рис. 4. Алгоритм математической модели контейнера для биопрепарата

Разработка программного обеспечения численного эксперимента

Для проведения численных экспериментов необходимо разработать программное обеспечение, в котором будет реализована систематическая модель объекта криоконсервации.

Для этих целей можно воспользоваться редактором "Exel"[1,3]. Программа оформлена в виде "Exel" книги, содержащей лист данных о теплофизических свойствах азота - "Свойства", лист на котором выполняется вычислительный эксперимент – "Эксперимент", а также лист для обработки результатов численного эксперимента – "Результат".

Данные листа "Свойства" составляются путём комбинирования значений ρ (плотности), ср (изобарной теплоемкости), μ (вязкости), β (коэффициент объемного расширения), λ (теплопроводности), и h (энтальпии) азота при температурах от 70 до 320 К и атмосферном давлении, из электронной базы данных "Азот"[8].

A

B

C

D

E

F

G

1

ρ

ср

h

µ

λ

β

2

70

837,866455

2023,857

110,510864

0,0002219836

0,149313

0,0143

3

71

833,697510

2065,252

112,554115

0,0002114583

0,147933

0,0141

4

72

829,443359

2091,495

114,635315

0,0002016271

0,146470

0,0139

5

73

825,108398

2107,342

116,733398

0,0001924382

0,144930

0,0137

6

74

820,698486

2115,186

118,847801

0,0001838496

0,143323

0,0135

7

75

816,217285

2118,043

120,963112

0,0001758174

0,141656

0,0133

8

76

811,673584

2117,216

123,081772

0,0001683045

0,139938

0,0132

9

77

807,068359

2114,184

125,198708

0,0001612723

0,138175

0,0130

Рис.5.1. Теплофизические свойства азота. Фрагмент листа "Свойства"

Формирование листа "Эксперимент" начинается с характеристики материалов подвергающихся охлаждению. Для этого в верхней части листа формируется и заполняется таблица (см. рис. 5.2.).

А

В

С

D

1

ХАРАКТЕРИСТИКА МАТЕРИАЛОВ

2

КОНТЕЙНЕР

ПРЕПАРАТ

3

Теплоёмкость, Дж/кг К

1640

4000

4

Теплопроводность, Вт/м К

0,41

0,5

5

Плотность, кг/м3

940

1010

6

Влажность

0

0,9

7

ТЕМПЕРАТУРА ДЕФРОСТАЦИИ, К

270

8

ТЕМПЕРАТУРА ЗАМОРАЖИВАНИЯ, К

240

9

ТЕПЛОТА ДЕФРОСТАЦИИ

0

300600

Рис.5.2. Характеристика материалов. Формирование листа "Эксперимент"

Затем определяются геометрические размеры объекта криоконсервации. Например, Dк=12,5 мм, Dпр=11 мм (данные значения используются для отладки программного обеспечения и впоследствии могут изменяться).

Вычисление осуществляется в радиальных координатах, поэтому по заданным диаметрам вычисляются радиусы, а полученные значения пересчитываются в метры. (рис. 5.3).

В

С

D

E

10

11

                                                                                        РАДИУС, мм   РАДИУС, м

12

ДИАМЕТР КОНТЕЙНЕРА, мм

12,5

6,25

0,00625

13

ДИАМЕРТ ПРЕПАРАТА, мм

11

5,5

0,0055

14

ВЫСОТА ОБЪЕКТА, м

0,1

15

ОБЪЕМ ПРЕПАРАТА, м3

9,50332Е-05