相场法模拟中取向与热扩散参数对Ni-Sn合金枝晶形貌的影响.pdf

收稿日期:2 0 2 2 1 2 0 8基金项目:国家自然科学基金资助项目(1 1 6 7 4 2 3 3)。作者简介:高 明(1 9 6 4),男,黑龙江伊春人,沈阳师范大学教授,博士。第4 1卷 第2期2023年 4月沈阳师范大学学报(自然科学版)J o u r n a l o fS h e n y a n gN o r m a lU n i v e r s i t y(N a t u r a l S c i e n c eE d i t i o n)V o l.4 1N o.2A p r.2 0 2 3文章编号:1 6 7 3 5 8 6 2(2 0 2 3)0 2 0 1 2 8 0 5相场法模拟中取向与热扩散参数对N i-S n合金枝晶形貌的影响高 明,焦艺伟(沈阳师范大学 物理科学与技术学院,沈阳 1 1 0 0 3 4)摘 要:利用相场模型与温度场进行耦合的控制方程来模拟N i-S n二元合金凝固过程中的枝晶生长情况。依据热力学定律、自由能泛函与熵泛函原理,建立N i-S n合金相场和温度场耦合的模型。用相场法模拟合金在非均匀温度下枝晶生长过程和枝晶生长形貌的情况。对形成的图像进行对比,研究各向异性强度参数与热扩散参数对枝晶生长形貌和生长速度的影响。模拟研究表明:在非均匀温度场下,当各向异性强度参数增大时,枝晶的生长速度会随之增大,并且二次支臂的生长数量也会增加。此外,当热扩散参数增大时,枝晶的生长速度会随之减小,枝晶主干变细,并且二次支臂的生长数量会有明显减少。关 键 词:相场法;枝晶生长;数值模拟;N i-S n合金中图分类号:T G 1 1 3 文献标志码:Ad o i:1 0.3 9 6 9/j.i s s n.1 6 7 3 5 8 6 2.2 0 2 3.0 2.0 0 6E f f e c t so fo r i e n t a t i o n a n dt h e r m a ld i f f u s i o n p a r a m e t e r so nd e n d r i t em o r p h o l o g yo fN i-S na l l o y i np h a s e f i e l ds i m u l a t i o nG A OM i n g,J I A OY i w e i(C o l l e g eo fP h y s i c a lS c i e n c ea n dT e c h n o l o g y,S h e n y a n gN o r m a lU n i v e r s i t y,S h e n y a n g1 1 0 0 3 4,C h i n a)A b s t r a c t:T h ep h a s ef i e l d m o d e lc o u p l e d w i t ht h et e m p e r a t u r ef i e l di su s e dt os i m u l a t et h ed e n d r i t eg r o w t ho fN i-S nb i n a r ya l l o yd u r i n gs o l i d i f i c a t i o n.B a s e do nt h et h e r m o d y n a m i c l a w,f r e ee n e r g yf u n c t i o n a la n de n t r o p yf u n c t i o n a l,t h ec o u p l i n g m o d e lo f N i-S na l l o y p h a s ef i e l da n dt e m p e r a t u r e f i e l di se s t a b l i s h e d.T h ep h a s ef i e l dm e t h o di su s e dt os i m u l a t et h ed e n d r i t eg r o w t hp r o c e s sa n d m o r p h o l o g yo ft h ea l l o ya tn o n-u n i f o r mt e m p e r a t u r e.T h ei m a g e sa r ec o m p a r e dt os t u d yt h ee f f e c t so fa n i s o t r o p i cs t r e n g t hp a r a m e t e r sa n dt h e r m a ld i f f u s i o np a r a m e t e r so nd e n d r i t eg r o w t hm o r p h o l o g ya n dg r o w t hr a t e.T h e s i m u l a t i o nr e s u l t s s h o wt h a t t h eg r o w t hr a t eo f d e n d r i t e sa n dt h en u m b e ro f s e c o n d a r ya r m sw i l l i n c r e a s ew i t ht h e i n c r e a s eo f a n i s o t r o p i cs t r e n g t hc o e f f i c i e n ti nt h en o n-u n i f o r mt e m p e r a t u r e f i e l d.I na d d i t i o n,w h e nt h e t h e r m a l d i f f u s i o nc o e f f i c i e n t i n c r e a s e s,t h eg r o w t hr a t eo fd e n d r i t e sw i l l d e c r e a s e,t h em a i nd e n d r i t ew i l lb e c o m e t h i n n e r,a n dt h en u m b e ro f s e c o n d a r ya r m sw i l lb es i g n i f i c a n t l yr e d u c e d.K e yw o r d s:p h a s e-f i e l dm e t h o d;d e n d r i t eg r o w t h;n u m e r i c a l s i m u l a t i o n;N i-S na l l o y材料的力学性能与其微观组织密切相关,枝晶是最常见的微观结构组织,所以它的形貌对材料的宏观特性有着至关重要的作用1。金属凝固过程中由于传热和动量传输的复杂性,需要控制大量参数和复杂条件。一般说来,枝晶的尺度大多在微米级,通过实验手段对其凝固过程中的生长进行研究较为困难,而采用相场法能够实现枝晶生长过程微观与宏观的统一23。由于相场法能够引入序参量,这就可以避免追踪固液界面,使模拟效果更接近真实情况4。相场法基于统计物理学,以G i n z b u r g-L a n d a u相变理论为基础,使用微分方程来表示复杂的扩散效应、有序化势和热力学驱动力5。引入相场变量描述固相与液相界面上的状态,并确定相场变量的控制方程与温度场的控制方程,从而更好地反映出金属凝固的真实情况67。本文采用相场与温度场耦合的方法,模拟了N i-S n二元合金枝晶生长过程,研究了各向异性强度参数和热扩散参数对枝晶生长的影响,并且分析了各参数对二元合金枝晶生长形貌及速度的影响,为合理选择控制参数提供了理论支持8。1 N i-S n二元合金相场模型1.1 相场控制方程基于热力学一致性,推导相场控制方程的过程可分为2种,一种为自由能泛函推导,另一种为熵泛函推导。本文以熵泛函理论为基础,结合热力学第二定律,假设系统密度均匀,液相无对流,并且在绝热的前提下推导相场方程91 0,过程如下:t=M22-(1-c)HA(,T)-c HB(,T)(1)=1+c o s(j(-0)(2)其中:为相场参量;M为相场参数;T为无量纲温度;为界面厚度有关的参数;上标A,B分别表示合金中组分N i和S n;c为合金中N i的组分;为各向异性强度参数;j为各向异性强度参数的模数。1.2 温度场控制方程温度场控制方程由热力学第一定律和热传导公式推导1 1,Tt=2T+T,可建立二元合金凝固过程温度场控制方程,如下:Tt=DT2T-(1-c)LA+c LBh()CPt+(LA-LB)h()CPct(3)其中:为无量纲热扩散参数,与潜热成正比;DT为合金的热扩散参数;CP为合金的比热容;LA,LB分别为合金中的组分N i,S n金属的熔化潜热;c为合金中N i的组分。2 数值模拟本文采用有限差分法,把模拟区域分成3 0 03 0 0网格,设初始晶核半径为r,并且使其位于计算区表1 N i-S n合金物性参数T a b l e1 P h y s i c a l p a r a m e t e r so fN i-S na l l o y物性参数N i(A)S n(B)TM/K17 2 6.1 55 0 5.0 4CP/(JKc m-3)4 6 02 2 0L/(Jc m-3)1 77 1 070 8 0DS/(c m2s-1)81 0-681 0-6DL/(c m2s-1)51 0-951 0-9域中心,并设置如下条件1 2:i-Nx22+j-Ny22r2,=1i-Nx22+j-Ny22r2,=0(4)在网格边界区域,选择Z e r o-N e u m a n n边界条件1 31 4。N i-S n合金的物性参数1 5见表1。3 模拟结果与分析3.1 各向异性强度参数对枝晶生长的影响3.1.1 各向异性强度参数对枝晶形貌的影响模拟时,只改变各向异性强度参数,即=0.0 0 5,=0.0 0 7和=0.0 0 9,并且选择了3个不同的时间步长进行计算,计算结果绘制为3组图片(图1)。其中,a,b,c对应着时间步长n=1 0 0 0,n=2 0 0 0,n=30 0 0;1,2,3则分别对应每个时间步长下,各向异性强度参数分别为0.0 0 5,0.0 0 7和0.0 0 9时的枝晶形貌。921 第2期 高 明,等:相场法模拟中取向与热扩散参数对N i-S n合金枝晶形貌的影响n=10 0 0n=20 0 0n=30 0 0(a 1)n=10 0 0,=0.0 0 5(a 2)n=20 0 0,=0.0 0 5(a 3)n=30 0 0,=0.0 0 5(b 1)n=10 0 0,=0.0 0 7(b 2)n=20 0 0,=0.0 0 7(b 3)n=30 0 0,=0.0 0 7(c 1)n=10 0 0,=0.0 0 9(c 2)n=20 0 0,=0.0 0 9(c 3)n=30 0 0,=0.0 0 9图1 不同各向异性强度参数下的枝晶形貌F i g.1 D e n d r i t em o r p h o l o g yu n d e rd i f f e r e n t a n i s o t r o p i cs t r e n g t hp a r a m e t e r s图2 各向异性强度参数对枝晶生长速度的影响F i g.2 T h e i n f l u e n c eo f a n i s o t r o p yc o e f f i c i e n to nd e n d r i t eg r o w t hr a t e从图1可以看出,当各向异性强度参数值为0.0 0 5时,晶粒主要沿着主轴方向生长,没有体现某一方向上的生长优势,二次支臂生长也不明显;当各向异性强度参数值为0.0 0 7时,随着各向异性强度增强,晶粒在主轴方向的支臂变得发达且细长,同时也能发现有二次支臂的生长;当各向异性强度参数值为0.0 0 9时,可以看出较为明显的二次支臂生长,其数量也有增加,这说明各向异性强度参数的增大会导致固相与液相接触面的不稳定扰动出现。图1从纵向对比,可以发现时间步长对枝晶生长的影响。当各向异性强度参数相同时,时间步长的增加会使主干支臂长度和二次支臂的数量都有所增加。3.1.2 各向异性强度参数对枝晶尖端生长速度的影响由图1可知,相同时间步长下枝晶尖端的生长速度也并不相同。为了探究枝晶尖端生长速度与各向异性强度参数之间的关系,绘制了不同各向异性强度参数下主干支臂长度随着时间步长的变化曲线图,如图2所示。从图2可以看出,当各向异性强度参数从0.0 0 5增加到0.0 0 9时,支臂长度的直线斜率会随之变大。这说明随着各向异性强度参数的增加,枝晶尖端的半径变小,加快了枝晶生长速度。031沈阳师范大学学报(自然科学版)第4 1卷 3.2 热扩散参数对枝晶生长的影响3.2.1 热扩散参数对枝晶形貌的影响为了探究热扩散参数与凝固过程中枝晶形貌之间的关系,选择3个热扩散参数分别为=1.4,=1.6,=1.8下的3个时间步长的模拟数据,并绘制了3组图片(图3)。其中,a,b,c对应着时间步长n=10 0 0,n=20 0 0,n=30 0 0;1,2,3则分别对应每个时间步长下,热扩散参数分别为1.4,1.6和1.8时的枝晶形貌。n=10 0 0n=20 0 0n=30 0 0(a 1)n=10 0 0,=1.4(a 2)n=20 0 0,=1.4(a 3)n=30 0 0,=1.4(b 1)n=10 0 0,=1.6(b 2)n=20 0 0,=1.6(b 3)n=30 0 0,=1.6(c 1)n=10 0 0,=1.8(c 2)n=20 0 0,=1.8(c 3)n=30 0 0,=1.8图3 不同热扩散参数下的枝晶形貌F i g.3 D e n d r i t em o r p h o l o g yu n d e rd i f f e r e n t t h e r m a l d i f f u s i o np a r a m e t e r s图4 不同热扩散参数下枝晶生长速度F i g.4 D e n d r i t eg r o w t h r a t eu n d e rd i f f e r e n tt h e r m a l d i f f u s i o np a r a m e t e r s从图3的横向对比可以看出,随着热扩散参数的增加,枝晶的主干支臂会变小,同时也会抑制二次支臂的生长。这是因为固液界面的前沿部分集中了大量的枝晶生长时释放的潜热,随着热扩散参数的增加,凝固潜热释放也会变快,导致枝晶的主干变得越来越光滑,二次支臂生长也不明显。从图3的纵向对比也可以看出时间步长对枝晶生长的影响。随着时间步长的增加,枝晶主干支臂长度也会随之减小,二次支臂的数量和长度也会减小。3.2.2 热扩散参数对枝晶尖端生长速度的影响图4为不同热扩散参数下枝晶主干支臂131 第2期 高 明,等:相场法模拟中取向与热扩散参数对N i-S n合金枝晶形貌的影响长度随着时间步长的变化曲线图。从图4可以看出,随着热扩散参数的增加,枝晶主干支臂长度随之减少,导致枝晶尖端生长速度变慢。这是由于凝固潜热的释放速度会随热扩散参数增加而变快,导致温度梯度变大,使枝晶尖端生长速度减慢。4 结 论通过对得到的模拟结果进行对比分析,可以得出以下结论:1)各向异性强度参数增大会导致枝晶的尖端生长速度变快,二次支臂更容易出现,其数量也会增加。各向异性强度参数相同时,运行的时间步长越多,枝晶的主干支臂生长越发达。2)热扩散参数增大会导致枝晶的尖端生长速度变慢。当热扩散参数相同时,运行的时间步长越多,枝晶的生长形貌越显著。参考文献:1KUR ZW,F I S HE RDJ.F u n d a m e n t a l o f s o l i d i f i c a t i o n4 t he d i t i o nM.S w i t z e r l a n d:T r a n sT e c hP u b l i c a t i o n,1 9 9 8:1 8 0 1 8 4.2 陈万春,简来成.相场方法及其在晶体生长中的应用J.人工晶体学报,2 0 0 2(3):2 4 5 2 5 9.3KO B AYA S H IR.M o d e l i n ga n dn u m e r i c a ls i m u l a t i o n so fd e n d r i t i cc r y s t a lg r o w t hJ.P h y s i c aD,1 9 9 3,6 3(3/4):4 1 0 4 2 3.4 高明,刘芳慧.相场法在组织模拟中的应用J.沈阳师范大学学报(自然科学版),2 0 1 3,3 1(3):3 1 8 3 2 5.5 龙文元.铝合金凝固过程枝晶生长的相场法数值模拟D.武汉:华中科技大学,2 0 0 4.6 冯力,王智平,路阳,等.二元合金多晶粒的枝晶生长的等温相场模型J.物理学报,2 0 0 8(2):1 0 8 4 1 0 9 0.7 刘芳慧.相场法模拟二元合金枝晶生长D.沈阳:沈阳师范大学,2 0 1 5.8L I U HS,WAN GJ,J I NZP.T h e r m o d y n a m i co p t i m i z a t i o no f t h eN i-S nb i n a r ys y s t e mJ.C a l p h a d,2 0 0 4,2 8(4):3 6 3 3 7 0.9 王智平,肖荣振,朱昌盛,等.N i-C u二元合金定向凝固的相场法数值模拟J.稀有金属材料与工程,2 0 0 6(S 2):3 6 93 7 3.1 0WHE E L E RAA,MUR R AYBT,S CHA E F E RRJ.C o m p u t a t i o no f d e n d r i t e su s i n gap h a s e f i e l dm o d e lJ.P h y s i c aD,1 9 9 3,6 6(1 0):2 4 3 2 6 2.1 1 冯力.多元合金凝固微观组织的相场法模拟研究D.兰州:兰州理工大学,2 0 0 9.1 2 佟乐乐.基于相场法的二元合金凝固过程三维模拟及研究D.太原:中北大学,2 0 1 6.1 3 袁训锋,丁 雨 田.相 场 法 模 拟 强 各 向 异 性 作 用 下 二 元 合 金 枝 晶 生 长 J.中 国 有 色 金 属 学 报,2 0 1 1,2 1(9):2 2 1 6 2 2 2 2.1 4J E ONGJH,G O L D E N F E L DN,D AN T Z I GJA.P h a s e f i e l dm o d e l f o r t h r e ed i m e n s i o n a l d e n d r i t i cg r o w t hw i t hf l u i df l o wJ.P h y sR e vE,2 0 0 1,6 4:0 4 1 6 0 2.1 5 张迪.N i-S n合金枝晶生长的相场法模拟D.太原:中北大学,2 0 2 1.231沈阳师范大学学报(自然科学版)第4 1卷