1、第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,2023文章编号:0253-4339(2023)05-0160-07doi:10.3969/j.issn.0253-4339.2023.05.160固体表面玻璃化冷冻方法降温结晶特性数值研究朱文欣1 赵芷慧1 黄永华1 程锦生2 陈 威3 李 铮3(1 上海交通大学机械与动力工程学院 上海 200240;2 上海东富龙医疗装备有限公司 上海 201108;3 上海交通大学附属第一人民医院泌尿中心男科 上海 200080)摘 要 固体表面玻璃化冷冻方法
2、(SSV)及其变式在细胞冻存领域的应用越来越广泛,但多基于实验和经验,缺乏仿真分析和物理机理解释,影响相关技术的通用化设计和研发效率。本文以质量分数 35%的丙二醇溶液为对象,考虑冷冻过程的非等温结晶特性,建立了耦合非等温结晶源项的 SSV 方法传热模型,研究了预冷温度、样品体积和形状对实施 SSV 方法的液体温度与结晶特性的影响。结果表明:该模型可以有效揭示冷冻过程中液滴样品内部的温度及结晶度变化规律。样品结晶危险区位于远离低温表面处,且最易结晶点不一定位于液滴顶点。预冷温度越低,冷冻样品体积越小,形状越扁平,样品降温速率越快,越容易实现玻璃化冷冻。关键词 细胞冻存;固体表面玻璃化冷冻;预冷
3、;快速降温;非等温结晶中图分类号:TB61+1;TQ025.3;R318.52文献标识码:ANumerical Study on Cooling and Crystallization Characteristics of Solid Surface Vitrification MethodZhu Wenxin1 Zhao Zhihui1 Huang Yonghua1 Cheng Jinsheng2 Chen Wei3 Li Zheng3(1.School of Mechanical Engineering,Shanghai Jiao Tong University,200240,China;
4、2.Shanghai Tofflon Medical Equipment Co.,Ltd.,Shanghai,201108,China;3.Department of Andrology,Urologic Medical Center,Shanghai General Hospital,Shanghai Jiao Tong University,Shanghai,200080,China)Abstract Solid surface vitrification(SSV)has been increasingly applied in the field of cell cryopreserva
5、tion.However,most studies based on empirical experimentations have not included simulation analysis and physical mechanism explanation,consequently affecting the universality of the designs and innovation efficiency in SSV-related technologies.This study focused on the non-isothermal crystallization
6、 characteristics and established a thermal diffusion model of SSV coupled with a non-isothermal crystallization source.A propylene glycol solution with the mass fraction of 35%was selected as the sample for analysis.Additionally,the effects of the pre-cooling temperature,sample volume,and sample sha
7、pe on the temperature and crystallization characteristics of the SSV method were investigated.The results show that the established model can directly reveal change patterns of the temperature and the crystallinity inside the sample during the freezing process.The region most likely to undergo cryst
8、allization is located at the upper part of the sample away from the precooled surface,and the most crystallizable region is not necessarily located at the top of the droplet sample.Furthermore,vitrification of the sample is easier to achieve because of the lower precooling temperature,smaller sample
9、 volume,and flatter shape.Keywords cryopreservation;solid surface vitrification;precool;rapid cooling;non-isothermal crystallization基金项目:国家重点研发计划(2022YFC2703000)资助项目。(The project was supported by the National Key Research and Development Program of China(No.2022YFC2703000).)收稿日期:2023-02-05;修回日期:2023
10、-03-22 细胞疗法被认为是未来医学的三大支柱之一,生物样本将成为重要战略资源1。由于生物的活性在低温环境中受到抑制,代谢等化学反应减慢近乎停滞2-3,使细胞等生物长期冻存成为可能4。然而,冷冻过程温跨大,造成细胞等生物在冷冻过程中产生不可逆损伤,如氧化性损伤、渗透性损伤和机械损伤等5-6。为减少冷冻损伤,多种冷冻方法在细胞等生物样本冻存领域获得应用7,如直接接触液氮冷却法8、环网状冷冻法9、微流控冷冻法10等。然而上述方法均易造成样品受热不均、回收困难并存在样品污染的风险。而固体表面玻璃化冷冻法(solid surface vitrification,SSV)由于降温速率快、降温速率控制性
11、好、样品不直接与液氮接触、易回收、无污染等优点,061第 44 卷 第 5 期2023 年 10 月固体表面玻璃化冷冻方法降温结晶特性数值研究Vol.44,No.5October,2023被越来越多地应用于低温生物冻存11-12。目前,针对 SSV 方法的研究多数为实验或经验性报道13-14,主要聚焦于降温速率的大小,对于 SSV 方法在冻存过程中造成的冷冻样品结晶机理与结晶度变化规律鲜有涉及。此外,由于冻存过程时间短,通过实验手段以高时间分辨率获得全冷冻过程样品全域温度、结晶度变化仍较为困难,如采用热电偶测温不仅侵入性大、易结露、漏热严重,而且无法测量结晶度;低温显微系统价格昂贵且无法定量分
12、析冷冻过程的结晶变化15等。鉴于此,本文采用数值计算方法探究 SSV 方法实现快速降温的核心要素,并通过建立冻存过程的非等温结晶预测模型,研究预冷温度、样品体积及形状对冷冻结晶效果的影响,为冷冻方法优化、冷冻对象体积形状等设计提供理论指导。1 基于非等温结晶模型的 SSV 方法建模 SSV 方法采用一定质量的高导热金属(铜、银或铝合金)块或片浸入液氮或其他低温介质预冷,将样品滴在预冷后的低温表面,通过瞬间接触导热实现快速降温冷冻,物理过程如图 1 所示。图 1 SSV 方法工作原理Fig.1 Process of the SSV method聚焦液滴碰触金属表面之后的传热过程,建立如图 2 所
13、示的数学物理模型。鉴于轴对称性,该模型采用二维轴对称形式,取一半,其中左上角凸出部分为液滴,其余长方形区域为金属基底载冷台。以 5 L液体冷冻为例,真实滴落液滴的几何尺寸可以近似为半径 1.82 mm,中心高度为 0.89 mm 的球冠,液滴的底部与金属块的上表面直接接触。此处考虑无氧铜材质的基底载冷台,其导热系数在低温下达 400 W/(m K)。为了消除边缘效应,取铜块的宽度为 50 mm,高度为 20 mm,此时相对液滴而言,可以认为铜块无限大。基于傅里叶导热定律,金属块内的导热可表示为:cpTt=(T)(1)式中:T 为温度,K;为密度,kg/m3;为导热系数,图 2 数学物理模型和边
14、界条件Fig.2 Mathematical physical model and boundary conditionsW/(mK);cp为比定压热容,J/(kgK);t 为时间,s。求解凝固问题的区域模型理论认为:液体的结晶过程发生在某一区域内,固液区间无清晰边界。因此,通过在能量方程中添加热源项,来描述结晶对传热过程的影响,将传热与结晶动力学相耦合。由于实际操作过程中样品体积小(一般为微升量级),传热时间短,忽略液体内部对流换热,则冷冻液滴内部含内热源 的热扩散方程为:cpTt=2T+(2)=ddt=LdXdt(3)式中:L 为熔融潜热,kJ/mol;X 为结晶度。忽略结晶界面处由于分子扩
15、散对冰晶前沿生长速率的影响,根据半经验非等温结晶模型16-20:dXdt=k1X2/3(1-X)(Tm-T)exp(-QRT)(4)式中:Tm为熔融温度,K;k1为结晶常数,与液体种类、质量浓度相关;Q 为活化能,kJ/mol;R 为气体常数 8.314 J/(molK)。由于实际采用的抛光铜表面粗糙度 Ra 为 0.8,视作光滑表面条件,且该铜表面与液滴直接接触,因此在进行建模时忽略了金属块与液滴样品间的接触热阻:TS=TC(5)Tn()S=Tn()C(6)式中:Tn为物体沿 n 方向的温度变化率;下标 C 表示铜块载体区域,S 表示被冷冻液滴样品。在边界条件方面,将金属块与液氮接触边界设为
16、第一类边界条件,由于金属块与液氮蒸气对流换热,因此其他边界设为第三类边界条件,对流换热表面传热系数 h 取为 100 W/(m2 K)21。传热过程的计算选取液滴假设为从低处自由落体并以规则的半球形态滴落在铜表面的瞬间为初始161第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,2023时刻。设置铜块初始温度已预冷至液氮温度 TC,0=78 K;样品液滴初始温度 TS,0=293.15 K。初始结晶度 X0=0,初始结晶度对时间的变化率dXdt=0。在热物性方面,铜块选用常物性参数,如表1 所示。表
17、 1 固体材料热物性参数Tab.1 Thermophysical properties of the copper材料导热系数/W/(m K)密度/(kg/m3)比定压热容/J/(kg K)铜4008 960385 需要说明的是,该模型侧重于了解被冷冻液滴所经历的传热过程,并实现与系统的传质过程解耦,其分析尺度大于冰晶生长的微结构,且不包含溶液传质浓度变化对冰晶微结构的影响。在 COMSOL 软件中,模型划分自由三角形网格,在样品与预冷表面接触处做细化处理,共划分了1 500、2 500、3 035、4 000 四种密度的网格,网格平均质量均优于 0.93。经网格无关性验证,取 3 035 个
18、网格进行计算在速度和精度方面均满足要求。对金属块用固体传热接口求解,对 CPA 液滴利用导热接口添加广义热源项,通过域常微分和微分代数方程接口定义非等温结晶源项,并耦合求解。添加温度判据以保证结晶在熔融温度下才发生。为保证对极其微小量数值计算正常,将结晶度初始值设为 110-16(可识别的接近 0的最小正值)。采用瞬态求解器,向后差分方法,时间步长不超过 0.01 s,相对容差控制为 110-6。图 3 模型计算值与文献22实验值对比Fig.3 Comparison of model calculated values with experimental values from the lit
19、erature22采用文献22的实测数据来验证模型的适用性,针对西林瓶(类薄层式 SSV 结构)中 1.9 mL 甘油水溶液置于液氮中冷冻过程温度变化进行对比验证。需要说明的是,验证时模型中固体材料选用玻璃,溶液选用甘油水溶液,与文献22的实验保持一致。预测结果和文献数据的对比如图 3 所示。由图 3 可知,耦合非等温结晶源项的传热区域模型可以较好地预测冷冻过程,最大误差约为 10%,平均误差小于2%。后期温度略高于实验值,可能是由于被测样品溶液的导热系数在低温时有所升高。2 计算结果与分析 模型验证采用的实验样本液滴为生物细胞冷冻操作常用的冷冻保护剂(CPA)丙二醇溶液(质量分数为 35%)
20、,其玻璃化转变温度 Tg为 160 K,冷冻结晶的危险温区介于 Tg与熔化温度 Tm之间(160 255.5 K)。由于 CPA 物性参数通常缺乏低温区的实测数据,采用常温区参数一般会低估物性参数的变化;而若采用水和冰的物性参数替代低温保护剂物性参数,会高估低温区物性参数的变化。在没有准确合适的物性参数数据时,有必要假设物性参数不随温度变化来参考计算23。一般而言,随着温度的降低,CPA 导热系数增加,而比热容和密度降低,总体上热扩散系数增加。因此利用常物性参数计算可能会低估冷却速率,即实际降温效果优于计算值。样品液滴(丙二醇溶液)物性参数来源于文献24-25,如表 2所示。利用上述模型和物性
21、参数,开展被冷冻样品滴落后的冷却和结晶特性研究。表 2 冷冻液滴热物性参数Tab.2 Thermophysical properties of the sample参数数值参数数值/W/(m K)0.23Q/(kJ/mol)51.1/(kg/m3)1 012.7L/(kJ/kg)289.9cp/J/(kg K)2 276k1/(min K)-14.561010 Tm/K255.52.1 非等温结晶特性 首先从液滴对称轴的高、中、低三个位点进行局部降温结晶特性分析,如图 4 所示。三个位点的冷冻降温过程均未呈现温度平台期,无明显的一级相变特征。由于低点最接近预冷表面,传热路径短、热阻小,因此液滴
22、温度在低点处最接近预冷表面温度,降温速率最快,相比于液滴中点及高点的降温速率有数量级优势。不同的降温特性导致液滴起始结晶时间和最终结晶度不同。根据 Uhlmann 理论,结晶度小于 10-6即为玻璃化26,低点的结晶度几乎不变,趋向初始设定值 10-16量级,实现了玻璃化转变。中点降温速度比高点更快,先于高点降温到熔融温度,所以结晶比高点更早产生。同时,由于中点比高点的传热路径更小,传热速率更大,能较快通过结晶危险温区,最终结晶度小于高点。只有以较高的温变速率快速通过结261第 44 卷 第 5 期2023 年 10 月固体表面玻璃化冷冻方法降温结晶特性数值研究Vol.44,No.5Octob
23、er,2023晶危险温区,才能实现被冷冻样品的玻璃化转变。对液滴样品选取切面进行整体降温结晶特性分析,如图 5 所示(注:图中只针对实现玻璃化冷冻及处于危险温区的部分显示云图)。液滴的降温最先从底部发生并实现玻璃化转变。样品中间结晶度增大,并向上扩散结晶,由于液滴边缘薄中间厚,结晶度与等温线均呈上凹形。但最大结晶度未出现在距离冷源最远、温度最高的液滴顶端,而出现于距顶端一定距离的中部。需要说明的是,由于本文结晶模型与传质解耦,因此,液滴中部存在最大结晶度这一现象与冷冻过程中样品界面处存在液气-固气相变无关,今后可进一步分析液滴内部复杂的温度变化情况对样品结晶度的影响。图 4 样品内低、中、高点
24、处温度与结晶度变化规律Fig.4 Variation of the temperature and the crystallinity at low,middle and high points within the sample图 5 液滴内部结晶度和温度云图Fig.5 The crystallinity and the temperature distribution inside the droplet图 6 液滴顶点、结晶度最大处的降温速率与结晶度变化速率Fig.6 The cooling rate and the crystallinity change rate at the ma
25、ximum crystallinity and the top of the droplet 进一步对冷冻过程中的液滴顶点处与结晶度最大点处开展降温速率和结晶度变化速率特性分析,如图 6 所示。顶点的降温速率和最大结晶处的降温速率趋势一致,但被冷冻液滴内部的降温速度存在不均匀性,不利于大体积样品的均匀冷冻。最大结晶度变化速率点比最大降温速率点滞后约 0.3 s,表明结晶过程存在滞后性。在整个液滴均出现结晶度后,顶点的降温速率始终大于最大结晶点处的降温速率,顶点的液体更快通过结晶危险温区,表明更大的降温速率可以抑制结晶,使被冷冻样品更容易实现玻璃化。需要说明的是,上文通过引入冷冻样品的结晶度X,
26、更直观地展示了样品在被冷冻过程中的内部状态,可预测样品的被冷冻效果,为设法实现无损冻存的各种尝试提供便捷高效的验证方法。2.2 固体表面玻璃化冷冻影响因素 针对不同的生物样本,通常需要设定不同的冷却预冷温度。因此,有必要探索预冷温度对样品在固体表面被玻璃化冷冻的影响。液氮是冷却金属载物台最常见的工质,其预冷极限温度 Tprec,min为 78 K,而35%丙二醇溶液的 Tm为 255.5 K,Tg为 160 K。因361第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,2023此,分别选取 Tpre
27、c为 80、100、120、140、160、200、250 K。图 7 所示为预冷温度对液滴结晶度和穿越危险温区时降温最慢处的平均降温速率 vcr的影响。预冷温度越低,降温速率越慢,结晶度越小,预冷温度与降温速率和结晶度均无线性关系。此外,当预冷温度低于样品玻璃化转变温度时,样品的降温速率均较快,达到 103 K/min 量级,其结晶度在 10-5量级;而当预冷温度高于 Tg时,降温速率在 102 K/min 量级,结晶度在 10-510-4量级;当预冷温度越接近熔融温度,降温速率越慢,极易结晶。因此,预冷温度越低,传热温差越大,越有利于实现样品的玻璃化冷冻。图 7 预冷温度对降温速率和结晶度
28、的影响Fig.7 Effect of pre-cooling temperature on the cooling rate and the crystallinity图 8 不同体积样品中结晶度和玻璃化转变率Fig.8 Crystallinity and glass transition rate in different volume samples前文非等温结晶模型计算表明,最大结晶处可能在被冷冻样品的中上部。因此,有必要探究样品体积对用 SSV 方法冷冻的样本的降温速率和结晶度的影响。基于现有实际应用和文献调研情况,样品体积 V分别选取 1、5、10、20、50、100、500 L,定
29、义玻璃化比例 VR 为样品中实现玻璃化冷冻体积与样品总体积之比,计算结果如图 8 所示。当体积为 1 L 时,样品完全玻璃化冷冻,因此,完全玻璃化转变的安全样品体积为 1 L。随着样本体积增加,样本内部导热路径变长,降温速率变慢,结晶显著增加,VR 不断下降,且呈非线性关系。在一定范围内,体积越大,结晶度变化越剧烈;而体积过大,结晶度将完全结晶。因此,样本体积是影响 SSV 方法的关键因素之一。本文提出的耦合热源项的非等温结晶模型可以预测不同体积的样品内部结晶程度和不同条件下玻璃化安全体积,为实际固体表面玻璃化冷冻设计提供一定依据。此外,根据傅里叶导热定律和牛顿冷却公式,有效传热面积是研究传热
30、问题的关键因素之一。因此,引入传热面积与总体积比值 A/V 参数,研究样品形状对SSV 方法冷冻的影响。样品体积为 5 L,改变样品形状,如图 9 所示。正球缺即为图 2 介绍的模型;反球缺与正球缺尺寸相同倒置在挖空的金属预冷表面;反半球形样品半径为 1.34 mm;反圆柱形半径为 1.51 mm,厚 0.7 mm;反薄层形半径为 2.3 mm,厚 0.3 mm。图 9 样品形状Fig.9 The sample shape不同形状液滴结晶度和玻璃化比例的计算结果如图 10 所示。由图 10 可知,相同体积下薄层形样品的 A/V 更大,玻璃化程度更高,而其他形状样品出现部分结晶。因此,A/V 一
31、定程度上可以反映传热拓扑结构的优劣。相同情况下,A/V 值越大,结晶度越小,样品越容易实现玻璃化冷冻。若以实现玻璃化冷冻为目标,样品形状应尽量薄层化。说明本文模型可以简便高效地研究不同形状样品在冷冻时内部结晶情况,可为实际冷冻样品形状设计提供依据。图 10 不同形状液滴结晶度和玻璃化比例Fig.10 Crystallinity and vitrification ratio of droplets with different shapes3 结论 本文对基于固体表面玻璃化冷冻方法的冷冻液滴建立了耦合非等温结晶源项的传热区域模型,研究了液滴降温特性和冷冻液滴内部的结晶情况,得到如下结论:461
32、第 44 卷 第 5 期2023 年 10 月固体表面玻璃化冷冻方法降温结晶特性数值研究Vol.44,No.5October,20231)固体表面玻璃化冷冻方法可以实现快速降温,降温速率可达到 103 K/min 量级。其可获得快速降温能力的核心要素是固体冷台得到了充分预冷并储存足够冷量。2)降温速率越快,结晶程度越小。结晶易出现在远离预冷表面的样品上部,结晶最严重区域不一定位于样品顶端。3)预冷温度越低、样品体积越小、样品面积/体积(即薄层化)越大,结晶程度越小,越容易实现玻璃化冷冻。4)耦合非等温结晶源项的传热区域模型可以通过结晶度 X 反映冷冻过程中样品内部结构变化,可作为设计开发新型冷
33、冻设备和低温保护剂的辅助手段。本文受上海市科委创新行动计划(20S31903400);上海交通大学“交大之星”计划医工交叉研究基金(YG2021QN86)项目资助。(The project was supported by Innovation Action Plan of Shanghai Science and Technology Commission(No.20S31903400),and Shanghai Jiao Tong University Med-Engineering Fund(No.YG2021QN86).)参考文献1 FISCHBACH M A,BLUESTONE J
34、A,LIM W A.Cell-based therapeutics:the next pillar of medicine J.Science Translational Medicine,2013,5(179):179ps7.2 华泽钊,任禾盛.低温生物医学技术M.北京:科学出版社,1994.(HUA Zezhao,REN Hesheng.Cryogenic biomedical technology M.Beijing:Science Press,1994.)3 MAZUR P.Freezing of living cells:mechanisms and implications J.A
35、merican Journal of Physiology-Cell Physiology,1984,247(3):C125-C142.4 POLGE C,SMITH A U,PARKES A S.Revival of spermatozoa after vitrification and dehydration at low temperaturesJ.Nature,1949,164(4172):666.5 VAJTA G,NAGY Z P.Are programmable freezers still needed in the embryo laboratory?Review on vi
36、trificationJ.Reproductive BioMedicine Online,2006,12(6):779-796.6 MAZUR P,LEIBO S P,CHU E H.A two-factor hypothesis of freezing injury:evidence from Chinese hamster tissue-culture cells J.Experimental Cell Research,1972,71(2):345-355.7 赵芷慧,黄永华.细胞冻融策略优化及设备研究进展J.制冷学报,2022,43(2):13-22.(ZHAO Zhihui,HUAN
37、G Yonghua.Advances on optimization strategies and devices for cell cryopreservation J.Journal of Refrigeration,2022,43(2):13-22.)8 LANDA V,TEPL O.Cryopreservation of mouse 8-cell embryos in microdropsJ.Folia Biologica,1990,36(3/4):153-158.9 TIERSCH C J,LIU Yue,TIERSCH T R,et al.3-D printed customiza
38、ble vitrification devices for preservation of genetic resources of aquatic species J.Aquacultural Engineering,2020,90:102097.10 周新丽,彭湉,王彪.PDMS-玻璃微流控芯片上附着态HepG2 细胞低温保存条件研究J.制冷学报,2020,41(5):159-166.(ZHOU Xinli,PENG Tian,WANG Biao.Research on cryopreservation conditions of attached HepG2 cells on PDMS-
39、glass microfluidic chipJ.Journal of Refrigeration,2020,41(5):159-166.)11 DOU R,SAUNDERS R E,MOHAMET L,et al.High throughput cryopreservation of cells by rapid freezing of subul drops using inkjet printing-cryoprinting J.Lab on a Chip,2015,15(17):3503-3513.12 EL ASSAL R,GUVEN S,GURKAN U A,et al.Bioin
40、spired cryoink preserves red blood cell phenotype and function during nanoliter vitrification J.Advanced Materials,2014,26(33):5815-5822.13 AKIYAMA Y,SHINOSE M,WATANABE H,et al.Cryoprotectant-free cryopreservation of mammalian cells by superflash freezing J.Proceedings of the National Academy of Sci
41、ences of the United States of America,2019,116(16):7738-7743.14 PRUSS D,YANG H,LUO X,et al.High-throughput droplet vitrification of stallion sperm using permeating cryoprotective agents J.Cryobiology,2021,101:67-77.15 DILLER K R,CRAVALHO E G.A cryomicroscopic investigation of intracellular ice forma
42、tion in frozen erythrocytesJ.Cryobiology,1971,8(4):398.16 BOUTRON P.Comparison with the theory of the kinetics and extent of ice crystallization and of the glass-forming tendency in aqueous cryoprotective solutions J.Cryobiology,1986,23(1):88-102.17 BAUDOT A,ALGER L,BOUTRON P.Glass-forming tendency
43、in the system water-dimethyl sulfoxide J.Cryobiology,2000,40(2):151-158.18 BOUTRON P,KAUFMANN A.Stability of the amorphous state in the system water-glycerol-ethylene glycol J.Cryobiology,1979,16(1):83-89.19 BOUTRON P,KAUFMANN A.Stability of the amorphous state in the system water-glycerol-dimethyls
44、ulfoxide J.Cryobiology,1978,15(1):93-108.20 BOUTRON P,MEHL P.Theoretical prediction of devitrification tendency:determination of critical warming 561第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,2023rates without using finite expansions J.Cryobiology,1990,27(4):359-377.21
45、ZHAN L,GUO S Z,KANGAS J,et al.Conduction cooling and plasmonic heating dramatically increase droplet vitrification volumes for cell cryopreservation J.Advanced Science,2021,8(11):2004605.22 马青.冻干保存中相变传热过程的数值仿真研究D.成都:电子科技大学,2018.(MA Qing.Simulation of the phase change heat transfer during freeze dryi
46、ngD.Chengdu:University of Electronic Science and Technology of China,2018.)23 OSTLER T,WOOLLEY T E,SWANN K,et al.Vitrifying multiple embryos in different arrangements does not alter the cooling rateJ.Cryobiology,2021,103:22-31.24 CHOI J,BISCHOF J C.Review of biomaterial thermal property measurements
47、 in the cryogenic regime and their use for prediction of equilibrium and non-equilibrium freezing applications in cryobiology J.Cryobiology,2010,60(1):52-70.25 任禾盛,华泽钊,俞国新,等.低温抗冻剂溶液玻璃化转变的临界冷却速率J.上海机械学院学报,1990,12(3):23-30.(REN Hesheng,HUA Zezhao,YU Guoxin,et al.The critical cooling rates for vitrific
48、ation of cryoprotective solutions J.Journal of University of Shanghai for Science and Technology,1990,12(3):23-30.)26 UHLMANN D R.A kinetic treatment of glass formationJ.Journal of Non-Crystalline Solids,1972,7(4):337-348.通信作者简介黄永华,男,研究员,博士生导师,上海交通大学机械与动力工程学院,021-34206295,E-mail:huangyh 。研究方向:低温绝热与传
49、热、低温流体与材料热物性质、低温冷冻保存。About the corresponding authorHuang Yonghua,male,researcher,Ph.D.supervisor,School of Mechanical Engineering,Shanghai Jiao Tong University,86-21-34206295,E-mail:huangyh .Research fields:cryogenic insulation and heat transfer,thermal physical properties of cryogenic fluids and ma
50、terials,cryo-preservation.(上接第 144 页)通信作者简介孙志高,男,博士,教授,苏州科技大学环境科学与工程学院,0512-63176109,E-mail:szg.yzu 。研究方向:储能与节能技术。About the corresponding authorSun Zhigao,male,Ph.D.,professor,School of Environmental Science and Engineering,Suzhou University of Science and Technology,86-512-63176109,E-mail:szg.yzu .
