求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

飞机着陆进近穿越风切变的轨迹响应

大气紊流和风切变对飞行活动影响较大,尤以低空风切变危害最大,严重威胁着民航飞行的安全。鉴于飞行器对抗大气扰动的实际能力有限,正确的飞行对策是回避或改出危险的大气扰动区。本文对某训练飞机着陆进近遭遇风切变时飞机轨迹的变化进行了模拟仿真计算。结果表明,采用恒定迎角策略,即采用最大允许迎角（抖动迎角）和最大油门时,改出效果最佳。同时计算出了某训练飞机在水平风为0时,所能抵抗的最大垂直阵风强度;以及在垂直阵风为0时,所能抵抗的最大水平风强度。     

Estimates of vertical land motion along the southwestern coasts of Turkey from coastal altimetry and tide gauge data

The differences between coastal altimetry and sea level time series of tide gauges in between March 1993 and December 2009 are used to estimate the rates of vertical land motion at three tide gauge locations along the southwestern coasts of Turkey. The CTOH/LEGOS along-track coastal altimetry retrieves altimetric sea level anomalies closer to the coast than the standard along-track altimetry products. However, the use of altimetry very close to the coast is not found to improve the results. On the contrary, the gridded and interpolated AVISO merged product exhibits the best agreement with tide gauge data as it provides the smoothest variability both in space and time compared with along track altimetry data. The Antalya gauge to the south (in the Mediterranean Sea) and the Mentes/Izmir gauge to the west (in the Aegean Sea) both show subsidence while the Bodrum tide gauge to the south (in the Aegean Sea) shows no significant vertical land motion. The results are compared and assessed with three independent geophysical vertical land motion estimates like from GPS. The GIA effect in the region is negligible. The VLM estimates from altimetry and tide gauge data are in good agreement both with GPS derived vertical velocity estimates and those inferred from geological and archaeological investigations.     

悬停共轴双旋翼干扰流动数值模拟

基于非结构动态嵌套网格方法,通过求解三维非定常Euler方程对悬停共轴双旋翼直升机的复杂流场进行了数值模拟研究.双旋翼之间的干扰使得两者的拉力性能下降,上旋翼对下旋翼的干扰程度远大于下旋翼对上旋翼的影响程度.下旋翼处的下洗速度比单旋翼的要大得多,而上旋翼的下洗速度比单旋翼的略大一些,几乎相同.桨叶总星巨角增加,则下洗速...     

基于欧拉方程的高超声速飞行器的壁面流线生成计算

提出了一种新的利用表面流函数法计算高超声速飞行器表面流线分布的方法.首先提出了表面流函数的概念,并通过理论推导,得到了表面流函数与表面流线的关系;然后运用结构化网格求解三维Eu-ler方程,计算得到高超声速钝头体的边界层外缘外部无粘流场气流参数;最后利用无粘流场气流参数和表面流函数的方法计算了高超声速飞行器的精确表面流线分布.结果表明,在无攻角和攻角小于20°的情况下均可以得到较好的壁面流线分布.高精度的表面流线的得到,为利用边界层内粘性主导区域的积分方程等方法进一步精确预测高超声速飞行器表面的气动加热奠定了基础.      

飞行器外挂投放数值模拟

采用非结构重叠网格Euler方程的CFD方法数值模拟了外挂投放。该方法在重叠网格"挖洞"时,以"离壁面距离"为挖洞判据;在建立重叠区域插值关系时,对Neighbor to Neighbor(N2N)和Alternating Digital Tree(ADT)两种方法的效率进行了比较。采用机翼/挂架/带舵外挂物模型开展了验证计算。在此基础上,完成了某飞行器外挂投放的初步计算。     

A series of small scientific satellite with flexible standard bus

Japan Aerospace Exploration Agency has a plan to develop the small satellite standard bus for various scientific missions and disaster monitoring missions. The satellite bus is a class of 250–400 kg mass with three-axis control capability of 0.02 accuracy. The science missions include X-ray astronomy missions, planetary telescope missions, and magnetosphere atmosphere missions. In order to adapt the wide range of mission requirements, the satellite bus has to be provided with flexibility. The concepts of modularization, reusability, and product line are applied to the standard bus system. This paper describes the characteristics of the small satellite standard bus which will be firstly launched in 2011.     

基于粒子群算法的翼型优化设计

采用粒子群算法（PSO）对层流翼型进行了以提高升阻比为目标的优化设计。翼型的设计达到了设计要求,优化设计后的翼型气动特性也有显著地改善,这表明了粒子群算法应用于翼型气动优化设计的可行性。在优化设计的过程中,粒子采用递减惯性权重,以加强粒子初期的全局搜索能力与后期的局部搜索能力。翼型由解析函数线性叠加法表示,目标函数和粒子的适应度由基于二维欧拉方程的流场数值解来提供。