循环流化床内颗粒团聚物流动特性的DSMC-LES模拟

采用直接模拟蒙特卡罗方法(DSMC)模拟颗粒间的碰撞,采用考虑颗粒脉动流动对气相湍流流动影响的大涡模拟(LES)研究气相湍流.单颗粒运动满足牛顿第二定律,颗粒相和气相相间作用的双向耦合由牛顿第三定律确定。数值模拟垂直管内气固两相上升流动,对管内气相速度和颗粒相速度、浓度以及聚团的流动特性进行分析.研究平均单个颗粒团聚物的存在时间、颗粒团聚物的时间份额和颗粒团聚物的生成频率分布特性,模拟结果与Manyele等(2002)和Sharma等(2000)试验结果吻合.     

LES-DSMC方法研究超细颗粒气固两相流动过程

采用考虑颗粒脉动流动对气相湍流流动影响的大涡模拟(LES)研究气相湍流,采用直接模拟蒙特卡罗方法(DSMC)模拟颗粒间的碰撞。单颗粒运动满足牛顿第二定律,颗粒相和气相相间作用的双向耦合由牛顿第三定律确定,考虑超细颗粒间的van der Waals作用力。数值模拟垂直管内超细颗粒气固两相流动,对颗粒相速度、浓度以及团聚物流动过程进行分析。     

三维射流中颗粒与流体的双向耦合作用

采用基于欧拉一拉格朗日的双向耦合模型对三维气固两相平面射流中颗粒与流体的双向耦合作用进行了直接数值模拟.在考虑颗粒相的反作用后,气相运动采用直接耦合求解,计算颗粒场时,选取Stokes数为0.1的较小颗粒,采用Lagrangian方法跟踪其运动.重点考察了颗粒相与流体相之间的相互作用,分析了不同固相载率的颗粒对流场特性以及对自身扩散的影响.模拟结果表明由于颗粒的影响,在射流入口处流场最初生成的两个大涡沿横向被拉伸,而在射流下游区域,涡结构则沿流向被拉伸;在射流的下游,颗粒降低了流场中心区域的流向平均速度,削弱了流场中心区域的湍流强度.此外,跟单向耦合相比,双向耦合情况下的颗粒分布更加均匀,并且均匀程度随固相载率的增加而增大.     

考虑双向耦合效应的格子Boltzmann气固两相湍流模型

在流体粒子概率密度函数输运方程中考虑颗粒对流体的反作用力,发展了考虑双向耦合效应的LB气固两相流模型,引入Smagorinsky亚格子模型模拟高雷诺数气相流场.对经典后台阶气固两相流动进行模拟,气相和颗粒相速度分布与实验结果进行比较,发现考虑双向耦合效应的LB气固两相流模型结果明显优于单向耦合结果.进一步研究不同惯性颗粒在流场中的弥散特性,小颗粒(St~O(0.1))对流体的跟随性较好,在流场中分布较为均匀;而St~O(1)的颗粒难被流场涡卷吸进入涡内,呈现倾向性弥散现象;大颗粒(St~O(10))由于自身惯性进入流场涡,在流场中分布较为均匀.     

三维气固两相混合层湍流拟序结构的直接数值模拟

本文对三维气固两相混合层湍流拟序结构进行了直接数值模拟。气相流场应用拟谱方法对Ｎ－Ｓ方程组进行直接求解,通过模拟时间模式的混合层流动,分析流场失稳后涡的卷起、配对、合并及撕裂过程,研究三维混合层湍流拟序结构的演变特征;计算颗粒场时,采用Ｌａｇｒａｎｇｉａｎ方法,针对颗粒不同的Ｓｔｏｋｅｓ数,模拟了颗粒场的瞬态分布,并分析流场三维大涡结构对颗粒分布的影响。     

基于SBES方法的压气机叶片动态载荷研究

针对RANS湍流模型对流场中分离涡的模拟精度不足,导致无法准确计算出分离涡引起的叶片动态载荷问题,本文采用混合RANS-LES湍流模型中的SBES方法,对带进口导叶的NASA Stage67压气机进行了流固耦合数值模拟,研究上游叶片尾迹脱落涡、叶尖泄漏涡引起的压气机叶片动态载荷及气弹响应特征,并将SBES方法和SST湍流模型的计算结果进行了对比分析。结果表明,SBES方法相比SST湍流模型能更准确地模拟出叶片尾迹脱落涡、叶尖泄漏涡的详细结构,计算的叶片表面动态载荷波动特征更加复杂,频谱中频率成分更加丰富。两种湍流模型计算的导叶尾迹扫掠引起的转子叶片高频气弹响应幅值相近,而在计算流场中大尺度分离涡引起的叶片低频气弹响应幅值上存在很大差异。     

空间发展湍流气粒两相平面混合层的大涡模拟与统计

