【摘 要】
:
线性方程组的求解开销往往是实际复杂应用在数值模拟时的主要开销。预处理迭代方法是求解大规模稀疏线性方程组的常用求解方法,常见的预处理方法和迭代方法往往聚焦于方法的通用性能而缺乏对于实际应用数值模拟特征的考虑。惯性约束聚变是一类强非线性、强间断、大变形、多介质的辐射流体应用,其在数值模拟时表现出各种特征:一方面,在模拟的一段时间内,一些物理量在局部计算区域内发生剧烈的变化,而在其他区域内变化不大;另一
【基金项目】
:
国防基础科研核科学挑战专题:百亿亿次计算科学的计算方法与高效能实现”项目(项目号:TZZT2016002); 国家自然科学基金:高阶并行可扩展 CFD 应用开发框架研究” 项目(项目号:No.61872380);
论文部分内容阅读
线性方程组的求解开销往往是实际复杂应用在数值模拟时的主要开销。预处理迭代方法是求解大规模稀疏线性方程组的常用求解方法,常见的预处理方法和迭代方法往往聚焦于方法的通用性能而缺乏对于实际应用数值模拟特征的考虑。惯性约束聚变是一类强非线性、强间断、大变形、多介质的辐射流体应用,其在数值模拟时表现出各种特征:一方面,在模拟的一段时间内,一些物理量在局部计算区域内发生剧烈的变化,而在其他区域内变化不大;另一方面,强间断和大变形等特点使得其离散所得到的系数矩阵元素大小存在量级上的差异。水下航行器则是一类不可压流体应用,在数值模拟中通常使用分离迭代方法求解不可压NS方程,分离迭代算法使得该应用在模拟时表现出数值震荡的特点。本文主要瞄准惯性约束聚变和水下航行器两类应用,利用其数值模拟特征,驱动模拟中所产生线性方程组高效求解。本文的主要贡献和创新如下:(1)针对数值模拟应用的局部特征,本文提出了一种基于局部特征的线性方程组求解方法。该方法是一种代数方法,主要包含三个步骤:首先,提取变化剧烈的局部区域;其次,求解局部区域对应的局部线性子系统;最后,求解整体线性方程组。本文对局部特征进行了数学抽象,给出了局部特征线性方程组定义和相关性质;并给出了基于梯度和基于残差两种局部区域提取方法。本文在二维热传导方程、多群辐射扩散方程、三温能量方程等问题中验证了该方法的有效性,该方法在多群辐射扩散方程和三温能量方程的线性方程组测试集中分别能达到1.61倍和1.65倍的加速比。(2)针对数值模拟应用的多尺度特征,本文提出了一种基于多尺度特征的预处理矩阵元素过滤方法。该方法首先根据元素的相对大小来建立求解变量之间的弱依赖关系矩阵,然后根据弱依赖矩阵以及一定的过滤策略来删除预处理矩阵中的元素。本文提供了四种过滤策略,包括双侧对称过滤、单侧非对称过滤以及两种相应的对角线修正策略。本文在泊松类方程、多群辐射扩散方程、三温能量方程等问题中验证了该方法的有效性,该方法在多群辐射扩散方程和三温能量方程中的线性方程组测试集中分别能达到1.47倍和1.55倍的加速比。(3)针对数值模拟应用的间断特征,本文提出了一种基于混合粗化策略的代数多重网格预处理方法。该方法的粗化算法主要包含两个步骤:首先,使用经典的粗化算法获得粗网格;其次,使用相容松弛迭代来衡量所得粗网格的质量,并挑选迭代中收敛较慢的细点作为粗点集合的补充。本文在间断系数的泊松类方程、三温能量方程以及三维翼身融合模拟等问题中验证了该方法的有效性。(4)针对分离迭代算法的数值震荡特征,本文提出了一种基于分离迭代算法特征的初值优化方法。该方法利用加权分组插值技术对分离迭代算法所产生的线性方程组迭代初值进行优化,主要包含三个步骤:首先,根据分离迭代算法配置特征,将所有线性方程组划分为若干泳道,并将各泳道中窗口范围内的已知解划分为若干小组;其次,每个小组内利用已知解进行初值预测;最后,将各小组的预测解进行加权平均并作为新的迭代初值。本文验证了该方法在pitz Daily、二维NACA0012、三维翼身融合等案例中的有效性,该方法在二维NACA0012和三维翼身融合等案例的数值模拟中分别能够获得2.58倍和1.87倍的加速比。
其他文献
蒎烯是一种自然界中存在的双环单萜类化合物,具有广泛的应用。蒎烯的两种同分异构体α-蒎烯和β-蒎烯二聚化后会形成具有极高能量密度的化合物,能够充当航天燃料。生产蒎烯的传统工业方法是作为造纸和提炼松节油的副产物,过程中需要消耗大量的原材料,且产量较低。为降低蒎烯的合成成本、提高产量,开发新的合成途径是十分必要的。目前,已有人构建了利用大肠杆菌合成蒎烯的方法,但所达到的滴度还只有32 mg/L[1]。对
量子多体系统在凝聚态物理和量子信息的研究中相当重要,许多神奇的现象出现在相互作用的多体系统。在拓扑相变的相关研究中,自旋系统扮演非常基础和重要的角色。例如,物理学家最开始发现的拓扑相变就是在XY自旋系统中的Berezinskii-Kosterlitz-Thouless(BKT)相变。在本论文中,我们主要研究一维Ising模型和XXZ模型的Berry相位和拓扑相变。本论文的主要框架如下:第一章,指出
一次飞行访问多个目标因其探测效率高、平均成本低而备受各大航天机构的青睐,是未来开展在轨服务和深空探测等航天任务的主要形式。多目标访问任务尤其是面向大规模目标访问任务的设计是一项极具挑战性的工作,其中访问序列的规划是重点与难点。为了提高复杂空间目标访问任务的设计水平,本文对面向大规模目标访问任务的飞行序列规划方法进行了研究,主要的研究内容和研究成果如下:提出了基于多神经网络的摄动多脉冲最优交会速度增
含能结构材料是一种兼具力学性能和能量释放特性的材料。以含能结构材料取代传统惰性金属材料,实现战斗部全能量化是现代高效毁伤武器的发展趋势。Al基复合材料和Zr基非晶合金因其优异的力学性能和能量密度,是当前含能结构材料的研究热点,但同样存在各自的技术和理论问题。对于Al-Ni等Al基复合材料,其冲击反应释能效果不理想;而对于Zr基非晶合金,则难以实现大尺寸成型。为提升Al-Ni复合材料冲击释能性质,并
TaxHf1-xC固溶陶瓷作为Ta C和Hf C的固溶体,是目前已知物质中熔点最高的一类材料,且相比于单一的Ta C、Hf C二元陶瓷而言,TaxHf1-xC固溶陶瓷还具备了更高的硬度和更优的抗氧化、耐烧蚀潜质,有望继Zr B2、Hf B2等热门超高温陶瓷材料之后成为未来热结构材料领域新的候选材料。但由于Ta C和Hf C陶瓷的熔点高且自扩散系数低,制备高致密且具有理想单相固溶结构特征的TaxHf
社会系统、信息系统、军事系统等典型复杂系统呈现出显著的层次性、差异性及动态性特征,传统的单层网络模型已无法充分描述以上复杂性,研究适用于层间耦合、结构差异、动态演化的多层网络模型尤为必要。网络鲁棒性是网络科学核心问题之一,多层网络鲁棒性更具挑战,已成为网络科学发展亟待解决的前沿课题。本文针对多层网络结构鲁棒性,重点开展多层网络的耦合作用机制、攻击级联失效以及结构状态恢复等研究。具体包括以下四个方面
本文以超燃冲压发动机中的液体燃料射流为研究对象,以实现超声速气流中液体燃料射流喷注、雾化、混合、蒸发及燃烧完整过程的数值模拟为基本目标,建立了两相燃烧大涡模拟方法,并对超声速气流中横向喷雾混合及燃烧过程中涉及的基本物理过程及机理进行了讨论分析。首先基于欧拉-拉格朗日框架建立起一套适用于超声速气流中两相流动的大涡模拟方法。通过基于网格控制体建立的动态数据链表实现了对液滴的高效管理以及液滴在变形网格条
在过去的四十多年中,热带气旋(TCs)的路径预报水平得到明显提高,但是TCs强度,尤其是快速增强过程的预报水平却提高缓慢。这主要是因为TCs的强度变化受到复杂的内部相互作用和外部应力的共同影响,如海气相互作用。作为中国近海中TCs最活跃的海域,南海的热动力结构复杂、观测不足,且经过的TCs发展和登陆时间很短,因此对其中快速变化TCs的预报造成了巨大挑战。所以,认识海洋与TCs之间的相互作用机制对提
信息技术的深入发展和网络应用的不断丰富,加上智能终端的广泛使用,使得社交网络成为人们发布和获取信息的重要媒介。借助于社交网络平台,人们在现实世界中的各种复杂关系在虚拟网络上都能得到体现与延伸,同时,在真实世界中发生的事件也会以网络信息的形式借由用户之间的线上交互进行扩散,进而对人们的线下行为产生影响。作为社交网络的主体要素,用户是网络信息的生成者与传播者,其表现出的社交影响力不仅是信息传播与关系结
随着互联网技术的飞速发展,人类社会可以利用的数据正以前所未有的速度增长,宣告了大数据时代的正式到来。作为数据分析的重要工具之一,复杂网络理论自诞生以来就受到了来自各个领域研究学者的广泛关注,特别是随着在线社交网络的兴起,利用复杂网络理论分析用户的行为模式已体现出其巨大的经济和社会价值。研究复杂网络的主要目的就是为了解决网络上的动力学问题,本文基于复杂网络传播动力学相关理论,分别开展了关于影响力极大