A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

Stationary flow fields prediction of variable physical domain based on proper orthogonal decomposition and kriging surrogate model

In this paper a new flow field prediction method which is independent of the governing equations, is developed to predict stationary flow fields of variable physical domain. Predicted flow fields come from linear superposition of selected basis modes generated by proper orthogonal decomposition(POD). Instead of traditional projection methods, kriging surrogate model is used to calculate the superposition coefficients through building approximate function relationships between profile geometry parameters of physical domain and these coefficients. In this context,the problem which troubles the traditional POD-projection method due to viscosity and compressibility has been avoided in the whole process. Moreover, there are no constraints for the inner product form, so two forms of simple ones are applied to improving computational efficiency and cope with variable physical domain problem. An iterative algorithm is developed to determine how many basis modes ranking front should be used in the prediction. Testing results prove the feasibility of this new method for subsonic flow field, but also prove that it is not proper for transonic flow field because of the poor predicted shock waves.     

人体热调节的有限元分析

本文根据人体的热调节机制建立了二维人体热调节的有限元模型，并设计了相应的实验对该模型进行了验证，实验结果与计算结果吻合较好，本模型对于以人体－环境热交换为内容的工程实践和理论研究有较高的应用价值。     

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

小型弹体转动惯量的测定

阐述了对称式三线摆的结构原理及该摆微角与大角摆动周期的计算方法 ,利用此摆可以测量构件的转动惯量。当 RL <130 、φ0 <6 0°时 ,该装置的测量准确度优于 1%。     

Self-similar solutions to stellar winds with cosmic rays

A two-fluid model is used to study the time evolution of stellar winds including the dynamical effect of cosmic rays. Neglecting the diffusion of cosmic rays, we seek self-similar solutions to spherically symmetric winds with a termination shock. The velocity upstream of the shock is taken to be zero. Physical solutions are those that can connect the shock to the star with the velocity approaching zero at the star. Two parameters govern the behaviour of the solutions, namely, the ratio of the upstream sound speed to the shock speed (in an inertial frame) and the gravitation potential of the central star. In some parameter regimes, no physical solution is possible.     

整体壳段典型特征高效模块化编程技术

整体壳段为航天产品核心零部件，具有结构复杂、尺寸精度高、材料去除率大、工艺过程复杂、加工周期长等特点。在壳段数控加工过程中，存在编程难度大、重复工作量大、效率低及编程质量与编程人员能力水平密切相关等问题。为了适应高效、高可靠性的任务需求，依托成熟商用CAM软件，通过对整体壳段典型特征的分析，基于UG NX软件平台研究并开发了典型特征自动化编程模块，有效提高了加工程序的质量稳定性，并极大降低了工艺人员的编程操作量，提高了工艺设计效率。     

一种冲击响应谱试验设备校准装置的研制

介绍了冲击响应谱基本概念及其试验设备的分类与工作原理，参考现行冲击试验台检定规程提出了针对该设备的校准方法，设计了一种基于虚拟仪器技术的新型冲击响应谱试验设备校准装置，给出了典型冲击波形（半正弦、矩形、后峰锯齿）的冲击响应谱仿真计算结果并进行了验证，计算结果与理论结果完全一致，证明算法准确可靠，能够实现对冲击响应谱试验设备的校准。     

How do interplanetary shock impact angles control the size of the geoeffective magnetosphere?

In this paper, we investigate temporal and spatial magnetosphere response to the impact of interplanetary (IP) shocks with different inclinations and speeds on the Earth’s magnetosphere. A data set with more than 500 IP shocks is used to identify positive sudden impulse (SI⁺) events as expressed by the SuperMAG partial ring current index. The SI⁺ rise time (RT), defined as the time interval between compression onset and maximum SI⁺ signature, is obtained for each event. We use RT and a model suggested by Takeuchi et al. (2002) to calculate the geoeffective magnetospheric distance (GMD) in the shock propagation direction as a function of shock impact angle and speed for each event. GMD is a generalization of the geoeffective magnetosphere length (GML) suggested by Takeuchi et al. (2002), defined from the subsolar point along the X line toward the tail. We estimate statistical GMD and GML values which are then reported for the first time. We also show that, similarly to well-known results for RT, the highest correlation coefficient for the GMD and impact angle is found for shocks with high speeds and small impact angles, and the faster and more frontal the shock, the smaller the GMD. This result indicates that the magnetospheric response depends heavily on shock impact angle. With these results, we argue that the prediction and forecasting of space weather events, such as those caused by coronal mass ejections, will not be accurately accomplished if the disturbances’ angles of impact are not considered as an important parameter within model and observation scheme capabilities.     

基于动网格技术的飞机舱弹分离过程 CFD/ RBD耦合数值仿真

飞机舱弹分离瞬间特性严重影响导弹/炸弹的打击精度。采用计算流体动力学(CFD)和刚体动力学(RBD)模型耦合求解的数值仿真方法,应用网格弹性光顺和局部重构相结合的动网格技术,有效地避免了刚体运动而引起的网格畸变,仿真计算了炸弹舱弹分离后俯仰角度变化过程 .仿真结果表明,舱弹分离时向炸弹施加 1个合适的抛 射力矩,有助于炸弹在下降过程中快速地调整其打击姿态。