CN106055887A - 基于相场模型下控制Ag2Ga纳米针长径比的研究方法 - Google Patents

基于相场模型下控制Ag2Ga纳米针长径比的研究方法 Download PDF

Info

Publication number
CN106055887A
CN106055887A CN201610361529.5A CN201610361529A CN106055887A CN 106055887 A CN106055887 A CN 106055887A CN 201610361529 A CN201610361529 A CN 201610361529A CN 106055887 A CN106055887 A CN 106055887A
Authority
CN
China
Prior art keywords
nanoneedle
draw ratio
microgranule
phase
reaction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610361529.5A
Other languages
English (en)
Other versions
CN106055887B (zh
Inventor
张俐楠
郑伟
陈超
吴立群
王洪成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiaxing Huaji Environmental Protection Technology Co., Ltd
Original Assignee
Hangzhou Electronic Science and Technology University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Electronic Science and Technology University filed Critical Hangzhou Electronic Science and Technology University
Priority to CN201610361529.5A priority Critical patent/CN106055887B/zh
Publication of CN106055887A publication Critical patent/CN106055887A/zh
Application granted granted Critical
Publication of CN106055887B publication Critical patent/CN106055887B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了基于相场模型下控制Ag2Ga纳米针长径比的研究方法,按如下步骤进行:一、建立Ag2Ga纳米针的三维计算模型;二、改变实验反应条件以及模拟时间来揭示纳米针长径比变化规律;在室温下,Ag微粒和液态Ga微粒之间会发生扩散形成Ag2Ga微粒,此时对其一端施加一个力拔出,形成Ag2Ga纳米针;在此过程中,实验反应的条件以及模拟时间都是影响Ag2Ga纳米针长径比的关键因素;进行如下操作:操作一:探究不同微粒浓度对Ag2Ga纳米针长径比变化的影响;操作二:探究不同模拟时间对Ag2Ga纳米针长径比变化的影响;操作三:探究不同拔出速度对Ag2Ga纳米针长径比变化的影响。

Description

基于相场模型下控制Ag2Ga纳米针长径比的研究方法
技术领域
本发明属于纳米材料研究技术领域,具体涉及一种基于相场模型下控制Ag2Ga纳米针长径比的研究方法。
背景技术
纳米针作为一维的纳米材料,具有独特的机械性、热稳定性、电子传输和光子传输性、光学性质、光导性和场发射效应等特点,因此,具有巨大的应用前景。例如,纳米针在力和质量传感器上的应用以及在纳米刻蚀上的应用。再如,纳米针在医学上也具有巨大得应用前景,如纳米针在医疗美容方面、利用穿刺功能向细胞注射药物等方面发挥作用。近十年来,国内外报道了大量的纳米针制备方法的研究。但是这些报道采用的都是一些传统的纳米针制备方法,这些方法具有加工成形难以控制,实验材料浪费的缺点。另外,利用传统的加工方法制备出的纳米针长径比较小,纳米针形态难以控制,应用价值较低,无法批量生产。因此,开发品质更优、直径更细长径比更高、制备更简单的纳米针具有更大的产业价值和商业价值。
在近几十年中,由于计算材料科学领域在各项计算方法上的改进使其计算成效已经获得了普遍的认可。在各种建模方法中,phase-field model(相场模型)已成为一个强大的并且被普遍认可的有效计算模型来模拟和预测材料微观结构的演化过程。Phase-fieldmodel的建立是基于diffuse-interface模型,这一模型概念的提出是在40年前由Cahn和Hilliard提出的。“phase-field”一词来源于diffuse-interface模型中用来定义空间和时间秩序的参数。换句话说,phase-field model是使用一组守恒和/或非守恒的变量在连续的界面区域里来描述材料的微观结构变化及其空间分布状况。与此同时,phase-fieldmodel能够直接链接到数据库使它可以预测任意形态复杂的微观结构的演化过程。这里指出的微观结构的长度尺度通常是在微米到纳米范围内,其大小、形状和空间排列结构特性在决定材料的物理和力学性能上起着至关重要的角色。微/纳米材料的各种运动演化过程在本质上是非线性的,因此,如何用数值方法表达出这种变化过程往往是急需解决的难题。其中,传统方法sharp-interface model是将二种或者多种成分或结构域之间的边界作为一个单一无过渡的界面。这种方法需要显式跟踪界面的位置,在一维系统中,此方法是可行的,但是发展到两维甚至于三维空间的状态下就会变得不切实际,此方法产生的计算难题将无法克服。
发明内容
基于上述现状,本发明公开一种基于相场模型下控制Ag2Ga纳米针长径比的研究方法。
本发明基于相场模型下控制Ag2Ga纳米针长径比的研究将基于相场模型建立Ag2Ga纳米针的三维计算模型;利用微观粒子Ag与Ga之间的化学反应并结合仿真模型,提出一种Ag2Ga纳米针成形机理研究的新方法;通过改变实验反应条件以及模拟时间来揭示纳米针长径比变化规律,提出一种控制纳米针长径比的新思路。
为实现上述目的,本发明采取如下技术方案:基于相场模型下控制Ag2Ga纳米针长径比的研究方法,按如下步骤进行:
一、建立Ag2Ga纳米针的三维计算模型;
二、改变实验反应条件以及模拟时间来揭示纳米针长径比变化规律;
在室温下,Ag微粒和液态Ga微粒之间会发生扩散形成Ag2Ga微粒,此时对其一端施加一个力拔出,形成Ag2Ga纳米针;在此过程中,实验反应的条件以及模拟时间都是影响Ag2Ga纳米针长径比的关键因素;进行如下操作:
操作一:探究不同微粒浓度对Ag2Ga纳米针长径比变化的影响;
操作二:探究不同模拟时间对Ag2Ga纳米针长径比变化的影响;
操作三:探究不同拔出速度对Ag2Ga纳米针长径比变化的影响。
所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,步骤一:Ag和Ga微粒的三维计算模型中包含四种成分相,定义Ag微粒为c1(x,y,z,t),Ag2Ga纳米针为c2(x,y,z,t),Ga微粒为c3(x,y,z,t),c1为Ag微粒的浓度占总的微粒浓度的百分比,c2为Ag2Ga微粒的浓度占总的微粒浓度的百分比,c3为Ga微粒的浓度占总的微粒浓度的百分比;模型的多种动力学包括Ag微粒、Ga微粒和Ag2Ga相互之间的扩散过程,Ag微粒、Ga微粒之间的化学反应过程,纳米针的形成过程和成形过程;模型的多种能量学包括自由能、化学能、梯度能;
定义在Ag微粒和Ga微粒之间发生的反应化学方程为:
2Ag+Ga=Ag2Ga (1)
对于上述的三维计算模型,根据Cahn-Hilliard模型,自由能采用如下方程式:
G = ∫ v { F c h + F g r a d } d V - - - ( 2 )
对于上述方程(2),Fch表示化学反应的能量变化,即在Ag和Ga发生化学反应生成Ag2Ga的过程中的化学能量变化,用如下方程式定义:
F c h = f 0 { c 1 2 ( 1 - c 1 ) 2 + c 2 2 ( 1 - c 2 ) 2 + c 3 2 ( 1 - c 3 ) 2 + ( 1 - c 1 - c 2 - c 3 ) 2 ( c 1 + c 2 + c 3 ) 2 } - - - ( 3 )
其中,f0表示反应过程的化学能量系数;
对于上述方程(2),Fgrad表示系统梯度能的能量变化,使用如下方程:
Fgrad=h11(▽c1)2+h22(▽c2)2+h33(▽c3)2+h12(▽c1c2)2+h13(▽c1c3)2+h23(▽c2c3)2 (4)
上式中,▽表示梯度计算,h11表示相c1内部梯度能量变化系数,h22表示相c2内部梯度能量变化系数,h33表示相c3内部梯度能量变化系数,h12表示相c1和相c2之间梯度能量变化系数,h13表示相c1和相c3之间梯度能量变化系数,h23表示相c2和相c3之间梯度能量变化系数;
对于上述模型,c1、c2、c3的密度场变量通过Cahn-Hilliard非线性扩散方程决定,并结合量守恒关系,得到如下控制方程:
∂ c 1 ∂ t = - ▿ ( - M 1 ▿ μ 1 + M 1 β ▿ φ ) - v 1 c 1 c 3 - - - ( 5 )
∂ c 3 ∂ t = - ▿ ( - M 3 ▿ μ 3 ) - v 2 c 1 c 3 - - - ( 7 )
其中,v1是c1的反应速率,v2是c3的反应速率,β是灵敏度系数;θ是关于Ag生成Ag2Ga的速率的系数;是关于Ga生成Ag2Ga的速率的系数;μ1、μ2、μ3分别表示Ag、Ag2Ga、Ga的化学势,分别用如下方程定义:
u 1 = ∂ G / ∂ c 1 - - - ( 8 )
u 2 = ∂ G / ∂ c 2 - - - ( 9 )
u 3 = ∂ G / ∂ c 3 - - - ( 10 )
对于上述方程(5)、(6)、(7),M1、M2、M3分别表示是微观粒子Ag、Ag2Ga、Ga的扩散率,分别用如下三个方程定义:
M 1 = M 0 [ { ∫ c 1 2 ( 1 - c 1 ) 2 dc 1 / ∫ 0 1 c 1 2 ( 1 - c 1 2 ) dc 1 } ( 1 - c 1 ) ] - - - ( 11 )
M 2 = M 0 [ { ∫ c 2 2 ( 1 - c 2 ) 2 dc 2 / ∫ 0 1 c 2 2 ( 1 - c 2 2 ) dc 2 } ( 1 - c 2 ) ] - - - ( 12 )
M 3 = M 0 [ { ∫ c 3 2 ( 1 - c 3 ) 2 dc 3 / ∫ 0 1 c 3 2 ( 1 - c 3 2 ) dc 3 } ( 1 - c 3 ) ] - - - ( 13 )
M0为材料的扩散系数。
所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,步骤二:操作一:通过探究不同的微粒浓度对纳米针长径比变化的影响,揭示纳米针长径比变化的规律,从而控制纳米针长径比变化;Ga的熔点为29.8摄氏度,因此实验温度控制在常温即可;Ag微粒与Ga微粒的反应用如下方程式表示
2Ag+Ga=Ag2Ga
取四组不同浓度比Ag:Ga分别为1.8:1、2:1、2.2:1、2.4:1,其它反应条件相同,其它反应条件包括温度、模拟时间、驱动力大小。
所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,步骤二:操作二:在常温下,取四组反应来模拟,模拟时间分别为t1、t2、t3、t4,t1<t2<t3<t4;在四组模拟时间下,分别对Ag2Ga纳米针一端施加相同大小的驱动力,并且将其缓慢的拔出;模拟时间t1时,时间较短,微粒之间的扩散和反应未充分完成,生成纳米针长径比小,直径较粗;模拟时间t2时,微粒之间的扩散和反应都较为完全,生成纳米针长径比较大,直径较细;模拟时间t3时,微粒之间的扩散和反应完全,生成纳米针长径比大,直径细;模拟时间t4时,微粒之间的扩散和反应过长,一些微粒吸附在纳米针表面,因此生成的纳米针较长,但直径较粗。
所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,步骤二:操作三:在常温下,通过不同的拔出速度V1、V2和V3,V1<V2<V3来探究,揭示拔出速度对纳米针的长径比变化影响;当拔出速度为V1且保持匀速,生成的纳米针的长径比较大,且纳米针的直径较细;当拔出速度为V2且保持匀速,生成的纳米针的长径比变小,纳米针直径变大;当拔出速度为V3且保持匀速,由于速度较大,纳米针一端很快脱离,生成的纳米针长度较短,直径较大,长径比小。
本发明基于相场模型下控制Ag2Ga纳米针长径比的研究通过基于相场模型下Ag2Ga纳米针三维计算模型,定义微观粒子的密度场变量,建立微观粒子的扩散方程,建立系统的自由能方程,建立Ag2Ga纳米针计算模型的控制方程。
与现有技术相比,本发明基于相场模型下控制Ag2Ga纳米针长径比的研究方法具有以下技术特点:
1、基于连续的相场模型,不需要显示跟踪界面的位置,大大降低了计算的复杂性。
2、建立基于相场模型的三维计算模型,提出了一种Ag2Ga纳米针成形制造过程研究的新方法。
3、利用仿真模型,改变实验反应条件以及模拟时间来揭示纳米针长径比变化规律,提出一种控制纳米针长径比的新思路。
附图说明
图1为Ag和Ga微粒的三维计算模型简单示意图。
图2为Ag和Ga形成Ag2Ga纳米针的简单示意图。
图3为不同微粒浓度纳米针长径比形态的简单示意图。
图4为不同模拟时间纳米针长径比形态的简单示意图。
具体实施方式
下面结合附图对本发明作进一步详细的说明。这些图均为简化的示意图,仅以示意方式说明本发明的基本结构,因此仅显示对本发明有关的构成。
本实施例基于相场模型下控制Ag2Ga纳米针长径比的研究方法,其按如下步骤进行:
建立Ag2Ga纳米针的三维计算模型
如图1所示,Ag和Ga微粒的三维计算模型中包含了多种能量学和多力学。
Ag和Ga微粒的三维计算模型中包含四种成分相,如图2所示,定义Ag微粒为c1(x,y,z,t),Ag2Ga纳米针为c2(x,y,z,t),Ga微粒为c3(x,y,z,t),空气相在此系统不予以考虑,c1为Ag微粒的浓度占系统总的微粒浓度的百分比,c2为Ag2Ga微粒的浓度占系统总的微粒浓度的百分比,c3为Ga微粒的浓度占系统总的微粒浓度的百分比。如图1所示的三维计算模型,此系统具有多种动力学和多种能量学特征;模型的多种动力学包括Ag微粒、Ga微粒和Ag2Ga相互之间的扩散过程,Ag微粒、Ga微粒之间的化学反应过程,纳米针的形成过程和成形过程;模型的多种能量学包括系统的自由能、化学能、梯度能等。
定义在Ag微粒和Ga微粒之间发生的反应化学方程为:
2Ag+Ga=Ag2Ga (1)
对于上述的三维计算模型,根据Cahn-Hilliard模型,系统的自由能可使用如下方程式:
G = ∫ v { F c h + F g r a d } d V - - - ( 2 )
对于上述方程(2),v为系统成分相的体积、dV求积分;Fch表示的是化学反应的能量变化,即在Ag和Ga发生化学反应生成Ag2Ga的过程中的化学能量变化,可使用如下方程式定义:
F c h = f 0 { c 1 2 ( 1 - c 1 ) 2 + c 2 2 ( 1 - c 2 ) 2 + c 3 2 ( 1 - c 3 ) 2 + ( 1 - c 1 - c 2 - c 3 ) 2 ( c 1 + c 2 + c 3 ) 2 } - - - ( 3 )
其中,f0表示反应过程的化学能量系数。
对于上述方程(2),Fgrad表示的是系统梯度能的能量变化,可使用如下方程:
Fgrad=h11(▽c1)2+h22(▽c2)2+h33(▽c3)2+h12(▽c1c2)2+h13(▽c1c3)2+h23(▽c2c3)2 (4)
▽表示的是梯度计算,h11表示的是相c1内部梯度能量变化系数,h22表示的是相c2内部梯度能量变化系数,h33表示的是相c3内部梯度能量变化系数,h12表示的是相c1和相c2之间梯度能量变化系数,h13表示的是相c1和相c3之间梯度能量变化系数,h23表示的是相c2和相c3之间梯度能量变化系数。
对于上述模型,c1、c2、c3的密度场变量通过Cahn-Hilliard非线性扩散方程决定,并结合量守恒关系,得到如下控制方程:
∂ c 1 ∂ t = - ▿ ( - M 1 ▿ μ 1 + M 1 β ▿ φ ) - v 1 c 1 c 3 - - - ( 5 )
∂ c 3 ∂ t = - ▿ ( - M 3 ▿ μ 3 ) - v 2 c 1 c 3 - - - ( 7 )
其中,v1是c1的反应速率,v2是c3的反应速率,β是灵敏度系数;θ是关于Ag生成Ag2Ga的速率的系数;是关于Ga生成Ag2Ga的速率的系数;对t求偏微分;μ1、μ2、μ3分别表示的是Ag、Ag2Ga、Ga的化学势,可分别用如下方程定义:
u 1 = ∂ G / ∂ c 1 - - - ( 8 )
u 2 = ∂ G / ∂ c 2 - - - ( 9 )
u 3 = ∂ G / ∂ c 3 - - - ( 10 )
对于上述方程(5)、(6)、(7),M1,M2,M3分别表示是微观粒子Ag,Ag2Ga,Ga的扩散率,可分别用如下三个方程定义:
M 1 = M 0 [ { ∫ c 1 2 ( 1 - c 1 ) 2 dc 1 / ∫ 0 1 c 1 2 ( 1 - c 1 2 ) dc 1 } ( 1 - c 1 ) ] - - - ( 11 )
M 2 = M 0 [ { ∫ c 2 2 ( 1 - c 2 ) 2 dc 2 / ∫ 0 1 c 2 2 ( 1 - c 2 2 ) dc 2 } ( 1 - c 2 ) ] - - - ( 12 )
M 3 = M 0 [ { ∫ c 3 2 ( 1 - c 3 ) 2 dc 3 / ∫ 0 1 c 3 2 ( 1 - c 3 2 ) dc 3 } ( 1 - c 3 ) ] - - - ( 13 )
M0为材料的扩散系数。
改变实验反应条件以及模拟时间来揭示纳米针长径比变化规律,提出一种控制纳米针长径比的新思路。
如图3所示,在室温下,Ag微粒和液态Ga微粒之间会发生扩散形成Ag2Ga微粒,这时候对其一端施加一个力拔出,形成Ag2Ga纳米针。在此过程中,实验反应的条件以及模拟时间都是影响Ag2Ga纳米针长径比的关键因素。
操作一:探究不同微粒浓度对Ag2Ga纳米针长径比变化的影响
根据控制方程(5)、(6)、(7),Ag2Ga纳米针的形成与Ag微粒的反应速率v1和Ga微粒的反应速率v2有着直接的影响,化学反应的速率又与反应物微粒的浓度以及反应的条件有关。因此,通过探究不同的微粒浓度对纳米针长径比变化的影响,揭示纳米针长径比变化的规律,从而控制纳米针长径比变化。Ga的熔点为29.8摄氏度,因此实验温度控制在常温即可。Ag微粒与Ga微粒的反应可用如下方程式表示
2Ag+Ga=Ag2Ga
取四组不同浓度比(Ag:Ga),分别为1.8:1、2:1、2.2:1、2.4:1,其它反应条件相同(包括温度,模拟时间,驱动力大小等)。整个操作过程可用图2简单表示。
操作二:探究不同模拟时间对Ag2Ga纳米针长径比变化的影响
根据控制方程(5)、(6)、(7),模拟时间长短与化学反应是否反应完全有着直接的联系,进而影响到生成Ag2Ga纳米针的形态。在常温下,取四组反应来模拟,模拟时间分别为t1、t2、t3、t4,t1<t2<t3<t4。在四组模拟时间下,分别对Ag2Ga纳米针一端施加相同大小的驱动力,并且将其缓慢的拔出。整个操作过程如图2简单示出。模拟时间t1时,由于时间较短,微粒之间的扩散和反应未充分完成,生成纳米针长径比小,直径较粗,如图3a所示;模拟时间t2时,微粒之间的扩散和反应都较为完全,生成纳米针长径比较大,直径较细,如图3b所示;模拟时间t3时,微粒之间的扩散和反应完全,生成纳米针长径比大,直径细,如图3c所示;模拟时间t4时,微粒之间的扩散和反应过长,一些微粒吸附在纳米针表面,因此生成的纳米针较长,但直径较粗,如图3d所示。
操作三:探究不同拔出速度对Ag2Ga纳米针长径比变化的影响
如图4所示。在Ag微粒和Ga微粒发生反应形成Ag2Ga纳米针后,需要对其一端施加一个驱动力将其拔出,拔出的速度对纳米针的成形也有着直接的关系。在常温下,其它反应条件相同,通过不同的拔出速度V1、V2和V3,V1<V2<V3来探究,揭示拔出速度对纳米针的长径比变化影响。如图4所示。当拔出速度为V1且保持匀速,生成的纳米针的长径比较大,且纳米针的直径较细,如图4a所示;当拔出速度为V2且保持匀速,生成的纳米针的长径比变小,纳米针直径变大,如图4b所示;当拔出速度为V3且保持匀速,由于速度较大,纳米针一端很快脱离,生成的纳米针长度较短,直径较大,长径比小,如图4c所示。
以上对本发明的优选实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。

Claims (5)

1.基于相场模型下控制Ag2Ga纳米针长径比的研究方法,其特征是按如下步骤进行:
一、建立Ag2Ga纳米针的三维计算模型;
二、改变实验反应条件以及模拟时间来揭示纳米针长径比变化规律;
在室温下,Ag微粒和液态Ga微粒之间会发生扩散形成Ag2Ga微粒,此时对其一端施加一个力拔出,形成Ag2Ga纳米针;在此过程中,实验反应的条件以及模拟时间都是影响Ag2Ga纳米针长径比的关键因素;进行如下操作:
操作一:探究不同微粒浓度对Ag2Ga纳米针长径比变化的影响;
操作二:探究不同模拟时间对Ag2Ga纳米针长径比变化的影响;
操作三:探究不同拔出速度对Ag2Ga纳米针长径比变化的影响。
2.如权利要求1所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,其特征是:
步骤一:Ag和Ga微粒的三维计算模型中包含四种成分相,定义Ag微粒为c1(x,y,z,t),Ag2Ga纳米针为c2(x,y,z,t),Ga微粒为c3(x,y,z,t),c1为Ag微粒的浓度占总的微粒浓度的百分比,c2为Ag2Ga微粒的浓度占总的微粒浓度的百分比,c3为Ga微粒的浓度占总的微粒浓度的百分比;模型的多种动力学包括Ag微粒、Ga微粒和Ag2Ga相互之间的扩散过程,Ag微粒、Ga微粒之间的化学反应过程,纳米针的形成过程和成形过程;模型的多种能量学包括自由能、化学能、梯度能;
定义在Ag微粒和Ga微粒之间发生的反应化学方程为:
2Ag+Ga=Ag2Ga (1)
对于上述的三维计算模型,根据Cahn-Hilliard模型,自由能采用如下方程式:
G = ∫ v { F c h + F g r a d } d V - - - ( 2 )
对于上述方程(2),Fch表示化学反应的能量变化,即在Ag和Ga发生化学反应生成Ag2Ga的过程中的化学能量变化,用如下方程式定义:
F c h = f 0 { c 1 2 ( 1 - c 1 ) 2 + c 2 2 ( 1 - c 2 ) 2 + c 3 2 ( 1 - c 3 ) 2 + ( 1 - c 1 - c 2 - c 3 ) 2 ( c 1 + c 2 + c 3 ) 2 } - - - ( 3 )
其中,f0表示反应过程的化学能量系数;
对于上述方程(2),Fgrad表示系统梯度能的能量变化,使用如下方程:
F g r a d = h 11 ( ▿ c 1 ) 2 + h 22 ( ▿ c 2 ) 2 + h 33 ( ▿ c 3 ) 2 + h 12 ( ▿ c 1 c 2 ) 2 + h 13 ( ▿ c 1 c 3 ) 2 + h 23 ( ▿ c 2 c 3 ) 2 - - - ( 4 )
上式中,表示梯度计算,h11表示相c1内部梯度能量变化系数,h22表示相c2内部梯度能量变化系数,h33表示相c3内部梯度能量变化系数,h12表示相c1和相c2之间梯度能量变化系数,h13表示相c1和相c3之间梯度能量变化系数,h23表示相c2和相c3之间梯度能量变化系数;
对于上述模型,c1、c2、c3的密度场变量通过Cahn-Hilliard非线性扩散方程决定,并结合量守恒关系,得到如下控制方程:
∂ c 1 ∂ t = - ▿ ( - M 1 ▿ μ 1 + M 1 β ▿ φ ) - v 1 c 1 c 3 - - - ( 5 )
∂ c 3 ∂ t = - ▿ ( - M 3 ▿ μ 3 ) - v 2 c 1 c 3 - - - ( 7 )
其中,v1是c1的反应速率,v2是c3的反应速率,β是灵敏度系数;θ是关于Ag生成Ag2Ga的速率的系数;是关于Ga生成Ag2Ga的速率的系数;μ1、μ2、μ3分别表示Ag、Ag2Ga、Ga的化学势,分别用如下方程定义:
u 1 = ∂ G / ∂ c 1 - - - ( 8 )
u 2 = ∂ G / ∂ c 2 - - - ( 9 )
u 3 = ∂ G / ∂ c 3 - - - ( 10 )
对于上述方程(5)、(6)、(7),M1、M2、M3分别表示是微观粒子Ag、Ag2Ga、Ga的扩散率,分别用如下三个方程定义:
M 1 = M 0 [ { ∫ c 1 2 ( 1 - c 1 ) 2 dc 1 / ∫ 0 1 c 1 2 ( 1 - c 1 2 ) dc 1 } ( 1 - c 1 ) ] - - - ( 11 )
M 2 = M 0 [ { ∫ c 2 2 ( 1 - c 2 ) 2 dc 2 / ∫ 0 1 c 2 2 ( 1 - c 2 2 ) dc 2 } ( 1 - c 2 ) ] - - - ( 12 )
M 3 = M 0 [ { ∫ c 3 2 ( 1 - c 3 ) 2 dc 3 / ∫ 0 1 c 3 2 ( 1 - c 3 2 ) dc 3 } ( 1 - c 3 ) ] - - - ( 13 )
M0为材料的扩散系数。
3.如权利要求2所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,其特征是:
步骤二:操作一:通过探究不同的微粒浓度对纳米针长径比变化的影响,揭示纳米针长径比变化的规律,从而控制纳米针长径比变化;Ga的熔点为29.8摄氏度,因此实验温度控制在常温即可;Ag微粒与Ga微粒的反应用如下方程式表示
2Ag+Ga=Ag2Ga
取四组不同浓度比Ag:Ga分别为1.8:1、2:1、2.2:1、2.4:1,其它反应条件相同,其它反应条件包括温度、模拟时间、驱动力大小。
4.如权利要求2或3所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,其特征是:
步骤二:操作二:在常温下,取四组反应来模拟,模拟时间分别为t1、t2、t3、t4,t1<t2<t3<t4;在四组模拟时间下,分别对Ag2Ga纳米针一端施加相同大小的驱动力,并且将其缓慢的拔出;模拟时间t1时,时间较短,微粒之间的扩散和反应未充分完成,生成纳米针长径比小,直径较粗;模拟时间t2时,微粒之间的扩散和反应都较为完全,生成纳米针长径比较大,直径较细;模拟时间t3时,微粒之间的扩散和反应完全,生成纳米针长径比大,直径细;模拟时间t4时,微粒之间的扩散和反应过长,一些微粒吸附在纳米针表面,因此生成的纳米针较长,但直径较粗。
5.如权利要求2或3所述基于相场模型下控制Ag2Ga纳米针长径比的研究方法,其特征是:
步骤二:操作三:在常温下,通过不同的拔出速度V1、V2和V3,V1<V2<V3来探究,揭示拔出速度对纳米针的长径比变化影响;当拔出速度为V1且保持匀速,生成的纳米针的长径比较大,且纳米针的直径较细;当拔出速度为V2且保持匀速,生成的纳米针的长径比变小,纳米针直径变大;当拔出速度为V3且保持匀速,由于速度较大,纳米针一端很快脱离,生成的纳米针长度较短,直径较大,长径比小。
CN201610361529.5A 2016-05-26 2016-05-26 基于相场模型下控制Ag2Ga纳米针长径比的研究方法 Active CN106055887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610361529.5A CN106055887B (zh) 2016-05-26 2016-05-26 基于相场模型下控制Ag2Ga纳米针长径比的研究方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610361529.5A CN106055887B (zh) 2016-05-26 2016-05-26 基于相场模型下控制Ag2Ga纳米针长径比的研究方法

Publications (2)

Publication Number Publication Date
CN106055887A true CN106055887A (zh) 2016-10-26
CN106055887B CN106055887B (zh) 2019-05-14

Family

ID=57174751

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610361529.5A Active CN106055887B (zh) 2016-05-26 2016-05-26 基于相场模型下控制Ag2Ga纳米针长径比的研究方法

Country Status (1)

Country Link
CN (1) CN106055887B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106771373A (zh) * 2017-01-10 2017-05-31 杭州电子科技大学 基于相场模型利用激光控制硅基表面形态的研究方法
CN113387323A (zh) * 2021-05-24 2021-09-14 杭州电子科技大学 一种基于电场控制Ag2Ga纳米针阵列成形的制备方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090233086A1 (en) * 2008-03-17 2009-09-17 Fujifilm Corporation Metal oxide microparticles, transparent conductive film, and dispersion
CN102311095A (zh) * 2011-08-09 2012-01-11 吉林大学 一种在微流控芯片中制备多级金属微纳结构的方法
CN103043597A (zh) * 2011-10-13 2013-04-17 聊城大学 以ZnO晶体为模型制备金属微纳米阵列电极的方法
KR20160012727A (ko) * 2014-07-25 2016-02-03 경희대학교 산학협력단 3차원 플라즈모닉 나노 구조체 및 그 제조방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090233086A1 (en) * 2008-03-17 2009-09-17 Fujifilm Corporation Metal oxide microparticles, transparent conductive film, and dispersion
CN102311095A (zh) * 2011-08-09 2012-01-11 吉林大学 一种在微流控芯片中制备多级金属微纳结构的方法
CN103043597A (zh) * 2011-10-13 2013-04-17 聊城大学 以ZnO晶体为模型制备金属微纳米阵列电极的方法
KR20160012727A (ko) * 2014-07-25 2016-02-03 경희대학교 산학협력단 3차원 플라즈모닉 나노 구조체 및 그 제조방법

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
曹保胜: "《纳米尺度下体系自由能对扩散模型的影响》", 《材料导报:研究篇》 *
王杰: "《材料微结构演化的相场模拟》", 《固体力学学报》 *
龙永强: "《Cu-Ni-Si合金微结构及相变的计算与模拟》", 《中国博士学位论文全文数据库》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106771373A (zh) * 2017-01-10 2017-05-31 杭州电子科技大学 基于相场模型利用激光控制硅基表面形态的研究方法
CN106771373B (zh) * 2017-01-10 2019-03-08 杭州电子科技大学 基于相场模型利用激光控制硅基表面形态的研究方法
CN113387323A (zh) * 2021-05-24 2021-09-14 杭州电子科技大学 一种基于电场控制Ag2Ga纳米针阵列成形的制备方法
CN113387323B (zh) * 2021-05-24 2024-04-05 杭州电子科技大学 一种基于电场控制Ag2Ga纳米针阵列成形的制备方法

Also Published As

Publication number Publication date
CN106055887B (zh) 2019-05-14

Similar Documents

Publication Publication Date Title
Wu et al. Cluster expansion method and its application in computational materials science
Li et al. Non-contact manipulation of nonmagnetic materials by using a uniform magnetic field: Experiment and simulation
Modather et al. An analytical study of MHD heat and mass transfer oscillatory flow of a micropolar fluid over a vertical permeable plate in a porous medium
Sarkar et al. Bioconvection in non-Newtonian nanofluid near a perforated Riga plate induced by haphazard motion of nanoparticles and gyrotactic microorganisms in the attendance of thermal radiation and Arrhenius chemical reaction: sensitivity analysis
Popa et al. Tutorial on the elastic theory of spin crossover materials
Xia et al. The role of computer simulation in nanoporous metals—A review
Ahmad et al. Influential study of novel microorganism and nanoparticles during heat and mass transport in Homann flow of visco-elastic materials
CN106055887A (zh) 基于相场模型下控制Ag2Ga纳米针长径比的研究方法
Liu et al. Discrete Element Modelling of the Mechanical Behavior of Sand–Rubber Mixtures under True Triaxial Tests
Aljuaydi et al. Numerical investigations of ion slip and Hall effects on Cattaneo-Christov heat and mass fluxes in Darcy-Forchheimer flow of Casson fluid within a porous medium, utilizing non-Fourier double diffusion theories through artificial neural networks ANNs
Gaikwad et al. Atomistic investigation of process parameter variations on material deformation behavior in nanoimprint lithography of gold
Reddy et al. Numerical study of TC4-NiCr/EG+ Water hybrid nanofluid over a porous cylinder with Thompson and Troian slip boundary condition: Artificial neural network model
Shoaib et al. A design of soft computing intelligent networks for MHD Carreau nanofluid model with thermal radiation
Kojima et al. Search strategy for rare microstructure to optimize material properties of filled rubber using machine learning based simulation
CN101853334A (zh) 化学反应扩散中的免方程多尺度模拟方法
Galal et al. Intelligent neural computing to investigate the heat and mass transmission in nanofluidic system between two rotating permeable disks: Supervised learning mechanism
CN112895408A (zh) 基于相场模型的聚合物表面微结构可控成形机理研究方法
Park et al. Mesoscale Simulations of Polymer Solution Self-Assembly: Selection of Model Parameters within an Implicit Solvent Approximation
Banerjee et al. Phase-diffusion pattern in quantum-nondemolition systems
CN106021938A (zh) 一种Ag2Ga纳米针成型机理及尖端形貌控制研究方法
CN104091020A (zh) 聚合物结构和性能的计算模拟方法
CN113205862A (zh) 一种诱导a(bc)3杂臂星形共聚物结构相变的方法
Kaswan et al. Stability Analysis of Dual Solutions for Nonlinear Radiative Magnetohydrodynamic Flow of Ag− TiO2/H2O Hybrid Nanofluid over a Nonlinearly Shrinking Surface
Li et al. Study of magnetic microspheres screening in microfluidic chip based on co-simulation method
Chan A continuum study of ionic layer analysis for single species ion transport in coaxial carbon nanotubes

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201112

Address after: 310012 room 2603, building 8, No. 2, xiyuanba Road, Sandun Town, Xihu District, Hangzhou City, Zhejiang Province

Patentee after: HANGZHOU ZHUILIE TECHNOLOGY Co.,Ltd.

Address before: 310018, No. 1, No. 2, Poplar Street, Hangzhou economic and Technological Development Zone, Hangzhou, Zhejiang

Patentee before: HANGZHOU DIANZI University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201204

Address after: Room 1004-3, building 8, 3333 Guangyi Road, Daqiao Town, Nanhu District, Jiaxing City, Zhejiang Province

Patentee after: Jiaxing Huaji Environmental Protection Technology Co., Ltd

Address before: 310012 room 2603, building 8, No. 2, xiyuanba Road, Sandun Town, Xihu District, Hangzhou City, Zhejiang Province

Patentee before: HANGZHOU ZHUILIE TECHNOLOGY Co.,Ltd.