基于离散元法的干湿颗粒系统仿真

介绍了基于离散元法的干湿颗粒系统仿真软件DEMSIM。对于干颗粒系统,DEMSIM可以分析二维和三维颗粒系统的弹性和塑性接触碰撞过程;对于湿颗粒系统,DEMSIM采用传统的液桥模型;对于颗粒-流体系统,DEMSIM采用CFD-DEM细观耦合模型模拟。一系列典型算例的模拟分析,验证了干湿颗粒系统仿真软件DEMSIM的精度和有效性。     

二维喷动床CFD-DEM模拟中曳力模型的影响

CFD-DEM已经广泛应用到喷动床的研究中,其模拟的准确性与用于处理颗粒-流体相互作用的曳力模型密切相关。为了探究不同曳力模型对喷动床CFD-DEM模拟结果的影响,基于非结构化网格的喷动床仿真,使用7个曳力模型分别对锥底喷动床内气固两相运动进行了数值模拟。综合床层压降、喷动高度和颗粒速度特性三个方面,Wen-Yu模型和Gibilaro模型预测的气固两相运动最剧烈,其次是Di Felice模型、Syamlal-O’Brien模型、Gidaspow模型和Huilin-Gidaspow模型,BVK模型预测的气固两相运动最平缓。由于模拟的气固两相体系属于密相体系,Huilin-Gidaspow模型的光滑过渡函数没有产生效果,所以Gidaspow模型和Huilin-Gidaspow模型在各个方面的预测结果基本一致。     

Entropy productions in granular materials

Granular materials display more abundant dissipation phenomena than ordinary materials. In this paper, a brief energy flow path with irreversible processes is illustrated, where the concept of granular temperature T_g, initially proposed for dilute systems, is extended to dense systems in order to quantify disordered force chain configurations. Additionally, we develop the concept of conjugate granular entropy s_g and its production equation. Our analyses find out that the granular entropy significantly undermined the elastic contact between particles, seriously affecting the transport coefficients in granular materials and creating new transport processes.     

Particle shape effect on hydrodynamics and heat transfer in spouted bed: A CFD–DEM study

Spouted bed has drawn much attention due to its good heat and mass transfer efficiency in many chemical units. Investigating the flow patterns and heat and mass transfer inside a spouted bed can help optimize the spouting process. Therefore, in this study, the effects of particle shape on the hydrodynamics and heat transfer in a spouted bed are investigated. This is done by using a validated computational fluid dynamics–discrete element method (CFD–DEM) model, considering volume–equivalent spheres and oblate and prolate spheroids. The results are analysed in detail in terms of the flow pattern, microstructure, and heat transfer characteristics. The numerical results show that the prolate spheroids (Ar = 2.4) form the largest bubble from the beginning of the spouting process and rise the highest because the fluid drag forces can overcome the interlocking and particle–particle frictional forces. Compared with spherical particles, ellipsoidal spheroids have better mobility because of the stronger rotational kinetic energy resulting from the rough surfaces and nonuniform torques. In addition, the oblate spheroid system exhibits better heat transfer performance benefiting from the larger surface area, while prolate spheroids have poor heat transfer efficiency because of their orientation distribution. These findings can serve as a reference for optimizing the design and operation of complex spouted beds.     

Correlation of data concerning resistance to flow of generalized Newtonian fluids through granular beds

Summary The paper is concerned with the pressure drop during the flow of rheologically complex fluids through granular beds. An approach is proposed which enables the correlation of the data for any generalized Newtonian fluid. As an example, a detailed derivation of a correlation equation is presented for Carreau fluids. The applicability of the derived equation was proved experimentally in the case of flow of molten poly(ethylene terephthalate) through granular beds.

---

Zusammenfassung Die Arbeit betrifft den Druckverlust bei der Strömung von Flüssigkeiten mit komplexen rheologischen Eigenschaften durch Kornschüttungen. Es wird ein Verfahren vorgeschlagen, das eine Korrelation der Daten für beliebige verallgemeinerte newtonsche Flüssigkeiten ermöglicht. Als Beispiel wird die Ableitung der Korrelationsgleichung für Carreau-Flüssigkeiten ausführlich dargestellt. Die Anwendbarkeit der abgeleiteten Gleichung wird für die Strömung von geschmolzenem Polyäthylenterephthalat durch Kornschüttungen experimentell bestätigt.

---

a _T temperature shift factor - A constant in eq. [1] - b _M molecular weight shift factor - d _p effective particle diameter, m - D capillary diameter, m - f _BK friction factor, defined by eq. [4] - K constant in eq. [22] - K ₀ constant depending on the shape of a conduit cross-section - K ₁ constant in eq. [11] - l bed height, m - l _e average length of the bed channels, m - L length of a conduit, m - M _n number-average molecular weight - M _w weight-average molecular weight - N constant in Carreau model - p pressure drop due to friction, Pa - r _h hydraulic radius, m - T absolute temperature, K - v mean linear velocity, m/s - v _e mean linear velocity in the bed channels, m/s - v ₀ superficial velocity, related to an empty cross-section of the column, m/s - dimensionless factor, defined by eq. [21] - shear rate, s^–1 - nominal shear rate at the wall of a circular pipe, s^–1 - average shear rate at the wall of a noncircular channel, s^–1 - bed porosity - shear dependent viscosity, Pa s - ₀ zero-shear rate viscosity, Pa s - infinite-shear rate viscosity, Pa s - [] intrinsic viscosity - time constant in Carreau model, s - µ Newtonian viscosity, Pa s - fluid density, kg/m³ - _w shear stress at the wall of a circular pipe, Pa - average shear stress at the wall of a noncircular channel, Pa - {Re} general form of Reynolds number in eq. [1] - Re _BK modified Reynolds number for Newtonian fluids, defined by eq. [25] - Re _BK ^* generalized Reynolds number for Carreau fluids, defined by eq. [24] With 5 figures and 2 tables     