本文对空间发展的湍流气固两相平面混合层流动进行了大涡模拟研究,其中气相亚网格尺度（SGS）使用结构函数模型,气相控制方程组采用SIMPLE方法求解,固体颗粒运动用拉格朗日方法计算。计算结果正确重现了流体涡结构的卷起、合并和破碎过程,以及小尺寸颗粒在涡边缘（低涡度区）的局部富集现象。对直径分别为42μm、72μm和135μm分别进行了模拟,并将统计结果和实验测量结果（Hishida　et　al[1]）比较,表明两者的平均速度吻合很好,但颗粒数密度和脉动速度存在较明显的差异,因此有必要对亚网格应力和颗粒之间的耦合作用以及拟序结构的三维性对颗粒运动的影响开展深入研究。     

三维格子涡方法模拟自由下落颗粒群

本文发展了基于双向耦合的三维格子涡方法,采用涡方法模拟流场中的涡量变化及涡元运动,使用双势法求解速度场,采用拉格朗日方法跟踪颗粒相.利用该模型模拟了颗粒群的自由下落及下落过程引起的气相流动,模拟结果与实验吻合良好。结果显示空气相速度径向分布满足高斯分布,颗粒相流量对颗粒群的扩散影响不显著,而流量增大会引起颗粒速度的增大.     

基于大涡模拟下风力机翼型非定常转捩流动的动态模态分解分析

翼型绕翼流动对风力机整机性能产生重要影响.本文基于大涡模拟方法,得到风力机翼型非定常转捩流动.通过对压力流场的动态模态分解(DMD)分析,发现翼型层流分离泡的生成和发展为流动主要非定常特征,且该特征具有主要频率.预估得到高增长率T-S扰动波频率与DMD模态频率接近,发现高频变化的分离泡由边界层分离点下游不远处的T-S扰动波诱导K-H不稳定性而主导,以及中低频的DMD模态表现为T-S扰动波所引起的湍流结构.对比不同攻角下的DMD分解结果,发现上表面分离泡会逐渐向前移动,长度变短;而下表面分离泡略微后移,长度变化不大.     

两相圆湍射流大颗粒调制作用的PDPA实验研究

采用相多普勒颗粒分析仪(PDPA)测量了雷诺数为8500的两相圆湍射流速度场,研究了两种大粒径颗粒对湍流的调制作用.发现在相同气载比条件下,大粒径颗粒加入后,在射流近场区域(X/D<10)气相湍流度增强,沿着射流中心轴线轴向和径向气相湍流度增强;较小粒径颗粒加入后削弱气相湍流度.在射流远场X/D=20截面,较大粒径颗粒的加入会削弱气相湍流度,但其削弱程度不如小粒径颗粒.随着粒径的增大,两相射流湍流的气相湍流度进一步增加.     

基于拟颗粒的气固两相曳力模型研究

为了研究气固两相流动大涡模拟中合适的曳力计算模型,本文引入拟颗粒和拟颗粒表面能的概念,通过拟颗粒表面能与外界输入能量之间的平衡关系来确定拟颗粒的粒径。根据拟颗粒粒径,得到运算量较小且考虑颗粒团聚效应的曳力计算模型。应用本文的曳力计算模型对二维竖直槽道内稠密气固两相流动进行了大涡模拟,结果表明颗粒的浓度分布具有上稀下浓,壁面附近浓中心稀及颗粒聚集等特点。这与实验结果在定性上是一致的。对气相和颗粒相的瞬时速度场进行了分析,发现气相和颗粒相速度场分布的非对称性是形成颗粒浓度分布壁面附近浓中心稀的重要原因之一。     

各向同性湍流中颗粒引起湍流变动的直接模拟

本文通过直接数值模拟对均匀各向同性湍流中颗粒对湍流的变动作用进行了研究.颗粒相的体积分数很小而质量载荷足够大,以至于颗粒之间的相互作用可以忽略不计,而重点考虑颗粒与湍流间能量的交换。颗粒对湍流的反向作用使得湍动能的耗散率增强,以至于湍动能的衰减速率增大.湍动能的衰减速率随颗粒惯性的增大而增大。三维湍动能谱显示,颗粒对湍动能的影响在不同的尺度上是不均匀的。在低波数段,流体带动颗粒,而高波数段则相反.     

Diffusion-broadened line shape near a turning point

The absorption spectrum of a gas with a large Doppler width and soft collisions between particles is studied. The particles are assumed to have a nonlinear dependence of the resonance frequency on velocity. The shape of the narrow peak in the spectrum resulting from an extremum of this dependence is calculated analytically for the first time. In the absence of collisions it has a characteristic asymmetric shape. Collisions are shown to broaden the line and change its shape. The profile of the probe-field spectrum is obtained for a three-level system with the strong field at an adjacent transition. The components of the Autler-Townes doublet are found to spread and repel each other because of collisions. Pis’ma Zh. éksp. Teor. Fiz. 68, No. 1, 27–32 (10 July 1998) Published in English in the original Russian journal. Edited by Steve Torstveit.     

