计算材料学 2 精选文档.pptVIP

FICK 第二定律 2 1 2 2 2 x u u u t u u x u t u t x x t x t x x t x t x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? v x C C C D t C C x C D t t x x t x t x x t x t x 2 1 2 2 2 C ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? Fick 第二定律 —— 稳态扩散解 REAL C0(0:1000+1),C(0:1000+1) !C(DISTANCE) C0=0.1 C0(0)=0.8 !BOUNDARY CONDITION C0(1001)=0.1 !BOUNDARY CONDITION C=C0 OPEN(1,FILE=F:\DIF.dat) DO JT=1,50000 ！ Time DO IX=1,1000 ! distance C(IX)=C0(IX)+0.45*(C0(IX+1)+C0(IX-1)-2.0*C0(IX)) ! C(0)=0.8 ! C(1001)=0.1 IF(JT==50000) WRITE(1,*) IX,C(IX) END DO C0=C END DO END 应用之二 生日问题 MC 模拟 假设有 n 个人在一起，各自的生日为 365 天之一，根据概率理论，与很多人的直 觉相反，只需 23 个人便有大于 50 ％的几 率人群中至少有 2 个人生日相同。 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 365 ) 1 ( 365 ...... 365 363 365 364 365 365 1 n n 理论几率 模拟几率 10 0.117 0.110 20 0.411 0.412 23 0.527 0.520 30 0.706 0.692 45 0.941 0.936 46 50 0.986 INTEGER M(1:10000), NUMBER1(0:364), NUMBER2 REAL X,Y ISEED=RTC() DO J=1,10000 NUMBER1=0 X=RAN(ISEED) NUMBER1(0)=INT(365*X+1) JJJ=1 DO I=1,365 Y=RAN(ISEED) NUMBER2=INT(365*Y+1) ETR=COUNT(NUMBER1.EQ.NUMBER2) IF (ETR==1) THEN EXIT ELSE JJJ=JJJ+1 M(J)=JJJ NUMBER1(I)=NUMBER2 MC 在材料学领域的应用 —— 随机行走 背景 如，布朗运动－－－最简单、无限制随机 行走 (Unrestricted randon walk,RW ) start end ? ? t R ? 平均平方端 - 端位移 : ?? v t t R 2 2 ? 2 1 ? v , 自然科学和社会生活中很多现象都与随机运动有关 可以模拟的 内容？ 扩散； 分子运动； 。。。。。。 如图所示，第 i 个分子在经过 N 步随机行走 后距原点距离为 R ，对 n 个分子每步的位移 平方求和后取平均值就得到了所有分子距 原点的方均距离 R2 ： 计算材料学 Computational Materials Science 第二讲 蒙特卡罗（ Monte Carlo ）方法 主讲：张晖 电话： Email: huizhangahut.edu 本堂课 主要内容