Large deformation simulations of geomaterials using moving particle semi-implicit method

Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.

Comparison of Construction Method for DEM Simulation of Ellipsoidal Particles

Discrete element model was developed to simulate the ellipsoidal particles moving in the moving bed. Multi-element model was used to describe a ellipsoidal particle, the contact detection algorithm of ellipsoidal particle was developed, and both contact force and gravity force were considered in the models. The simulation results were validated by our experiment. Three algorithms for representing an ellipsoidal particle were compared in macro and micro aspects. The results show that there exists big difference in the microscopic parameters such as kinetic energy, rotational kinetic energy, deformation, contact force and collision number which leads to the difference of macroscopic parameters. The relative error in the discharge rate and tracer particle position is the largest between 3-tangent-element representation and experimental results. The flow pattern is similar for the 5-element and 3-intersection representations. The only difference is the discharge rate of 5-element representation is larger than the experimental value and that of the 3-intersection representation has the contrary result. Finally the 3-intersection- element reoresentation is chosen in the simulation due to less comouting time than that of the 5-element renresentation.

Effect of moving baffle on average velocity and mixing of binary particles in rotating drums

Adding a moving baffle to the drum is a new way to enhance the motion and mixing of particles in rotating drums.To obtain its influence on binary particles,horizontal rotating drums provided with a moving baffle were investigated by discrete element method(DEM).AtΩ=15 r/min,increasing the length of moving baffle can increase the fluctuation amplitude of average particle velocity.AtΩ=60 r/min,the influence of the moving baffle on the average velocity fluctuation tends to be more random.At both rotational speeds,the moving baffle causes the average particle velocity to fluctuate more sharply.The moving baffle can enhance particle mixing.AtΩ=15 r/min,the moving baffle with length ofδ=1/3 can best enhance particle mixing.However,atΩ=60 r/min,only the moving baffle with a specific length(δ=1/4)can enhance mixing.This basic research has a positive reference value for the application of the moving baffle in industry.

Moving particle semi-implicit simulation on the molten Wood's metal downward relocation process

In the case of a severe accident involving nuclear reactors,an important aspect that should be considered is the leakage of molten material from the inside of the reactor into the environment.These molten materials damage other reactor components,such as electrical tubes,grid plates and core catchers.In this study,the moving particle semi-implicit(MPS)method is adopted and improved to analyze the twodimensional downward relocation process of molten Wood’s metal as a representation of molten material in a nuclear reactor.The molten material impinges the Wood’s metal plate(WMP),which is mounted on a rigid dummy stainless steel in a cylindrical test vessel.The breaching process occurs because of heat transfer between the molten material and WMP.The formed breach areas were in good agreement with the experimental results,and they showed that the molten Wood’s metal spread above the WMP.The solid WMP fraction decreased with time until it reached the termination time of the simulation.The present results show that the MPS method can be applied to simulate and analyze the downward relocation process of molten material in the grid plate of a nuclear reactor.

Simulation of Solid Particle Interactions Including Segregated Lamination by Using MPS Method

A new MPS(Moving Particle Semi-implicit)method is developed to simulate the behaviors and interactions of multiple fine solid particles as a continuum.As fluid particles are affected by viscosity,so solid particles are affected by friction.The solid particle dynamics for landslides,dumping,and gravity sorting etc.which can be difficult to simulate using conventional MPS methods,are modeled in this paper using the developed multi-solid-particle MPS method that benefits from drawing comparisons with the corresponding fluid particle behaviors.The present MPS results for dumping solid particles are verified against the corresponding DEM(Discrete Element Method)results.The shape and angle of repose for solid particles are shown to be highly dependent on the friction coefficient between grains.The peculiar phenomenon of segregated lamination(gravity sorting)among grains of different densities has been successfully reproduced using the multi-solid-particle MPS method.Lamination quality is found to be dependent on the densities and frictional coefficients of the constituent particles.The behavior of heterogeneous mixtures of multiple solid and liquid particles are also compared and discussed.This newly developed tool offers a window into the physical dynamics of sedimentology that the broader geoscience community might benefit from.

