菲克定律应用

1、菲克定律1扩散动力学方程菲克定律1.1菲克第一定律1.1.1宏观表达式1858年，菲克（Fick）参照了傅里叶（Fourier）于1822年建立的导热方程，建立定量公式。在也时间内，沿X方向通过X处截面所迁移的物质的量侦与X处的浓度梯度成正比:有：Am 氐 CAAtAXdCdm即有：Am 氐 CAAtAXdCdm即dxAdtdCdx（7-1）图7-1扩散过程中溶质原子的根据上式引入扩散通量概念分布式（7-1）即菲克第一定律。式中J称为扩散通量，常用单位是mol/（ cm2 s）;竺浓度梯度；dx.D扩散系数，它表示单位浓度梯度下的 :：通量，单位为cm2/s或m2 /s ；气负号表示扩散方向与

2、浓度梯度方向相 工 反见图7-2。图7-2溶质原子流动1.1.2微观表达式的方向与浓度降低的方 微观模型：向相一致设任选的参考平面1、平面2上扩 散原子面密度分别为n1和n2,若n1=n2，则无净扩散流。假定原子在平衡位置的振动周期为t则一个原子单位时间内离 开相对平衡位置跃迁次数的平均值，即跃迁频率为r =1（7-2）T由于每个坐标轴有正、负两个方向，所以向给定坐标轴正向跃迁的几率是-r。6设由平面l向平面2的跳动原子通量为J量为J12，由平面2向平面1的跳动原图7-3 一维扩散的微观模型子通量为J21 TOC o 1-5 h z J 12= 6 Y(7-3)J 21= 6 n2 r(7-4

3、)注意到正、反两个方向，则通过平面1沿x方向的扩散通量为(7-5)J = J J = r(n n )(7-5)11221612而浓度可表示为(7-6)式(7-6)中的1表示取代单位面积计算e表示沿扩散方向的跳动距离(见图7-3)，则由式(7-5)、式(7-6)得j = -r(c c )5 = -r(c c )5 = -rs2竺=d竺(7-7)1 6126216 dx dx式(7-7 )即菲克第一定律的微观表达式，其中d =1 rs 2(7-8)6式(7-8)反映了扩散系数与晶体结构微观参量之间的关系，是扩散 系数的微观表达式。三维情况下，对于各向同性材料(D相同)，则(7-9)J= J + J

4、 + J = D(i 竺 + j 竺 + k 竺)=DV- C x y zdx dx (7-9)式中：V = iJj 1 + k区为梯度算符。dx dx dx对于各向异性材料，扩散系数D为二阶张量，这时，f D D f D D 气气气D23ID D D ) 3132 33 /办竺办竺dx ）(7-10)对于菲克第一定律，有以下三点值得注意：（1）式（7-1）是唯象的关系式，其中并不涉及扩散系统内部原 子运动的微观过程。（2）扩散系数反映了扩散系统的特性，并不仅仅取决于某一种组元的特性。（3）式（7-1）不仅适用于扩散系统的任何位置，而且适用于扩 散过程的任一时刻。其中，j、d、竺可以是常量，也

5、可以是变量， dx即式（7-1）既可适用于稳态扩散，也可适用于非稳态扩散。1.2菲克第二定律当扩散处于非稳态，即各点的浓度随时间而改变时，利用（7-1）不容易求出C（x/）。但通常的扩散过程大都是非稳态扩散，为便于求 出C（x/），菲克从物质的平衡关系着手，建立了第二个微分方程式。1.2.1 一维扩散则有Am = (J A - J AxA)At图7-4扩散流通过微小体Am _ J - 则有Am = (J A - J AxA)At图7-4扩散流通过微小体Am _ J - J =xx I AxAxAAtAx积的情况当Ax、At时，有告=-J将式（7-1）代入上式得竺工(D竺)dt dx dx(7-

6、11)(7-12(7-12)(7-13)如果扩散系数D与浓度无关，则式（7-11）可写成dC d 2 Cdtdx 2般称式（7-11）、式（7-12）为菲克第二定律。1.2.2三维扩散（1）直角坐标系中竺工(D竺)+ 1 (D竺)+ 1 (D竺)dt dx dx dy dy dz dz当扩散系数与浓度无关，即与空间位置无关时，dC = D( d 2 C + d 2 C + d 2 C) dtdx 2dy 2dz 2(7-14)或简记为：竺=访2C dt(7-15)式中： v2 =F + f 为 Laplace 算符。柱坐标系中通过坐标变换 x = cos，体积元各边为dr，湖，dz，则有:y

7、= r sin 0 TOC o 1-5 h z dC1f aded n ded8C、=(rD )() +(rD )dtr drdrd0 r d0dzdz7-16)对柱对称扩散，且D与浓度无关时有dCD ddC=(r)dtr drdr(7-17)球坐标系中x = r sin 0 cos 甲通过坐标变换y = r sin0 sin平，体积元各边为dr，rd0，z = rcos0r sin0 d中，则有: TOC o 1-5 h z dC1ddC 1 ddC 0d 2 CdC =上任(r2D 竺)+- 与(D sin0 生)+ .dtr 2drdrsin 0 d0d0sin 2 0d甲 2（7-18

