离散型湍流多相流动的研究进展和需求

离散型多相流动,指气体-颗粒(气-固)、液体-颗粒(液-固)、液体-气泡、气体-液雾以及气泡-液体-颗粒等两相或三相流动.这种类型的多相流动广泛存在于能源, 航天和航空, 化工和冶金,交通运输, 水利, 核能等领域.本文阐述了离散型多相流动的国内外基础研究,包括颗粒/液滴/气泡在流场中受流体动力作用力的研究, 颗粒-颗粒,液滴-液滴,气泡-气泡之间以及颗粒/液滴和壁面之间碰撞和聚集规律的研究,颗粒-气体和气泡-液体湍流相互作用的研究, 和数值模拟的研究,包括多相流动的雷诺平均模拟、大涡模拟和直接数值模拟的研究进展.最后, 归纳了目前尚待研究的需求.      

高浓度固-液两相流紊流的动理学模型

采用分子动理学方法,基于固-液两相流液相分子或颗粒相颗粒的Boltzmann方程,对Boltzmann方程分别取零矩和一次矩,则得到高浓度固-液两相流紊流的连续方程和动量方程,再和较成熟的低浓度两相流连续方程和动量方程比较,取低浓度两相流控制方程中较成熟合理的有关项和高浓度时由动理学方法推导出的颗粒间碰撞项,则得到高浓度固-液两相流紊流的最终控制方程:连续方程和动量方程.     

血液入口段两相流动的层流NS方程数值模拟

本文采用球形刚性颗粒悬浮于牛顿流体的血液模型，基于有相间滑移的两相流动层流Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，采用相间滑移算法（Ｉｎｔｅｒ－ＰｈａｓｅＳｌｉｐＡｌｇｏｒｉｔｈｍ）对圆管内定常轴对称血液入口两相流动问题进行了计算，获得了与实验结果吻合的计算结果，较好地模拟血液流动中红细胞的径向迁移现象，结果表明，采用二相流动模型研究血液流动是一种有前途的方法。     

气液两相流动与固壁相互作用耦合求解的研究

气液两相流动与固壁相互作用的研究是液滴撞击壁面运动研究的重要基础.以结合了VOF和Level Set两种方法优点的用于气液相界面追踪的复合Level Set-VOF方法和利用唯象分析方法建立的能够反映接触角滞后性及壁面性质对润湿过程影响的壁面润湿模型为基础,提出了气液两相流动与固壁相互作用耦合求解流程,给出了气液两相流动与固壁相互作用耦合求解过程中接触线速度的计算方法及边界条件的确定方法.通过与已有实验结果的对比,对提出的气液两相流动与固壁相互作用耦合求解方法的有效性进行了验证.     

航空发动机轴承腔润滑的气液两相均匀流研究

基于轴承腔中润滑油的两相均匀流动模型,采用湍流模式和有限差分数值方法计算轴承腔内三维定常N-S方程,对腔内润滑油的气液两相均匀流动特性进行研究,以获得气液两相均匀流条件下润滑油流场、压力场和速度场在轴承腔内的分布情况,分析转子转速、含气率和润滑油进口速度对润滑油出口压力以及润滑油与壁面之间剪切力的影响,同时对单相流和两相均匀流润滑性能差异进行比较.结果表明,转子转速、含气率和润滑油进口速度对润滑油出口压力和腔内壁面与润滑油间的平均剪切力具有不同影响,而采用2种流动模型计算出的轴承腔润滑油出口压力的差异较大,同时支持了开展航空发动机轴承腔润滑两相流动分析的必要性.     

气体速度对液膜在预膜板表面流动形态影响的LBM模拟

为了探究气流剪切作用对航空燃油在气动雾化喷嘴预膜板上流动形态的影响,首先对基于相场理论的两相流格子Boltzmann模型进行修正,并通过经典算例验证了修正后模型的准确性和可靠性.随后利用该模型模拟了同向气流驱动下液膜在水平预膜板表面上的流动,分析了气流速度对液膜流动形态的影响规律.研究表明,该模型可准确追踪具有大密度比的气液相界面的形态变化;气液剪切速度差会诱发两相界面出现Kelvin-Helm-holtz不稳定性现象,因而当气体速度升高时,气液剪切速度差增大,不仅液膜流动速度随之增高,且在铺展阶段液膜会产生较高振幅的波浪面并加快铺展过程,自由表面波动增强.     

可压缩气固混合层中离散相与连续相的相互作用研究

