An improved KGF‐SPH with a novel discrete scheme of Laplacian operator for viscous incompressible fluid flows

The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

Coupling finite difference method with finite particle method for modeling viscous incompressible flows

A hybrid approach to couple finite difference method (FDM) with finite particle method (FPM) (ie, FDM-FPM) is developed to simulate viscous incompressible flows. FDM is a grid-based method that is convenient for implementing multiple or adaptive resolutions and is computationally efficient. FPM is an improved smoothed particle hydrodynamics (SPH), which is widely used in modeling fluid flows with free surfaces and complex boundaries. The proposed FDM-FPM leverages their advantages and is appealing in modeling viscous incompressible flows to balance accuracy and efficiency. In order to exchange the interface information between FDM and FPM for achieving consistency, stability, and convergence, a transition region is created in the particle region to maintain the stability of the interface between two methods. The mass flux algorithm is defined to control the particle creation and deletion. The mass is updated by N-S equations instead of the interpolation. In order to allow information exchange, an overlapping zone is defined near the interface. The information of overlapping zone is obtained by an FPM-type interpolation. Taylor-Green vortices and lid-driven shear cavity flows are simulated to test the accuracy and the conservation of the FDM-FPM hybrid approach. The standing waves and flows around NACA airfoils are further simulated to test the ability to deal with free surfaces and complex boundaries. The results show that FDM-FPM retains not only the high efficiency of FDM with multiple resolutions but also the ability of FPM in modeling free surfaces and complex boundaries.     

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.     

SPH方法Delaunay三角刨分与自由液面重构

光滑粒子法(SPH)作为一种拉格朗日型无网格方法,兼具欧拉网格方法和拉格朗日网格方法的优势,已经成功应用于科学和工程的众多领域。SPH方法后处理一般基于无规则分布的粒子,不如网格类方法后处理简便、直接。另外,SPH方法模拟自由液面流动等问题时,通过粒子位置难以重构自由液面的准确位置。发展一种基于Delaunay三角刨分的SPH后处理方法,即先基于SPH粒子位置利用Delaunay三角刨分建立三角网格,然后将粒子信息转化成网格单元/节点信息,从而可以在三角网格上进行后处理,实现基于网格方法的后处理功能,并可以在三角网格上直接提取或重构自由液面。将本文的方法应用到液滴碰撞和溃坝流SPH模拟结果的后处理中,得到了非常好的结果,表明本文的方法有效可靠。     

A corrected symmetric SPH method to simulate viscoelastic free surface flows based on the PTT model

In this work, a corrected symmetric and periodic density reinitialized SPH (CSPDR‐SPH) method is proposed and extended to simulate the viscoelastic free surface flows based on the Phan–Thien–Tanner model. The improvements mainly lie in deriving a corrected symmetric kernel gradient, and combining it with a periodic density reinitialization procedure. In addition, a simple artificial viscosity and a simple artificial stress form are adopted. Thus, the CSPDR‐SPH method has higher accuracy and better stability than the SPH method, and conserves both linear and angular momentums. The consistency and convergence of the CSPDR‐SPH method are justified by approximating a function in one and two dimensions. The merits of CSPDR‐SPH method are demonstrated by several benchmarks. The simple flow in a two‐dimensional channel is investigated to show the capability of the CSPDR‐SPH method to simulate the viscoelastic free surface flow. Then the CSPDR‐SPH method is extended to simulate the impacting drop problem. Numerical results show that the CSPDR‐SPH method can precisely capture the viscoelastic free surface. The Reynolds number, Weissenberg number and elongation parameter have remarkable effect on the flows. Copyright © 2012 John Wiley & Sons, Ltd.     

一种SPH应力修正算法及自由表面流中的应用

提出了一种SPH应力修正算法,即模型中的拉应力和压应力分别采用不同的插值核函数和状态方程来处理,改善应力稳定性问题。介绍了一种改进的Quintic核函数,用于改善模型中压应力的稳定性。通过增加钟型核函数的光滑长度,改善模型中拉应力的稳定性。采用该应力修正算法模拟了无重力条件下方形液滴的震荡变形过程,对比分析了不同算法的模拟结果。此外,为进一步验证算法的适用性,模拟了溃坝算例。研究表明,改进的Quintic型核函数明显改善了粒子聚集现象,该SPH应力修正方法可以使液滴具有更均匀的粒子分布以及更光滑的自由表面,有效改善了SPH方法中的压应力不稳定作用以及自由表面流的模拟精度。     

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.     

SPH方法中常数一致性核函数的建立及公式化

施加边界条件的不足是SPH方法的一个棘手问题,因为在结构边界外没有颗粒的存在,使得在边界处核函数的单位特性不能得到满足。施加"伪"颗粒是目前通用的一种方法,但是对于不规则结构和复杂几何边界,确定这些"伪"颗粒非常困难。本文讨论通过使用满足常数一致性的核函数来改善边界的不足。文章首先通过三种方法推导了满足常数一致性条件的核函数及其函数梯度的表达式,发现了两个不同分母式的表达,分析了满足常数一致性的修正核函数的数学特性。开展了二维和三维的算例比较,结果发现使用修正的核函数对边界条件有明显改善,对计算精度和稳定性也有显著提高。     

一种改进的可阻止SPH数值断裂的方法

