Accelerating CFD–DEM simulation of processes with wide particle size distributions

Size-reduction systems have been extensively used in industry for many years. Nevertheless, reliable engineering tools to be used to predict the comminution of particles are scarce. Computational fluid dynamics(CFD)–discrete element model(DEM) numerical simulation may be used to predict such a complex phenomenon and therefore establish a proper design and optimization model for comminution systems.They may also be used to predict attrition in systems where particle attrition is significant. Therefore,empirical comminution functions(which are applicable for any attrition/comminution process), such as:strength distribution, selection, equivalence, breakage, and fatigue, have been integrated into the threedimensional CFD–DEM simulation tool. The main drawback of such a design tool is the long computational time required owing to the large number of particles and the minute time-step required to maintain a steady solution while simulating the flow of particulate materials with very fine particles.The present study developed several methods to accelerate CFD–DEM simulations: reducing the number of operations carried out at the single-particle level, constructing a DEM grid detached from the CFD grid enabling a no binary search, generating a sub-grid within the DEM grid to enable a no binary search for fine particles, and increasing the computational time-step and eliminating the finest particles in the simulation while still tracking their contribution to the process.The total speedup of the simulation process without the elimination of the finest particles was a factor of about 17. The elimination of the finest particles gave additional speedup of a factor of at least 18.Therefore, the simulation of a grinding process can run at least 300 times faster than the conventional method in which a standard no binary search is employed and the smallest particles are tracked.     

Comparative study of two fluid model and dense discrete phase model for simulations of gas-solid hydrodynamics in circulating fluidized beds

Computational fluid dynamics(CFD)has become a valuable tool to study the complex gas-solid hydrodynamics in the circulating fluidized bed(CFB).Based on the two fluid model(TFM)under the Eulerian-Eulerian framework and the dense discrete phase model(DDPM)under the Eulerian-Lagrangian framework,this work conducts the comparative study of the gas-solid hydrodynamics in a CFB riser by these two different models.Results show that DDPM could be used to predict gas-solid hydrodynamics in the circulating fluidized bed,and there are differences between TFM and DDPM,especially in the radial distribution profiles of solid phase.Sensitivity analysis results show that the gas-solid drag model exhibits significant effects on the results for both the two models.The specularity coefficient and the restitution coefficient in the TFM,as well as the reflection coefficient and the parcel number in the DDPM,exhibit less impact on the simulated results.     

Anomalous behavior of fluid flow through thin carbon nanotubes

Theoretical and Computational Fluid Dynamics - Molecular dynamics simulation is used to study flow rate behavior of monoatomic fluid through carbon nanotube (CNT) against the pore diameter. All...     

Conservation of circulation in SPH for 2D free‐surface flows

In this paper, we study how accurately the Smoothed Particle Hydrodynamics (SPH) scheme accounts for the conservation and the generation of vorticity and circulation, in a low viscosity, weakly compressible, barotropic fluid in the context of free‐surface flows. We consider a number of simple examples to clarify the processes involved and the accuracy of the simulations. The first example is a differentially rotating fluid where the integration path for the circulation becomes progressively more complicated, whereas the structure of the velocity field remains simple. The second example is the collision of two rectangular regions of fluid. We show that SPH accurately predicts the time variation of the circulation as well as the total vorticity for selected domains advected by the fluid. Finally, a breaking wave is considered. For such a problem we show how the dynamics of the vorticity generated by the breaking process is captured by the SPH model. Copyright © 2012 John Wiley & Sons, Ltd.     

Periodic sticking motion in a two-degree-of-freedom impact oscillator

