带自由表面水流与离散颗粒耦合模型研究

基于考察泥沙运动的细观行为特征,采用离散单元法(DEM)模拟泥沙颗粒运动,结合带自由表面的水动力学计算模型,建立了CFD-DEM耦合数值模型。计算程序开发基于Fortran语言来实现。耦合模型中实现了硬球模型和软球模型两种颗粒碰撞模型,应用范围较广。作为自由表面水流与泥沙颗粒流数值模型的初步研究,在模型建立的基础上,对模型做了基本的验证。分别通过单颗粒静水沉降和混合颗粒群分选两个计算工况,验证了模型的正确性及模拟精度。该耦合模型可进一步丰富带自由表面水流条件下泥沙运动的研究手段。     

含刚性线夹杂及裂纹的各向异性压电材料耦合场分析

采用各向异性弹性力学中Ｓｔｒｏｈ方法对含刚性线夹杂及裂纹的无限大各向异性压电材料耦合的弹性场和电场进行了分析。并得到夹杂和基体界面间耦合场的实型显函表达式及夹杂尖端的１／２阶奇异性。     

流固耦合破坏分析的多分辨率PD-SPH方法

流固耦合破坏是一类涉及结构变形与破坏以及复杂自由表面现象的强非线性力学问题.结合近场动力学(peridynamics, PD)与光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)各自的优势并考虑其计算效率问题,提出一种适用于分析流-固耦合破坏问题的多分辨率PD-SPH混合方法.分别采用SPH和PD方法以不同的空间和时间分辨率对流体和结构进行离散与求解,利用具有与流体粒子相同光滑长度的虚粒子处理流-固界面,以高精度满足界面边界条件.通过两个经典算例:液柱静压力下弹性板的变形和溃坝流体冲击弹性闸门的变形问题,表明提出的多分辨率PD-SPH方法兼具较高的计算精度和计算效率;对含裂缝的Koyna重力坝水力劈裂问题进行模拟计算,所得裂缝扩展路径与文献结果吻合,说明该方法适用于涉及结构破坏的流固耦合问题仿真.最后尝试采用该方法进行流体冲击作用下含裂纹混凝土板崩塌过程数值仿真,准确描述混凝土板的断裂破坏和全过程中的流体运动.多分辨率PD-SPH混合方法或可为流-固耦合作用下的结构损伤破坏仿真提供一种新选择.     

颗粒体流-固特性的实验研究方法评述

颗粒材料是一种常见的复杂多体物质,其形态各不相同,颗粒之间通过相互作用使其表现出十分复杂的流体或固体特性。因此结合理论、实验及模拟的方法对颗粒材料进行研究,可以进一步发现颗粒的物理规律并对其进行表征,同时也有助于学者对颗粒材料的特性形成更加系统的认识。本文对颗粒材料在不同物质状态下表现出的不同流态特性与多尺度特点进行简要介绍,同时回顾了颗粒材料流-固特性的研究现状与2种状态下的实验研究方法,概括了相关测量方法的优缺点。     

颗粒流动力学及其离散模型评述

颗粒流是由众多颗粒组成的具有内在相互作用的非经典介质流动. 自然界常见颗粒流都是密集流, 颗粒间接触形成力链, 诸多力链相互交接构成支撑整个颗粒流重量和外载荷的网络, 其局部构型及强度在外载荷下演化, 是颗粒流摩擦特性和接触应力的来源.本文介绍球形颗粒间无粘连作用时的Hertz法向接触理论和Mindlin-Deresiewicz切向接触理论. Campbell依据是否生成较为稳定的力链把颗粒流分为弹性流和惯性流两大类, 其中弹性-准静态流和惯性-碰撞流分别对应准静态流和快速流, 作为两种极端流动情况通常处理成连续体, 分别采用摩擦塑性模型和动理论予以描述, 但是表征接触力链的颗粒弹性参数并不出现这两个模型和理论框架中, 如何进一步考虑颗粒弹性参数将非常困难. 目前离散动力学方法逐渐成为复现其复杂颗粒流动现象、提取实验不可能获得的内部流动信息进而综合起来探索颗粒流问题的一种有效工具, 其真实性强于连续介质理论的描述. 软球模型对颗粒间接触力简化处理, 忽略了切向接触力对法向接触力及其加载历史的依赖, 带来了法向和切向刚度系数如何标度等更艰难的物理问题, 但由于计算强度小而广泛应用于工程问题中. 硬球模型不考虑颗粒接触变形, 因而不能描述颗粒流内在接触应变等物理机理, 仅适用于快速颗粒流, 这不仅仅是由于两体碰撞的限制. 因此基于颗粒接触力学的离散颗粒动力学模型是崭新的模型,适用于准静态流到快速流整个颗粒流态的模拟, 可以细致考虑接触形变及接触力的细节,建立更为合理的颗粒流本构关系, 进而有力的促进颗粒流这一非经典介质流动的研究.      

