超声速流动中非线性EASM湍流模式应用研究

针对超声速复杂流动区域精确模拟的需要,发展了基于k-ω可压缩修正形式的非线性显式代数雷诺应力模式(EASM),提高了该模式对超声速复杂流动的数值模拟精度。通过对二维超声速凹槽和三维双椭球的数值计算表明,与SA和SST常规线性涡黏性湍流模式比较,非线性的EASM模式对大分离以及剪切层流动结构的刻画能力更精细,对剪切层再附区的压力及摩擦系数分布模拟更加精确;EASM模式能够准确地模拟二次激波引起的压强和热流分布情况。     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

反应堆三角形通道棒束间二次流动的数值研究

采用非线性Ｋ－ε湍流模式数值模拟三角形通道棒束中的二次流动，并考察其对流动和传热的影响。数值方法采用非正交曲线坐标系下求解控制方程的非交错网格方法。计算结果表明该模式能够较为有效地反映棒束中的二次流动，进一步分析表明二次流动有利于改善棒束中的流动和传热特性。     

非线性湍流模式研究及进展

现代湍流模式研究已经超出了经典的Ｂｏｕｓｓｉｎｅｓｑ涡粘性概念和线性的雷诺应力输运范畴,湍流运动过程中的非线性本质已成为模式研究人员所关心的中心问题。其目的在于使湍流模式能更加真实地再现湍流运动的复杂性,提高模式的适用范围,使复杂湍流能够得到合理的模拟,非线性湍流模式在解决复杂湍流运动的计算中已经取得可喜进展,正逐步应用于工程湍流的计算。同时,工程中的湍流问题计算也已走出了简单剪切流动类型及传统的ｋ－ε（及其它形式的）二方程模式框架,二阶矩封闭模式在先进的工程计算中已被用来解决诸如可压缩的空气动力学、发动机气缸及三维复杂几何场内等具有重要应用背景的流动问题,并逐步进入计算流体力学商业软件包。      

超声速平面剪切层声辐射涡模态数值分析

对Mc = 1.2二维超声速空间发展平面自由剪切层, 进行了扰动模态及流动结构的数值分析. 采用时空三阶改进MacCormack格式, 差分求解可压缩扰动Navier-Stokes方程, 直接数值模拟入口不同基频谐波扰动的非线性演化特征. 采用空间线性稳定性理论证明, 计算所促发的扰动波是声辐射涡模态. 扰动参数及特征函数分析显示, 声辐射涡模态是弱色散的快／慢两种外部模态, 在扰动对流Mach数为超声速一侧呈膨胀／压缩状辐射. 单频受迫扰动可无相差地促发多模态混合扰动波, 而在自然扰动条件下, 剪切层的稳定性受慢模态主导.     

基于格子Boltzmann方法的低雷诺数四圆柱绕流流动模式与受力特性模拟

利用格子Boltzmann方法模拟了雷诺数为100时,均匀来流条件下的二维菱形排布的四柱绕流现象,得到了不同柱间距比下的绕流流动模式及阻力变化规律。结果表明:圆柱互扰效应与柱间距比有关,当L/D≤1.2时为单钝体模式,圆柱互扰效应以临近效应为主;当1.2相似文献   

改进虚拟边界算法在超声速流动问题求解中的应用

提出了一种改进的虚拟单元浸没边界法, 并与一种高阶格式的有限差分算法相结合, 运用于求解超声速复杂几何绕流问题.该算法的核心思想在于在固体边界的内部和外部分别施加满足边界关系的作用点, 使得几何边界离散更加细化, 起到了壁面附近网格局部加密的作用.采用源空间内流体点作为反距离插值算法的重构点, 有效避免了插值点数目过少而与作用点相重合情况.通过对二维激波反射现象 (马赫数为 2.81) 和三维超声速球体绕流问题 (马赫数为 1.2) 的数值模拟, 与实验结果对比表明, 本文改进算法相对一般的虚拟边界法来说能显著提高数值精度, 减小计算误差.计算结果揭示了球体绕流中剪切层、压缩波系和尾迹的相互作用导致自由剪切层失稳的机理.剪切层厚度和湍流雷诺脉动经历了线性增长、大幅度震荡和小幅度波动三个阶段, 导致剪切层表面褶皱因子变化呈指数规律增长.其湍流结构表现出明显的各向异性, 具体在流向雷诺正应力在湍流脉动中占主导地位, 激波的压缩作用对不同方向雷诺正应力的影响存在空间迟滞效应.      

基于非结构网格的TTGC有限元格式的实现及在超声速流动中的应用

