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.

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.

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.

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.

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.

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.

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.

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.

Numerical Simulation of Unidirectional Stratified Flow by Moving Particle Semi Implicit Method

Numerical simulation of stratified flow of two fluids between two infinite parallel plates using the Moving Particle Semi-implicit(MPS)method is presented.The developing process from entrance to fully development flow is captured.In the simulation,the computational domain is represented by various types of particles.Governing equations are described based on particles and their interactions.Grids are not necessary in any calculation steps of the simulation.The particle number density is implicitly required to be constant to satisfy incompressibility.The weight function is used to describe the interaction between different particles.The particle is considered to constitute the free interface if the particle number density is below a set point.Results for various combinations of density,viscosity,mass flow rates,and distance between the two parallel plates are presented.The proposed procedure is validated using the derived exact solution and the earlier numerical results from the Level-Set method.Furthermore,the evolution of the interface in the developing region is captured and compares well with the derived exact solutions in the developed region.

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.

Analysis of shock wave reflection from fixed and moving boundaries using a stabilized particle method

In the present paper, the efficiency of an enhanced formulation of the stabilized corrective smoothed particle method （CSPM） for simulation of shock wave propagation and reflection from fixed and moving solid boundaries in compressible fluids is investigated. The Lagrangian nature and its accuracy for imposing the boundary conditions are the two main reasons for adoption of CSPM. The governing equations are further modified for imposition of moving solid boundary conditions. In addition to the traditional artificial viscosity, which can remove numerically induced abnormal jumps in the field values, a velocity field smoothing technique is introduced as an efficient method for stabilizing the solution. The method has been implemented for one- and two-dimensional shock wave propagation and reflection from fixed and moving boundaries and the results have been compared with other available solutions. The method has also been adopted for simulation of shock wave propagation and reflection from infinite and finite solid boundaries.

浸入运动边界-格子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.分析油液黏度和储能弹簧刚度对该快速机构切换特性的影响,得到快速机构在油液环境下切换过程中飞轮受到的油液阻力矩的变化规律,通过流场分析揭示了飞轮受到的阻力矩变化和波动的原因.

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

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

考虑负荷特性与工况的公交充电站光储配置

随着光储充一体化充电站的发展,如何在现有公交充电站基础上合理确定光储容量成为了焦点。因此提出了考虑负荷特性与运行工况的光储公交充电站容量配置模型。首先考虑公交充电站充电负荷特性,通过移动边界法确定储能充放电状态的阈值,构建更加最优的光储公交充电站运行工况及运行条件。其次,以极大化光伏利用率和极小化系统投资、运行成本为目标,采用粒子群算法进行优化求解,确定出光储公交充电站最优的光伏以及储能容量。基于实际运营数据,对考虑负荷特性和运行工况的光储公交充电站和未考虑的光储公交充电站进行经济性对比。结果表明,在相同配置情况下考虑负荷特性和运行工况能够得到更加合理地光伏和储能容量。

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 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.

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.