Research on stability of moving jet containing dense suspended solid particles

Ｎｏｍｅｎｃｌａｔｕｒｅαｊ:ｐｈａｓｅｊｖｏｌｕｍｅｆｒａｃｔｉｏｎ;ｕｊ:ｐｈａｓｅｊｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;Ｕ :ｍｅａｎｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;Ｕ :ｍｅａｎａｘｉａｌｖｅｌｏｃｉｔｙ;Ｕｏ:ｍｅａｎｘ＿ａｘｉａｌｖｅｌｏｃｉｔｙｏｆｔｈｅｅｘｉｔ;Ｕ∞ :ｖｅｌｏｃｉｔｙｏｆｊｅｔｔｉｎｇｄｅｖｉｃｅ;ΔＵ :ｊｕｍｐｉｎｇｖｅｌｏｃｉｔｙｏｆｍｏｖｉｎｇｊｅｔ,ｃｈａｒａｃｔｅｒｉｓｔｉｃｖｅｌｏｃｉｔｙｏｆｆｌｏｗ＿ｆｉｅｌｄ ,ΔＵ =Ｕｏ-Ｕ∞ ;ｐ:ｐｒｅｓｓｕｒｅ;Ｐ :ｅｉｇｅｎｐｒｅｓ…     

Research on stability of interface of jet containing suspended solid particles

Ｎｏｍｅｎｃｌａｔｕｒｅ　αｓ:ｐａｒｔｉｃｌｅｖｏｌｕｍｅｆｒａｃｔｉｏｎ;　ｕｇ:ｇａｓｉｎｓｔａｎｔａｎｅｏｕｓｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;　ｕｓ:ｐａｒｔｉｃｌｅｉｎｓｔａｎｔａｎｅｏｕｓｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;　Ｕ :ｍｅａｎａｘｉａｌｖｅｌｏｃｉｔｙｏｆｔｗｏ＿ｐｈａｓｅｊｅｔ;　Ｕ0 :ｘ＿ａｘｉａｌｖｅｌｏｃｉｔｙｏｆｔｈｅｅｘｉｔ,ｃｈａｒａｃｔｅｒｉｓｔｉｃｖｅ ｌｏｃｉｔｙｏｆｆｌｏｗ＿ｆｉｅｌｄ;　ｐ:ｉｎｓｔａｎｔａｎｅｏｕｓｆｌｏｗ＿ｆｉｅｌｄｐｒｅｓｓｕｒｅ;　 ｐ:ｍｅａｎ…     

Numerical simulation of supersonic reacting mixing layer

In this paper, the supersonic chemically reacting mixing layer is simulated with the third order ENN scheme, based on the Navier-Stokes equations, containing transport equations of all species. The numerical results show that the thickness of mixing layer increases gradually along the flow direction, and that the Kelvin-Helmholtz instabilities may not exist in mixing layer flows. The project supported by the National Natural Science Foundation of China     

含悬浮固粒的旋转射流剪切层稳定性研究

本文在理想不可压旋转圆射流的运动方程中添加了固粒作用项，由此推得了时间增长率的表达式，进而得关于含悬浮固粒放置射流稳定性的修正瑞利稳定性准则，求出了不同固粒质量密度固－气脉动速度比值，固气脉动速度相位差及Ｓｔｏｋｅｓ数情况下旋转射流场的增长率与径向空间波数的关系曲线，在比较这些曲线的基础上，给出了关于固粒属性对旋转射流场稳定性影响的几个重要结论为控制旋转射流场和后续发展提供了依据。     

The distribution of solid particles suspended in a turbulent flow: A stochastic approach

Basic fluid mechanics and stochastic theories are applied to show that the concentration distribution of suspended solid particles in a direction normal to the mean streamlines of a two-dimensional turbulent flow is greatly influenced by the lift force exerted on them in the vicinity of the wall. Analytic solution shows that, when the direction of the mean flow is horizontal, the probability density functionp (y, t) for random displacements of the particles will have a maximum value at a point from the wall where the perpendicular component of the lift force precisely balances particle gravity. Interpretation of experimental observations is presented using this theory.     

超音速混合层稳定性分析及增强混合的研究

利用流动稳定性提高超音速混合层的混合效率,对于提高超音速流的高效混合是一个有效途径。研究结果表明,有展向曲率的三维混合层中,三维扰动的增长率很大,且法向的掺混能力也较强,可以有效地增强混合。对于高马赫数来流的超音速混合层,这一特性依然存在,这将有利于提高高超音速混合层的混合能力。     

Wavelet analysis of coherent structures in a three-dimensional mixing layer

