双流体颗粒-壁面碰撞模型用于旋流流动

对旋流两相流动,用考虑壁面粗糙度的双流体颗粒-壁面碰撞模型,结合二阶矩两相湍流模型和不同颗粒相边界条件,进行了数值模拟。结果表明,由于考虑了各方向雷诺应力之间的相互转化,雷诺应力从平均运动中得到能量,以及壁面对颗粒运动的衰减作用等因素,该颗粒壁面碰撞模型给出的模拟结果与实验符合较好,而广泛应用的零梯度颗粒壁面边界条件给出的模拟结果最差。     

两相可压缩边界层的颗粒运动特性分析

本文采用高精度直接数值模拟方法,研究了马赫数为2的可压缩湍流边界层中,两种不同粒径颗粒的运动特性.研究发现近壁面的颗粒流向平均速度要小于流体;颗粒具有法向指向壁面的平均速度;近壁面颗粒的流向雷诺正应力要略大于流体:其中大颗粒具有更加显著的上述趋势.小颗粒的法向和展向雷诺正应力基本等于流体,而大颗粒则都要远小于流体.条件平均统计显示,大颗粒在低速区域的富集导致了这一现象.     

下降管中稠密两相湍流的数值模拟

本文将统一二阶矩两相湍流模型和颗粒动力学理论结合,推导并封闭了稠密两相流考虑颗粒间碰撞的统一二阶矩两相湍流模型。该模型用颗粒动力学理论模拟颗粒之间的碰撞,用各向异性的统一二阶矩模型考虑气相和颗粒相的湍流脉动,并用输运方程描述气固两相湍流之间的相互作用。最后用该模型对下降管中的气粒两相流动进行了模拟,模拟所得的沿径向颗粒浓度分布和速度分布与实验吻合较好,预报结果也反应出了颗粒雷诺应力的各向异性。     

稠密气固流动二阶矩摩擦应力模型和模拟

基于统一二阶矩两相湍流应力模型和稠密颗粒动力学理论,建立了欧拉-欧拉双流体二阶矩颗粒摩擦应力模型。模型充分反映各向异性的气固两相相间雷诺应力相互作用,并引入有效颗粒弹性恢复系数,考虑了因颗粒表面不光滑产生的摩擦力对湍流流动结构和颗粒弥散特性的影响,对于下降管的计算结果与实验吻合较好。因颗粒摩擦产生能量耗散降低了颗粒温度和导致颗粒雷诺应力再分布。在入口和出口区域因颗粒碰撞频率较高而产生的能量耗散对流动结构影响明显。     

两相壁面射流PDF方程有限分析方法及数值模拟

提出求解位置-速度相空间中高维两相流PDF(probability density function)方程的有限分析方法,将位置-速度相空间颗粒PDF方程约化到速度空间,并解析求解,颗粒的位置PDF用轨道方法求解.对壁面射流两相流动进行数值模拟,并与颗粒雷诺应力轨道方法进行比较计算,结果优于颗粒雷诺应力轨道方法.     

突扩湍流气粒流动两相脉动速度关联和颗粒速度概率密度分布的实验研究

本文用PDPA对突扩气粒两相流动中的两相脉动速度关联进行了测量,得到了30和50微米颗粒的两相脉动速度关联,其空间分布和气体及颗粒的雷诺应力分布一致．同时发现,气粒两相速度概率密度分布并不都是高斯分布。     

颗粒碰撞率的实验研究

本文采用一种新的研究思路,以高速摄像作为手段,结合各种新颖的PTV算法,在毫秒级时间尺度和毫米级空间尺度上(颗粒尺寸)来对颗粒碰撞率做实验研究,对长期以来一直沿用的描述颗粒碰撞率的关系进行检验与修正。通过对实际颗粒碰撞率与颗粒浓度、平均相对速度的关系进行测量研究,发现理论关系与实验存在较大的差别,理论关系较大地高估了碰撞率。指出了需要进一步深入和完善描述颗粒碰撞率的实验规律,这些规律可直接应用于气固两相流动的数值计算与理论分析中。     

气固两相圆柱绕流的直接数值模拟

本文运用谱元方法对气固两相圆柱绕流进行了直接数值模拟。在获得高精度计算流场信息基础上,进行颗粒扩散运动的研究。通过研究颗粒大小对颗粒扩散运动的影响来刻画气固两相圆柱绕流的物理特征。研究时发现,由于圆柱尾流独特的涡结构特征,颗粒在流场中心轴线附近区域的扩散机理主要是决定于吸力作用,该作用会驱使小颗粒卷吸进入圆柱尾部近壁区,甚至与圆柱后壁发生碰撞,这种扩散机理与混合层中颗粒扩散机理明显不同。     

一种新的壁面函数

