气粒两相流传热问题的光滑离散颗粒流体动力学方法数值模拟

在气粒两相流动问题中,颗粒间以及气体与颗粒间的传热问题不可忽略.光滑离散颗粒流体动力学(SDPH)模型作为一种新的求解气粒两相流动问题的方法,已经成功应用于模拟风沙运动等问题.在此基础上,提出了SDPH方法的热传导模型,模拟了气粒两相流动问题中的热传导过程以及颗粒蒸发过程.首先引入各相的能量方程,利用有限差分与光滑粒子流体动力学一阶导数相结合的方法,处理各相内部热传导项中的二阶导数问题,基于气粒两相间温度差及对流换热系数计算颗粒与气体间的热传导量,推导得到了含热传导模型的气粒两相流SDPH计算方程组,模拟计算了圆盘形颗粒团算例及鼓泡流化床内部热传导算例,并与双流体模型计算结果进行对比,结果基本符合;其次利用离散液滴模型中的颗粒蒸发传质传热定律计算颗粒的蒸发过程,数值模拟了颗粒射流蒸发过程,并与离散颗粒模型结果进行对比,两者符合得较好,验证了该方法的准确性及实用性.     

直接模拟中不同边界条件的实施及对沉降规律的影响

在牛顿流体中, 对颗粒在4种不同边界的垂直通道中的沉降运动进行了直接数值模拟. 计算结果表明:通过计算区域随颗粒运动而移动构建的无限长通道能准确模拟颗粒自由下落到稳定沉降的发展过程; 周期性边界条件由于流场变化, 对颗粒沉降产生了影响, 不能模拟颗粒的自由沉降过程; 底部封闭边界适合模拟封闭容器内颗粒与固壁的相互作用过程, 若颗粒达到稳定沉降, 也能模拟无限长通道内的沉降过程; 流化边界适合模拟流化床内气固两相流动. 计算结果有助于更好地理解和使用不同边界条件. 关键词： 直接数值模拟 边界条件 沉降 任意拉格朗日-欧拉方法     

突扩圆管内液-固两相流固体颗粒运动特性的DPM数值模拟

采用考虑颗粒碰撞的欧拉一拉格朗日数值模拟方法(DPM),对水平突扩圆管中液固两相流固体颗粒的碰撞过程进行了数值计算.在模型中,对液相采用欧拉法建立控制方程,对离散颗粒采用拉格朗日方法模拟.采用硬球模型描述颗粒间的碰撞作用.计算结果表明,该模型可以真实地模拟液固两相流中固体颗粒运动的动态变化过程以及颗粒的非均匀分布特征,从单颗粒层次上提供颗粒的运动信息,这有助于深入研究液固两相流中固体颗粒的运动规律.     

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

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

基于间断有限元方法的焦炭燃烧直接模拟

焦炭燃烧是煤燃烧的重要过程,准确模拟焦炭燃烧,深入研究焦炭燃烧的物理规律,对于预测及优化煤燃烧过程有很大意义。间断有限元方法具有良好几何适应性和稳定性,是准确模拟气固两相流动问题的较好方法。然而该方法尚未应用于求解气固两相燃烧,且存在通量重构格式无法应用于动网格计算的问题。本文利用间断有限元方法进行空间离散,结合任意拉格朗日-欧拉方法处理动网格问题,进行了二维流场下焦炭颗粒燃烧过程的直接模拟,并将该方法应用于多个焦炭颗粒运动燃烧的计算,得到了较准确的结果,揭示了单颗粒及运动多颗粒焦炭燃烧过程中流场的变化规律。     

鼓泡床内气固两相流的变尺度格子气模拟

基于格子气自动机的建模思想,建立气固两相流的变尺度格子气模型,用微观层次的气固作用规则描述两相流的宏观行为.二维流场空间被离散成两层不同尺度的正六边形网格,固体颗粒和气体粒子分别在各自的网格上运动;在格子气自动机基本规则基础上设计气固相作用规则和附加演化规则;确定模型宏观物理量的统计计算方法,并根据相似原理建立模型物理量与真实物理量之间的转换关系.用建立的模型对鼓泡床内气固流动行为进行模拟,通过比较发现,变尺度格子气模型的模拟结果与文献中的实验结果和双流体模型的模拟结果吻合较好,且平均相对误差更小,表明建模方法的正确性和模型的有效性.     

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

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

喷动流化床的数值模拟研究

