雾层两相耦合流场及气溶胶颗粒物的动力学性质

雾层气溶胶系统涉及复杂的动力学演变过程:碰撞、凝并、破碎、冷凝/蒸发、成核、沉积、表面化学反应等.因此,发展雾层与气相流场耦合的Eulerian-Lagrangian两相流模型、颗粒动力学及随机轨道模型,考虑重力、曳力、布朗力、Basset力等对颗粒相的作用.基于SIMPLE和多重Monte Carlo算法求解颗粒群平衡方程,自行开发了FAD程序首先对室内燃烧源细微颗粒物的扩散实验展开数值模拟,计算结果与实验数据吻合较好.将建立的模型和方法数值研究气溶胶污染物在雾环境中的输运过程,分析雾消散阶段颗粒相浓度、平均尺度的时空分布.结果显示:当时间演化至60 min,雾滴的平均尺度减小到初始的65.67％,而气溶胶颗粒最大数目对应的尺度为0.006 μm.     

气固两相流中颗粒碰撞的Monte-Carlo数值模拟

利用颗粒碰撞动力学模型和颗粒几何碰撞率模型,采用Monte—Carlo算法来模拟颗粒之间碰撞,把该算法与求解雷诺应力-概率密度函数模型的有限差分-Monte Carlo算法耦合起来,对轴对称突扩通道内的两相旋流场进行了数值模拟,模拟结果表明,由于颗粒碰撞使颗粒的动能和湍动能在三个坐标方向上进行了再分配,从而导致颗粒的动能和湍动能在三个坐标方向上趋于各向同性;另外,由于颗粒碰撞破坏了颗粒-颗粒、颗粒-流体微团之间的速度关联,从而造成颗粒湍动能及两相速度脉动关联的降低。     

硬球加强模型在CFD-DEM耦合计算中的验证与分析

针对CFD-DEM耦合计算中,颗粒计算时间步的选取影响颗粒碰撞计算精度和效率的问题。本文引入插值算法,将动量定理求解颗粒碰撞前后速度进行加权平均;根据弹性理论计算得到颗粒碰撞力,进行动力学方程求解;通过速度收敛准则修正初值速度并自动调整迭代求解次数,提出一种计算精度不受计算时间步长影响,无需对碰撞过程进行精细描述的高效率和高精度的加强硬球模型。对两个颗粒匀和变速碰撞算例进行数值模拟,碰撞后速度、碰撞力和碰撞时间与理论计算误差小于4%,与采用软球碰撞模型的DEM方法相比,颗粒碰撞计算精度不受计算时间步长影响,计算效率提高36.3%和36.8%。对单个颗粒在静水中沉降进行数值模拟,计算步长取10 s～5 s,颗粒与壁面即可得到精确解,计算效率提高33.5%。通过压力损失实验验证了该模型能够准确计算颗粒体积分数小于12%条件下两相流的压力损失。     

颗粒群碰撞搜索及CFD-DEM耦合分域求解的推进算法研究

在采用计算流体力学-离散元耦合方法(computational fluiddynamics-discrete element method, CFD-DEM)进行固液两相耦合分析时, 颗粒计算时间步的选取直接影响到耦合计算精度和计算效率. 为此, 本文选取每个目标颗粒为研究对象, 引入插值函数计算时间步的运动位移, 构建可变空间搜索网格; 通过筛选可能碰撞颗粒建立搜索列表, 采用逆向搜索方式判断碰撞颗粒, 从而提出一种改进的DEM方法(modified discreteelement method, MDEM). 该算法在颗粒群与流体耦合计算中, 颗粒计算初始时间步选取不受颗粒碰撞时间限制, 通过自动调整和修正实现大步长, 由颗粒和流体耦合条件实时更新流体计算时间步, 使颗粒计算时间步选取过小导致计算效率低、选取过大导致颗粒碰撞漏判的问题得以解决, 为颗粒与流体耦合的数值模拟提供了行之有效的计算方法. 通过两个颗粒和多个颗粒的数值模拟, 得到的颗粒间碰撞力、碰撞位置及次数, 与理论计算结果的相对误差均低于2%, 与传统的DEM碰撞搜索算法相比, 在选取的3种计算时间步均不会影响计算精度, 且有较高的计算效率. 通过多个颗粒与流体的耦合数值模拟, 采用传统的CFD-DEM方法, 只有颗粒计算时间步选取10$^{-6}$ s或更小才能得到精确解, 而采用本文方法取10$^{-4}$ s也能够得到精确解, 避免了颗粒碰撞随时间步增大而出现的漏判问题, 且计算耗时降低了16.7%.      

