高精度PPM格式在湍流燃烧中的应用

利用处理三维可压缩粘性流体流动问题中的沉浸边界法,并结合基于PPM方法的高精度TVD格式,对三维方形管道中部的圆柱火焰绕流及惰性气体绕流问题进行了数值模拟。计算湍流时采用大涡模拟(LES),化学反应速率采用EBU漩涡破碎模型。通过计算结果与实验结果的比较,发现高精度PPM格式能精确模拟两类圆柱绕流问题。计算中还发现,火焰圆柱绕流算例中,在火焰翻越圆柱前,由于燃烧的膨胀作用,使得火焰正面前的未燃气体流动并形成惰性气体绕流,这与无燃烧时的惰性气体绕流类似。但当火焰翻越圆柱过程中及完全翻越圆柱后,两种算例绕流流场出现明显变化。     

三维大钝头体简谐振动绕流的数值模拟

用数值模拟方法研究了大钝头体外形（双子星座飞船）在超声速下俯仰振动的绕流特性，数值模拟的出发方程为三维全Ｎ－Ｓ方程，格式为二阶迎风ＴＶＤ格式，研究结果给出了物体在俯仰振动过程中激波、分离结构和物面压力分布等非定常特征。     

三角翼的双襟翼控涡作用的数值模拟研究

对装有“前端襟翼”和“前缘襟翼”的７４°后掠三角翼的不可压缩流场作了数值模拟,以研究襟翼的旋涡控制作用．数值模拟是用拟压缩性方法求解一般曲线坐标系下的三维不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,时间离散用向后Ｅｕｌｅｒ差分,空间无粘项的离散用二阶迎风ＴＶＤ格式,所得的离散方程用对角化形式的近似隐式因子分解格式求解．湍流模型用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模式．计算了三种平面形状的机翼在迎角范围为１０°～５０°的绕流和气动特性．计算和实验的比较表明,襟翼向下偏转可以推迟旋涡破裂,且对提高机翼的减阻能力、升阻比和改善失速前后的气动特性有明显效果,双襟翼具有更佳的控涡效果．     

基于反馈力浸入边界法模拟复杂动边界流动

浸入边界法是模拟流固耦合的重要数值方法之一。本文采用反馈力浸入边界方法,对旋转圆柱和水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟。其中,固体边界采用一系列离散的点近似代替,流体为不可压缩牛顿流体,使用笛卡尔自适应加密网格,利用有限差分法进行求解。固体对流场的作用通过构造适宜的反馈力函数实现。本文首先通过旋转圆柱绕流的计算结果同实验结果进行对比,吻合较好,验证了该计算方法的可靠性。然后针对水电站水力过渡过程中水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟,得到导叶动态绕流后的流场分布特性和涡结构的演化特性。     

翼型冰增长和结冰影响的数值模拟研究

翼面结冰是威胁飞行安全的重要因素之一.本文对于霜状结冰,在给定的结冰气候条件下,采用数值方法,预测整个结冰过程,建立冰增长模型,求解水滴轨迹运动方程,分析结冰过程中翼型绕流流场的变化,以及冰层的发展和形成过程.针对给定的三种不同冰型,即钝头型、尖头型和双角型,进一步采用结构化网格生成技术,结合中心有限体积法和LU-SGS隐式算法,利用B-L代数湍流模型,完成了绕流流场的N-S方程数值模拟,分析不同形状的冰型对翼型绕流及气动特性的影响.计算结果与试验数据进行了对比,表明采用的方法是正确、可行的.     

绕椭球的低速流动研究

