Fast spacecraft solar radiation pressure modeling by ray tracing on graphics processing unit

A novel method is presented to evaluate on the graphics processing unit (GPU) the force and torque on a spacecraft due to solar radiation pressure. The method employs efficient ray tracing techniques, developed in the graphics rendering discipline, to resolve spacecraft self-shadowing and reflections at faster than real-time computation speed. The primary algorithmic components of the ray tracing process which contribute to the method’s computational efficiency are described. These components include two-level bounding volume hierarchy acceleration data structures, fast ray to bounding box intersection testing using the slab intersection algorithm and fast triangle intersection testing using the Möller-Trumbore algorithm. Spacecraft material optical properties are represented as a combination of Lambertian diffuse and ideal specular reflections. Both diffuse and specular ray-surface interactions are modeled. The approach is implemented using C++ and OpenCL and executed on a consumer grade GPU. Model validation is presented comparing ray traced force and torque values to the same quantities produce by a faceted analytic model. Numerical results illustrate the impact of self-shadowing on the force and torque calculation, and demonstrate the fast computational speed that is enabled with this implementation.     

Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.     

基于ALE法的水陆两栖飞机船体着水性能分析

水陆两栖飞机的着水性能是体现其综合性能的重要因素之一.应用ALE法对三种水陆两栖飞机船体截面进行入水仿真计算,从着水载荷、喷溅性能及运动响应三个方面对着水性能进行综合评估.对比分析了直线型、带舭弯弧形、复合型三种构型截面在不同速度、不同重量下的压力分布、喷溅性能以及运动响应,得到了三种截面横向压力分布规律.计算结果显示相同状态下,喷溅性能复合型优于带舭弯弧形、带舭弯弧形优于直线型;过载复合型远大于带舭弯弧形,带舭弯弧形稍大于直线型,综合看来带舭弯弧形截面着水性能最佳.     

长板件数控加工技术的探索与研究

结合某型机板件在数控加工中的研制生产状况和存在的问题,论述了实现各种长板件类结构零件装夹定位、刀具及其参数、加工方案及程编策略.通过这些技术措施的实施,在一定程度上解决了在加工长板类零件装夹和定位时出现的零件变形、让刀效率低、质量不稳定等问题.     

动力失效对起飞纵向气动特性影响数值研究

研究民用飞机动力失效对飞机起飞阶段纵向气动特性的影响规律及其机理,对保证飞机有效操纵和安全飞行具有重要意义.采用在点对点多块结构化网格系统上求解三维可压缩雷诺平均N-S方程的数值方法,研究发动机动力失效对某民用飞机起飞构型纵向气动特性的影响.通过DLR-F11模型验证研究方法对民用飞机高升力构型气动特性的预测能力;针对安装动力短舱的某翼吊涡扇发动机民用飞机起飞构型,通过对比其发动机在正常工况和失效时飞机气动特性的差异,得出动力失效对飞机纵向气动特性的影响规律及机理.结果表明:动力失效后,不但溢流效应会使飞机阻力系数增大;还会导致发动机进、排气特性较正常工作状态明显不同,恶化短舱附近流场,对飞机的升力、失速特性带来不利影响.     

典型工况下低排放燃烧室的压力振荡特性

为了研究低排放燃烧室在典型工况下的压力振荡特性,针对模型燃烧室进行了燃烧自激振荡特性试验.在试验中测量了采用贫油预混预蒸发(LPP)燃烧技术的低排放燃烧室在典型工况下的压力振荡频率和幅值,在燃烧室进口压力为1.10～2.77MPa、燃烧室进口温度为656～845 K、燃烧室压降为3.41％～4.35％范围内,分析了燃油粒径变化对振荡特性的影响.分析结果表明:局部当量比脉动是引发燃烧不稳定的因素之一.通过计算燃油二次雾化状态下的液滴最大粒径,发现燃油液滴粒径的变化对主燃级出口处的局部当量比脉动有直接影响,从而引起燃烧室压力振荡幅值和频率的变化.     

