基于蒙特卡罗法辐射传递系数的统计分析

文章基于蒙特卡罗方法定义了表面热辐射换热辐射传递系数的表达式,确定其概率分布和置信度水平并进行统计分析。推导封闭腔内辐射传递系数的相对和绝对误差公式,为求解辐射传递系数所需追踪射线数量提供参考依据。     

发射特性对线性各向异性散射介质内瞬态温度分布的影响

考虑各向异性散射介质内的辐射传递方向，本文考察了辐射能的重新分配过程。利用射线踪迹法结合Hottel和Sarofim区域法推导了辐射传递系数(RTCs)。考虑一维各向异性参与性半透明灰介质。两边界表面为镜反射，一个表面半透明，另一个表面不透明。辐射传递系数用于计算辐射热源项。瞬态能量方程由全隐式控制容积法求解。外界辐射加热与对流冷却同时作用于边界表面。考虑四个线性散射相函数。以各向同性散射情况比较，研究不透明边界两侧发射特性和线性散射相函数对瞬态温度分布的影响。结果表明，各向异性散射与各向同性散射介质间的相对温度差异与散射不对称参数近似成正比。     

流线曲率法在多级跨声速轴流压气机特性预测中的应用

为研究多级跨声速压气机的分析问题,以通流理论为基础,采用了一系列适用于跨声速压气机的攻角、落后角和损失等经验模型,发展了一套基于流线曲率法的通流计算程序来预测跨声速压气机流场及其工作特性。为提高经验模型的预测精度,考虑到真实压气机中复杂的三维流动效应,针对部分早期模型进行了合理改进,包括改进了落后角模型使其适用于更大弯度范围叶型,以及采用一种更为合理的可变结构激波损失预测模型。针对两台跨声速压气机算例进行了计算校验,并将校验结果与实验值和三维数值计算进行对比。对比表明,设计工况下总压比最大计算误差为4.1%,效率误差为1.1%,在非设计工况特性预测和展向流场参数计算中也能得到和实验值相符的变化趋势,该通流计算方法可为现代跨声速轴流多级压气机特性分析提供具有参考价值的预测结果。     

进口水平投影可控的流线追踪内收缩进气道设计

为了满足两侧进气布局飞行器的乘波前体与进气道一体化设计要求,提出了一种进口水平投影可控的流线追踪内收缩进气道设计方法。基于马赫数分布可控的轴对称基准流场,在指定进口水平投影为椭圆的条件下,采用该方法设计了内收缩进气道并在设计点(Ma=5.4)和接力点(Ma=4.0)对其进行数值研究。结果表明,设计点时进气道都能保持基准流场的波系结构和沿程压力分布,无粘时可以全捕获自由来流,喉道性能与基准流场几乎相等。有粘条件下,设计点和接力点时进气道具有较高的压缩效率和良好的流量捕获能力,接力点的流量系数高达0.85。该设计方法为内收缩进气道与乘波前体的一体化设计提供了新途径。     

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.     

前后缘型线同时可控的乘波体设计

提出了一种前后缘型线同时可控的乘波体设计方法,在马赫数可控的外锥形曲面基准流场中,结合流线追踪技术和混合函数,实现了椭圆前缘转椭圆后缘的乘波体设计,并在设计点(Ma=6.0)和接力点(Ma=4.0)对其进行数值仿真研究。前后缘同时可控的乘波体在型面剧烈过渡处产生了较弱的激波,出口两侧存在高温高压区,后部对称面附近的激波形状由圆弧变为平直线且出口处流场基本均匀,非常有利于与进气道匹配设计。另外,该乘波体具有较高的容积率和预压缩效率,附面层修正后的容积率为0.24,设计点时乘波特性较好,接力点时前部完全乘波,具有较高的升阻比,有黏条件下设计点和接力点的升阻比分别为2.54和2.41。此外,与给定前缘的乘波体相比,其升力、阻力、俯仰力矩和出口增压比都有明显增加,但是升阻比和出口总压恢复系数有所降低,在设计点无黏升阻比由3.56降为3.00。以上研究表明,本文的设计方法可行且更加灵活,拓宽了乘波体的选择范围。     

航空紧固件生产车间流水线平衡量化分析及改善