Wavelet analysis is applied to the results obtained by the direct numerical simulation of a three-dimensional (3D) mixing layer in order to investigate coherent structures in dimension of scale. First, 3D orthonormal wavelet bases are constructed, and the corresponding decomposition algorithm is developed. Then the Navier-Stokes equations are transformed into the wavelet space and the architecture for multi-scale analysis is established. From this architecture, the coarse field images in different scales are obtained and some local statistical quantities are calculated. The results show that, with the development of a mixing layer, the energy spectrum densities for different wavenumbers increase and the energy is transferred from the average flow to vortex structures in different scales. Due to the non-linear interactions between different scales, cascade processes of energy are very complex. Because vortices always roll and pair at special areas, for a definite scale, the energy is obtained from other scales at some areas while it is transferred to other scales at other areas. In addition, energy dissipation and transfer always occur where an intense interaction between vortices exists. The project supported by the Research Fund for the Doctoral Program of Higher Education and the National Natural Science Foundation for Outstanding Youth of China (19925210)     

空间模式下可压缩气固两相混合层流动特性研究

采用Euler-Lagrange颗粒-轨道双向耦合模型对空间模式下含有固粒的二维可压缩混合层流场进行了研究。气相流场采用具有空间三阶精度的WNND格式进行数值模拟,固相方程采用单边三点差分离散。在考虑流场对固粒作用的同时,也计及了固粒对流场的反作用。在对流马赫数为0.5时,研究了颗粒相对密度、颗粒尺寸、Stokes数等因素对粒子运动和流场结构的影响。研究结果表明:在可压缩空间模式混合层中,固粒的Stokes数仍然是主要影响参数之一;相同Stokes数下不同密度的固粒对流场的干扰不同,轻固粒对流场的干扰明显要小。     

Mode decomposition of a noise suppressed mixing layer

Noise is generated in atwo-dimensional mixing layer due to the growing of instability waves and vortex pairings. The adjoint-based control methodology has shown to be arobust tool to suppress noise radiation. The mode decomposition algorithms such as the compressible versionof proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are employed toanalyze thespatial/spatial-temporal coherent structures for a consecutive data sets of the controlled mixing layer and itsuncontrolled counterpart. The analyses of POD indicate that the y-direction body forcecontrol mainly modify themost energetic spatialstructures, and increase the uniformity of the flow. The analyses of DMD show us prevalent frequencies andcorresponding mode structures, and the stability characteristics of each mode can be obtained fromDMD-spectrum. The spectral signatures illustrate that a lot of neutral/slightly damping modesemerging in uncontrolled flow within the frequency range (ω < 0.4) are suppressed due to control, relevant spatial-temporal structures are also varied, which iscoincident with the change of far-field noise spectra. From the view of mode decomposition, the action of control redistribute the energy forfrequency components of ω < 0.4 by weakening nonlinearities and regularizing corresponding dynamicstructures in streamwise direction, and thus suppress the noise radiation. Moreover, the POD- and DMD-analysis in this studydemonstrate that DMD can serve as an important supplement for POD in analyzing a time-resolved physicalprocess.     

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

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

Research on particle dispersion in a plane mixing layer with coherent structures

The numerical simulation with two-way coupling was performed in a liquid -particle mixing layer and the corresponding experiment study was made. In the process of vortex rolling up and vortices pairing, the particles with different St number have a very different pattern of dispersion. The mean velocity of particle with St = 1 is higher than that of the fluid phase on the low-speed side, and lower than that of the fluid phase on the high-speed side. The RMS of particle approaches that of the fluid phase with decreasing particle St number. The RMS in the transverse direction is smaller than that in the streamwise direction. The velocity fluctuation correlation of particle is smaller than the Reynolds shear stress, the “overshoot“ phenomenon that the velocity fluctuation correlation of particle is larger than the Reynolds shear stress does not appear. The larger the St number of particle is, the wider the range of the particle dispersion will be. The computed results are in agreement with the experimental ones.     

The distribution of particles in a shock-induced boundary layer of a dusty gas over a solid surface

The laminar boundary layer behind a constant-speed shock wave moving through a dusty gas along a solid surface is studied. The Saffman lift force acting on a spherical particle in a gas boundary layer is taken into account. A method for calculating the density profile of dispersed phase near the wall is proposed and some numerical results are given. It is shown that behind the shock wave, there exists a curved thin layer where the density of particles is many times higher than the original one. This dust collection effect may be of essential importance to the problem of dust explosion in industry.     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.     

LOAD CAPACITY OF OILFILM——LUBRICATED BALL-BEARING WITH SPIRAL GROOVES

In order to improve load capacity and stability of a fluid moving-pressure bearing, at present, measure of cutting grooves on a turning-shaft is generally adopted.The performances of spiral grooves are better in various groovetypes. Therefore, a theoretical analysis and an experimental investigation of the bearing with spiral grooves attract general attention. But the results of most articles published in china are numerical solutions of the computer. Using the methods of parameter perturbation in the paper, an approaching analytical solution of Reynolds equation of the moving-pressure oilfilm-lubricated ball-bearing with spiral grooves has been obtained. We have combined the numerical solutions of the computer, and the influence of the bearing parameters upon the load capacity is discussed; the optimum values of the groove parameter are obtained. The theoretical conclusion of this paper is in good agreement with the experimental data obtained by 651 Institute.     