Periodic sticking motions can occur in vibro-impact systems for certain parameter ranges. When the coefficient of restitution is low (or zero), the range of periodic sticking motions can become large. In this work the dynamics of periodic sticking orbits with both zero and non-zero coefficient of restitution are considered. The dynamics of the periodic orbit is simulated as the forcing frequency of the system is varied. In particular, the loci of Poincaré fixed points in the sticking plane are computed as the forcing frequency of the system is varied. For zero coefficient of restitution, the size of the sticking region for a particular choice of parameters appears to be maximized. We consider this idea by computing the sticking region for zero and non-zero coefficient of restitution values. It has been shown that periodic sticking orbits can bifurcate via the rising/multi-sliding bifurcation. In the final part of this paper, we describe three types of post-bifurcation behavior which occur for the zero coefficient of restitution case. This includes two types of rising bifurcation and a border orbit crossing event.     

VIBRO-IMPACT DYNAMICS OF PIPE CONVEYING FLUID SUBJECTED TO RIGID CLEARANCE CONSTRAINT 1)

管道与间隙约束间的碰撞振动是工程输流管结构的一个重要动力学现象. 迄今,人们通常采用光滑的非线性弹簧来模拟管道与间隙约束之间的碰撞力,但这种光滑的碰撞力无法真实反映碰撞前后管道状态向量的非光滑传递特征. 本文基于非光滑理论建立了具有刚性间隙约束简支输流管的非线性碰撞振动模型,利用 Galerkin 法离散了无穷维的管道模型, 并引入恢复系数构造了碰撞前后管道各处状态向量的传递矩阵,运用四阶龙格库塔法分析了脉动内流激励下管道与刚性间隙约束的非光滑碰撞振动现象,着重讨论了刚性间隙约束参数对管道动态响应随流速脉动频率变化的影响规律,特别是碰撞振动的周期性运动规律. 研究结果表明,刚性约束间隙值对管道碰撞振动行为的影响较大,在某些脉动频率下管道会出现多周期和非周期性的运动形态,还可出现非光滑系统特有的黏滑现象. 此外,碰撞恢复系数对管道振动的影响也比较显著,较小的恢复系数值更容易使管道在大范围脉动频率区间出现复杂的非周期碰撞振动.     

具有刚性间隙约束输流管的碰撞振动

管道与间隙约束间的碰撞振动是工程输流管结构的一个重要动力学现象. 迄今,人们通常采用光滑的非线性弹簧来模拟管道与间隙约束之间的碰撞力,但这种光滑的碰撞力无法真实反映碰撞前后管道状态向量的非光滑传递特征. 本文基于非光滑理论建立了具有刚性间隙约束简支输流管的非线性碰撞振动模型,利用 Galerkin 法离散了无穷维的管道模型, 并引入恢复系数构造了碰撞前后管道各处状态向量的传递矩阵,运用四阶龙格库塔法分析了脉动内流激励下管道与刚性间隙约束的非光滑碰撞振动现象,着重讨论了刚性间隙约束参数对管道动态响应随流速脉动频率变化的影响规律,特别是碰撞振动的周期性运动规律. 研究结果表明,刚性约束间隙值对管道碰撞振动行为的影响较大,在某些脉动频率下管道会出现多周期和非周期性的运动形态,还可出现非光滑系统特有的黏滑现象. 此外,碰撞恢复系数对管道振动的影响也比较显著,较小的恢复系数值更容易使管道在大范围脉动频率区间出现复杂的非周期碰撞振动.      

Direct and indirect methods of calculating entropy generation rates in turbulent convective heat transfer problems

Computational fluid dynamics (CFD) solutions of turbulent convective heat transfer problems based on the mass, momentum and energy conservation principle provide all information to calculate the entropy production rate in such a transfer process. It can be determined in the post processing phase of a CFD calculation. Two methods are discussed in detail which can provide the information about the entropy production with different degrees of accuracy.     

Contributions of Computational Aeroacoustics to Jet Noise Research and Prediction

