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

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

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

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

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

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

采用多重网格法进行电子结构的计算

采用多重网格法和对称性理论以及密度泛函理论中著名的Ｋｏｈｎ－Ｓｈａｍ方程,开发了一种计算具有一定对称性的原子簇的电子结构的方法,并对Ｈ２分子的基态能进行了的计算。计算结果与用其他理论计算得到的结果符合的很好。     

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

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

低雷诺数下方柱和圆柱涡致振动的数值分析

通过数值求解二维不可压粘性流体的N-S方程,研究了方柱和圆柱在低雷诺数（Re＝90～130）下的涡致振动．柱体振动对流体的影响用动网格法模拟,N-S方程采用时步分裂技术先求出中间速度场,再用多重网格法求解泊松方程得到压力场,最后根据压力场对中间速度场进行修正．柱体的振动方程用Newmark-β法求解．通过计算柱体和流体的相互作用,成功地捕捉到了“锁定”和“拍”现象．计算结果与实验结果相当吻合,并和已有的用任意拉格朗日-欧拉公式有限元法计算的结果进行了比较．     

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

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

低雷诺数下方柱和圆柱涡致振动的数值分析

通过数值求解二维不可压粘性流体的N－S方程研究了方柱和圆柱在低雷诺数（Re一叨一130）下的涡致振动．柱体振动对流体的影响用动网格法模拟．N－S方程采用时步分裂技术先求出中间速度场,再用多重网格法求解泊松方程得到压力场,最后根据压力场对中间速度场进行修正．柱体的振动方程用Newlllaix一月法求解．通过计算柱体和流体的相互作用,笔者成功地捕捉到了“销定”和“拍”现象．计算结果与实验结果相当吻合,并和已有的用任意拉格朗日一欧拉公式的有限元法计算的结果进行了比较．低雷诺数下方柱和圆柱涡致振动的数值分析@曹丰产@项…     

求解二次Lagrangian有限元方程的瀑布型多重网格法

针对二次Lagrangian有限元方程,通过将新外推公式和二次插值技巧相结合,为细层提供好的初始值,设计了新的瀑布型多重网格法.数值实验表明,与基于部分几何信息的代数多重网格法相比,新算法有更好的精度和效率.     

Mortar型旋转Q1元的瀑布型多重网格方法

展现了mortar型旋转Q1元的瀑布型多重网格方法.证明了采用共轭梯度作为光滑子的瀑布型多重网格法是最优的,而采用其它传统迭代作光滑子的瀑布型多重网格法是拟最优的.并通过数值试验验证了我们的理论结果.     

圆柱绕流的三维数值模拟

利用计算流体力学软件CFX－4,对粘性不可压缩流体的圆柱绕流进行了三维数值模拟,采用有限体积法和SIMPLE计算程式,利用不可压缩Navier-Stokes方程,模拟雷诺数在亚临界区内的绕流流动,并计算了流体的水动力特性。为克服数值模拟高雷诺数时的数值不稳定性,计算中采用了QUICK迎风格式,其对流项为三阶精度,其余项如扩散项等为二阶精度,圆柱两端边界采用周期性边界条件。计算结果表明,高雷诺数时圆柱周围的流动具有明显的三维特性,且沿柱长方向不同断面的升力和阻力系数并不相同。同时,对圆柱绕流进行了二维数值模拟,并与三维数值结果进行比较,发现三维模拟的升力和阻力系数均小于二维模拟。     

水下圆柱体带空泡轴向绕流研究(英文)

研究了水下圆柱体带空泡轴向绕流问题.用数值方法研究了不同长细比圆柱体水下运动产生空泡的形态及圆柱体水下带空泡绕流的阻力特性,圆柱体绕流空泡形态以及阻力系数数值计算结果与试验对比,符合较好.圆柱体长细比对空泡形态影响可以忽略.随长细比增大:空化数较高时,阻力系数先减小后增大;空化数较小时,其阻力系数基本保持不变.小空化数条件下,给定长细比,空化数越小,阻力系数越小.     