壁面函数是处理近壁区湍流的一种工程方法。常用的一种壁面函数是以混合长度模型为基础的,求出壁面应力后采用雷诺比拟求壁面热流。这种壁面函数在求分离点或再附点附近的壁面热流时遇到了困难,因为实验表明此时壁面应力极小,而壁面热流达到极大,二者不成比例。     

风沙运动的DPM数值模拟

采用离散颗粒模型对风沙气固两相流动进行了数值计算．在该模型中,采用体平均的Navier-Stokes方程来描述气相的运动,对离散颗粒采用拉格朗日方法模拟,求解颗粒运动方程．采用硬球模型描述颗粒问碰撞作用．计算结果表明,沙粒平均水平速度廓线在0．02 m以上高度可以表示为按对数函数或幂函数规律增加,在0．02 m以下则发生偏离;在大于0．02 m的高度,输沙通量随高度按指数规律衰减,而在地面附近由于颗粒蠕移的影响发生偏离,这与已有文献结果一致。本文的模拟有助于风沙运动规律的研究与掌握。     

Surface roughness effects on the turbulence statistics in a low Reynolds number channel flow

A high-resolution particle image velocimetry was used to characterize a low Reynolds number turbulent flow in a channel. Experiments were conducted over a sand grain-coated surface of large relative roughness, and the results were compared with measurements over a smooth surface. The roughness perturbation significantly modified the outer layer. Even though the streamwise Reynolds stress shows less sensitivity in the outer layer to the boundary condition, significant enhancements were observed in the wall-normal Reynolds stress and the Reynolds shear stress. These modifications were considered as footprints of the larger-scale eddies transporting intense wall-normal motions away from the rough wall. A quadrant decomposition shows that strong and more frequent ejections are responsible for the larger values of the mean Reynolds shear stress over the rough wall. The results also indicate that spanwise vortex cores with mean vorticity of the same sign as the mean shear are the dominant smaller-scale vortical structures over the smooth and rough walls. A linear stochastic estimation-based analysis shows that the average larger-scale structure associated with these vortices is a shear layer that strongly connects the outer layer flow to the near-wall flow. A proper orthogonal decomposition of the flow suggests that the large-scale eddy is more energetic for the rough wall, and contributes more significantly to the resolved turbulent kinetic energy and the Reynolds shear stress than the smooth wall.     

Long-time behavior of sand ripples induced by water shear flow

The formation of sand ripples under water shear flow in a narrow annular channel and the approach of the ripple pattern towards a steady state were studied experimentally. Four results are obtained: i) The mean amplitude, the average drift velocity and the mean sediment transport rate of the evolving bed shape are strongly related. A quantitative characterization of this relation is given. ii) The ripple pattern reaches a stationary state with a finite ripple amplitude and wavelength. The time needed to reach the state depends on the shear stress and may be several days. iii) The onset of ripple formation is determined by the bed shear stress, but it seems neither to depend on the grain diameter nor on the depth of the water layer. iv) The ripple amplitude, drift velocity and sediment transport in this stationary state depend on the grain size. This dependency is neither captured by the particle Reynolds number nor by the Shields parameter: an empirical scaling law is presented instead.     

Formation of ripples over a sand bed submitted to a turbulent shear flow

We investigate the process of ripple formation when a sand bed is submitted to a steady and turbulent liquid flow. The sand transport dynamics is described in terms of a simple relaxation law which accounts for the fact that the transport rate does not adapt instantaneously to its equilibrium value. The equilibrium sand flux is evaluated using a standard law based on the estimation of the flow shear stress calculated at the sand bed surface. The latter is estimated from an analytical resolution of the flow over a deformed sand bed which is based on the Jackson and Hunt calculation [J.C.R. Hunt, Quart. J. R. Met. Soc. 101, 929 (1975)]. Within this model, we investigate the stability of the sand bed and are able to derive analytical scaling laws for the wavelength and phase velocity of the most dangerous mode. In the deep flow limit, the model predicts the occurrence of a single mode of instability corresponding to the formation of ripples. Predictions of our model are compared with previous models and available experimental data.     

Effects of inclined impermeable bed on the turbulent characteristics of the flow using particle image velocimetry