CFRP加固RC梁的模糊随机可靠度数值模拟

CFRP加固RC结构可靠度的不确定性,不仅要考虑随机性,也要考虑模糊性.本文根据模糊数学与可靠性理论,研编了引入模糊示性函数的Monte Carlo计算程序,通过试验讨论模糊示性函数的确定方法.利用Matlab软件,直接产生随机变量数组,解决Monte Carlo 抽样困难.CFRP加固梁破坏试验获得全部随机变量的统计参量.编制基于Matlab的Monte Carlo计算程序,完成了对CFRP加固RC梁抗弯模糊随机可靠度数值模拟,模拟计算结果得到传统FOSM法验证.     

离散系统通用动力学方程求解算法的研究进展

离散系统中的颗粒物在凝并、破碎、冷凝/蒸发、成核、沉积等事件作用下颗粒尺度分布的时间演变由通用动力学方程所描述.该方程为一典型的部分积分微分方程,普通数值方法难以求解.本文详细介绍了求解通用动力学方程的矩方法、分区法、离散法、离散-分区法、MonteCarlo方法等几种算法的原理、优缺点和最新的研究进展,并着重介绍了MonteCarlo算法,包括基于时间驱动Monte Carlo方法、基于事件驱动MonteCarlo方法、常数目法、常体积法以及多重Monte Carlo算法.      

基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

Fokker-Planck方程有限解析/Monte Carlo数值模拟方法

对白噪声驱动随机系统的Fokker-Planck方程进行约化，求得约化方程的解析解，使 用局部解析解和Monte Carlo结合方法求解常系数Fokker-Planck方程，并与常系数Fokker-Planck方程的精确解 进行对比，之后求解了变驱动力系统的行为. 数值模拟结果表明，有限解析/Monte Carlo结合的方法，能成功求解一维Fokker-Planck方程，求解粒子数为10$^{5}$个，能获得 十分光滑的PDF分布曲线，计算颗粒在300个时，就能获得较好的均值. 其研究为两相 湍流PDF模型新计算方法研究提供基础.     

三方程线性弹性-阻尼DEM模型及碰撞参数确定

建立了一种考虑法向接触力、切向接触力(含静滑动摩擦力及动滑动摩擦力)和力矩(含由切向力产生的力矩及静滚动摩擦力矩和动滚动摩擦力矩)的三方程线性弹性-阻尼离散单元模型,并将该模型应用到颗粒物料的三维数值模拟中,讨论了模型中几个重要碰撞参数--刚性系数、阻尼系数及摩擦系数的选择及其对计算结果的影响,同时也探讨了时间步长等计算参数对模拟结果的影响.为了验证算法和参数选择的正确性,本文对几个有代表性的颗粒系统进行了数值试验研究,并对计算结果进行了细致的分析,验证了新模型和参数选择的正确性.     

随机结构正交展开分析的Ritz动力聚缩法

针对随机结果正交展开理论计算上的弱点 ,本文在分析扩阶矩阵特性的基础上 ,于 Ritz模态向量子空间中对扩阶方程实现动力聚缩 ,大大提高了正交展开理论对实际工程问题的分析能力。分析实例表明 :即使结构参数具有很大变异性 (如δ =0 .4 )时 ,该算法依然能理想地与 Monte Carlo法模拟结果相吻合 ,计算时间则远远小于 Monte Carlo模拟法。同时 ,分析例证再一次强化了在结构动力分析中考虑结构参数随机性的必要性     

