HTR-10平衡态运行方式研究

为了使１０ＭＷ高温气冷实验堆（ＨＴＲ－１０）运行在安全、经济的状态下，研究了５次通过、８次通过和１０次通过三种运行方式下平衡态ＨＴＲ－１０堆芯的特性，利用高温气冷堆物理设计程序 ＶＳＯＰ对所选方案进行分析计算。结论表明：在最大燃耗不超过１０１ ＧＷｄ／ｔ的条件下，增大燃料球通过堆芯的次数并缩短每次通过堆芯所需的时间，将会使乏燃料平均燃耗提高，使ＨＴＲ－１０的燃料得到更有效的利用。     

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

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

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

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

计算流体力学与传热中的多重网格块修正算法

在流体力学的数值模拟计算中经常会遇到大型代数方程组的求解,而求解代数方程组的代价随着方程个数的增加而迅速增大。 主要将块修正多重网格算法引入ＳＩＭＰＬＥ通用程序中,并以流体力学中五个典型的实际问题为例,分别采用四种不同的迭代方法进行计算比较,展示了块修正多重网格法加速收敛速度的特性。     

HTR-10 大破口事故下结构特性的瞬态分析

１０ＭＷ模块式高温气冷反应堆（ＨＴＲ－１０）的热气导管压力壳的双端断裂事故属假想极限事故。当热气导管断裂后,高压的一回路冷却剂气体（３ＭＰａ）通过破口向堆舱猛烈喷放,一回路迅速卸压,并形成强卸压冲击波。利用了改进后的Ｋ－ＦＩＸ（ＦＬＸ）程序,对该事故下冷却气体喷放过程中堆体内主要结构部件（压力壳左支承、堆芯壳支承、堆芯壳）的安全特性进行了瞬态分析。通过计算,给出了破口处的压力瞬变、流量瞬变和堆芯壳上下两端的压差瞬变,以及在卸压冲击波作用下堆芯壳的膜应力和弯曲应力。最后,利用计算分析定量地给出了大破口极限事故下ＨＴＲ－１０的堆体主要结构部件的安全系数。结果表明,即使在极限事故下,ＨＴＲ－１０堆体结构仍有良好的安全特性。     

最小二乘等几何分析的多重网格方法

针对最小二乘等几何分析得到的代数方程系数矩阵的条件数大、迭代求解成本高的问题,提出了求解该方程的多重网格法。该方法在密网格上进行误差光顺,使高频误差快速衰减,在疏网格上进行误差修正,使低频误差快速衰减。通过节点插入算法自动生成不同尺寸的网格,根据离散B样条建立网格转换矩阵。采用该方法求解了泊松方程,对比了多重网格迭代与Gauss-Seidel迭代、PCG迭代的收敛性,结果表明Gauss-Seidel迭代收敛速度最慢,PCG迭代收敛速度随着代数方程自由度的增加而变慢,多重网格的收敛速度最快,能够有效求解最小二乘等几何分析得到的代数方程,解决了矩阵条件数过大的问题,并且收敛速度与网格尺寸无关。     

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

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

薄平板气动导数的数值计算

用基于一般曲线坐标系和交错网格的差分法求解原始变量二维不可压粘性流体的Ｎ－Ｓ方程,计算了宽高比为２２．５、两端圆滑的薄平板的气动导数．将薄平板当作刚性物体,首先分别计算物体做竖向强迫振动和扭转强迫振动的气动力,雷诺数取为１０５,物体运动和气流的耦合作用用动网格法考虑,Ｎ－Ｓ方程用二阶投影法和多重网格法求解．然后用最小二乘法确定Ｓｉｍｉｕ和Ｓｃａｎｌａｎ给出的气动力表达式中的８个气动导数．绝大部分气动导数的计算结果和Ｔｈｅｏｄｏｒｓｅｎ的理论值及实验值都十分吻合     

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

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

HTR-10初装堆芯及过渡过程物理计算分析

选取石墨球与燃料球均匀混合作为１０ＭＷ高温气冷堆（ＨＴＲ－１０）的初始装料方案，利用高温堆物理模拟程序ＶＳＯＰ及二维ＳＮ程序，分析计算了初始装料时ＨＴＲ－１０堆芯进水反应性效应、控制棒及第二停堆系统反应性当量，研究了初装堆向平衡态过渡过程中的临界性、单球最大功率、最大比燃耗等变化情况。结果表明：ＨＴＲ－１０初装堆的进水反应性效应比平衡态小；控制棒及第二停堆系统反应性当量比平衡态的大。但是，初装堆冷态下反应性控制系统当量裕量比平衡态小；过渡过程中有效增殖因数在很小范围内变化，燃料最大比燃耗不超过１００ＧＷｄ／ｔ。     

