Prediction of local cooling rates and cell survival during the freezing of a cylindrical specimen

A finite element numerical model was implemented to simulate the freezing process of an aqueous salt solution in a cylindrical container. Local cooling rates within the container were computed for several defined cooling protocols applied at the boundary. Characteristic cell survival signatures were used to predict the associated local survival rates throughout the system. These calculations show that there are two definite time domains during a typical freezing process: (1) while the surface temperature is changing and (2) after the surface temperature reaches a constant storage value. The calculations also show significant spatial variations in the local cooling rates within the container and considerable local deviation from the volumetric average survival for various simulated freezing protocols.