Multi-monte-carlo method for general dynamic equation considering particle coagulation

Introduction ParticulateMatter(PM10,aerodynamicdiameterlessthan10μm)hasbeenwidely investigated.Somekindsofcombustionincludingcoalcombustion,gasoline/dieseloil combustionofvehicle,municipalsolidwastecombustion,etc,areoneofthemainsourcesof PM10[1].Particle…     

Effect of collision‐partner selection schemes on the accuracy and efficiency of the direct simulation Monte Carlo method

The effect of the collision‐partner selection scheme on the accuracy and the efficiency of the Direct Simulation Monte Carlo method is investigated. Several schemes that reduce the mean collision separation, including the fixed sub‐cell scheme, the transient adaptive sub‐cell scheme, and the virtual sub‐cell scheme, are evaluated. Additionally, a new scheme is proposed that limits the population from which collision partners are selected based on the distance traveled by a simulator and performs near‐neighbor collisions using this population. These collision‐partner selection schemes are assessed for Fourier flow (heat conduction between parallel plates) and a standard hypersonic benchmark problem (Mach 15.6 nitrogen flow over a 25–55° biconic). The new limited‐selection near‐neighbor scheme has superior performance compared to the other schemes for both flows and reduces both the spatial and temporal discretization errors relative to random‐selection and nearest‐neighbor collision‐partner selection schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

Modeling wall effects in a micro-scale shock tube using hybrid MD–DSMC algorithm

Wall effects in a micro-scale shock tube are investigated using the Direct Simulation Monte Carlo method as well as a hybrid Molecular Dynamics–Direct Simulation Monte Carlo algorithm. In the Direct Simulation Monte Carlo simulations, the Cercignani–Lampis–Lord model of gas–surface interactions is employed to incorporate the wall effects, and it is shown that the shock attenuation is significantly affected by the choice of the values of tangential momentum accommodation coefficient. A loosely coupled Molecular Dynamics–Direct Simulation Monte Carlo approach is then employed to demonstrate incomplete accommodation in micro-scale shock tube flows. This approach uses fixed values of the accommodation coefficients in the gas–surface interaction model, with their values determined from a separate dynamically similar Molecular Dynamics simulation. Finally, a completely coupled Molecular Dynamics–Direct Simulation Monte Carlo algorithm is used, wherein the bulk of the flow is modeled using Direct Simulation Monte Carlo, while the interaction of gas molecules with the shock tube walls is modeled using Molecular Dynamics. The two regions are separate and coupled both ways using buffer zones and a bootstrap coupling algorithm that accounts for the mismatch of the number of molecules in both regions. It is shown that the hybrid method captures the effect of local properties that cannot be captured using a single value of accommodation coefficient for the entire domain.     

Differentially weighted direct simulation Monte Carlo method for particle collision in gas-solid flows

