AUSM格式在二维浅水方程求解中的应用

引入空气动力学中发展起来的AUSM格式,在三角形网格上采用有限体积法,对二维浅水控制方程进行数值求解;引入过程中,对浅水问题中的源项及露滩问题等采用了相应的特殊处理。通过对二维计算常用的溃坝算例及扎龙湿地实际算例的计算结果分析可知:二维浅水问题中引入该计算格式后,计算结果较好且保持了原格式的优良特性,同时也验证了该格式在二维浅水方程中应用的可行性及可靠性。

基于Van Leer+AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer+AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明:发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。

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

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

AUSM+格式的改进

为了提高AUSM+格式对流场计算中激波和剪切流的分辨率,减小数值振荡和激波后数值过冲问题,本文在原有AUSM+格式的基础上经过数值研究分析,对AUSM+格式进行了声速、特征马赫数的处理方式,对网格界面马赫数等进行了改进,并用改进后的AUSM+格式求解了含有脱体激波、斜激波、剪切流、膨胀扇、及湍流边界层等流动特征的复杂流动。经典型算例检验证实,改进后的AUSM+格式对激波和剪切流的分辨率得到了加强,求解的准确性和收敛性有一定的提高。改进后的AUSM+格式继承了原AUSM+格式优点,计算处理问题能力增强,可以有效地应用于亚声、跨声和高超声流动的工程实际问题中。

基于Van Leer+AUSM混合格式超声速流场的并行数值算法研究

基于CLUSTER集群并行体系结构与消息传递库MPI并行环境,开展了超声速流场的并行数值算法研究。计算网格采用多块对接粘性结构网格,控制方程为积分形式的Navier-Stokes(N-S)方程,空间离散采用基于Van Leer+AUSM混合格式的有限体积法,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,应用一方程S-A模型计算湍流的粘性。并行计算采用分区并行的思想,设计了满足负载平衡要求的区域分解和通讯方法,编程上采用单控制流多数据流(SPMD)模型以及Master/Slave并行模式。最后,通过对某弹丸和某制导火箭弹超声速流动的数值模拟,验证了并行算法的准确性、高效性和处理复杂几何外形的能力。

AUSM类格式在热流数值模拟中的比较研究

系统分析了AUSM类空间离散格式中AUSM+、AUSMDV及AUSMPW三种格式的构造方法及特点,并分别对钝锥、钝双锥的高超粘性绕流流场进行了数值模拟,重点分析了格式对热流数值模拟结果的影响。数值结果表明:AUSM类格式在高超音速流动模拟中具有高间断分辨率和数值稳定性,其中AUSM+和AUSMPW格式更适用于热流数值模拟,而AUSMDV格式可能会出现"红玉"现象。

AUSM+格式在亚音、跨音和超音速流动数值模拟中的应用

在任意曲线坐标系下,采用有限体积法计算可压缩流动。利用AUSM+格式替换传统的Rhie-Chow动量插值方法,给出有限体积界面的速度和压力求解方式。为了适用于亚音速流动,通过预处理的方法给出了速度-压力耦合的表达式。算例考核结果表明,该算法在不可压流动到超音速的全范围内都有良好的适用性,而AUSM+格式的计算结果均优于Rhie-Chow动量插值方法。

上风格式的若干性能分析

采用理论分析和数值试验相结合的方法,对上风格式的三个典型代表:FDS、FVS及AUSM格式的性能进行了分析比较:FDS格式具有高的间断分辨率但可能出现非物理解;FVS格式捕捉强非线性波的能力很强,计算效率高,但耗散比较大;AUSM格式兼有高间断分辨率和高计算效率且克服了两者的弱点。以一维无粘激波管流场为算例考察了这些格式性能与间断分辨率。

一种基于AUSM思想的通量分裂方法

根据对流迎风分裂(AUSM)思想提出一种通量分裂方法,称为K-CUSP格式.它与传统H-CUSP和E-CUSP格式的最大差异在于总能量的分裂:K-CUSP格式将无粘守恒通量中所有的运动学量分裂到对流项,所有的热力学量分裂到压力项,即总能量被分裂成动能和静焓.对于压力项的数值通量,采用一种新的界面构造方法.数值测试表明:1K-CUSP格式继承了FVS格式的简单性和稳健性.在激波后不易出现压力过冲,在膨胀区域没有振荡,优于AUSM和WPS格式;2K-CUSP格式继承了FDS格式的分辨率.激波间断的分辨率和H-CUSP、E-CUSP格式基本相同,接触间断的分辨率高于FVS格式,低于Roe、AUSM和WPS格式.AUSM和WPS格式在计算运动接触间断时,速度存在很大振荡,而新格式不存在振荡.

广义拉格朗日乘子MHD方程组的两类逆风分裂格式

