建立两相流方程的动力论方法

本文从Boltzmann方程出发建立两相流基本方程。其中碰撞项的处理采用了类似于经典动力论在处理稠密气体输运性质时采用过的方法。动力论方法表明,在气固两相流中除通常的分压和分热流外,还存在碰撞压强和碰撞热流,它们同通常的分压和分热流有本质的差别。本文还讨论了二体碰撞假设和混沌假设在气固两相流中的适用性。     

On pressure and thermal flux in gas-mixture flows and in two-phase flows

To begin with, two different definitions of pressure, thermal flux, etc. in the diffusion model and two-fluid model are given. Then the physical interpretations of the pressure and the thermal flux are provided by introducing the momentum and energy fluxes,M and ε, through a surface dS in the flow field. The quantities defined in the diffusion model are suggested when the motion of the mixture is studied as a whole, while the quantities defined in the two-fluid model are suggested when the motion of individual species is studied. The collision pressure and thermal flux in dense gas-mixtures are also discussed in detail, i.e. their origin, their expressions in the momentum and energy equations, and their distinctions from the normal partial pressure and thermal flux. A gas-particle flow can be treated as a flow of dense gas-mixtures. The long-standing controversy whether the “inertial coupling term” should exist in the momentum equation can be clarified by the two different definitions of pressure.     

基于颗粒湍流和颗粒碰撞相互作用规律的湍流模型的研究

颗粒湍流和颗粒碰撞的相互作用规律是两相流动中的核心问题。用颗粒湍流模型和颗粒碰撞的动力论模型叠加的方法在研究两相湍流流动方面取得了一定的成效,但是还有待改进。本文基于颗粒湍流形成大尺度脉动和颗粒间碰撞引起小尺度脉动的概念,从双流体模型出发,建立了两相流动的双尺度kp-pε两相湍流模型。利用该模型对下行床和突扩室内的气固...     

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

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

扩展重要抽样法及其在平尾转轴可靠性分析中的应用

提出了扩展重要抽样法，用以计算结构系统的多个失效模式中含有不全相同基本随机变量时的系统失效概率。通过构造扩展重要抽样法的抽样密度函数，给出该方法失效概率的估计值，以及其方差和变异系数的计算公式。并将此方法用于某型飞机平尾转轴的可靠性分析，算例结果表明其优越性。     

Improved subgrid scale model for dense turbulent solid-liquid two-phase flows

The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules. Assuming that the solid-phase velocity distributions obey the Maxwell equations, the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow, the solid-particle‘s governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations. Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale (SGS) model, in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor, is proposed to model the two-phase governing equations by applying dimension analyses. Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls, the velocity and pressure fields, and the volumetric concentration are calculated. The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.     

Two-phase turbulence models for simulating dense gas-particle flows

The two-fluid model is widely adopted in simulations of dense gas-particle flows in engineering facili- ties. Present two-phase turbulence models for two-fluid modeling are isotropic. However, turbulence in actual gas-particle flows is not isotropic. Moreover, in these models the two-phase velocity correlation is closed using dimensional analysis, leading to discrepancies between the numerical results, theoretical analysis and experiments. To rectify this problem, some two-phase turbulence models were proposed by the authors and are applied to simulate dense gas-particle flows in downers, risers, and horizontal channels; Experimental results validate the simulation results. Among these models the USM-O and the two-scale USM models are shown to give a better account of both anisotropic particle turbulence and particle-particle collision using the transport equation model for the two-phase velocity correlation.     

On basic equations of turbulent swirling gas-solid flows and their application in cyclones

The basic equations of turbulent gas-solid flows are derived by using the pseudo-fluid model of particle phase with a refined two-phase turbulence model. These equations are then applied to swirling gas-particle flows for analyzing the collection efficiency in cyclone separators.     

Measurement and simulation of the two-phase velocity correlation in sudden-expansion gas-particle flows