理解和预测绕椭球的流动对指导飞行器和潜艇等交通工具的设计具有很强的工程意义. 近年来, 针对椭球绕流开展了大量的实验和数值模拟研究. 对有攻角下椭球绕流分离的定性描述和定量研究, 促进了对三维分离的辨识和拓扑研究. 文章对流场特性进行了分析, 介绍了分离对气动力、噪声、尾迹的影响, 以及实验条件对流动的影响. 上述定常流动与非定常机动过程之间存在明显差异, 非定常机动过程不能作为定常或准定常问题处理, 在机动过程中, 分离出现明显延迟, 气动力出现明显变化. 随后介绍了数值模拟在求解绕椭球流动中的进展, 当前求解雷诺平均的N-S方程湍流模式仍然是解决绕椭球大范围分离流动的主要工程方法, 大涡模拟和分离涡模拟等也逐渐得到了广泛应用. 受限于计算能力, 直接数据模拟只能用于较低雷诺数, 在高雷诺数流动中还不适用. 非定常机动过程的数值模拟较定常状态, 与实验结果的差距要大一些. 最后, 介绍了对椭球绕流场转捩的研究进展, 对T-S转捩与横流转捩的机理和辨识已经较为准确, 数值模拟结果与实验结果基本相符, 但对再附转捩的认识还不够清晰, 尤其是迎风面, 因此椭球绕流转捩的研究还需要依靠实验.      

超声速平板边界层斜波失稳转捩过程研究

以5阶迎风和6阶对称紧致格式混合差分求解三维可压缩滤波Navier-Stokes方程，对Mach 数为4.5, Reynolds数为10000的空间发展平板边界层湍流进行了大涡模拟. 时间推进采用 紧致存储3阶Runge-Kutta方法，亚格子尺度模型为修正Smagorinsky涡黏性模型. 通过在 入口边界叠加一对线性最不稳定第一模态斜波扰动，数值模拟得到了平板层流边界层失稳转 捩直至湍流的演化过程. 对流场转捩过程中瞬时量及统计平均量的分析表明，数值模拟结果 与理论吻合，得到的Y型剪切层、交替\Lambda涡结构以及转捩后期的发卡涡结构的发展 变化与相关文献结果一致，湍流流谱定性合理.     

高温真实气体底部流动的NS方程数值计算

本文数值模拟了高超音速飞行时钝锥的底部流动。采用轴对称NS方程并考虑真实气体效应。湍流模型采用修正的Baldwin-Lomax涡粘性代数模型,数值方法空间离散对流项采用显式NND格式,粘性项采用中心差。时间离散采用三阶的龙格-库塔法。真实气体模型采用考虑七种组分四种反应的汉森模型。给出了底部流场的压力和温度分布及各组分的浓度分布。可以看出在近底部区域高速流-绕过拐角就产生一回流旋涡区。由于温度变化很大,气体的热力学特性受气体离解、复合和振动能激发的影响。所以整个流动过程变得十分复杂。     

基于车身绕流的低雷诺数湍流模型改进研究

本文将汽车绕流模块化为各典型局部流动,通过常用湍流模型对各典型局部流动进行数值模拟,结果验证了湍流模型对转捩的捕捉能力是准确模拟汽车绕流的关键. 在分析汽车绕流分离及转捩机理的基础上,优化了稳态和瞬态求解方法,改进了湍流模型对转捩的预测能力,进而提高了湍流模型在汽车流场模拟上的精度. 针对汽车绕流的稳态问题,将流线曲率因子及 响应阈值引入 LRN $k$-$\varepsilon $ (low Reynolds number $k$-$\varepsilon $) 模型,获得了一种能够更准确预 测转捩的改进低雷诺数湍流模型 (modified LRN $k$-$\varepsilon $),改善了原模型对湍流耗散率的过强依赖性及全应力发展预测不足等问题;针对汽车绕流瞬态求解,通过分析 RANS/LES 混合湍流模型的构造思想及特点,引入约束大涡模拟方法,结合本文提出的改进的 LRN $k$-$\varepsilon $ 湍流模型,提出了一种能准确捕捉转捩现象 的转捩 LRN CLES 模型. 分别将改进的模型用于某实车外流场和风振噪声仿真中,通过 Ansys Fluent 求解器计算,并将计算结果与常用湍流模型的仿真结果、HD-2 风洞试验结果和实车道路实验结果进行对比,表明改进后的湍流模型能够更准确模拟复杂实车的稳态和瞬态特性,为汽车气动特性的研究提供了可靠理论依据及有效数值解决方法.      