基于流体力学方程的二维烟雾模拟

介绍了一种基于流体力学方程并采取欧拉法实时模拟二维真实感烟雾的算法.为了获得视觉真实感,烟雾流体场的物理模型是根据Navier-Stokes方程定义;采取欧拉法将烟雾流体场划分为二维网格空间,并将烟雾速度场和浓度场分布定义于每个网格中心点;最后,通过双线性插值获得整幅图象中每个象素的烟雾浓度.试验证明该算法能实时高效地模拟具有真实感的二维烟雾效果.     

P_1非协调四边形元解Stokes问题的多重网格算法

研究了用P1-Q0元(其中P1表示P1非协调四边形元)解Stokes问题的多重网格算法.由于P1-Q0元不满足LBB条件,因此其不能直接用来求解Stokes问题.本文基于曾提出的一种P1-Q0元解Stokes问题的非协调混合有限元稳定化逼近方法,提出了W循环多重网格方法,证明了该方法的最优收敛性.最后给出的数值算例验证了该理论结果.     

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

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

基于有限分析法和高阶紧致格式的新计算方法

将有限分析法和高阶紧致格式两种计算方法相结合,并提出了一种新的压力方程修正方法,发展了一种新的流体力学计算方法,并用它求解三维N-S方程组,计算结果表明,此算法较适应于复杂流动问题的计算。     

N S方程全离散非线性Galerkin格式异步并行算法的计算稳定性分析

给出了Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程全离散非线性Ｇａｌｅｒｋｉｎ格式,提出了一种新的异步并行算法,并对其计算稳定性进行分析。文中证明了此格式的计算稳定性。     

SIMPLE算法中压力修正方程边界条件确定新探

不可压缩粘性流体Ｎａｖｅｒ－Ｓｔｏｋｅｓ方程ＳＩＭＰＬＥ算法求解中压力修正方程起着很关键的作用。但其边界处理的特殊性常给正确应用带来的一定困难。提出从连续性方程本质出发确定压力修正方程边界条件的新方法。     

用代数应力紊流模型预估水轮机转轮内部三维流场

以三峡水轮机组转轮内部流动研究为背景,基于雷诺时均方程和代数应力紊流模型,采用贴体坐标和交错网络系统,用ＳＩＭＰＬＥＣ算法对一比转速为２４０的混流式水轮机模型转轮内部三维紊流进行了数值模拟,得出了设计工况的流场特性,并初步分析了转轮的汽蚀和能量特性。结果表明,研制的程序可作为水轮机组转轮内部流场分析、数值试验和性能预估的工具。计算结果可为改进和优化转轮设计提供有益的参考资料。     

格子Boltzmann方法中的几个问题

简单介绍了在格子Ｂｏｌｚｔｚｍａｎｎ方法的最新动态,分析了该方法目前存在的问题。给出了几种解决问题的途径,通过数值分析和理论分析得到,用格子Ｂｏｌｔｚｍａｎｎ方法求解流体力学问题,有三方面的优点：（１）格子Ｂｏｌｔｚｍａｎｎ方程是 维的,对流项是线性的,简单的松驰方程可以得到复杂的非线性的流体力学方程,已表现出数值稳定,结果准确等特点。（２）在不可压的极限下可以得到不可压缩Ｎ－Ｓ方程,压力由状态方     

基于多重网格法的相位解缠算法

指出了最小二乘相位解缠算法是求解二维相位解缠问题最稳健的方法之一,并可等效为求解一大型的稀疏线性方程系统.求解大型线性方程组通常采用迭代法,然而其收敛速度非常慢.为了改善收敛特性,提出了一种新的相位解缠算法——多重网格法,该方法通过在疏密不同的网格层上进行迭代,以平滑不同频率的误差分量,从而加快系统的收敛速度.仿真实验表明:该方法能够很好地恢复真实相位,具有解缠精度高,收敛速度快等优点.     

全三维粘性流场计算的快速收敛方法

以新的高收敛率的 L U型隐式格式和高精度、高分辨率的 MU SCL TVD迎风格式为基础 ,提出了一种快速求解三维粘性流场的计算方法。为有效减少计算工作量 ,近固壁区采用了壁面函数方法。湍流模型采用简单的混合长度模型。为提高迭代收敛速度 ,从多重网格理论出发 ,结合改良型L U隐式和迎风格式 ,对 Reynolds平均三维可压缩 Navier-Stokes方程进行求解。用此方法数值模拟了 NASA L ewis37跨音速压气机转子流场。计算得到的流场与实验结果进行了对比。结果表明 ,此方法可快速得到三维粘性流场的流动特性且计算结果精度较高 ,可望在叶轮机械多叶片排全三维粘性流场校核计算中发挥作用。