脉动压力脉动速度变形平均项的表达式

脉动压力脉动速度变形平均项,最初是Rotta[21]改写压力梯度做功项得来的.但是,这项处理起来都很困难,从Rotta开始,以后Launder等人都对这项做过一些假定.本文根据脉动速度所满足的方程解出脉动压力,然后进而求出脉动压力乘上脉动速度变形的平均值,得到了脉动压力脉动速度变形平均项的完整表达式.这个表达式说明了Rotta和Launder等人的有限表达式是有一定道理的.本文所得的完整表达式分为两种情形加以讨论.一种是几种涡旋不分开的情形,另一种是三种涡旋分开考虑的情形.由此,本文为雷诺应力模式和三涡旋模式等湍流模式提供了完整的脉动压力脉动速度变形平均项的表达式.     

层状二维流动的基本方程式

在很多海洋、大气等二维流动问题中所用的动力学方程往往沿用推广后的河流水力学方程或"纳维-斯托克斯方程"其中把湍流阻力项写成这样的方程式和湍流阻力项用到实际问题上去,无疑是存在着极大的局限性,而将导致矛盾百出.本文则从雷诺方程出发,把所有的物理量沿深度加以平均,求出平均以后的物理量所满足的运动方程,连续方程和扩散方程.     

多孔介质两相驱动问题的特征差分和Schwarz区域分裂并行算法

本文考虑多孔介质两相驱动问题的数值解法。用混合元方法求解压力方程，可同时得到速度和压力的近似；对浓度方程，给出了两类特征差分与Ｓｃｈｗａｒｚ型区域分裂引结合的数值格式，以减小对流项产生的数值弥散，减小所处理问题的规模和实行并行计算。     

可混溶驱动问题的超收敛性

本文讨论多孔介质中两相可混溶渗流驱动问题的有限元方法，采用一致网格剖分、指标为k的Raviart－Thomas空间对压力作混合有限元逼近，用正则剖分、逼近阶为l的标准有限元方法处理浓度方程，通过核函数对有限元解作卷积进行局部平均确定非线性项的系数，得到了浓度误差H1范数的超收敛估计，经高阶插值，得到了整体高精度的逼近．     

延拓Chaplygin气体黎曼解的流逼近极限过程中的集中性研究(英文)

本文研究了延拓Chaplygin气体的黎曼解在流逼近极限过程中的集中现象和狄拉克激波的两种形成机制问题.利用相平面分析法和广义特征分析法,构造出了延拓Chaplygin气体的整体黎曼解,并获得了两个结果:当压力消失时,延拓Chaplygin气体的包含两个激波的解收敛到输运方程的狄拉克激波解;当压力项趋近于广义Chaplygin压力项时,延拓Chaplygin气体的包含两个激波的解收敛到广义Chaplygin气体的狄拉克激波解.结论推广到了延拓Chaplygin气体.     

螺旋型旋风分离器两相流场的数值模拟

对螺旋型旋风分离器进行了两相流场的三维数值模拟．气体流场通过求解三维N-S方程得到,湍流模型采用了雷诺应力模型．计算结果表明,旋风分离器内部的流场分为两部分:螺旋通道内比较稳定的流场和筒体中心区域的复合涡结构流场．对颗粒运动轨迹的计算表明,颗粒在入口处的初始位置对颗粒分离有比较显著的影响．同时得到了不同入口速度下颗粒的分级效率曲线,并给出了气体流量对旋风分离器性能的影响,结果显示:气体流量的增加会提高分离效率,但同时导致压力损失的急剧增加．     

基于网格变动的油水驱动半正定问题迎风混合元方法

在网格随时间变动的有限元空间上研究了不可压缩的两相渗流驱动问题.分别对饱和度方程扩散矩阵正定和半正定的情形,提出了基于网格变动的迎风混合元方法混合元逼进压力方程,饱和度方程的对流项采用迎风格式来处理,扩散项则采用推广的混合元来逼进.在网格任意变动的情形下得到几乎最优的误差估计;对正定问题的格式进行改进,即在两个网格之间投影变化时采取近似解的线性构造,可以得到与固定网格时相同的最优收敛阶.     

凝析气田压力曲线特征分析

本文采用一个带汇项的流动方程描述了凝析气井在压力大于最大凝析压力时的压力时间变化特征.根据Duhamel叠加原理建立了反映凝析作用的汇项;引入了描述这一汇项的两个新的物理参数:凝析强度R_D和凝析弛豫时间λ_D.由汇项在流动方程中的作用,给出了凝析作用影响井口压力的定量表述,阐明了凝析作用影响井口压力的机理.本文给出了这一流动方程外边界条件为无限大、内边界以定产量生产情况下的解析解.本文的结论可用于试井理论.     