尽管已有许多文献采用数值模拟方法研究两相流问题,但主要是集中不可压流动方面.本文采用Eul-er-Lagrange颗粒-轨道双向耦合模型对时间模式下含有固粒的二维可压缩混合层流动进行了研究.气相流场采用非定常全Navier-Stokes方程描述,并应用具有空间三阶精度的WNND(Weighted Non-Oscillatory, Contai-ning No Free Parameters and Dissipative)格式进行数值高散.固相方程采用二阶单边三点差分离散.在考虑流场对固粒作用的同时,也计及颗粒对流场的反作用.主要研究混合层大尺度涡对颗粒扩散特性的影响及颗粒对流场结构的影响问题.在对流马赫数为0.5时,研究不同Stokes数颗粒在连续流场中的扩散特性,而在对流马赫数为0.8时研究了不同Stokes数颗粒对流场小激波结构的影响.     

压力梯度对边界层两相流动稳定性的影响

分析了有压力梯度的边界层两相流动稳定性,推导出类似于Saffman理论的修正的稳定性方程,数值计算采用高精度的谱方法。结果说明,压力梯度对边界层两相流动稳定性有显著的影响,顺压梯度增强流动稳定性,而逆压梯度则促进流动失稳。在不同的压力梯度和浓度下,Stokes数对流动稳定性的影响是一致的,存在一个临界Stokes数,小Stokes数促进流动失稳,而大Stokes数则提高临界雷诺数,抑制流动失稳的最佳Stokes数为10的量级。     

喷雾二相流动的初步分析与实验

Quandt、Muench等曾初步考虑了应用二相流动的一般理论来分析喷雾推进系统的性能。他们的工作与实验相结合,表明在船速高于75节时,喷雾推进可能优于螺旋桨推进和喷水推进。 本文就喷雾状二相喷管流动的基本方程进行了研讨。并和实验数据进行了比较。 一、基本假定和流动方程 假定气-液二相流动过程为一维定常流动,气相为均匀连续介质,且满足理想气体状     

旋转直管内气固或液固两相流动浓度分布

本文针对旋转直管内气固或液固两相流动，建立了经过适当简经处理的颗粒无量纲运动方程，得到了稳定状态下直管内颗粒的浓度分布公式，从而为气固两相流风机或液固两相流泵的颗粒浓度分布研究打下了基础。     

高超声速层流尾迹的数值模拟

本文利用无波动、无自由参数、耗散的差分格式(NND格式),通过求解NS方程,数值模拟了高超声速层流尾迹的流动,清晰地给出了主激波、拐角膨胀波、迹激波及自由剪切层,所得流场物理量的分布与实验结果甚为一致。计算发现了底部迴流区由起始向定常的发展中,在瞬时流线图上经历了极限环形成、胀大、缩小、再胀大最后消失的演变过程。     

Gas-particle two-phase turbulent flow in a vertical duct

Two-phase gas-phase turbulent flows at various loadings between the two vertical parallel plates are analyzed. A thermodynamically consistent turbulent two-phase flow model that accounts for the phase fluctuation energy transport and interaction is used. The governing equation of the gas-phase is upgraded to a two-equation low Reynolds number turbulence closure model that can be integrated directly to the wall. A no-slip boundary condition for the gas-phase and slip-boundary condition for the particulate phase are used. The computational model is first applied to dilute gas-particle turbulent flow between two parallel vertical walls. The predicted mean velocity and turbulence intensity profiles are compared with the experimental data of Tsuji et al. (1984) for vertical pipe flows, and good agreement is observed. Examples of additional flow properties such as the phasic fluctuation energy, phasic fluctuation energy production and dissipation, as well as interaction momentum and energy supply terms are also presented and discussed.

Applications to the relatively dense gas-particle turbulent flows in a vertical channel are also studied. The model predictions are compared with the experimental data of Miller & Gidaspow and reasonable agreement is observed. It is shown that flow behavior is strongly affected by the phasic fluctuation energy, and the momentum and energy transfer between the particulate and the fluid constituents.     

双圆诱导不可压缩势流的几个精确解

Lagally(1929)给出了两个静止圆在均匀来流中的扰动速度势精确解。本文通过使用 Apollonius 保角变换,将两圆外部区域变为圆环内部区域,并在变换平面上用 Fo-urier 级数方法求解 Laplace 方程,从而得出双圆扰动势流的另外三个精确解。它们分别对应于静止流体中两圆膨胀(收缩)、沿连心线相向(相背)运动和垂直连心线反向运动(错动)所诱导的速势场。     

Numerical solution of the turbulent flow of an incompressible fluid in curved planar and axially symmetrical channels

