平面湍流混合层的准相似性理论

本文应用周培源的理论,给出了不可压缩平面湍流混合层的一级近似解,我们考虑连续性方程,平均运动方程和二阶脉动速度关联方程,忽略二阶脉动速度关联方程中的三阶关联项,并引进自模拟假设,得到了和实验符合的平均速度,二阶脉动速度关联项和湍流微尺度等的理论计算结果。     

湍流两相流的脉动速度联合PDF输运方程

概率密度函数（PDF）的方法是构造两相湍流模型的一种重要的方法.构建气体-颗粒速度联合PDF输运方程的关键是颗粒所见气体微团速度的Langevin方程.首先由N-S方程出发,精确推导出颗粒所见气体微团脉动速度的Langevin方程,进而通过理论分析表明,对比通常采用的颗粒所见气体微团瞬时速度的Langevin方程而言,采用前者能有效地减少关联量的统计偏差.最后,给出了颗粒-气体脉动速度的联合PDF输运方程.     

视密度加权的时平均统一二阶矩两相湍流模型

在气粒两相湍流的双流体模型中,颗粒相的视（表观）密度是有脉动的,在时平均的统一二阶矩（USM）模型中出现了和颗粒数密度或视密度脉动有关的项和方程,使模型方程比较复杂。实际上,用LDV或PDPA测量的流体（用小颗粒代表）和颗粒速度都是颗粒数加权平均的结果。因此,在视密度加权平均基础上推导两相湍流模型更为合理。通过推导和封闭了视密度加权平均的统一二阶矩模型（MUSM）方程组,改进了两相速度脉动关联的封闭,并引入了颗粒遇到的气体脉动速度及其输运方程。MUSM模型可以减少所用方程数,节省计算量。视密度加权平均的统一二阶矩两相湍流模型是一种对原有时间平均的统一二阶矩模型和改进和发展。     

用PDF方程法分析悬沙垂线浓度分布

从颗粒运动的PDF(概率密度分布函数)输运方程出发。建立颗粒相的质量、动量和脉动速度二阶矩方程．对于明渠二维恒定均匀流。利用垂向动量方程导出了新的泥沙扩散方程。表明颗粒脉动强度梯度、升力、重力沉降和紊动扩散都影响悬沙运动。说明了传统扩散方程的不足．理论分析了水沙两相物理属性和水流条件对泥沙扩散系数和浓度分布特征的影响。并通过细颗粒试验资料的分析进行了定量研究．     

可压缩流场中气泡脉动数值模拟

在应用边界元方法对气泡动力学的研究中, 绝大多数模型是建立在不压缩势流理论基础之上, 针对可压缩流场中气泡运动特性的研究很少. 从波动方程出发, 分别在气泡运动前期和后期对波动方程进行简化, 得到气泡运动局部和全局简化方程, 采用双渐进方法对简化方程进行匹配, 提出了考虑流场可压缩性的非球状气泡运动模型. 该模型的计算结果与Prospertti 等的解析结果吻合很好, 气泡脉动最大半径和内部最大压力随气泡脉动逐渐减小. 基于该模型对比了自由场中药包爆炸考虑可压缩性与不考虑可压缩性的计算结果, 发现考虑可压缩性气泡射流速度较小, 随后基于该模型计算了刚性边界下气泡的运动特性.     

汶川地震滑坡与地震参数及地质地貌因素之间的相关关系

采用大涡模拟方法,模拟了槽道湍流,得到了不同雷诺数下槽道湍流的结果. 在此基 础上,研究了平均速度、雷诺应力、脉动动能和脉动速度均方根的分布;讨论了平均速度的 壁面律问题;给出了雷诺应力、脉动动能和脉动速度均方根随雷诺数的变化规律,其中雷诺 应力、脉动动能给出了定量公式.     

稠密气固两相流各向异性颗粒相矩方法

基于气体分子动力学和颗粒动理学方法,考虑颗粒速度脉动各向异性,建立颗粒相二阶矩模型.应用初等输运理论,对三阶关联项进行模化和封闭.考虑颗粒与壁面之间的能量传递和交换,建立颗粒相边界条件模型.采用Koch等计算方法模拟气固脉动速度关联矩.考虑气体-颗粒间相互作用,建立稠密气体-颗粒流动模型.数值模拟提升管内气固两相流动特性,模拟结果表明提升管内颗粒相湍流脉动具有明显的各向异性.预测颗粒速度、浓度和颗粒脉动速度二阶矩与Tartan等实测结果相吻合.模拟结果表明轴向颗粒速度脉动强度约为平均颗粒相脉动强度的1.5倍,轴向颗粒脉动能大约是径向颗粒脉动能3.0倍.     

近壁湍流统计性质研究

采用大涡模拟方法,模拟了槽道湍流,得到了不同雷诺数下槽道湍流的结果. 在此基础上,研究了平均速度、雷诺应力、脉动动能和脉动速度均方根的分布;讨论了平均速度的壁面律问题;给出了雷诺应力、脉动动能和脉动速度均方根随雷诺数的变化规律,其中雷诺应力、脉动动能给出了定量公式.      

湍流边界层相干结构剩余脉动雷诺应力项的实验测量

用IFA300恒温热线风速仪和X形二分量热线探针以采样间隔小于湍流耗散时间尺度的分辨率精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、法向速度分量的时间序列信号。用条件采样和相位平均技术提取了相干结构猝发过程中相干结构剩余脉动雷诺应力和随机脉动对相干结构贡献的雷诺应力的条件相位平均波形。基于理论上对湍流相干结构动量方程中随机脉动对相干结构贡献的雷诺应力和相干结构剩余脉动雷诺应力项的分析,对两种雷诺应力项进行了对比研究。研究发现,相干结构剩余脉动雷诺应力项在数值上具有和随机脉动对相干结构贡献的雷诺应力相同的数量级,表明在相干结构动力学模型方程中,相干结构剩余脉动雷诺应力项并不像以前估计的那样可以忽略不计。     