A brief review of recent progress in the field of computational aeroacoustics (CAA) is proposed. This paper is complementary to the previous reviews of Tam [(1995a) “Computational aeroacoustics: issues and methods”, AIAA J. 33(10), 1788–1796], Lele [(1997) “Computational Aeroacoustics: a review”, AIAA Paper 97–0018, 35th Aerospace Sciences Meeting and Exhibit, Reno, Nevada] and Glegg [(1999) “Recent advances aeroacoustics: the influence of computational fluid dynamics”, 6th International Congress on Sound and Vibration, Copenhagen, Danemark, 5–8 July, 43–58] on advances in CAA. After a short introduction concerning the current motivations of jet noise studies, connections between computational fluid dynamics (CFD) and CAA using hybrid approaches are discussed in the first part. The most spectacular advances are probably provided by the direct computation of jet noise, and some recent results are shown in the second part.     

计算气动弹性力学中的界面映射方法研究

非线性气动弹性体振动研究中,涉及到非线性的结构动力学和非线性的流体动力学耦合问题,在耦合边界上要满足两个系统的连续性相容条件,必须在边界处进行数据的交换。本文针对非线性气动弹性问题的计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法,在常体积转换法(CVT)的基础上,发展了一种耦合界面的数据映射矩阵(IMM)。该方法仅需要局部的网格信息,将耦合边界上载荷信息和位移信息的转换放在同一个映射矩阵中来处理,并且该矩阵可以通用求解CFD/CSD的耦合问题,克服了占用大量CPU时间和内存的弊端。最后将该界面映射方法应用于柔性大展弦比机翼的气动弹性计算和AGARD445.6机翼的颤振预测中,结果表明该方法能够高效、高精度地处理不同网格体系间的数据交换,并具有处理复杂非规则几何体信息转换的能力。     

Mathematical modeling of salt transport by coupled subsurface and surface water flows

A mathematical model is examined to describe the transport of salts by coupled flows of surface, soil, and subsoil waters for large-scale objects characterized by complex hydrogeologic conditions. Computational algorithms and computer programs developed for realization of the model are based on the use of finite-difference methods to distinguish between different physical processes and modeling regions. Results are presented from examples of calculations to illustrate characteristic features of the problem of mass transfer by coupled flows. Lavrent'ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 4, pp. 118–126, July–August, 1998.     

Godunov型CSPM方法

CSPM方法是传统SPH方法的改进,它有效地解决了传统SPH方法在求解含边界问题时的困难。然而最近研究表明,CSPM方法在求解含间断算例时却存在计算不稳定的问题。本文基于Godunov间断分解思想和CSPM方法提出了Godunov型CSPM方法。该方法成功地解决了CSPM在求解含间断流体学问题中遇到的困难,且不用引入人工粘性。为了验证算法的有效性,本文采用一维激波管算例对不同方法进行了对比计算。计算结果表明,本文提出的方法比以往SPH方法有更高的间断捕获精度。     

Recent advances of computational aeroacoustics

Computational aeroacoustics (CAA) is an interdiscipline of aeroacoustics and computational fluid dynamics (CFD) for the investigation of sound generation and propagation from various aeroacoustics problems. In this review, the foundation and research scope of CAA are introduced firstly. A review of the early advances and applications of CAA is then briefly surveyed, focusing on two key issues, namely, high order finite difference scheme and non-reflecting boundary condition. Furthermore, the advances of CAA during the past five years are highlighted. Finally, the future prospective of CAA is briefly discussed.     

基于流动仿真大数据应对旋涡-湍流的研究进展

文章以流体科学进入二十一世纪后，在大规模超算、云存储、数据通信和人工智能为支撑的大数据时代背景下，结合目前在复旦大学航空航天系所构建的热流体湍流直接数值仿真数据库，以及复旦大学团队近期与美国德州大学刘超群教授、上海理工大学蔡小舒教授以及国内水动力学杂志编辑部所合作开展的第三代涡识别技术研究，初步概念性地展示旋涡和湍流，特别是针对有工程实际背景和直接应用价值的壁湍流，在这两个流体力学关键基础议题上的最新认知，和基于大数据深度学习的相关湍流工程模拟实践成果．这些成果包括：(1) 基于第三代涡识别技术的尾迹湍流中的涡运动学和动力学探索；(2) 流-热统一完整的类-1、类-2湍流边界层壁面律构建；(3) 基于第三代涡识别量对Kolmogorov 幂次律的再认知；(4) 基于DNS统计数据和神经网络深度学习构建新型湍流封闭模型及RANS计算实践．通过这些成果展示，论证解决这两个基础流体科学议题的技术路径，进而促进流体及相关学科研究在现代大数据背景下取得实质性进展和突破，并惠及现代流体、气动、水利、动力和化工等工程领域．     