叶片反扭对跨声速大涵道比风扇性能的影响

使用基于流固耦合算法的叶片反扭程序,考虑了非定常气动力对叶片变形的非线性作用,研究了叶片反扭对跨声速大涵道比风扇性能的影响.以冷态叶型为起点,先计算离心力作用下的叶片变形,在此基础上使用流固耦合程序获得非定常气动力作用下的变形.考察了3个转速下叶片的动态变形对大涵道比风扇气动性能的影响.结果表明:在跨声速工况下,叶片表面激波位置的变化对叶片反扭角有很大作用,在考察的转速范围内,堵塞点使用设计叶型计算的流量大于动态叶型下的流量,数值达7%,将导致发动机起飞推力小于预测值.结果表明在大涵道比风扇设计阶段,预测气动性能使用准确叶型的重要性.      

不同进口截面下液力透平非定常压力脉动计算

为了研究不同进口截面对液力透平内压力脉动的影响,利用CFX软件对一单级液力透平进行非定常数值计算,通过设置监测点,计算各监测点在不同进口截面下的压力脉动,通过快速傅里叶变换将其压力脉动计算结果做相应分析,分析各监测点处压力脉动的时域和频域分布.结果表明:蜗壳内在大蜗壳进口下,随着蜗壳进口直径的增加,每个监测点处的压力逐渐增加,而在小蜗壳进口下,每个监测点处的压力逐渐减小.在距离蜗壳收缩管较远处,大蜗壳进口下的压力脉动较小,而在距离收缩管较近处,小蜗壳进口下的压力脉动较小.在叶轮内不同进口截面压力脉动的差异在同一时刻从进口到出口逐渐减小.尾水管内在小蜗壳进口下尾水管进口处的压力脉动最大.      

翅片位置对复合式减涡器减阻性能影响数值模拟

为了探索翅片-管复合式减涡器的翅片安装位置对共转盘腔径向内流压力损失的影响规律,对不同转速、翅片周向位置及安装角度下的去旋系统开展了数值研究,得到了不同工况下共转盘腔径向内流的流场结构及压力损失分布曲线。研究结果表明:减涡管能引导流体径向流入,并降低流体的旋流比;相比于管式减涡器,翅片-管复合式减涡器能明显降低盘腔内的总压损失;在不同旋转雷诺数下,翅片的周向安装位置α及安装角β均存在最佳值;在中、高旋转雷诺数下,最佳值分别为α=9°,β=30°,最佳结构下总压损失较基础模型低40%左右;改变翅片周向位置及安装角度可以明显改变气流进入减涡管的角度,在较优情况下,可以减小流体流入减涡管的阻力及在减涡管内的流动阻力,整体上减小了盘腔内总压损失。     

Al/Ni含能多层膜在Al2O3陶瓷/不锈钢焊接中的应用

采用磁控溅射设备,生长AuSn合金做焊料层、Al/Ni含能多层膜做热量提供层,实现了不锈钢和Al_2O_3间的异质材料自蔓延高温扩散焊。利用SEM、XRD和DSC等测试手段表征AuSn合金和Al/Ni含能多层膜的微观形貌、相成分和放热量;用万能试验机测试焊接接头的力学性能。结果表明,AuSn合金的质量比基本达到80∶20,而多层膜的层状结构清晰,反应热达到1 239 J/g。焊接实验结果表明,仅使用AuSn焊料时,剪切强度仅为46 MPa,在增加Al/Ni含能多层膜后,其剪切强度可达90 MPa,强度提高了约一倍。焊接接头的界面显微形貌和相结构研究表明,剪切强度的增强主要是Al/Ni多层膜提供了额外能量使得界面处的反应剧烈,陶瓷金属化层与中间层的反应加剧,形成了新的反应生成物。