抛物问题的时空局部网格加密方法

来源 :山东大学 | 被引量 : 2次 | 上传用户:minifeng
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  一些大规模的物理过程往往在空间和时间上具有很强的局部性质,由井、裂缝、障碍、区域边界等引起的局部性质在空间上是固定的,还有一些情况下的局部性质是随着时间而变化的,如油田开发过程中的流体驱替前缘。这些物理过程都可归结为求解大型偏微分方程组。如果使用均匀的网格系统,要得到关于这类问题的较精确的解,需要使用足够小的空间步长和时间步长,这就大大增加了计算成本,从而限制了所解决的问题的规模。自1982年起,国外文献中陆续提出了解决这类问题的局部加密方法[1 5],提出的局部网格加密方法仍然使用原来的网格系统,但在需要加密的地方,即局部性质较强的区域,对网格进行加密,这样就可以在不太增加计算量的前提下,大大提高计算精度。目前,局部网格加密方法已被广泛的应用到科学工程计算领域,同时在理论研究方面,许多学者进行了大量的研究工作,提出了多种局部网格加密离散格式。  对时空局部网格复合网格的形成,我们采用规则的局部网格加密方法[6,7]:首先在整个区域引入总的时间步长和空间步长,然后利用自适应网格加密方法或根据解的先验信息,选择一些子区域,在这些子区域上引入局部时间步长和局部空间步长,且总的时间步长和空间步长分别为局部时间步长和空间步长的整数倍,这样就形成了时空复合网格。在这种情形下的第一个问题是构造稳定和精确的近似格式,困难主要在粗细网格区域的交界面;第二,我们需要构造求解复合网格上代数系统的有效的求解方法,文献[8]中Dawson和Du基于Galerkin离散,给出了区域分解算法,[9]中Ewing、Lazarov、Pasciak和Vassilevski利用间断Galerkin方法,构造了离散格式,并利用迭代方法求解相应的复合网格系统,这部分的内容不在本文讨论范围之内,可参考文献[8 10]。  本文针对局部网格加密方法,作了部分研究工作。首先,针对几类问题,在矩形剖分或三角剖分时空局部加密复合网格上,建立了有限差分格式,从理论上研究了离散格式的稳定性和收敛性,并通过数值试验验证了理论分析结果;然后,把局部加密方法应用到多轮次调剖数值模拟器中,通过数值模拟结果证实了局部网格加密方法的有效性。这些研究结果对于油气资源的勘探、开发,环境科学的数值模拟和计算,以及半导体器件问题的数值模拟,均具有重要的理论和实用价值[11 13]。  全文共分三章。  本文第一章,研究了半导体器件问题,提出了在矩形剖分时空局部加密复合网
其他文献
本文首先回顾了常微分方程理论,特别是可积性理论的发展历史及研究现状,介绍了Liouville可积的相关概念与定义,以及形如dy/dx=∑ni=0fi(x)yi的“方程在线性变换群下的不变量”等
  两相流体动态边界在许多科学应用中都发挥着重要的作用[1],例如LOX/H2火箭发动机。动态边界决定了原子化过程,比如液体喷射层打破变成小水滴,这个动态过程可以分为两个连续
  数学物理及工程问题,如油气藏的勘探与开发,大型结构工程,天气预报等,无不归结为求解大型偏微分方程,面临着如何有效解决大规模科学计算的困难。众所周知,区域分裂能够降
  本文主要考虑的是:若随机变量ξ,η∈L2(Ω,FT,P;R),且ξ和η同P-分布,那么对泛函εg[·]而言,εg[ξ]=εg[η]是否总成立呢?事实上,当且仅当εg[·]退化为E[·]时,上式总成
  广义大系统的稳定性和镇定问题是广义大系统理论的基本问题。由于广义系统的特殊性,对其稳定性和镇定问题的研究不仅要考虑稳定性,还要考虑正则性。而广义系统的鲁棒稳定性
分形插值替代传统的插值技术,给出了一个更广泛的插值函数集,它为理解现实世界的现象提供了一种很好的确定性方法。用这种方法,我们不仅能构造非整数维的插值函数,而且也能够构造
度量空间中的不动点理论是非线性泛函分析的重要组成部分,一直是人们研究的热门领域.自锥度量空间被提出以来,其空间结构与性质,尤其是锥度量空间中的不动点定理,更成为近年来国