非线性椭圆问题的非精确牛顿代数多重网格法

采用基于矩阵图集的粗化算法形成粗点集,构造改进的插值算子,结合V型多重网格法和瀑布型多重网格法的算法结构,提出了一种改进的代数多重网格（IAMG）法,并估计了该算法的计算量。将IAMG法运用于求解牛顿算法中线性校正方程,提出了求解非线性椭圆型问题的非精确牛顿代数多重网格（IN-AMG）法。数值实验表明与对比算法相比,IN-AMG法在求解线性校正方程方面的整体计算量更少、计算时间更短。     

求解变分不等式的多重网格算法

多重网格法是求解椭圆型偏微分方程边值问题的一种快速、有效的数值方法.本文将多重网格算法应用于变分不等式问题的数值求解.将不动点法与多重网格过程相结合提出了求解变分不等式问题的一种多重网格算法.以障碍问题及其特例—弹、塑性杆的自由扭转问题为例,给出了求解所得的数值结果,讨论了这种算法的收敛性情况.实例表明,文中提出的算法保持了一般多重网格过程的主要特点.它具有远小于1的收敛比率;松弛因子的改变对收敛速率的影响很不灵敏;求解变分不等式问题的计算量接近或略小于相应的变分问题.     

三维弹性力学问题中有限元方程的预处理方法

针对三维弹性问题中有限元方程的数值求解，建立了一类简单且实用的代数多重网格预处理共轭梯度法(AMG-CG法)，详细描述了相应代数多重网格方法的粗化技术及网格转移算子的构造．由于该预处理方法能有效地降低刚度矩阵的条件数，使刚度矩阵的谱分布更集中，从而大大提高了计算效率．数值结果表明，AMGCG法对求解三维弹性问题有限元方程是十分有效和健壮的。     

二维非结构网格Euler方程高效解法

在非结构网格上应用多重网格技术加速 Euler 方程的收敛,在多重网格中通过聚合法进行粗网格生成,并对粗网格中的多边形网格做了等价面处理.在空间离散上采用 Roe 格式,在时间推进上分别采用了显式和隐式算法.通过对 NA-CA0012 翼型和 RAE2822 翼型的流场模拟,比较了显式多重网格法和隐式多重网格法的计算效率.     

一维抛物方程初边值问题拟多重网格预处理迭代法

本文将求解椭圆方程边值问题的拟多重网格预处理迭代法推广到求解抛物方程初边值问题,将多重网格法的优点和预处理方法很好的结合到一起,加快迭代的收敛速度,从而减少解抛物方程的计算量。     

多重网格法求解原始变量形式的Navier-Stokes方程

采用多重网格法 ,求解原始变量形式的Navier Stokes方程 .控制方程在交错网格上离散 ,并采用SIMPLE算法计算 .为加速收敛 ,结合使用适合非线性方程的全近似格式的多重网格法FAS .文章计算了雷诺数Re =10 0 ,10 0 0时不可压粘性流体绕圆柱的流动 .计算结果与实验结果和国外的计算结果基本吻合 ,说明了该方法的准确性和有效性     

基于多核异构的代数多重网格的并行算法实现

近年来,受GPU其高浮点峰值性能的提高和应用领域中大规模科学计算问题的驱动,高性能领域中利用代数多重网格(AMG)求解稀疏线性方程组成为研究热点。针对经典的AMG算法,探究建立阶段(网格粗化)和求解阶段的并行计算结构,提出基于多核异构的AMG并行计算模式。数值实验表明,并行计算模式计算效率相对于串行提高了3~4倍,加速效果明显。     

基于多重网格法的MT正演模型边界截取

为了验证基于多重网格法边界截取的效果,通过数值模拟,分别对大地电磁法问题的均匀半空间模型和三层模型进行基于多重网格法的边界截取研究和分析.采用加权平均法产生粗网格上系数矩阵,并与传统的方法对比.研究结果表明:粗网格加权平均法的渐进收敛性明显比Galerkin法(采用简单的插值矩阵)的好,但比粗网格近似的几何法的渐进收敛性略差;收敛性的广义Fourier谱分析在大地电磁法中对介质的非连续性不敏感,不能区分模型参数变化所引起的收敛性变化;随着边界截取的增加,满足特定收敛标准的数值解与理论解偏差增大;当截取到的深度约为趋附深度的2倍时,截取后的计算结果接近理论值,说明基于多重网格法的边界截取有效.     

一类新的瀑布型代数多重网格方法

对瀑布型多重网格(CMG)法和代数多重网格(AMG)法进行组合,提出一种新的求解二维椭圆型边值问题的瀑布型代数多重网格(CAMG)法,并进行数值实验.结果表明,CAMG法所得解的误差小于10-6,并且每层的迭代次数都少于AMG法,特别在最细层上的迭代次数远远少于AMG法.CAMG法是收敛,高效的迭代算法.     

多重网格法及在HTR-10堆芯动态仿真

研究了ＨＴＲ－１０堆芯的动态仿真,为反应堆控制系统的分析和设计服务。由于堆芯模型方程是典型的非定常可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,要求快速准确的求解。求解算法采用了结合多重网格的ＳＩＭＰＬＥＣ算法,并且讨论了多重网格法的基本思想、算法步骤和计算流程。通过对求解域中的不规则区域的处理、合理的进行粗细网格的划分和校正、设置出口缓冲区等措施,克服了在ＨＴＲ－１０堆芯动态仿真中运用多重网格法的困难,并且所得的结果与单网格算法进行了比较。仿真结果表明：达到同样精度,采用多重网格方法,求解过程比单网格算法的迭代次数减少了约５／６,计算时间减少了约一半,为快速进行堆芯动态仿真提供了条件。