A modified SPH method to model the coalescence of colliding non-Newtonian liquid droplets

This paper develops a modified smoothed particle hydrodynamics (SPH) method to model the coalescence of colliding non-Newtonian liquid droplets. In the present SPH, a van der Waals (vdW) equation of state is particularly used to represent the gas-to-liquid phase transition similar to that of a real fluid. To remove the unphysical behavior of the particle clustering, also known as tensile instability, an optimized particle shifting technique is implemented in the simulations. To validate the numerical method, the formation of a Newtonian vdW droplet is first tested, and it clearly demonstrates that the tensile instability can be effectively removed. The method is then extended to simulate the head-on binary collision of vdW liquid droplets. Both Newtonian and non-Newtonian fluid flows are considered. The effect of Reynolds number on the coalescence process of droplets is analyzed. It is observed that the time up to the completion of the first oscillation period does not always increase as the Reynolds number increases. Results for the off-center binary collision of non-Newtonian vdW liquid droplets are lastly presented. All the results enrich the simulations of the droplet dynamics and deepen understandings of flow physics. Also, the present SPH is able to model the coalescence of colliding non-Newtonian liquid droplets without tensile instability.     

基于SPH/FEM方法的浅水冲击喷溅及其抑制结构研究

物体浅水冲击会造成结构响应和自由液面的剧烈喷溅,以往的研究多集中在物体受水的冲击响应方面,而对浅水层受物体冲击引起的剧烈喷溅问题及其喷溅抑制结构的相关研究很少。本文基于SPH/FEM耦合方法对浅水冲击的液体喷溅特性开展研究。首先,计算了刚性长方体结构的冲击喷溅,与文献中准二维溅水试验结果进行了比较,将浅水喷溅分为初始喷溅、二次喷溅和底部射流三部分。分别研究了各部分喷溅的形成机理,指出二次喷溅由结构对水层的挤压和底部射流的能量传递共同引起,并分析了水深对二次喷溅的影响规律。其次,对喷溅抑制结构的机理做了数值研究,发现拱形翻边构型可以抑制初始喷溅,但不能有效降低二次喷溅。最后,提出了一种带凹槽和阻水边条的新型翻边构型,对二次喷溅具有更好的抑制效果。     

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

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

Sloshing and energy dissipation in an egg: SPH simulations and experiments

Tall structures, such as towers and bridges, can oscillate at excessive magnitudes when subjected to wind and earthquake loads. Liquid sloshing absorbers can be used to suppress these excessive oscillations by tuning the frequency of the sloshing to the critical frequency of the structure. Sloshing absorbers are simple structures consisting of a partially full container of liquid with a free surface. Tuning ensures that significant amounts of harmful energy can be extracted from the structure to the sloshing liquid. However, there needs to be a rapid means of dissipating this energy to avoid its returning back to the structure (then back to the liquid periodically).A hen׳s egg seems to have evolved to efficiently dissipate energy to protect its embryo using sloshing of its liquid content. Hence, the potential to implement the egg׳s unique properties as a sloshing absorber for structural control, is the main focus of this study. Numerical simulations, using Smoothed Particle Hydrodynamics (SPH), and experimental comparisons are presented in this paper. One objective is to demonstrate the ability of SPH to simulate complex free surface behaviour in three dimensions. Such a tool is then useful to identify different dissipation modes. Effects of fill volume and viscosity on the rate of dissipation, are also investigated.     

Simulating bubble dynamics in a buoyant system

Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two-dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three-dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.     

An improved particle shifting technique for incompressible smoothed particle hydrodynamics methods

The smoothed particle hydrodynamics (SPH) method is one of the powerful Lagrangian tools for modeling free surface flows. However, it suffers from particle disorder, which leads to interpolation and numerical errors. To overcome this problem, several techniques have been introduced until now, among which the particle shifting technique (PST) based on Fick's law is an efficient one. The current form of this method needs tuning parameters to fulfill numerical stability criteria. In this study, to eliminate calibration factors, a new shifting coefficient is derived theoretically based on particle positions before and after shifting, regardless of other parameters such as velocity, pressure, time step intervals, etc. The only required input is particle positions, and the main concern is conserving particle densities in their updated positions. In addition to the proposed PST, a new distribution index (DI) is introduced for measuring the spatial uniformity of particles. Furthering the research, some novel treatments are also studied to improve particle movements near free surface boundary. The proposed idea is only assessed for ISPH method in this study, and its performance in other SPH schemes needs more investigations. Following this innovative method, it is validated by modeling different cases including dam break flow, paddle movement, and elliptical water drop. In all cases, particle arrangements have been improved by means of the modified shifting method. In that sense, good agreements between simulation results with experimental data, analytical solutions, and other numerical methods approve the ability of the developed method in simulating free surface flows.     

弹塑性波计算中的光滑粒子法

采用改进的光滑粒子法 ,对脉冲应力载荷下板中弹塑性一维应变波的传播进行了数值计算 ,比较了人工粘性法和通量修正法处理冲击波间断面的效果。结果表明 ,改进的光滑粒子法在应力波数值计算中有良好的精度。     

NUMERICAL SIMULATION FOR FORMED PROJECTILE OF DEPLETED URANIUM ALLOY

IntroductionThecharactersoftheprocessesofformedprojectilebydetonationarelargestrainandhighstrainrate .InanalysisoftheprocesseswithLagrangianfinite_elementorfinite_differenceapproachestherearemajorproblemsthatsomemeshesalwaysbecomeseverelydistortiontocauselargecomputationalerrors .Toovercomethesedifficulties,themeshrezoningtechniquewasused[1] butrezoningnumberwaslimited .Smoothedparticlehydrodynamics (SPH)isuniqueincomputationalcontinuumdynamicsinthatitusesnospatialmesh .ItisapureLagrangianp…     

海冰动力学数值模拟中改进的PIC方法

为了准确地模拟海冰的动力过程,需要建立精确有效的数值方法。本文结合质点网格法(PIC)和光滑质点流体动力学方法(SPH)发展了一种改进的PIC方法。该方法在欧拉坐标下对海冰动量方程进行差分计算,在拉格朗日坐标下进行海冰质点位移、厚度和密集度计算,并采用Gauss函数进行欧拉网格点与拉格朗日质点间海冰参数的交互插值。采用改进的PIC方法对规则区域内的海冰堆积过程进行了数值试验,对渤海海冰的动力过程进行了72小时数值模拟。计算结果均表明改进的PIC方法具有计算量小,计算结果平稳精确的优点,可很好地适用于海冰动力作用过程的数值模拟。