论文部分内容阅读
生物大分子在微尺度下的流动特性，是近年来生物大分子研究的热点之一。DNA是典型的生物大分子。DNA的流动特性，即DNA在流场的不同条件下与时间有关的变形和流动规律，
从19世纪60年代开始，人们把经典力学的方法和分子生物的实验和统计分析技术相结合对DNA的杆状螺旋结构进行研究，取得了大量成果。近年来，随着计算机水平的大幅提高，布朗动力学的计算机模拟作为一种数值研究手段在流变学领域得到广泛的利用。对给定的模型可以很方便的设计与之相应的计算机模拟程序。
本文采用了带有粘性的FENE模型来实现对DNA分子动力学的模拟，在模拟时运用了布朗动力学方法，所得到的主要结论如下：
（1）导出了带有粘性的FENE模型的动力学方程和Euler迭代公式，运用Euler方法模拟了该模型在延展流中的动力学，将教值模拟结果与Smith的实验数据进行了比较，发现结果与实验数据比较吻合。
（2）给出了强1阶收敛的Milstein方法的迭代公式，模拟了带有粘性的珠-簧链模型在剪切流中的动力学问题，计算出不同Weissenberg数下的分子分数伸长，统计出分子在不同Weissenberg数下的伸长分布，将数值模拟结果与Smith的实验数据进行了比较分析，得到了较好的结论。
（3）由带有粘性的FENE模型的动力学方程，推出了其所受应力的微分方程，可以用来研究应力的性质。