基于非结构网格系统,实现了时空三阶精度的TTGC有限元格式,并在三阶TTGC格式上发展了基于人工粘性的激波捕捉技术。在非结构网格下,采用这种方法对若干典型的超声速流动问题(SOD激波管、马赫数为3的前台阶流动以及马赫数为8的高超声速圆柱流动)进行了验证计算。结果表明,TTGC格式分辨率高,在粗糙网格下能够准确的模拟超声速流场中的激波、接触间断等复杂流动现象,并且能有效的控制间断附近的数值色散现象。与传统的有限体积方法相比,本文实现的TTGC有限元格式在模拟超声速流动问题方面具有格式精度高、数值耗散小等优点。     

对二维分离流涡黏性系数非线性分布的新认识

以弱非线性涡黏性模型为出发点,对Delery分离流动实验结果进行分析并获得了非平衡态分离区涡黏性系数与形状因子J之间的非线性关系. 该非线性关系显示在分离起始阶段,涡黏性系数较平衡态先减小,后增大;再附阶段,涡黏性系数较平衡态数值逐渐增大,并在再附点位置接近最大,而后又逐渐减小,恢复到平衡态水平. 总结涡黏性系数的这种非线性发展数学关系式,并将它应用于BL模型,在不添加微分方程的情况下发展出一种适用于分离流动的改进代数湍流模型. 对低速平板流动,跨声速,超声速以及高超声速分离流动的计算结果表明,该改进湍流模型可以较准确地模拟各类复杂分离流动,计算精度明显优于传统代数模型以及一些两方程模型,而计算工作量仍与BL模型相当. 这表明所提出的涡黏性系数非线性发展规律是正确的,且应用在二维分离流动中具有一定的普适性.      

逆压梯度边界层未分离情况下的湍流模式适用性研究Turbulence model applicability for boundary layer flow with adverse pressure gradient in attached case

边界层逆压梯度作用下的流动是许多工程中的一个基础问题,由于逆压梯度作用,流动形态复杂,使得数值模拟有很大的难度。基于雷诺平均纳维‐斯托克斯RANS（Reynolds Averaged Navier‐Stokes）方程对二维平板逆压梯度边界层作数值计算研究,选取6种代表性的湍流模式,得到局部摩擦系数的数值解,与实验值比较,发现k‐ω模式具有很好的精度。基于该湍流模式,给出了湍动能分布,该结果有助于认识逆压梯度边界层流动的复杂特征。     

关于二相流、多相流、多流体模型和非牛顿流等概念的探讨

本文分析了单相流、二相流和多相流等概念上的差异,也分析了单流体模型、双流体模型和多流体模型等概念上的差异,指出前面三种概念是按流动介质的客观物理构成划分的,而后者是按主观采用的研究方法划分的．目前这些概念在使用中存在一些混乱,如二相流与多相流,多相流与多流体模型等．本文还研究了扩散模型、非牛顿流模型和颗粒流模型等,指出前两种模型在分类上属于单流体模型,分析了非牛顿流模型、扩散模型和双（多）流体模型的特点和应用范围,最后,以泥石流为例讨论了以上概念的应用．      

A THREE-FLUID MODEL OF THE SAND-DRIVEN FLOW

ＡＴＨＲＥＥ－ＦＬＵＩＤＭＯＤＥＬＯＦＴＨＥＳＡＮＤ－ＤＲＩＶＥＮＦＬＯＷ￥(刘大有，董飞)ＬｉｕＤａｙｏｕ；ＤｏｎｇＦｅｉ（ＩｎｓｔｉｔｕｔｅｏｆＭｅｃｈａｎｉｃｓ，ＡｃａｄｅｍｉａＳｉｎｉｃａ，Ｂｅｉｊｉｎｇ１０００８０，Ｐ．Ｒ．Ｃｈｉｎａ）Ａｂｓ...     

The study on three-dimensional mathematical model of river bed erosion for water-sediment two-phase flow

Based on the tensor analysis of water-sediment two-phase flow, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent flow. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.     

Numerical predictions of the bifurcation of confined swirling flows

The bifurcation of confined swirling flows was numerically investigated by employing both the k-? and algebraic stress turbulence models. Depending upon the branch solution examined, dual flow patterns were predicted at certain swirl levels. In the lower-branch solution which is obtained by gradually increasing the swirl level from a low-swirl flow, the flow changes with increasing swirl number from the low-swirl flow pattern to a high-swirl flow pattern. In the upper-branch solution which is acquired by gradually decreasing the swirl level from a high-swirl flow, on the other hand, the flow can maintain itself in the high-swirl flow pattern at the swirl levels where it exhibits the low-swirl flow pattern in the lower branch. The bifurcation of confined swirling flows was predicted with either the k-? model or the algebraic stress model being employed. Both the k-? and algebraic stress models result in comparable and sufficiently good predictions for confined swirling flows if high-order numerical schemes are used. The reported poor performance of the k-? model was clarified to be mainly attributable to the occurrence of the bifurcation and the use of low-order numerical schemes.     

On Modelling Periodic Motion with Turbulence Closures

