Computational fluid dynamics-discrete element method simulation of stirred tank reactor for graphene production

Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.

A review of methods,applications and limitations for incorporating fluid flow in the discrete element method

The past decade has witnessed the substantial growth in research interests and progress on the subject of coupled hydro-mechanical processes in rocks and soils,driven mainly by the surge of research in unconventional hydrocarbon reservoirs and associated hazards.Many coupling techniques have been developed to include the effects of fluid flow in the discrete element method(DEM),and the techniques have been applied to a variety of geomechanical problems.Although these coupling methods have been successfully applied in various engineering fields,no single fluid/DEM coupling method is universal due to the complexity of engineering problems and the limitations of the numerical methods.For researchers and engineers,the key to solve a specific problem is to select the most appropriate fluid/DEM coupling method among these modeling technologies.The purpose of this paper is to give a comprehensive review of fluid flow/DEM coupling methods and relevant research.Given their importance,the availability or unavailability of best practice guidelines is outlined.The theoretical background and current status of DEM are introduced first,and the principles,applications,and advantages and disadvantages of different fluid flow/DEM coupling methods are discussed.Finally,a summary with speculation on future development trends is given.

Stabilization for Equal-Order Polygonal Finite Element Method for High Fluid Velocity and Pressure Gradient

This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that are governed by Stokes equations system.This technique is constructed by a local pressure projection which is extremely simple,yet effective,to eliminate the poor or even non-convergence as well as the instability of equal-order mixed polygonal technique.In this research,some numerical examples of incompressible Stokes fluid flow that is coded and programmed by MATLAB will be presented to examine the effectiveness of the proposed stabilised method.

MIXED COMPATIBLE ELEMENT AND MIXED HYBRID INCOMPATIBLE ELEMENT VARIATIONAL METHODS IN DYNAMICS OF VISCOUS BAROTROPIC FLUIDS

This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress can be naturally satisfied. The gene-rallzed variational principles with mixed hybrid incompatible finite elements are alsopresented and proved, and they can reduce the computation of incompatible elements indynamics of viscous barotropic flows.

An explicit finite element method for dynamic analysis in three-medium coupling system and its application

In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.

Methodology for Comparing Coupling Algorithms for Fluid-Structure Interaction Problems

The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.

NUMERICAL SIMULATION OF UNSTEADY-STATE UNDEREXPANDED JET USING DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD

A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.

Computational fluid dynamics simulations of respiratory airflow in human nasal cavity and its characteristic dimension study

To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.

RESEARCH ON SOLID-LIQUID COUPLING DYNAMICSOF PIPE CONVEYING FLUID

On the basis of Hamilton principle. the equation of sonlid-liquid coupling vibration of pipe conveying fluid is deduced. An asymmetrical sonlid-liquid coupling damp matrix and a symmetrical solid-liquid coupling Stiffness matrix are obtained. Using QR method , pipe’s nature frequencies are calculated. The curves of the first four orders of natural frequency-flow velocity of pipe waw given .The influence of flowing velocity ,pressure, solid-liquid coupling damp and solid-liquid coupling stiffness on natural frequency are discussed respectively.The dynamic respondence of the pipes for stepload with different flow velocity are calculated by Newmark method .It is found that,with the flow velocity increased, the nature frequency of the pipes reduced, increased,reduced again and so on.

Structural Parameter Analyses on Rotor Airloads with New Type Blade-Tip Based on CFD/CSD Coupling Method

For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.

Application of the Hierarchical Functions Expansion Method for the Solution of the Two Dimensional Navier-Stokes Equations for Compressible Fluids in High Velocity

This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as SUPG (Streamline Upwind Petrov-Galerkin), applied with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed simulations are performed for some cases whose theoretical solutions are known. These cases are the following: continuity test, stability and convergence test, temperature step problem, and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations with the proposed method were good both qualitatively and quantitatively when compared with the theoretical solutions. This allows concluding that the objectives of this work are reached.

CVBEM and FVM Computational Model Comparison for Solving Ideal Fluid Flow in a 90-Degree Bend

While finite volume methodologies (FVM) have predominated in fluid flow computations, many flow problems, including groundwater models, would benefit from the use of boundary methods, such as the Complex Variable Boundary Element Method (CVBEM). However, to date, there has been no reporting of a comparison of computational results between the FVM and the CVBEM in the assessment of flow field characteristics. In this work, the CVBEM is used to develop a flow field vector outcome of ideal fluid flow in a 90-degree bend which is then compared to the computational results from a finite volume model of the same situation. The focus of the modelling comparison in the current work is flow field trajectory vectors of the fluid flow, with respect to vector magnitude and direction. Such a comparison is necessary to validate the development of flow field vectors from the CVBEM and is of interest to many engineering flow problems, specifically groundwater modelling. Comparison of the CVBEM and FVM flow field trajectory vectors for the target problem of ideal flow in a 90-degree bend shows good agreement between the considered methodologies.