本文借鉴 Ｃｕｎｄａｌｌ＆ Ｓｔｒａｃｋ提出的碰撞模型,将欧拉方法和拉格朗日法结合起来,建立了喷动流化床气固两相流动的数学模型,在计算机上对试验台进行了模拟,探讨了床内颗粒的运动特征。按１：１的比例建立了喷动流化床的模型试验台并进行了试验,对模拟结果和试验结果进行了比较与分析。颗粒在喷动床内的运动特征,模拟结果都能够再现出来,模拟结果与试验结果取得了很好的一致性。     

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

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

六分离器循环流化床内气固流动特性的LES-DEM研究

工业生产中通常采用多分离器布置的方案增大循环流化床的容量.本文采用大涡模拟-离散单元法(Large eddy simulation-Discrete element method,LES-DEM)对六分离器循环流化床内气固流动特性进行了全循环数值模拟研究。该方法在欧拉框架下求解流体的运动,在拉格朗日框架下追踪每个颗粒的运动。结果表明:六个旋风分离器回路中存在着气固分配不均现象,旋风分离器的轴对称布置方式优于中心对称布置方式.提升管内呈现强烈的颗粒聚团行为,颗粒团的体积大致为正态分布频率,其纵横比都大于1,并且随着提升管高度先增加后减小.     

Time-resolved tracking of a sound scatterer in a complex flow: nonstationary signal analysis and applications

It is known that ultrasound techniques yield nonintrusive measurements of hydrodynamic flows. For example, the study of the echoes produced by a large number of particles insonified by pulsed wavetrains has led to a now-standard velocimetry device. In this paper, a new technique for the measurement of the velocity of individual solid particles moving in fluid flows is proposed. It relies on the ability to resolve in time the Doppler shift of the sound scattered by the continuously insonified particle. For this signal-processing problem two classes of approaches can be used: time-frequency analysis and parametric high-resolution methods. In the first class the spectrogram and reassigned spectrogram is considered, and applied to detect the motion of a small bead settling in a fluid at rest. In nonstationary flows, methods in the second class are more robust. An approximated maximum likelihood (AML) technique has been adapted, coupled with a generalized Kalman filter. This method allows for the estimation of rapidly varying frequencies; the parametric nature of the algorithm also provides an estimate of the variance of the identified frequency parameters.     

Parallel computation and numerical visualization for non-uniform particles included in gas and liquid flows

A parallel computation method has been proposed for the mixing and segregation of granular mixture included in gas and liquid flows. In this method, a three-dimensional (3D) computational volume is decomposed into multiple sub-blocks and their geometries are represented by 3D body-fitted coordinates. The fluid-particle interactions are treated by two types of models: a two-way model for liquid-solid flows and a one-way model in case of gas-solid flows. The computations of the particle motions in the multiple sub-blocks are executed simultaneously on the basis of the distinct element method (DEM). Since a graphic process is also executed as one of the parallel jobs, the particle distributions can be visualized during the computations. The computational method was applied to the gas-solid flows consisting of different diameters and densities in the horizontal and inclined cylinders rotating around their axes. From the comparison with the experimental results, nearly uniform mixing and particle segregation are successfully predicted in the oscillating liquid flows. In addition, it has been indicated that the particle pathline is very effective to visualize and understand the flow patterns of the particles with different properties. The result of the computations for the liquid-solid flows demonstrated that the vertical segregation of the non-uniform particles is reasonably reproduced.     

Flows with suspended and floating particles

The evolution of the configuration of a set of particles dispersed in a flowing liquid is crucial in many applications such as sedimentation, slurry transport, rheology and structured arrays of micro- and nano-particles. Direct simulation based on what is called fictitious domain method coupled with finite element method has been used to study particulate flows and sedimentation process. Here we extend the previously proposed formulations to naturally include buoyancy force and the capillary driven attraction or repulsion of particles located at fluid interfaces. The set of differential equations is discretized using a fully implicit-fully coupled fictitious domain/finite element approach, avoiding numerical instabilities that may arise from explicit integration. The proposed formulation and implementation are validated by comparing the predictions of simple 2D flows to available numerical or analytical solutions. The method is then used to analyze the flotation of 2D particles and capillary driven aggregation at fluid interfaces.     

A lattice Boltzmann method for a two-phase flow containing solid bodies with viscoelastic membranes

A lattice Boltzmann method (LBM) for two-phase flows containing solid bodies with viscoelastic membranes is proposed. The method is based on the two-phase LBM, in which one phase is regarded as the solid phase. In the present model, the membrane is assumed to be composed of identical particles that are connected to their neighboring particles by elastic springs to take account of stretching and compression effects. The method is applied to two representative problems, namely the behavior of a viscoelastic body under shear flow and the motion of a viscoelastic body in a Poiseuille flow. Tank-tread motion and axial migration, which are both characteristic of the motion of viscoelastic bodies, are simulated by using the method. These results indicate that the method is capable of simulating the complex behavior of viscoelastic bodies in capillaries, such as the motion of red blood cells in blood flows.     