In periodic flows near walls transport effects may be considerably larger than in a steady turbulent boundary layer. The question explored in this contribution is, therefore, whether providing separate transport equations for each of the Reynolds stresses consequently leads to a better modelling of a periodic flow's evolution than an eddy viscosity scheme whose constitutive equation is inherently linked to the generation and dissipation terms being in balance (local equilibrium). Our conclusion is that, while a stress-transport scheme is indeed better equipped to reproduce the phenomena, it does not consistently out-perform the EVM over the range of flows studied. In some cases it is suggested that more attention must be paid to modelling diffusive transport in order to secure the benefits of second-moment closure. To illustrate sensitivity to diffusive transport, two different diffusion models are tested, one of which leads to different effective transport coefficients in each stress component. This revised version was published online in July 2006 with corrections to the Cover Date.     

高温流体与变物性固体耦合系统温度场及热变形的数值计算

在腔体内部有高温流体,其外壳等固体材料的物性温度而为化的瞬态非线性问题中,其边界条件的确定十分困难。本文采用流一固“混合流”模式进行温度场的计算,由于方程计及了瞬态、复杂形状、浮动外边界条件等因系,采用变步长的控制容积法及时域有限元进行计算。对某个高温燃气配气阀的温度场及热变形进行了数值计算,得到与实测相符的结果。     

Vertical integration of three-phase flow equations for analysis of light hydrocarbon plume movement

A mathematical model is derived for areal flow of water and light hydrocarbon in the presence of gas at atmospheric pressure. Vertical integration of the governing three-dimensional, three-phase flow equations is performed under the assumption of local vertical equilibrium to reduce the dimensionality of the problem to two orthogonal horizontal directions. Independent variables in the coupled water and hydrocarbon areal flow equations are specified as the elevation of zero gauge hydrocarbon pressure (air-oil table) and the elevation of zero gauge water pressure (air-water table). Constitutive relations required in the areal flow model are vertically integrated fluid saturations and vertically integrated fluid conductivities as functions of air-oil and air-water table elevations. Closed-form expressions for the vertically integrated constitutive relations are derived based on a three-phase extension of the Brooks-Corey saturation-capillary pressure function. Closed-form Brooks-Corey relations are compared with numerically computed analogs based on the Van Genuchten retention function. Close agreement between the two constitutive models is observed except at low oil volumes when the Brooks-Corey model predicts lower oil volumes and transmissivities owing to the assumption of a distinct fluid entry pressure. Nonlinearity in the vertically integrated constitutive relations is much less severe than in the unintegrated relations. Reduction in dimensionality combined with diminished nonlinearity, makes the vertically integrated water and hydrocarbon model an efficient formulation for analyzing field-scale problems involving hydrocarbon spreading or recovery under conditions for which the vertical equilibrium assumption is expected to be a satisfactory approximation.     

坡面流及土壤侵蚀动力学(I)——————坡面流

坡面流是坡面土壤侵蚀的根本动力,是研究土壤侵蚀动力学过程的基础.本文简要介绍了坡面流的基本特征和水力学特性,综述了坡面入渗产流过程、坡面流阻力,以及坡面流运动的数学描述和预报模型等方面的研究进展.指出复杂坡面条件下的坡面流运动,以及从小尺度过渡到大尺度将是今后坡面流研究的重要发展方向.      

低浓度固液两相流的颗粒相动理学模型

用广义Fokker-Planck扩散模型描述液相湍动对颗粒的挟带作用，用修正的BGK模型描述粒间碰撞效应，建立了封闭的颗粒相PDF输运方程．运用Chapman-Enskog迭代法求得方程的二阶近似解，获得颗粒相脉动速度二阶矩和三阶矩闭合关系．模型与颗粒流模型相容，与液相湍流闭合模型是否相容依赖于扩散模型的具体形式，并据此比较了不同的涡一颗粒作用模型．模型与二维明渠流轻质沙和天然沙试验资料符合很好．表明细小粒径颗粒能够充分跟随水流运动；大粒径颗粒的相间平均速度差和壁面滑移速度明显，近壁区内的颗粒沿流向和垂向脉动强度都可能大于水流，并存在一定程度的颗粒碰撞效应．     

Effect of collision and velocity model of lattice Boltzmann model on three-dimensional turbulent flow simulation

Various collision and velocity models of the lattice Boltzmann model (LBM) were compared to determine their effects on the efficiency of a three-dimensional homogeneous isotropic decaying turbulent flow simulation. We determined that a decrease in the number of velocities, in particular, 13-velocities, which can be used in the quasi-equilibrium lattice Boltzmann and in the multiple-relaxation time models (MRT), could considerably decrease the computational effort. However, decreasing the number of velocities deteriorates the stability and the accuracy of the results. By comparing the collision models, we also determined that the stability of the entropic lattice Boltzmann model (ELBM), and 19- and 27- velocity MRT is much higher than in other models. However, the numerical viscosity introduced by the ELBM underestimates the enstrophy, and the computational effort increases because of the calculation overhead required to solve the additional equations if special care is not given to the calculation.