Study of fluid cell coarsening for CFD-DEM simulations of polydisperse gas–solid flows

Particle polydispersity is ubiquitous in industrial fluidized beds, which possesses a significant impact on hydrodynamics of gas–solid flow. Computational fluid dynamics-discrete element method (CFD-DEM) is promising to adequately simulate gas–solid flows with continuous particle size distribution (PSD) while it still suffers from high computational cost. Corresponding coarsening models are thereby desired. This work extends the coarse-grid model to polydisperse systems. Well-resolved simulations with different PSDs are processed through a filtering procedure to modify the gas–particle drag force in coarse-grid simulations. We reveal that the drag correction of individual particle exhibits a dependence on filtered solid volume fraction and filtered slip velocity for both monodisperse and polydisperse systems. Subsequently, the effect of particle size and surrounding PSD is quantified by the ratio of particle size to Sauter mean diameter. Drag correction models for systems with monodisperse and continuous PSD are developed. A priori analysis demonstrates that the developed models exhibit reliable prediction accuracy.     

On boundary potential energies in deformational and configurational mechanics

This contribution deals with the implications of boundary potential energies, i.e. in short surface, curve and point potentials, on deformational and configurational mechanics. Within the realm of deformational mechanics the surface/curve potentials are allowed in the most general case to depend on the deformation, the surface/curve deformation gradient and the spatial surface normal/curve tangent and are parametrised in the material placement and the material surface normal/curve tangent. The point potentials depend on the deformation and are parametrised in the material placement. From the configurational mechanics perspective the roles of fields and parametrisations are reversed. By considering variational arguments based on the kinematics of deforming surfaces/curves, in particular the relevant surface/curve stresses and distributed forces contributing to (localized) deformational and configurational force balances at surfaces/curves/points, which extend the common traction boundary conditions, are derived. Thereby, dissipative distributed configurational forces that are energetically conjugate to configurational changes are introduced as definitions. The (localized) force balances at surfaces/curves/points together with the contributing stresses and distributed forces within deformational and configurational mechanics display an intriguing duality. The resulting dissipative configurational tractions at the boundary are exemplified for some illustrative cases of boundary potentials.     

Numerical modeling of fluid-particle interaction in granular media

Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular media as a whole. This paper presents a coupled computational fluid dynamics and discrete element method (CFD-DEM) approach for this purpose. The granular particle system is modeled by DEM, while the fluid flow is simulated by solving the locally averaged Navier–Stokes equation with CFD. The coupling is considered by exchanging such interaction forces as drag force and buoyancy force between the DEM and CFD. The approach is benchmarked by two classic geomechanics problems for which analytical solutions are available, and is further applied to the prediction of sand heap formation in water through hopper flow. It is demonstrated that the key characteristic of granular materials interacting with pore water can be successfully captured by the proposed method.     

颗粒介质结构与力学特征研究综述

颗粒介质由大量离散的粗颗粒聚集而成,如自然界中的粗砂和碎屑堆积体等. 在工程实践中,人们依据经验和实验数据建立了许多模型,虽然可以满意地描述某些力学现象,但是对颗粒介质力学性质全貌的认识以及颗粒介质物理本质的理解仍远远不够. 颗粒介质长程无序、短程有序的结构和复杂的能量转化过程,注定了其独特的力学性质. 该文综述了颗粒介质结构探测和表征技术、热力学理论和固态-流态转变方面的新进展,特别介绍了清华大学近5 年来开展的颗粒介质结构模型化方法和双颗粒温度热力学理论. 最后,提出了开展结构分析-热力学理论的联合研究思路,以期更加深入认识颗粒介质的力学特性,探究颗粒介质的热力学根源,改善现有唯象研究现状.     

Aeroelastic stability consideration of supersonic flight vehicle using nonlinear aerodynamic response surfaces

Aeroelastic stability of a flexible supersonic flight vehicle is considered using nonlinear dynamics, nonlinear aerodynamics, and a linear structural model. Response surfaces including global multivariate orthogonal modeling functions are invoked to derive applied nonlinear aerodynamic coefficients. A modified Gram–Schmidt method is utilized to orthogonalize the produced polynomial multivariate functions, selected and ranked by predicted squared error metric. Local variation of angle-of-attack and side-slip angle is applied to the analytical model. Identification of nonlinear aerodynamic coefficients of the flight vehicle is conducted employing a CFD code and the required analytical model for simulation purposes is constructed. The method is used to determine the aeroelastic instability and response of a selected flight vehicle.