DOI:10.19476/j.ysxb.1004.0609.2001.s2.009
液态金属快速凝固过程的PFM模型
王君 孙加林
昆明理工大学材料与冶金学院
昆明理工大学材料与冶金学院 昆明650093
摘 要:
近年来发展起来的PFM模型根据“溶质捕捉”和“溶质牵引”两种机制对金属及合金快速凝固过程进行数学描述及模拟 , 较好地反映了金属凝固的物理本质。介绍了PFM模型发展过程、研究现状及理论模型 , 并对发展前景进行了展望。
关键词:
相区场模型 ;溶质捕捉 ;合金凝固 ;枝晶生长 ;
中图分类号: TG111.4
收稿日期: 2001-03-06
Phase-field model of rapid solidification process of liquid alloy
Abstract:
In recent years, phase field model (PFM) of alloy solidification based on solute trapping and solute drag was introduced as a powerful tool to simulate the process of solidification. An overview introduction about its development was given and an emphasis on its theoretical frame was put.
Keyword:
phase field model; solute trapping; alloy solidification; dendritic growth;
Received: 2001-03-06
凝固过程的物理本质及数学模拟一直是凝聚态物理和材料科学领域的研究重点之一。 在平衡凝固和近平衡凝固条件下, 经典热力学和扩散动力学可以较好地解决问题
[1 ,2 ]
。 但在快速凝固过程中, 金属液的冷却速度很快, 不仅固相和液相内部的溶质来不及充分扩散, 凝固界面上的溶质迁移也偏离平衡。 合理地反映界面迁移过程中的物质、 能量传输过程是数学模型必须考虑的问题。
对凝固前沿固液界面的处理大体可以归为两类
[3 ,4 ,5 ,6 ,7 ,8 ,9 ]
: 不连续界面模型 (Sharp-interface model)
[10 ]
和扩散界面模型 (Diffusive-interface model) 。 不连续界面模型针对凝固界面两侧的微小区域进行研究, 认为该区域内由独立的固相和液相组成, 两者之间的界面为简单的几何界面。 该类模型的界面生长速度和液相浓度之间的关系, 较好地反映了溶质传输过程。 但该类模型存在两个明显的缺点: 1) 忽略了界面前沿实际存在的固液共存区, 使结果和实际情况不能相符合; 2) 将界面温度和液相浓度用不同的控制方程来表示, 而不是有机的结合起来。
扩散界面模型认为液相和固相之间存在具有一定厚度的界面, 该界面区域在浓度、 密度和温度等性质上和毗邻的液相和固相区都有着一定的差别, 凝固组织前沿的生长和该界面向液相内部的推移是同时进行的。 同时, 固相和液相之间的传输过程也必须经过这个区域。 通过Ginzburg-Landau相变理论
[11 ,12 ]
发展起来的相区场模型 (Phase-field model) 是其中重要的代表
[13 ]
。
根据对界面微区自由能密度表达方式的不同, 相区场模型又可以分为3类: 由Wheeler, Bottinger和Mcfadden建立的WBM模型
[13 ,14 ,15 ,16 ]
, Steinbach等人建立的适用于稀溶液的模型
[17 ,18 ]
, 以及Losert等人建立的模型
[19 ,20 ]
。 其中以WBM模型的应用最为广泛
[13 ,14 ,21 ,22 ]
, 另外两种则由于有过多的假设条件而缺乏普遍性
[23 ]
。 本文将以WBM模型为例说明相区场模型中的核心问题。
1数学模型
1.1PFM模型
基于Helmholtz自由能公式, Caginalp和Xie给出的WBM1
[13 ,24 ,25 ,26 ,27 ,28 ]
模型具有如下形式:
F
=
∫
Ω
[
f
(
?
,
c
,
Τ
)
+
ε
2
2
|
?
?
2
|
]
d
v
?
?
?
(
1
)
式中 Ω 为研究对象所处的微区, 梯度项反映了界面迁移过程中能量的变化, 常数ε 为梯度能量系数。 二元合金自由能密度函数f (?, c , T ) 由相区场函数?、 浓度场函数c 以及温度场函数T 共同决定:
f
(
?
,
c
,
Τ
)
=
p
(
?
)
f
S
(
c
,
Τ
)
+
[
1
-
p
(
?
)
]
f
L
(
c
,
Τ
)
+
W
(
c
)
4
g
(
?
)
?
?
?
(
2
)
相区场函数? (x , t ) 用以表征凝固界面中每一点相的特征
[21 ]
。 其中的双阱函数g (?) =?2 (1-?) 2 , 在?=0和?=1处有极值, 分别对应于纯液相和固相。 p (?) 满足p (0) =0和p (1) =1, 用来表明凝固界面的微区中由于温度和浓度不同而引起的差异, 具有不同的多项式表达方式
[29 ,30 ]
。 W (c ) 是和表面能以及界面厚度有关的能垒大小。
为使Helmholtz自由能函数随时间变化单调递减, 可令??/?t =-M 1 ?F /??, ?c /?t =?·M 2 ·[c (1-c ) ]? (?f /?c ) ]。 式中常数M 1 为界面动力学系数, 迁移系数M 2 为和溶质扩散系数有关的常数。
通过上述方程, 可以写出如下的界面生长速度v 的控制方程表达式, 这里的z 满足-l /2<z <l /2。
-
v
Μ
1
d
?
d
z
=
ε
2
d
2
?
d
z
2
-
?
f
(
?
,
c
,
Τ
)
?
?
?
?
?
(
3
)
-
v
V
m
d
c
d
z
=
d
J
d
z
?
?
?
(
4
)
1.1.1 溶质分配驱动力
对式 (4) 进行积分可得v (c -c S ) V m =J , 扩散通量J =-c (1-c ) M 2 V -1 m dμ /dz , 界面层化学位具有如下的形式, 积分后便可以得到溶质分配驱动力的表达方式:
d
μ
d
z
=
-
v
m
J
c
(
1
-
c
)
Μ
2
=
-
v
(
c
-
c
S
)
c
(
1
-
c
)
Μ
2
?
?
?
(
5
)
μ
(
c
L
)
-
μ
(
c
S
)
=
-
V
∫
-
l
/
2
l
/
2
c
-
c
S
c
(
1
-
c
)
Μ
2
d
z
?
?
?
(
6
)
1.1.2 相变驱动力
对式 (3) 两边同乘d?/dz 并进行积分, 得到
-
v
Μ
1
∫
-
l
/
2
l
/
2
?
z
2
d
z
=
ε
2
[
?
z
2
]
-
l
/
2
l
/
2
-
∫
-
l
/
2
l
/
2
f
?
?
z
d
z
?
?
?
(
7
)
分步积分得到
∫
-
l
/
2
l
/
2
f ? ?z dz =[f (?, c ) ]
-
l
/
2
l
/
2
-∫
-
l
/
2
l
/
2
f c ?z dz
=[f (?, c ) ]
-
l
/
2
l
/
2
-∫
-
l
/
2
l
/
2
f c (c -c S ) z dz
=[f - (c -c S ) f c ]
-
l
/
2
l
/
2
+∫
-
l
/
2
l
/
2
(c -c S ) μ z dz
由μ =?f /?c 结合 (6) 式上式可以重写为
∫
-
l
/
2
l
/
2
f
?
?
z
d
z
=
[
f
-
(
c
-
c
S
)
f
c
]
-
l
/
2
l
/
2
+
V
m
v
∫
-
l
/
2
l
/
2
J
μ
z
d
z
?
?
?
(
8
)
令a =∫
-
l
/
2
l
/
2
?
z
2
d z, 同时
ε
2
[
?
z
2
]
-
l
/
2
l
/
2
≈
0
, 则可以得到相变时自由能变化的表达式, 亦及相变的驱动力:
ΔF S =f (c S ) -[f (c L ) + (c S -c L ) f c (c L ) ] (9)
1.2溶质捕捉模型 [21,31,32,33]
二元合金的溶质模型是通过对生长前沿浓度变化的分析而建立起来的。 为此, 需要对扩散和凝固所需要的能量进行分析。 设固相和液相内的自由能密度分别为f S (c S , T ) 和f L (c L , T ) , 其中的T 代表温度, c S 和c L 分别代表固相和液相中溶质的摩尔分数。 在液相和固相中A组元和B组元的化学位可以表示如下:
μ
A
L
=
f
L
(
c
L
,
Τ
)
-
c
L
?
f
L
?
c
L
(
c
L
,
Τ
)
?
?
?
(
1
0
)
μ
B
L
=
f
L
(
c
L
,
Τ
)
+
(
1
-
c
L
)
?
f
L
?
c
L
(
c
L
,
Τ
)
?
?
?
(
1
1
)
μ
A
S
=
f
S
(
c
S
,
Τ
)
-
c
S
?
f
S
?
c
S
(
c
S
,
Τ
)
?
?
?
(
1
2
)
μ
B
S
=
f
S
(
c
S
,
Τ
)
-
(
1
-
c
S
)
?
f
S
?
c
S
(
c
S
,
Τ
)
?
?
?
(
1
3
)
热力学平衡时界面处的化学位相等, μ
A
S
=μ
A
L
且μ
B
S
=μ
B
L
。 可以通过互扩散化学势来表示平衡条件: μ L =μ
B
L
-μ
A
L
=?f L /?c L , μ S =μ
B
S
-μ
A
S
=?f S /?c S 。 凝固时的自由能变化可以表示为
Δ
F
S
=
f
S
(
c
S
,
Τ
)
-
[
f
L
(
c
L
,
Τ
)
+
(
c
S
-
c
L
)
?
f
L
?
c
L
(
c
L
,
Τ
)
]
?
?
?
(
1
4
)
该式表明了液相向固相传质的能量变化情况。 平衡时, μ L -μ S =0, ΔF S =0。 在稀溶液条件下, 通过μ L -μ S =0可以推导出c S =k E c L , 并得到平衡分配系数k E 的定义式。 通过ΔF S =0推导出液相温度的表达式T =T M +m L c L 。 在非平衡的条件下, μ L -μ S ≠0, ΔF S ≠0, 两式分别可以理解为溶质分配、 相变的驱动力。
1.2.1 分配系数
假设凝固界面厚度的尺度为原子尺寸, 同时为简化表达, 添加热力学因子c L (1-c L ) , 则扩散通量J 可以表示为如下形式:
J
=
-
Μ
0
V
m
c
L
(
1
-
c
L
)
μ
L
-
μ
S
l
=
v
V
m
(
c
L
-
c
S
)
?
?
?
(
1
5
)
式中 V m 为摩尔体积, M 0 为与浓度相关的迁移系数, v 为界面生长速率。 则界面生长速度和溶质分配驱动力有如下关系:
μ
L
-
μ
S
=
-
v
(
c
L
-
c
S
)
l
Μ
0
c
L
(
1
-
c
L
)
?
?
?
(
1
6
)
在稀溶液中, 分配系数可以表达为ln (k /k E ) =v (1-k ) /v D 。
1.2.2 界面温度
通过凝固自由能变化ΔF S 可以获得界面温度的表达式。 自由能变化和界面迁移速率具有如下的关系
[21 ]
,
v
=
v
c
[
1
-
exp
(
V
m
Δ
F
S
R
Τ
)
]
≈
-
v
c
V
m
Δ
F
S
R
Τ
?
?
?
(
1
7
)
该模型假设凝固过程中不存在能量损失, 其中v c 为可能的最大生长速度。 在溶质捕获模型中, 由于是开放系统, 存在着能量的耗散ΔF D , 则用于凝固过程的自由能为ΔF C =ΔF S -ΔF D , 式 (17) 相应地变为
v
=
v
c
[
1
-
exp
(
V
m
Δ
F
C
R
Τ
)
]
≈
-
v
c
V
m
(
Δ
F
S
-
Δ
F
D
)
R
Τ
?
?
?
(
1
8
)
Aziz和Kaplan专门提出耗散自由能的表达式, 并由此推出稀溶液条件下界面温度的表达式
[34 ,35 ]
。 其中的参数α 可以为溶质从液相向固相浓度转变的度量, α =1时表示存在溶质牵引效应, α =0则不存在溶质牵引。 m L 为平衡液相线斜率, 动力学系数u =v C L / (RT 2 M ) 。
ΔF D = (c L -c S ) {[μ
A
S
(c S , T ) -μ L A (c L , T ) ]-
[μ
B
S
(c S , T ) -μ
B
L
(c L , T ) ]} (19)
Τ
=
Τ
Μ
+
m
L
c
∞
k
?
{
1
-
k
+
[
k
+
(
1
-
k
)
α
]
ln
(
k
/
k
E
)
1
-
k
E
}
-
v
u
?
?
?
(
2
0
)
2应用及展望
针对纯金属材料的平衡及非平衡过程应用PFM模型已经有多年的历史, 随着计算机计算能力的不断提高, 已经能够模拟诸如树枝晶生长等复杂的动态过程
[36 ,37 ,38 ,39 ,40 ,41 ,42 ]
。
针对二元合金凝固过程的研究是WBM1模型建立以后的事情
[13 ,21 ,37 ]
, Xie和Caginalp将PFM模型应用于二元合金的非平衡凝固过程。 此后, 针对二元合金非平衡凝固组织生长的研究越来越多
[14 ,37 ,38 ,39 ,40 ,41 ,42 ,43 ,44 ,45 ,46 ]
, 并和实验结果取得了较好的一致
[21 ,47 ]
。
随着研究的不断深入, 针对更深层次现象的研究也越来越多。 Beckermann等人研究了凝固过程的对流条件影响以及微观偏聚现象
[48 ,49 ,50 ,51 ]
, Conti等人针对PFM模型中的热现象进行了研究
[19 ,52 ]
。 合理借鉴物理学和其它领域的一些理论, 可以帮助完善PFM模型。 如Cahn-Hilliard方程对形核 (调幅分解前期) 的处理非常成功
[53 ,54 ,55 ,56 ,57 ,58 ,59 ]
, 而PFM模型对形核的描述并不充分。 PFM模型还不能解决由于晶格类型不同而导致的枝晶生长各向异性问题, 而固态相变弹性场理论对各向异性弹性场的处理则较为成熟
[60 ,61 ]
。 此外, 枝晶生长作为一种典型的自组织现象, 自组织生长的DLA动力学模型
[62 ,63 ,64 ,65 ]
, 也有一定的参考价值。
参考文献
[1] ZHOUYaoHe (周尧和) , HUZhuangQi (胡壮麒) .So lidificationTechnology (凝固技术) [M ].Beijing:Me chanicalIndustryPress, 1998.10.
[2] FifePC .Modelsforphaseseparationandtheirmathe matics[J].ElectronicJournalofDifferentialEquations, 2000, 2000 (48) :1-26.
[3] CaginalpG , XieW .Phase fieldandsharp interfacealloymodels[J].PhysicalReviewE , 1993, 48:1897-1909.
[4] CaginalpG .Phasefieldmodelsandsharpinterfacelim its:somedifferencesinsubtlesituations[J].RockyMountainJournalofMathematics, 1991, 21 (2) :603-616.
[5] FabbriM , VollerVR .Thephase fieldinthesharpinter facelimit:acomparisonbetweenmodelparameters[J].JCompPhys, 1997, 130:256, .
[6] KesslerDA , KoplikJ , LevineH .Patternselectioninfingeredgrowthphenomena[J].AdvancesinPhysics, 1988, 37:255-339.
[7] LangerJS .Instabilitiesandpatternformationincrystalgrowth[J].RevModPhys, 1980, 52:1-28.
[8] TaylorJE , CahnJW .Linkanisotropicsharpanddiffusesurfacemotionlawsviagradientflows[J].JStatPhys, 1992, 77:183-197.
[9] VladimrovaN , MalagoliA , MauriR .Diffusiondrivenphaseseparationofdeeplyquenchedmixtures[J].PhysRewE , 1998, 58 (6) :7691-7699.
[10] CharachC .ZaltzmanB .Interfacialkineticseffectinplanarsolidificationproblemswithoutinitialundercool ing[J].MathematicalModelsandMethodsinApplSci, 1994, 4:331-354.
[11] ColletP .Thermodynamiclimitoftheginzburg landauequation[J].Nonlinearity, 1994, 7:1175-1190.
[12] LevermoreCD , OliverM .ThecomplexGinzburgLan dauequationasamodelproblem[J].LecturesinApplMath, 1996, 31:141-189.
[13] WheelerAA , BottingerWJ , McFaddenGB .Phase fieldmodelforisothermalphasetransitionsinbinaryal loys[J].PhysRevA , 1992, 45 (10) :7424-7439.
[14] WheelerAA , McFaddenGB , BottingerWJ.Phasefieldmodelforsolidificationofaeutecticalloy[J].ProcRoySocLond, 1996, A452:495-525.
[15] BatesP , FifeP .SpectralcomparisonprinciplesfortheCahnHilliardandphasefieldequations, andtimescalesforcoarsening[J].PhysicaD , 1990, 43:335-348.
[16] WarrenJA , BoettingerWJ , CaginalpG .Aphase fieldmodelofdendriticgrowthinabinaryalloy[J].ActaMetallMater, 1995, A43:689-698.
[17] SteinbachI, PezzollaF , NestlerB .A phase fieldcon ceptformultiphasesystem[J].PhysicaD , 1996, 94:135-147.
[18] ContiM .Thermalandchemicaldiffusionintherapidsolidificationofbinaryalloys[J].PhysRevE , 2000, 61 (1) :642-650.
[19] ContiM .Interfacedynamics, instabilitiesandsolutebandsinrapiddirectionalsolidification[J].PhysRevE , 1998, 58 (2) :2071-2084.
[20] LosertW , StillmanDA , CumminsHZ , etal.Selec tionofdoubletcellularpatternsindirectionalsolidifica tionthroughspatiallyperiodicperturbations[J].PhysRevE , 1998, 58 (6) :7492-7506.
[21] AhmadNA , WheelerAA , BoettingerWJ, etal.So lutetrappingandsolutedraginaphase fieldmodelofrapidsolidification[J].PhysRevE , 1998, 58:3436-3450.
[22] KimSG , KimWT , SuzukiT .Phase fieldmodelforbinaryalloys[J].PhysRevE , 1999, 60:7186-7197.
[23] ContiM .Planarisothermalsolidificationfromanunder cooledmelt:unsteadysolutesegregationstudiedwiththephase fieldmodel[J].PhysRevE , 1997, 55:701-707.
[24] CaginalpG .Ananalysisofaphasefieldmodelofafreeboundary[J].ArchRationalMechAnal, 1986, 92:205-245.
[25] CaginalpG .StefanandHeleShawtypemodelsasasymptoticlimitsofphasefieldequations[J].PhysRevA , 1989, 39:887-896.
[26] CaginalpG , ChenX .Convergenceofthephasefieldmodeltoitssharpinterfacelimits[J].EurJApplMath, 1998, 9:417-445.
[27] CaginalpG , JonesJ.Aderivationandanalysisofphasefieldmodelsofthermalalloys[J].AnnPhys, 1995, 237:66-107.
[28] CaginalpG , FifePC .Phase fieldmethodsforinterfa cialboundaries[J].PhysRevB , 1986, 33 (11) :7792-7795.
[29] WangSL , SekerkaRF , WheelerAA , etal.Thermo dynamically consistentphase fieldmodelsforsolidifica tion[J].PhysicaD , 1993, 69:189-200.
[30] FriedE , GutinME .Continuumtheoryofthermallyin ducedphasetransitionsbasedonanorder parameter[J].PhysicaD , 1993, 68:326-343.
[31] ContiM .Solutetrappingindirectionalsolidificationathighspeed:a1Dstudywithaphase fieldmodel[J].PhysRevE , 1997, 56:3717-3720.
[32] CharachC , KeizmanY .Solutetrappingeffectsinpla narisothermalsolidificationofdilutebinaryalloys[J].PhysRevE , 1997, 54:588-599.
[33] WeelerAA , BottingerWJ, Mcfadden.Phase fieldmodelofsolutetrappingduringsolidification[J].PhysRevE , 1993, 47 (3) :1893-1909.
[34] KaplanT , AzizMJ, GrayLJ .ApplicationofOnsager’sreciprocityrelationstointerfacemotiondur ingphasetransformations[J].JChemPhys, 1989, 90:1133-1140.
[35] AzizMJ , BottingerWJ.Onthetransitionfromshort rangediffusion limitedtocollision limitedgrowthinal loysolidification[J].ActaMetallMater, 1994, 42 (2) :527-537.
[36] CollinsJB , LevineH .Diffuseinterfacemodelofdiffu sion limitedcrystalgrowth[J].PhysRevB , 1985, 31:6119-6122.
[37] BechhoeferJ, LowenH , TuckermanLS .Dynamicalmechanismfortheformationofmetastablephases[J].PhysRevLett, 1991, 67:1266-1269.
[38] RyoKobayshi.Modelingandnumericalsimulationofdendriticcrystalgrowth[J].PhysicaD , 1993, 63:410-423.
[39] WangSL , SekerkaRF .Computationofthedendriticoperatingstateatlargesupercoolingsbythephasefieldmodel[J].PhysRevE , 1996, 53 (4) :3760-3776.
[40] KarmaA , RappelWJ .Phase fieldmethodforcompu tationallyefficientmodelingofsolidificationwitharbi traryinterfacekinetics[J].PhysRevE , 1996, 53:3017.
[41] WangSL , SekerkaRF .Algorithmsforphasefieldcomputationofthedentriticoperatingstateatlargesu percoolings[J].JComputPhys, 1996, 127:110-117.
[42] WheelerAA , MurrayBT , SchaeferRJ .Computationofdendritesusingaphasefieldmodel[J].PhysicaD , 1993, 66:243-262.
[43] KarmaA .Phase fieldmodelofeutecticgrowth[J].PhysRevE , 1994, 49 (3) :2245-2250.
[44] ElderKR , DroletF , KosterlitzJM , etal.Directionalsolidificationintwoandthreedimensions[J].PhysRevLett, 1994, 72:677.
[45] ElderKR , GuntonJD , GrantM .Non isothermaleutecticcrystallization[J].PhysRevE , 1996, 54:6476-6484.
[46] DroletF , ElderKR , GrantM , etal.Phase fieldmod elingofeutecticgrowth[J].PhysRevE , 2000, 61 (6) :6705-6720.
[47] McffaddenGB , WheelerAA , BraunRJ, etal.Phase fieldmodelsforanistropicinterface[J].PhysRevE , 1993, 48 (3) :2016-2024.
[48] BeckermannC , DiepersHJ, SteinbachI, etal.Model ingmeltconvectioninphase fieldsimulationsofsolidifi cation[J].JComputPhys, 1999, 154:468-496.
[49] DiepersHJ, BeckermannC , SteinbachI .Simulationofconvectionandripeninginabinaryalloymushusingthephase fieldmethod[J].ActaMater, 1999, 47:3663-3678.
[50] TongX , BeckermannC , KarmaA .Velocityandshapeselectionsofdendriticcrystalsinaforcedflow[J].PhysRevE , 2000, 61:R49-R52.
[51] VollerVR , BeckermannC .Approximatemodelsofmi crosegregationwithcoarsening[J].MetandMatTransA , 1999, 30A :3016-3019.
[52] ContiM .Solidificationofbinaryalloys:thermaleffectsstudiedwiththephase fieldmodel[J].PhysRevE , 1997, 55:765-771.
[53] BatesPW , FifePC .DynamicsofnucleationfortheCahnHilliardequation[J].SIAMJApplMath, 1993, 53:990-1008.
[54] LangerJS .Theoryofspinodaldecompositioninalloys[J].AnnalsofPhysics, 1971, 65:53-86.
[55] CahnJW .Phaseseparationbyspinodaldecompositioninisotropicsystems[J].JChemPhys, 1965, 42:93-99.
[56] KogaT , KawasakiK .Spinodaldecompositioninbinaryfluids:Effectsofhydrodynamicinteractions[J].PhysRevA , 1991, 44:R817.
[57] GrantCP .SpinodaldecompositionfortheCahnHilliardequation[J].CommunicationsinPartialDifferentialE quations, 1993, 18:453-490.
[58] EyreDJ.SystemsofCahnHilliardequations[J].SIAMJApplMath, 1993, 53:1686-1712.
[59] PegoRL .FrontmigrationinthenonlinearCahnHilliardequation[J].ProcRoySocLondon, 1989, A422:261-278.
[60] BrushLN , SekerkaRF .Anumericalstudyoftwo di mensionalgrowthformsinthepresenceofanisotropicgrowthkinetics[J].JCrystalGrowth, 1989, 96:419-441.
[61] AkaiwaN , MeironDI.Two dimensionallate stagecoarseningfornucleationandgrowthathigh areafrac tions[J].PhysRevE , 1996, 54 (1) :R13-R16.
[62] MeakinP .Formationoffractalclustersandnetworksbyirreversiblediffusion limitedaggregation[J].PhysRevLett, 1983, 51:1119-1124.
[63] AlsunaidiA , Lach hubM .Cluster clusteraggregationinbinarymixtures[J].PhysRevE , 2000, 61 (1) :550-556.
[64] BogoyavlenskiyVA , ChernovaNA .Diffusion limitedaggregation:Arevisedmean fieldapproach[J].PhysRevE , 2000, 61 (5) :5422-5428.
[65] GadeDM , JoyMP .Self organizedcriticalityindy namicswithoutbranch[J].PhysRevE , 1998, 57 (5) :5019-5022.