具有表面弛豫的液体自扩散的Monte Carlo模拟

计算物理 ›› 2009, Vol. 26 ›› Issue (1): 42-48.

具有表面弛豫的液体自扩散的Monte Carlo模拟

蔡淑惠, 陈巧龙, 蔡聪波

厦门大学物理系和通信工程系, 福建厦门 361005

收稿日期:2007-07-24 修回日期:2008-01-23 出版日期:2009-01-25 发布日期:2009-01-25
作者简介:蔡淑惠(1968-),女,福建,教授,博士,从事磁共振相关研究．
基金资助:
国家自然科学基金(批准号10575085和10605019);福建省自然科学基金(项目编号2008J0028)资助项目

Monte Carlo Simulation of Liquid Self-diffusion with Surface Relaxation

CAI Shuhui, CHEN Qiaolong, CAI Congbo

Departments of Physics and Communications Engineering, Xiamen University, Xiamen 361005, China

Received:2007-07-24 Revised:2008-01-23 Online:2009-01-25 Published:2009-01-25

摘要/Abstract

摘要： 研究多孔介质中液体的自扩散行为能获得介质的微观结构信息,有助于了解介质中液体的传输性质.以具有不同孔隙大小的无限长圆柱体模型中的液体为对象,采用Monte Carlo随机行走方法,模拟存在表面弛豫时液体的自扩散系数和核自旋磁化强度随时间的变化,导出将NMR弛豫参数和随机行走参数联系在一起的表达式.结果表明:在快扩散区,液体的核自旋磁化强度随弛豫时间呈单指数衰减,且自扩散系数在短时情况下独立于表面弛豫率;在慢扩散区,液体的核自旋磁化强度衰减和自扩散系数在短时情况下均与表面弛豫率无关.模拟结果与理论分析相吻合,可用于求解介质的表面积与体积之比及孔径大小等结构信息.

关键词: 自扩散, 表面弛豫, 磁化强度衰减, Monte Carlo模拟

Abstract: Time-dependent self-diffusion coefficients and magnetization attenuation of liquid with surface relaxation in porous media in an infinitely long cylinder are calculated with Monte Carlo random walk method.Relationship between parameters of NMR relaxation and random walk is obtained.In fast-diffusion area, magnetization of nuclear spins attenuates exponentially,and self-diffusion coefficients are independent of surface relaxivity within short diffusion time.In slow-diffusion area,both magnetization attenuation and self-diffusion coefficients are independent of surface relaxivity within short diffusion time.These results agree well with theoretical analysis.Surfaceto-volume ratio and pore size can be determined.

Key words: self-diffusion, surface relaxation, magnetization attenuation, Monte Carlo simulation

中图分类号:

O411.3

蔡淑惠, 陈巧龙, 蔡聪波. 具有表面弛豫的液体自扩散的Monte Carlo模拟[J]. 计算物理, 2009, 26(1): 42-48.

CAI Shuhui, CHEN Qiaolong, CAI Congbo. Monte Carlo Simulation of Liquid Self-diffusion with Surface Relaxation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 26(1): 42-48.

[1]	龙光波, 郑永刚, 欧建文, 王兴华, 黄霞, 王艳芳. 逆Compton散射的MonteCarlo模拟[J]. 计算物理, 2016, 33(2): 136-146.
[2]	周璇, 张志东. 基于修正的Gruhn-Hess两体势模型研究弯曲形变向列相液晶盒[J]. 计算物理, 2010, 27(4): 598-602.
[3]	许晓军, 王凤飞. 不同温度下团簇生长的Monte Carlo模拟[J]. 计算物理, 2009, 26(5): 758-762.
[4]	寻之朋, 唐刚, 韩奎, 夏辉, 郝大鹏, 周伟, 杨细全. 1+1维Wolf-Villain模型奇异标度行为的数值模拟[J]. 计算物理, 2009, 26(2): 287-292.
[5]	刘艳辉, 胡林. DNA微环拓扑性质的Monte Carlo模拟[J]. 计算物理, 2009, 26(1): 152-158.
[6]	张连珠, 于富强, 孔令宜, 何琪. 等离子体金属表面氮化中放电温度对氮原子态粒子的影响[J]. 计算物理, 2004, 21(2): 173-178.
[7]	邓辉球, 胡望宇, 舒小林. 二元合金表面偏析的Monte Carlo方法模拟[J]. 计算物理, 2003, 20(2): 137-141.
[8]	李之杰, 王泽辉, 川村孝一, 山村泰道. 采用壳层效应屏蔽长度Monte Carlo方法计算溅射产额[J]. 计算物理, 2002, 19(4): 290-292.
[9]	贺朝会, 陈晓华, 李国政. 高能质子单粒子翻转效应的模拟计算[J]. 计算物理, 2002, 19(4): 367-371.
[10]	王瑜, 周洪楠. 海尔-波普彗星旋喷流结构的数值模拟[J]. 计算物理, 2001, 18(1): 72-74.
[11]	王治虎, 罗孟波. 一端附壁高分子链的Monte Carlo研究[J]. 计算物理, 2000, 17(6): 645-648.
[12]	孔繁敏, 冯存峰, 张学尧, 傅宇, 张乃健, 何瑁, 王承瑞, 谭有恒. 羊八井ARGO实验中原初γ和质子的区分[J]. 计算物理, 2000, 17(3): 280-285.
[13]	冯存峰, 孔繁敏, 张学尧, 何瑁, 戴志强, 张乃健. 用神经网络分辨“膝”区原初质子成分[J]. 计算物理, 1999, 16(6): 568-572.
[14]	查敏, 朱清棋, 丁林恺. 对宇宙线广延大气簇射Čerenkov光的计算机模拟[J]. 计算物理, 1999, 16(6): 573-579.
[15]	张春, 吴克勇, 胡翠英. 电子注入在介质中输运的数值模拟[J]. 计算物理, 1999, 16(2): 183-191.