Numerical Study of Heat and Mass Transfer during Cryopreservation Process with Application of Directed Interval Arithmetic

In the present paper, numerical modelling of heat and mass transfer proceeding in a two-dimensional axially symmetrical articular cartilage sample subjected to a cryopreservation process is presented. In the model under consideration, interval parameters were assumed. The heat transfer process is described using the Fourier interval equation, while the cryoprotectant transport (DMSO) across the cell membrane is analyzed using a two-parameter model taking into account the simulation of the water volume in the chondrocytes and the change in DMSO concentration over time. The liquidus tracking (LT) protocol introduced by Pegg et al. was used to model the cryopreservation process. This procedure divides the heating and cooling phases into eight and seven steps, respectively, allowing precise regulation of temperature and cryoprotectant (CPA) concentration of bathing solutions. This protocol protects chondrocytes from ice crystal, osmotic stress, and electrolyte damage. The obtained interval concentrations of cryoprotectant in chondrocytes were compared with previous simulations obtained using the deterministic model and they are mostly in agreement with the simulation data.