双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

燃气自由射流的正格式数值模拟

采用二阶正格式方法对复杂燃气自由射流进行了数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对不同欠膨胀压力比下的燃气射流进行了数值计算。计算结果表明,该方法能够较好捕捉到波胞、滑移线、射流边界以及三波点等复杂射流流场波系结构,与实验照片反映的流动特征及高精度、高分辨率的三阶间断有限元方法的计算结果吻合较好,马赫盘位置与理论预估值和实验测量结果误差较小。证明该方法对激波具有较强的捕捉能力,为此类流动的数值研究提供了一种新手段。     

带副翼的翼身组合体绕流的 Euler方程解

提出了一种求解带副翼偏转的翼身组合体绕流的Euler 方程计算方法.将对接分区网格与分区求解算法相结合,有效地求解了绕此外形的复杂流动.提出了一种满足通量守恒的内边界耦合条件.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种限制器(Limiter)函数以保证TVD性质.数值算例表明本文方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

消除粘性项高阶离散数值振荡的半结点-结点交错方法

通常的粘性项处理方法是连续两次采用结点型中心差分格式求得.但是,通过两次采用高阶结点型中心格式求得的粘性项易在流场中的"间断"附近产生数值振荡.本文采用多种格式,通过求解一维Burgers方程来充分展示了这种振荡现象.为了消除这种数值振荡,一种半结点型的高阶紧致格式被用来求解粘性项.Fourier频谱分析表明,这种格式具有非常优越的保频谱性能.一维和二维的数值计算结果表明,通过"半结点-结点"交错采用这种半结点型格式的方法可以非常有效地避免粘性项高阶离散可能导致的数值振荡.     

深空组阵迭代FX相关合成算法及性能分析

针对火星距离下大规模天线组阵中极微弱接收信号的高效合成问题,在VLBI(Very Long Baseline Interferometry,甚长基线干涉测量)的FX相关器开环处理架构基础上提出了一种新的迭代FX相关合成算法,对算法原理进行了详细理论推导,对数据处理流程、算法合成性能等进行深入分析,与FSC(Full-Spectrum Combining,全频谱合成)、BC(Baseband Combining,基带合成)、SSC(Symbol-Stream Combining,符号流合成)等组阵方案进行性能仿真对比,并给出该算法在36接收天线组阵的实测试验结果。仿真和试验结果表明,迭代FX相关合成算法相比于FSC、BC、SSC等组阵方案具有较好的合成性能,在输入带内信噪比不低于-20d B时,合成效率达到了95%以上,非常适用于火星天线组阵微弱接收信号的组阵合成。     

数值模拟中的边界处理方法研究

为适应于模拟超声速流场中的一些特殊物理现象,研究了不同物理现象对边界处理的要求。针对激波、旋涡和声波等流场中典型的流动现象,采用常用的外推、特征和缓冲三类边界处理方法分别进行数值模拟。计算结果表明超声速流场边界采用外推技术,只适于计算带有激波的边界。特征方法还能适于旋涡运动的计算。声波的计算只能采用缓冲层技术。     

基于加速响应面逼近的改进CSSO优化过程研究

实现基于响应面的并行子空间优化过程(CSSO-RS),关键在于如何高效率高精度的构造近似响应面.本文针对如何快速提高响应面近似精度的问题,提出了方案数据库探测加点策略和限制记忆法构造响应面策略两种加速响应面逼近的措施.用经典测试函数和飞机总体设计优化两个算例研究表明,这两种策略能够有效加快响应面的逼近速度,使CSSO-RS过程的优化效率得到较大改善.     

一类新的二阶精度非交错形式差分格式的构造

本文对二维标量非线性双曲型守恒律，以非交错形式Ｌａｘ－Ｆｒｉｅｄｒｉｃｈｓ格式为基本模块，构造了一类新的不须解Ｒｉｅｍａｎｎ问题的差分格式，并证明其具有时空二阶精度。数值实验表明，格式具有较高的分辨激波、稀疏波及接触不连续的能力     

Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor

A robust unsteady rotor flowfield solver CLORNS code is established to predict the complex unsteady aerodynamic characteristics of rotor flowfield. In order to handle the difficult problem about grid generation around rotor with complex aerodynamic shape in this CFD code,a parameterized grid generated method is established, and the moving-embedded grids are constructed by several proposed universal methods. In this work, the unsteady Reynolds-Averaged Navier-Stokes(RANS) equations with Spalart-Allmaras are selected as the governing equations to predict the unsteady flowfield of helicopter rotor. The discretization of convective fluxes is accomplished by employing the second-order central difference scheme, third-order MUSCL-Roe scheme, and fifth-order WENO-Roe scheme. Aimed at simulating the unsteady aerodynamic characteristics of helicopter rotor, the dual-time scheme with implicit LU-SGS scheme is employed to accomplish the temporal discretization. In order to improve the computational efficiency of holecells and donor elements searching of the moving-embedded grid technology, the ‘‘disturbance diffraction method" and ‘‘minimum distance scheme of donor elements method" are established in this work. To improve the computational efficiency, Message Passing Interface(MPI) parallel method based on subdivision of grid, local preconditioning method and Full Approximation Storage(FAS) multi-grid method are combined in this code. By comparison of the numerical results simulated by CLORNS code with test data, it is illustrated that the present code could simulate the aerodynamic loads and aerodynamic noise characteristics of helicopter rotor accurately.     

高超音速双锥绕流数值模拟的格式效应分析

通过模拟双锥外形高超音速绕流,对各种计算流体力学空间离散格式的性能进行分析.研究了Roe格式、van Leer格式、AUSM(Advection Up-stream Splitting Method)格式、AUSMDV格式、AUSM+格式、上风TVD(Total Variation Diminishing)格式等对双锥高超音速流场特征的模拟能力.结果表明:各空间格式激波分辨率都较高,但粘性分辨率相差很大.其特性与格式耗散性有关,耗散性越小,粘性分辨率越高;格式耗散性对热流极值影响显著,耗散性越小,热流极值越大,其位置越靠近前缘.