金属学报  2014 , 50 (11): 1403-1412 https://doi.org/10.11900/0412.1961.2014.00300

定向凝固藕状多孔Al生长过程的模拟仿真

杨倩倩1, 刘源12, 李言祥12

1 清华大学材料学院, 北京 100084
2 清华大学先进成形制造教育部重点实验室, 北京 100084

MODELING AND SIMULATION OF STRUCTURAL FORMATION OF POROUS ALUMINUM IN GASAR SOLIDIFICATION

YANG Qianqian1, LIU Yuan12, LI Yanxiang12

1 School of Materials Science and Engineering, Tsinghua University, Beijing 100084
2 Key Laboratory for Advanced Materials Processing Technology (Ministry of Education), Tsinghua University, Beijing 100084

中图分类号:  TG146, TG249

通讯作者:  Correspondent: LIU Yuan, associate professor, Tel: (010)62789328, E-mail: yuanliu@tsinghua.edu.cn

收稿日期: 2014-09-4

修回日期:  2014-06-9

网络出版日期:  --

版权声明:  2014 《金属学报》编辑部 版权所有 2014, 金属学报编辑部。使用时,请务必标明出处。

基金资助:  * 国家自然科学基金项目51271096和教育部新世纪优秀人才支持计划项目NCET-12-0310资助

作者简介:

杨倩倩, 女, 1986年生, 博士生

展开

摘要

通过对Gasar凝固过程中的传质、气泡形核、气孔生长、中断及脱离等的理论分析, 建立了一个描述单气孔演变过程的非稳态三维模型, 并采用有限差分的方法模拟了不同凝固速率下定向凝固多孔Al的气孔形貌. 基于Al-H2系的研究结果表明: 当凝固速率在0.15~0.005 mm/s范围内时, 固/气两相能够维持协同生长, 气孔的平均孔径分布在100~1100 mm之间, 但随凝固速率的降低会逐渐增加, 气孔长度亦随凝固速率的降低逐渐增加, 长径比则基本保持在40左右; 当凝固速率为0.015 mm/s时, 气孔孔径的模拟值与实验值吻合较好, 之后随凝固速率的降低, 模拟孔径略低于实验值, 分析认为, 实际凝固过程中熔体上方的H2向熔体内的不断扩散是导致该差异的主要原因; 随着熔体过热度和H2分压的逐渐增大, 对应藕状多孔Al固/气两相协同生长凝固速率范围的最大值由不足0.01 mm/s先逐渐增加之后稳定在0.15 mm/s, 最小值则由0.0001 mm/s左右逐渐增加至约0.01 mm/s; 对Al-H2系和Cu-H2系相关参数的比较分析表明, H2在金属熔体中的溶解度是决定Al-H2系Gasar结构中固/气两相协同生长凝固速率范围的主要参数.

关键词: 定向凝固 ; 协同生长 ; 模拟仿真 ; 多孔Al ; 凝固速率

Abstract

The solid/gas eutectic unidirectional solidification process is a new kind of technology fabricating the lotus-type porous structure. Besides having the properties of large specific surface area, excellent sound absorption, penetrating performance etc. of traditional porous materials, the particularity of the lotus-type porous structure makes it has extraordinary mechanical and thermal properties. There is a great potential for lotus-type porous aluminum in the field of lightweight engineering and heat dissipation of chip owing to its low density, outstanding corrosion resistance and high thermal conductivity. However, the fabrication of lotus-type porous aluminum has always been more difficulty than other metals. Until porous aluminum with excellent pore structure was fabricated under very small solidification rates (0.008~0.015 mm/s) and high superheat degrees of melt (240~340 K), the proper processing parameters were recognized to be essential for the coupled growth of solid/gas phases for Al-H2 system, especially the solidification rate. In this work, a three-dimensional time-dependent model describing the evolution of single pore was established based on the theoretical analysis of mass transfer, gas bubble nucleation, pore growth, interruption and detachment. The morphology of single pore under different solidification rates were simulated during the Gasar process by using the finite difference method for Al-H2 system. The research reveals: coupled growth of solid/gas phases can be achieved under the solidification rates from 0.15 mm/s to 0.005 mm/s. The average pore diameter which ranges from 100 mm to 1100 mm increases with decreasing the solidification rate. The pore length also increases while the pore aspect ratio is nearly a constant about 40 with decreasing the solidification rate. The simulated average pore diameter is in good agreement with the experimental values when solidification rate equals to 0.015 mm/s, and then being slightly lower than the experimental ones with decreasing the solidification rate. The diffusion of hydrogen into the melt during the solidification process is regarded as the main reason of the discrepancy between simulated and experimental average pore diameters. The maximum value of the simulated solidification rates for coupled growth of solid/gas phases in Al-H2 system first increases from less than 0.01 mm/s to 0.15 mm/s and then being a constant, while the minimum ones increase from about 0.0001 mm/s to 0.01 mm/s with improving the overheat degree of melt and hydrogen partial pressure. By comparing the relative parameters of Al-H2 and Cu-H2 systems, the solubility of hydrogen is regarded to be the main parameter which determines solidification rates of coupled growth of solid/gas phases for Al-H2 system.