In gas-solid flows,particle-particle interaction(typical,particle collision) is highly significant,despite the small particles fractional volume.Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles(e.g.,agglomeration and fragmentation) in spite of their initial monodisperse particle distribution.The conventional direct simulation Monte Carlo(DSMC)method for particle collision tracks equally weighted simulation particles,which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions.In this study,a new differentially weighted DSMC(DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian-Lagrangian models for hydrodynamics.Three schemes(mass,momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass,momentum or energy of the whole system unchanged respectively.A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency.The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best.Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC,compared with the equally weighted DSMC.Meanwhile,computational cost can be largely reduced in contrast with direct numerical simulation.     

Differentially weighted direct simulation Monte Carlo method for particle collision in gas–solid flows

In gas–solid flows, particle–particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles (e.g., agglomeration and fragmentation) in spite of their initial monodisperse particle distribution. The conventional direct simulation Monte Carlo (DSMC) method for particle collision tracks equally weighted simulation particles, which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions. In this study, a new differentially weighted DSMC (DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian–Lagrangian models for hydrodynamics. Three schemes (mass, momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass, momentum or energy of the whole system unchanged respectively. A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency. The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best. Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC, compared with the equally weighted DSMC. Meanwhile, computational cost can be largely reduced in contrast with direct numerical simulation.     

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

An analytical approach for the closure equations of gas–solid flows with inter-particle collisions

This work examines the effect of inter-particle collisions on the motion of solid particles in two-phase turbulent pipe and channel flows. Two mechanisms for the particle–particle collisions are considered, with and without friction sliding. Based on these collision mechanisms, the correlations of the various velocity components of colliding particles are obtained analytically by using an averaging procedure. This takes into account three collision coordinates, two angles and the distance between the centers of colliding particles. The various stress tensor components are obtained and then introduced in the mass, linear momentum and angular momentum equations of the dispersed phase. The current approach applies to particle–particle collisions that result from both the average velocity difference and the turbulent velocity fluctuations. In order to close the governing equations of the dispersed phase, the pseudo-viscosity coefficients are defined and determined by the time of duration of the inter-particle collision process. The model is general enough to apply to both polydisperse and monodisperse particulate systems and has been validated by comparisons with experimental data.     

Role of particle stiffness and inter-particle sliding friction in milling of particles

Discrete element method （DEM） has been used to investigate the effects of particle elastic modulus and coefficient of inter-particle sliding friction on milling of mineral particles. An autogeneous mill of 600 mm diameter and 320 mm length with 14,500 particles has been selected for the simulation. Various mill performance parameters, for example, particle trajectories, collision frequency, collision energy and mill power have been evaluated to understand the effects of particle elastic modulus and inter-particle sliding friction during milling of particles. For the given model, it has been concluded that at high energy range, as the elastic modulus and particle sliding friction increase the energy dissipated among the particles increases. The collision frequency increases with the increase in elastic modulus, however, this trend is not clearly observed with increasing inter-particle sliding friction. The power draw of the mill increases with the increase in fraction of mill critical speed.     

A novel efficient hybrid DSMC–dynamic collision limiter algorithm for multiscale transitional flows

A new 2D parallel multispecies polyatomic particle–based hybrid flow solver is developed by coupling the Direct Simulation Monte Carlo (DSMC) method with a novel Dynamic Collision Limiter (DCL) approach to solve multiscale transitional flows. The hybrid DSMC‐DCL solver can solve nonequilibrium multiscale flows with length scales ranging from continuum to rarefied. The DCL method, developed in this work, dynamically assigns different number of collisions in cells, which is based on the local value of K‐S parameter such that the number of collisions per time step is limited in near‐equilibrium flow regions. Present hybrid solver uses the Kolmogorov‐Smirnov statistical test as the continuum breakdown parameter, based on which, the solution domain is decomposed into near‐equilibrium and nonequilibrium flow regions. Direct Simulation Monte Carlo is used where nonequilibrium flow regions are encountered, while the DCL method is used where flow regions are found to be in near‐equilibrium state. In this work, we have studied hypersonic flow of nitrogen over a blunt body with an aerospike and supersonic flow of argon through a micronozzle. The results obtained by the hybrid DSMC‐DCL solver are compared and shown to agree well with the experimental data and with those obtained from DSMC, with significant savings in the computational cost.     

Coupled Particle-In-Cell and Direct Simulation Monte Carlo method for simulating reactive plasma flows

Plasma flows with high Knudsen numbers cannot be treated with classic continuum methods, as represented for example by the Navier–Stokes or the magnetohydrodynamic equations. Instead, the more fundamental Boltzmann equation has to be solved, which is done here approximately by particle based methods that also allow for thermal and chemical non-equilibrium. The Particle-In-Cell method is used to treat the collisionless Vlasov–Maxwell system, while neutral reactive flows are treated by the Direct Simulation Monte Carlo method. In this article, a combined approach is presented that allows the simulation of reactive, partially or fully ionized plasma flows. Both particle methods are briefly outlined and the coupling and parallelization strategies are described. As an example, the results of a streamer discharge simulation are presented and discussed in order to demonstrate the capabilities of the coupled method.