应力历史对颗粒材料强度和变形特征影响的离散元模拟

主要关注了颗粒材料前期所受的应力历史对其后期宏观力学响应的影响。该应力历史由一段等比例加载应力路径以及卸载垂直方向应力至与水平方向围压相同的卸载段描述。具体工作为:基于PFC2D双轴压缩数值实验,调查了应力历史对颗粒样本的强度、变形特征、细观参量如组构的影响,得出颗粒样本的偏应力-应变、体积应变曲线及其发生破坏时的名义应变云图。数值结果表明:高低围压下样本分别发生剪切破坏和弥散破坏,高围压下弹性阶段的刚度受应力历史影响较大,而低围压下样本在刚进入塑性至应力峰值点阶段的弹塑性刚度变化较为明显。随着应力历史中等比例加载系数的增大,剪胀加快,变形局部化范围有所不同;另一方面,颗粒形状的不规则性会增强颗粒材料的各向异性,导致样本强度更高。     

基于Level Set方法的气-液-固三相流动模型与模拟

采用基于Level Set方法与离散颗粒模型相结合的方法,建立了一个用于描述气液固三相流动的新模型.模型耦合了颗粒与气泡、颗粒与液相以及气泡与液相之间的相互作用.应用该模型对液固悬浮液中的典型现象--气泡的单孔及多孔形成过程以及颗粒夹带现象进行了三维模拟,检验了其可行性.并进一步研究了颗粒的存在对气泡的形成与上升过程的影响以及气泡诱导的液相流动对颗粒行为的影响.研究结果表明,所提出的模型能够真实地预测三相流中气泡与颗粒分散相的特征,为研究多尺度的三相流动提供了一种新途径.     

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

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

气液两相瞬变流的流固耦合模型研究

传统的气液两相流瞬态分析和管道动力响应计算是分开的，存在一定的缺陷．针对石油工业中常见的多相混输问题，介绍了常见的气液两相流瞬态模型和流固之间存在的耦合机理，在不作薄壁管假设的前提下推导出了气液两相瞬变流的流固耦合模型．与现有相似模型的对比分析表明，这一模型比较全面地考虑了流体和管道的特性以及不同的耦合形式，可以适应实验和／或仿真研究的需要．     

基于细观方向平均模型的颗粒材料宏观Cosserat连续体本构关系

提出了基于细观微-方向模型（Micro—Directional Model）的宏观Cosserat连续体本构关系。在细观尺度上考虑颗粒旋转自由度及接触力矩,微结构的影响通过接触分布函数体现。给出均质各向同性Cosserat连续体模型弹性常数的细观参数表达式,并建议了二维情况下内尺度参数的细观力学表达式。对颗粒材料宏观行...     

A novel method to calculate the pressure interaction between dust and fluid using SPH

Unstable behavior of smoothed particle hydrodynamics (SPH) dust particles, such as clumping or fingering under certain conditions, has been reported by several researchers who have conducted studies on dusty fluid SPH. The simulation results in this study show that this instability is numerical, and the instability is mainly attributable to the ill‐interpolated pressure gradient in the interaction term between 2 phases. In this paper, we introduce a new method to calculate the pressure force interaction term between dust and fluid particles. The key idea is to first interpolate the pressure gradient at SPH fluid particles and then use the values to calculate the pressure gradient at SPH dust particles, in a consecutive manner. To compare the new method with the existing method, we first conducted an interpolation of pressure gradient at hydrostatic equilibrium under gravity to estimate any error. The results show that the new method is more accurate. We then conducted additional numerical tests, namely, dust‐liquid counterflow, sedimentation in a confined tank, and sedimentation in the presence of turbulence. The unphysical unstable behavior of SPH dust particles such as clumping or fingering was significantly reduced in the new method. The results also show that the instability becomes more significant when using the existing method especially for the case when simulating a flow with relatively high concentration of dust or for the case in which inertia dominates the dynamics of dust particles. Especially, in those cases, the existing method should be avoided, and the newly proposed method is highly recommended.     

