Estimation of sea surface temperature, wind speed and water vapor with microwave radiometer data based on simulated annealing

A method for estimation of sea surface temperature, ocean wind speed and water vapor with microwave radiometer data based on simulated annealing is proposed. The proposed method shows about 60% improvement of sea surface temperature estimation accuracy in comparison to the existing method using Newton’s iterative algorithm.     

三维非结构网格的生成及优化

叙述了一种新的生成三维非结构网格的Delaunay方法,该方法能够高效地生成优质四面体网格。由于采用了局部剖分的方法来保证边界的完整性,因而不需要添加新的节点到原始表面三角形中。该方法在四面体外接球的球心处插入新的节点,在插入全部内点后,采用删除边、边/面交换、网格光顺等方法对生成的网格进行了优化,基本能删除质量差的四面体。     

在非结构网格上求解非线性航空声学问题的高精度有限体积法

介绍了有限体积数值方法,该方法适用于在任意非结构网格上求解线性和非线性的航空声学问题.本方法基于角点-中心的多参数格式,可在笛卡尔网格上达到六阶精度,对于可能的不连续性采用了自适应耗散.通过一系列算例研究了该方法的特性,结果表明:在模拟谐振型管中的噪声抑制中,所提出的方法是很有效的.     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

基于扫掠法的六面体网格生成算法及实现

提出了一种基于扫掠法的六面体网格生成算法,该算法通过源和目标曲面参数域间的近似最小二乘仿射映射将源曲面的网格投影到目标曲面,并利用背景网格插值的方法计算内结点,从而将2.5维实体生成六面体网格。该算法避免了空间网格投影过程中耗时的“寻根”过程,并且不需要导致效率降低的网格优化,能够处理各种2.5维实体,生成的网格质量较好。     

Automatic differentiation based discrete adjoint method for aerodynamic design optimization on unstructured meshes

A systematic methodology for formulating, implementing, solving and verifying discrete adjoint of the compressible Reynolds-averaged Navier-Stokes (RANS) equations for aerodynamic design optimization on unstructured meshes is proposed. First, a general adjoint formulation is con-structed for the entire optimization problem, including parameterization, mesh deformation, flow solution and computation of the objective function, which is followed by detailed formulations of matrix-vector products arising in the adjoint model. According to this formulation, procedural components of implementing the required matrix-vector products are generated by means of auto-matic differentiation (AD) in a structured and modular manner. Furthermore, a duality-preserving iterative algorithm is employed to solve flow adjoint equations arising in the adjoint model, ensur-ing identical convergence rates for the tangent and the adjoint models. A three-step strategy is adopted to verify the adjoint computation. The proposed method has several remarkable features:the use of AD techniques avoids tedious and error-prone manual derivation and programming;duality is strictly preserved so that consistent and highly accurate discrete sensitivities can be obtained; and comparable efficiency to hand-coded implementation can be achieved. Upon the cur-rent discrete adjoint method, a gradient-based optimization framework has been developed and applied to a drag reduction problem.     

Atmospheric response to a Hot SST Event in November 2006 as observed by the AIRS instrument

Using Atmospheric Infrared Sounder (AIRS) products of atmospheric temperature and geopotential height, we investigate the atmospheric response to HE0611, which was found and investigated by [Qin, H., Kawamura, H., Sakaida, F., Ando, K. A case study of the tropical Hot Event in November 2006 (HE0611) using a geostationary meteorological satellite and the TAO/TRITON mooring array. J. Geophys. Res. 113, C08045, doi: 10.1029/2007JC004640, 2008]. HE0611 was formed by connecting two very high SST areas, HE0611-East and HE0611-West. The period-mean atmosphere temperatures at levels of 925 and 850 hPa in HE0611-West are higher, by about 0.5 K, than those in WE0611-East while the atmospheric temperatures at middle to high levels (700–300 hPa) are higher in HE0611-East. The period-mean geopotential heights HE0611-East are much lower than those in HE0611-West for the levels from the surface to 400 hPa. The mean geopotential heights from 400 hPa to 200 hPa are higher in HE0611-East. In the middle and high layers over HE0611-West, the atmosphere temperatures gradually decrease from 7th to 17th, and then increase significantly. The increase in HE0611-East starts from 15th November, which is earlier than that of HE0611-West. The geopotential heights in the high layer of both the areas also show corresponding behaviors. The lagged atmospheric response in the western part is confirmed by the correlation analysis. It emerges that the atmospheric response to HE0611 is well organized and associated with deep convention in HE0611-East and subsidence in HE0611-West. These are also consistent with the HE0611 features and evolution revealed by earlier HE studies.     

厚覆冰导线升力突变机理及数值模拟研究

针对新月形厚覆冰导线的升力系数在风攻角15°附近存在突变的问题,分别采用基于k-ωSST湍流模型的雷诺时均法和大涡模拟(LES)的数值方法对新月形厚覆冰导线在风攻角10°~20°范围进行了模拟。通过对比两种数值方法计算得到的覆冰导线气动力系数、流场结构和表面风压,发现LES方法能够更好地捕捉新月形覆冰导线表面的小尺度涡结构,得到的覆冰导线气动力参数计算结果与风洞试验数据高度吻合;而k-ωSST湍流模型难以模拟壁面上小尺度涡,捕捉不到升力系数的突变。根据覆冰导线不同壁面区域的压力分布,发现上侧壁面处的涡结构影响整体流场,并在下侧壁面曲率、来流夹角和壁面切线方向共同作用下导致升力系数突变。LES的气动力参数模拟结果可为覆冰导线防舞提供参考。     

Gurney Flap绕流气动特性研究

采用C-H型多块结构网格和SST k-ω湍流模型求解二维定常雷诺平均NS方程,对NACA4412单段翼型后缘附近和NACA63(2)-215B多段翼型襟翼后缘附近的Gurney Flap绕流进行了数值模拟。讨论了局部网格对计算结果的影响,重点研究了Gurney Flap位置和高度对气动特性的影响。计算结果表明:对于NACA4412单段翼型,Gurney Flap位置对气动性能的影响较小,但随着高度的增加,零升攻角和失速攻角减小,最大升力系数增大,同时阻力系数和力矩系数增加;对于NACA63(2)-215B多段翼型,Gurney Flap对气动性能的影响类似于NACA4412单段翼型情形。