高超声速热流高精度数值模拟方法

通过高超声速下钝锥热流的计算,对影响热流计算精度的因素进行了综合研究.获得了壁面法向网格Re数对热流计算的影响规律.采用边界高阶插值,提高了热流计算精度.用Roe和AUSM(Advection Upstream Splitting Method)+格式搭配minmod和混合limiter两种限制器,探讨了空间离散方法及限制器的耗散性对热流结果的影响.应用五阶WENO(Weighted Essentially Non-Oscillatory)方法,对高阶格式计算热流的性能进行了细致研究,最终确立了网格依赖小,热流精度高的计算方法.研究认为:高阶格式不但可以放宽热流计算对网格Re数的限制,同时使得格式对本身计算方法的依赖性也变小.      

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM+ (AdvectionUpstreamSplittingMethod)格式 ,推广应用于二维非结构网格中 .利用格林 高斯公式对控制体内的变量进行线性重构 ,获得空间二阶精度 ,并采用Barth型限制器以抑制数值解的振荡 .为提高计算效率 ,采用了自适应多重网格法 .最后给出了NACA0 0 1 2翼型和平板激波反射的几个Euler方程的算例 ,并进行了分析比较 .     

高超声速气动热数值计算壁面网格准则

首先对高超声速气动热数值模拟中的3种壁面法向网格准则(平均自由程(MFP)、自由来流参数网格雷诺数和壁面参数网格雷诺数准则)进行了分析.随后,提出了壁面参数预估方法,使得壁面参数网格雷诺数和MFP两种准则无需进行试算即可直接预估壁面法向网格尺度.其次,运用提出的壁面参数预估方法结合不同网格准则,在同一实验条件下确定多个壁面法向网格尺度,通过与实验数据对比研究各网格准则确定的壁面法向网格尺度是否满足热环境模拟精度要求,并对提出的壁面参数预估方法进行验证.最后,开展了高超声速二维完全气体及真实气体效应算例、三维钝双锥算例的数值模拟.研究表明,壁面参数预估方法是可靠的;自由来流参数网格雷诺数准则随来流静温变化所确定的网格尺度变化趋势与物理推理相悖;基于壁面参数网格雷诺数和MFP两种准则主要受壁面温度影响,其在同一条件下所确定的网格尺度基本一致,并满足热环境模拟的精度要求.      

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中.利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡.为提高计算效率,采用了自适应多重网格法.最后给出了NACA0012翼型和平板激波反射的几个Euler方程的算例,并进行了分析比较.     

高超音速热流数值计算中的误差匹配原则

通过对粘性边界层的数值模拟和误差分析,提出了差分格式的误差匹配原则.在边界层内的计算网格一般都具有很大的长宽比,而计算误差与网格长宽比密切相关.为了使不同坐标方向上的误差保持在同一量级上,需根据误差匹配原则确定计算网格的尺寸.对绕钝锥的高超音速流动进行了数值模拟,进而计算了钝锥表面的热流分布.数值计算表明,参照误差匹配原则进行网格划分是有益的,可在保持流场数值计算精度和提高热流计算准确性的同时,使网格的数量和分布更为合理.     

限制器对高超声速体表面热流数值模拟的影响

采用Roe的FDS（Flux Difference Splitting）差分格式结合min mod,van Leer和Osher-C限制器,对高超声速体外部流场进行数值模拟.研究了限制器对高超声速体气动力/热数值模拟的区别,对比分析了3种限制器对三维模型表面热流数值模拟结果的影响.研究表明,限制器对数值模拟精度有重要影响,且对气动热数值模拟的影响大于对气动力的影响;对简单或较复杂的气动外形,采用Roe的FDS格式结合min mod限制器获得的表面热特性比它结合van Leer和Osher-C限制器时要好,可为热防护设计提供可靠结果.      

有限体积法中面积分离散格式的精度分析

讨论了应用顶点中心格式对欧拉方程进行有限体积离散时的精度,详细推导了顶点中心格式在平滑网格和扭曲网格上的截断误差,同时与中心格式和顶点格式进行了比较.格式精度分析采用泰勒级数展开的方法.分析结果表明:顶点中心格式和顶点格式类似,总能达到至少是一阶的精度,并且可以通过加密网格的办法减少误差,而且顶点中心格式更适合于分区计算.通过圆弧亚音速绕流和NACA0012翼型绕流流场的数值计算,验证了分析的结论.     

使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