The governing equations for axially symmetric flow, where the Reynolds stresses are expressed by scalar turbulent viscosity, are the Reynolds equations. The turbulence model k, ? is used in the well-known form for fully developed turbulent flow.The numerical method, a continuation of the MAC system¹, is adapted so that even for high Reynolds cell numbers precision (δx²) can be achieved for the steady flow. Irregular cells join the rectangular network on the curved surface. Von Neumann's stability condition of the linearised numerical system is investigated. Special problems concerning the numerical solution of the turbulence model equations are stated and a special procedure is worked out to ensure that the fields k, ? do not converge to physically meaningless values. The program for the computer is universal in that the boundary problems can be assigned by input data.As an example, an axially symmetrical diffuser with an area ratio of widening 1.40 is computed. Fields of velocity and pressure at the wall as well as fields v_T and k are assessed. The results are compared with an experiment. The conclusion is that this method is suitable for the problems mentioned in this study as well as for unsteady flow.     

雷诺数对三角翼绕流的影响

应用染色液流动显示技术研究了雷诺数对６０°尖前缘三角翼前缘涡破裂位置、背风面流动结构等的影响，并详细分析了背风面流动随攻角的变化．     

Theoretical analysis of using airflow to purge residual water in an inclined pipe

ＩｎｔｒｏｄｕｃｔｉｏｎＰｕｒｇｉｎｇｒｅｓｉｄｕａｌｌｉｑｕｉｄｉｎａＵ＿ｓｈａｐｅｄｐｉｐｅｌｉｎｅｉｓｃｈａｌｌｅｎｇｅｆｏｒｃｈｅｍｉｃａｌａｎｄｔｒａｎｓｐｏｒｔａｔｉｏｎｉｎｄｕｓｔｒｙ[1].Ｍａｎｙｔｅｃｈｎｉｑｕｅｓｈａｖｅｂｅｅｎｅｍｐｌｏｙｅｄｔｏｐｕｒｇｅｏｕｔｒｅｓｉｄｕａｌｗａｔｅｒ[2 ],ｂｕｔｎｏｎｅｈａｖｅｂｅｅｎｃｏｍｐｌｅｔｅｌｙｓａｔｉｓｆａｃｔｏｒｙ .Ｏｎｅｐｏｓｓｉｂｌｅｍｅｔｈｏｄｆｏｒｄｏｉｎｇｔｈｉｓｈａｓｂｅｅｎｔｏｂｌｏｗａｇｒｅａｔｑｕａｎｔｉｔ…     

Sidewall flow stability in the MHD channel flows

The velocity distribution between two sidewalls is M-shaped for the MHD channel, flows with rectangular cross section and thin conducting walls in a strong transverse magnetic field. Assume that the dimensionless numbersR _m ?1,M, N? 1, and σ* and that the distance between two perpendicular walls is very long in comparison with the distance between two sidewalls. First, the equation for steady flow is established, and the solution of M-shaped velocity distribution is given. Then, an equation for stability of small disturbances is derived based on the velocity distribution obtained. Finally, it is proved that the stability equation for sidewall flow can be transformed into the famous Orr-Sommerfeld equation, in addition, the following theorems are also proved, namely, the analogy theorem, the generalized Rayleigh's theorem, the generalized Fjørtoft's theorem and the generalized Joseph's theorems.     

Solution of the Incompressible Unsteady Navier–Stokes Equations Only in Terms of the Velocity Components

A method which uses only the velocity components as primitive variables is described for solution of the incompressible unsteady Navier–Stokes equations. The method involves the multiplication of the primitive variable-based Navier–Stokes equations with the unit normal vector of finite volume elements and the integration of the resulting equations along the boundaries of four-node quadrilateral finite volume elements. Therefore, the pressure term is eliminated from the governing equations and any difficulty associated with pressure or vorticity boundary conditions is avoided. The equations are discretized on four-node quadrilateral finite volume elements by using the second-order-accurate central finite differences with the mid-point integral rule in space and the first-order-accurate backward finite differences in time. The resulting system of algebraic equations is solved in coupled form using a direct solver. As a test case, an impulsively accelerated lid-driven cavity flow in a square enclosure is solved in order to verify the accuracy of the present method.     

A collocated finite volume method for predicting flows at all speeds

An existing two-dimensional method for the prediction of steady-state incompressible flows in complex geometry is extended to treat also compressible flows at all speeds. The primary variables are the Cartesian velocity components, pressure and temperature. Density is linked to pressure via an equation of state. The influence of pressure on density in the case of compressible flows is implicitly incorporated into the extended SIMPLE algorithm, which in the limit of incompressible flow reduces to its well-known form. Special attention is paid to the numerical treatment of boundary conditions. The method is verified on a number of test cases (inviscid and viscous flows), and both the results and convergence properties compare favourably with other numerical results available in the literature.     

弯曲渠道中的振荡流和二次流

本文以弯道理想不可压缩流动的流线和等势线为曲线坐标系给出了弯道内不可压缩粘性振荡流动的一阶解以及边界层流型下的二阶解,以直角弯道为例对弯道內的粘性振荡流动以及二次流的影响进行了讨论。