Investigation of Numerical Conditions of Moving Particle Semi‑implicit for Two‑Dimensional Wedge Slamming

The sensitivity of moving particle semi-implicit(MPS)simulations to numerical parameters is investigated in this study.Although the verifcation and validation(V&V)are important to ensure accurate numerical results,the MPS has poor perfor-mance in convergences with a time step size.Therefore,users of the MPS need to tune numerical parameters to ft results into benchmarks.However,such tuning parameters are not always valid for other simulations.We propose a practical numerical condition for the MPS simulation of a two-dimensional wedge slamming problem(i.e.,an MPS-slamming condition).The MPS-slamming condition is represented by an MPS-slamming number,which provides the optimum time step size once the MPS-slamming number,slamming velocity,deadrise angle of the wedge,and particle size are decided.The simulation study shows that the MPS results can be characterized by the proposed MPS-slamming condition,and the use of the same MPS-slamming number provides a similar fow.

On the limitations of linear beams for the problems of moving mass-beam interaction using a meshfree method

This paper deals with the capabilities of linear and nonlinear beam theories in predicting the dynamic response of an elastically supported thin beam traversed by a moving mass. To this end, the discrete equations of motion are developed based on Lagrange＇s equations via reproducing kernel particle method （RKPM）. For a particular case of a simply supported beam, Galerkin method is also employed to verify the results obtained by RKPM, and a reasonably good agreement is achieved. Variations of the maximum dynamic deflection and bending moment associated with the linear and nonlinear beam theories are investigated in terms of moving mass weight and velocity for various beam boundary conditions. It is demonstrated that for majority of the moving mass velocities, the differences between the results of linear and nonlinear analyses become remarkable as the moving mass weight increases, particularly for high levels of moving mass velocity. Except for the cantilever beam, the nonlinear beam theory predicts higher possibility of moving mass separation from the base beam compared to the linear one. Furthermore, the accuracy levels of the linear beam theory are determined for thin beams under large deflections and small rotations as a function of moving mass weight and velocity in various boundary conditions.

Simulation of mould filling process using smoothed particle hydrodynamics

The implementation of high pressure die casting （HPDC） filling process modeling based on smoothed particle hydrodynamics （SPH） was discussed. A new treatment of inlet boundary was established by discriminating fluid particles from inlet particles. The roles of artificial viscosity and moving least squares method in the present model were compared in the handling pressure oscillation. The final model was substantiated by simulating filling process in HPDC in both two and three dimensions. The simulated results from SPH and finite difference method （FDM） were compared with the experiments. The results show the former is in a better agreement with experiments. It demonstrates the efficiency and precision of this SPH model in describing flow pattern in filling process.

Numerical Simulation of Sloshing with Large Deforming Free Surface by MPS-LES Method

Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

Efficient rendering of breaking waves using MPS method

This paper proposes an approach for rendering breaking waves out of large-scale of particle-based simulation. Moving particle semi-implicit (MPS) is used to solve the governing equation, and 2D simulation is expanded to 3D representation by giving motion variation using fractional Brownian motion (fBm). The waterbody surface is reconstructed from the outlines of 2D simulation. The splashing effect is computed according to the properties of the particles. Realistic features of the wave are ren-dered on GPU, including the reflective and refractive effect and the effect of splash. Experiments showed that the proposed method can simulate large scale breaking waves efficiently.

Numerical Simulation of the Solitary Wave Interacting with an Elastic Structure Using MPS-FEM Coupled Method

Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.

Implementation of the moving particle semi-implicit method on GPU