绕横向振动圆柱分离流动的数值研究

本文用由作者提出并发展的一种基于区域分解思想.综合了解N—S方程的有限差分法及涡法各自优点的新数值方法,研究了均匀流绕横向振动圆柱的分离流动,考察了圆柱振动对涡旋脱落,升力,阻力等的影响.成功地模拟了“锁定”这一重要的液固相互作用的非线性现象。计算结果与国外实验结果和其他计算结果符合较好。     

二维圆柱层流绕流及其控制数值模拟

采用有限体积法求解二维N-S方程,通过改变时间、空间计算参数,对雷诺数100的二维圆柱非定常流场进行了数值模拟,将所得结果与实验数据以及已有的计算结果进行了对比.并对在尾部多个位置加入不同长度阻隔板控制流动的圆柱流场进行了考察.结果发现:在合适的位置放置阻隔板,可以有效抑制尾迹中涡脱落的形成.     

用速度-涡量方法解具有表面抽吸的圆柱绕流问题

用有限差分法计算了不可压缩粘性流体绕具有表面抽吸圆柱的流动，控制方程采用涡量－速度形式的N－S方程，并引入对数有坐标变换，以使近壁处的网格加密，其中涡量输运方程采用二步R－K方法求解，速度泊松方程采用LSOR方法求解，并对速度进行了无散修正，计算表明合理选择抽吸的位置与强度，可有效减少阻力和升力，得到较好的流动控制效果。     

圆柱绕流的边界层及其电磁控制

利用电介质溶液中圆柱体侧表面附近分布的电磁产生电磁力,作用于流体边界层,达到控制圆柱绕流的目的。在层流模式下(Re<150),采用极坐标系中NavierStokes方程的涡量-流函数形式,对置于弱电介质溶液中包覆电磁场激活板的圆柱绕流的尾迹和边界层的速度分布进行了数值研究。涡量传输方程采用ADI(隐式交替)格式,流函数方程采用FFT(快速傅里叶变换)格式,计算结果具有二阶空间精度和一阶时间精度,作了相应的实验,对实验结果进行了分析和讨论。     

用格子Boltzmann方法计算椭圆柱绕流的阻力

用格子Boltzmann方法模拟椭圆柱绕流, 研究椭圆柱形状对阻力的影响. 对圆柱绕流问题进行了数值模拟, 阻力系数的数值计算结果与相关文献数值相符. 计算了当Re=200, 椭圆柱纵轴长度不变、 横轴长度逐渐变大时几种不同形状的椭圆柱绕流, 并用插值方法处理了曲线边界, 用动量转换法计算了曲线边界受力. 计算得到了不同形状椭圆柱绕流的流线、 涡线以及阻力系数随横轴/纵轴长度比的变化趋势. 通过分析流线和涡线的变化, 给出了阻力变化的机理.     

圆柱绕流低维Galerkin方法的推广

把求解静止圆柱绕流问题的低维Ｇａｌｅｒｋｉｎ方法推广至可以求解绕横向或流向振动圆柱、旋转或旋转振荡圆柱等多种流动问题．方法的关键在于对以上四种流动分别选取不同的基本模态，使之满足相应的边界条件．对于振动或旋转振动圆柱绕流问题，它们与时间有关．于是扰动模态满足齐次边界条件，可以由此构造标准的Ｇａｌｅｒｋｉｎ算法．算例表明，用推广了的低维Ｇａｌｅｒｋｉｎ方法定性研究上述流动问题简便实用．     

聚焦波浪对直立圆柱作用的试验研究

为研究聚焦波浪对直立圆柱墩的作用情况,根据波的传播原理,采用JONSWAP谱生成不同参数下的聚焦波浪,用波峰波谷聚焦相加减的方法分析出聚焦波浪对直立圆柱墩的二阶绕射作用强于一阶波浪的绕射作用.波浪涌高数据统计分析表明:波浪涌高随波陡及结构物相对尺度的增大而增大;圆柱周围不同受浪角度的涌高也是变化的.