Biphasic investigation of tissue mechanical response during freezing front propagation

Cryopreservation of engineered tissue (ET) has achieved limited success due to limited understanding of freezing-induced biophysical phenomena in ETs, especially fluid-matrix interaction within ETs. To further our understanding of the freezing-induced fluid-matrix interaction, we have developed a biphasic model formulation that simulates the transient heat transfer and volumetric expansion during freezing, its resulting fluid movement in the ET, elastic deformation of the solid matrix, and the corresponding pressure redistribution within. Treated as a biphasic material, the ET consists of a porous solid matrix fully saturated with interstitial fluid. Temperature-dependent material properties were employed, and phase change was included by incorporating the latent heat of phase change into an effective specific heat term. Model-predicted temperature distribution, the location of the moving freezing front, and the ET deformation rates through the time course compare reasonably well with experiments reported previously. Results from our theoretical model show that behind the marching freezing front, the ET undergoes expansion due to phase change of its fluid contents. It compresses the region preceding the freezing front leading to its fluid expulsion and reduced regional fluid volume fractions. The expelled fluid is forced forward and upward into the region further ahead of the compression zone causing a secondary expansion zone, which then compresses the region further downstream with much reduced intensity. Overall, it forms an alternating expansion-compression pattern, which moves with the marching freezing front. The present biphasic model helps us to gain insights into some facets of the freezing process and cryopreservation treatment that could not be gleaned experimentally. Its resulting understanding will ultimately be useful to design and improve cryopreservation protocols for ETs.