用格子Boltzmann方法数值模拟混合层中旋涡的合并过程

用格子Boltzmann方法计算混合层中的流动问题。在流场的入口处加不同频率、振幅和相位的小扰动,观察混合层中旋涡的演进机理,模拟二维混合层中旋涡合并现象。在基本扰动波的基础上,又加入频率为基本波频率一半的亚谐波,得到了两个涡合并的计算结果,当加入的亚谐波频率为基本波频率的三分之一时,得到了三个涡合并的计算结果。这些计算结果与已有文献的结果基本一致,显示用格子Boltzmann方法模拟混合层问题是可行的。     

Mixing and segregation assessment of bi-disperse solid particles in a double paddle mixer

A double paddle blender's flow patterns and mixing mechanisms were analyzed using discrete element method (DEM) and experiments. The mixing performance of this type of the blender containing bi-disperse particles has been rarely studied in the literature. Plackett-Burman design of experiments (DoE) methodology was used to calibrate the DEM input parameters. Subsequently, the impact of the particle number ratio, vessel fill level, and paddle rotational speed on mixing performance was investigated using the calibrated DEM model. The mixing performance was assessed using relative standard deviation and segregation intensity. Mixing performance was significantly affected by the paddle rotational speed and particle number ratio. Moreover, the Peclet number and diffusivity coefficient were used to evaluate the mixing mechanism in the blender. Results revealed that the diffusion was the predominant mixing mechanism, and the best mixing performance was observed when the diffusivity coefficients of 3 mm and 5 mm particles were almost equal.     

Eulerian-Lagranigan simulation of aerosol evolution in turbulent mixing layer

The formation and evolution of aerosol in turbulent flows are ubiquitous in both industrial processes and nature. The intricate interaction of turbulent mixing and aerosol evolution in a canonical turbulent mixing layer was investigated by a direct numerical simulation (DNS) in a recent study (Zhou, K., Attili, A., Alshaarawi, A., and Bisetti, F. Simulation of aerosol nucleation and growth in a turbulent mixing layer. Physics of Fluids, 26, 065106 (2014)). In this work, Monte Carlo (MC) simulation of aerosol evolution is carried out along Lagrangian trajectories obtained in the previous simulation, in order to quantify the error of the moment method used in the previous simulation. Moreover, the particle size distribution (PSD), not available in the previous works, is also investigated. Along a fluid parcel moving through the turbulent flow, temperature and vapor concentration exhibit complex fluctuations, triggering complicate aerosol processes and rendering complex PSD. However, the mean PSD is found to be bi-modal in most of the mixing layer except that a tri-modal distribution is found in the turbulent transition region. The simulated PSDs agree with the experiment observations available in the literature. A different explanation on the formation of such PSDs is provided.     

Spanwise domain effects on the evolution of the plane turbulent mixing layer

Large Eddy Simulation is used to simulate a series of plane mixing layers. The influence of the spanwise domain on the development of the mixing layer, and the evolution of the coherent structures, are considered. The mixing layers originate from laminar conditions, and an idealised inflow condition is found to produce accurate flow predictions when the spanwise computational domain extent is sufficient to avoid confinement effects. Spanwise domain confinement of the flow occurs when the ratio of spanwise domain extent to local momentum thickness reaches a value of ten. Flow confinement results in changes to both the growth mechanism of the turbulent coherent structures, and the nature of the interactions that occur between them. The results demonstrate that simulations of the two-dimensional mixing layer flow requires a three-dimensional computational domain in order that the flow will evolve in a manner that is free from restraints imposed by the spanwise domain.     

Transport of particles in an atmospheric turbulent boundary layer

A program incorporating the parallel code of large eddy simulation (LES) and particle transportation model is developed to simulate the motion of particles in an atmospheric turbulent boundary layer (ATBL). A model of particles of 100-micrometer order coupling with large scale ATBL is proposed. Two typical cases are studied, one focuses on the evolution of particle profile in the ATBL and the landing displacement of particles, whereas the other on the motion of particle stream.The English text was polished by Yunming Chen.     

Effect of coarse particles on the heat transfer in a particle-laden turbulent boundary layer

The temperature distribution in particle-laden turbulent flow, in a flume, was investigated both by DNS and experimentally. Simulations were performed at Re=171 and Pr=5.4 in order to study the interaction between the particle motion and flow turbulence. Two-way coupling was used to obtain various turbulence statistics, the grid resolution was sufficiently fine to resolve all essential turbulent scales. The effect of particle diameter on momentum, heat transfer and particle deposition was considered. The details of particle-turbulence interaction depend on the particle Stokes number and the particle Reynolds number.

The spatial structures of instantaneous flow and temperature fields were visualized. Low frequency small oscillations of deposited particles were observed. It was found that these small deviations from the initial position, caused strong changes in the instantaneous temperature field near the particle.

The experiments provided details of the temperature field on the heated wall close to the particle. In the front of the particle, a sharp increase in heat transfer coefficient was observed. The experimental results agree well with the computational predictions.