菲克定律应用

1、1 扩散动力学方程菲克定律1.1 菲克第一定律宏观表达式1858年，菲克（Fick）参照了傅里叶（Fourier）于1822年建立的导热方程，建立定量公式。在时间内，沿方向通过处截面所迁移的物质的量与处的浓度梯度成正比：图7-1 扩散过程中溶质原子的分布即 根据上式引入扩散通量概念，则有： （7-1）图7-2 溶质原子流动的方向与浓度降低的方向相一致式（7-1）即菲克第一定律。式中称为扩散通量，常用单位是/（；浓度梯度；D扩散系数，它表示单位浓度梯度下的通量，单位为/或；负号表示扩散方向与浓度梯度方向相反见图7-2。1.1.2微观表达式微观模型：设任选的参考平面1、平面2上扩散原子面密度分别为

2、n1和n2，若n1n2，则无净扩散流。假定原子在平衡位置的振动周期为，则一个原子单位时间内离开相对平衡位置跃迁次数的平均值，即跃迁频率G为图7-3 一维扩散的微观模型 (7-2) 由于每个坐标轴有正、负两个方向，所以向给定坐标轴正向跃迁的几率是。设由平面l向平面2的跳动原子通量为J12，由平面2向平面1的跳动原子通量为J21 (7-3) (7-4)注意到正、反两个方向，则通过平面1沿x方向的扩散通量为 (7-5)而浓度可表示为 (7-6)式(7-6)中的1表示取代单位面积计算，表示沿扩散方向的跳动距离（见图7-3），则由式(7-5)、式(7-6)得 (7-7)式(7-7)即菲克第一定律的微观表

3、达式，其中 (7-8)式（7-8）反映了扩散系数与晶体结构微观参量之间的关系，是扩散系数的微观表达式。三维情况下，对于各向同性材料（D相同），则 （7-9）式中：为梯度算符。对于各向异性材料，扩散系数D为二阶张量，这时， （7-10）对于菲克第一定律，有以下三点值得注意：（1）式（7-1）是唯象的关系式，其中并不涉及扩散系统内部原子运动的微观过程。（2）扩散系数反映了扩散系统的特性，并不仅仅取决于某一种组元的特性。（3）式（7-1）不仅适用于扩散系统的任何位置，而且适用于扩散过程的任一时刻。其中，、可以是常量，也可以是变量，即式（7-1）既可适用于稳态扩散，也可适用于非稳态扩散。1.2 菲克第

4、二定律当扩散处于非稳态，即各点的浓度随时间而改变时，利用式（7-1）不容易求出C(x,t)。但通常的扩散过程大都是非稳态扩散，为便于求出C(x,t)，菲克从物质的平衡关系着手，建立了第二个微分方程式。 一维扩散图7-4 扩散流通过微小体积的情况如图7-4所示，在扩散方向上取体积元，Jx和分别表示流入体积元及流出体积元的扩散通量，则在时间内，体积元中扩散物质的积累量为则有 当、0时，有 将式（7-1）代入上式得 （7-11）如果扩散系数与浓度无关，则式（7-11）可写成 （7-12）一般称式（7-11）、式（7-12）为菲克第二定律。1.2.2 三维扩散（1）直角坐标系中 （7-13）当扩散系数

5、与浓度无关，即与空间位置无关时， （7-14）或简记为： （7-15）式中：为Laplace算符。（2）柱坐标系中通过坐标变换 ，体积元各边为，则有： （7-16）对柱对称扩散，且与浓度无关时有 （7-17）（3）球坐标系中通过坐标变换 ，体积元各边为，则有： （7-18）对球对称扩散，且与浓度无关时有： （7-19）图7-5 菲克第一、第二定律的关系从形式上看，菲克第二定律表示，在扩散过程中某点浓度随时间的变化率与浓度分布曲线在该点的二阶导数成正比。如图7-5所示，若曲线在该点的二阶导数大于0，即曲线为凹形，则该点的浓度会随时间的增加而增加，即0；若曲线在该点的二阶导数小于0，即曲线为凸形，