8、）对球对称扩散，且D与浓度无关时有：dC D a6C、=(r 2)6t r 2 drdr(7-19)从形式上看，菲克第二定律表 示，在扩散过程中某点浓度随时间 的变化率与浓度分布曲线在该点的 二阶导数成正比。如图7-5从形式上看，菲克第二定律表 示，在扩散过程中某点浓度随时间 的变化率与浓度分布曲线在该点的 二阶导数成正比。如图7-5所示， 若曲线在该点的二阶导数竺大于 dx 20，即曲线为凹形，则该点的浓度会图7-5菲克第一、第二定律的关系随时间的增加而增加，即竺0;若曲线在该点的二阶导数竺小于 dtdx 20，即曲线为凸形，则该点的浓度会随时间的增加而降低，即竺0。 dt而菲克第一定律表示

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

10、低而恒定的压力P/扩散一定时间以 后，金属膜中建立起稳定的浓度分布。图7-6氢对金属膜的一维稳态扩散图7-6氢对金属膜的一维稳态扩散达到稳态扩散时的边界条件：c | x=0=c2C| x=8 WCC2可由热解反应 H2tH+H的平衡常数K确定，根据K的定 义K=产物活度积反应物活度积设氢原子的浓度为C,则艮 口C = JKp = S.Jp（7-20）式（7-20）中S为西佛特（Sievert）定律常数，其物理意义是，当空 间压力p=1MPa时金属表面的溶解浓度。式（7-20）表明，金属表面 气体的溶解浓度与空间压力的平方根成正比。因此，边界条件为：C C 点=S死（7-21）根据稳定扩散条件，

11、有dedt= dedt= （D 竺）=0dxdx所以积分得de _=const = a dxC = ax + b（7-22）式（7-22）表明金属膜中氢原子的浓度为直线分布，其中积分常数a、b由边界条件式（7-21）确定C - C SJa = = 5 侦P1-、p2）Lb = C 2 = S 互将常数a、b值代入式（7-22）得C （x） = S（v P1- 互）x + S.p 2（7-23）单位时间透过面积为A的金属膜的氢气量虹=JA = -DA竺=-DAa = -DAS（、- ,p ）（7-24）dtdx5 i 2由式（7-24）可知，在本例所示一维扩散的情况下，只要保持p1、 p2恒定，

12、膜中任意点的浓度就会保持不变，而且通过任何截面的流 量四、通量J均为相等的常数。dt引入金属的透气率P表示单位厚度金属在单位压差（以MPa为单位）下、单位面积透过的气体流量P = DS（7-25）式中：D为扩散系数，S为气体在金属中的溶解度，则有J = P （汉-互）（7-26）在实际应用中，为了减少氢气的渗漏现象，多采用球形容器、选 用氢的扩散系数及溶解度较小的金属、以及尽量增加容器壁厚等。2.1.2柱对称稳态扩散史密斯（Smith）利用柱对称稳态扩散测定了碳在y铁中的扩散系 数。将长度为L、半径为尸的薄壁铁管在1000C退火，管内及管外 分别通以压力保持恒定的渗碳及脱碳气氛，当时间足够长，

