1控制方程.doc

CSTAM2014-B01-0250 混合网格上的并行高阶间断有限元 黏流数值方法研究 秦望龙* ，吕宏强 ，伍贻兆 南京航空航天大学航空宇航学院空气动力学系 第八届全国流体力学学术会议 2014年9月18—21日 甘肃·兰州 Copyright ? 2012 版权所有 中国力学学会地址: 北京市北四环西路15号　邮政编码:100190 　Address: No.15 Beisihuanxi Road, Beijing 10间断Galerkin有限元（discontinous Galerkin, DG）方法由于精度高，模板小，易于实现网格和阶数自适应。前高阶间断有限元方法求解流体力学问题存在几个难点：复杂流动模型求解，激波捕捉和计算效率。相比于有限体积方法，间断有限元方法要求解更多的单元自由度，且随着阶数的提高，单元自由度增加。二维方程，有限体积方法单元自由度为4，间断有限元方法自由度为，为，为，为60这使得一些计算在串行计算机上。单一的网格单元常常不能满足效率和精度的双重要求流场区域和特征。 单流多数据流single program multiple data, SPMD）编程模式和信息传递界面（Message Passing Interface, MPI）数据交换方法[13]，文章有效地求解了Navier-Stokes(N-S)方程和方程，。本文的网格分区采用METIS黏性流动的求解一直是计算流体力学中的热点。升阻力预测，尤其是阻力预测一直是计算流体力学领域的一大难题。为了提高CFD阻力预测的可信度，增强工业界对CFD结果的信心，美国航空航天学会（AIAA）专门成立了阻力预测研讨会（DPW）[15-17]，我国也启动了航空高可信CFD开放式专题研究（HiCFD）活动，同时考虑升力和阻力预测问题。本文采用高阶间断有限元方法对黏流问题进行数值求解，考量该方法对于阻力的预测结果1.1 Navier-Stokes方程 守恒形式可压缩的Navier-Stokes方程 （1.1） 其中，是有界域。守恒变量，无粘数值通量，黏性数值通量可以写成如下矢量形式 ，,（1.2）式中，，，，分别是密度，压强和单位总能和单位总焓，是笛卡尔坐标系下的速度分量。是黏性应力分量，是热通量分量。若不考虑 ，则式（1.1）退化为守恒形式的欧拉方程。 (1.3) (1.4) (1.5) (1.6) 是定压比热容与定容比热容的比值，文中取1.4。是层流普朗特数，空气取0.72，是流体的温度。是动力黏性系数，可以通过Sutherland公式得出。 (1.7) 是温度等于288.15K时空气的粘度= 288.15K，=110.4K 1.2 Spalart-Allmaras一方程模型 文章采用修正的Spalart-Allmaras（S-A） （1.8）来替换传统S-A模型[18]方程中黏性通量和源项中的项，使得变量为负数时原方程退化为简单的对流方程，提高了方程稳定性[19-21]。 （1.9） （1.10） （1.11）是湍流涡粘系数 （1.12） (1.13) (1.14) (1.15) 其中是旋转张量，d是积分点到物面的实际距离。公式中常数的选取见文献[18] 式（）分段函数保证了为正数时,湍流生成项是正数,提高了数值求解的鲁棒性[20]。 2.高阶间断有限元数值离散 方程（1.1）两边同乘测试函数V离散，有 是对流通量和粘性通量的和。将计算域表示成网格单元的集合，对单个单元则有 （2.2）是单元，是单元的边界。在单元内，和有如下高阶表达形式 （2.3） （2.4） 其中，是阶数p对应的个基函数。将式(2.3),(2.4)代入公式（2.2），得： （2.）BR2方法[23](the second Bassi_Rebay method)进行求解，即在一阶导数的求解中额外添加单元交界面间断引起的梯度贡献 （2.6） 其中，和分别是面上左单元和右单元的变量值。是单元的一阶导数额外项，对于面的一阶导数项，有 （2.7） 可以看出， （2.8）来代替式（2.5）项 （2.9） 通量函数包括无粘数值通量与粘性数值通量函数。文中无粘数值通量采用LLF数值通量函数[24],粘性数值通