Keywords: unidirectional solidification ; coupled growth ; modeling and simulation ; porous aluminum ; solidification rate

0

PDF (2900KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

杨倩倩, 刘源, 李言祥. 定向凝固藕状多孔Al生长过程的模拟仿真[J]. , 2014, 50(11): 1403-1412 https://doi.org/10.11900/0412.1961.2014.00300

YANG Qianqian, LIU Yuan, LI Yanxiang. MODELING AND SIMULATION OF STRUCTURAL FORMATION OF POROUS ALUMINUM IN GASAR SOLIDIFICATION[J]. 金属学报, 2014, 50(11): 1403-1412 https://doi.org/10.11900/0412.1961.2014.00300

多孔金属材料具有许多不同于其致密结构的优异的物理和力学性能, 如轻质、大比表面积、良好的吸声隔音以及渗透性能等, 正成为一种重要的功能材料. 作为一种新型多孔材料, 采用定向凝固方法制备的藕状多孔金属(又名Gasar), 其圆柱状气孔沿凝固方向规则排布, 且可以实现对孔径和孔隙率的有效控制[1-5]. 藕状多孔金属的特殊结构使其在兼具传统多孔金属特征性能的同时, 还具有特殊的力学[6,7]和热学性能, 尤其是作为一种特殊的微通道结构表现出优异的散热性能, 在高热流密度散热领域展现出重要的应用前景[8-10].

在诸多金属当中, Al因其密度低、耐腐蚀性好和导热性能高等特点, 其藕状多孔结构有望在轻质工程和芯片散热领域得到应用[11]. 然而, 与制备工艺成熟的藕状多孔Cu[12-14], Mg[15-18]等相比, 藕状多孔Al的制备一直比较困难. 张华伟等[19]最先采用与藕状多孔Cu和Mg相同的模铸法Gasar制备工艺来制备藕状多孔Al, 得到孔隙率不仅很低(小于5%), 而且气孔为近球形的多孔Al, 并不是所期望的藕状多孔结构, 这是因为H2在Al中的溶解度太小(0.1 MPa下, H2在Al中溶解度比在Mg和Cu中溶解度低约1~2个数量级), 导致气泡生长过程中没有足够H2扩散进入气泡而发生生长中断. Komissarchuk等[20]也采用同样的普通模铸法, 但在凝固开始时通过减小H2压力来降低气泡生长的阻力, 从而减小气泡生长所需H2扩散供应量, 希望能得到长条形气孔, 结果得到的是孔径在0.2~0.3 mm, 孔长在1.5~2.5 mm的近椭球形气孔(长径比7~8). 虽然孔长有所增加, 但也还不是理想的藕状多孔结构. Ide等[21]采用连铸工艺, 在比普通模铸法低1~2个数量级的极低速(凝固速率v=0.008~0.015 mm/s)和高过热度(240~340 K)的条件下成功制备出了气孔长径比大于10, 孔隙率高达40%的藕状规则多孔Al. 以上实验结果表明, 藕状多孔结构的形成受凝固速率、气体压力、熔体过热度等的共同影响. 只有当试样的凝固速率与其他工艺参数相匹配时, 才能够得到孔分布均匀、孔径一致且具有较大长径比(不小于15)的理想藕状多孔结构, 即实现固/气两相的协同生长. 否则就会导致气孔的迅速中断, 甚至得到近球形气孔. 由于影响因素较多, 且相互耦合, 仅通过实验研究无法获得各种工艺参数以及不同体系的物性参数差异对制备工艺的影响. 与之相比, 模拟仿真能够综合考虑各物性参数和工艺参数的影响, 是一种研究非稳态凝固过程的有效手段.

本工作通过对Gasar凝固过程中的传质、气泡形核、气孔生长、中断及脱离等的理论分析, 建立了一个描述单气孔演变过程的非稳态三维模型. 以Al-H2系为研究对象, 采用有限差分的方法模拟其定向凝固过程中的气孔形貌变化, 得到形成藕状多孔Al的凝固速率范围, 并与实验结果进行了比较. 通过模拟仿真系统研究了熔体过热度和H2分压对固/气两相协同生长凝固速率的影响, 并得出了不同熔体过热度和H2分压条件下藕状多孔Al固/气两相协同生长的凝固速率范围. 通过对Al-H2系和Cu-H2系相关参数的比较分析, 认为H2在金属熔体中的溶解度是决定Al-H2系固/气两相协同生长凝固速率范围的主要参数.

1 仿真理论模型

在铸锭定向凝固过程中, 凝固速率取决于温度场, 进而影响界面前沿熔体中的溶质分布和扩散. 为便于研究凝固速率对气孔结构的影响, 在模型中忽略温度场的计算, 将凝固速率作为初始条件直接输入, 同时对模型做出如下简化: (1) 忽略熔体的充型过程以及熔体在凝固过程中的流动; (2) 假设凝固界面上方熔体各高度处的温度梯度GTL相等, 为熔体过热度DT与熔体高度h的比值; (3) 已凝固固相中的温度梯度GTS为定值, 在模拟过程中作为初始条件直接给定.

图1   理想藕状多孔结构横截面示意图

Fig.1   Schematic of ideal lotus-type porous structure (rG—pore radius at the solid/liquid interface, rS—equivalent radius of circle which is approximately equal to half of the inter-pore spacing)

1.1 传质过程

由于圆柱状气孔的存在, Gasar凝固过程中, 凝固界面前沿H2的浓度在水平方向分布不均匀并产生传质. Gasar的理想横截面结构可以抽象为类似于固/固共晶的正六边形排布, 如图1所示. 为了计算的方便, 模型中将该六边形等效为面积相等的圆形区域[22]. 考虑到H2在固相中的扩散系数比在熔体中小约3个数量级, 对孔结构的影响基本可以忽略[23], 因此, 本工作计算中只考虑熔体中的传质过程, 同时忽略熔体中的对流. 如图2所示, 在凝固界面上建立柱状坐标系, H2在凝固界面前沿熔体中的浓度CL (r, z, t)满足如下微分方程:

图2   熔体中的传质边界条件和坐标系

Fig.2   Solute diffusion boundary conditions and the corresponding coordinate selection (zl/g—gas/liquid interface, z0—height of remaining melt)

CLt=DL2CLr2+1rCLr+2CLz2

式中, DL为H2在熔体中的扩散系数, r和z分别为径向和轴向坐标, t为凝固时间. 熔体中的初始H2浓度CL(r, z, 0)可由Sieverts定律[24]得到:

CLr,z,0=C0=ξ(Tp)pH2

式中, pH2为熔体上方气氛中的H2分压; Tp为熔体初始温度; x(Tp)=exp(-A/Tp+B)为系数项, A和B为对应体系的溶解度常数. 边界条件如图2所示.

(1) 边界1 (自由液面):

CLr,z0,t=C0=ξTppH2

(2) 边界2 (侧向边界和中心处):

CL0,z,tr=CLrS,z,tr=0

(3) 边界3 (气/液边界):

CLr,zl/g,t=ξTmpbpb=pat+ρLgh+2σlg/RG

式中, rS为图1所示圆的等效半径; zl/g为气/液界面; Tm为金属熔点; pb和pat分别为界面前沿气泡内的气体压力和凝固时熔体上方的气室压力( pat=pH2+pAr, pAr为熔体上方的Ar分压); ρL, g和h分别为金属熔体的密度、重力加速度以及气泡上方的熔体高度; σlg和RG分别为金属液/气界面张力和界面处球缺的曲率半径.

(4) 边界4 (固/液边界):

CLr,0,tz=-vDL1-k0CLr,0,tr[rG,rS]

式中, v和k0分别为试样的凝固速率以及H2在金属固/液两相中的分配系数.

1.2 气泡生长模型

随着凝固界面的推进, 新凝固固相析出的H2在界面前沿的熔体中富集并部分扩散进入气孔内, 使得气孔与固相一起生长. 气孔在生长过程中会形成凸入熔体中的气泡. 本模型中, 假设初始状态下气泡顶部的凸起为半球形, 之后气泡形状随着H2扩散量的变化作出相应调整. 图3为一个时间步长Dt前后气孔的生长示意图. 气孔顶部球缺的曲率半径RG1, 固/液界面处的气孔半径rG1以及界面角θ1经过Dt时间的生长后分别变为RG2, rG2和θ2. 这3个变量之间始终满足如下关系式:

图3   Gasar凝固中Dt时间前后气孔生长示意图

Fig.3   Schematic describing the gas pore growth in a short time Dt for Gasar process (v—solidificationrate; Dt—time increment; rG1, rG2—pore radius at the solid/liquid interface before and after Dt time; RG1, RG2—pore radius of curvature before and after Dt time; θ1, θ2—interfacial angle before and after Dt time; DV—increased volume of pore in Dt time; θ′—interfacial angle between the old hemispherical before Dt time and the new solidification height after Dt time)

sinθ=rG/RG

当熔体中的传质以及温度、浓度分布已知后, 可以通过以下方法计算Dt时间后气孔的结构参数:

(1) Dt时间内通过气/液界面进入气泡内的H2量DN与Dt时间前后气孔内的H2总量之间满足:

ΔN=ΔNr+ΔNz=-ΔtDLzl/gCLrdSr+zl/gCLzdSzNbt+Δt=Nbt+ΔN

式中, NbtNbt+Δt分别为Dt时间前后气孔内的H2总量, DNr和DNz分别为Dt时间内在气/液界面处通过径向和轴向扩散进入气泡内的H2量, Sr和Sz分别表示H2在气/液界面处沿径向和轴向扩散的面积.

(2) Dt时间后凝固界面处气孔的半径为:

rG2=RG12-RG1cosθ1-vΔt2RG1cosθ1-vΔt<RG10RG1cosθ1-vΔtRG1

(3) Dt时间内由于H2往气泡内的扩散所增加的气孔体积DV可以表示为:

ΔV=V2-VV2=13πRG232+3cosθ2-cos3θ2V=13πRG132+3cosθ-cos3θsinθ=rG2/RG1

如图3所示, V2为Dt时间后气孔顶部新球缺的体积, V′为Dt时间前气孔顶部的球缺在Dt时间后界面前沿剩余部分的体积, θ′为Dt时间前气孔顶部的球缺在Dt时间后的剩余部分与新凝固界面间的界面角.

(4) 根据理想气体定律, DV可由下式求得:

ΔV=Vbt+Δt-Vbt=Nbt+ΔtRgTbt+Δtpbt+Δt-Vbt

式中, VbtVbt+Δt分别为Dt时间前后气孔的总体积, Tbt+Δtpbt+Δt分别为Dt时间后气孔内的平均温度和气体压力, Rg为理想气体常数.

在已知t时刻气孔各参数的情况下, 联立式(7)~(11)可以计算得到Dt时间后气孔的各参数值.

1.3 气泡形核、气孔中断与脱离

当熔体中以均质形核的方式出现球形胚胎气泡时, 体系Gibbs自由能的变化由体积自由能的降低( -4πRG3ΔGv/3)和界面自由能的升高( 4πRG2σlg) 2部分组成[25]. 在经典形核理论中, DGv为液、固相单位体积的Gibbs自由能, 而对于Gasar凝固中的气泡形核, 张华伟等[26]认为其为气泡反抗气泡外阻力、排开单位体积熔体所作的体积功, 并给出气泡在熔体中的临界形核半径Rc表达式为:

Rc=2σlg/3pat+ρLgh

由上式求得本工作对应工艺参数条件(pat=0.5 MPa, h=0.15 m)下多孔Al的临界气泡形核半径为1.2 μm.

根据理想气体定律 V=nRgT/pH2, 由于固相温度梯度的存在, 气孔内气体的平均温度T随着凝固界面的升高而逐渐降低, 由于气孔内的H2压力 pH2几乎保持不变, 势必会导致新增体积V随凝固高度的增加而减小. 当扩散进入气孔的H2量n所引起的体积增加V不足以维持气孔的正常生长时, 气泡半径减小, 界面角q会相应增大, 当q180°时, 气孔发生中断. 与之相反, 若气孔内H2的扩散量远大于维持气孔正常生长所需的量时, 气泡半径增大, 界面角相应减小, 界面前沿的气泡大于半球形, 当q0°时, 界面处的气泡发生脱离. 本计算中, 设q≥160°时气孔发生中断; 与此相应, 设q≤20°时,界面处的近球形气泡与界面间形成缩颈并脱离.

一旦气泡脱离凝固界面, 就会以一定的速率上浮逸出. 一般认为半径1 mm以下球形气泡的上浮速率可由Stokes公式计算得到. 然而, 气泡自由表面存在的表面张力梯度会使气泡受到向上的表面张力, 从而产生Marangoni运动. 陈自发[27]认为熔体中气泡的实际上浮速率vf应由Stokes运动和Marangoni运动共同决定, 并对气泡的上浮速率进行修正, 给出如下计算公式:

图4   浓度场网格划分示意图

Fig.4   Schematic of mesh devision of concentration field

vf=2gρL9ηRG2+GTLdσlg/dT2ηRG

式中, h为熔体黏度, dslg/dT为气/液界面张力对温度的导数. 当vf大于v时, 脱离后的气泡不能被凝固界面重新捕获, 最终上浮逸出.

2 仿真结果与分析

在上述理论模型的基础之上编写计算机程序, 对不同凝固速率下藕状多孔Al单气孔的生长过程进行模拟, 并与实验结果进行比较. 该过程中浓度场的计算采用有限差分法. 如图4所示, 将柱状坐标系原点固定在凝固界面中心处, 实时划分凝固界面附件的熔体网格, 把熔体划分为Nr×Nz个单元(Nr和Nz分别为径向和横向的单元数量), (i, k)代表分别沿r和z方向第i和k (0≤i≤Nr-1, 0≤k≤Nz-1)个单元, 单元尺寸为p×[(i+1)2-i2]×Dr2×Dz (Dr和Dz分别为径向和轴向的单元尺寸). k=0时单元位于凝固界面上, i=0时圆柱形单元体的中心位于熔体轴对称中心线上, 单元属性由单元体中心的位置决定. 对于i≠0的圆环柱单元体(i, k), 其属性则取决于该单元内[(i+0.5)×Dr, (k+0.5)×Dz]的坐标位置与气泡球心位置间的距离与RG的相对大小. 如图5所示, 若该距离大于RG (虚线所示为气泡边界), 则为熔体单元. 反之, 则认为该单元体为气相单元(蓝色区域所示), 气体压力为pb, 气泡内等效H2浓度Cb为气/液边界处的H2浓度CL(r, zl/g, t).

图5   浓度场气相单元确定示意图

Fig.5   Schematic of determination of gas phase in concentration field(a) cross-section(b) longitudinal-section

分别计算Dt时间后熔体单元的H2浓度和该时间段内通过气/液界面由熔体单元扩散进入气泡内H2的物质的量, 进而对整个浓度场进行求解, 最终可得到藕状多孔Al的气孔结构.

2.1 凝固速率对气孔结构的影响

Liu等[28]和李再久等[29]的研究结果认为, 藕状多孔金属稳态生长时的凝固速率v和气孔间距L之间满足类似于Jackson-Hunt关系的关系式: vL2=F (F是由合金系、熔体保温温度和气体压力共同决定的常数). 根据Ide等[21]采用普通连铸法对Al-H2系的实验研究, v为0.015 mm/s时, 藕状多孔Al的气孔平均孔径2rG和孔隙率分别约为700 mm和10%, 与之对应的气孔间距约为2.2 mm. 以此为基准, 本工作采用与上述实验相同的工艺参数( pH2= pAr=0.25 MPa, DT=290 K, GTS=10 K/mm), 并分别模拟v为0.5, 0.15, 0.05, 0.015, 0.01, 0.005和0.0015 mm/s时, Al-H2系单气孔的纵截面生长形貌, 其孔间距(L≈2rS)分别取为0.3, 0.7, 1.2, 2.2, 2.7, 3.8和7.0 mm.

图6为不同凝固速率下模拟出的Al-H2系纵截面生长形貌图. 可以看出, v为0.5 mm/s时, 由于H2扩散量不足, 气泡形核后孔径迅速减小并发生中断, 只得到近球形气泡. 当v降低至0.15 mm/s, 圆柱状气孔开始形成, 但孔径很小(不足200 mm); 随着v不断降低至0.005 mm/s, 气孔孔径逐渐增大, 固/气两相能够保持协同生长. 然而, 当v继续降低到0.0015 mm/s时, 由于孔径增长过快, 界面前沿的气泡达到脱离条件而上浮逸出. 因此, 在v为0.15~0.005 mm/s范围内, Al-H2系固/气两相能够维持协同生长. 在此范围内, 在凝固初始阶段, 单位时间内H2的扩散量大于气孔生长所需的量, 气孔孔径随凝固高度的增加迅速增大至最大值, 之后由于固相温度梯度的存在缓慢降低, 直至最后发生中断.

图6   不同凝固速率下气孔生长纵截面形貌的仿真图

Fig.6   Cross-sections parallel to the solidification direction of porous Al under solidification rate v= 0.5 mm/s (a), 0.15 mm/s (b), 0.05 mm/s (c), 0.015 mm/s (d), 0.01 mm/s (e), 0.005 mm/s (f) and 0.0015 mm/s (g) (pH2=pAr=0.25 MPa, superheat degree of melt DT=290 K, temperature gradient of solid phase GTS=10 K/mm, H—the solidification height, pH2—hydrogen pressure, pAr—argon pressure )

图7给出了不同凝固速率下气孔平均孔径的模拟值与实验值[21]的比较. 可以看出, 气孔的平均孔径模拟值分布在100~1100 mm之间, 且随v的降低逐渐增加. 当v=0.015 mm/s时, 模拟值与实验值吻合较好, 之后随v的降低, 模拟孔径略低于实验值. 图8给出了不同凝固速率下气孔长度和长径比的模拟结果. 可以看出, 气孔长度随凝固速率的升高而显著降低, 而长径比则基本保持在40左右.

图7   不同凝固速率下气孔平均孔径的模拟值与实验值[21]的比较

Fig.7   Comparisons of the simulated and experimental[21] average pore diameter results

2.2 凝固速率对浓度场分布的影响

图9显示了凝固速率为0.15, 0.05和0.015 mm/s时, 气孔稳定生长时界面前沿熔体中的浓度场分布(对应凝固高度为该速率下气孔总长度一半左右). 可以看出, 沿水平方向(r轴方向), 凝固界面处熔体中的H2浓度逐渐增大, 随着v的降低, 熔体中的H2扩散更加充分, 径向浓度梯度逐渐减小. 沿高度方向(z轴方向), v较大时(v=0.15 mm/s), 气泡上方熔体中的H2浓度随高度增加先升高后降低至初始浓度值. 在此情况下, 对应的气泡内的H2由来自于界面前沿的熔体向固相转变时排出的H2以及气泡上方一定范围内熔体中向气泡内扩散的H2组成. 随着v的降低, 气泡上方熔体中的H2浓度随高度的增加逐渐增大, 且v越小, H2浓度随高度增加越显著. 这是因为随着v的降低, 气泡周围影响其生长的熔体区域相应增大, 相当于更大范围熔体中的H2扩散进入气孔内, 产生单一的随高度增加而逐渐增大的浓度变化, 与之对应的气泡内的H2主要来自于气泡上方更大范围熔体中向气泡内扩散的H2.

图8   不同凝固速率下的气孔长度和长径比的模拟结果

Fig.8   Simulated results of pore length and pore aspect ratio

可见, 随着v的降低, 气孔内的H2由主要依靠界面前沿的熔体向固相转变时排出的H2提供转变为主要依靠气泡上方更大范围熔体中的H2向气泡内的扩散供应. 而在采用普通模铸法制备多孔Al时, 由于v太大(约为0.5~1.5 mm/s), 且H2在Al熔体中本身的溶解度较低, 气泡内的H2主要依靠界面前沿熔体向固相转变时排出的H2来提供, 其扩散量不足以维持气孔的正常生长. 而通过降低v至0.15~0.005 mm/s, 气泡周围一定范围熔体中的H2向气泡内的扩散量能够满足多孔Al中气孔正常生长的需要, 因此固/气两相能够维持协同生长.

图9   不同凝固速率下气孔生长界面前沿的浓度场分布模拟结果

Fig.9   Simulated hydrogen concentration distributions in the melt ahead of the front of solidification under v=0.15 mm/s (a), 0.05 mm/s (b) and 0.015 mm/s (c) (pH2=pAr=0.25 MPa, DT=290 K, GTS=10 K/mm)

从图9中还可以看出, v越低, 气泡上方熔体中的浓度梯度越高(熔体中的H2浓度均低于熔体上表面的H2浓度), 熔体中向气泡提供溶质的影响区域越大. 然而, 实际的实验过程中, 随着凝固界面的不断升高, 气泡上方剩余熔体的高度h不断变小, 熔体上表面的H2由于向气泡内的扩散会不断减少. 根据气液两相的浓度平衡原理, 气室内的H2会源源不断地通过界面扩散到熔体中, 以保持熔体上表面始终为C0, 如图10所示. 最终, 导致扩散进入气孔内的H2量增加, 孔径增大. 而模拟过程中, 浓度场计算高度始终为定值, 与实际实验制备过程的不断变化的浓度场扩散高度存在差异. 因此, 认为实际凝固过程中熔体上方的H2向熔体内的不断扩散是造成图7中v<0.015 mm/s时气孔的模拟平均孔径低于实验值的主要原因. 与此同时, 文献[21]中并未给出界面前沿熔体中的温度梯度GTL, 模拟过程中则始终假定剩余熔体中各高度处的GTL相同, 并取值为熔体过热度DT与剩余熔体高度h的比值. 因此, 真实GTL和模拟值之间的差异也是导致气孔的模拟平均孔径低于实验值的重要原因.

图10   H2往熔体内扩散的示意图

Fig.10   Schematic of diffusion of hydrogen into the melt

2.3 工艺参数对协同生长凝固速率范围的影响

图11给出了不同DT和 pH2条件下Al-H2系固/气两相协同生长的凝固速率范围. 可以看出, 当DT约在100 K和 pH2约在0.1 MPa以下时, 对应凝固速率范围的最大值随DT和 pH2的升高由约0.01 mm/s逐渐增大至约0.15 mm/s. 然而, 一旦超过上述临界值, 最大凝固速率随DT和 pH2的升高保持不变. 与之不同的是, 该范围的最小凝固速率在本工作所示参数范围内随着DT和 pH2的升高由约0.0001 mm/s逐渐增加至约0.01 mm/s.

根据式(13), vf由Stokes运动速度和Marangoni迁移速度共同决定. 与此同时, 张华伟[30]给出饱和溶解有H2的熔体在凝固开始前发生H2逸出的工艺判据为:

图11   不同熔体过热度DT和H2分压pH2条件下Al-H2系固/气两相协同生长的凝固速率范围

Fig.11   Simulated range of solidification rate for coupled growth of solid/gas phases under different superheat degrees of melt (pH2=pAr=0.25 MPa, GTS=10 K/mm) (a) and hydrogen partial pressures (pH2+pAr=0.5 MPa, GTS=10 K/mm, DT=290 K) (b) for Al-H2 system

C0Cbvfν

式中, Cb为气泡的生长浓度( Cb=ξTmpb). 当vf=v时, 将式(13)代入上式, 求得H2逸出的临界半径 RG表达式为[26]:

RG=-9GTLdσlg/dT+81GTLdσlg/dT2+288ηgρLv8gρL

根据式(15)可求出不同凝固速率对应的 RG和毛细压力(pc=2slg/ RG), 结合式(14)进而得到凝固开始前H2逸出的临界熔体过热度DT′和H2分压 pH2分别为[30]:

ΔT=11Tm-12Alnpat+ρLgh+pc/pH2-TmpH2=ξ2Tmpat+ρLgh+pcξ2Tm+ΔT

图12给出了不同凝固速率下, 多孔Al凝固前熔体中H2逸出的DT′和 pH2. 不同凝固速率下 ΔTpH2均随着凝固速率的增大而逐渐降低(即v越大, H2逸出越严重). v=0.15 mm/s时, DT′和 pH2分别为108 K和0.07 MPa. 即DT和 pH2在上述临界值以下且v≤0.15 mm/s时, 提高熔体过热度和H2分压均能有效增加熔体中的H2浓度进而提高固/气两相协同生长的最大凝固速率; 一旦超过上述临界值且v≥0.15 mm/s时, 由于更高过热和H2分压, 使多溶入熔体中的H2在降温到凝固温度之前就逸出了, 熔体开始凝固时的H2浓度分别为相应凝固速率下临界熔体过热度和H2分压的浓度, 固/气两相协同生长的最大凝固速率也不会继续升高.

在v分别低至0.01和0.013 mm/s时, 对应的DT′和 pH2分别约为500 K和0.5 MPa. 图11中, 不同熔体过热度和H2分压下固/气两相协同生长的最小凝固速率基本位于上述临界值以下, 即在该速率范围内, 熔体中的H2不会逸出, H2浓度随DT和 pH2的升高而逐渐增大, 相同凝固速率下形核后的气泡更容易达到脱离条件而上浮, 以至于不能实现与固相的协同生长, 因此, 固/气两相协同生长的最小凝固速率随DT和 pH2的升高而逐渐增大.

图12   不同凝固速率下Al-H2系中H2逸出的临界熔体过热度ΔT′和H2分压pH2

Fig.12   Critical superheat degree of melt ΔT(pH2=pAr=0.25 MPa, GTS=10 K/mm) and hydrogen partial pressure pH2(pH2+pAr=0.5 MPa, GTS=10 K/mm, DT=290 K) under different solidification rates for Al-H2 system

根据上述分析, Gasar凝固中固/气两相的协同生长主要取决于单位时间内H2的扩散量Q1与气孔生长所需的H2量Q2是否匹配. 在总压一定的情况下, 前者主要取决于熔体中H2的扩散系数DL以及熔体中的H2浓度CL, 后者则主要由v决定. 即H2的扩散系数和H2浓度决定了该体系Gasar凝固中固/气两相协同生长的凝固速率. 其中, H2在熔体中的扩散系数由体系本身决定, 熔体中H2浓度的提升虽然可以通过在一定范围内增大熔体过热度DT和H2分压 pH2来实现, 但由于H2逸出的影响, 调整上述工艺参数对提升固/气两相协同生长凝固速率的影响有限.

因此, 文献[19]采用普通模铸法制备多孔Al, 其对应的凝固速率过大, 导致Q1始终小于Q2, 是造成气孔迅速中断的主要原因. 而Komissarchuk等[20]采用普通模铸法在凝固开始时采用降低H2压力的方法, 虽然一定程度上减小了气孔生长所需的H2量, 但因其孔径较小且对应凝固速率下的固相温度梯度较大, 气孔内气体发生较为严重的冷却收缩, 导致其藕状多孔结构不理想.

另外, 通过对Al-H2和Cu-H2的比较发现, H2在金属熔点温度附近熔体中的扩散系数相差并不大(分别为0.070[31]和0.058[32] mm2/s), 而0.1 MPa下H2在熔体中的溶解度则相差数十倍(分别为0.7[24]和18.8[24] mol/m3). 因此, 分析认为H2在熔体中的溶解度是决定Al-H2系固/气两相协同生长凝固速率范围的主要参数.

3 结论

(1) 通过建立描述Gasar凝固中单气孔演变过程的非稳态三维模型, 模拟了不同凝固速率下Al-H2单气孔的生长形貌, 并得出了藕状多孔Al固/气两相协同生长的凝固速率范围为0.15~0.005 mm/s. 该范围内气孔的平均孔径分布在100~1100 mm之间, 且随凝固速率的降低逐渐增加, 气孔长度亦随凝固速率的降低逐渐增加, 长径比则基本保持在40左右; 当凝固速率为0.015 mm/s时, 气孔孔径的模拟值与实验值吻合较好, 之后随凝固速率的降低, 模拟孔径略低于实验值, 分析认为实际凝固过程中熔体上方的H2向熔体内的不断扩散是导致该差异的主要原因.

(2) 得出了不同熔体过热度和H2分压条件下Al-H2系固/气两相协同生长的凝固速率范围. 随着熔体过热度和H2分压的增加, 该范围对应凝固速率的最大值由不足0.01 mm/s先逐渐增加之后稳定在0.15 mm/s, 最小值则由0.0001 mm/s左右逐渐增加至约0.01 mm/s. 通过对Al-H2系和Cu-H2系相关参数的比较分析表明, H2在金属熔体中的溶解度是决定Al-H2系固/气两相协同生长凝固速率范围的主要参数.

The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。


参考文献

[1] Banhart J. Prog Mater Sci, 2001; 46: 559

[本文引用: 1]     

[2] Nakajima H. Prog Mater Sci, 2007; 52: 1091

[3] Jiang G R, Li Y X, Liu Y. Metall Mater Trans, 2010; 41A: 3405

[4] Li Z J, Jin Q L, Yang T W, Jiang Y H, Zhou R. Acta Metall Sin, 2013; 49: 757

(李再久, 金青林, 杨天武, 蒋业华, 周 荣. 金属学报, 2013; 49: 757)

[5] Shapovalov V I. US Pat, 5181549, 1993

[本文引用: 1]     

[6] Simone A E, Gibson L J. Acta Mater, 1996; 44: 1437

[本文引用: 1]     

[7] Hyun S K, Nakajima H. Mater Sci Eng, 2003; A340: 258

[本文引用: 1]     

[8] Muramatsu K, Ide T, Nakajima H, Eaton J K. J Heat Transfer, 2013; 135: 072601

[本文引用: 1]     

[9] Chen L T, Zhang H W, Liu Y, Li Y X. Acta Metall Sin, 2012; 48: 329

(陈刘涛, 张华伟, 刘 源, 李言祥. 金属学报, 2012; 48: 329)

[10] Zhang H W, Chen L T, Liu Y, Li Y X. Int J Heat Mass Transfer, 2013; 56: 172

[本文引用: 1]     

[11] Nakajima H. Proc Jpn Acad Ser B, 2010; 86: 884

[本文引用: 1]     

[12] Zhang H W, Li Y X, Liu Y. Acta Metall Sin, 2006; 42: 1165

[本文引用: 1]     

(张华伟, 李言祥, 刘 源. 金属学报, 2006; 42: 1165)

[本文引用: 1]     

[13] Park J S, Hyun S K, Suzuki S, Nakajima H. Acta Mater, 2007; 55: 5646

[14] Zhou R, Li Z H, Jiang Y H, Zhou R F, Yang T W, Jin Q L. Chin Pat, 200910094262.8, 2009

[本文引用: 1]     

(周 荣, 黎振华, 蒋业华, 周荣锋, 杨天武, 金青林. 中国专利, 200910094262.8, 2009)

[本文引用: 1]     

[15] Zahrani M M, Meratian M, Kabiri Y. Mater Lett, 2012; 85: 14

[本文引用: 1]     

[16] Zhou R, Jiang Y H, Li Z H, Jin Q L, Yang T W, Zhou R F. Chin Pat, 200910094261.3, 2009

(周 荣, 蒋业华, 黎振华, 金青林, 杨天武, 周荣锋. 中国专利, 200910094261.3, 2009)

[17] Tane M, Nakajima H. J Mater Res, 2008; 23: 849

[18] Liu Y, Li Y X, Zhang H W. Acta Metall Sin, 2004; 40: 1121

[本文引用: 1]     

(刘 源, 李言祥, 张华伟. 金属学报, 2004; 40: 1121)

[本文引用: 1]     

[19] Zhang H W, Li Y X, Liu Y. Acta Metall Sin, 2007; 43: 11

[本文引用: 1]     

(张华伟, 李言祥, 刘 源. 金属学报, 2007; 43: 11)

[本文引用: 1]     

[20] Komissarchuk O, Xu Z B, Hao H, Zhang X L, Vladimir K. China Foundry, 2014; 11: 1

[本文引用: 2]     

[21] Ide T, Iio Y, Nakajima H. Metall Mater Trans, 2012; 43A: 5140

[本文引用: 5]     

[22] Drenchev L, Sobczak J, Malinov S, Sha W. Modell Simul Mater Sci Eng, 2006; 14: 663

[本文引用: 1]     

[23] Li Y X,Wu A P. Principle of Materials Processing. Beijing: Tsinghua University Press, 2005: 83

[本文引用: 1]     

(李言祥,吴爱萍. 材料加工原理. 北京: 清华大学出版社, 2005: 83)

[本文引用: 1]     

[24] Zhang H W, Li Y X, Liu Y. Acta Metall Sin, 2007; 43: 113

[本文引用: 3]     

(张华伟, 李言祥, 刘 源. 金属学报, 2007; 43: 113)

[本文引用: 3]     

[25] Kurz W, Fisher D J. Fundamental of Solidification. Switzerland: Trans Tech Publications, 1998: 22

[本文引用: 1]     

[26] Zhang H W, Li Y X. Acta Phys Sin, 2007; 56: 4864

[本文引用: 2]     

(张华伟, 李言祥. 物理学报, 2007; 56: 4864)

[本文引用: 2]     

[27] Chen Z F. Bachelor Thesis, Tsinghua University, Beijing, 2011

[本文引用: 1]     

(陈自发. 清华大学学士学位论文, 北京, 2011)

[本文引用: 1]     

[28] Liu Y, Li Y X. Scr Mater, 2003; 49: 379

[本文引用: 1]     

[29] Li Z J, Jin Q L, Yang T W, Zhou R, Jiang Y H. Acta Metall Sin, 2014; 50: 507

[本文引用: 1]     

(李再久, 金青林, 杨天武, 周 荣, 蒋业华. 金属学报, 2014; 50: 507)

[本文引用: 1]     

[30] Zhang H W. PhD Dissertation, Tsinghua University, Beijing, 2006

[本文引用: 2]     

(张华伟. 清华大学博士学位论文, 北京, 2006)

[本文引用: 2]     

[31] Fisher D J. Hydrogen Diffusion in Metals—A 30 Year Retrospective. Switzerland: Scitec Publications Ltd, 1999: 37

[本文引用: 1]     

[32] Sigrist P, Feichtinger H K, Marincek B. Z Phys Chem-Frankfurt, 1977; 107: 211

[本文引用: 1]     

/