移动有限元方法在流体中的应用

来源 :北京大学 | 被引量 : 0次 | 上传用户:kongguoying
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在过去的二十年里,自适应方法的研究发展非常迅速。在相当多的实际问题中,由于解的奇异性质的存在,传统的一致网格计算会极大地浪费计算机资源,并使得数值求解无法实现。相反地,通过自适应的方法来分布网格,可以在很大程度上提高计算的性能。[R.Li,T.Tang,andP.W.Zhang,J.Comput.Phys.,170(2001),pp.562-588]中提出了一种基于调和映射的移动有限元方法,将方程求解和网格移动完全分开,从而使得移动网格方法在各类不同问题中的应用归结为控制函数的构造问题。在它的基础上,本文首次尝试针对不可压缩流体的原始变量Navier-Stokes方程,设计移动网格方法。对于不可压缩流体,设计的主要难点在于如何在网格移动的过程中仍然保持不可压缩条件,所以设计一个保持对速度场散度为零条件的迭代算法是本算法的主要特点之一。另外,不可压缩流体解的奇异性不像间断问题那么强,流场在几乎全局具有非平凡解。这就使得一方面要求精确求解带有奇性的区域,另一方面在其他区域也要保持足够的精度。我们通过对速度场的后验误差估计来构造控制函数。用这种移动网格方法对标准的周期双剪切层问题进行数值模拟,在相当稀疏的网格上就不再有虚假涡出现。 我们把上述针对不可压缩流体的移动网格方法和水平集方法结合起来模拟带表面张力的不可压缩二相流。根据二相流控制方程的自身特点,对不同网格之间的插值格式作了相应的修改。并设计了一个简单易行的控制函数,仅仅考虑了水平集函数和分界面曲率的影响。我们分别对单个气泡的上浮、两气泡的融合、液滴的下落等物理问题进行了模拟,数值结果充分说明了该方法的有效性。 本文还将此移动网格方法推广到球面问题上,但困难在于球面的曲率为正,使得球面到球面的调和映射不再具有唯一性。我们通过引进摄动(perturbed)能量泛函,得到了一个类调和映射,解决了这个问题。对于这个类调和映射的非线性EulerLagrange方程的求解,采用类似Gauss-Seidel迭代的技巧设计数值格式。由于对初值的合理猜测和类调和映射作为泛函极小的性质,这样的迭代算法被证明是收敛的。在液晶聚合物的模拟中,Doi模型被普遍采用来描述方向分布函数的时间发展,它能够预测几种不同的聚合物现象。当聚合物浓度增加时,方向分布函数变得越来越趋向于Delta函数,要刻画它在很小的区域上就需要很多的计算网格。用这套移动网格方法,对于Fokker-Planck方程,我们可以求解球调和展开方法无法处理的大聚合物浓度O(1000)的情形。另外本文还根据方向分布函数自身所具有的独特性质设计了一种特殊的移动网格方法,可以大幅度地提高数值模拟的速度,数值例子中给出了对应较大聚合物浓度的相图。
其他文献
本文对多目标规划问题中函数的凸性和对偶问题进行了讨论。第一章介绍了我研究的思路和对这方面工作的一些看法。在第二章中,我给出了多目标规划问题有效解和弱有效解的定义,并
对角占优矩阵是一类有着广泛应用背景的特殊矩阵,它在数学、物理和工程技术等实际问题中出现的常微分方程、偏微分方程和大型线性系统的算法研究中有着十分重要的作用.尤其是
所谓函数的唯一性理论主要是探讨在什么情况下只存在一个函数满足给定的条件.近几十年来,它倍受关注,已成为国际上较为活跃的研究课题.而且,随着研究的不断深入和发展,它被赋予了
代数K-理论与代数数论有着密切的联系。假设F是一个数域,O〈,F〉是F的整数环。对于Tame核K〈,2〉O〈,F〉的结构的研究是热门的前沿课题之一,许多数学家对此进行了大量的研究。
本文研究了一个倒向随机微分方程的理论问题和一个数理金融中的实践问题。 在对一般鞅驱动的倒向随机微分方程的研究中,通过对驱动鞅过程和σ-代数流进行停时化处理,改善了
在经典小波分析的基础上,本文提出了n维广义拟实数进制小波分析的理论框架,增加了构造小波的自由度。同时给出了正交和双正交情形下的相关定理,证明了在此理论框架下,广义的Malla
非正则哈密尔顿系统广泛应用于物理力学诸多领域,譬如等离子体中的导心系统就是一个典型的例子.对于这类系统,传统意义上的标准辛算法不再适用.以往的处理方法有两种:一种是针
在数学中,图论是研究图的理论.图是离散数学中最重要的研究对象之一,它是用来描述集合元素间二元关系的一种数学结构.一个图是由图中的顶点以及连接它们的边所组成.具体来说,一
  文章分为四个部分。第一部分在文献阅读基础上回顾了协整理论的发展历程,并且总结了时间序列分析与协整分析的研究成果。第二部分首先阐述了维纳过程和基于其上的泛函中心
本论文给出了在一个较弱条件下的一种精确Gross-Zagier公式.特别的,我们还推广了Gross在文章”X0(N)上的Heegner点”中的工作,考察更为广泛的模曲线上的Heegner点的性质.通过运