非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ＥＮＯ（ＥｓｅｎｔｉａｌｙＮｏｎｏｓｃｉｌａｔｏｒｙＳｃｈｅｍｅ）型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶ＴＶＤＲｕｎｇｅＫｕｔａ方法．最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算     

三单元细化准则在随机振动分析中的应用

从弹性波波动理论出发,提出了一种新的有限元网格细化准则——三单元细化准则,可以解释模态数细化准则实施过程中出现的模态数突变现象。对于以封闭式盒段组合结构为特征的卫星电子设备,讨论了其有限元模型的简化建模方法。利用三单元细化准则,验证了某卫星电子设备随机振动加速度响应的NASTRAN软件分析结果满足工程计算精度的要求。     

用于立方体卫星的高增益天线技术与发展趋势

随着MarCO与RainCube的成功在轨验证,立方体卫星已经在深空探测与遥感领域展现出独特的优势。立方体卫星搭载的天线不仅需要提供良好的辐射特性,还必须满足立方体卫星苛刻的体积和重量限制,这种限制条件在设计大口径高增益天线时尤为重要。首先简要介绍了立方体卫星的起源、发展及其应用前景,其次重点介绍了立方体卫星高增益天线的技术特点、研究现状及其在轨验证情况,最后归纳总结了立方体卫星高增益天线技术的未来发展趋势,为我国后续开展立方体卫星研究与应用提供参考。     

Models for estimating daily photosynthetically active radiation in oceanic and mediterranean climates and their improvement by site adaptation techniques

In this work Photosynthetically Active Radiation (PAR) in oceanic and mediterranean climates is modeled. Twenty-two different models have been developed and tested: eleven Multilinear Regression (MR) models and eleven Artificial Neuron Network (ANN) models, using combinations of variables such as Global Horizontal Irradiance (GHI), Global Extraterrestrial Irradiance (G₀), Temperature (T) and Relative Humidity (RH). Data provided by Satellite Application Facility on Climate Monitoring (CM SAF) are used to develop and train the models, while the models have been validated using field data from four stations located in Spain, covering the different study climates.According to the results, zones with different climate conditions need different models, both for the case of MR and ANN. The results show the need of including the GHI in all models in order to obtain accurate estimates; in fact, the presence of more variables only improves slightly the results in mediterranean climate, while in oceanic climate no improvement is observed.On the other hand, comparing MR and ANN models, ANN models did not show better results than those of MR models in no one of the cases studied. Regarding the climate, both types of models are clearly better for the mediterranean case than for the oceanic one. In order to improve the performance of the model for oceanic climate a correction based on the site adaptation technique was carried out. The good results obtained by this technique fully justify its use.The best proposed models provide better performance than other models which are restricted to certain locations. Besides, the clustering technique based on the PAR variable, used in this work, allows obtaining useful models for a whole region. Finally, another advantage of this methodology is that there is no need of ground measurements for its development, except for the site adaptation technique.     

Experimental and numerical study on pressure oscillation in a combustor with rotary valve

A rotary valve coupled in a combustor assembly can generate periodic pressure oscillations inside the combustor and can be used to study the combustion instability in the combustion chamber of a rocket engine. This paper proposes a cold gas flow experimental system based on a rotary valve and a corresponding rotation model. A 3D numerical model is proposed to obtain the transient flow inside the rotary valve, and the dynamic mesh technique and User-Defined Functions (UDFs) are adopted to implement a swing motion instead of a rotary motion. Several cold gas flow experiments are carried out at rotating speeds of 75 225, and 375 rpm to verify the validity of the numerical model. The effects of rotating speed, stroke length ratio and radius of the RED (Rotor Exhaust Duct) on the pressure oscillations are studied using this numerical model. The results show that the maximum and peak-to-peak values of the pressure oscillations gradually decrease with increasing rotating speed. The ratio of the corresponding peak-to-peak value to the maximum pressure (pressure amplitude ratio) is reduced from 1.81 to 0.6%. The stroke length ratio affects the pressure waveform because it leads to a change in the time spent in the non-exhaust stage. When the stroke length ratio is 0, the waveform closely resembles a sine wave. With the increase in the stroke length ratio, the pressure waveform exhibits a more square or a triangular wave shape. Finally, a high-frequency and high-amplitude pressure oscillation can be obtained by appropriately increasing the radius of the RED.     

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

一种求解轨迹优化问题的改进多分辨率技术

针对轨迹优化多分辨率技术可能发生细化遗漏或者细化失败的缺陷，设计了一种改进多分辨率技术。该方法通过在每次网格细化迭代过程中增加两层高分辨率节点并且引入节点检测算法，确保在轨迹的非光滑区域能够高效、持续地细化网格。针对原始多分辨率技术的初始网格的节点数目存在二进限制的弊端，定义了一种新的广义二分网格，使得多分辨率技术的初始网格节点数目可以为任意奇数。采用多个轨迹优化算例验证了本文方法的有效性并且与其它方法进行了对比。此外，仿真结果还表明，网格充分细化对目标函数影响不大，但是能够更加准确地满足路径约束和终端状态约束。     

霜状冰结冰数值模拟研究

提出了一种三维霜冰结冰的数值模拟方法。基于欧拉-拉格朗日法计算空气-水滴两相流流场,获得了表面水滴撞击特性以及撞击量,求解质量守恒和能量守恒方程获得网格单元内的结冰量,三维霜冰的结冰计算程序采用用户自定义（UDF）编写,FLUENT提供的动网格功能实现了结冰后的网格重构。计算了4°攻角,液态水含量为1 g/m3,水滴直径20μm,结冰时间6 min,温度分别为-19.4°和-28.3°时NACA0012翼型的结冰,并与相关文献结果进行了对比,计算结果与文献实验结果吻合较好,验证了方法的合理性和可靠性。分析了不同的粒径、液态水含量、来流速度等对撞击特性以及结冰的影响,随着粒径、来流速度的增加,水滴收集系数、撞击区和结冰厚度增大,随着液态水含量的增加,水滴收集系数和撞击区是不变的,但是结冰厚度增加。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

应用动网格技术计算着陆过程中缓冲发动机喷流流场

利用Fluent软件对着陆过程中缓冲发动机的喷流流场进行数值计算,发动机与地面相对位置的不断变化使计算域形状不断改变.为了保证计算的正确性,必须采用动网格技术.介绍了实现移动网格的基本方法,并结合使用网格弹性平滑移动和网格重构方法实现了着陆过程中缓冲发动机喷流流场的非定常数值计算;得到了随发动机着陆产生的喷流流场的压强...