Large surface deformation due to thermo-mechanical effects during cryopreservation by vitrification - mathematical model and experimental validation

PLoS One. 2023 Mar 9;18(3):e0282613. doi: 10.1371/journal.pone.0282613. eCollection 2023.

Abstract

This study presents a simplified thermal-fluids (TF) mathematical model to analyze large surface deformations in cryoprotective agents (CPA) during cryopreservation by vitrification. The CPA deforms during vitrification due to material flow caused by the combined effects of thermal gradients within the domain, thermal contraction due to temperature, and exponential increase in the viscosity of the CPA as it is cooled towards glass transition. While it is well understood that vitrification is associated with thermo-mechanical stress, which might lead to structural damage, those large deformations can lead to stress concentration, further intensifying the probability to structural failure. The results of the TF model are experimentally validated by means of cryomacroscopy on a cuvette containing 7.05M dimethyl sulfoxide (DMSO) as a representative CPA. The TF model presented in this study is a simplified version of a previously presented thermo-mechanics (TM) model, where the TM model is set to solve the coupled heat transfer, fluid mechanics and solid mechanics problems, while the TF model omits further deformations in the solid state. It is demonstrated in this study that the TF model alone is sufficient to capture large-body deformations during vitrification. However, the TF model alone cannot be used to estimate mechanical stresses, which become significant only when the deformation rates become so small that the deformed body practically behaves as an amorphous solid. This study demonstrates the high sensitivity of deformation predictions to variation in material properties, chief among which are the variations of density and viscosity with temperature. Finally, this study includes a discussion on the possibility of turning on and off the TF and TM models in respective parts of the domain, in order to solve the multiphysics problem in a computationally cost-effective manner.

Publication types

  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Cryopreservation* / methods
  • Cryoprotective Agents / chemistry
  • Cryoprotective Agents / pharmacology
  • Dimethyl Sulfoxide / pharmacology
  • Models, Theoretical
  • Vitrification*

Substances

  • Cryoprotective Agents
  • Dimethyl Sulfoxide

Grants and funding

Research reported in this paper was supported by the National Science Foundation (NSF) award number EEC 1941543. The NSF had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The content of this paper is solely the responsibility of the authors and does not necessarily represent the official views of the NSF.