Newton迭代法的P．C．格式

基于Newton迭代法 ,提出了一种有效的预估校正 (P .C .)迭代格式。本方法把牛顿法中的微分以差商替代 ,因而 ,既可求解具离散根的方程 ,也适用于有重根的方程求解 ,大量的试验结果表明 ,本文算法不仅可有效解决重根问题 ,而且有较高的收敛速度。     

用高效矢通量分裂格式求解二维叶栅的跨音流动

提出了一种高效隐式矢通量分裂格式。该格式通过在显格式上添加合适的高阶小量, 避免了近似因子分解和系数矩阵运算, 大大地减少了计算量。同时该格式在收敛性及激波捕捉精度等方面均保持了矢通量分裂格式所具有的优势。将该格式用到了二维叶栅的无粘跨音流场的数值模拟中去, 结果令人满意, 证实了无矩阵运算、无近似因子分解的隐格式构造有一定的推广意义。      

双曲型守恒律的无振荡中心差分格式

讨论了双曲型守恒律的一类无振荡中心差分格式。Ｈ．Ｎｅｓｙａｈｕ和Ｅ．Ｔａｄｍｏｒ研究了以交错型ＬａｘＦｒｉｅｄｒｉｃｈ格式（ＬｘＦ）为模块的无振荡中心格式的构造与熵不等式。此类格式利用高阶的ＭＵＳＣＬ型插值替代一阶分片常数逼近,减少了ＬａｘＦｒｉｅｄｒｉｃｈ格式的过多数值粘性,建立了一维标量非线性双曲型守恒律的一类高分辨格式。讨论以非交错ＬａｘＦｒｉｅｄｒｉｃｈ格式为模块建立起的差分格式。证明了此格式具有二阶精度、ＴＶＤ性质并在一定条件下满足熵条件。     

空气动力学数值模拟的一种分量形式TVD格式

本文选取一个精度较高的ＴＶＤ限制子,把著名的二阶Ｍａｃ Ｃｏｒｍａｃｋ格式和Ｗａｒｍｉｎｇ－Ｂｅａｍ格式结合,构造了一种分量形成的二阶ＴＶＤ差分格式。通过两个典型算例的计算,并与ＥＮＯ格式的结果比较,显示了该格式精度高和计算逻辑简单的特点。     

一类求解Hamilton-Jacobi方程的交错网格差分格式

Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程在控制论和微分对策中有广泛的应用,由于其表达形式与双曲守恒律方程极为相近,这有利于借助于求解双曲守恒律方程的差分格式来构造求解Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程的差分格式。文中将Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程变化为双曲守恒律方程,利用求解双曲守恒律方程的交错网格的Ｇａｕｓｓ型差分格式,构造了一类求解Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程的交错网格的Ｇａｕｓｓ型差分     

基于速度空间非结构网格和守恒修正的改进离散速度方法

传统离散速度方法在求解跨流域流动问题时,通常只求解动理学模型方程,即Boltzmann-BGK方程。与传统方法不同,改进离散速度方法同步求解了动理学模型方程和相应的宏观伴随方程。通过这种方式,可以将分子碰撞影响考虑到宏观伴随方程的通量计算中,同时宏观方程预估得到的结果可以用于预估平衡态,从而实现Boltzmann-BGK方程的全隐式离散。这两点改进可以有效克服传统方法在连续和近连续流区域计算效率低、计算精度差的缺陷。为了进一步减少速度空间的网格量和避免数值求积时的Runge现象,采用了非结构网格结合矩形律来离散速度空间并引进守恒修正来强制满足相容性条件。算例测试表明,采用速度空间非结构网格和守恒修正可以有效减少改进离散速度方法的计算量和内存花销。     

湍流方程扩散项离散对某径向叶轮流场计算的影响

在相对旋转坐标系下采用Harten,Lax and van Leer contact （HLLC） 格式离散对流项,自行开发了基于多块结构化网格的有限体积程序,实现了对叶轮机械内部流场的数值求解.分别对半开式径向叶轮、闭式后弯叶轮展开数值模拟,程序和商业软件计算得到的不同叶高处表面压力数据,其相对差异不超过1%,验证了算法的正确性.针对湍流方程的扩散项,分别使用完全离散和略去交叉导数项离散,通过对湍流黏度等值线、气动轴向力和力矩的比较表明:在网格正交性较好的情况下,略去交叉导数项的离散对计算结果的影响小于1%,显著地减小湍流方程离散的计算量.      

GPU平台上的叶轮机械CFD加速计算

通过数据并行的方式对一个成熟的叶轮机多块网格气动计算程序(MAP)进行了并行化处理,利用计算统一设备架构(CUDA)技术实现了在图形处理单元(GPU)上的并行计算.保留了原程序中的2阶空间迎风格式和隐式时间离散格式,并采用了隐式迭代对线性系统进行求解.经过2个叶轮机械算例的测试,与在传统的中央处理器(CPU)上运行的原程序相比,在计算结果完全一致的前提下,单GPU的计算速度最高可达单CPU计算速度的8.89倍,与四核并行的CPU计算相比可以得到2.39倍的加速.     

捕捉激波的群速度控制方法

回顾及综述了群速度控制方法基本原理、发展过程及最新进展。群速度控制方法是通过控制色散误差以捕捉间断的一种数值方法。该方法通过在间断两侧采用不同色散类型的格式,将数值波的传播压制在间断附近,避免数值振荡向周围扩散,从而达到抑制数值振荡的效果。该方法主体上采用数值色散抑制振荡,实际使用中仅需补充较低水平的数值耗散就可将振荡压制在可接受的范围内。因而群速度控制格式的总体耗散水平较低,易于实现小尺度波的高分辨率捕捉,是复杂流动高分辨率数值模拟的有效方法。     

一类中心型三阶格式及应用

通过分析格式的耗散关系和重线性化技术,发展相应格式构造技术,提出了一类中心型的三阶非线性格式,格式具有与NND格式相当的网格模板.由于格式的线性形式具有很好的耗散特性,同时重线性化技术的使用使格式在光滑区尽量采用线性形式,从而理论上提升了计算对旋涡等结构的分辨和保持能力;另一方面,格式的非线性机制使格式在非光滑区退化为类NND格式,实现耗散和稳定机制.验证计算表明,格式能够开展粘性和旋涡及脉动流场的计算应用,算法形式简单,便于工程应用,容易进行有限体积推广.