Entanglement entropy and mutual information production rates in acoustic black holes

A method to investigate acoustic Hawking radiation is proposed, where entanglement entropy and mutual information are measured from the fluctuations of the number of particles. The rate of entropy radiated per one-dimensional (1D) channel is given by S=κ/12, where κ is the sound acceleration on the sonic horizon. This entropy production is accompanied by a corresponding formation of mutual information to ensure the overall conservation of information. The predictions are confirmed using an ab initio analytical approach in transonic flows of 1D degenerate ideal Fermi fluids.     

改进粒子水平集法及其应用

针对传统粒子水平集方法的缺陷提出一种虚拟粒子对距离函数的修正方法.改进算法针对虚拟粒子与待修正网格点的各种相对位置进行误差修正,可以避免虚拟粒子在非界面法线方向移动而产生的误差.数值求解涡旋流场问题显示:改进算法的界面捕捉性能得到显著提高.在此基础上,结合投影法求解Navier-Stokes方程组,数值模拟存在高密度比(10³量级)及粘性比(10²量级)的不可压缩两相流问题(气泡上升及自由面不稳定性问题).计算结果与已有结果吻合良好,详细分析主要控制参数对瑞利-泰勒不稳定性现象的影响.     

Numerical investigation of two-phase flows in plug nozzles

The results of numerical investigations of gas dynamics of two-phase (including vortex) flows in plug nozzles are presented. It is shown that most particles settle on the nozzle walls, which are contoured for ideal-gas flow. A method is studied for decreasing with the help of vorticity the precipitation of particles on the plug contour. Methods are proposed for contouring plug nozzles for two-phase flows.Scientific-Research Institute of Applied Mathematics and Mechanics at Tomsk State University. Academician V. P. Makeev Design Office of Machine Building. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 4, pp. 81–91, April, 1993.     

An energy-consistency-preserving large eddy simulation-scalar filtered mass density function (LES-SFMDF) method for high-speed flows

An energy-consistency-preserving large eddy simulation-scalar filtered mass density function (LES-SFMDF) method is developed to improve the existing LES-SFMDF method in high-speed flows, especially supersonic flows. The high-speed source term in the SFMDF transport equation is analysed theoretically from a new point of view, and then several primary principles are proposed for the LES-SFMDF to achieve a good consistency even along discontinuities. Based on these principles and further theoretical analysis, the high-speed source term of the enthalpy in the SFMDF is modelled and computed from both MC and LES variables rather than the usually used solely LES variables. This new LES-SFMDF method is used for simulating the flows in a shock tube and in a subsonic temporally developing mixing layer. This method shows a better particle energy consistency than the existing method when applied across discontinuities in supersonic laminar flows. Unlike the existing method, with this energy-consistency-preserving LES-SFMDF method, particle energy consistency is consistent with particle mass consistency so that particle energy consistency can benefit from particle velocity correction. This method also demonstrates robustness for various numbers of particles.     

A Multi-Speed Compressible Lattice-Boltzmann Model

An extension of the lattice Boltzmann BGK method to compressible flows is presented that combines three novel additions: (1) particles move density and energy weights in multiple velocity bins (11 for 1-D flow) to nearby cell centers. (2) the equilibrium distribution remains an unexpanded Maxwellian; and (3) transport and relaxation to equilibrium are performed implicitly at each node. These advances allow for the parallel modeling of high Mach number shocks and high Reynolds number flows, while avoiding advective numerical diffusion, the need for Riemann solvers, and non-linear limiters. A 1D shock tube application is shown. Generalization to higher dimensions and multi-materials are discussed.     

自由剪切流动中颗粒扩散的自相似特性

为了研究自由剪切流动中颗粒扩散的统计特性,对空间发展模式的三维气固两相射流和时间发展模式的三维气固两相混合层进行了直接数值模拟.其中对气相不可压缩Navier-Stokes方程的求解分别采用有限容积方法和拟谱方法,对离散惯性颗粒的跟踪采用单向耦合的拉格朗日方法.统计结果显示,当流场进入自相似状态后,以流场速度梯度二阶不变量Q为自变量的颗粒数的概率分布函数也呈现出与颗粒尺寸和时间发展无关的自相似特性.