The purpose of this work is to demonstrate the existing possibilities and methodology of mathematical modeling of the freezing phase in the technology of lyophilization of vials filled with a solution containing a pharmaceutical substance and auxiliary substances. The freezing process is crucial for the subsequent stages of primary sublimation and the drying stage. So it is at the stage of freezing that ice crystals are formed and a certain microstructure forming the pores for subsequent stages and influencing the speed of all stages in the future. The developed methods allow estimating the kinetics of cooling, freezing, calculating the sizes of ice crystals and determining the permeability of the sublimation layer. Of course, for these calculations it is necessary to know the conditions and values of the variables, calculated and measured during the pilot freezing. However, modeling allows to reduce the time of development and optimization of the process, increases the speed of transfer of the technological process to other equipment. In this review, formulas are analyzed, the thermal conductivity equations for each cooling zone, equations for estimating the crystal size and pore size in the dried layer, and calculating the permeability of the sublimated layer. The article draws conclusions about the perspective of modeling methods for the freezing phase, examines the equations used in modeling and demonstrates the model of supercooling, crystal formation and other mass and heat exchange processes during the cooling and freezing stage. The calculations presented in this paper are confirmed by references to experimental data and have great practical value.
Цель работы – продемонстрировать существующие возможности и методику математического моделирования этапа заморозки в технологии лиофилизации флаконов, наполненных раствором, содержащим фармацевтическую субстанцию и вспомогательные вещества. Процесс заморозки является определяющим для последующих этапов первичной сублимации и этапа досушки. Так, именно на этапе замораживания формируются кристаллы льда и определенная микроструктура, формирующая поры для последующих этапов и влияющая в дальнейшем на скорость всех этапов. Разработанные методики позволяют оценить кинетику охлаждения, замерзания, просчитать размеры ледяных кристаллов и определить проницаемость сублимационного слоя. Для этих расчетов необходимо знать условия и значения переменных, рассчитываемые и измеряемые во время пилотной заморозки. Однако моделирование позволяет сократить время разработки и оптимизации процесса, повышает скорость переноса технологического процесса на другое оборудование. В данном обзоре разобраны формулы, уравнения теплопроводности для каждой зоны охлаждения, уравнения для оценки размеров кристаллов и размера пор в высушенном слое, а также для расчета проницаемости сублимированного слоя. В статье сделаны выводы о перспективе методов моделирования этапа заморозки, рассмотрены уравнения, используемые в моделировании, и наглядно продемонстрированы моде ли переохлаждения, кристаллообразования и других массо- и теплообменных процессов во время этапа охлаждения и замораживания. Расчеты, представленные в этой работе, подтверждаются ссылками на экспериментальные данные и имеют большую практическую ценность.