A kernel gradient free (KGF) SPH method

The finite particle method (FPM) is a modified SPH method with high order accuracy while retaining the advantages of SPH in modeling problems with free surfaces, moving interfaces, and large deformations. In both SPH and FPM, kernel gradient is necessary in kernel and particle approximation of a field function and its derivatives. In this paper, a new FPM is presented, which only involves kernel function itself in kernel and particle approximation. The kernel gradient is not necessary in the whole computation, and this approach is thus referred to as a kernel gradient free (KGF) SPH method. This is helpful when a kernel function is not differentiable or the resultant kernel gradient is not sufficiently smooth, and thus it is more general in selecting a kernel function. Moreover, different from the original FPM with an asymmetric corrective matrix, in the new FPM, the resultant corrective matrix is symmetric, and this is advantageous in particle approximations. A series of numerical examples have been conducted to show the efficiencies of KGF‐SPH including one‐dimensional mathematical tests of polynomial functions with equal or variable smoothing length and two‐dimensional incompressible fluid flow of shear cavity. It is found that KGF‐SPH is comparable with FPM in accuracy and is flexible as SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

Normalized SPH without boundary deficiency and its application to transient solid mechanics problems

Smoothed particle hydrodynamics (SPH) method is a powerful tool for modeling solid mechanics problems, especially for large deformation problems. However, it suffers from boundary deficiency and difficulty of boundary condition treatment. In this work, a normalized SPH method is proposed to overcome these problems. The method is based on a newly developed normalized particle approximation. To derive this particle approximation, a normalized kernel approximation which is accurate for derivatives of linear functions everywhere in a problem domain is constructed, and all integral terms of the normalized kernel approximation including boundary terms are discretized by particle summations. The normalized particle approximation is free of matrix inversion, consequently attractive in computational stability and simplicity compared with other corrective particle approximations. Its approximation accuracy is demonstrated by calculating derivatives of test functions. Based on this particle approximation, the formulation of the normalized SPH method for transient solid mechanics problems is derived. Moreover, a direct method of treating traction boundary conditions is presented by making use of the boundary term of the normalized particle approximation. The accuracy and capability of the normalized SPH method are validated by the calculation of elastic wave propagation in solids and compared with commonly used SPH method.

     

固体介质中SPH方法的拉伸不稳定性问题研究进展

光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种基于核估计的无网格Lagrange数值方法.它用粒子方程离散流体动力学的连续方程, 既可以处理有限元难于处理的大变形和严重扭曲问题, 又可以处理有限差分法不易处理的自由边界和材料界面的问题, 在固体力学中的冲击、爆炸和裂纹模拟中具有广阔的发展前景.但是, 该算法的拉伸不稳定性(tensile instability)问题是它在固体力学领域中应用的最大障碍.对SPH稳定性分析表明, 算法不稳定性的条件仅与应力状态和核函数的2阶导数有关.目前, 应力点法(stress points)、Lagrange核函数法、人工应力法(artificialstress)、修正光滑粒子法(corrective smoothed particle method, CSPM)和守恒光滑法(conservativesmoothing)以及其他一些方法成功地改善了SPH的拉伸不稳定性, 但是每一种方法都不能彻底解决SPH的拉伸不稳定性问题.本文介绍了SPH法的方程和Von Neumann稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.      

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

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

Adaptive smoothed particle hydrodynamics for high strain hydrodynamics with material strength

This paper presents the implementation of an adaptive smoothed particle hydrodynamics (ASPH) method for high strain Lagrangian hydrodynamics with material strength. In ASPH, the isotropic kernel in the standard SPH is replaced with an anisotropic kernel whose axes evolve automatically to follow the mean particle spacing as it varies in time, space, and direction around each particle. Except for the features inherited from the standard SPH, ASPH can capture dimension-dependent features such as anisotropic deformations with a more generalized elliptical or ellipsoidal influence domain. Two numerical examples, the impact of a plate against a rigid surface and the penetration of a cylinder through a plate, are investigated using both SPH and ASPH. The comparative studies show that ASPH has better accuracy than the standard SPH when being used for high strain hydrodynamic problems with inherent anisotropic deformations. PACS 46.15.-x, 83.10.Rs, 83.50.-v     

RESEARCH DEVELOPMENTS OF SMOOTHED PARTICLE HYDRODYNAMICS METHOD AND ITS COUPLING WITH FINITE ELEMENT METHOD