Dynamo transformation of the collisional R–T in a weakly ionized plasma

Theoretical prediction of a new kind of normal mode behaviour of electro-mechanical nature was first time reported by Dwivedi and Das in 1992 in the context of mesospheric modeling of observed neutral induced turbulence. Local dynamo action (due to relative neutral flow) governs the basic physical principle for linear excitation of the neutral induced low frequency instability (NILF) in mesospheric plasma, which comprises of weakly ionized inhomogeneous gas confined by the external gravity and ambient magnetic field. The present contribution offers physical explanation in terms of dynamo transformation of neutral drag effect as a source to understand complete suppression of the usual collisional R-T and in turn linear driving of the NILF. It is therefore emphasized, worth calling it as the dynamo instability.     

Flow visualization and laser measurement on particle modulation to gas-phase turbulence

Particle modulation to turbulence is investigated experimentally by means of PDPA, PIV and flow visualization for a gas-particle two-phase jet flow. Large particles can enhance the small-scale vortex, so that gas-phase turbulence intensity is increased, while small particles may delay the rolling up of the gas vortex, so that gas-phase turbulence intensity is attenuated. The critical particle size range for such different effects is between 150 σm and 200 σm, corresponding to the Stokes number is between 88 and 157 under the present flow conditions. The PIV results show small particles can retain the gas-phase vortex structure, while large particles can break large vortex structure. The particle Stokes number is not the only judgment standard whether particles enhance or attenuate gas-phase turbulence. The CTI (Change of Turbulence Intensity) number can mark off particle modulation on turbulence in two-phase flow, but more studies are needed.     

3-D numerical analysis of EHD turbulent flow and mono-disperse charged particle transport and collection in a wire-plate ESP

The present study attempts to develop a detailed numerical approach and a simulation procedure to predict the motion of gas, ions and particles inside a simple parallel plate channel containing a single corona wire. A hybrid Finite Element (FEM)-Flux Corrected Transport (FCT)-Finite Volume (FVM) method is used: the FEM–FCT numerical algorithm is applied for modeling the steady-state corona discharge, while the turbulent gas flow and the particle motion under electrostatic forces are modeled using the commercial CFD code FLUENT. Calculations for the gas flow are carried out by solving the Reynolds-averaged Navier–Stokes equations and turbulence is modeled using the k–? turbulence model. An additional source term is added to the gas flow equation to include the effect of the electric field, obtained by solving a coupled system of the electric field and charge transport equations using User-Defined Functions (UDFs). The particle phase is simulated based on the Lagrangian approach, where a large number of particles is traced with their motion affected by the gas flow and electrostatic forces using the Discrete Phase Model (DPM) in FLUENT. The developed model is useful to gain insight into the particle collection phenomena that take place inside an ESP.     

Calculations of the flow resistance and heat emission of a sphere in the laminar and high-turbulent gas flows

An early drag crisis can occur at high turbulence of incoming gas flow to a sphere. To study the influence of a crisis on heat transfer from a sphere to gas, a numerical experiment was carried out in which the free gas flow around a sphere with a temperature lower than the sphere temperature was simulated for two cases. The flow was laminar in the first case and highly turbulent in the second case. To take into account turbulence, the kinematic coefficient of turbulent viscosity with a value, which is much higher (up to 2000 times) than that for physical viscosity, was introduced. The results of calculations show that the early drag crisis occurs at Reynolds numbers of about 100 and results in considerable (by four to seven times) decrease in the hydrodynamic force and sphere drag coefficient C _d . The early drag crisis is also accompanied by the crisis of heat transfer from a sphere to gas with a decrease in Nusselt numbers Nu by three to six times.     

Lagrangian chaos and the effect of drag on the enstrophy cascade in two-dimensional turbulence

We investigate the effect of drag force on the enstrophy cascade of two-dimensional Navier-Stokes turbulence. We find a power law decrease of the energy wave number (k) spectrum that is faster than the classical (no-drag) prediction of k(-3). It is shown that the enstrophy cascade with drag can be analyzed by making use of a previous theory for finite lifetime passive scalars advected by a Lagrangian chaotic fluid flow. Using this we relate the power law exponent of the energy wave number spectrum to the distribution of finite time Lyapunov exponents and the drag coefficient.     

Fluctuations and Correlations of Pure Quantum Turbulence in Superfluid 3He-B

We describe the first measurements of line-density fluctuations and spatial correlations of quantum turbulence in superfluid 3He-B. All of the measurements are performed in the low-temperature regime, where the normal-fluid density is negligible. The quantum turbulence is generated by a vibrating grid. The vortex-line density is found to have large length-scale correlations, indicating large-scale collective motion of vortices. Furthermore, we find that the power spectrum of fluctuations versus frequency obeys a -5/3 power law which verifies recent speculations that this behavior is a generic feature of fully developed quantum turbulence, reminiscent of the Kolmogorov spectrum for velocity fluctuations in classical turbulence.