高速旋转冲压增程弹用进气道复杂流场数值模拟

采用块结构网格与二阶精度流场分区求解技术,对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟,得到了高速旋转工况下对应于不同来流马赫数和攻角,以及临界工况时超音速进气道内外流场复杂的波系结构。在旋转工况下,进气道流场结构与不旋转时相似,但旋转速度以及攻角的存在对冲压发动机进气道的总体性能产生了负面影响,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大;冲压发动机进气道在较低马赫数工作时的综合性能优于在较高马赫数工作时的综合性能。     

数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

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

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

转速对高速旋转冲压增程弹用进气道影响数值模拟

结合分区对接网格技术以及二阶精度区域分解算法，对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟。得到了高速旋转工况下对应于不同来流攻角和旋转角速度，临界工况时，超声速进气道内外流场复杂的波系结构。随着旋转角速度的提高，进气道总压恢复系数和动能效率均有所降低，而流场畸变指数则显著增大。特别是当转速达到20kr／min，进气道总压恢复系数和动能效率下降趋势以及流场畸变指数增大趋势更明显。攻角的存在对冲压发动机进气道的总体性能产生了负面影响。     

高超声速三维碳—碳烧浊流场的数值研究

本文针对高超声速再体的烧蚀现象 ,利用简单隐式TVD差分格式和激波捕捉法 ,数值求解三维化学非平衡Navier Stokes方程 ,其中化学模型是碳 碳 (C C)空气化学模型 ,考虑 12个化学组分和 31个化学反应过程 ,研究了C C烧蚀对再入体头部区域的壁面温度和热流分布的影响。为了计算效率和稳定性提出壁面条件显式处理的方法。对再入高度为 6 5km和速度为 8km/s的再入体头部区域烧蚀流场进行了数值模拟 ,用飞行迎角α =0°的计算结果与国外文献进行了比较 ,符合得较好。同时给出了三维小迎角α =5°的计算结果     

TVD格式数值模拟斜腔内旋涡运动

 采用ＴＶＤ有限体积格式,通过求解非定常可压缩Ｎ－Ｓ方程,对恒定出入口压力比驱动下的斜腔内部流动进行了数值模拟。计算后得到清晰的流场涡系结构及其发展过程；显示了空腔流动中的自维持振荡现象；并分析了剪切层不稳定对周期性振荡的作用机制。     

预估校正的高分辨率格式及其应用

本文从基于通量分裂的一阶迎风格式出发,构造了一类校正型具有高分辨瘁的二阶全变差递减(TVD)格式,并用此格式进行了一维激波管问题及二维问题的计算。结果表明格式的解为二阶精度,激波处无明显的振荡,激波过渡区为1～2个点。     

考虑粘性时激波反射及壁面热流计算

本文采用显式 TVD 格式,求解薄层 Navier-Stokes 方程,计算了激波与楔面相互干扰的流场情况,得出了马赫反射的各种等强度线及固壁面上详细的物理量分布,根据计算得出的流场,计算了壁面摩擦阻力及热流分布。     

A Class of TVD Type Combined Numerical Scheme for MHD Equations With a Survey About Numerical Methods in Solar Wind Simulations

It has been believed that three-dimensional, numerical, magnetohydrodynamic (MHD) modelling must play a crucial role in a seamless forecasting system. This system refers to space weather originating on the sun; propagation of disturbances through the solar wind and interplanetary magnetic field (IMF), and thence, transmission into the magnetosphere, ionosphere, and thermosphere. This role comes as no surprise to numerical modelers that participate in the numerical modelling of atmospheric environments as well as the meteorological conditions at Earth. Space scientists have paid great attention to operational numerical space weather prediction models. To this purpose practical progress has been made in the past years. Here first is reviewed the progress of the numerical methods in solar wind modelling. Then, based on our discussion, a new numerical scheme of total variation diminishing (TVD) type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account convergence, stability and resolution. This new MHD model is established by solving the fluid equations of MHD system with a modified Lax-Friedrichs scheme and the magnetic induction equations with MacCormack II scheme for the purpose of developing a combined scheme of quick convergence as well as of TVD property. To verify the validation of the scheme, the propagation of one-dimensional MHD fast and slow shock problem is discussed with the numerical results conforming to the existing results obtained by the piece-wise parabolic method (PPM). Finally, some conclusions are made. This revised version was published online in August 2006 with corrections to the Cover Date.     

自身场磁等离子体推力器数值仿真流场特性分析

磁等离子体推力器以其推力大、比冲高等特点,成为未来深空探测、星际航行任务首选的电推力器类型,研究推力器内部等离子体流场特性,有利于解释推力器出现的物理现象。针对特定自身场磁等离子体推力器,建立磁流体模型,使用TVD Lax-Friedrich格式以及ADI方法对推力器内部及羽流进行数值求解,得到等离子体流场及等离子体参数分布情况,仿真结果显示在高电流工况下,等离子体羽流更加集中,轴向加速效果更加显著,但高电流模式下阴极温度较高,不利于阴极的使用寿命。