射线跟踪法的改进及其应用

研究复杂组合体目标后向散射场问题，并对射线跟踪法的典型算法进行改进，以提高计算精度和运算速度。推导精确的积分公式，用典型体目标对改进的算法进行验证，进而计算复杂组合体目标的后向散射场。矩形截面直进气道的计算结果不仅比学者Hao Ling的计算结果更为精确，而且计算速度提高了两倍。三面角反射器的计算结果与实验结果基本吻合。最后，计算了某型航母甲板上复杂目标组合体目标的RCS，并与近似法和像素法的计算结果进行了比较分析。     

后向台阶层流分离剪切层特性研究

应用激光测速仪对生向台阶层流边界层分离产生的剪切层在再附前的发展进行了测量，得到了时均速度、湍流强度及剪切层厚度等的分布特性，并研究了来流湍流度对台阶后流场湍强度分布及剪切层发展的影响，实验中发现台阶后有一个流速降低的区域，对此从涡动力学的角度作了解释。     

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

为了满足两侧进气布局飞行器的乘波前体与进气道一体化设计要求,提出了一种进口水平投影可控的流线追踪内收缩进气道设计方法。基于马赫数分布可控的轴对称基准流场,在指定进口水平投影为椭圆的条件下,采用该方法设计了内收缩进气道并在设计点(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.     

Sunspot activity and cosmic ray modulation at 1 a.u. for 1900–2013

The descent of sunspot cycle 23 to an unprecedented minimum of long duration in 2006–2009 led to a prolonged galactic cosmic ray (GCR) recovery to the highest level observed in the instrumental era for a variety of energetic charged particle species on Earth, over a wide range of rigidities. The remarkable GCR increase measured by several ground-based, balloon-borne, and detectors on a satellite is described and discussed. It is accompanied by a decrease in solar wind velocity and interplanetary magnetic field at 1 a.u., reaching the lowest values since measurements of the solar wind began in October 1963; the solar polar field strength (μT) measured at the Wilcox Solar Observatory (WSO) is also significantly reduced compared to prior cycles since the start of the program in 1976, the polar field in the northern hemisphere reversed in June 2012 and again in February 2014, that in the southern hemisphere reversed in July 2013. If updates of WSO data confirm the second reversal in northern solar hemisphere, it would pose a serious challenge to the Dynamo Theory. The long-term change in solar behavior may have begun in 1992, perhaps earlier. The physical underpinnings of these solar changes need to be understood and their effect on GCR modulation processes clarified. The study discusses the recent phenomena in the context of GCR modulation since 1900. These happenings affected our empirical predictions for the key parameters for the next two sunspot cycles (they may be progressively less active than sunspot cycle 24) but it enhanced support for our prediction that solar activity is descending into a Dalton-like grand minimum in the middle of the twentyfirst century, reducing the frequency of the coronal mass ejections; they determine the space weather affecting the quality of life on Earth, radiation dose for hardware and human activities in space as well as the frequency of large Forbush decreases at 1 a.u.     

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

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

GH4169合金反挤压成形模具磨损

基于修正 Archard 磨损模型,采用数值模拟方法系统分析了 GH4169合金反挤压成形过程中各挤压工艺参数对模具磨损的影响规律。结果表明：在选取的参数范围内,挤压凸模最易产生磨损失效的区域为凸模圆角处,模具最大磨损深度随凸模圆角半径及坯料预热温度的增大而降低,随摩擦系数的增大而增大;当挤压速率小于100 mm ／s时,模具最大磨损深度随挤压速率的增大而减小,当挤压速率大于100mm ／s 时,模具最大磨损深度随挤压速率的增大先增大后减小。最佳工艺参数坯料预热温度1020℃,摩擦系数0．05,变形速率100mm ／s,模具预热温度300℃时模具磨损量最小,为9．28×10－3 mm。     

一种基于实数域特征值分解的前后向平滑算法

为减小算法的运算量,对常规前后向空间平滑(FBSS)算法进行改进,提出了一种基于实数域特征值分解的前后向平滑(RDFBSS)算法。算法根据前后向平滑算法协方差矩阵的中心-厄米特对称特点,通过实值化变换,使FBSS算法的特征值分解从复数域转为实数域,大幅减少了算法的运算量,且不影响解相干和统计性能,提高了实用性。仿真结果证实了算法的有效性。     

光线追踪求交算法的研究及其实现

本文对光线追踪求交算法进行了深入讨论。文中结合形体设计实例说明了此算法的正确性和实用性。     

Mean energy of response to galactic cosmic ray spectrum: IMP 8 and Climax neutron monitor

Solar modulations of galactic cosmic ray (GCR) intensity contain a wealth of information about their transport in the heliosphere. To extract this information from the data one studies the dependence of the observed modulations on the mean energy of response of detectors providing data for the analyses. There is a great deal of confusion about the detector energy response to GCR spectrum in the literature. We present a preliminary report on the computations of the mean energy of response for the Climax neutron monitor (CL/NM) and IMP 8 cosmic ray nuclear composition instrument to GCR protons for 1973–1998, covering the solar cycles 21 and 22. We find that for penetrating proton channel on IMP 8 the mean energy changes by a factor of over two whereas for the neutron monitor the change is only 21%. However, the corresponding change for the computed modulation function is a factor of about 3.5.