首先通过人为设定数值断裂,从数值上分析数值断裂对计算精度的影响。然后引进一种可适用于大变形的并保证系统质量守恒、动量守恒和能量基本守恒的简单加点技术,用以阻止数值断裂。并用改进的smoothed particle hydrodynamics(SPH)方法对钢弹撞击飞机蒙皮进行数值模拟以验证其阻止数值断裂的有效性。在130~250 m/s的碰撞速度下,用改进的SPH方法得到的变形区直径、剩余速度以及弹坑深度与试验结果吻合较好。     

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.     

Numerical simulation of soil–water interaction using smoothed particle hydrodynamics (SPH) method

An application of smoothed particle hydrodynamics (SPH) to simulation of soil–water interaction is presented. In this calculation, water is modeled as a viscous fluid with week compressibility and soil is modeled as an elastic–perfectly plastic material. The Mohr–Coulomb failure criterion is applied to describe the stress states of soil in the plastic flow regime. Dry soil is modeled by one-phase flow while saturated soil is modeled by separate water and soil phases. Interaction between soil and water is taken into account by means of pore water pressure and seepage force. Simulation tests of soil excavation by a water jet are calculated as a challenging example to verify the broad applicability of the SPH method. The excavations are carried out in two different soil models, one is dry soil and the other is fully saturated soil. Numerical results obtained in this paper have shown that the gross discontinuities of soil failure can be simulated without any difficulties. This supports the feasibility and attractiveness of this a new approach in geomechanics applications. Advantages of the method are robustness, conceptual simplicity and relative ease of incorporating new physics.     

应力波数值计算中的SPH方法

对一维波动方程的SPH(smoothed particle hydrodynamics)格式和有限差分格式进行比较,并采用SPH法模拟了一维应力/应变波, 获得1个可衡量SPH法模拟应力波准确性的重要指标。结果表明,SPH法模拟应力波传播中采用的光滑长度必须不小于粒子间距;采用B-样条核函数和高斯型核函数能够获得良好的应力波图像,而二次型核函数不能,因此二次型核函数不适用于冲击动力学的数值计算。     

高速碰撞数值计算中的SPH自适应粒子分布法

针对传统光滑粒子法在计算高速碰撞问题时会出现近邻粒子逸出核函数影响域而产生数值破坏这一缺陷,提出了一种根据粒子间距变化自动添加、合并粒子的SPH自适应粒子分布算法。采用该方法对Taylor碰撞和超高速碰撞问题进行了数值模拟,结果表明,该方法可以有效地消除计算中出现的数值破坏,提高计算精度。     

光滑粒子法中的一种新的核函数

分析了传统的核函数产生压缩失稳现象的原因,提出了消除这种压缩失稳现象的一种新的核函数。采用改进的光滑粒子法,对几种常用的核函数进行了一维应变波的对比计算。结果表明:所提出新的核函数在应力波计算中不但保证了计算精度,还能有效地消除压缩失稳。     

避免核函数导数的二阶导数核近似方法的修正

采用核近似光滑粒子流体力学SPH（Smooth Particle Hydrodynamics）方法计算一元函数二阶导数和多元函数二阶偏导数,对避免核函数导数算法进行了改进,提出了对核函数光滑长度处处具有o（h）精度的算法,并分别推导出一元和多元函数的修正公式。用不同的粒子间距和不同的光滑长度进行了二阶导数的计算和热传导问题的模拟,进行了修正方法与原方法的误差分析。结果表明,本文提出的修正措施在提高精度、减少误差及加快收敛速度等方面起到很大的作用。     

Investigation on charge parameters of underwater contact explosion based on axisymmetric SPH method

This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.     

SPH computation of plunging waves using a 2‐D sub‐particle scale (SPS) turbulence model

The paper presents a 2‐D large eddy simulation (LES) modelling approach to investigate the properties of the plunging waves. The numerical model is based on the smoothed particle hydrodynamics (SPH) method. SPH is a mesh‐free Lagrangian particle approach which is capable of tracking the free surfaces of large deformation in an easy and accurate way. The Smagorinsky model is used as the turbulence model due to its simplicity and effectiveness. The proposed 2‐D SPH–LES model is applied to a cnoidal wave breaking and plunging over a mild slope. The computations are in good agreement with the documented data. Especially the computed turbulence quantities under the breaking waves agree better with the experiments as compared with the numerical results obtained by using the k–ε model. The sensitivity analyses of the SPH–LES computations indicate that both the turbulence model and the spatial resolution play an important role in the model predictions and the contributions from the sub‐particle scale (SPS) turbulence decrease with the particle size refinement. Copyright © 2006 John Wiley & Sons, Ltd.     

基于核重构思想的配点型无网格方法的研究--一维问题

无网格方法按其离散原理可分为Galerkin型、配点型等。其中Galerkin型无网格方法的实施需要背景网格,不属于真正的无网格法;配点型无网格方法的实施不需要背景网格,是真正的无网格法。本文首先介绍了重构核点法的基本原理,然后基于核重构思想,与配点法相结合,以一维问题为例,研究了配点型无网格方法,对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨。并结合若干典型算例,检验了其计算精度与收敛姓。     

h-ADAPTIVITY ANALYSIS BASED ON MULTIPLE SCALE REPRODUCING KERNEL PARTICLE METHOD

IntroductionAs one of the meshfree methods, multiple scale reproducing kernel particle method(RKPM)[1,2]bears, besides the common features of all meshfree methods, a uniqueinherent feature of multiresolution analysis by which the structural response can b…     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.