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.

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.

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.

A spawning particle filter for defocused moving target detection in GNSS-based passive radar

Global Navigation Satellite System(GNSS)-based passive radar(GBPR)has been widely used in remote sensing applications.However,for moving target detection(MTD),the quadratic phase error(QPE)introduced by the non-cooperative target motion is usually difficult to be compensated,as the low power level of the GBPR echo signal renders the estimation of the Doppler rate less effective.Consequently,the moving target in GBPR image is usually defocused,which aggravates the difficulty of target detection even further.In this paper,a spawning particle filter(SPF)is proposed for defocused MTD.Firstly,the measurement model and the likelihood ratio function(LRF)of the defocused point-like target image are deduced.Then,a spawning particle set is generated for subsequent target detection,with reference to traditional particles in particle filter(PF)as their parent.After that,based on the PF estimator,the SPF algorithm and its sequential Monte Carlo(SMC)implementation are proposed with a novel amplitude estimation method to decrease the target state dimension.Finally,the effectiveness of the proposed SPF is demonstrated by numerical simulations and pre-liminary experimental results,showing that the target range and Doppler can be estimated accurately.

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.

RESEARCH ON STABILITY OF MOVING JET CONTAINING DENSE SUSPENDED SOLID PARTICLES

The spatial stability equation of moving jet containing dense suspended solid particles was derived out by means of the continuum phase-coupled model. The stability curves of moving jet far different downstream distances, Reynolds number of flow-field, particle properties and velocities of jetting device are got by the finite difference method based on the asymptotic method and the Eulerian conservative difference scheme. Founded on the analysis of the obtained stability curves it is found that the positive velocity of jetting device widens the unstable frequency range of flow-field hut the effect of the negative one is contrary. In addition, particles existing in the flow-field curb the instability of flow-field and the effect enhances with the decrease of Reynolds number of flow-field. These conclusions benefit learning the development of moving two-phase jet.

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.

Sedimentation motion of sand particles in moving water(Ⅰ):The resistance on a small sphere moving in non-uniform flow

In hydraulics,when we deal with the problem of sand particles moving relative to the surrounding water,Stokes'formula of resistance has usually been used to render the velocity of sedimentation of the particles.But such an approach has not been proved rigorously,and its accuracy must be carefully considered.In this paper,we discuss the problem of a sphere moving in a non-uniform flow field,on the basis of the fundamental theory of hydrodynamics.We introduce two assumptions:i)the diameter of the sphere is much smaller than the linear dimension of the flow field,and ii)the velocity of the sphere relative to the surrounding water is very small.Using these two assumptions,we solve the linearized Navier-Stokes equations and equations of continuity by the method of Laplace transform,and finally we obtain a formula for the resistance acting on a sphere moving in a non-uniform flow field.

Nearfield acoustic holography in a moving medium based on particle velocity input using nonsingular propagator

Nearfield acoustic holography in a moving medium is a technique which is typically suitable for sound sources identification in a flow.In the process of sound field reconstruction,sound pressure is usually used as the input,but it may contain considerable background noise due to the interactions between microphones and flow moving at a high velocity.To avoid this problem,particle velocity is an alternative input,which can be obtained by using laser Doppler velocimetry in a non-intrusive way.However,there is a singular problem in the conventional propagator relating the particle velocity to the pressure,and it could lead to significant errors or even false results.In view of this,in this paper,nonsingular propagators are deduced to realize accurate reconstruction in both cases that the hologram is parallel to and perpendicular to the flow direction.The advantages of the proposed method are analyzed,and simulations are conducted to verify the validation.The results show that the method can overcome the singular problem effectively,and the reconstruction errors are at a low level for different flow velocities,frequencies,and signal-to-noise ratios.

Effect of Ultrafine Particles on Flow Field and Transport Properties near the Interface Around a Moving Bubble