RESEARCH ON IMPROVEMENTS OF LRN TURBULENCE MODEL BASED ON FLOW AROUND AUTOMOBILE BODY 1)

本文将汽车绕流模块化为各典型局部流动,通过常用湍流模型对各典型局部流动进行数值模拟,结果验证了湍流模型对转捩的捕捉能力是准确模拟汽车绕流的关键. 在分析汽车绕流分离及转捩机理的基础上,优化了稳态和瞬态求解方法,改进了湍流模型对转捩的预测能力,进而提高了湍流模型在汽车流场模拟上的精度. 针对汽车绕流的稳态问题,将流线曲率因子及 响应阈值引入 LRN $k$-$\varepsilon $ (low Reynolds number $k$-$\varepsilon $) 模型,获得了一种能够更准确预 测转捩的改进低雷诺数湍流模型 (modified LRN $k$-$\varepsilon $),改善了原模型对湍流耗散率的过强依赖性及全应力发展预测不足等问题;针对汽车绕流瞬态求解,通过分析 RANS/LES 混合湍流模型的构造思想及特点,引入约束大涡模拟方法,结合本文提出的改进的 LRN $k$-$\varepsilon $ 湍流模型,提出了一种能准确捕捉转捩现象 的转捩 LRN CLES 模型. 分别将改进的模型用于某实车外流场和风振噪声仿真中,通过 Ansys Fluent 求解器计算,并将计算结果与常用湍流模型的仿真结果、HD-2 风洞试验结果和实车道路实验结果进行对比,表明改进后的湍流模型能够更准确模拟复杂实车的稳态和瞬态特性,为汽车气动特性的研究提供了可靠理论依据及有效数值解决方法.     

Secondary flow coefficient of overbank flow

This paper presents a 2D analytical solution for the transverse velocity distribution in compound open channels based on the Shiono and Knight method （SKM）, in which the secondary flow coefficient （K-value） is introduced to take into account the effect of the secondary flow. The modeling results agree well with the experimental results from the Science and Engineering Research Council-Flood Channel Facility （SERC-FCF）. Based on the SERC-FCF, the effects of geography on the secondary flow coefficient and the reason for such effects are analyzed. The modeling results show that the intensity of the secondary flow is related to the geometry of the section of the compound channel, and the sign of the K-value is related to the rotating direction of the secondary flow cell. This study provides a scientific reference to the selection of theK-value.     

The presence of slug flow in horizontal two-phase flow

One of the flow regimes occurring in horizontal two-phase flows is characterized by periodic large waves “surging” along the tube. This flow, called “slug” flow, has been frequently observed in low and high pressure gas liquid systems, but it has been noticed that slugging is absent in certain liquid-liquid two-phase systems. A method is developed giving the necessary conditions for the presence of slug flow. This method quantitatively explains the observed absence of slugging in certain liquid-liquid flows.     

Fabrication of rigid and flexible refractive-index-matched flow phantoms for flow visualisation and optical flow measurements

A method for the construction of both rigid and compliant (flexible) transparent flow phantoms of biological flow structures, suitable for PIV and other optical flow methods with refractive-index-matched working fluid is described in detail. Methods for matching the in vivo compliance and elastic wave propagation wavelength are presented. The manipulation of MRI and CT scan data through an investment casting mould is described. A method for the casting of bubble-free phantoms in silicone elastomer is given. The method is applied to fabricate flexible phantoms of the carotid artery (with and without stenosis), the carotid artery bifurcation (idealised and patient-specific) and the human upper airway (nasal cavity). The fidelity of the phantoms to the original scan data is measured, and it is shown that the cross-sectional error is less than 5% for phantoms of simple shape but up to 16% for complex cross-sectional shapes such as the nasal cavity. This error is mainly due to the application of a PVA coating to the inner mould and can be reduced by shrinking the digital model. Sixteen per cent variation in area is less than the natural patient to patient variation of the physiological geometries. The compliance of the phantom walls is controlled within physiologically realistic ranges, by choice of the wall thickness, transmural pressure and Young’s modulus of the elastomer. Data for the dependence of Young’s modulus on curing temperature are given for Sylgard 184. Data for the temperature dependence of density, viscosity and refractive index of the refractive-index-matched working liquid (i.e. water–glycerol mixtures) are also presented.     

