Numerical investigations of the effect of oblique impact on particle deformation in cold spraying by the SPH method

In this study, a systematic examination of the oblique impacting of copper particles in cold spraying was conducted by using the smoothed particle hydrodynamics (SPH) method compared to the Lagrangian method. 3D models were employed owing to the asymmetric characteristic of the oblique impacting. It is found that in the oblique impact, the additional tangential component of particle velocity along the substrate surface could create a tensile force and decrease the total contact area and bonding strength between the particle and the substrate. The simulation results compare fairly well to the experiment results. Meanwhile, the asymmetric deformation can result in the focus of the shear friction on a small contact zone at one side, which may rise the interfacial temperature and thus facilitate the occurrence of the possible shear instability. Therefore, there probably exists an angle range, where the deposition efficiency may be promoted rather than the normal angle. Moreover, the particle deformation behavior simulated by the SPH method is well comparable to that simulated by the Lagrangian method and the experimental results, which indicates the applicability of the SPH method for simulating the impact process in cold spraying besides the previously used Arbitrary Lagrangian Eulerian (ALE) method.     

Group and phase velocities in the free-surface visco-potential flow: New kind of boundary layer induced instability

Water wave propagation can be attenuated by various physical mechanisms. One of the main sources of wave energy dissipation lies in boundary layers. The present work is entirely devoted to thorough analysis of the dispersion relation of the novel visco-potential formulation. Namely, in this study we relax all assumptions of the weak dependence of the wave frequency on time. As a result, we have to deal with complex integro-differential equations that describe transient behaviour of the phase and group velocities. Using numerical computations, we show several snapshots of these important quantities at different times as functions of the wave number. Good qualitative agreement with previous study [D. Dutykh, Eur. J. Mech. B/Fluids 28 (2009) 430] is obtained. Thus, we validate in some sense approximations made anteriorly. There is an unexpected conclusion of this study. According to our computations, the bottom boundary layer creates disintegrating modes in the group velocity. In the same time, the imaginary part of the phase velocity remains negative for all times. This result can be interpreted as a new kind of instability which is induced by the bottom boundary layer effect.     

改进的物理粘性SPH方法及其在溃坝问题中的应用

在低雷诺数物理粘性SPH方法基础上引入再生核粒子法进行密度重构,既避免了用人工粘性所导致的数值耗散问题,又提高了低雷诺数物理粘性SPH方法的数值稳定性;以溃坝问题为例,对比分析低雷诺数物理粘性SPH方法和本文方法的仿真结果表明,本文方法可有效消除数值不稳定,压强和速度分布更加光滑,粒子秩序更好,可应用于雷诺数较高或粘性不可忽略的流动问题.     

Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method

The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503–1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710–1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.     

高速碰撞SPH方法模拟中的初始光滑长度和粒子间距

采用光滑粒子流体动力学方法对高速碰撞问题作数值模拟,分析初始光滑长度和非一致粒子间距对计算结果的影响,提出修正光滑长度法,并引入XSPH速度纠错公式.数值结果表明,初始光滑长度越大,弹丸的刚性越小,h₀的合理取值范围应为d₀<h₀≤1.5d₀;粒子间距越小,材料的刚性越小.非一致粒子间距和均匀粒子分布的计算结果吻合的较好.     

Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach

The stability and accuracy of three methods which enforce either a divergence-free velocity field, density invariance, or their combination are tested here through the standard Taylor–Green and spin-down vortex problems. While various approaches to incompressible SPH (ISPH) have been proposed in the past decade, the present paper is restricted to the projection method for the pressure and velocity coupling. It is shown that the divergence-free ISPH method cannot maintain stability in certain situations although it is accurate before instability sets in. The density-invariant ISPH method is stable but inaccurate with random-noise like disturbances. The combined ISPH, combining advantages in divergence-free ISPH and density-invariant ISPH, can maintain accuracy and stability although at a higher computational cost. Redistribution of particles on a fixed uniform mesh is also shown to be effective but the attraction of a mesh-free method is lost. A new divergence-free ISPH approach is proposed here which maintains accuracy and stability while remaining mesh free without increasing computational cost by slightly shifting particles away from streamlines, although the necessary interpolation of hydrodynamic characteristics means the formulation ceases to be strictly conservative. This avoids the highly anisotropic particle spacing which eventually triggers instability. Importantly pressure fields are free from spurious oscillations, up to the highest Reynolds numbers tested.     

SPH二阶粒子近似光滑函数的充分条件及数值验证

为实现SPH二阶粒子近似,讨论光滑函数应满足的充分条件,给出对应的光滑函数形式,分析二阶粒子近似解的唯一性,并进行数值验证.结果表明：和标准SPH方法相比,在周期边界条件下,随着粒子数加密,SPH方法二阶粒子近似的L₁误差减小且L₁阶趋近2.0.     

Dissipative particle dynamics study of velocity autocorrelation function and self-diffusion coefficient in terms of interaction potential strength

This research focuses on numerically investigating the self-diffusion coefficient and velocity autocorrelation function (VACF) of a dissipative particle dynamics (DPD) fluid as a function of the conservative interaction strength. Analytic solutions to VACF and self-diffusion coefficients in DPD were obtained by many researchers in some restricted cases including ideal gases, without the account of conservative force. As departure from the ideal gas conditions are accentuated with increasing the relative proportion of conservative force, it is anticipated that the VACF should gradually deviate from its normally expected exponentially decay. This trend is confirmed through numerical simulations and an expression in terms of the conservative force parameter, density and temperature is proposed for the self-diffusion coefficient. As it concerned the VACF, the equivalent Langevin equation describing Brownian motion of particles with a harmonic potential is adapted to the problem and reveals an exponentially decaying oscillatory pattern influenced by the conservative force parameter, dissipative parameter and temperature. Although the proposed model for obtaining the self-diffusion coefficient with consideration of the conservative force could not be verified due to computational complexities, nonetheless the Arrhenius dependency of the self-diffusion coefficient to temperature and pressure permits to certify our model over a definite range of DPD parameters.     

波面斜率结构函数及其在自适应光学扩展信标波前探测中的应用

从结构函数的定义出发 ,推导出了斜率结构函数表达式 ,并且在模拟大气随机波前的仿真实验中得以验证。同时 ,将斜率结构函数应用于扩展信标波前探测中 ,可作为选择最佳基准子波面的依据。