13、管壁内10 / 25L4E 9 8 6 4 2 LLO.O.0-0,氽ttL4E 9 8 6 4 2 LLO.O.0-0,氽tt悒SU图7-7在1000C碳通过溥壁铁从图7-7还可以引出一个重要管的稳态扩散中，碳的浓度各点的碳浓度不再随时间而变，即竺=0时，单位时间内通过管壁的dt碳量m/t为常数，其中m是t时间内流入或流出管壁的碳量，按照通量的定义mJ =2 兀 rLt（7-27） 由菲克第一定律式（7-1）有mdC=D 2r兀 Ltdr或m = D（2兀 Lt ）-Cd In r（7-28）式中m、L、t以及碳沿管壁的径向分布都可以测量，D可以由C对lnr图的斜率确定（见图7-7）。分布的

14、概念：由于m/t为常数，如果D分布不随浓度而变，则dC也应是常d ln r数，C对lnr作图应当是一直线。但实验指出，在浓度高的区域，q d ln r小，D大；而浓度低的区域，E大，D小。由图7-7算出，在1000C， d ln r碳在y铁中的扩散系数为：当碳的质量分数为0.15 %时，D=2.5x10-7cm2/s；当质量分数为 1.4%时，D=7.7x10-7cm2/s。可见 D是浓度的函数，只有当浓度很小时、或浓度差很小时，D才近似为 常数。2.1.3球对称稳态扩散如图7-8所示，有内径为r外径为r2的球壳，若分别维持内表面、外表面的浓度C面、外表面的浓度CC2保持不变则可实现球对称稳态

15、扩散。边界条件E 1边界条件E 1 r=广 C 2由稳态扩散，并利用式（7-19D 8 , 8C、D 8 , 8C、八(r2 一)= 08tr 2 8r8r徂8C得r 2 =const = a8rdC图7-8球壳中可实现球对称稳态扩散解得（解得（7-29）C=-+b r代入边界条件，确定待定常数a, b_ ri r2(C2 C1)r 一 r b = C 2C11 r r求得浓度分布rr (C - C ) C r - CrC (r) = 2+ I1r(r r )r r（7-30）在实际中，往往需要求出单位时间内通过球壳的扩散量业，并利用 dtr2竺=a的关系dr TOC o 1-5 h z dm

16、dC /=JA = 一 D -4兀 r2 = 4 兀Da dtdr/ C C=4 兀Dr r 21(7-31)而不同球面上的扩散通量21(7-31)而不同球面上的扩散通量j = dm = 1 dm = jrrC- CAdt 4兀 r2 dtr 2 r - r（7-32）图7-9过饱和固溶体的析出可见，对球对称稳态扩散来说， 在不同的球面上，剪相同，但j并不相同。dt上述球对称稳态扩散的分析方法对处理固态相变过程中球形晶核的生长速率是很重要的。匚f了_封严如图7-9中的二元相图所示，局以成分为C0的单相a固溶体从高温学爹爹在V-T冷却，进入双相区并在保温。此时会在过饱和固溶体a 中析出成菲克定律

17、应用分为a的P相，与之平衡的a相成分为Cap。在晶核生长初期，设相晶核半径为r，母相在半径为r的球体中成分由C逐渐降为C120a随着时间由七,匕,七变化，浓度分布曲线逐渐变化，相变过程中各相成 分分布如图7-10所示。般说来，这种相变速度较慢，而且涉及的范围较广，因此可将晶核生长过程当作准稳态扩散处理，即在晶核生长初期任何时刻， 浓度分布曲线保持不变。由球对称稳态扩散的分析结果式（7-31）， 并利用r1r2，即新相晶核很小、扩散范围很大的条件。应特别注意 分析的对象是内径为、外径为r2的球壳，由扩散通过球壳的流量 dm，其负值即为新相晶核的生长速率。dt也=-D 4兀rr 二D 4兀r 2