流体与结构相互作用问题的FE-SPH耦合模拟

应用有限元(FE)-光滑粒子流体动力学(SPH)耦合法模拟了具有自由表面的不可压流体与结构的相互作用问题.流体和结构分别采用SPH法和有限元法同时求解,两者在交界面处的相互作用通过接触算法进行处理.为了避免隐式计算压力,通过引入人工压缩率,将不可压流体近似为人工可压缩流体.采用FE-SPH耦合法对弹性板在随时间变化的水压作用下的变形以及倒塌水柱冲击弹性结构两个问题进行了模拟.模拟结果与实验结果以及其他已有数值结果符合良好,说明FE-SPH耦合法用于流体与结构相互作用问题的模拟是可行和有效的.     

攒尖四坡屋面风致雪漂移的数值模拟

为获得攒尖四坡屋面的风致雪漂移规律,基于欧拉-欧拉方法和风雪单向耦合假定,运用计算流体动力学软件,选取Mixture模型分别对立方体周边和高低屋面上的风致雪漂移运动进行数值模拟,将模拟结果与两者的实地观测数据对比,探讨分析数值风洞的关键技术和参数设置,验证数值模拟方法的合理性与可靠性。依据攒尖四坡房屋的使用功能要求,设计分析模型与分析工况,在试算的基础上对屋面进行分区。以风速5 m/s,7 m/s,9 m/s,11 m/s,13 m/s和15 m/s,风向角0°,15°,30°和45°以及屋面坡度25°,30°,34°,40°和45°为分析参数,对攒尖四坡房屋的120种工况进行数值模拟,得到屋面各分区侵蚀沉积的基本规律,提出可用于抗雪设计的屋面积雪分布系数。研究表明,风向角的改变会使屋面积雪分布状态发生较大程度的变化,风速和屋面坡度的变化对屋面整体积雪量有较大影响。     

A multiblock/multilevel mesh refinement procedure for CFD computations

A multiblock/multilevel algorithm with local refinement for general two‐ and three‐dimensional fluid flow is presented. The patched‐based local refinement procedure is presented in detail and algorithmic implementations are also presented. The multiblock implementation is essentially block‐unstructured, i.e. each block having its own local curvilinear co‐ordinate system. Refined grid patches can be put anywhere in the computational domain and can extend across block boundaries. To simplify the implementation, while still maintaining sufficient generality, the refinement is restricted to a refinement of the grid successively halving the grid size within a selected patch. The multiblock approach is implemented within the framework of the well‐known SIMPLE solution strategy. Computational experiments showing the effect of using the multilevel solution procedure are presented for a sample elliptic problem and a few benchmark problems of computational fluid dynamics (CFD). Copyright © 2001 John Wiley & Sons, Ltd.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics (CFD) modeling of the gas–solids two-phase flow in a circulating fluidized bed (CFB) riser is carried out. The Eularian–Eularian method with the kinetic theory of granular flow is used to solve the gas–solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory (Johnson & Jackson, 1987) with specularity coefficient and particle–wall restitution coefficient. The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operating conditions of the CFB riser.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics （CFD） modeling of the gas-solids two-phase flow in a circulating fluidized bed （CFB） riser is carried out. The Eularian-Eularian method with the kinetic theory of granular flow is used to solve the gas-solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory （Johnson ＆ Jackson, 1987） with specularity coefficient and particle-wall restitution coefficient.The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operatin~ conditions of the CFB riser.