基于微观流固耦合的超深层致密砂岩气藏应力敏感性分析

超深层致密砂岩气藏具有强应力敏感性,目前常用的研究方法包括压力脉冲实验法和实时在线CT扫描法两种,压力脉冲实验法不能揭示产生应力敏感性的微观机理,而实时在线CT扫描法也无法模拟深部地层高压、高应力的条件。为解决应力敏感性实验研究的不足,基于离散单元法与管道网络模型建立了微观流固耦合算法,编制了模拟器,并对模拟器力学计算和流固耦合模块的正确性进行了验证,分析了应力大小、加载方向和孔隙压力对岩心渗透率的影响,最后从微观上揭示了超深层致密砂岩气藏的应力敏感性机理。研究结果表明:①应力通过增加与之垂直方向上喉道两侧的法向压力,减小喉道的水力半径,进而降低储层的渗透率;②较高的孔隙压力能够阻碍岩石颗粒在应力作用下的移动,从而减缓了孔隙和喉道的变形,使模型保持较高的渗透率;③致密砂岩气藏的渗透率受到应力和地层压力的共同控制,并且具有各向异性,在垂直于最小主应力方向上形成渗透率较大的优势通道;④异常高压阻碍了地应力的压实作用,有利于保护储层孔隙,使地层有较好的储集性能和较高的渗透率。结论认为,根据离散元法结合孔隙网络模型建立的流固耦合方法可为理解超深层致密砂岩应力敏感性提供理论参考,并为超深层致密砂岩气藏的科学高效开发提供指导。

基于CFD-DEM耦合的梯级溜槽的设计与分析

针对传统物料转运过程中溜槽和输送带磨损严重、出口处粉尘浓度过高的问题,建立含臂架的梯级溜槽几何模型,采用基于计算流体力学与离散单元法(CFD-DEM)耦合的数值模拟方法,分析了臂架对转运溜槽的磨损以及对其出口处粉尘排放浓度的影响。仿真结果表明:含臂架的梯级溜槽可以有效控制物料流的速度和方向,降低对溜槽内表面的冲击磨损,降低出口处的粉尘量。

氦冷固态包层氚增殖球床气体和粉末流动特性的数值研究

基于计算流体力学和离散元耦合分析,研究了氦冷固态增殖包层球床中吹扫氦气的速度分布和破碎的正硅酸锂粉末速度随时间的变化及粉末数量随氦气入口流速变化。模拟结果表明,在吹扫氦气流动方向上,用入口流速归一化的氦气速度分布与入口流速无关,在靠近球床边缘速度变化较大,在球床中部速度变化较小;用入口流速归一化的粉末平均速度与入口流速关系不大,会随吹扫时间逐渐趋于稳定;对小粒径粉末,吹扫气体流速越大,越容易被吹出球床,对大粒径粉末,因其本身体积过大,易在球床中形成堵塞。

多级水力压裂应力阴影效应的数值分析

多级水力压裂作为储层改造的主要技术手段,广泛用于提高致密油气藏的气体产量。裂缝在扩展过程中的应力阴影效应对于有效连通储层的天然裂缝并形成有利于油、气体流动的复杂裂缝网络至关重要。本研究采用了自适应的有限元-离散元方法,通过改进网格自动细化和识别多个裂缝扩展模拟含天然裂缝储层的多级水力压裂。数值模型中考虑了水力裂缝、天然裂缝和微观孔隙结构的相互作用,整合了非线性的Carter滤失模型来描述多级压裂过程中压裂液的滤失效应。引入了理想平行板流的支撑剂运输方程,并采用了达西定律来分析裂缝网络中的流体渗透效应。通过比较无裂缝均质模型和天然裂缝模型的裂缝网络和流体流动,评估了天然裂缝对多级压裂裂缝行为和气体产量的影响。研究为确定和优化致密气藏中压裂簇间距提供了一种新的方法。

恶劣海况下船舶砰击颤振响应特性数值计算与试验研究