泥石流固液分相流速计算方法研究

泥石流固液分相流速是泥石流对岸坡、防治结构冲击、磨损机理的核心问题．将泥石流体简化为具有相同粒径的固相和具有相同力学性质的液相,基于泥石流体为沿流动方向的一维两相流体,运用两相流理论建立了泥石流固液分相流速控制方程．构建了泥石流平均压力、彻体力及平均表面力的计算方法,尤其通过浆体的Binhanm体流变方程、Bagnold颗粒相互作用试验成果建立了控制体平均表面力计算方法;建立了固液两相流速比例系数,以及理论固相流速与实际流速的比例系数．据此求解控制方程得到了固液分相流速计算方法,该方法既可同时适用于粘性泥石流和稀性泥石流,也可在泥石流爆发以后通过现场采集沉积物分析反求泥石流爆发期间的分相流速．工程实例分析显示,该方法计算结果与实测结果吻合较好．     

伴有排水的两相饱和介质动力问题的LAMB积分公式

饱和土在动力作用下伴有排水的反应是工程中常见的现象.由于饱和土以两相饱和介质模拟,在已经求得的集中力作用下,两相饱和介质基本解Green函数和三项流相Green函数基础上,由Betti定理出发,推出排水问题的Lamb积分公式和应力公式,并代入相关参数,通过计算得到了中心扩张源问题的位移场、排水状态、孔隙压力等解答.最后给出了当两相介质蜕化为单相介质后,与经典解答比较的结果.     

固液两相流中的K-ε双方程湍流模式*

建立了固液两相流中更一般的K-ε双方程湍流模式。模化了固相和液相的连续方程、动量方程及K方程和ε方程。该湍流模型考虑了固液两相间速度的滑移,颗粒间的作用及相间作用。使用本文所建立的湍流模型,数值预测了一管湍流中的沙水混合流动,其预测结果与实验结果比较一致。     

Numerical modeling of turbulence-induced interfacial instability in two-phase flow with moving interface

The present paper introduces a new interfacial marker-level set method (IMLS) which is coupled with the Reynolds averaged Navier–Stokes (RANS) equations to predict the turbulence-induced interfacial instability of two-phase flow with moving interface. The governing RANS equations for time-dependent, axisymmetric and incompressible two-phase flow are described in both phases and solved separately using the control volume approach on structured cell-centered collocated grids. The transition from one phase to another is performed through a consistent balance of kinematic and dynamic conditions on the interface separating the two phases. The topological changes of the interface are predicted by applying the level set approach. By fitting a number of interfacial markers on the intersection points of the computational grids with the interface, the interfacial stresses and consequently, the interfacial driving forces are easily estimated. Moreover, the normal interface velocity, calculated at the interfacial markers positions, can be extended to the higher dimensional level set function and used for the interface advection process. The performance of linear and non-linear two-equation k–ε turbulence models is investigated in the context of the considered two-phase flow impinging problem, where a turbulent gas jet impinging on a free liquid surface. The numerical results obtained are evaluated through the comparison with the available experimental and analytical data. The nonlinear turbulence model showed superiority in predicting the interface deformation resulting from turbulent normal stresses. However, both linear and nonlinear turbulence models showed a similar behavior in predicting the interface deformation due to turbulent tangential stresses. In general, the developed IMLS numerical method showed a remarkable capability in predicting the dynamics of the considered two-phase immiscible flow problems and therefore it can be applied to quite a number of interface stability problems.     

A 3-dimensional Eulerian–Eulerian CFD simulation of a hydrocyclone

Hydrocyclones are used in mineral industries for classification and separation of solid particles of different sizes and densities suspended in water medium. In the present study an Eulerian–Eulerian CFD simulation of a solid–liquid hydrocyclone has been carried out taking into account two solid phases and one liquid phase. The average size of the larger particle was 0.6117 and that of the smaller particle was 0.09875 mm. Three separate momentum balance equations for the three phases have been considered unlike that in the mixture model where a single momentum equation is solved for the three phases. Two turbulent models i.e. the Reynolds stress model (RSM) and the standard k–ε model were studied. Comparison of the two turbulence models showed slight variation in prediction of the velocity profile and the separation efficiency. The maximum deviation between the two models was observed near the wall where the stress was maximum for larger size particles.     

Influence of Mass Loading on Particle-Laden Turbulent Pipe Flow

