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

AW40

fx = ЕСЛИ(AW38>$C$19;”                                “$D$33)

AV

AW

AX

AY

AZ

BA

37

38

11

11

12

12

12

13

39

Ti

hi

qi+1

Ti

hi

qi+1

40

277

2800

0

41

277,00

2800

0

42

277,00

2800

0

Рис.5.13. Определение начальных значений энтальпий для всех точек моделируемого массива

После того как все ячейки строки τ=0 заполнены, вычисления переходят на следующий временной слой. Первым определяется значение энтальпии для слоя i=1.

AV

AW

AX

AY

AZ

BA

37

38

11

11

12

12

12

13

39

Ti

hi

qi+1

Ti

hi

qi+1

40

=ЕСЛИ($K40=1;ЕСЛИ(AW$38<$C$19;(AY40-AV40)*

СМЕЩ($B$27;0;AW$38);ЕСЛИ(AW$38=$C$19;0;”        “));”

41

Рис.5.14. Определение значений энтальпии i-та

Сумма тепловых потоков на границе элемента умножается на шаг по времени и делится на массу текущего элемента, полученный результат имеет отрицательный знак. При его сложении с предыдущим значением энтальпии получается новое значение hi для слоя τ=Δτ. Полученное значение hiпозволяет определить новое значение Ti. Возможные три состояния элемента: талое, замерзающее и мёрзлое. В состоянии замерзания r*< hi<0, тогда температура равна температуре дефростатции Ti =270 К. Температура для мёрзлого и талого состояния рассчитывается. Формула для определения температуры содержит пять операторов "Если".

P

Q

R

S

T

U

38

0                                                                                                       

1

1

1

2

2

39

qi+1

Ti

hi

qi+1

Ti

40

200

0

277

2800

0

277

41

=ЕСЛИ(S$38<$C$19+1;ЕСЛИ($K41=1;S40+(Q40+T40)*$B$38/СМЕЩ($C$26;0;S$38);”       “);”             “)

42

44

ЕСЛИ(лог_выражение;[значение_если_истина];[значение_если_ложь])

 

Рис.5.15. Определение температуры элемента на новом временном слое

Для вычисления энтальпии для элемента i=2 формула корректируется, так как для всех участков кроме i=1, значения qi-1 не вычисляется, так как уже рассчитано при составлении теплового баланса точки i-1. Для второго элемента (см. рис. 5.16.) тепловой поток к предстоящему элементу численно равен значению, содержащемуся в ячейки Т41, но учтённому с отрицательным знаком. В остальном формулы идентичны.

T

U

V

W

X

Y

38

2

2

2

3

3

3

39

qi+1

Ti

hi

qi+1

Ti

hi

40

0

277

28000

0

277

28000

41

0

277,00

28000

0

277,00

28000

42

=ЕСЛИ(V$38<$C$19+1;ЕСЛИ($K42=1;V41+(-T41+W41)*$B$38/ СМЕЩ($C$26;0;V$38);”       “);”             “)

43

44

ЕСЛИ(лог_выражение;[значение_если_истина];[значение_если_ложь])