6、则该点的浓度会随时间的增加而降低，即0。而菲克第一定律表示扩散方向与浓度降低的方向相一致。从上述意义讲菲克第一、第二定律本质上是一个定律，均表明扩散的结果总是使不均匀体系均匀化，由非平衡逐渐达到平衡。2 菲克定律的应用涉及扩散的实际问题有两类:其一是求解通过某一曲面（如平面、柱面、球面等）的通量J，以解决单位时间通过该面的物质流量；其二是求解浓度分布C(x,t)，以解决材料的组分及显微结构控制，为此需要分别求解菲克第一定律及菲克第二定律。2.1 稳态扩散及其应用 一维稳态扩散 考虑氢通过金属膜的扩散。如图7-6所示，金属膜的厚度为d，取x轴垂直于膜面。考虑金属膜两边供气与抽气同时进行，一面保持

7、高而恒定的压力p2，另一面保持低而恒定的压力p1。扩散一定时间以后，金属膜中建立起稳定的浓度分布。图7-6 氢对金属膜的一维稳态扩散 氢的扩散包括氢气吸附于金属膜表面，氢分子分解为原子、离子，以及氢离子在金属膜中的扩散等过程。达到稳态扩散时的边界条件： C| x=0 =C2C| x=d =C1C1、C2可由热解反应 H2®H+H的平衡常数K确定，根据K的定义K=设氢原子的浓度为C，则K=即 （7-20）式（7-20）中S为西佛特（Sievert）定律常数，其物理意义是，当空间压力p=1MPa时金属表面的溶解浓度。式（7-20）表明，金属表面气体的溶解浓度与空间压力的平方根成正比。因此

8、，边界条件为： C| x=0 =S C| x=d =S （7-21）根据稳定扩散条件，有=(D)=0所以 =积分得 （7-22）式（7-22）表明金属膜中氢原子的浓度为直线分布，其中积分常数a、b由边界条件式（7-21）确定 将常数a、b值代入式（7-22）得 （7-23）单位时间透过面积为A的金属膜的氢气量 （7-24）由式（7-24）可知，在本例所示一维扩散的情况下，只要保持p1、p2恒定，膜中任意点的浓度就会保持不变，而且通过任何截面的流量、通量J均为相等的常数。 引入金属的透气率P表示单位厚度金属在单位压差（以MPa为单位）下、单位面积透过的气体流量 （7-25）式中：D为扩散系数，S

9、为气体在金属中的溶解度，则有 （7-26）在实际应用中，为了减少氢气的渗漏现象，多采用球形容器、选用氢的扩散系数及溶解度较小的金属、以及尽量增加容器壁厚等。2.1.2 柱对称稳态扩散史密斯（Smith）利用柱对称稳态扩散测定了碳在g铁中的扩散系数。将长度为L、半径为r的薄壁铁管在1000退火，管内及管外分别通以压力保持恒定的渗碳及脱碳气氛，当时间足够长，管壁内各点的碳浓度不再随时间而变，即时，单位时间内通过管壁的碳量m/t为常数，其中m是t时间内流入或流出管壁的碳量，按照通量的定义 （7-27）由菲克第一定律式（7-1）有图7-7 在1000碳通过薄壁铁管的稳态扩散中，碳的浓度分布 或 （7-

10、28）式中m、L、t以及碳沿管壁的径向分布都可以测量，D可以由C对lnr图的斜率确定（见图7-7）。从图7-7还可以引出一个重要的概念：由于m/t为常数，如果D不随浓度而变，则也应是常数，C对lnr作图应当是一直线。但实验指出，在浓度高的区域，小，D大；而浓度低的区域，大，D小。由图7-7算出，在1000，碳在g铁中的扩散系数为：当碳的质量分数为0.15时，D=2.5´10-7cm2/s；当质量分数为1.4时，D=7.7´10-7cm2/s。可见D是浓度的函数，只有当浓度很小时、或浓度差很小时，D才近似为常数。2.1.3 球对称稳态扩散如图7-8所示，有内径为r1、外径为r

