Numerical simulations of turbulent flows in aeroramp injector/gas-pilot flame scramjet

To uncover the internal flow characteristics in an ethylene-fueled aeroramp injector/gaspilot(ARI/G-P)flame scramjet,a Reynolds-averaged Navier-Stokes(RANS)solver is constructed under a hybrid polyhedral cell finite volume frame.The shear stress transport(SST)k-x model is used to predict the turbulence,while the Overmann’s compressibility corrected laminar flamelet model is adopted to simulate the turbulent combustion.Nonreactive computations for Case 1(G-P jet on),Case 2(ARI jets on),and Case 3(both ARI and G-P jets on)were conducted to analyze the mixing mechanism,while reactive Cases 4–7 at equivalent ratios of 0.380,0.278,0.199 and0.167 respectively were calculated to investigate the flame structure and combustion modes.The numerical results are compared well to those of the experiments.It is shown that the G-P jet plays significant role in both the fuel/air mixing and flame holding processes;the combustion for the four reactive cases takes place intensively in the regions downstream of the ARI/G-P unit;Cases 4 and 5are under subsonic combustion mode,whereas Cases 6 and 7 are mode transition critical and supersonic combustion cases,respectively;the mode transition equivalent ratio is approximately 0.20.     

Atmospheric response to a Hot SST Event in November 2006 as observed by the AIRS instrument

Using Atmospheric Infrared Sounder (AIRS) products of atmospheric temperature and geopotential height, we investigate the atmospheric response to HE0611, which was found and investigated by [Qin, H., Kawamura, H., Sakaida, F., Ando, K. A case study of the tropical Hot Event in November 2006 (HE0611) using a geostationary meteorological satellite and the TAO/TRITON mooring array. J. Geophys. Res. 113, C08045, doi: 10.1029/2007JC004640, 2008]. HE0611 was formed by connecting two very high SST areas, HE0611-East and HE0611-West. The period-mean atmosphere temperatures at levels of 925 and 850 hPa in HE0611-West are higher, by about 0.5 K, than those in WE0611-East while the atmospheric temperatures at middle to high levels (700–300 hPa) are higher in HE0611-East. The period-mean geopotential heights HE0611-East are much lower than those in HE0611-West for the levels from the surface to 400 hPa. The mean geopotential heights from 400 hPa to 200 hPa are higher in HE0611-East. In the middle and high layers over HE0611-West, the atmosphere temperatures gradually decrease from 7th to 17th, and then increase significantly. The increase in HE0611-East starts from 15th November, which is earlier than that of HE0611-West. The geopotential heights in the high layer of both the areas also show corresponding behaviors. The lagged atmospheric response in the western part is confirmed by the correlation analysis. It emerges that the atmospheric response to HE0611 is well organized and associated with deep convention in HE0611-East and subsidence in HE0611-West. These are also consistent with the HE0611 features and evolution revealed by earlier HE studies.     

Estimation of sea surface temperature, wind speed and water vapor with microwave radiometer data based on simulated annealing

A method for estimation of sea surface temperature, ocean wind speed and water vapor with microwave radiometer data based on simulated annealing is proposed. The proposed method shows about 60% improvement of sea surface temperature estimation accuracy in comparison to the existing method using Newton’s iterative algorithm.     

厚覆冰导线升力突变机理及数值模拟研究

针对新月形厚覆冰导线的升力系数在风攻角15°附近存在突变的问题,分别采用基于k-ωSST湍流模型的雷诺时均法和大涡模拟(LES)的数值方法对新月形厚覆冰导线在风攻角10°~20°范围进行了模拟。通过对比两种数值方法计算得到的覆冰导线气动力系数、流场结构和表面风压,发现LES方法能够更好地捕捉新月形覆冰导线表面的小尺度涡结构,得到的覆冰导线气动力参数计算结果与风洞试验数据高度吻合;而k-ωSST湍流模型难以模拟壁面上小尺度涡,捕捉不到升力系数的突变。根据覆冰导线不同壁面区域的压力分布,发现上侧壁面处的涡结构影响整体流场,并在下侧壁面曲率、来流夹角和壁面切线方向共同作用下导致升力系数突变。LES的气动力参数模拟结果可为覆冰导线防舞提供参考。     

Gurney Flap绕流气动特性研究

采用C-H型多块结构网格和SST k-ω湍流模型求解二维定常雷诺平均NS方程,对NACA4412单段翼型后缘附近和NACA63(2)-215B多段翼型襟翼后缘附近的Gurney Flap绕流进行了数值模拟。讨论了局部网格对计算结果的影响,重点研究了Gurney Flap位置和高度对气动特性的影响。计算结果表明:对于NACA4412单段翼型,Gurney Flap位置对气动性能的影响较小,但随着高度的增加,零升攻角和失速攻角减小,最大升力系数增大,同时阻力系数和力矩系数增加;对于NACA63(2)-215B多段翼型,Gurney Flap对气动性能的影响类似于NACA4412单段翼型情形。     

用转捩模型预测转捩及确定最佳粗糙带高度

采用SST k-ω二方程湍流模型,附加γ-Reθ转捩模型对不同迎角下的NACA0012翼型绕流进行了转捩预测,计算预测的转捩位置与实验结果符合较好,说明γ-Reθ转捩模型对转捩位置具有较好的预测能力。进一步用模型对带有粗糙带的翼型进行转捩预测,获得了不同雷诺数下刚好激发转捩又不引起过大局部压力扰动和附加阻力的最低粗糙带高度,对固定转捩风洞实验有指导意义。     

卫星星座整网轨道确定分析与仿真

编队飞行的卫星或卫星星座对轨道确定自主性和精度提出了较高要求。针对这个问题,通过建立星座的轨道动力学模型和星间观测的测量模型,将星座中的星间观测数据和地面观测数据融合起来,将待估的卫星轨道参数和部分动力学参数进行适当的分类,研究卫星星座整网轨道确定的新方法,并在理论上分析了整网定轨方法能提高定轨精度的原因;最后采用自主开发的卫星星座整网轨道确定软件进行了仿真计算。计算表明,该方法能有效地减少对地面站的依赖,并较大幅度提高定轨时卫星绝对位置和相对位置精度。
     

湍流模型源项的数值处理方法对计算结果的影响研究

涡粘性湍流模型广泛应用于各种工程湍流问题的计算。对于大多数湍流模型,在湍流控制方程的右端会出现源项,其刚性给数值计算带来很大的影响。从源项的物理意义出发分析了源项导致求解困难的原因,并在通用的求解湍流模型方程的点隐法基础上,以源项弱刚性的S-A湍流模型和源项强刚性的k-ωSST湍流模型为例,论述了生成项和耗散项之间的平衡关系对数值模拟的重要意义,给出了针对不同源项的具体处理方法。对RAE2882翼型跨声速流动算例的模拟结果表明,这些处理方法有效的提高了控制方程组的计算稳定性。     

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。