The Moving Particle Semi-implicit (MPS) method performs well in simulating violent free surface flow and hence becomes popular in the area of fluid flow simulation. However, the implementations of searching neighbouring particles and solving the large sparse matrix equations (Poisson-type equation) are very time-consuming. In order to utilize the tremendous power of parallel computation of Graphics Processing Units (GPU), this study has developed a GPU-based MPS model employing the Compute Unified Device Architecture (CUDA) on NVIDIA GTX 280. The efficient neighbourhood particle searching is done through an indirect method and the Poisson-type pressure equation is solved by the Bi-Conjugate Gradient (BiCG) method. Four different optimization levels for the present general parallel GPU-based MPS model are demonstrated. In addition, the elaborate optimization of GPU code is also discussed. A benchmark problem of dam-breaking flow is simulated using both codes of the present GPU-based MPS and the original CPU-based MPS. The comparisons between them show that the GPU-based MPS model outperforms 26 times the traditional CPU model.

Comparison of Improved Meshless Interpolation Schemes for SPH Method and Accuracy Analysis

In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been used for this purpose.Traditional particle interpolation (TPI) is simple and easy to do, but its low accuracy has become an obstacle to its wider application.This can be seen in the cases of particle disorder arrangements and derivative calculations.There are many different methods to improve accuracy, with the moving least square (MLS) method one of the most important meshless interpolation methods.Unfortunately, it requires complex matrix computing and so is quite time-consuming.The authors developed a simpler scheme, called higher-order particle interpolation (HPI).This scheme can get more accurate derivatives than the MLS method, and its function value and derivatives can be obtained simultaneously.Although this scheme was developed for the SPH method, it has been found useful for other meshless methods.

Two-phase micro-and macro-time scales in particle-laden turbulent channel flows

The micro- and macro-time scales in two-phase turbulent channel flows are investigated using the direct nu- merical simulation and the Lagrangian particle trajectory methods for the fluid- and the particle-phases, respectively. Lagrangian and Eulerian time scales of both phases are cal- culated using velocity correlation functions. Due to flow anisotropy, micro-time scales are not the same with the theo- retical estimations in large Reynolds number （isotropic） tur- bulence. Lagrangian macro-time scales of particle-phase and of fluid-phase seen by particles are both dependent on particle Stokes number. The fluid-phase Lagrangian inte- gral time scales increase with distance from the wall, longer than those time scales seen by particles. The Eulerian inte- gral macro-time scales increase in near-wall regions but de- crease in out-layer regions. The moving Eulerian time scales are also investigated and compared with Lagrangian integral time scales, and in good agreement with previous measure- ments and numerical predictions. For the fluid particles the micro Eulerian time scales are longer than the Lagrangian ones in the near wall regions, while away from the walls the micro Lagrangian time scales are longer. The Lagrangian integral time scales are longer than the Eulerian ones. The results are useful for further understanding two-phase flow physics and especially for constructing accurate prediction models of inertial particle dispersion.

Research on Numerical Wave Tank Based on the Improved Moving-Particle Semi-Implicit Method with Large Eddy Simulation

Moving-particle semi-implicit(MPS) method is a new mesh-free numerical method based on Lagrangian particle. In this paper, MPS method is applied to the study on numerical wave tank. For the purpose of simulating numerical wave, we combine the MPS method with large eddy simulation(LES) which can simulate the turbulence in the flow. The intense pressure fluctuation is a significant shortcoming in MPS method. So, we improve the original MPS method by using a new pressure Poisson equation to ease the pressure fluctuation. Divergencefree condition representing fluid incompressible is used to calculate pressure smoothly. Then, area-time average technique is used to deal with the calculation. With these improvements, the modified MPS-LES method is applied to the simulation of numerical wave. As a contrast, we also use the original MPS-LES method to simulate the wave in a numerical wave tank. The result shows that the new method is better than the original MPS-LES method.

浸入运动边界-格子Boltzmann方法4种固含率计算方法对比研究