开发高性能的磁流体力学数值模拟方法是提高空间天气数值预报研究的一个重要方面.有限体积法的逆风分裂格式具有良好的间断捕获能力,Steger-Warming和AUSM(Advection Upstream Splitting Method)是逆风分裂格式FVS(Flux Vector Splitting)方法中具有代表性的两种格式.采用这两种格式求解具有伽利略不变性的扩展型广义拉格朗日乘子磁流体力学(EGLM-MHD)方程组,对Orszag-Tang涡流问题和三维爆炸波问题进行数值模拟,结果表明两种格式均能得到稳定精确的数值结果.与StegerWarming格式相比,AUSM格式产生的磁场散度误差更小,计算速度更快.

几种典型迎风格式的分析与比较

重点分析了AUSM+格式的构造方法,同时将其与Van Leer,Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM+格式,Van-Leer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率。获取了对于基础研究和工程应用很有意义的结论。

高超声速稀薄流场AUSM分裂式DSMC-IP方法研究

针对直接模拟蒙特卡洛(DSMC)方法统计耗散较大而传统的信息保存(DSMC-IP)方法难以有效模拟强激波的问题,采用对流迎风分裂(AUSM)通量计算格式对IP方法进行改造。以局部马赫数为标准重构控制方程中的关联项通量,使计算更加准确的符合激波两侧的流动特征,从而形成一种具有较高统计精度和高超声速流动模拟能力的新型DSMC-IP方法。采用该方法对超声速圆柱和带扩张角的喷管高超声速流动进行数值模拟,结果显示,该方法与DSMC方法的计算结果相比,其流场结构基本相符,而前者的统计耗散明显更低,且流动特征更清晰,表面特征系数与实验值或参考值相差在5.5%以内,证明了该方法在高超声速稀薄流场模拟中的准确性和有效性。

先进上风格式跨声速CFD计算研究

介绍了空气动力学跨音速计算中应用广泛的几种上风格式,同时引入了JAMSON中心差分格式作为比较,针对激波分辨率,格式精度,计算效率等方面的问题对ONERAM6机翼进行了对比性计算,讨论了ROE的FDS格式,AUSM+格式和AUSMDV格式的计算性能,得出了很有意义的结论。

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。

稠密气固两相流的流场实验显示和粒化温度模型

在水平激波管中 ,研究激波与堆积粉尘的相互作用 ,并通过带示踪粒子的X射线和阴影摄影 ,记录堆积粉尘内外流场的波系结构 .引入“粒化温度” ,基于稠密气体分子动力学的思路导出堆积粉尘的守恒方程和本构方程 ,利用AUSM+ 格式对激波与堆积粉尘相互作用现象进行了数值模拟 .

叶栅激波流场的数值模拟

给出了ＡＵＳＭ 格式在任意曲线坐标下的形式，将ＡＵＳＭ 格式和三阶精度ＭＵＳＣＬ格式融合，并结合ＬＵ－ＳＧＳ隐式时间推进法，对典型双圆弧叶栅无粘跨音速和超音速流动进行了数值模拟．

Spalart-Allmaras湍流模型在高超声速气动加热计算中的应用

将AUSM+格式与LU SGS隐式迭代相结合,采用一方程的Spalart Allmaras湍流模型和二阶迎风MUSCL格式,研究了火星"探路者号"实验模型在10马赫来流条件下的气动热问题。采用量热完全和热完全两种气体模型的计算结果与实验数据进行了比较。热完全气体模型模拟的热流分布规律与实验吻合较好,但是对分离区的热流估计过高。在网格合理划分,尤其是满足近壁面的网格分辨率要求的前提下,该方法能够合理反映高超声速流场除分离区以外的气动热特性。

激波碰撞干扰流动非定常效应的数值研究

将AUSM +上风格式应用于数值模拟半球绕流激波碰撞干扰的流动中 ,揭示了该流动在定常来流条件下的非定常本质 .分析弓形激波后超音速“喷流”的流动结构及其非定常效应的形成机制 ,提出一种新的激波碰撞干扰流动的非定常形成机理 .

后效期火药气体流空过程数值模拟

后效期火药气体流空规律与火药气体作用系数历来是各国兵器研究的重要内容。现有的火药气体作用系数的计算方法多是在准定常假设下得到的理论或经验关系式,适应面较窄,误差普遍偏高。本文采用AUSM+格式和双时间步法,对65式37mm高炮(无膛口装置和带不同锥形消焰器)共5个算例的火药气体流空过程进行了数值模拟,根据数值结果计算了火药气体作用系数,与实验值相比,误差能控制在6%以内,符合较好。为进一步研究复杂条件下的后效期流空规律打下了基础。

近空间高超声速弹箭气动数值模拟

针对近空间高超声速飞行弹箭气动力、气动热问题,采用Spalart-Allmaras(SA)湍流模型和热完全气体模型,运用时间相关有限体积方法及AUSM+格式对三维非定常可压缩流体Navier-Stokes方程进行了数值模拟。再现了超高速条件下弹箭周围的复杂流动现象,得到热流分布规律,为进一步近空间超高速弹箭研究奠定了基础。计算结果表明阻力系数随着攻角的变化显著变大;弹箭尾翼气动加热很严重,弹箭尾翼变大时气动加热量明显增大,在设计弹箭时应引起足够重视。