为了量化地评价航空紧固件生产流水线的平衡性并对其进行改善,对一种典型的航空紧固件生产流水线进行了分析,结合流水线平衡的概念和航空紧固件生产流水线的特点,通过公式推导建立了一套可以从不同角度量化评价流水线平衡性的指标体系,阐述了其应用的方法和可行性;通过对比几种当前主要的研究流水线平衡的方法,指出航空紧固件生产流水线的平衡更适合采用工业工程的方法进行分析;研究了平衡性改进问题的解决思路,指出瓶颈工序对流水线的平衡性有重要影响,通过缩短瓶颈工序的加工时间可使整条流水线的平衡性得到有效改善,同时也使流水线上各工序的加工效率获得显著提升。最后通过实例计算对以上理论和方法进行了应用和验证。     

Numerical study of traveling ionosphere disturbances with vertical incidence data

High frequency ionosphere vertical sounding, as an important and representative application for detecting the ionosphere and studying the characteristics of radio propagation, can be utilized to monitor the ionosphere continuously variation and to acquire the ionosphere asymmetrical features of diverse scale above the ionosphere vertical sounding stations. This is a first article on real time application of numerical methods to obtain the parameters of traveling ionosphere disturbances (TIDs) using vertical incident ionograms. In this paper, the distribution of ionosphere electron density with TIDs is constructed using a background ionosphere model superimposed a perturbation theory model. The background ionosphere electron density is modelled by the inversion of vertical incident ionograms which are observed before the appearance of the disturbance. Based on the fourth order Adams-Bashforth-Moulton (the so-called ABM) predictor corrector method, instead of Runge-Kutta method, the fast digital ray tracing method is established. According to process of the disturbed trace simulation and parameters inversion, the characteristic parameters of ionosphere disturbance at different detection time can be obtained in real time. The numerical analysis of TIDs is then captured completely.     

基于射线追踪的快速雷达图像仿真方法

为了解决目标识别、隐身目标设计等研究领域对雷达图像需求日益增加的问题,提出了一种基于射线追踪的快速雷达图像仿真方法。该方法利用射线追踪计算每一个射线管的散射电场及其在雷达成像平面的投影坐标,通过将所有射线管的贡献在像域进行累加得到1幅二维雷达图像。通过对不同路径射线的聚类,可实现对散射贡献的分离,从而建立目标部件与强散射源的对应关系,并对复杂目标雷达图像进行解释。给出了Slicy目标的快速雷达图像仿真与算例分析,分析结果表明:基于射线追踪的快速雷达图像仿真方法可对复杂雷达图像进行有效预测,并能对强散射源进行成因分析。     

A comparative study of ionospheric IRIEup and ISP assimilative models during some intense and severe geomagnetic storms

Three-dimensional (3-D) electron density matrices, computed in the Mediterranean area by the IRI climatological model and IRIEup and ISP nowcasting models, during some intense and severe geomagnetic-ionospheric storms, were ingested by the ray tracing software tool IONORT, to synthesize quasi-vertical ionograms. IRIEup model was run in different operational modes: (1) assimilating validated autoscaled electron density profiles only from a limited area which, in our case, is the Mediterranean sector (IRIEup_re(V) mode); (2) assimilating electron density profiles from a larger region including several stations spread across Europe: (a) without taking care of validating the autoscaled data in the assimilation process (IRIEup(NV)); (b) validating carefully the autoscaled electron density profiles before their assimilation (IRIEup(V)).The comparative analysis was carried out comparing IRI, IRIEup_re(V), ISP, IRIEup(NV), and IRIEup(V) foF2 synthesized values, with corresponding foF2 measurements autoscaled by ARTIST, and then validated, at the truth sites of Roquetes (40.80°N, 0.50°E, Spain), San Vito (40.60°N, 17.80°E, Italy), Athens (38.00°N, 23.50°E, Greece), and Nicosia, (35.03°N, 33.16°E, Cyprus). The outcomes demonstrate that: (1) IRIEup_re(V), performs better than ISP in the western Mediterranean (around Roquetes); (2) ISP performs slightly better than IRIEup_re(V) in the central part of Mediterranean (around Athens and San Vito); (3) ISP performance is better than the IRIEup_re(V) one in the eastern Mediterranean (around Nicosia); (4) IRIEup(NV) performance is worse than the IRIEup(V) one; (5) in the central Mediterranean area, IRIEup(V) performance is better than the IRIEup_re(V) one, and it is practically the same for the western and eastern sectors.Concerning the overall performance, nowcasting models proved to be considerably more reliable than the climatological IRI model to represent the ionosphere behaviour during geomagnetic-ionospheric storm conditions; ISP and IRIEup(V) provided the best performance, but neither of them has clearly prevailed over the other one.