拉格朗日型的有限元法和光滑粒子法在模拟材料大变形问题时各存优缺点, 而有限元与光滑粒子耦合算法实现了在小变形区域采用有限元法计算, 在局部的大变形区域采用光滑粒子法计算, 有效地综合了有限元法计算效率高和光滑粒子法能够自然地模拟材料大变形问题的特点.重点论述了有限元法、光滑粒子法以及有限元与光滑粒子耦合算法的研究现状及应用进展, 并讨论了各方法中需要进一步解决的问题.最后通过算例对3种方法的计算精度和计算效率进行了分析, 供研究人员参考.     

A compatibly differenced total energy conserving form of SPH

We describe a modified form of smoothed particle hydrodynamics (SPH) in which the specific thermal energy equation is based on a compatibly differenced formalism, guaranteeing exact conservation of the total energy. We compare the errors and convergence rates of the standard and compatible SPH formalisms on a variety of shock test problems with analytic answers. We find that the new compatible formalism reliably achieves the expected first‐order convergence for these analytic shock tests and, in all cases, improves the accuracy of the numerical solution over the standard formalism. We also examine the performance of our new formalism on a more complicated applied problem: the diversion of an asteroid by a kinetic impactor. We find the compatible discretization demonstrates measurable improvement in the convergence of properties such as the deflection velocity in this kind of applied problem as well. Copyright © 2014 John Wiley & Sons, Ltd.     

A new version of smoothed point method to simulate linear elastic wave propagation

By using the kernel function of the smoothed particle hydrodynamics (SPH) and modification of statistical volumes of the boundary points and their kernel functions, a new version of smoothed point method is established for simulating elastic waves in solid. With the simplicity of SPH kept, the method is easy to handle stress boundary conditions, especially for the transmitting boundary condition. A result improving by de-convolution is also proposed to achieve high accuracy under a relatively large smooth length. A numerical example is given and compared favorably with the analytical solution.     

SPH 方法在聚能装药射流三维数值模拟中的应用

采用有限元软件LS-DYNA中的SPH 算法实现聚能装药射流形成过程的三维数值模拟。将SPH方法和FEM 方法与理论计算结果相比较。研究结果表明:SPH 数值模拟计算过程十分稳定,避免了有限元数值模拟过程中的网格扭曲、缠绕和物质穿透等问题,而且计算精度比有限元方法更高;采用SPH 方法计算的射流速度、射流长度与有限元计算结果和理论计算结果的误差均小于5%,因此该方法对于模拟诸如聚能装药等多介质、大变形问题十分有效;初始相邻粒子数N 的最佳取值范围为600~1800,在该范围内计算结果的误差小于3%,且计算效率高。     

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

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

Meshfree particle simulation of the detonation process for high explosives in shaped charge unlined cavity configurations

The numerical simulation of the detonation of a high explosive (HE) is generally not an easy task for traditional grid based methods. Smoothed particle hydrodynamics (SPH) method, as a meshfree, Lagrangian and particle method, provides a very attractive approach in dealing with large deformations and large inhomogeneities in the extremely transient high explosive detonation and later expansion process. This paper presents the application of SPH to simulate and analyze the detonation process of high explosive in shaped charge. A three-dimensional SPH code is developed and applied to simulate the shaped charge detonation process in different scenarios. It is observed that for high explosive in a shaped charge, the detonation produced gaseous products experience strong convergence that forms an extremely high-pressure gas jet. Factors such as different charge cavity shapes and different detonation models lead to quite different behavior of the gas jet convergence and later divergence. Further analyses reveal that a critical value for the charge head length exists. Beyond this critical value, increasing the charge head length will not result in improvement on the gas jet convergence performance.Received: 11 March 2002, Accepted: 9 December 2002, Published online: 28 April 2003     

SPH-FEM接触算法在冲击动力学数值计算中的应用

为了充分发挥光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)易于处理大变形以及有限元(Finite Element Method,FEM)计算精度和效率高的优势,论文基于无网格粒子接触算法,在有限元节点处设置背景粒子,通过接触力的方式计算SPH粒子和有限单元之间的接触问题.使用SPH-FEM接触算法分别对球头钢弹斜冲击钢板和平头钢弹正冲击钢板的情况进行了三维数值计算,SPH采用完全变光滑长度算法,FEM采用矩阵向量积的EBE(Element-By-Element)算法.将SPH-FEM接触算法的计算结果与LS-DYNA的计算结果以及实验结果进行了对比验证.     

多相流动的光滑粒子流体动力学方法研究综述

