coeffd.m

来自「用于计算脉冲激光照射下物体的温度分布」· M 代码 · 共 17 行

function f=CoeffD(m,n,l)
%  计算控制容积界面上的扩散系数，由相邻结点的扩散系数值插值得到，是温度和浓度的函数
%  输入结点的序号m和界面标记号n,时层序号l
%  返回该结点该界面上的扩散系数值
global R Z r z rn zn L1 M1;
i=mod(m,L1);
j=(m-i)/L1+1;
switch n
        case 1%aE
        f=((Z(i)-z(i))/(z(i+1)-z(i))/CoeffP(i,j,l)+(z(i+1)-Z(i))/(z(i+1)-z(i))/CoeffP(i+1,j,l) )^(-1);
        case 2%aW
        f=((Z(i-1)-z(i-1))/(z(i)-z(i-1))/CoeffP(i-1,j,l)+(z(i)-Z(i-1))/(z(i)-z(i-1))/CoeffP(i,j,l) )^(-1);
        case 3%aN
        f=((R(j)-r(j))/(r(j+1)-r(j))/CoeffP(i,j,l)+(r(j+1)-R(j))/(r(j+1)-r(j))/CoeffP(i+1,j,l) )^(-1);
        case 4
        f=((R(j-1)-r(j-1))/(r(j)-r(j-1))/CoeffP(i,j-1,l)+(r(j)-R(j-1))/(r(j)-r(j-1))/CoeffP(i,j,l) )^(-1);
end