快速大时间步长熵条件格式的分辨率研究

提出了高分辨率快速大时间步长熵条件格式的构造方法.用激波管问题对一族熵条件格式进行研究.在精度、步长、限制器方面进行了详细的数值实验,研究了同样计算量下各种格式的表现品质.从理论上保证了大时间步长格式的无振荡性质,从具体的数值实验分析中确定了大时间步长格式的分辨率问题.      

气动热CFD计算的格式效应研究

通过对N-S方程的求解计算了二维圆柱绕流及三维绕钝双锥高超声速流动过程中的气动热问题.对于求解方程,在空间离散的格式上,选用了3种上风格式——Roe的FDS,Van Leer的FVS和AUSM+格式以及一种中心差分格式,将各种格式所得的热流结果与实验结果进行了比较.对于AUSM+格式,还应用了minmod limiter 和 双 minmod limiter 两种限制器,并对结果进行了对比.      

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

绕AOTV高超音速流场的数值计算

本文采用实质上为二阶精度的隐式无波动、无自由参数耗散(NND)差分格式结合Steger-Warming矢通量分裂方法求解N-S方程,对AOTV外型的高超音速绕流进行了模拟。在计算中,对脱体激波采用了捕获法。作为验证,首先计算了二维圆柱绕流,并同文献[1]的数据进行了比较,取得了满意的结果,算例的计算结果还表明,采用本文的方法对激波的模拟具有较高的分辨能力。     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

新型三阶TVD限制器性能分析

在计算流体力学(CFD)方法中,限制技术是影响计算精度和计算稳定性的重要因素,目前应用较广的经典二阶总变差衰减(TVD)限制器虽能较好地满足计算要求,但性能差异大且分辨率和耗散的性能间并未得到良好权衡。对一种新型的三阶TVD插值限制器(T-3限制器)进行了研究并将其与3种经典限制器进行对比。首先通过一维黎曼问题,得出T-3限制器兼顾较高间断分辨率和良好稳定性的特点;接着通过高超声速双锥绕流和X-33外形飞行器的数值实验,得到T-3限制器具有刻画复杂流动的能力以及较优的气动热计算性能。     

非定常N－S方程的线化及其差分求解

在微幅简谐振荡的假设下，对计算平面内非定常Ｎ－Ｓ方程进行了线比，并在ＭａｃＣｏｒｍａｒｋ于１９８０年提出的半隐格式的基础上，发展了一种适用于求解计算平面内线化Ｎ－Ｓ方程的差分格式．用该差分格式求解计算平面内的线化Ｎ－Ｓ方程时，可以避免求解五对角方程；且时间步长可以取得较大，计算效率高．     

Newton-Kantorovich/Radau pseudospectral solution to rigid spacecraft time-optimal three-axis reorientation

The Newton-Kantorovich/Radau pseudospectral method is implemented in a homotopy scheme to solve the rigid spacecraft time-optimal three-axis reorientation problem. First, the reorientation differential equations associated with the current homotopic variable are linearized about the nominal (previous homotopic variable) solution. Then, the Radau pseudospectral method is used to solve the resulting linearized differential equations. By successively updating the nominal solution with the Newton-Kantorovich approach, the linearized solving procedure is iterated to achieve high-order precision. Because no numerical integration is needed and only linear algebraic equations are solved, the proposed method is computationally efficient. Finally, performance of the proposed method is illustrated through numerical examples.     

环境压强对气动塞式喷管底部压强的影响

为了了解环境压强对塞式喷管底部压强的影响,利用隐式NND格式(Non-oscillation, No free coefficients Dissipative scheme)求解二维层流N-S方程,计算了在没有底部二次流和二次流质量流率为1.4%的情况下,不同环境压强时的底部压强.结果表明,随着环境压强的升高,底部压强也在增加,但在闭尾迹时增加缓慢,在开尾迹时,增加迅速,并可能出现波动.      

超声速磁流体管道流动的数值模拟

应用可同时适用于描述磁流体MHD (MagnetoHydroDynamic)强化超燃冲压发动机磁流体发电器与磁流体加速器的,包含磁流体N-S(Navier-Stokes)方程与磁感应方程的可压缩磁流体方程组对超声速磁流体管道流动进行了数值模拟研究,以验证本项工作所提出的对2组方程进行交替迭代的分裂式算法.给出了算例的计算结果,并与无磁场及恒磁场下相同管道的流动计算结果进行了对比分析.