Effects of flow dispersers upstream of two-phase flow monitoring instruments

The effects of flow dispersing screens on the accuracy of two-phase flow monitoring instruments was studied with an air-water two-phase flow loop in which air and water flow rates were monitored separately. Flow dispersers were installed at the inlet of an instrumented spool piece containing a drag disk and a turbine meter. For the present tests, the spool piece was mounted vertically and flow was always down.When the turbine was upstream of the drag disk, two-phase mass flow rates indicated by the instruments were approximately one-third actual mass flows both in the presence and absence of flow dispersers. However, when the drag disk was upstream of the turbine, the use of flow dispersing screens upstream of the drag disk resulted in flow measurements within $\text{± 10\%}$ of actual flow rates over the range 0–160 kg/m²· s.     

Two-phase flow distribution of air–water annular flow in a parallel flow heat exchanger

The air and water flow distribution are experimentally studied for a round header – flat tube geometry simulating a parallel flow heat exchanger. The number of branch flat tube is 30. The effects of tube outlet direction, tube protrusion depth as well as mass flux, and quality are investigated. The flow at the header inlet is identified as annular. For the downward flow configuration, the water flow distribution is significantly affected by the tube protrusion depth. For flush-mounted configuration, most of the water flows through frontal part of the header. As the protrusion depth increases, more water is forced to the rear part of the header. The effect of mass flux or quality is qualitatively the same as that of the protrusion depth. Increase of the mass flux or quality forces the water to rear part of the header. For the upward flow configuration, however, most of the water flows through rear part of the header. The protrusion depth, mass flux, or quality does not significantly alter the flow pattern. Possible explanations are provided based on the flow visualization results. Negligible difference on the water flow distribution was observed between the parallel and the reverse flow configuration.     

Transient flow redistribution in annular two-phase flow

A model is described for the prediction of transient flow redistribution in vertical annular two-phase flow. The model is based on an analysis of the local parameters controlling the flow and takes account of the diffusive motion of entrained droplets and the delay time for change in the wave structure on the film. Comparisons are made with experimental results on inlet effects and it is shown that the wall injection experimental results can be described by the model. The jet injection results are not fitted by the model and it is shown that some additional deposition mechanism must be present.     

Analysis of entry flow to determine elongation flow properties revisited

The minimisation technique proposed by Binding (J. Non-Newtonian Fluid Mech., 27 (1988) 173) was used in our Generalised Engineering Bernoulli Equation framework to relate the entry pressure and stress power. We arrived at a final result similar to Binding's using assumed kinematics. Through subsequent assumptions to the kinematics we finally arrive at a result exactly equivalent to Cogswell's technique (Trans. Soc. Rheol., 16 (1972) 383). Thus, these two techniques are related in this general framework. The techniques were used to predict elongation flow properties of a polymer melt and polymer solution. The results for the polymer melt clearly show Cogswell's technique is adequate at high elongation rates. All these techniques require minimisation of the stress power with respect to the flow volume and discussion is given as to the validity of this minimisation technique. In addition, the approximate variational technique we propose gives clears limits as to when a technique, such as Cogswell's, can be applied.     

Prediction of the slug-to-churn flow transition in vertical two-phase flow

An assessment is made of the various viewpoints on the slug-to-churn flow transition in vertical upward flow in the light of recent experimental results obtained at Harwell Laboratory. It is found that the flooding model of McQuillan & Whalley and the bubble entrainment model of Barnea & Brauner give satisfactory results at low and high liquid flow rates, respectively. An improved model for flooding, which takes account of the effect of the falling film, has been proposed. It is shown that this new model is in good agreement with experimental results at both low and high liquid flow rates.