前缘对进气道性能影响的数值模拟

应用质量平均的N-S方程和B-L代数湍流模型，在不同压缩面前缘半径条件下，以高超声速进气道的流场进行了数值模拟，对比研究了进气道压前缘半径对流场结构及性能的影响，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散。     

计算跨音速N─S方程组的二阶隐式矢通量分裂格式

提出了一种二阶精度隐式矢通量分裂单步差分格式并用于求解时间平均Ｎ－Ｓ方程组。与传统方法相比,该方法是Ｊａｃｏｂｉ谱半径格式和无近似因子分解。数值试验表明：该方法具有精度高、计算量小和收敛快等优点。在平面叶栅和扩压器进气道的跨音速流场的求解中,较好地捕获了激波、分离等现象,与实验数据比较结果令人满意     

二维超／高超声速进气道流场数值模拟

对超／高超声速三级压缩进气道流场进行了数值模拟,来流马赫数为4、6,进气道内流动为层流状态,根据二维Navier-Stokes方程,采用二阶精度Roe格式进行离散。按照流场特点,合理地设计网格分布及调整不同黏性范围的熵修正,防止了壁面附近过大的数值耗散,使计算结果更加合理。在进气道模型的各级压缩折转角处,获得了清晰的激波结构,在进气道内部的各种波系的相交、反射和激波诱导的边界层分离等现象都得到合理的描述。计算得到的压力分布,在各级压缩斜板上同简单波理论结果十分接近。用本文方法计算了另一个二级压缩进气道,沿上、下壁面的压力分布与试验比较符合得较好。     

基于激波拟合与柔性壁的反设计数值方法研究

将激波和待求激波压缩面作为运动间断边界拟合,研究了超声速来流下的型面反设计数值计算方法.算法秉承拟合法思想,采用二阶格式时间推进欧拉方程组,在运动边界和动网格条件下,按预估的波系结构划分网格块,反求可产生指定压力分布的三维型面,并获得精确的激波型面.通过包括凸包进气道在内的三个反设计算例验证了本文算法的合理性及其对网格质量的适应性,证明了在型面反设计中,采用激波拟合法可大大提高型面精度.本算法可用于超声速进气道压缩面和唇口完全三维反设计.     

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

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

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

高超声速飞行器加速度导数数值模拟

在Etkin的非定常气动力模型基础上,采用非定常Navier-Stokes方程描述飞行器的非定常强迫沉浮振动流场,采用空间二阶精度的NND格式和"双时间"时间离散法进行数值求解,研究飞行器加速度导数的数值计算方法。并给出了返回舱外形和HBS(Hyper Ballistic Shape)标模外形的加速度导数、旋转导数和俯仰阻尼导数的计算结果及分析。     

非定常流动计算的混合时间推进方法研究

结合ADI方法的二阶时间精度特性和高斯-赛德尔迭代法的良好收敛性特性,提出了一种ADI-高斯-赛德尔混合迭代格式,应用该格式,对NACA0012翼型的跨声速俯仰振荡进行了EULER方程求解,空间离散格式采用ROE’S FDS格式,计算结果与实验及成熟的LUSGS双时间步迭代方法计算结果进行了分析比较。计算结果表明,该格式精度能达到LUSGS双时间步迭代方法完全收敛的二阶精度,且该格式精度对内迭代次数的依赖性较低。分析表明,本文提出的方法具有较高的精度和较好的收敛性。     

多重网格加速的LUSGS算法用于喷管流场数值模拟

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用独特的显式算法克服了两方程湍流模型刚性对流场计算收敛速度的影响,并给出了该方法的理论依据。通过轴对称喷管内流场、二元喷管内外流场的数值模拟验证了采用以上方法在计算效率上的提高。最后将RNG模型计算结果与SA模型计算结果及实验数据做对比分析。     

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

高精度两步TVD格式的构造及数值检验

基于通量分裂的思想,利用Taylor级数理论和第二步限制器函数的控制,构造了在光滑区域空间为三阶、时间为一阶精度的高分辨率全变差递减(TVD)格式.该格式在激波过渡点降为一阶迎风格式.并从单个线形方程推广到非线性方程及方程组情形.通过几个典型算例的计算,并与二阶TVD格式作了比较,表明该方法对流场中的激波有较高的分辨率,且是无波动的.     

更准确、更精确、更高效——高超声速流动数值模拟研究进展

从准度、精度和效率3方面回顾了近几十年来高超声速流动数值模拟研究的进展。在物理模型方面,介绍了高超声速数值模拟中高温气体效应、稀薄气体效应以及湍流效应的建模与模拟,基于雷诺平均Navier-Stokes(RANS)方程重点对现阶段较为关注的高超声速边界层转捩的模式理论研究进行了介绍。在空间离散算法方面,主要介绍了高超声速数值模拟中常用的二阶精度迎风格式以及高阶精度格式的发展及其应用。在时间推进方面,主要回顾了隐式时间推进方法的发展及其应用。在误差和不确定度估计方面,主要介绍了其概念、来源以及常用的分析方法,同时给出了迭代误差估计、Richardson外插法以及敏感性导数方法等初步研究结果。最后,讨论了高超声速流动数值模拟中下一步需关注的问题。     

WNND格式在稳态问题中的应用

基于二阶NND格式,引入Jiang和Shu的空间模板加权思想,构造了一种空间守恒变量型三阶精度的Weighted NND格式(WNND格式),通过组合波、球头以及Stardust返回舱绕流等多个算例对WNND格式进行了验证,数值结果表明,WNND格式引起的耗散和波动较小,并且能够高精度地分辨间断.在返回舱数值模拟中,引入了ADI隐式时间离散方法,加快了收敛速度,提高了计算效率.     

应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

某小型涡扇发动机全流道准三维数值法

采用高阶精度Godunov格式和时间相关法,用多部件模型在小型涡扇发动机的全流道S2流面上求解带黏性力项的非定常Euler方程组,并将整机模拟的结果与试验结果进行了对比,在合理给出损失的情况下,计算结果与试验结果基本吻合.     

有限体积时间推进法在跨音速扩压器湍流流场中的应用

本文提出了一种有限体积时间推进法,将其应用于均化N-S方程的求解,模拟跨音扩压器的湍流流场,并采用空间变时间步长和多重网格技术加快收敛速度。计算结果与有关文献的实验数据吻合良好。     

Magnetic Reconnection Phenomena In Interplanetary Space

Interplanetary magnetic reconnection(IMR) phenomena are explored based on the observational data with various time resolutions from Helios, IMP-8, ISEE3, Wind, etc. We discover that the observational evidence of the magnetic reconnection may be found in the various solar wind structures, such as at the boundary of magnetic cloud, near the current sheet, and small-scale turbulence structures, etc. We have developed a third order accuracy upwind compact difference scheme to numerically study the magnetic reconnection phenomena with high-magnetic Reynolds number (R _M=2000–10000) in interplanetary space. The simulated results show that the magnetic reconnection process could occur under the typical interplanetary conditions. These obtained magnetic reconnection processes own basic characteristics of the high R _M reconnection in interplanetary space, including multiple X-line reconnection, vortex velocity structures, filament current systems, splitting, collapse of plasma bulk, merging and evolving of magnetic islands, and lifetime in the range from minutes to hours, etc. These results could be helpful for further understanding the interplanetary basic physical processes. This revised version was published online in August 2006 with corrections to the Cover Date.