A CFD code in the framework of OpenFOAM was validated for simulations of particle-laden pipe and channel flows at low to intermediate mass loadings. The code is based on an Eulerian two-fluid approach with Reynolds-averaged conservation equations, including turbulence modeling and four-way coupling. Pipe flow simulations of particles in air against gravity were conducted at Reynolds numbers up to 50000. The particle mass loading was varied and its effect on the mean velocities and turbulent fluctuations of the two phases was studied. Special attention was paid to the influence of mass loading on the centerline velocity and the wall shear velocity of the fluid phase for various flow parameters and particle properties. Empirical correlations were established between these two quantities and the flow Reynolds number, particle Reynolds number, Stokes number and particle to fluid density ratio for a range of particle mass loadings. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Direct numerical simulation of turbulent flows in eccentric pipes

A numerical algorithm was developed for solving the incompressible Navier-Stokes equations in curvilinear orthogonal coordinates. The algorithm is based on a central-difference discretization in space and on a third-order accurate semi-implicit Runge-Kutta scheme for time integration. The discrete equations inherit some properties of the original differential equations, in particular, the neutrality of the convective terms and the pressure gradient in the kinetic energy production. The method was applied to the direct numerical simulation of turbulent flows between two eccentric cylinders. Numerical computations were performed at Re = 4000 (where the Reynolds number Re was defined in terms of the mean velocity and the hydraulic diameter). It was found that two types of flow develop depending on the geometric parameters. In the flow of one type, turbulent fluctuations were observed over the entire cross section of the pipe, including the narrowest gap, where the local Reynolds number was only about 500. The flow of the other type was divided into turbulent and laminar regions (in the wide and narrow parts of the gap, respectively).     

Computation of dilute two-phase flow in a pump

This paper is a report on a joint project between academia and industry which is concerned with computation of dilute two-phase flow through a pump in turbulent condition. The flow field for the continuous phase is computed using the Reynolds averaged Navier–Stokes equations together with mixing length turbulence modeling. The dispersed phase is treated using the Lagrangian approach by tracking it's trajectory along which the information is passed. It is found that the bubbles and small solid particles flow out of the chamber (between the rotating impeller and the casing wall) with the conveying fluid. The solid particles of relatively bigger sizes accumulate at the low pressure zones near the cashing wall or the rotating shaft.     

明渠中变密度变粘度牛顿流体紊动基本方程式

本文利用Navier－Stokes方程及雷诺时均法则,导出了变密度变粘度牛顿流体的紊动微分方程式,并进一步导出了变密度变粘度牛顿流体在明渠中紊流流动时的运动微分方程式.文中首次提出了密度紊动应力与粘度紊动应力的概念.     

宾汉流体稠密两相湍流流动中的二阶矩模型

建立的Bingham流体稠密两相流动的二阶矩-颗粒动力论湍流模型（USM-theta模型）既体现了两相的作用,又体现了屈服应力所引起的附加项,并提出了USM-theta模型下考虑浓度修正值影响的两相湍流流动的算法．利用该模型对圆管内Bingham流体的单相湍流流动、稠密液固两相的湍流流动进行了计算,并和五方程湍流模型进行了比较,结果表明该模型的预测效果更好．利用USM-theta模型对含颗粒的Bingham流体的两相湍流流动进行了模拟,随着屈服应力的增加,Bingham流体相与颗粒相在管道中心附近的主流速度减小．液固两相湍流和Bingham流体两相湍流的计算结果表明屈服应力引起的附加项对流动有很重要的影响．     

Simulation of free turbulent particle-laden jet using Reynolds-stress gas turbulence model

Free two-phase flows occur in many practical applications, such as sprays or particle drying and combustion. This paper deals with mathematical modelling of a free turbulent two-phase jet. A steady, axisymmetric, dilute, monodisperse, particle-laden, turbulent jet injected into a still environment, has been considered. The model treats the gas-phase from an Eulerian standpoint and the motion of particles from a Lagrangian one. Closure of the system of time averaged transport equations has been accomplished by using a Reynolds-stress turbulence model. The particles–fluid interaction has been considered by the PSI-Cell concept. Both the effect of interphase slip and the effect of particle dispersion have been taken into account.     

密相两相流基本方程组的封闭性和分散相在壁面处的剪应力

本文对将连续相和分散相分别作连续介质考虑的两相流基本方程组,通过分析密相下两相流动的物理机制,提出了分散相的应力满足宾汉流体的本构方程,并用分散相在壁面处剪应力作为分散相的一个边界条件,从而得到了一个封闭的处理密相系统的两相流数学模型.最后,还讨论了分散相在壁面处剪切力的表达式.