In this study, the turbulent characteristics of the flow in an open channel with horizontal and inclined impermeable beds were studied experimentally using two-dimensional particle image velocimetry (PIV). The experiments were conducted in a channel of 6.5 m in length, 7.5 cm in width and 25 cm in height. The slope of the channel was S = 0 for the horizontal impermeable bed and for S = ?0.002, S = ?0.004 and S = ?0.006 for the inclined impermeable bed. Hydraulic characteristics such as distributions of velocities, turbulent intensities and Reynolds stress are investigated at a fine resolution using the PIV. Velocity is measured above the horizontal and inclined impermeable bed for the same different heights (h = 5, 7, 9, 11 and 13 cm) and for the same different discharges (Q = 0.735, 0.845 and 0.970 lt/s). Results show that the channel slope influences significantly near the impermeable bed but not near the free surface the variation of turbulent characteristics of the flow and also the alteration of the channel slope from ?0.002 to ?0.006 doesn't influence the variation of turbulent characteristics of the flow, which are the longitudinal turbulent intensity u′U_*, the vertical turbulent intensity v′/U_* and the turbulent kinetic energy. The channel slope doesn't influence the Reynolds stress.     

Wall Shear Stress and Heat Transfer of Downward Bubbly Flow at Low Flow Rates of Liquid and Gas

The effect of suppression of turbulence in a downward bubbly flow and its impact on the wall shear stress and heat transfer are discussed. Measurements were carried out for Reynolds numbers Re = 5000–10000, which were calculated from the velocity of the liquid phase and with the gas volumetric flow rate ratio β = 0–0.05. Data on the size of bubbles detaching from the edges of an array of capillaries in a liquid flow are given. The influence of the disperse phase dimensions on the wall shear stress and heat transfer is discussed. It is shown that change in the size of the dispersed phase can lead to both intensification and deterioration of heat transfer as compared with a single-phase flow at constant flow rates of liquid and gas at the channel inlet. The cause of the heat transfer deterioration is “laminarization” of the flow in the near-wall region. An analysis of the spectral power of signals is given.     

Flow characteristics of an offset jet over a surface mounted square rib

A three-dimensional offset jet flow over a surface mounted square rib was investigated using particle image velocimetry at Reynolds number of 7600. The square rib was mounted at three different locations downstream of the nozzle exit. Contour and profile plots in the symmetry plane were used to investigate the effect of the rib on the mean flow distribution and turbulence statistics. Changing the rib location enhanced the mixing and entrainment characteristics of the jet thereby leading to larger wall-normal spread rates. The quadrant analysis revealed a significant contribution of all events towards the Reynolds shear stress. The joint probability density function analysis showed that the structures within the outer shear layer of the jet differed from those of boundary layers. Results from the two-point correlation analysis indicated that the introduction of a surface mounted rib modified the flow structures within the developing region. However, there was no significant structural difference between the flow with mounted rib and the no rib case within the self-similar region. The proper orthogonal decomposition was used to examine how changes in rib location modified the dynamics of the energetic modes towards the turbulent kinetic energy and Reynolds shear stress within the developing and self-similar regions.     

Experimental study of vertical stress profiles of a confined granular bed under static and dynamic conditions

In a wet agglomeration process inside a low shear mixer, the blade function is to induce i) homogenization of the liquid sprayed on the powder surface and ii) a stress field able to transfer the mechanical energy at the particle scale. In this work we study the mechanical state of a confined powder bed through the analysis of stress distributions (by force measurements) in a rectangular cell in two cases: for a classical model powder (i.e. glass beads) and a complex powder (i.e. wheat semolina). Two types of vertical stress profiles are obtained according to the type of measurements carried out in the powder bed, either locally (at different positions in the cell) or globally (at the entire base). The global vertical stress profile follows Janssen's model and the local vertical stress profile highlights a critical length, identified as the percolation threshold of the force network, and a shielding length near the bottom, which is similar to an influence length of the side walls. In the context of wet agglomeration, the results allow to consider the role of the characteristic lengths in the mixing bowl under vertical mechanical solicitation.     

A model of Barchan dunes including lateral shear stress

Barchan dunes are found where sand availability is low and wind direction quite constant. The two dimensional shear stress of the wind field and the sand movement by saltation and avalanches over a barchan dune are simulated. The model with one dimensional shear stress is extended including surface diffusion and lateral shear stress. The resulting final shape is compared to the results of the model with a one dimensional shear stress and confirmed by comparison to measurements. We found agreement and improvements with respect to the model with one dimensional shear stress. Additionally, a characteristic edge at the center of the windward side is discovered which is also observed for big barchans. Diffusion effects reduce this effect for small dunes.     

微结构超疏水表面减阻特性数值研究

通过湍流状态下超疏水表面的数值模拟,研究超疏水表面湍流流动的速度和剪应力分布,分析其与减阻的内在关系,初步总结出微观结构与减阻率之间的基本规律.数值模拟采用非定常雷诺平均模型,气液两相流则采用VOF模型.结果表明:超疏水表面存在滑移流动和周期分布的剪应力,其不仅可以实现减阻也有可能导致增阻;微结构尺寸对减阻率有显著影响,为了尽可能增大减阻率,矩形微结构的深宽比应大于3:2,且凹槽间距应尽可能小,凹槽宽度应小于200微米.