Laser Doppler Anemometer has been used to measure the flow field characteristics near the interface around a moving bubble in the presence of ultrafine particles. In order to model a moving bubble, the bubble was fixed into the counter flow liquid by a metal mesh. Experimental materials are air and water, and the particles are complex oxidate powder. Experiments were carried out under the operating conditions: the liquid flow velocity u 0 is 12.6 cm/s, the equivalent diameter d e is 0.6 cm, the mass concentration of particle is 0.2 0 0 ,the average particle diameter is about 10 nm and the density is 2 g/cm 3. The velocity profiles of both frontal and tail vortex areas were measured respectively. The experimental results show that the velocity fields are obviously changed in the existence of particles. In the frontal area of the bubble, both tangential and normal velocities decrease due to the presence of particles, but in tail vortex area, the tangential velocities increase remarkably, and normal velocities rise gradually from the center towards the fringe in the opposite tendency to that of no particles. The influences of flow field change in the presence of particles on gas liquid mass transfer are analyzed and discussed.

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.

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.

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.

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.

基于IPSO-LSTM的井下动目标位置预测实验研究

提升井下人员定位精度能够加强矿山安全监测,最大程度保障井下人员的生命安全。针对现有测距类算法受现场环境影响致使定位精度不足的问题,提出一种基于IPSO-LSTM的定位模型,应用于井下动目标的位置预测。采用LSTM构建指纹定位模型,通过UWB无线模块采集距离信息以构建距离-位置指纹关系数据库,利用数据库对PSO-LSTM模型进行训练,最后将训练好的模型进行目标轨迹预测。为比较不同改进策略对PSO的提升效果,对比了混沌映射随机初始化种群位置、非线性惯性权重递减、非对称优化学习因子和适应度函数优化4种改进策略,实验证明改进的PSO优化算法收敛速度快、鲁棒性好。为验证IPSO-LSTM的定位效果,以平均定位误差作为评价指标,将IPSO-LSTM模型与Chan算法、PSO-LSTM模型、LSTM神经网络、SSA-LSTM模型和GWO-LSTM进行对比,结果显示,IPSO-LSTM定位模型的平均定位误差为30 mm,相对传统Chan算法、LSTM、PSO-LSTM模型分别提升了76%、49%、24%。为降低局部误差偏大的现象,采用中值滤波对输入信息处理,进一步提升了定位精度。研究对进一步提高现有井下动目标定位系统的精度和稳定性具有重要意义和参考价值。

浸入运动边界-格子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耦合模型中固含率计算方法优选提供借鉴.

改进海洋捕食者算法的机器人路径规划研究

为解决海洋捕食者算法(MPA)收敛速度慢、收敛精度低、易陷入局部最优的问题,提出一种多策略改进海洋捕食者算法(IMPA)。引入Logistic混沌映射初始化种群,增加捕食者种群多样性;基于当前迭代次数t的自适应移动步长动态调整策略,增强算法逃离局部最优的能力;在IMPA迭代后期,加入中垂线算法(MA),基于中垂线策略的游离粒子位置更新方法,能够加快更新捕食者的位置,增强算法的寻优速度和寻优精度,避免算法陷入局部最优。最后通过改变IMPA阶段转换寻优过程,进一步平衡搜索过程,加强全局与局部适应性。选用6个基准测试函数对算法性能进行测试,测试结果显示:IMPA收敛速度更快,收敛精度更高;最后将改进算法应用于移动机器人路径规划,仿真结果表明:该算法规划的路径长度更短,搜索效率更高。

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

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

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

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

明渠湍流结构对粒子运动跟随性影响综述

近年来大量学者通过室内试验、现场观测以及数值模拟等方法对明渠湍流中粒子的运动规律进行了深入细致的研究,逐步揭示了明渠湍流中粒子运动复杂性的本质和规律,并将成果应用于河流泥沙输移、河流地貌演变以及生态环境保护等多个领域。鉴于明渠湍流具有明显的分区效应,其湍流流速分布、分区结构组成以及相干结构分布特性等对水中粒子的运动过程均会产生显著影响,对水中粒子在明渠湍流结构中的整体运动规律、湍流结构内区以及外区的运动过程进行了总结;重点阐述了近壁区粒子运动规律的反常现象和原因,以及湍流结构中粒径、密度和Stokes数这3个影响粒子运动过程的主要因素。研究成果可为泥沙运动以及生态保护等研究及应用提供参考。