偏微分方程迭代并行解法与网格优化方法

来源 :中国工程物理研究院 | 被引量 : 4次 | 上传用户:zybp821
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
科学计算对计算规模的要求是无止境的,随着计算机能力的飞速的发展,计算的规模已经从几千几万到百万甚至千万上亿的规模。算法也相应地由串行的算法发展到适合大规模并行计算机的并行算法。计算规模的扩大带来了许多新问题。首先,并行计算机的体系结构要求相应的算法能充分挖掘计算能力,促进了并行计算方法的研究;其次,大规模的数值模拟计算对于计算方法的稳定性和计算网格的品质提出了更高的要求,推动了算法的稳定性和流体计算中计算网格的优化研究。因此,本文的主要研究目标是针对偏微分方程大规模的数值模拟中出现的这些问题进行研究,提出解决的方法和进行理论分析。 现代并行计算面临许多问题。从并行机的体系结构来说,大规模分布式内存并行计算机(MPP,或PC Cluster)是发展的主导方向,其主要特征是采用消息传递进行不同处理机之间的通信。现代大规模并行计算面临的主要瓶颈是全局通信问题和同步问题。全局通信体现了整个问题的全局综合,这是种全局性无法避免,但可以通过别的途径使之达到相同的功能,比如通过迭代方法。在构造偏微分方程离散计算格式时,应尽可能考虑避免整体数据的强相关和无需全局通讯的具有高并行度的格式。研究突破代数的层面,从差分格式乃至偏微分方程,甚至物理问题的层面来提出新的算法。偏微分方程的差分格式是离散逼近的重要的一步,其并行性方面的研究最近二十年来得到广泛的重视和发展。但离散格式的构造和格式的代数求解没有充分地结合起来,这两者相结合的问题没有受到足够的重视。实际上,差分格式的成功和代数方程组的求解息息相关,只有把差分格式和方程组的迭代求解结合起来研究,才可能使代数方程组的一些比较有效的迭代方法发挥更大的作用,构造真正适合应用的差分格式。大规模科学计算的另外一个问题是对于高速流动的流体,当用Lagrange方法或ALE方法进行计算时,容易出现网格折叠的情况而中断计算,这是精密物理的大规模数值模拟中目前遇到的瓶颈问题之一。 针对以上的问题,本论文对偏微分方程的差分数值求解进行了讨论。针对在现代并行计算机条件下如何充分有效地利用大量的计算机资源,高效稳定地数值求解偏微分方程,进行了详细的讨论。本论文的特点是结合迭代算法研究差分格式和结合差分格式研究迭代算法,从而使两者有机地结合到一起,为更好的解决辐射流体热传导和粒子输运计算等问题,探索新的研究途径。对区域分解的并行算法进行了理论分析和数值试验比较研究,将这些方法应用于抛物型方程和中子输运方程,得到了一些有创新性的结果;对代数方程组的求解的红黑排序混合算法进行了详细的理论分析和数值试验,并将其应用于辐射流体力学计算和油藏模拟中,取得了满意的计算结果;结合并行差分格式对Krylov子空间方法提出了一种并行预条件技术,并对其谱性质进行了较深入的分析;对于流体计算网格的优化,提出了一种组合优化的方法,这种优化方法的计算量小,并且能优化凹网格的问题,同时能较好的保持优化过的网格的重心能匹配原来网格的重心,能比较有效的解决网格的优化问题。 全文的主要内容可以大致分为三个部分:并行差分方法,线性代数方程组的并行高效求解和姆格优化这花个部分既相互独立,同时又相互影响和相互促进。木文结构如下:首先是绪论,概括介绍了本文的研究背萦和主要内容;第一章、第二章、式的钩造技术,主要研究了迭代并行差分格式的稳定性和收敛性,第说章研究并行差分格以及结合迭代技术的迭代并行差分格式的构造和理论分析;第四章研究中子输运方程离散纵标方法的迭代并行差分格式:第五章讨论了红黑排序混合算法的收敛速度问题,以及该算法在实际工程计算中的应娜:第六章是墓于界荡预估和修正的并行预条件技术。第七章讨论了计算流体的A LE方法中的卿格重构和优化技术。本文的杰要_f作如一F:口}给出DFF差分格式的一个2范数的严格的无条件稳定性的证明,并估计出稳定性的常数, 同时证明了用OFF差分格式构造的区域分解并行差分格式的稳定性的必要条件。阴对显式迭代并行差分格式证明了其对于模型问题的稳定性条件比纯显式差分格式放宽了 封仍以_七。睁}从空间区域分解出发,综合考虑迭代算法和差分格式的相互关系,使得两者相互触合,取 长补辣,设计了高效的迭代并行差分格式,这些新构造的并行格式其有良好的性质,主要 表现在如下豹几个方面: 给一)差分格式豹稳定性非常好。 (二)差分格式有较高的精度。 (蕊)差分格式构造比较简单,对区域的剂分没有任何限制。 (洲)差分格式不需要全局通信,有很好的扩展性。 (五)便于对已有的程序进行改造。拼〕将构造抛物型方程并行差分格式的方法推广到中子输运方程,进行了数值试验。吓!对线性代数方程组的迭代算法,红黑排序混合算法,分析了其算法特性和算法的收敛加速 悄况。证明了红黑排序混合算法比自然排序的算法,!讨样用共辘梯度法和Jacob玄迭代方法 求解,混合算法收敛速度是自然排序的2倍;用GMRES算法,混合算法收敛速度是自然 排序算法的2到4倍。大金的数徽试验验证了理沦?
其他文献
在海河流域废污水处理及再生利用现状调查的基础上,分析了海河流域废污水再生利用存在问题。
送股又称为股票股利,是股利政策的重要组成。中国市场上送股概念又分为送红股和转股两种,前者来源于公司的留存收益,后者来源于公司资本公积,美国市场送股则按比例划分。本文
1979年至2010年,中国实际利用外商直接投资10483.86亿美元,平均每年327.62亿美元,年均增长26.69%。2010年中国实际吸收外资金额突破千亿美元,达到1057.4亿美元,连续18年居发
随着飞行器控制及导航技术的迅猛发展,饱和问题越来越受到人们的关注。其主要体现在两个方面,一方面是控制系统的执行器饱和,其表现在执行器工作过程中的输出有界性。另一方
镁合金具有密度小、比强度和比刚度高、尺寸稳定性好等特性,其产品在电池、交通工具、航空航天、武器等领域备受青睐。工业挤压态镁合金具有密排六方晶体结构,而且晶粒较为粗大
随着汽车工业的高速发展以及对汽车零部件减重的迫切需求,开发工艺简单、价格低廉的高强度钢铁材料是降低油耗、减少排放、提高安全性的必然趋势和根本途径。目前重型商用车普
衡量各发电商的市场力与投机行为 ,研究市场的电价变化 ,是电力市场的深层次问题。在总结传统市场力理论的基础上 ,基于厂商所面对的需求曲线 ,提出了反映机组市场力的市场电
随着当前各种新型技术的迅速发展,企业人力资源管理的信息化发生了根本的变化,]--h资源管理信息化属于企业的一项重要课题。本文分析了企业人力资源管理中信息化所出现的问题,简
为研究现浇和预制装配式混凝土扇形铅黏弹性阻尼减震框架抗震性能的差异性,对二层二跨预制装配式和现浇混凝土框架扇形铅黏弹性阻尼减震框架结构进行了低周反复加载试验,对比
面对日益突出的环境和资源问题,各国逐渐加强技术革新以实现节能减排的目标。以火力发电为主的能源结构迫使技术突破的重心转向如何有效地提高燃煤发电机组的效率。燃煤的转