非结构网格上一种新型的高分辨率高精度无波动的有限元格式

在非结构网格上构造出无波动无自由参数耗散性有限元格式 ,即NND有限元格式 .通过若干个典型二维跨音速和超音速可压缩无粘定常流动的算例证明这确是一个高精度的 ,对激波具有高分辨率的无波动的新型有限元格式 .特别与网格自适应相结合 ,可得到十分满意的结果     

用于高速可压缩流体分析的带多维耗散格式的自适应Delaunay三角剖分

利用自适应Delaunay三角剖分并结合胞格中心迎风算法,分析非粘滞高速可压缩流体问题.推导了多维耗散格式,并采用非结构化三角网格的迎风算法,改善了激波的计算结果.解精度评价中引入误差估计,在网格重划分算法中,解梯度变化大的区域生成小单元格,解梯度变化小的区域使用大单元格.该格式能进一步推广到高阶时空的解精度分析中.通过稳态和不稳态的高速可压缩流体超音速激波和激波传播特性的分析,可以评估该算法的效率.     

一类TVD型的迎风紧致差分格式

给出一种迎风型TVD(total variation diminishing)格式的构造方法,该方法通过限制器来抑制线性紧致格式在模拟间断流场时的非物理波动,可构造出非线性TVD型紧致格式(CTVD)．然后采用该法构造出了3阶和5阶的TVD型紧致格式,并通过模拟一维组合波和Riemann问题,二维激波-涡相互干扰和激波-边界层相互作用等来考察它们的性能．数值实验表明了该类格式的高阶精度和分辨率,且过间断基本无振荡．     

计算激波的高精度数值方法

在分析了数值解在激波附近产生非物理振荡的原因后,构造了一个三阶迎风紧致格式以及激波的捕捉技术,并且,提出一种称为准装配法的新的激波装配方法.一维流动的数值试验表明,新方法是非常令人满意的.     

群速度直接控制四阶迎风紧致格式

为求解多尺度复杂流动在有限网格点的情况下高精度方法是可供选择的方法之一,其中紧致型格式数值解的分辨率更好.用迎风型紧致格式计算激波时,数值解中仍有数值振荡产生,这是由于对应于不同波数之群速度不均一所致.采用群速度直接控制方法重构紧致型格式以提高捕捉激波的能力.该方法简单,精度高,网格基架点少.算例表明,这一新的方法是可行的.     

守恒型双曲方程的中心差分TVD格式研究

研究了如何利用迎风格式的耗散性构造中心差分ＴＶＤ格式的方法，给 相应的定理，构造出新的耗散表达式。新格式既保留了二阶中心差分格式灵活方便的优点，又吸收了迎风格式耗散项比较精细的特点，同时具有ＴＶＤ性质，使得新格式具有较同的激波分辨率。     

基于AUSM分裂的二维通量分裂格式

基于对流迎风分裂思想构造的AUSM类格式具有简单、高效、分辨率高等优点,在计算流体力学中得到了广泛的应用.传统的AUSM类格式在计算界面数值通量时只考虑网格界面法向的波系,忽略了网格界面横向波系的影响.使用Liou-Steffen通量分裂方法将二维Euler方程的通量分裂成对流通量和压力通量,采用AUSM格式来分别计算对流数值通量和压力数值通量.通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的AUSM通量分裂格式.在计算一维算例时,该格式保留了精确捕捉激波和接触间断的优点.在计算二维算例时,该格式不仅具有更高的分辨率而且表现出更好的鲁棒性,可以消除强激波波后的不稳定现象.此外,在多维问题的数值模拟中,该格式大大地提高了稳定性CFL数,具有更高的计算效率.因此,它是一种精确、高效并且强鲁棒性的数值方法.     

基于有限体积法的非结构网格大涡模拟离散方法研究

非结构网格下的大涡模拟是解决复杂几何体高Reynolds(雷诺)数流动的有效途径.首先,基于有限体积法,研究了对流项和扩散项非结构网格下的离散方法.研究结果表明:基于TVD(total variation diminishing)限制器的限制中心差分格式保证了对流项的二阶精度并抑制了非物理振荡,同时,线性迎风格式虽然稳定,但数值耗散过大,且不能保证有界,中心差分格式引起了周期性非物理振荡;扩散项的超松弛非正交修正减小了网格非正交带来的离散误差,但修正系数须根据网格非正交的程度进行合理选取.为验证所述离散方法对大涡模拟的适用性,数值计算了Re=1.14×10~6下的非定常三维小球绕流,计算方法包括:计算网格用基于Delaunay三角剖分和Netgen前沿推进算法的四面体非结构网格;湍流模型用改进的延迟分离涡大涡模型;在离散格式的选取上,对流项用限制中心差分,扩散项加入非正交修正,插值格式用最小二乘法,时间项用二阶后向差分.计算结果表明,所用离散方法稳定收敛并且与实验数据基本吻合.     

激波捕捉差分方法研究

在迎风型格式和矢通量分裂技术的基础之上,对捕捉激波方法进行一种新的尝试．该方法首先对原始格式在特征方向上进行投影,然后用限制器对这些特征分量的变化幅值进行限制以抑止非物理波动,最后再把它转换成守恒形式,得到了基本上无振荡的激波捕捉格式．用该方法对两种迎风显示格式(二阶和三阶)和3种迎风紧致格式(三阶、五阶和七阶)进行处理,并在一维和二维的情况下进行了应用测试．通过与高阶WENO、MP、Compact-WENO等格式的比较,表明该方法在光滑捕捉激波的前提下仍有较高精度和分辨率．     

一种基于TV分裂的真正多维Riemann解法器

给出了一种真正多维的HLL Riemann解法器.采用TV(Toro-Vázquez)分裂将通量分裂成对流通量和压力通量,其中对流通量的计算采用类似于AUSM格式的迎风方法,压力通量的计算采用波速基于压力系统特征值的HLL格式,并将HLL格式耗散项中的密度差用压力差代替,来克服传统的HLL格式不能分辨接触间断的缺点.为了实现数值格式真正多维的特性,分别计算网格界面中点和角点上的数值通量,并且采用Simpson公式加权中点和角点上的数值通量来得到网格界面上的数值通量.采用基于SDWLS(solution dependent weighted least squares)梯度的线性重构来获得空间的二阶精度,时间离散采用二阶Runge-Kutta格式.数值实验表明,相比于传统的一维HLL格式,该文的真正多维HLL格式具有能够分辨接触间断,消除慢行激波波后振荡以及更大的时间步长等优点.并且,与其他能够分辨接触间断的格式（例如HLLC格式）不同的是,真正多维的HLL格式在计算二维问题时不会出现数值激波不稳定现象.