[目的]针对恶劣海况下船舶所受砰击颤振响应现象,探究船舶非线性波浪载荷与瞬态高幅值砰击载荷的耦合作用。[方法]采用计算流体动力学与有限元方法(CFD-FEM)相结合的双向流固耦合方法对S175集装箱船进行数值仿真计算,并与试验结果及切片理论计算结果进行对比验证;采用分段变截面弹性龙骨梁模型开展船舶的砰击颤振特性模型试验,基于CFD-FEM双向流固耦合方法开展船艏砰击载荷及高频非线性砰击颤振响应特性分析,并与模型试验结果进行对比验证。[结果]结果显示,波浪砰击载荷对船艏颤振响应的影响不可忽视,6级海况下由砰击颤振诱发的二阶高频成分分量占低频波浪弯矩的59.86%。[结论]采用基于CFD-FEM的双向流固耦合方法可准确计算船首砰击颤振响应;在高海况下船舶所受非线性波浪载荷及结构动态响应易受船首瞬态砰击载荷的影响,在船舶结构设计与安全评估中需考虑高频砰击颤振的情况。

调谐液柱阻尼器-结构系统风致振动响应的CFD/CSD耦合分析方法

针对调谐液柱阻尼器(tuned liquid column damper, TLCD)难以建立其精确的非线性理论分析模型,且其力学性能试验成本高和耗时长等问题,首先采用计算流体动力学(CFD)数值模拟方法,对TLCD系统的力学性能和动力特征进行仿真模拟,在此基础上进一步提出了基于计算流体动力学/计算结构动力学(CFD/CSD)耦合分析方法,求解带TLCD系统的高层建筑结构的风致动力响应。通过开展某一TLCD系统在特定底部激励下的力学性能和动力特性试验,得到其内液体晃荡的自由液面波高和晃动力时程,验证了CFD数值模拟方法可以准确地分析TLCD水箱内液体的非线性晃动特征。随后对风工程领域广泛采用的76层建筑结构振动控制Benchmark模型,假设其顶部设置TLCD系统时主体结构在三种风速重现期(10、50和100年)风速对应的横风向动力风荷载激励下的风致控制效率,采用提出的CFD/CSD耦合分析方法,进行了数值仿真模拟分析。耦合分析结果表明,TLCD系统对Benchmark模型的风致加速度、速度和位移响应均有一定的控制效果,对加速度响应的控制效果要优于对位移响应的控制效果。该研究方法可为复杂TLCD系统对高层建筑的风振控制分析提供有效的参考。

龙羊峡重力拱坝模态识别与材料参数反演

基于龙羊峡重力拱坝实测振动数据,利用随机子空间法进行了大坝模态参数识别,建立了坝体-地基-库水系统的三维流固耦合有限元模型。根据识别出的一阶频率,分别考虑均质地基及分层地基,反演了坝体与基岩材料参数,并通过有限元软件ABAQUS进行了大坝的动力特征分析。结果表明:龙羊峡重力拱坝的一阶频率为2.778~3.196 Hz,反演得到坝体的弹性模量约为50 GPa~52 GPa;对于均值地基,其弹性模量的反演值为30 GPa;对于分层地基,2 490 m高程以下区域的弹性模量为34 GPa;实测振动数据反演得到的坝体弹性模量和地基弹性模量高于其设计值;龙羊峡重力拱坝的前3阶振型分别为正对称、反对称与正对称,其1阶频率和2阶频率相差较大,与重力坝自振频率分布特征相似。

大变幅加卸载下特厚煤层底板断层突水机理模拟研究

特厚煤层采场空间大、扰动范围广,强烈大变幅荷载易导致底板断层破裂加剧并诱发水害。数值模拟是揭示特厚煤层底板断层活化与突水机理的重要方法,准确反映大变幅加卸载下岩体破裂与裂隙水耦合特征是其合理性的关键。构建损伤变量与塑性应变、应力之间的相关关系,得到完整岩块的拉、加压卸载损伤演化方程;以平方拉剪应力与Benzeggagh-Kenane为初始、完全断裂准则,建立塑性位移与强度劣化关系,建立加卸载韧性断裂本构关系;基于实验数据建立贯通裂隙加卸载剪切本构关系。以基本方程与状态方程为基础,结合浸没边界方法,形成裂隙岩体水力学模拟理论。由此编制流体动力学-有限离散元CFD-FDEM耦合程序,模拟研究特厚煤层底板断层活化突水过程。结果表明:CFD-FDEM耦合程序可数值实现特厚煤层底板断层从(准)连续体到离散体转化,以及断层带裂隙水运移过程。底板断层采动破坏包络线呈W形,最深位于断层及其上盘(48.6m),最浅位于断层下盘(23m)。特厚煤层采场底板断层及其上盘受到较大超前集中应力,而后在采空区内大幅卸载,导致该位置出现显著二次破坏,并形成主要导水通道。研究成果为特厚煤层工作面底板断层水害防治提供理论支撑。