In this paper the present authors measured the gas-particle two-phase velocity correlation in sudden expansion gas-particle flows with a phase Doppler particle anemometer(PDPA) and simulated the system behavior by using both a Reynolds-averaged Navier-Stokes(RANS) model and a large-eddy simulation(LES).The results of the measurements yield the axial and radial time-averaged velocities as well as the fluctuation velocities of gas and three particle-size groups(30 μ m,50 μ m,and 95 μ m) and the gas-particle velocity correlation for 30 μ m and 50 μ m particles.From the measurements,theoretical analysis,and simulation,it is found that the two-phase velocity correlation of sudden-expansion flows,like that of jet flows,is less than the gas and particle Reynolds stresses.What distinguishes the two-phase velocity correlations of sudden-expansion flow from those of jet and channel flows is the absence of a clear relationship between the two-phase velocity correlation and particle size in sudden-expansion flows.The measurements,theoretical analysis,and numerical simulation all lead to the above-stated conclusions.Quantitatively,the results of the LES are better than those of the RANS model.     

Evaluation of a Modified Reynolds Stress Model for Turbulent Dispersed Two-Phase Flows Including Two-Way Coupling

A modified Reynolds stress turbulence model for the pressure rate of strain can be derived for dispersed two-phase flows taking into account gas-particle interaction. The transport equations for the Reynolds stresses as well as the equation for the fluctuating pressure can be derived starting from the multiphase Navier–Stokes equations. The unknown pressure rate of strain correlation in the Reynolds stress equations is then modelled by considering the multiphase equation for the fluctuating pressure. This leads to a multiphase pressure rate of strain model. The extra particle interaction source terms occurring in the model for the pressure rate of strain can be constructed easily, with no noticeable extra computational cost. Eulerian–Lagrangian simulation results of a turbulent dispersed two-phase jet are presented to show the differences in results with and without the new two-way coupling terms.     

Lie symmetry solutions for two-phase mass flows

We apply Lie symmetry method to a set of non-linear partial differential equations, which describes a two-phase rapid gravity mass flow as a mixture of solid particles and viscous fluid down a slope (Pudasaini, J. Geophys. Res. 117 (2012) F03010, 28 pp [1]). In order to systematically explore the mathematical structure and underlying physics of the two-phase mixture flow, we generate several similarity forms in general form and construct self-similar solutions. Our analysis generalizes the results, obtained by applying the Lie symmetry method to relatively simple single-phase pressure-driven gravity mass flows, to the two-phase mass flows that include several dominant driving forces and strong phase-interactions. Analytical and numerical solutions are presented for the symmetry-reduced homogeneous and non-homogeneous systems of equations. Analytical and numerical results show that the new models presented here can adequately describe the dynamics of two-phase debris flows, and produce observable phenomena that are consistent with the physics of the flow. The solutions are strongly dependent on the choice of the symmetry-reduced model, as characterized by the group parameters, and the physical parameters of the flows. These solutions reveal strong non-linear and distinct dynamic evolutions, and phase-interactions between the solid and fluid phases, namely the phase-heights and phase-velocities.     

A two-scale second-order moment two-phase turbulence model for simulating dense gas-particle flows

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the single-scale second-order moment two-phase turbulence model superposed with a particle collision model (USM-θ model) in most regions. The project supported by the Special Funds for Major State Basic Research, China (G-1999-0222-08), and the Postdoctoral Science Foundation (2004036239) The English text was polished by Keren Wang     

非对称射流形成的几何理论(二)

将平板非对称正碰撞形成射流满足的定常方程组封闭条件推广到非正碰情况，提出了平面二维条件下非对称射流形成的一般理论，给出了射流形成的临界条件和出流射流的理论公式。对斜碰撞情形下非对称射流的出流方向、厚度随碰撞角和来流初始位形的变化规律进行了理论预测和数值模拟，理论解与二维欧拉程序的仿真结果符合较好。     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 2. Flows with heat transfer and combustion

The objective of this part of the paper is to summarize the information concerning the authors' works in the field of simulation of two-phase gas-particle turbulent flows with heat transfer and combustion. A kinetic equation had been derived for the probability density function (PDF) of the particle velocity, temperature, and mass distributions in turbulent flows. This PDF equation is used for the construction of the governing conservation equations of mass, momentum, and heat transfer in the dispersed particle phase.The numerical scheme incorporates two-phase fluid dynamics, convective and radiative heat transfer, and combustion. The proposed models have been applied to the calculation of various particle-laden turbulent flows in jets, combustion and gasification chambers, and furnaces.     

激波主导流动下壁板的热气动弹性稳定性理论分析

