砌体结构破坏倒塌过程的离散单元模拟

基于传统颗粒离散单元模型（DEM）,开发了一种绑定式离散单元模型（BTDEM）用于砌体结构崩塌机制的研究,并对一个典型砌体结构的破坏倒塌过程进行数值模拟。结果证实BTDEM方法能够详细地模拟出砌体建筑结构的崩塌过程,同时还可以模拟出砌块的破碎,是一种普适性较强的模拟砌体结构建筑的高效数值方法。BTDEM方法具有计算速度...     

基于刚（粘）塑性流动理论的自然单元法研究

将自然单元法与刚（粘）塑性流动理论相结合,对自然单元法在金属塑性成形过程数值模拟中的应用进行了研究。采用基于Voronoi图和Delaunay三角化结构的Non-Sibsonian插值方法构造近似速度场向量,实现无网格方法中速度边界条件的直接精确施加,提出了基于刚（粘）塑性流动理论的无网格自然单元法。运用不完全广义变分...     

入水后爆炸问题的数值模拟技术

应用各种数值计算方法及计算处理技术,编制程序实现了对弹体入水后爆炸问题全过程的数值模拟.其中弹体与水,爆轰产物气体与弹壳之间的相互作用通过流固耦合技术来描述;水面与空气、爆轰产物与水、空气之间的相互作用采用VOF方法(Volume-of-Fluid)来描述;采用了刚体-柔体转换、单元失效删除等计算技术以更高效、更好地模...     

磁流体压缩管道热电磁流动特性研究

利用磁流体五波模型对低磁雷诺数下压缩管道中磁流体流动进行数值模拟。该模型由带有电磁作用强制项的Navier-Stokes方程组与电势Poisson方程组成,数值格式分别采用严格保证熵条件的熵条件格式和中心差分格式。数值模拟对不同磁作用数下的不同几何外形管道进行数值模拟研究,结果表明在磁流体压缩管道中,由于发生器模式提取...     

多尺度模拟与计算研究进展

简要介绍了多尺问题与研究方法.重点论述了两类常见多尺度问题的模拟计算方法与研究进展,分析了各自的优缺点和使用范围.对现有研究的局限性和存在的问题进行分析,指出了进一步研究多尺度模拟与计算的必要性.介绍了求解含有孤立缺陷问题的非局部准连续体法,MAAD等方法以及求解基于微观模型本构模拟问题的局部连续体法、HMM等方法.文章对多尺度模拟与计算的前景进行展望,提出了一些亟待解决的问题.     

金属-碳管复合结构的计算力学研究

系统地研究了金属-碳纳米管复合结构的力学行为,考察了单轴压缩载荷作用下填充管的临界屈曲应变对管内金属原子数目的依赖性,分析了管的几何特征,包括管径、管长及手性,对填充管变形与力学行为的影响,并与连续体力学模型的预测进+行了对比分析.本文的研究结果对金属-碳管复合结构的理论研究和工程应用都具有较好的指导意义.     

钝感炸药爆轰波传播和驱动的数值模拟研究

为实现对非理想爆轰波传播和爆轰驱动问题的精确计算,研究了与贴体坐标下Level Set方法相结合的爆轰冲击波动力学(DSD)计算方法;根据Hamilton-Jacobi方程的Godunov差分格式,提出了非正交的贴体坐标系中Level Set函数方程的差分格式及其相应的数值方法;并将该DSD方法与二维流体动力学程序TD...     

重大工程引起的珠江河口潮波逆向不稳定传播现象

采用模式ccost-2d对50年前白藤堵海工程及10年之后泥湾门区域的水动力特征变化进行了模拟和分析.结果可见,堵海工程造成泥湾门水道白蕉站潮差变小,高低潮时刻滞后2h以上,低潮水位抬升超过0.5m,引起严重的内涝问题.通过数值解对整个河网区的流动过程分析发现,由于堵海工程堵塞了泥湾门水道直接通向大海的通道,仅留的鸡啼...     

Two-component Brownian coagulation: Monte Carlo simulation and process characterization

The compositional distribution within aggregates of a given size is essential to the functionality of composite aggregates that are usually enlarged by rapid Brownian coagulation. There is no analytical solution for the process of such two-component systems. Monte Carlo method is an effective numerical approach for two-component coagulation. In this paper, the differentially weighted Monte Carlo method is used to investigate two-component Brownian coagulation, respectively, in the continuum regime, the free-molecular regime and the transition regime. It is found that (1) for Brownian coagulation in the continuum regime and in the free-molecular regime, the mono-variate compositional distribution, i.e., the number density distribution function of one component amount (in the form of volume of the component in aggregates) satisfies self-preserving form the same as particle size distribution in mono-component Brownian coagulation; (2) however, for Brownian coagulation in the transition regime the mono-variate compositional distribution cannot reach self-similarity; and (3) the bivariate compositional distribution, i.e., the combined number density distribution function of two component amounts in the three regimes satisfies a semi self-preserving form. Moreover, other new features inherent to aggregative mixing are also demonstrated; e.g., the degree of mixing between components, which is largely controlled by the initial compositional mass fraction, improves as aggregate size increases.     

Validation of EMMS-based drag model using lattice Boltzmann simulations on GPUs

Interphase momentum transport in heterogeneous gas–solid systems with multi-scale structure is of great importance in process engineering. In this article, lattice Boltzmann simulations are performed on graphics processing units (GPUs), the computational power of which exceeds that of CPUs by more than one order of magnitude, to investigate incompressible Newtonian flow in idealized multi-scale particle–fluid systems. The structure consists of a periodic array of clusters, each constructed by a bundle of cylinders. Fixed pressure boundary condition is implemented by applying a constant body force to the flow through the medium. The bounce-back scheme is adopted on the fluid–solid interfaces, which ensures the no-slip boundary condition. The structure is studied under a wide range of particle diameters and packing fractions, and the drag coefficient of the structure is found to be a function of voidages and fractions of the clusters, besides the traditional Reynolds number and the solid volume fractions. Parameters reflecting multi-scale characters are, therefore, demonstrated to be necessary in quantifying the drag force of heterogeneous gas–solid system. The numerical results in the range 0.1 ≤ Re ≤ 10 and 0 < ? < 0.25 are compared with Wen and Yu's correlation, Gibilaro equation, EMMS-based drag model, the Beetstra correlation and the Benyahia correlation, and good agreement is found between the simulations and the EMMS-based drag model for heterogeneous systems.