FV/MC混合算法求解轴对称钝体后湍流流场

介绍一种有限容积／Monte Carlo结合求解湍流流场的相容的混合算法．有限容积法求解Reynolds平均的动量方程和能量方程，Monte Carlo方法求解模化的脉动速度—频率—标量联合的PDF方程．将该算法发展到无结构网格，探讨了在无结构网格中实现两种方法的耦合，包括颗粒定位，颗粒场和平均场之间数据交换等问题．并以二维轴对称钝体后湍流流场作为算例，比较了计算结果与实验结果．     

The new method of solving velocity correlation functions and application to the plane turbulent wake

The results of plane turbulent wake given by Chou Peiyuan^[1] are considered as the first order approximation and put into the equations of turbulent fluctuation. The equations are solved numerically within the range of micro-scale by means of spectrum method. The double, triple and quadruple fluctuating velocity correlations are obtained by computation. They are in good agreement with experimental results.     

Computation of Compressible Flows using a Two-equation Turbulence Model on Unstructured Grids

This paper presents a numerical method for solving compressible turbulent flows using a k - l turbulence model on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multi-stage Runge-Kutta time stepping scheme, while the turbulence equations are advanced using a multi-stage point-implicit scheme. The positivity of turbulence variables is achieved using a simple change of dependent variables. The developed method is used to compute a variety of turbulent flow problems. The results obtained are in good agreement with theoretical and experimental data, indicating that the present method provides a viable and robust algorithm for computing turbulent flows on unstructured meshes.     

Pulsatory Turbulent Compressible Fluid Flow and Propagation of Pressure Waves in a Channel

The properties of the damping coefficient and phase velocity of propagation of small-amplitude pressure waves as functions of the oscillation frequency are investigated for the turbulent flow of a weakly compressible fluid in a circular pipe. The wall friction is found by solving numerically the equation of motion and the relaxation equations for the turbulent shear stress and viscosity which provide the basis for a turbulent transfer model developed for unsteady conditions. The properties are explained in terms of an analysis of the calculated data on turbulent transfer. The results obtained are compared with experiments.     

Numerical simulation of swirling flow in diffusers

A reduced form of Navier–Stokes equations is developed which does not have the usual minimum axial step size restriction. The equations are able to predict accurately turbulent swirling flow in diffusers. An efficient single sweep implicit scheme is developed in conjunction with a variable grid size domain-conforming co-ordinate system. The present scheme indicates good agreement with experimental results for (1) turbulent pipe flow, (2) turbulent diffuser flow, (3) turbulent swirling diffuser flow. The strong coupling between the swirl and the axial velocity profiles outside of the boundary layer region is demonstrated.     

Numerical viscosity reduction in the resolution of the shallow water equations with turbulent term

A first‐order finite volume model for the resolution of the 2D shallow water equations with turbulent term is presented. An upwind discretization of the equations that include the turbulent term is carried out. A method to reduce the excess of numerical viscosity (or diffusion) produced by the upwinding of the flux term is proposed. Two different discretizations of the turbulent term are compared, and results for uniform distributions of the viscosity are presented. Finally, two discretizations of the time derivative which are more efficient than Euler's are proposed and compared. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical modeling of a two-dimensional vertical turbulent two-phase jet

The results of numerically modeling two-dimensional two-phase flow of the “gas-solid particles” type in a vertical turbulent jet are presented for three cases of its configuration, namely, descending, ascending, and without account of gravity. Both flow phases are modeled on the basis of the Navier-Stokes equations averaged within the framework of the Reynolds approximation and closed by an extended k-? turbulence model. The averaged two-phase flow parameters (particle and gas velocities, particle concentration, turbulent kinetic energy, and its dissipation) are described using the model of mutually-penetrating continua. The model developed allows for both the direct effect of turbulence on the motion of disperse-phase particles and the inverse effect of the particles on turbulence leading to either an increase or a decrease in the turbulent kinetic energy of the gas. The model takes account for gravity, viscous drag, and the Saffman lift. The system of equations is solved using a difference method. The calculated results are in good agreement with the corresponding experimental data which confirms the effect of solid particles on the mean and turbulent characteristics of gas jets.     

DNS and scaling laws from new symmetry groups of ZPG turbulent boundary layer flow

The Lie group, or symmetry approach, developed by Oberlack (see e.g. Oberlack [26] and references therein) is used to derive new scaling laws for various quantities of a zero pressure gradient turbulent boundary layer flow. The approach unifies and extends the work done by Oberlack for the mean velocity of stationary parallel turbulent shear flows. From the two-point correlation (TPC) equations the knowledge of the symmetries allows us to derive a variety of invariant solutions (scaling laws) for turbulent flows, one of which is the new exponential mean velocity profile that is found in the mid-wake region of flat-plate boundary layers. Further, a third scaling group was found in the TPC equations for the one-dimensional turbulent boundary layer. This is in contrast to the Navier–Stokes and Euler equations, which have one and two scaling groups, respectively. The present focus is on the exponential law in the outer region of turbulent boundary layer corresponding new scaling laws for one- and two-point correlation functions. A direct numerical simulation (DNS) of a flat plate turbulent boundary layer with zero pressure gradient was performed at two different Reynolds numbers Re=750,2240. The Navier–Stokes equations were numerically solved using a spectral method with up to 140 million grid points. The results of the numerical simulations are compared with the new scaling laws. TPC functions are presented. The numerical simulation shows good agreement with the theoretical results, however only for a limited range of applicability. PACS 02.20.-a, 47.11.+j, 47.27.Nz, 47.27.Eq