粘性流体换热问题用直角坐标网格的贴体解法

采用直角坐标网格但通过在边界附近保留不规则控制体的非结构化直角坐标网格方法求解粘性不可压流体流动和换热问题.在边界附近切割单元的求解采用常规的方法处理,而不采用插值方法处理,并在同位网格上对控制方程组进行离散.因而该算法是贴体的,并可适用于解决复杂边界流动换热问题.通过对同轴环行空间和二维倾斜空腔内的对流换热问题的计算对该数值方法在传热问题中的应用进行验证.对比结果表明,计算结果与精确解和STAR-CD的结果在一定Re数和网格数时是很接近的,可以满足一定的精度要求,说明该数值计算方法可用于传热问题.      

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

过冷水滴撞击三维机翼的数值模拟

基于欧拉两相流理论,提出一种数值模拟过冷水滴撞击三维机翼的方法.根据水滴拟流体模型建立三维水滴控制方程;提出一种可穿透型壁面边界模拟水滴对机翼表面的撞击;由水滴流场的求解结果得到机翼表面的水滴收集特性.对ONERA M6机翼表面的水滴收集特性进行了研究,计算结果表明所提出的数值模拟方法是有效、可行的,可以为三维积冰的数值模拟研究奠定良好的基础.      

单次扫描空间推进算法在高超声速流动中的适用性研究

通过高超声速流动的三个典型算例,研究单次扫描空间推进(SSPNS)算法计算高超声速流场气动力、气动热的精度以及适用范围.结果表明:①在楔角不大于15°,流动无流向分离的情形下,用SSPNS算法计算压缩楔板流动是可行的;②对绕圆锥流动的计算表明,SSPNS算法对横向分离的模拟是值得信赖的;③对于12°攻角以内的小攻角绕翼流动,SSPNS算法的精度与时间迭代法基本一致;④空间推进法计算速度较时间迭代CFD(计算流体动力学)方法提高1～2个数量级.      

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

Oscillatory thermocapillary flows in simulated floating-zones with time-dependent temperature boundary conditions

This study deals with numerical simulations of the Maxus sounding rocket experiment on oscillatory Marangoni convection in liquid bridges. The problem is investigated through direct numerical solution of the non-linear, time-dependent, three-dimensional Navier-Stokes equations. In particular, a liquid bridge of silicon oil 2[cs] with a length L = 20 [mm] and a diameter D = 20 (mm) is considered. A temperature difference ΔT = 30 [K] is imposed between the supporting disks, by heating the top disk and cooling the bottom one with different rates of ramping. The results show that the oscillatory flow starts as an ‘axially running wave’, but after a transient time the instability is described by the dynamic model of a ‘standing wave’, with an azimuthal spatial distribution corresponding to m = 1 (where m is the critical wave number). After the transition, the disturbances become larger and the azimuthal velocity plays a more important role and the oscillatory field is characterized by a travelling wave. The characteristic times for the onset of the different flow regimes are computed for different rates of ramping.     

物面突变引起的二维分离流动的数值模拟

本文用一个单步二阶精度格式,对由物体截面突变引起的二维超声速分离流动进行了数值模拟。计算了近尾迹流场,并详细研究了来流M数和Re数对绕后台阶分离流场结构的影响,得到了与理论和前人计算一致的结果。另外,对超声速绕竖舨流动进行了数值计算,给出了细致的流场结构图画。