11、2的球壳，若分别维持内表面、外表面的浓度C1、C2保持不变，则可实现球对称稳态扩散。图7-8 球壳中可实现球对称稳态扩散边界条件 |由稳态扩散，并利用式（7-19）得 解得 （7-29）代入边界条件，确定待定常数求得浓度分布 （7-30）在实际中，往往需要求出单位时间内通过球壳的扩散量，并利用的关系图7-9 过饱和固溶体的析出 (7-31)而不同球面上的扩散通量 （7-32）可见，对球对称稳态扩散来说，在不同的球面上，相同，但并不相同。图7-10 球形晶核的生长过程上述球对称稳态扩散的分析方法对处理固态相变过程中球形晶核的生长速率是很重要的。如图7-9中的二元相图所示，成分为的单相固溶体从高温

12、冷却，进入双相区并在保温。此时会在过饱和固溶体中析出成分为的相，与之平衡的相成分为。在晶核生长初期，设相晶核半径为，母相在半径为的球体中成分由逐渐降为，随着时间由变化，浓度分布曲线逐渐变化，相变过程中各相成分分布如图7-10所示。一般说来，这种相变速度较慢，而且涉及的范围较广，因此可将晶核生长过程当作准稳态扩散处理，即在晶核生长初期任何时刻，浓度分布曲线保持不变。由球对称稳态扩散的分析结果式（7-31），并利用>>，即新相晶核很小、扩散范围很大的条件。应特别注意分析的对象是内径为、外径为的球壳，由扩散通过球壳的流量，其负值即为新相晶核的生长速率。 = （7-33）应注意式（7-33

13、）与菲克第一定律的区别，因为式中的并不是浓度梯度。2.2 非稳态扩散非稳态扩散方程的解，只能根据所讨论的初始条件和边界条件而定，过程的条件不同方程的解也不同，下面分几种情况加以讨论。.1 一维无穷长物体的扩散无穷长的意义是相对于扩散区长度而言，若一维扩散物体的长度大于，则可按一维无穷长处理。由于固体的扩散系数D在10-210-12cm2×s-1很大的范围内变化，因此这里所说的无穷并不等同于表观无穷长。设A，B是两根成分均匀的等截面金属棒，长度符合上述无穷长的要求。A的成分是C2，B的成分是C1。将两根金属棒加压焊上，形成扩散偶。取焊接面为坐标原点，扩散方向沿X方向，扩散偶成分随时间的

14、变化如图7-11所示。求解的扩散方程为式(7-12) 初始条件 t=0时，C=C1，(x0) C=C2，(x0) (7-35)边界条件 t时，C=C1，(x=) C=C2，(x=) （7-36）图7-11 扩散偶成分随时间的变化求解扩散方程的目的在于求出任何时刻的浓度分布C（x,t）可采用分离变量法，拉氏变换法，但在式（7-12），式（7-35），式（7-36）的特定条件下，采用波耳兹曼变换更为方便，即令 （7-37）代入式（7-12）左边右边故式（7-12）变成了一个常微分方程 （7-38）令，代入式（7-38）得 （7-39）解得 （7-40）式（7-40）代入到中，有将上式积分， （7-

15、41）再令，则式（7-41）可改写为 （7-42）图7-12 用定积分表示浓度注意式（7-42）是用定积分，即图7-12中斜线所示的面积来表示的，被积函数为高斯函数，积分上限为。根据高斯误差积分 （7-43）因为)，利用边界条件式（7-36）在t0时，分别有故 ，求出积分常数a，b分别为， （7-44）将式（7-44）代入式（7-428-32）有 （7-45）式（7-45）中的积分函数称为高斯误差函数，用表示（见图7-12），定义为= （7-46）值对应的值列于表7-1。这样式（7-45）可改写成 （7-47）式（7-47）即为扩散偶在扩散过程中，溶质浓度随，即随的变化关系式。（1）式（7-4

