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

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

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

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

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

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

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