壳核组织形成过程的数值模拟研究.pdfVIP

中国科学G辑：物理学 力学 天文学 2007年第37卷第3期：409—416 《中国科学》杂志社 http：／／www．scichina．com SCIENCEINCHINAPRESS 壳核组织形成过程的数值模拟研究 秦涛王海鹏魏炳波4 (西北工业大学应用物理系，西安710072) 摘要 对空间环境条件下偏晶凝固过程中壳核组织的形成机制进行了数值模 拟研究．在总自由能中引入表面自由能项对H模型进行了修正，使之能够描述表 面偏析作用下的相分离．选取Fe—Cu合金作为模型体系，系统模拟了微重力状态 下深过冷合金液滴的亚稳相分离过程，揭示了表面偏析效应和温度梯度引起的 织从弥散结构演化为3层壳核组织；Fe50Cu50合金经历了“双连续相_÷4层壳核_3 层壳核”的组织演化；Fe35Cu65合金的组织由弥散结构向两层壳核组织演化．偏析 作用使得最外层总是富cu相，温度梯度引起的Marangoni对流决定了内部体相的 组织演化过程． 关键词 空间环境Fe．Cu合金相分离壳核组织计算模拟 偏晶合金的液态相分离一直是材料物理和空间科学领域中的重要研究课题【l~3】．目前的研 离过程中母相分解和组织演化的重要影响因素拉’4J．因为金属的不透明性和苛刻的实验要求， 致使深入的理论研究困难较大，所以对其演化过程和组织形成的微观机制还不十分清楚．数 值模拟作为相分离研究的一种有效方法，近年来被广泛应用于有机聚合物体系相分离的研究 中15,61，但在金属体系中应用的报道还不多，而且当前相分离的模拟工作仅局限于矩形区域和 均匀温度场．在空间环境条件下的液滴具有高度的球对称性和非均匀的温度场【1，7】，采用传统 方法进行模拟将受到限制． H模型是基于黏流体系扩散过程并考虑流场的作用，将扩散方程和N—S方程相耦合而得 到鲥81．它成功地描述了二元液态体系的体相分离过程，但用以研究表面偏析作用下相分离的 报道还很少． Fe—Cu包晶合金是一个重要的亚稳偏晶体系，其平衡相图一1如图1．当过冷至不混溶温度 以下时，合金熔体分离为富Fe液相和富Cu液相．由于相分离发生于热力学亚稳态，过冷熔体 中很容易发生7Fe的形核和生长，从而形成不同特征的组织形貌{7，m】．虽然对这一包晶合金体 收稿日期：2007．02．12；接受日期：2007—05．28 十联系人，E．mail：bbwei@nwpu．edu．cn 410 中国科学G辑物理学力学天文学 第37卷 系的亚稳相分离已进行初步研究，其相分离和组织演化机制还有待进一步探索． 本文采用数值模拟方法研究了空间环境条件下深过冷Fe—Cu合金液滴中的组织演化规律， 分析了表面偏析和温度梯度的存在对组织形貌的作用机制．模拟过程所选取的3种合金成分在 图1中标出． 图1三种合金成分在Fe—Cu二兀相图中的位置 1数学模型与计算方法 H模型描述了浓度场和流场耦合情况下的相分离组织演化，其无量纲化形式为…1 芸一v．㈣叫则卅V等|+v．孝， ㈣ 变分算子；f为随机Gauss白噪声，其涨落幅值为F；F为体系的自由能泛函． 以往的模拟研究仅考虑了体相自由能凡和界面浓度梯度引起的自由能，盯。。两项对总自由 能的贡献，没有涉及表面自由能凡．事实上，表面自由能砖在相分离过程中起着非常重要的 作用．这里将氏引入到总自由能中，可得 F=民+‰+Fs， (2) 为了简化计算而又不失一般性，这里选取R为正规溶液自由能的形式： 吒=瓦o+0[0tn痧+(1一妒)lIl(1一妒)]+29t(1一痧)，(3) 其中Ro为常量部分；9为约化温度，定义为丌孔．心。d表示为 ％ad=寺lV∥，