18、MS dt1 2 r - r1 r（7-33）=-D 4兀 r 2 C （7-33）1r1应注意式（7-33）与菲克第一定律的区别，因为式中的C0 -命应注意式（7-33）是浓度梯度。2.2非稳态扩散非稳态扩散方程的解，只能根据所讨论的初始条件和边界条件而定，过程的条件不同方程的解也不同，下面分几种情况加以讨论。7.2.2.1 一维无穷长物体的扩散菲克定律应用无穷长的意义是相对于扩散区长度而言，若一维扩散物体的长度 大于顷，则可按一维无穷长处理。由于固体的扩散系数 D在 10-210-12cm2.s-1很大的范围内变化，因此这里所说的无穷并不等同于 表观无穷长。设A，B是两根成分均匀的等截面金

19、属棒，长度符合上述无穷长 的要求。A的成分是C2, B的成分是C】。将两根金属棒加压焊上， 形成扩散偶。取焊接面为坐标原点，扩散方向沿X方向，扩散偶成 分随时间的变化如图7-11所示。求解的扩散方程为式(7-12)dC 8 2 C 8t 8x 2初始条件 t=0 时，C=C1，(x0)C=C2，(x V 0)(7-35)边界条件 t N 时，C=C1，(x= 8 )C=C2，(x= 8 )(7-36)求解扩散方程的目的在于求出任 何时刻的浓度分布C （x,t）可采用分 离变量法，拉氏变换法，但在式 （7-12），式（7-35），式（7-36）的特 定条件下，采用波耳兹曼变换更为方图7-11扩散

