磁性高分子链的磁性质和构象性质的Monte Carlo模拟

采用退火 (Annealing)MonteCarlo方法 ,从高温到低温顺序模拟了简立方格点上考虑最近邻Ising相互作用的磁性高分子链在不同温度的磁性质和构象性质 .磁性高分子链在低温下存在自发磁矩 ,无限长链的临界温度Tc=1 77± 0 0 5J kB.在临界温度附近 ,高分子链经历了从伸展的无规线团到紧缩球体的塌缩相变 .对链的尺寸、形状、近邻数及能量的分析表明 ,高分子链的构象性质从温度Tc=1 77开始发生较明显的变化 ,这表明高分子Ising链的相变是Ising相互作用和链节运动协同作用的结果 .     

吸波材料与微波暗室问题的数学模型

为分析尖劈形材料的吸波性能,建立了半圆镜像模型将辐射波在尖劈空间的传播等效为射线对等分半圆的切割,确定了反射波线和吸波材料几何参数之间的解析关系.对于微波暗室问题,假设在余弦辐射体上随机选取一个方向的能量波束集中了其所有能量,分别建立了加权能量分布模型和Monte Carlo模型,仿真验证了两模型结论的一致性.     

磁性纳米颗粒系统偶极相互作用的Monte Carlo研究

采用局域Monte Carlo方法模拟不同易轴分布的简单立方排列单分散单畴Fe纳米颗粒系统的ZFC-FC曲线及磁滞回线.结果表明：随着偶极相互作用的增强，系统的阻塞温度T_B逐渐增大，且ZFC曲线的峰变宽.说明偶极相互作用使得系统的有效能垒提高，分布宽度增加.研究FC曲线磁化强度的倒数与温度关系，发现偶极相互作用系统中存在反铁磁有序.系统的阻塞态及超顺磁态的磁滞回线表明，极低低温下，随着偶极相互作用的增强，系统的矫顽力和剩磁减小，偶极相互作用阻碍系统的磁化；系统处于超顺磁态，各向异性作用及偶极相互作用使得系统的磁化曲线偏离Langevin曲线且偶极相互作用展现出退磁相互作用效应.偶极相互作用增强，系统磁化曲线与Langevin曲线偏差量的最大值向低场移动.在偶极相互作用下，易轴与外场夹角为45°的磁性纳米颗粒系统的平均有效能垒和有效能垒分布宽度较易轴随机分布系统的大.     

随机波动模型下算术亚式期权的Monte Carlo模拟定价

在随机波动模型下,研究亚式期权的定价问题.推导出了标的资产及其随机波动模型的路径,利用对偶变量法对亚式期权进行数值模拟计算,并对随机波动模型下与B-S模型下的欧式期权和亚式期权定价结果进行比较,最后给出了具有固定敲定价格和浮动敲定价格的算术亚式期权的数值计算结果.     

ABC三嵌段共聚物胶束从多间隔结构到多核结构转变的Monte Carlo模拟

采用Monte Carlo模拟方法研究了疏水-亲水-疏水(H-P-H)型ABC三嵌段共聚物在B嵌段的选择性溶剂中的自组装行为. 模拟结果表明, 通过调节A嵌段和C嵌段的疏水性和二者之间的不相容性, 体系中可以形成多种形貌各异的胶束. 根据胶束中疏水核结构的特点, 这些胶束大体上可以被分为多核型胶束和多间隔型胶束两种类型. 通过增强疏水嵌段的疏水性或降低A嵌段和C嵌段间的不相容性, H-P-H型ABC三嵌段共聚物胶束能够发生从多核型胶束向多间隔型胶束的转变. 进一步分析胶束中聚合物的链构象等微观结构信息发现, A嵌段和C嵌段间的排斥作用和疏水作用之间存在竞争关系, 而这种竞争关系是影响体系中形成多核型胶束还是多间隔型胶束的决定性因素.     

快中子照相模拟分析与实验验证

首次系统地推导出快中子照相像素值形成解析式, 建立图像反差不等式, 并利用该不等式首次对图像对比度与源强、照射时间和散射之间的关系进行说明. 并在像素值解析式基础上编制快中子照相模拟程序, 利用该程序对空间分辨率和图像对比度进行模拟, 并与实验对照, 研究结果表明空间分辨率模拟效果好于实验, 图像对比度模拟效果与实验相当. 最后通过对狭缝、方孔以及多材质组成的复杂样品模拟并与实验对照, 结果显示模拟效果与实验照片在反差灵敏度效果上非常一致, 该模拟计算方法可为实验设计和工程应用提供参考.     

扩散模型和凝聚模型耦合作用下胶体凝聚动力学的Monte Carlo模拟研究

以扩散模型(Ds(γ)=D0×sγ)和凝聚模型(Pij(σ)=P0×(i×j)σ)为基础,对胶体体系随时间的演变、团簇大小分布及其标度关系、团簇的重均大小S(t)的变化规律以及模型对最终分形维数的影响四个角度进行了比较研究,发现扩散指数γ0和凝聚概率指数σ0对胶体的凝聚动力学过程有相似的影响.本文在较宽的γ和σ取值范围内,对胶体的凝聚动力学进行了模拟研究,对慢速凝聚向快速凝聚的转化机理作了定量分析,并进一步分析了在团簇-团簇凝聚(CCA)模型下,得到类似扩散置限凝聚(DLA)模型的凝聚体的物理意义,结果表明:(1)γ0代表了体系中团簇或单粒做"定向运动"而非无规则的布朗运动的情况.这种"定向运动"的推动力可能来自于大团簇产生的强"长程范德华力"、"电场力"等,或来自于体系边界处的外力场的作用.(2)当σ0时,体系成为先快后慢的慢速凝聚,这可能对应大团簇为一排斥中心,即胶体颗粒存在"排斥力场"的现象.(3)证实了团簇的重均大小在凝聚过程的早期按指数规律增长,而后期按幂函数规律增长的实验现象.模拟研究还表明,胶体体系的凝聚动力学过程,在σ0时是一个存在正反馈机制的非线性动力学过程,而在σ0时则体现出负反馈的特征.     

运用Mathematica软件包模拟M/M/1队列

运用Mathematica软件包模拟一个M/M/1队列问题,显现该软件包对于数学研究与教学的重要价值.通过对问题的分析制定模拟方案,利用计算机仿真形象地展示排队过程.其方法适用于一般队列问题.     

随机介质背散射二维Mueller矩阵的数值模拟与分析

推导了随机介质背散射Mueller矩阵的直接计算公式,并运用矢量Monte Carlo方法进行了数值模拟.结果表明随机介质背散射二维Mueller矩阵方位关系随散射系数的减小而增强,而与微粒大小关系不大;Mueller矩阵元素绝对值的空间分布随径向呈近似指数规律衰减,矩阵元素的方位变化具有周期性.对称系统的二维Mueller矩阵的花样图中仅有7幅独立,其余9幅可通过对称、旋转变换得到.