为了达到流固耦合,格子Boltzmann方法(LBM)可采用浸入运动边界法(IMB)实现移动颗粒边界上的无滑移条件.该耦合方式(IMB-LBM)中固含率计算方法对流固耦合计算精度和效率有影响.对常用的固含率4种计算方法,即蒙特卡洛法(MCM)、单元分解法(UDM)、近似多边形法(APM)和闭合边界法(CBM),分别阐述其具体算法,对比了它们的计算精度和计算效率;最后通过圆盘颗粒非连续变形分析方法(DDDA)与IMB-LBM耦合模型下的一个多颗粒沉降流固耦合算例,对比分析了它们在流固耦合计算过程中的耗时.结果表明:1)CBM无误差,MCM和UDM在随机点数取1000,子单元数取100时误差稳定在1%以下,APM在颗粒直径大于格子长度10倍时,误差小于0.44%;2)MCM和UDM的计算精度及耗时分别与随机点数和子单元数相关,它们的计算耗时大于APM和CBM;3)计算效率上,APM>CBM>UDM>MCM,其中CBM计算耗时略微大于APM,APM和UDM计算耗时分别比MCM少2个和1个数量级.该结果可为IMB-LBM耦合模型中固含率计算方法优选提供借鉴.

油液环境中有载分接开关快速机构切换特性

利用基于移动粒子半隐式方法的Particleworks,分析有载分接开关快速机构在油液环境中的切换特性.介绍有载分接开关快速机构的工作原理,设计有载分接开关快速机构的机械结构,给出快速机构在油液环境下的动力学计算基础.建立该有载分接开关快速机构在油液环境中的动力学仿真模型,对照了该快速机构在有无变压器油液环境中的切换特性差异,得到驱动电机的驱动扭矩要求大于40 N·m.分析油液黏度和储能弹簧刚度对该快速机构切换特性的影响,得到快速机构在油液环境下切换过程中飞轮受到的油液阻力矩的变化规律,通过流场分析揭示了飞轮受到的阻力矩变化和波动的原因.

液滴撞击弹性悬臂梁行为的全拉格朗日数值方法

针对以往研究容易忽略弹性壁面变形对液滴撞击壁面的影响,以及传统网格方法在模拟自由表面大变形、弹性体耦合上易出现网格畸变等问题,提出了耦合移动粒子半隐式法(MPS)与离散单元法(DEM)的全拉格朗日液滴-弹性体相互作用数值方法。基于该方法,模拟了跳板效应、液滴反弹和铺展、板变形和振动等非线性现象,并通过与实验或理论解的对比,验证了方法的准确性。进一步探究了液滴撞击悬臂薄板的动态过程,分析了撞击速度和角度对液滴与薄板相互作用的影响。数值结果表明:悬臂梁末端最大偏移量随法向撞击速度呈线性增加,当液滴撞击角度θ≥90°时,接触时间随撞击角度的增加呈现先减小后增大然后再减小的趋势,当撞击角度θ<90°时,接触时间随撞击角度的增加呈现先增大后减小的趋势;液滴重心投影点的位置是影响液滴铺展系数与悬臂梁最大偏移量的主要因素,法向韦伯数相同的情况下,θ≥90°时的液滴铺展系数始终大于θ<90°的,而θ<90°时的悬臂梁最大偏移量始终大于θ≥90°的。该研究加深了对液滴撞击动力学的理解,可为涂层工艺和液滴传感技术等提供参考。

基于改进对流粒子域插值物质点法的隧道大变形分析

物质点法(MPM)在模拟大变形问题时具有很好的效果,然而传统的MPM在粒子穿越网格边界时存在单元穿越误差,导致精度降低。为克服传统MPM的单元穿越误差,基于对流粒子域插值物质点法(CPDI)理论框架,采用自适应正交改进插值移动最小二乘法(AOIIMLS),提出了改进CPDI方法。AOIIMLS通过构造加权正交基函数,并且忽略了新对角矩阵中的零元素或极小元素的贡献,以避免求解逆矩阵,增强了鲁棒性。改进CPDI采用速度梯度计算粒子域的速度场,粒子速度和粒子域角点速度用于重构背景网格速度函数。通过一维柱在自重作用下的压缩、砂柱坍塌和隧道坍塌离心机试验验证了改进CPDI方法的准确性和适用性,结果表明改进CPDI降低了单元穿越误差,得到了更高的精度。最后,采用改进CPDI方法模拟了青岛地铁4号线静沙区间地面塌陷全过程,验证了改进CPDI方法在岩土工程大变形领域的适用性及优势。