20、偶成分随时间便，即令的变化X E（7-37） 代入式（7-12）左边dCdC凯dC xdC人=, = , = , dt8人dt8人2t3/2d人2t8 2 8 2 CD=D8x 22 +竺在=d竺.1故式（7-12）变成了一个常微 8K 8x 2dK t分方程（7-38）dC令dC令*- = U血代入式（7-38）得（7-39（7-39）（7-40）人 n dU一一u = D-2 肽解得u=亦硕-务） 式（7-40）代入到咎U中，有将上式积分，(7-41)dC d=a、exp(-C 将上式积分，(7-41)dC d=a、exp(-C = a j * exp( -)dk + bo4 D再令p=X

21、 /（2侦万），则式（7-41）可改写为C = a - 2：D jp exp（一。2）d。+ b = ajp exp（一。2）d。+ b（7-42）oo注意式（7-42）是用定积分，即图7-12中斜线所示的面积来表 示的，被积函数为高斯函数 exp（一。2），积分上限为。根据高斯误差积分j%xp（-。2）d。= -2（7-43）图7-12用定积分表示浓度因为p=X /（2,切）=x /（2板瓦），利用边界条件式（7-36）在tNO时，分别有C = C = a j+8 e 邛 2 dp + b C = C2 = a J- e 邛2 dp + b故 C = a + b， C =一。 + b（7-4

22、4）求出积分常数a, （7-44）a三，b =工2 麝2将式（7-44）代入式（7-428-32）有（7-45）C =堂工 -.史 j6 exp（-P 2（7-45）22 部0式（7-45）中的积分函数称为高斯误差函数，用erf （p ）表示（见图7-12），定义为erf (erf (P )= pexp(-P2)dP 疽冗0（7-46）p值对应的erf （P）值列于表7-1。这样式（7-45）可改写成C = Ci ； C2 - C2 一 Ci erf （p ）（ 7-47）式（7-47）即为扩散偶在扩散过程中，溶质浓度随p，即随erf （P）的 变化关系式。（1）式（7-47）的用法给定扩散系

23、统，已知扩散时间t，可求出浓度分布曲线C（x,t）。具体的方法是，查表求出扩散系数D，由D、t以及确定的，求出0=x /(2D)，查表 7-1 求出 erf (P)，代入式(7-47 )求出 C(x,t)。已知某一时刻C(x,t)的曲线，可求出不同浓度下的扩散系数。 具体的方法是，由C(x,t)计算出erf (P )，查表7-1求出P，t、X已知， 利用P f /(2.、而)可求出扩散系数D。任一时刻C(x,t)曲线的特点对于x=0的平面，即原始接触面，有P=0，即erf(P)，因此该 平面的浓度c =。1 ； C2恒定不变；在X = 3，即边界处浓度，有 C = CC =C2，即边界处浓度也

24、恒定不变。曲线斜率竺=竺.亚= _ = e* . ._!(7-48)办邓办 22、：Dt 郁由式(7-47)，式(7-48)可以看出，浓度曲线关于中心(x=0,C = C藐)是对称的。随着时间增加，曲线斜率变小，当t T8时， 2各点浓度都达到C1 + C2，实现了均匀化。2抛物线扩散规律由图7-12及式(7-47)可知，浓度C(x,t)与p有 对应的关系， 由于p= x/(2板瓦),因此C(x,t)与x/5之间也存在一一对应的关系，设 K(C)是决定于浓度C的常数，必有(7-49)x2=K（C）t(7-49)式（7-49 ）称为抛物线扩散规律，其应用范围为不发生相变的扩散。散时间之比为1：

25、4：图散时间之比为1： 4：图7-13抛物线扩散规律如图7-13所示，若等浓度q的扩散等距离之比为1： 2： 3： 4,则所用的扩9： 16。（4）式（7-47）的恒等变形式5可以写成:耘切:匿mtC = C2 + C1 -C2 + C1 -C rf(P) = C 1 -：rf(P)+ Cerf(p)221J0 J1 7图7-14 一维无穷长物体扩散的两种(7-50)特殊情况式中:（a）式中:（a）镀层的扩散、异种金属的扩散当g=0时（镀层的焊扩散，异种金属的扩散焊），（b）扩散，异种金属的扩散焊），（b）真空除气、表面脱碳如图7-14（a），有C = C01 - erf （p ）（7-51）

26、当C0=0时（除气初期，真空除气以及板材的表面脱碳等）， 如图7-14（b），有菲克定律应用C = Cerf (p )(7-52)近似估算由查表 7-1 可知，当 p =0.5 时，erf (p)=0.52040.5,亦即当 x2=Dt 时，根据式(7-51)有60.5CO。由于扩散，如果某处的浓度达到初 始浓度的一半，一般称该处发生了显著扩散。关于显著扩散，利用 X2=Dt，给出x可求t，给出t可求X。2.2.2半无穷长物体的扩散半无穷长物体扩散的特点是，表面浓度保持恒定，而物体的长 度大于顷。对于金属表面的渗碳、渗氮处理来说，金属外表面的 气体浓度就是该温度下相应气体在金属中的饱和溶解度C

27、。，它是恒 定不变的；而对于真空除气来说，表面浓度为0，也是恒定不变的。钢铁渗碳是半无穷长物体扩散的典型实例。例如将工业纯铁在 927 C进行渗碳处理，假定在渗碳炉内工件表面很快就达到碳的饱和 浓度(1.3%C)，而后保持不变，同时碳原子不断地向里扩散。这样， 渗碳层的厚度、渗碳层中的碳浓度和渗碳时间的关系，便可由式 (7-51)求得。初始条件，t=0，xo，C=0；边界条件，tN0，x=8，C=0； x=0，C0=1.3927 C时的碳在铁中扩散系数D=1.5 X 10-7cm2s-1，所以C = 1.31 - erf,x 1.31 - erf1.29 x 103 .=I 2J1.5 x 1

28、0 -7 t )1&)渗碳10（3.6X 104s）后渗碳层中的碳分布C = 1.31 - erf （6.8 尤）在实际生产中，渗碳处理常用于低碳钢，如含碳量为0.25%的钢。这时为了计算的方便，可将碳的浓度坐标移到0.25为原点，这样就 可以采用与工业纯铁同样的计算方法。2.2.3瞬时平面源在单位面积的纯金属表面涂上扩散元素组成平面源，然后对接成扩散偶进行扩散。若扩散系数为常数，其扩散方程为式（7-12）：注意到涂层的厚度为0，因此方程式（7-12）的初始、边界条件为（7-53）当 t 0时，C| = 8，c| （7-53）当t 0时，C 0X=3由微分知识可知，满足方程式（7-12）及上述初始、边界条件的解具有下述形式C-M1C-M11/2exp 4Dt)（7-54）式中a是待定常数。可以利用扩散物质的总量M来求积分常数s(7-55)(7-55)(7-56)(7-57)M =M Cdx3如果浓度分布由式（7-54 ）表示，并令X 2=B 24 Dt则有dx = 2（