针对激波主导流动下弹性壁板的热气动弹性稳定性分析问题,建立了基于当地活塞流理论的分析模型,并用数值仿真方法来验证其正确性. 首先基于Hamilton原理和Von-Karman大变形理论,建立壁板的热气动弹性运动方程,其中假设壁板受热后温度均匀分布,激波前后区域的气动力模型采用当地一阶活塞流理论;利用Galerkin方法将具有连续参数系统的偏微分颤振方程离散为有限个自由度的常微分方程;基于李雅普诺夫间接法将非线性颤振方程组在平衡位置处进行线化,再用Routh-Hurwits判据来判断线性系统的稳定性,从而来推论出非线性颤振系统的气动弹性稳定性. 在时域中采用龙格--库塔法对非线性颤振方程进行数值积分,得到壁板非线性颤振响应的时间历程,与理论分析结果进行对比. 研究结果表明,壁板受到斜激波冲击时,更容易发生颤振失稳,并且激波强度越大,极限环幅值和频率越大;激波主导流场中的壁板失稳边界不同于传统单纯超声速气流中壁板颤振的失稳边界;只有在斜激波前后不同的动压值都满足颤振稳定性边界的条件下,壁板才可能保持其气动弹性稳定性.      

A two-fluid model for critical vapour-liquid flow

One-dimensional two-fluid equations are used to calculate the mass flux of initially saturated or subcooled water discharging from a pipe in critical flow. The model allows in a general way for thermal non-equilibrium between the liquid and vapour bubbles, and for interphase relative motion. The theory is shown to be in good agreement with the measured critical flow-rates of pressurised water over a wide pressure range for a single choice of parameters characterising (i) the density of nucleation sites in the liquid, (ii) the liquid superheat required to cause bubble nucleation. Predictions are made of the critical mass flux of initially saturated water in pipes of the range of sizes of interest in water-reactor blowdown safety analysis. Results indicate that for pipes up to ten diameters in length flows will be significantly higher than values obtained from conventional homogeneous thermal equilibrium flow theory.     

Modeling the evolution of droplet size distribution in two-phase flows

A theoretical model is developed in the present study to simulate droplet motion and the evolution of droplet size distribution (DSD) in two-phase air/dispersed water spray flows. The model takes into account several processes which influence DSD and droplet trajectory: droplet collision and coalescence, evaporation and cooling, gravitational settling, and turbulent dispersion of dispersed phase. The DSDs determined by the model at different locations in a two-phase flow are evaluated by comparing them to experimental observations obtained in an icing wind tunnel. The satisfactory coincidence between simulation and experimental results proves that the model is reliable when modeling two-phase flows under icing conditions. The model is applied for two particular examples in which the modification of DSD is calculated in two-phase flows under conditions describing in-cloud icing and freezing drizzle.     

Simulation of sudden-expansion and swirling gas-particle flows using a two-fluid particle-wall collision model with consideration of the wall roughness

A two-fluid particle-wall collision model with consideration of wall roughness is proposed. It takes into account the effects of the friction, restitution and in particular the wall roughness, and hence the redistribution of Reynolds stress in different directions, the absorption of turbulent energy from the mean motion and the attenuation of particle motion by the wall. The proposed model is used to simulate sudden-expansion and swirling gas-particle flows and is validated by comparing with expermental results. The results show that the proposed model gives better results than those obtained by the presently used zero-gradient condition. Hence, it is suggested that the proposed model should be used as the wall boundary condition for the particle phase in place of the presently used boundary condition. The project supported by the Special Funds for the Major State Basic Research, China (G-1999-0222-08)     

Effect of drag force on backward-facing step gas-particle turbulent flows

An improved drag force coefficient of gas-particle interaction based on the traditional Wen’s 1966 model is proposed. In this model, a two-stage continuous function is used to correct the discontinuous switch when porosity less than 0.2. Using this proposed correlation and the Wen’s 1966 model, a gas-particle kinetic energy and particle temperature model is developed to predict the hydrodynamic characteristics in backward-facing step gas-particle two-phase turbulent flows. Numerically results showed that they are in good agreement with experiment measurements and presented model are better due to a improvement of momentum transport between gas and particle phases. Particle dispersions take on the distinctively anisotropic behaviors at every directions and gas phase fluctuation velocity are about twice larger than particle phases. Particle phase has a unique transportation mechanism and completely different from the gas phase due to different density. Furthermore, the correlation values of axial–axial gas-particle are always greater than the radial–radial values at fully flow regions. The gas-particle two-phase interactions will make influence on two-phase turbulent flow behaviors.