用于有限体积法的中心WENO混合格式

通过理论推导得到6阶中心格式与5阶WENO格式相结合的混合格式(CWENO6格式),为了提高有限体积法中高精度格式的计算效率,引入3阶MUSCL反卷积方法以取代传统有限体积法中的反卷积方法。通过不同的算例比较了CWENO6格式与WENO5格式的数值表现。结果表明,相比于传统的WENO5格式,在网格量相同的前提下,由于数值耗散较小,CWENO6格式对流动细节的描述更为准确,为有限体积法模拟三维多尺度复杂提供了一种更为有效的方法。     

WENO格式权重分析与改进

为了优化WENO格式计算性能,在对Jiang和Shu的经典WENO格式(记为WENO-JS)加权方法分析的基础上,通过引入间接光滑指数,构造出一种新的WENO格式——WENO-E格式,取得减小间断区耗散的效果。理论分析表明,该格式与WENO-JS格式计算效率基本相同,可达到相同阶的计算精度;但在相同网格下,较之WENO-JS格式,该格式对光滑区域的求解有更小的截断误差,对间断的捕捉有更高的分辨率。与WENO-JS格式相比,采用WENO-E格式进行线性迁移方程、非线性Burgers方程、欧拉方程等相关问题的数值实验,均能取得更好的数值结果。     

采用重组模板的权重优化WENO-Z格式

针对精确模拟含激波等复杂流动结构的流场对高精度格式的低耗散低色散要求,基于5阶有限差分WENO-Z格式,提出一种模板重组技术。在计算WENO非线性权时,引入一个由3点模板重新组合的4点模板,优化原格式中各模板的权重分配,进而提出了两种改进WENO-Z格式。采用近似色散关系分析方法对改进前后格式色散与耗散特性进行了对比与分析。分析表明：两种改进格式耗散有不同程度的降低。数值实验表明:改进格式具有更优越的激波捕捉性能,对小尺度流场结构具有更高的分辨率。     

用高阶高精度 WENO 格式求解二维激波附面层干扰流场

采用时间相关法求解二维 Navier-Stokes方程 ,数值模拟二维平板层流附面层与激波干扰流场 ,给出了物面压力分布和应力分布。计算中 ,对流项空间导数的差分离散采用高阶高精度 WENO格式 ,时间方向采用具有 TVD性质的 Runge-Kutta方法 ,粘性项采用二阶中心差分。所得压力分布和应力分布与国外实验结果吻合较好 ,计算实践表明高阶 WENO格式具有优异的性能 ,应用前景广阔     

高精度高分辨率WENO格式分析与改进

高精度、高分辨率计算格式对提高超声速复杂流场的计算水平有重要的意义。WENO格式是在ENO格式基础上发展起来的一类新型高精度、高分辨率格式,本文在格式加权算法理论和数值分析的基础,对加权算法进行了改进,通过模型问题应用可以看出:改进后格式权值的计算更为光滑,仍具有高精度、高分辨率的性质,并克服了随网格变化而出现的波动现象。     

高精度有限体积格式在三维曲线坐标系下的应用

为了构造在光滑区具有较高分辨率并且可以无振荡捕捉激波的高精度有限体积格式,同时降低格式在模板选择上遇到的困难,基于逐维重构方法,发展了结构网格下的高精度有限体积格式,并将这一格式推广到三维曲线坐标系下,从而可以适应相对复杂外形下的计算。为充分验证格式的有效性,选取一系列典型算例进行计算:在等熵涡输运和二维喷管流动中验证了格式的精度可以达到设计精度;在双马赫反射问题中格式也表现出良好的捕捉激波的能力。数值计算表明,上述格式在曲线网格上具有较高的数值精度和鲁棒的激波捕捉能力,适用于流体力学方程的计算。     

一种基于低扩散通量分裂方法的WENO格式研究

将低扩散通量分裂格式(LDFSS)和加权基本无振荡格式(WENO)相结合,构造出一种混合格式,其中WENO格式用于物理量重构,而LDFSS用于通量分裂。采用这种格式对Riemann问题、钝头体高超声速无粘绕流流场进行了计算,并对超声速平板湍流边界层进行了混合LES/RANS模拟,计算结果表明:相对于采用Lax-Fridrichs分裂的WENO格式来说,这种混合格式对于激波和接触间断的分辨率更高,并且在标量保正性方面更优,收敛性更好;而相对于采用带有Minmod限制器的MUSCL方法进行物理量高阶重构的LDFSS格式来说,这种混合格式在混合模拟的计算中能够更好地反映湍流流场的脉动特性,计算得到的湍流速度脉动量的统计值更加准确。     

