间断Galerkin有限元方法应用及生物分子模拟中的方法研究

来源 :中国科学院数学与系统科学研究院 | 被引量 : 0次 | 上传用户:B08050402
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文分为两部分:第一部分主要研究了间断Galerkin有限元方法(DG)在特征值问题中的应用;第_部分主要研究了分子力学模拟中如何构造模型才能够更快速高效的模拟分子系统,如蛋白质分子的模拟,提出三层模型,并导出了半解析镜像电荷方法和最小二乘镜像电荷方法。   DG方法由于其易于并行、精度高等特点被广泛应用。本文将一类DG方法应用到椭圆型特征值问题中,并对非光滑区域,如区域有裂缝、区域含有凹角等问题,计算了Laplace特征值问题的特征值。通过比较DG方法和协调有限元、非协调有限元及其他有限元校正方法发现,DG方法在处理区域有裂缝或者某些凹角等问题的特征值问题中有更高的精度。同时,针对这些不规则的区域,本文给出了在这些区域上构造DG方法的详细过程,通过修改可以应用到其他类似的特征值问题中。   静电场在生物分子的计算过程中发挥着重要作用。静电相互作用是长程作用力,强烈地依赖丁待研究的生物分子周崮的溶液环境。因此,模拟该生物分子系统的计算量非常大。对连续模型,即考虑靠近生物分子的所有溶液中的分子或离子,而对于较远的溶液视为高介电常数的连续介质,本文针对一个点电荷的模型问题,提出三层分子模型,并且可以将其应用到大型生物分子计算中。对一个点电荷处于球形真空中,外而被水溶液或者离子溶液包围,如果介电常数在界面上不连续,反应电场会在球形区域内部靠近界面的时候出现问题。因此,本文在真空和溶液交界面构造了一个过渡层,其介电常数跟该处到球心的距离相关,并且使得介电常数在界面上也连续。这就克服了以前的模型中由于介电常数的不连续性导致的没有物理意义的反应场。   针对该三层模型,本文分别构造了半解析镜像电荷方法和最小二乘镜像电荷方法对其进行计算。通过跟序列形式的电势解析解比较发现,甚至只用一个镜像电荷,两种方法都能够很好的近似电势的精确解。特别地,通过选取多个镜像电荷,最小二乘镜像电荷方法可以得到非常高的精度。而且,在对离子溶液的计算中,最小二乘镜像方法不依赖于溶液强度。对于外层是水溶液的情形,针对不同位置的点电荷,本文给出了其相应的第一个和第二个镜像电荷,在以后的应用中可以直接通过表中的数据插值得到任意位置点上电荷的镜像电荷。
其他文献
最小支撑树问题是一类经典的组合优化问题,它已经有很好的解决方法.本论文重点研究在无向图上支撑树扩容问题,这是对最小支撑树问题的推广。   给定连通图G=(V,E;l,b,c)及正常数
由于许多现实世界的物理系统可以由分数阶状态方程更好地描述,分数阶微积分越来越被广泛地应用于信号处理与控制、图像处理、流体力学、分形理论、电力分形网络、分数阶正弦振
作为一种处理模糊和不确定性问题的数据分析工具,粗糙集理论是非常具有优势的。它能有效的处理不完备的、不精确的、不确定的那些问题,现在已经和数据挖掘、模式识别成功的结合
本文从向量丛E及其上的Ehresmann联络出发,在一阶射流丛()1E上构造了不依赖于度量的协变经典场论。通过对规范不变性的分析讨论了从主丛出发构造规范场论的原因。在此基础上介绍
本文主要研究欧式外汇期权定价模型。由于两个外汇欧式期权的定价和外汇的波动率微笑以及其间的相关系数紧密相关,仅利用Black-Scholes框架无法进行正确定价。因此本文以Black
近年来,鲁棒主成分分析模型(Robust Principal Component Analysis,RPCA)在图像去噪、图像修复、图像分类以及运动目标分割等领域取得了成功运用。其主旨思想认为一系列含噪音
一阶拟线性双曲组的节点状态的精确边界能控性在等温气体传输等问题中有广泛的应用([26]).本文基于具有非齐次项的一阶拟线性双曲组的半整体经典解理论,利用统一的构造性方法
本文认为,诚信、团结合作教育可以很好地解决我国教育工作中广泛存在的高分低能问题,所以把诚信、团结合作教育上升到一定高度,并针对农村小学生的特点,探索了一些可行的方法
计算机代数方法和工具,如Grobner基计算、QEPCAD、DISCOVERER等已经应用到不变量生成,寻找秩函数,可达集计算等程序验证领域。本文结合用DISCOVERER寻找秩函数及浮点误差分析方
本论文主要分为两大部分。   第一部分主要研究半空间Zd×Z+上的Bernoulli边渗流模型生成的无穷开簇上的接触过程.得到的结论是:当p大于Zd×Z+上的Bernoulli边渗流模型的临