斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。     

Experimental investigation and correlation development of jet impingement heat transfer with two rows of aligned jet holes on an internal surface of a wing leading edge

Extensive experimental studies on the heat transfer characteristics of two rows of aligned jet holes impinging on a concave surface in a wing leading edge were conducted, where50000 Rej 90000, 1.74 H/d 27.5, 66° a 90°, and 13.2 r/d 42.03. The finding was that the heat transfer performance at the jet-impingement stagnation point with two rows of aligned jet holes was the same as that with a single row of jet holes or the middle row of three-row configurations when the circumferential angle of the two jet holes was larger than 30°. The attenuation coefficient distribution of the jet impingement heat transfer in the chordwise direction was so complicated that two zones were divided for a better analysis. It indicated that: the attenuation coefficient curve in the jet impingement zone exhibited an approximate upside-down bell shape with double peaks and a single valley; the attenuation coefficient curve in the non-jet impingement zone was like a half-bell shape, which was similar to that with three rows of aligned jet holes; the factors,including Rej, H/d and r/d, affected the attenuation coefficient value at the valley significantly.When r/d was increased from 30.75 to 42.03, the attenuation rates of attenuation coefficient increased only by 1.8%. Consequently, experimental data-based correlation equations of the Nusselt number for the heat transfer at the jet-impingement stagnation point and the distributionof the attenuation coefficient in the chordwise direction were acquired, which play an important role in designing the wing leading edge anti-icing system with two rows of aligned jet holes.     

Rio de la Plata estuarine system: Relationship between river flow and frontal variability

This paper deals with the application of SeaWIFS images to characterize spatial and temporal variability of fronts in the Rio de la Plata estuarine system over the period 2000–2003. We aim to depict the relationship between river outflow and variability of fronts’ loci on monthly to ENSO-related timescales and the influence of the winds along Rio de la Plata (axial winds) on the abrupt changes in frontal dynamics over synoptic timescales. During the studied period both La Niña (July 1999–June 2000) and El Niño (April 2002–May 2003) events induced significant displacements of fronts. Three distinct fronts were analyzed between river, estuarine, coastal and marine waters of the Rio de la Plata: Main Turbidity Front, Main Marine Front, and Secondary Marine Front. Their number, location and separation seem to be mainly related to river outflow and second, to fresh (>8 m/s) axial winds. During low discharge periods (i.e. summer time and/or La Niña events) these winds induce abrupt changes in the location of fronts (100–200 km) and greater separation between them over synoptic timescale, whereas during high river discharge or ENSO years some of the variability of fronts location is explained by the river’s outflow fluctuations, especially by the high variability of the River Uruguay discharge.     

水平表面气流剪切作用下的水膜厚度

飞机结冰表面上的液态水受气流吹拂作用会发生向后溢流,从而影响结冰区域范围及防冰系统设计;为了获得水膜流动规律,对水平平板表面上气流剪切驱动的水膜流动进行了实验测量和建模分析。通过水膜流动风洞试验台产生高速气流驱动水膜的流动,使用色散共焦位移计测量同一位置的水膜在不同时刻的厚度变化,结果表明气-液界面由底层薄水膜和多种尺度的波动组成,具有变化速度快随机性强的特点。通过水膜厚度随气流速度及水膜雷诺数的变化规律,发现平均水膜厚度与两者均呈现出单调非线性的依赖关系。基于薄水膜流动理论和平均水膜厚度实验结果,提出了高速气流剪切作用下的气-液波动界面剪切因子计算式,适用于风速17.8~52.2m/s,水膜雷诺数26~128之间的平板水膜流动计算。     

圆截面曲线管道内二次流动的Galerkin解

本文推导了曲线直角坐标系内的Ｎ－Ｓ方程。结果表明Ｇａｌｅｒｋｉｎ方法可以克服摄动法的小参数局限，得到了纯数值法无法求得的增 解析解。     

亚跨超CFD软件平台研制

亚跨ＣＦＤ软件平台研制的宗旨是综合国内外计算流体动力学的最新研究成果，用计算机软件技术的新发展，建立一个面向对象的计算流体动力学的研究开发环境，加快计算流体动力学最新成果向工程应用的转换速度。本文介绍初级版本的研究进展和研究成果，主要匝数值网格生成技术，流场解算器，面向对象的界面设计等三部分组成。     

圆球绕流场的尾涡分析和升阻力研究

采用自编的分块算法程序模拟了雷诺数在20～1000之间的圆球绕流场.在Re=25时捕捉到流动分离,与试验结果一致.对Re=250时稳态非对称尾流区内流体的输运情况进行研究,证实圆球绕流场中纵向对称面的存在.对于非稳态流动,主要研究了阻力和横向力的时间变化规律,并分析了频谱特征.计算发现横向合力主要在横截面上两个很窄的区间内偏移.经过几种不同涡定义方法的优劣测试,本文采用Hunt等人介绍的Q-定义方法对非稳态绕流场的尾涡结构进行了直观描述.     

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.     

针栓式喷注单元膜束撞击雾化混合过程数值模拟

为了全面认识针栓式喷注器喷雾场结构，基于自适应网格加密技术和分三相计算的PLIC VOF（Piecewise Linear Interface Calculation Volume of Fluid）方法对针栓式喷注单元膜束撞击雾化混合过程进行了仿真分析，通过对两路推进剂分别进行界面追踪，获得了膜束撞击雾化混合过程的详细结构特征，与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。以此为基础对膜束撞击的喷雾场结构、撞击变形过程、流场涡结构、雾化破碎典型特征及破碎后的雾化混合分布特征进行了识别分析，结果表明：膜束撞击形成了液束未穿透液膜和液束穿透液膜2种不同的喷雾扇结构。膜束撞击形成的喷雾扇呈"Ω"形，膜束同时发生弯曲变形和横截面变形。另外，膜束撞击同时受到正压和剪切应力作用，导致了一系列复杂涡流现象，使得相互作用增强，雾化混合均增强，这也是膜束撞击喷注构型优于膜膜撞击的本质原因。最后，还发现膜束撞击喷雾场液滴分布呈现分区结构特征，分别是液束控制主导的上雾化区、液膜控制主导的下雾化区及夹在中间的混合区，实际中应兼顾雾化特性和混合特性，选取中等动量比膜束撞击，这可为针栓式喷注器的理论研究和工程设计提供重要参考。     

多孔介质空心圆柱绕流的摄动解

本文主要讨论了绕多孔介质空心圆柱的有势流动问题。运用摄动方法,给出了速度势的四阶近似,然后得到了圆柱表面的无量纲的压力系数和阻力系数。结果发现在势流的情况下多孔介质圆柱受到附加阻力作用,并且在圆柱的喷射或吸收流体的特征速度与来流特征速度的比值是一个小量的时候,阻力系数与特征速度比成简单的线性关系。