矢通量分裂-TVD杂交新格式及其应用

在贴体曲线坐标系中,本文提出了一个高分辨率的新格式,它是Steger-Warming的矢通量分裂与Harten的TVD(Total Variation Diminishing)格式的杂交,用于计算跨音速定常流动和捕获激波。典型的跨音速内流和叶栅流的算例表明:在60×15网格下捕获的激波过渡区不超过两个网格,证明其分辨率较高;分析激波附近的数值结果,没出现“低亏、过跳”、伪振荡现象;与实验结果的比较表明,对定常流在无人为附加耗散项的条件下它能给出较高质量无数值波动的激波流场解。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

用于二维无粘跨音速流计算的隐式TVD有限体积法(英文)

本文结合Harten的TVD(Total variation diminishing)概念和有限体公式建立一种求解守恒型Euler方程组的有效方法,采用通量分裂的迎风格式对方程组的隐式算子部分离散化,得到系数矩阵是主对角线占优的线性方程组,采用线Gauss-Seidel迭代推进法求解,提高了数值收敛速度,给出了二维跨音速内流的计算结果。     

三维湍流高速进气道内外流场的高效高分辨率解

本文提出一种高效、高分辨率求解三维复杂流场的隐式算法，并成功地用于求解高速进气道内外流场的三维、可压缩、雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，获得了三个工况五种流场的三维数值解；与有关实验数据比较，结果令人满意。这种算法由Ｊａｍｅｓｏｎ－Ｔｕｒｋｅｌ的ＬＵ分解法和张涵信ＮＮＤ格式的数值通量耦合而成，湍流采用ＢａｌｄｗｉｎＬｏｍａｘ模型。大量数值计算表明，该算法能捕捉到复杂的激波系，有较好的激波     

烧蚀外形三维粘性流场数值模拟

本文通过求解层流Navier-Stokes方程,数值模拟了三维奶头型烧蚀外形粘性流场。差分格式为Yee和Harten的迎风隐式TVD格式,在空间方向具有二阶精度。计算得到了高质量的激波和合理的流场结构。物面压力分布与实验结果进行比较,符合较好。     

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。     

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

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

自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

捕捉间断的高精度数值方法

为发展适用于捕捉超声速流场中各种间断的高精度算法,将通量限制的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的通量限制型差分格式.通过在时间方向上利用一阶精度格式计算的一维定常激波,以及在时间方向采用多步Runge-Kutta方法计算的一维非定常激波管问题上的数值试验与二阶精度的TVD格式所计算的结果比较,表明新方法比二阶精度方法在间断的捕捉上具有明显的优势.通过新方法的计算结果与精确解的比较,表明新方法的准度也是非常令人满意的.     

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

高速粘性内流的高分辨率高精度迎风型杂交格式

本文在有限体积离散和ＬＵ分解的基础上，构造出一个新的隐式迎风型杂交格式，并用于求解定常流动的稳态解。     

数值求解可压缩N—S方程的增量迭代法

本文在Beam-Warming格式的基础上,引入通量分裂的概念,应用行高斯迭代技术提出了一个二阶隐式的求解BFC坐标变换下的N-S方程的增量迭代法。它消除了因式分解带来的寄生误差,合理地处理了控制方程中的粘性项,具有较好的收敛性和稳定性。二维Poissuill流动和Couette流动的计算结果和精确解吻合良好。S形进气道的数值解反映出的现象和实验结果是一致的。     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

旋转坐标系N-S 方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明:所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

多级风扇三维粘性流场的数值模拟

 进行了多级跨音速风扇粘性流动的数值模拟研究。求解全三维N-S方程,通过高分辨率格式提高跨音速流场计算的准确性,利用LU-SGS隐式解法加快计算收敛速度。给出了一双级高负荷跨音速风扇设计和非设计工况下的详细计算结果,各工况下的计算与实验数据均吻合较好。     

PARALLEL COMPUTATION OF SUPER-SONIC BLUNT BODY VISCOUS FLOWFIELDS IN PVM

ＰＡＲＡＬＬＥＬＣＯＭＰＵＴＡＴＩＯＮＯＦＳＵＰＥＲ┐ＳＯＮＩＣＢＬＵＮＴＢＯＤＹＶＩＳＣＯＵＳＦＬＯＷＦＩＥＬＤＳＩＮＰＶＭＹａｎｇＸｉａｏｈｕｉ，ＷａｎｇＺｈｅｎｇｈｕａ，ＷａｎｇＣｈｅｎｇｙａｏ（Ｄｅｐｔ．１，ＮａｔｉｏｎａｌＵｎｉｖｅｒｓｉｔ...