高精度有限体积法与间断有限元法的比较研究

本文通过数值算例,比较了高精度有限体积法以及间断有限元法在求解不同问题时表现。研究发现：在精度相同的条件下,间断有限元法的计算误差要明显小于有限体积法;间断有限元法的重构过程与高精度有限体积法相比较为简单,但高阶情形下解多项式的自由度较多并且需要计算体积分,因此整个求解时间较长;在同样精度条件下,高精度有限体积法的计算效率要明显高于间断有限元法。降低时间积分时解多项式的自由度数目是实现高精度算法在实际问题中应用的重要手段。     

高精度有限体积法与间断有限元法的比较

通过数值算例,比较了高精度有限体积法和间断有限元法在求解不同问题时的表现。研究发现:在精度相同的条件下,间断有限元法的计算误差要明显小于有限体积法;间断有限元法的重构过程与高精度有限体积法相比较为简单,但高阶情形下解多项式的自由度较多并且需要计算体积分,因此整个求解时间较长。降低时间积分时解多项式的自由度数目是实现高精度算法在实际问题中应用的重要手段。     

测试投稿，大家不要管哈~

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

超声速飞行器流场的并行数值计算及效率分析

基于有限体积方法、TVD差分格式和显式Runge-Kutta迭代方法的框架,针对超声速/高超声速飞行器绕流流场,在超级并行计算机上完成了2~64个CPU并行数值计算工作。通过测试程序在超级计算机上的并行效率,并将并行程序应用于航天飞机绕流流场计算,检验了计算程序进行大规模并行计算的性能。结果表明,在负载平衡的条件下,程序在该超级并行计算机上达到了不同程度的超线性加速比,并行效率最高达到了126%,远远高于微机Cluster并行平台上的结果,适合复杂流场的大规模并行计算。     

一个二阶精度的大粒子有限体积方法

对二维Euler方程给出一个二阶精度的非结构网格有限体积方法,方法的主要思想是把Euler方程的一个时间步分两步计算:第一步只考虑压力加速效应,第二步再考虑输运效应,离散过程中采用一种具有最小二乘思想的线性重构函数以计算交界面的流通量,数值实验表明该方法是非常有效的。     

高超声速横向喷流的三维数值模拟

采用有限体积法,结合半隐的Harten TVD格式和Baldwin Lomax代数湍流模型,求解了三维全N S方程,计算了超声速气流干扰流场,并与实验数据进行了比较,模拟清晰地捕捉到喷流附近复杂波系结构和旋涡结构,并且显示出了弹体表面受到喷流影响的区域,说明横向喷流干扰对弹体影响非常大。     

高超声速乘波飞行器三维流场的并行数值模拟

针对高超声速乘波飞行器三维绕流流场,在基于LINUX+MPI系统的分布式并行计算平台上,并行求解了三维雷诺平均的N-S方程。并行数值方法采用的是有限体积方法(FVM)、OC-TVD差分格式、B-L代数湍流模型及流场分区的并行方法。计算结果表明,所采用的并行数值模拟方法能够求解包含强激波的流场,激波穿越区域边界时无断层、错位等通量不守恒的现象。并行计算效率高,8个处理机计算时的并行加速比达到了6 8。     

等截面隔离段中激波串结构的数值模拟

采用二阶TVD格式的有限体积法偶合Baldwin Lomax代数湍流模型求解雷诺平均Navier stokes方程 ,数值模拟了二维矩形等截面隔离段中的流动现象 ,很好地模拟出由激波 /附面层干扰所形成的复杂的激波串流场结构。计算结果与国外有关文献的结果进行了比较。     

非定常Euler方程数值计算中高精度格式应用

采用五阶精度加权紧致非线性格式(WCNS)和非定常“双时间步”方法求解非定常Euler方程,模拟NACA0012翼型强迫俯仰振动流场,研究了高精度格式应用到非定常计算时“双时间步”方法物理时间步长、子迭代收敛判据、子迭代步数以及物理时间导数离散方法对计算精度和计算效率的影响.     

激波诱导惰性堆积粉尘运动的建模及数值模拟

基于空间平均的双流体模型 ,引入描述颗粒脉动速度的粒化温度 ,建立颗粒相的本构方程 ,将粒化温度模型推广到解决高速流动问题。采用AUSM+ 有限差分法 ,数值模拟激波在一定厚度的惰性粉尘床中传播及诱导粉尘颗粒运动的过程。结果表明 ,粒化温度模型较好地描述了这一过程 ,证实了此模型在处理稠密颗粒高速流动和流场存在强间断中的可行性