【物理应用】基于相场法模拟金属镍的晶粒的长大过程附matlab完整代码
1 内容介绍
金属镍是支撑信息技术产业和光伏产业发展的基础材料,它的性能直接影响和决定着相关高技术产业的发展.随着直拉硅单晶的直径不断增大以及集成电路特性线宽持续减小,金属镍的缺陷控制及其质量保证尤为重要,特别是新型的基于快速热处理(RTP)的"内吸杂"结构在热处理过程中硅片内部演变行为的研究.近些年来,计算机数值模拟技术越来越多地被用于材料微缺陷演变的研究中.其中,以吉布斯-朗道自由能理论为基础的相场法,已经成为研究材料介-微观缺陷的一种重要数值模拟方法.
2 完整代码
%模拟镍的相场计算(无任何耦合、粗网格、变边界条件)-----------------------------
clear;clc;clf;tic
close all
%参数----------------------------------------------------------------------
m=0.05; eps=0.005; r=0.05; %动力学系数、界面相关参数、各向异性系数
Tm=1728; L=2350; DT=0.155; Cp=5.42; %熔点、潜热、热扩散率、比热
Delta=-0.6; R=0.05; alfa=400; %无量纲过冷度、形核半径、耦合相关参数alfa
EPS=eps^2;A=m/EPS;B=(m/EPS)*EPS;C=eps*alfa*(-Delta);D=30/(-Delta);
%网格划分-------------------------------------------------------------------
N=100; NTimeSteps=5000; Dx=0.002; Dy=Dx; Dt=1e-5; E=(m*Dt)/eps^2; DXX=Dx^2;
[x,y]=