不同舵翼操纵面形式的水动力及艇后流场特性数值模拟研究

为研究共翼型舵和非共翼型舵两种操纵面形式的艇后流场特性及螺旋桨推进特性，基于SST(Menter)湍流模型建立了SUBOFF标准潜艇模型尾流场数值预报模型。经试验结果验证，所建立数值模型进行潜艇尾流场及螺旋桨推进性能预报具有较高精度。将SUBOFF潜艇模型的水平舵改进为共翼型舵及非共翼型舵，对艇后流场及水动力性能进行了预报。数值结果表明：在舵角小于10°时，共翼型舵使潜艇的俯仰力矩和垂向力相对非共翼型舵提升20%以上；在舵角超过10°时，共翼型舵的水动力优势随着舵角增大而减小。尾部流场预报显示：共翼型舵在小舵角时可以有效的消除舵翼结合处的涡流，同时共翼型舵可以有效的降低桨盘面伴流的不均匀性，对尾流品质的改善效果优于非共翼型舵。螺旋桨计算结果显示：共翼型舵螺旋桨推进系数在大部分舵角下都小于非共翼型舵，在舵角为20°时，共翼型舵相对非共翼型舵推力系数下降3.5%，扭矩系数下降2.4%；同时共翼型舵的桨盘面流场均匀度要优于非共翼型舵，舵角为5°时，共翼型舵桨盘面处流场不均匀度相对于非共翼型舵要降低7.1%，舵角为25°时则降低25.1%。     

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

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

Preliminary study on heat flux measurement data of TT-0 flight test

With the explosive development of aerospace science, the design of the new generation airliner at higher speeds is attracting more attentions. To achieve this goal, it is necessary to achieve accurate prediction of the aerodynamic heating / force loads and successful reduction of drag and heat flux. As a remedy for the existing studies which are based upon the CFD and wind tunnel tests, this study presents a flight test for the drag and heat reduction spike technology. The principal goals of this flight test were to provide reference for verifying the accuracy of the prediction technology on ground and promote the development of the drag and heat reduction technology. By adopting the OS-X rocket, the TT-0 test vehicle designed by Shenyang Aircraft Design & Research Institute reached a maximum Mach number of 5.8 and a maximum altitude of 38 km. Hypersonic and supersonic pressure data by pressure scanning valves and heat fluxes by gauges at different locations were obtained successfully. Also, heat fluxes obtained by in-house CFD code are illustrated in comparison with the flight data. The results indicate that the numerical errors are large in most cases. More technologies, such as more CFD codes and more numerical procedures, should be adopted to conduct studies on this issue in the future.     

Effect of scavenge port angles on flow distribution and performance of swirl-loop scavenging in 2-stroke aircraft diesel engine

Swirl-Loop Scavenging (SLS) improves the performance of 2-stroke aircraft diesel engine because the involved swirl may not only benefit the scavenging process, but also facilitate the fuel atomization and combustion. The arrangement of scavenge port angles greatly influences in-cylinder flow distribution and swirl intensity, as well as the performance of the SLS engine. However, the mechanism of the effect and visualization experiment are rarely mentioned in the literature. To further investigate the SLS, Particle Image Velocimetry (PIV) experiment and Computational Fluid Dynamics (CFD) simulation are adopted to obtain its swirl distribution characteristics, and the effect of port angles on scavenging performance is discussed based on engine fired cycle simulation. The results illustrate that Reynolds Stress Turbulence model is accurate enough for in-cylinder flow simulation. Tangential and axial velocity distribution of the flow, as well as the scavenging performance, are mainly determined by geometric scavenge port angles α_geom and β_geom. For reinforcement of scavenging on cross-sections and meridian planes, α_geom value of 27° and β_geom value of 60° are preferred, under which the scavenging efficiency reaches as high as 73.7%. Excessive swirl intensity has a negative effect on SLS performance, which should be controlled to a proper extent.     

Computational engineering analysis of external geometrical modifications on MQ-1 unmanned combat aerial vehicle

