基于模拟退火算法与响应面模型的三维气动优化设计方法

提出了一套适用于高耗时三维气动设计问题的优化设计体系.其主要思想是采用改进拉丁超立方体试验设计选取样本点,应用自行开发的三维粘性流场求解程序进行流场计算建立数据库,采用二次响应面方法建立近似模型,再应用高效模拟退火算法进行全局寻优.以NASA rotor57为对象,在详细进行流场计算基础上采用所提出的优化体系对其进行了三维积叠优化设计.在对流量、效率加以严格约束的条件下,总压比可提高1.8%,对流场结构进行了分析.优化结果表明本优化方法省时,适于三维气动设计的特点.     

亚跨超CFD软件平台研制

亚跨ＣＦＤ软件平台研制的宗旨是综合国内外计算流体动力学的最新研究成果，用计算机软件技术的新发展，建立一个面向对象的计算流体动力学的研究开发环境，加快计算流体动力学最新成果向工程应用的转换速度。本文介绍初级版本的研究进展和研究成果，主要匝数值网格生成技术，流场解算器，面向对象的界面设计等三部分组成。     

求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

A hybrid lattice Boltzmann flux solver for integrated hypersonic fluid-thermal-structural analysis

In this paper, a hybrid Lattice Boltzmann Flux Solver (LBFS) with an improved switch function is proposed for simulation of integrated hypersonic fluid-thermal-structural problems. In the solver, the macroscopic Navier–Stokes equations and structural heat transfer equation are discretized by the finite volume method, and the numerical fluxes at the cell interface are reconstructed by the local solution of the Boltzmann equation. To compute the numerical fluxes, two equilibrium distribution functions are introduced. One is the D1Q4 discrete velocity model for calculating the inviscid flux across the cell interface of Navier–Stokes equations, and the other is the D2Q4 model for evaluating the flux of structural energy equation. In this work, a new dual thermal resistance model is proposed to calculate the thermal properties at the fluid–solid interface. The accuracy and stability of the present hybrid solver are validated by simulating several numerical examples, including the fluid-thermal-structural problem of cylindrical leading edge. Numerical results show that the present solver can accurately predict the thermal properties of hypersonic fluid-thermal-structural problems and has the great potential for solving fluid-thermal-structural problems of long-endurance high-speed vehicles.     

基于OpenFOAM的超临界压力下低温射流大涡模拟

为了深入认识低温射流在超临界压力下的流动特性,将真实流体热物性模型写入开源CFD程序OpenFOAM平台中,开发出了用于模拟超临界压力下流动过程的均相求解器。针对液氮跨临界、超临界射流进行大涡模拟,研究了环境条件变化对于低温射流流动特性的影响。结果表明：采用的真实流体模型能够在广泛的温度与压力范围内准确计算流体的热物理性质,在此基础上开发的均相求解器能够准确描述超临界条件下的低温射流,其流动特性主要表现为稠密液体与环境气体之间的剪切层不稳定以及湍流混合,此类变密度射流在湍流充分发展区域同常密度射流一样具有自相似特性;超临界射流相比于跨临界射流具有更高的混合效率以及更短的液核穿透长度;环境温度升高以及环境压力降低均使得低温射流与周边气体之间的密度梯度增加,对剪切层中不稳定波动的发展起到抑制作用,导致了射流液核长度的增加以及扩张角度的减小。     

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed, the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows. However, it needs numerous logical operations in supersonic cases, which limit the method’s applicability in engineering problems greatly. In order to overcome this defect, a hybrid multidimensional Riemann solver, called HMTHS (Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures), is proposed. It simulates the strongly interacting zone by adopting the MHLLES (Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures) scheme at subsonic speeds, which is with a high resolution by considering the second moment in the similarity variables. Also, it adopts the MULTV (Multidimensional Toro and Vasquez) scheme, which is with a high resolution in capturing discontinuities, to simulate the flux at supersonic speeds. Systematic numerical experiments, including both one-dimensional cases and two-dimensional cases, are conducted. One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves, shock waves, and linear contact discontinuities. Two-dimensional cases, such as the double Mach reflection case, the supersonic shock / boundary layer interaction case, the hypersonic flow over the cylinder case, and the hypersonic viscous flow over the double-ellipsoid case, indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows. Therefore, it is promising to be widely applied in both scholar and engineering areas.     

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

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

High-order implicit discontinuous Galerkin schemes for unsteady compressible Navier–Stokes equations

Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.     

二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

基于可满足性问题求解器的星上FPGA永久损伤容错技术研究

现代卫星广泛使用的FPGA在空间高能粒子的影响下,会产生门电路的永久性损伤。而传统的三模冗余等容错方法不但成倍增加了系统硬件开销,还存在因冗余器件耗尽而失效的风险。因此,提出一种利用FPGA自身冗余资源,修复永久性损伤的容错方案。该方案通过建立FPGA内部资源的功能模型,将容错问题转化为数学上的可满足性问题。并且利用经过改进的GSAT算法对该问题求解,可以获得在功能上与损伤前完全相同的电路结构,及其所对应的FPGA配置文件。将该文件重新下载到FPGA中,可以屏蔽损伤带来的影响,从而达到利用FPGA自身冗余资源容错的目的。通过实验和分析可以看出,本文方案具有对损伤修复成功率高、计算量小和需要内存空间少的特点,因此符合星上计算能力和硬件资源十分有限的实际情况。