16、7）的用法 给定扩散系统，已知扩散时间t，可求出浓度分布曲线C(x,t)。具体的方法是，查表求出扩散系数D，由D、t以及确定的，求出，查表7-1求出，代入式（7-47）求出C(x,t)。 已知某一时刻C(x,t)的曲线，可求出不同浓度下的扩散系数。具体的方法是，由C(x,t)计算出，查表7-1求出，t、x已知，利用可求出扩散系数D。（2）任一时刻C(x,t)曲线的特点 对于x=0的平面，即原始接触面，有=0，即,因此该平面的浓度恒定不变；在，即边界处浓度，有，即边界处浓度也恒定不变。 曲线斜率 （7-48）由式（7-47），式（7-48）可以看出，浓度曲线关于中心（x=0, ）是对称的。随着时

17、间增加，曲线斜率变小，当时，各点浓度都达到，实现了均匀化。（3）抛物线扩散规律由图7-12及式（7-47）可知，浓度C(x,t)与有一一对应的关系，由于,因此C(x,t)与之间也存在一一对应的关系，设K(C)是决定于浓度C的常数，必有x2=K(C)t (7-49)图7-13 抛物线扩散规律式（7-49）称为抛物线扩散规律，其应用范围为不发生相变的扩散。如图7-13所示，若等浓度C1的扩散等距离之比为1：2：3：4，则所用的扩散时间之比为1：4：9：16。图7-14 一维无穷长物体扩散的两种特殊情况（a）镀层的扩散、异种金属的扩散焊 (b)真空除气、表面脱碳（4）式（7-47）的恒等变形式（7-

18、47）可以写成 （7-50）式中：。 当C1=0时（镀层的扩散，异种金属的扩散焊），如图7-14（a），有 （7-51） 当C0=0时（除气初期，真空除气以及板材的表面脱碳等），如图7-14（b），有 （7-52）（5）近似估算由查表7-1可知，当=0.5时，=0.52040.5，亦即当x2=Dt时，根据式（7-51）有C0.5C0。由于扩散，如果某处的浓度达到初始浓度的一半，一般称该处发生了显著扩散。关于显著扩散，利用x2=Dt，给出x可求t，给出t可求x。2.2.2 半无穷长物体的扩散半无穷长物体扩散的特点是，表面浓度保持恒定，而物体的长度大于。对于金属表面的渗碳、渗氮处理来说，金属外表面

19、的气体浓度就是该温度下相应气体在金属中的饱和溶解度C0，它是恒定不变的；而对于真空除气来说，表面浓度为0，也是恒定不变的。钢铁渗碳是半无穷长物体扩散的典型实例。例如将工业纯铁在927进行渗碳处理，假定在渗碳炉内工件表面很快就达到碳的饱和浓度（1.3%C），而后保持不变，同时碳原子不断地向里扩散。这样，渗碳层的厚度、渗碳层中的碳浓度和渗碳时间的关系，便可由式（7-51）求得。初始条件，t=0，xo，C=0；边界条件，t0，x=，C=0；x=0，C0=1.3927时的碳在铁中扩散系数D=1.5×10-7cm2×s-1，所以渗碳10h(3.6×104s)后渗碳层中的碳分

20、布在实际生产中，渗碳处理常用于低碳钢，如含碳量为0.25%的钢。这时为了计算的方便，可将碳的浓度坐标移到0.25为原点，这样就可以采用与工业纯铁同样的计算方法。2.2.3 瞬时平面源在单位面积的纯金属表面涂上扩散元素组成平面源，然后对接成扩散偶进行扩散。若扩散系数为常数，其扩散方程为式（7-12）：注意到涂层的厚度为0，因此方程式（7-12）的初始、边界条件为 （7-53）由微分知识可知，满足方程式（7-12）及上述初始、边界条件的解具有下述形式 （7-54）式中a是待定常数。可以利用扩散物质的总量M来求积分常数a，有 （7-55）如果浓度分布由式（7-54）表示，并令 （7-56）则有，将其代入式（7-55）得将上式代入式（7-54）可得 （7-57）图7-15示出了不同Dt值时由式（7-57）确定的浓度分布曲线