This paper focuses on the effects of external geometrical modifications on the aerodynamic characteristics of the MQ-1 predator Unmanned Combat Aerial Vehicle (UCAV) using computational fluid dynamics. The investigations are performed for 16 flight conditions at an altitude of 7.6 km and at a constant speed of 56.32 m/s. Two models are analysed, namely the baseline model and the model with external geometrical modifications installed on it. Both the models are investigated for various angles of attack from −4° to 16°, angles of bank from 0° to 6° and angles of yaw from 0° to 4°. Due to the unavailability of any experimental (wind tunnel or flight test) data for this UCAV in the literature, a thorough verification of calculations process is presented to demonstrate confidence level in the numerical simulations. The analysis quantifies the loss of lift and increase in drag for the modified version of the MQ-1 predator UCAV along with the identification of stall conditions. Local improvement (in drag) of up to 96% has been obtained by relocating external modifications, whereas global drag force reduction of roughly 0.5% is observed. The effects of external geometrical modifications on the control surfaces indicate the blanking phenomenon and reduction in forces on the control surfaces that can reduce the aerodynamic performance of the UCAV.     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

火星大气来流模拟装置CFD仿真与试验

根据我国火星着陆巡视器工作过程,其着陆发动机需要在相对火星大气高速迎风运动中可靠点火。由于巡视器着陆时发动机喷管出口气流与火星稀薄气流方向相反,目前无法通过理论计算准确获得着陆过程的动态流场对发动机起动过程的影响量值。为验证火星着陆环境下发动机点火的适应性,需要建立发动机的火星大气来流试验环境模拟条件。为模拟发动机在火星大气条件下的相对运动,在真空舱内发动机保持固定,前端设置环形来流形成装置,该装置在发动机喷管周围形成一定速度的逆向来流包络。采用数值模拟技术结合试验验证方法,在火星着陆器巡视器主发动机性能考核试验中,针对来流的形成装置开展了设计研究工作。来流模拟试验测试数据表明:在确保贮箱供应压力稳定的条件下,来流模拟系统能够形成100~200 m/s速度的稳定来流,发动机在来流下能稳定启动工作,真空舱压力满足试验要求。     

针对旋翼流场模拟的自适应网格算法

针对基于格心格式求解器的旋翼流场模拟,提出了相应的自适应笛卡尔网格的数据存储结构及自适应算法。给出了相应的单元处理策略,简化了对自适应笛卡尔网格的处理;对于频繁的自适应加密过程中产生的大量重复点,采用了高效的交替数字树算法(Alternating digital tree,ADT)予以删除;对于自适应疏化过程中产生的大量无用点,提出了标记-删除-移动(Mark,delete,move,MDM)算法予以快速地删除,减少了不必要的计算资源消耗。对CaradonnaTung旋翼在不同悬停状态下进行了模拟验证,对比了压力分布系数与桨尖涡位置。之后对HELISHAPE 7A旋翼在前飞情况下进行了模拟验证,计算值与实验值吻合。此外,求解器对桨尖涡的捕捉效果得到了明显的提高,表明本文方法具有良好的有效性与鲁棒性。     

风-雨耦合环境超高层三塔连体建筑雨压分布研究

以中国东南沿海某超高层三塔连体建筑为研究对象，以风-雨双向耦合算法为核心，基于计算流体动力学技术采用连续相和离散相模型进行风场和雨场的迭代模拟。首先基于9种风雨组合工况进行三塔连体建筑非定常脉动风场模拟，探讨超高层三塔连体建筑平均风压分布、表面速度流线和流场干扰机理。然后对比研究不同风雨组合工况下主塔表面雨滴附着数量、雨滴冲击力和雨压系数的分布规律，揭示风-雨耦合场中结构表面速度流线、雨滴运行轨迹和最终速度的作用机理。最后提炼出超高层三塔连体建筑最不利风-雨组合工况，并给出对应的雨压系数取值建议。研究表明：风-雨耦合环境下超高层连体建筑迎风面雨荷载作用最为显著，此时雨荷载与风荷载最大比值可达23.81%，局部测点最大雨压系数达到0.301，100 a重现期风速和强大暴雨组合为风-雨耦合作用的最不利组合工况。     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。