Superquadric DEM-SPH coupling method for interaction between non-spherical granular materials and fluids

The research on the coupling method of non-spherical granular materials and fluids aims to predict the particle–fluid interaction in this study. A coupling method based on superquadric elements is developed to describe the interaction between non-spherical solid particles and fluids. The discrete element method (DEM) and the smoothed particle hydrodynamics (SPH) are adopted to simulate granular materials and fluids. The repulsive force model is adopted to calculate the coupling force and then a contact detection method is established for the interaction between the superquadric element and the fluid particle. The contact detection method captures the shape of superquadric element and calculates the distance from the fluid particle to the surface of superquadric element. Simulation cases focusing on the coupling force model, energy transfer, and large-scale calculations have been implemented to verify the validity of the proposed coupling method. The coupling force model accurately represents the water entry process of a spherical solid particle, and reasonably reflects the difference of solid particles with different shapes. In the water entry process of multiple solid particles, the total energy of the water entry process of multiple solid particles tends to be stable. The collapse process of the partially submerged granular column is simulated and analyzed under different parameters. Therefore, this coupling method is suitable to simulate fluid–particle systems containing solid particles with multiple shapes.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

A coupled WC-TL SPH method for simulation of hydroelastic problems

A coupled weakly compressible (WC) and total Lagrangian (TL) smoothed particle hydrodynamics (SPH) method is developed for simulating hydroelastic problems. The fluid phase is simulated using WCSPH method, while the structural dynamics are solved using TLSPH method. Fluid and solid components of the method are validated separately. A sloshing water tank problem is solved to test the WCSPH method while oscillation of a thin plate and large deformation of a cantilever beam are simulated to test the TLSPH method. After validating each component, the coupled WC-TL SPH scheme is used to simulate two benchmark hydroelastic problems. The first test case shows the evolution of water column with an elastic boundary gate, and the second one investigates the breaking water column impact on elastic structures. The agreement between WC-TL SPH results and literature data shows the ability of the proposed method in simulating hydroelastic phenomena.     

考虑破碎影响的颗粒材料亚塑性模型及应变局部化模拟

在基于2ndP-K应力率的亚塑性模型基础上,通过引入一个能够考虑颗粒破碎影响的孔隙比-平均压力临界状态方程,形成了一个能够模拟颗粒破碎影响的颗粒材料亚塑性模型,数值算例考查了颗粒破碎对应变局部化模式及位移-承载曲线的影响,结果表明,所建议模型具有模拟破碎对颗粒材料应变局部影响的良好性能。     

A consistent incompressible SPH method for internal flows with fixed and moving boundaries

An improved incompressible smoothed particle hydrodynamics (ISPH) method is presented, which employs first‐order consistent discretization schemes both for the first‐order and second‐order spatial derivatives. A recently introduced wall boundary condition is implemented in the context of ISPH method, which does not rely on using dummy particles and, as a result, can be applied more efficiently and with less computational complexity. To assess the accuracy and computational efficiency of this improved ISPH method, a number of two‐dimensional incompressible laminar internal flow benchmark problems are solved and the results are compared with available analytical solutions and numerical data. It is shown that using smaller smoothing lengths, the proposed method can provide desirable accuracies with relatively less computational cost for two‐dimensional problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.     

Coupling of discrete-element method and smoothed particle hydrodynamics for liquid-solid flows

Particle based methods can be used for both the simulations of solid and fluid phases in multiphase medium, such as the discrete-element method for solid phase and the smoothed particle hydrodynamics for fluid phase. This paper presents a computational method combining these two methods for solid-liquid medium. The two phases are coupled by using an improved model from a reported Lagrangian-Eulerian method. The technique is verified by simulating liquid-solid flows in a two-dimensional lid-driven cavity.     

一种改进的MPS粒子作用模型

从分析角度出发,为移动粒子半隐式法(MPS)提供了一组新的粒子作用模型,该模型在粒子局部坐标系下进行泰勒展开,同步获取各算子的粒子近似式,确保了离散精度的一致性.方法对于随机分布的粒子和边界附近的粒子同样适用.文中利用三个算例对该方法进行了检验,计算结果与分析解符合较好.最后,运用该方法成功模拟了自由面流动的二维溃坝问题.