光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)具有粒子方法的无网格和全拉格朗日特征, 适用于具有界面大变形、不连续性和多物理场的多相流的高精度模拟. SPH方法模拟多相流已有大量报道, 具体的实现方式也大不相同. 本文首先阐述了采用SPH方法模拟流体的基本控制方程, 以及求解过程中需要考虑的流体压力求解、表面张力、固体边界等问题. 整理和总结了基于SPH方法进行多相流模拟的主要实现方式: (1)双流体模型的拉格朗日求解器: 两相离散为两组独立SPH粒子, 并用显式相间作用耦合两相; (2)多相SPH方法: SPH方法对多相流模拟的自然延伸, 相间作用由SPH参数隐式描述; (3) SPH与其他离散方法的耦合: 差异较大的两相各自采用不同离散方法, 发挥不同拉格朗日方法的优点; (4) SPH和基于网格方法的耦合: 网格方法处理简单的单相流动主体, 获得精度和效率间的平衡. 另外, 还在模拟参数物理化等方面论述了与SPH方法模拟多相流相关的一些改进和修正方法, 并在最后讨论和建议了提高多相流SPH模拟效率和精度的措施.      

12.7 mm弹侵彻不同强度钢靶的数值模拟

针对12.7 mm弹侵彻不同强度钢靶时可能出现子弹保持完整或发生破碎的情况,过去的数值模拟仅限于模拟单一模式的子弹侵彻行为。为了克服这种数值模拟的局限性,开展了模型算法、网格尺寸对模拟结果影响的研究,并将模拟结果与实验结果进行了对比,提出了一种能够用于模拟子弹保持完整或破碎的弹靶模型。研究结果表明,为模拟子弹保持完整状态,子弹和靶板应分别采用基于Lagrange算法的有限元法和光滑粒子算法,而且子弹网格尺寸和靶板粒子间距之比应至少保持在5.3左右,否则弹头会产生与实验结果不符合的异常变形。但是,在模拟子弹发生破碎侵蚀时,该比例的网格/粒子尺寸比会引起计算中止。为了克服该问题,进一步建立了一种弹体表面采用大尺寸网格、内部采用细化小尺寸网格的有限元/光滑粒子法耦合弹靶模型。计算结果表明,改进的弹靶模型可模拟子弹保持完整或者发生破碎的情况。

     

Simulation of container filling process with two inlets by improved smoothed particle hydrodynamics (SPH) method

In this article, an improved smoothed particle hydrodynamics (SPH) method is proposed to simulate the filling process with two inlets. Improvements are achieved by deriving a corrected kernel gradient of SPH and a density re-initialisation. In addition, a new treatment of solid wall boundaries is presented. Thus, the improved SPH method has higher accuracy and better stability, and conserves both linear and angular momentums. The validity of the new boundary treatment is shown by simulating the spin-down problem. The bench tests are also presented to demonstrate the performance of the improved SPH method. Then the filling process with a single inlet is simulated to show the ability to capture complex-free surface of the proposed method. Finally, the filling process with two inlets is numerically investigated. The numerical results show that the filling patterns are affected significantly by Reynolds number, aspect ratio of the container and the location of the inlets.     

无网格稳定配点法及其在弹性力学中的应用

伽辽金型无网格法具有精度高、稳定性好的优点,但是实现高阶准确积分过程复杂,计算效率低.配点型无网格法的计算效率高,但是其在求解复杂问题时往往会出现精度和稳定性较差的结果.本文介绍一种新的无网格法-无网格稳定配点法,采用重构核近似作为近似函数,在规则子域内非常容易实现高阶准确积分,既保留了配点型无网格法效率高的特点,又具...     

A monotone finite volume scheme for advection–diffusion equations on distorted meshes

A new monotone finite volume method with second‐order accuracy is presented for the steady‐state advection–diffusion equation. The method uses a nonlinear approximation for both diffusive and advective fluxes that guarantee the positivity of the numerical solution. The approximation of the diffusive flux is based on nonlinear two‐point approximation, and the approximation of the advective flux is based on the second‐order upwind method with proper slope limiter. The second‐order convergence rate for concentration and the monotonicity of the nonlinear finite volume method are verified with numerical experiments. Copyright © 2011 John Wiley & Sons, Ltd.     

A robust shell element in meshfree SPH method

With the incorporation of total Lagrangian smoothed particle hydrodynamics(SPH) method equation and moving least square(MLS) function,the traditional SPH method is improved regarding the stability and consistency.Based on Mindlin-Ressiner plate theory,the SPH method simulating dynamic behavior via one layer of particles is applied to plate’s mid-plane,i.e.,a SPH shell model is constructed.Finally,through comparative analyses on the dynamic response of square,stiffened shells and cylindrical shells under various strong impact loads with common finite element software,the feasibility,validity and numerical accuracy of the SPH shell method are verified.Consequently,further researches on SPH shell may well pave the way towards solving problems involving dynamic plastic damage,tearing or even crushing.