基于SPH方法的低韦伯数下三维液滴碰撞聚合与反弹数值模拟研究

光滑粒子流体动力学方法（SPH方法）作为纯拉格朗日粒子方法,可以有效避免网格法在模拟大变形过程中带来的网格扭曲等缺陷,适合模拟含大变形的液滴碰撞聚合与反弹过程。该文基于Ott和Schnetter提出的修正SPH方法,利用有限差分与SPH一阶导数相结合的方法处理粘性项中的二阶导数问题,进行Couette流算例验证,数值解...     

Smoothed particle hydrodynamics simulation of shear-induced powder migration in injection moulding

We present the application of the smoothed particle hydrodynamics (SPH) discretization scheme to Phillips' model for shear-induced particle migration in concentrated suspensions. This model provides an evolution equation for the scalar mean volume fraction of idealized spherical solid particles of equal diameter which is discretized by the SPH formalism. In order to obtain a discrete evolution equation with exact conservation properties we treat in fact the occupied volume of the solid particles as the degree of freedom for the fluid particles. We present simulation results in two- and three-dimensional channel flow. The two-dimensional results serve as a verification by a comparison to analytic solutions. The three-dimensional results are used for a comparison with experimental measurements obtained from computer tomography of injection moulded ceramic microparts. We observe the best agreement of measurements with snapshots of the transient simulation for a ratio D(c)/D(η)=0.1 of the two model parameters.     

修正的SPH方法模拟强制旋转动边界问题研究

采用移动最小二乘(movingleastsquares,MLS)修正的光滑粒子流体动力学(smoothedparticlehydrodynamics,SPH)算法模拟了一个强制旋转动边界问题模型。提出了一种施加强制旋转粒子动边界方案。阐述了修正方法的原理并给出了具体的修正操作方法。同时,还构建了静边界模型,并分别使用修正的SPH以及商业软件FLUENT计算。对计算结果的比较分析表明:该修正SPH方法,能够消除压力振荡;该方法及动边界处理方案可以有效计算该强制动边界问题,为进一步计算更加复杂模型奠定了理论基础。同时,该方法对于不同领域中的动边界问题也具有一定参考价值。     

Prediction of powder particle size during centrifugal atomisation using a rotating disk

The centrifugal atomisation of metallic melts using a rotating disk is an important process for powder production and spray deposition. The theoretical prediction of powder particle size is desirable for the design of atomisers. In this paper, wave theory was applied to analyse the disintegration of metallic melts in the film disintegration regime during centrifugal atomisation using a rotating disk. A mathematical model was proposed to predict the spray parameters. The governing equation for the fastest-growing wave number was developed and solved numerically. The effect of the variation in film thickness during film extension was taken into account. Film length and powder particle size were calculated and compared with available experimental data in the literature, and a good agreement was achieved. The influence of the break-up parameter was studied, and it is shown that the break-up parameter is not sensitive to the predicted powder particle size. Both simulated results and experimental data showed that fine powders can be produced by increasing disk speed.     

Development of computationally efficient augmented Lagrangian SPH for incompressible flows and its quantitative comparison with WCSPH simulating flow past a circular cylinder

In Lagrangian particle-based methods such as smoothed particle hydrodynamics (SPH), computing totally divergence-free velocity field in a flow domain with the smallest error possible is the most critical issue, which might be achieved through solving pressure Poisson equation implicitly with higher particle resolutions. However, implicit solutions are computationally expensive and may be particularly challenging in the solution of multiphase flows with highly nonlinear deformations as well as fluid-structure interaction problems. Augmented Lagrangian SPH (ALSPH) method is a new alternative algorithm as a prevalent pressure solver where the divergence-free velocity field is achieved by iterative calculation of velocity and pressure fields. This study investigates the performance of the ALSPH technique by solving a challenging flow problem such as two-dimensional flow around a cylinder within the Reynolds number range of 50 to 500 in terms of improved robustness, accuracy, and computational efficiency. The same flow conditions are also simulated using the conventional weakly compressible SPH (WCSPH) method. The results of ALSPH and WCSPH solutions are not only compared in terms of numerical validation/ verification studies, but also rigorous investigations are performed for all related physical flow characteristics, namely, hydrodynamic coefficients, frequency domain analyses, and velocity divergence fields.     

Prediction of powder particle size during centrifugal atomisation using a rotating disk

The centrifugal atomisation of metallic melts using a rotating disk is an important process for powder production and spray deposition. The theoretical prediction of powder particle size is desirable for the design of atomisers. In this paper, wave theory was applied to analyse the disintegration of metallic melts in the film disintegration regime during centrifugal atomisation using a rotating disk. A mathematical model was proposed to predict the spray parameters. The governing equation for the fastest-growing wave number was developed and solved numerically. The effect of the variation in film thickness during film extension was taken into account. Film length and powder particle size were calculated and compared with available experimental data in the literature, and a good agreement was achieved. The influence of the break-up parameter was studied, and it is shown that the break-up parameter is not sensitive to the predicted powder particle size. Both simulated results and experimental data showed that fine powders can be produced by increasing disk speed.     

SPH with the multiple boundary tangent method

In this article, we present an improved solid boundary treatment formulation for the smoothed particle hydrodynamics (SPH) method. Benchmark simulations using previously reported boundary treatments can suffer from particle penetration and may produce results that numerically blow up near solid boundaries. As well, current SPH boundary approaches do not properly treat curved boundaries in complicated flow domains. These drawbacks have been remedied in a new boundary treatment method presented in this article, called the multiple boundary tangent (MBT) approach. In this article we present two important benchmark problems to validate the developed algorithm and show that the multiple boundary tangent treatment produces results that agree with known numerical and experimental solutions. The two benchmark problems chosen are the lid‐driven cavity problem, and flow over a cylinder. The SPH solutions using the MBT approach and the results from literature are in very good agreement. These solutions involved solid boundaries, but the approach presented herein should be extendable to time‐evolving, free‐surface boundaries. Copyright © 2008 John Wiley & Sons, Ltd.     

基于SPH方法的油箱抗冲击有限元分析

目的验证SPH方法在流固耦合问题上的适用性。方法采用SPH方法(光滑粒子流体动力学方法)模拟无人机油箱抗冲击过程,并与采用CEL方法(耦合欧拉-拉格朗日方法)模拟得到的结果进行对比。结果 2种模拟方法得到的燃油晃荡情况、油箱应力分布及地面支反力结果接近。结论 SPH方法在流固耦合问题上是适用的。     

基于SPH方法的剪切流驱动液滴在固体表面变形运动数值模拟研究

光滑粒子流体动力学(SPH)方法是纯拉格朗日粒子方法,可以有效避免网格法在模拟大变形过程中带来的网格扭曲等缺陷,适合模拟含大变形的剪切流驱动液滴在固体表面变形运动过程。在基于CSF模型的表面张力SPH方法基础上,采用新的边界处理方式和界面法向修正方法,引入Brackbill提出的壁面附着力边界条件处理方法,得到了含壁面附着力边界条件的表面张力算法。基于新方法模拟了剪切流驱动液滴在固体表面变形运动过程并与实验结果和VOF方法模拟结果进行了对比验证。结果表明:该方法在处理壁面附着力问题时精度较高,稳定性较好,适合处理工程中剪切流驱动液滴在固体表面变形运动问题。     

Simulation of concrete flow in V-funnel test and the proper range of viscosity and yield stress for SCC

The present paper highlights the flow simulation of self consolidating concrete (SCC) in V-funnel test that is used to determine the concrete filling ability and its resistance against segregation. Simulations were performed using a two-dimensional smoothed particle hydrodynamic (SPH) method to determine the discharge time where SCC was considered as a homogeneous Bingham fluid. The numerical predictions are lower than experimental data because of the assumptions of two-dimensional and homogeneous flow. Having the SPH method employed, SCCs with different viscosities and yield stresses were simulated to compare the discharge time with the suggested criteria in EFNARC (2002) and (2005) guidelines. Based on simulations results, the appropriate range of viscosities and yield stresses as well as a relation between rheological properties and discharge time for SCC taking into account EFNARC (2002) and (2005) guidelines are suggested. Using the suggested relations, one can assess the proper SCC filling ability without conducting the V-funnel test.     

POWDER COATING PROCESS PARAMETERS FOR A TRANSFER EFFICIENCY MODEL

The trajectories of charged powder particles in a powder coating system are governed by the electrostatic, gravitational and aerodynamic forces acting on the particles. A mathematical model of particle trajectories inside a powder coating booth must consider (1) the aerodynamic flow field, (2) particle size and charge distributions, (3) the electrostatic field distribution, and (4)the geometry of the target. Our approach is to employ a grid generation and flow solver to examine the air flow pattern and an iterative technique where the Charge Simulation Method can be used to compute the electric field strength and the Method of Characteristics can be used to compute the charge density in the gun-to-target region. The electrostatic forces due to the deposited powder layer and image charge are to be taken into account to determine if the particle will deposit on the substrate or not. The model is applied to the geometry of a high-voltage electrode consisting of a long thin rod with a hemispherical end cap and a grounded flat disk substrate. An experimental system to measure transfer efficiency, with the ability to control various parameters effecting transfer efficiency, has been developed to verify the theoretical model. The simulation results can provide valuable information concerning particle deposition and optimization of transfer efficiency. This paper describes (1) system parameters involved in modeling the transfer efficiency, (2) an approach to develop such a model with preliminary data on the simulation of particle track, and (3) experimental data on the real-time measurements of first pass transfer efficiency.     

SPH原理、发展现状及热传导问题模型

对光滑粒子流体动力学方法（smoothed particle hydrodynamics,SPH）的基本技术原理及发展现状进行了综述,利用SPH方法对两个典型二维非线性动力学算例进行了数值模拟。同时,利用国外学者提出的方法初步探索了有关热传导问题的SPH模型。     

The numerical investigation of spreading process of two viscoelastic droplets impact problem by using an improved SPH scheme

In this paper, the spreading process of two XPP model droplets impacting on a plate in sequence at low Reynolds number is numerically simulated by using an improved smoothed particle hydrodynamics (I-SPH) method. The I-SPH method is a coupled approach which uses the traditional SPH (TSPH) method near the boundary domain and uses a kernel-gradient-corrected SPH method in the interior of fluid flow for the reason of remedying the accuracy and stability of TSPH. Meanwhile, an artificial stress term and a periodic density re-initialization technique are presented to eliminate the tensile instability and restrain pressure oscillation, respectively. A new boundary treatment is also adopted. The ability and merit of proposed I-SPH method combined with other techniques are first illustrated by simulating three typical examples. Subsequently, the deformation phenomena of two viscoelastic droplets impacting and spreading on a plate in sequence are numerically investigated. Particularly, the influences of the falling time interval, Weissenberg number and other rheological parameters on the deformation process are studied respectively. All numerical results agree well with the available data.     

Dynamic flow‐based particle splitting in smoothed particle hydrodynamics

In this work, a dynamic procedure for local particle refinement to be used in smoothed particle hydrodynamics (SPH) is presented. The algorithm is able to consistently produce successive levels of particle splitting in accordance to a flow‐based criterion. It has been applied together with accurate and robust formulations for variable spatial resolution in the framework of a semi‐implicit, truly incompressible scheme for SPH. Different test cases have been considered to assess the capabilities and advantages of the proposed procedure, namely, the laminar flow around circular and square obstacles in a plane channel for various regimes. Such flow cases entail the simulation of attached and separated shear layers, recirculating flow, vortex shedding and surface discontinuities. The results obtained for two levels of particle splitting have demonstrated that significant improvements may be obtained with respect to uniform particle spacing solutions in a variety of situations, thus presenting an excellent trade‐off between accuracy and computational cost. Copyright © 2015 John Wiley & Sons, Ltd.     

粉末高速压制成形密度分布的数值模拟及影响因素分析

将研究不连续体力学行为的离散单元法应用于粉末高速压制致密化过程的研究,将粉末视为黏弹性的离散颗粒,建立粉末高速压制过程颗粒接触模型及每个颗粒的基本运动方程,推导了力与位移表达的粉末高速压制黏弹性本构关系。基于PFC软件实现了铁粉高速压制过程中粉末颗粒二维流动情况及压坯密度分布的数值模拟,模拟结果的密度分布规律与实际压制的密度分布规律较为一致;利用数值模拟结果对影响压坯密度分布的摩擦因数、高径比、双向压制因素进行了具体分析。     

Numerical simulation of reactive polymer flow during rotational molding using smoothed particle hydrodynamics method and experimental verification

A numerical study of flow behavior and the heat transfer process of reactive polymer in reactive rotational molding systems is carried out using the Smoothed Particle Hydrodynamics method (SPH). The flow of reactive polymer during rotational molding inside a horizontal rotating cylinder is modeled as the slightly compressible viscous fluid flow and with free surface. These simulations show the influence of process parameters on the flow of polymer. Especially, the influence of the change of viscosity on the flow, due to the chemical reactions, is simulated by using an adapted viscosity expression. For a constant rotational speed, several flow regimes are observed as follows: the fluid remains at the bottom of the mold with the formation of a thin layer on the surface if the viscosity of reactive polymer is very low (less than 850 Pa.s), the presence of cascades is shown (after the point of maximum height) if the viscosity is higher (850 to 980 Pa.s), and the fluid moves at the same speed as the mold if the viscosity is sufficiently high (approximately 1230 Pa.s). The heat transfer between the mold and reactive polymer is also simulated. The results of SPH simulations are then compared to the experimental results conducted on polyurethane (TPU) system based on 1,6-hexamethylene diisocyanate (HDI, Sigma Aldrich, France), polyethylene glycol (PEG 1000) as macrodiol, 1,3-propanediol (PDO), and dibutyl tin dilaurate (DBTDL).     

A Lagrangian gradient smoothing method for solid‐flow problems using simplicial mesh

A novel Lagrangian gradient smoothing method (L‐GSM) is developed to solve “solid‐flow” (flow media with material strength) problems governed by Lagrangian form of Navier‐Stokes equations. It is a particle‐like method, similar to the smoothed particle hydrodynamics (SPH) method but without the so‐called tensile instability that exists in the SPH since its birth. The L‐GSM uses gradient smoothing technique to approximate the gradient of the field variables, based on the standard GSM that was found working well with Euler grids for general fluids. The Delaunay triangulation algorithm is adopted to update the connectivity of the particles, so that supporting neighboring particles can be determined for accurate gradient approximations. Special techniques are also devised for treatments of 3 types of boundaries: no‐slip solid boundary, free‐surface boundary, and periodical boundary. An advanced GSM operation for better consistency condition is then developed. Tensile stability condition of L‐GSM is investigated through the von Neumann stability analysis as well as numerical tests. The proposed L‐GSM is validated by using benchmarking examples of incompressible flows, including the Couette flow, Poiseuille flow, and 2D shear‐driven cavity. It is then applied to solve a practical problem of solid flows: the natural failure process of soil and the resultant soil flows. The numerical results are compared with theoretical solutions, experimental data, and other numerical results by SPH and FDM to evaluate further L‐GSM performance. It shows that the L‐GSM scheme can give a very accurate result for all these examples. Both the theoretical analysis and the numerical testing results demonstrate that the proposed L‐GSM approach restores first‐order accuracy unconditionally and does not suffer from the tensile instability. It is also shown that the L‐GSM is much more computational efficient compared with SPH, especially when a large number of particles are employed in simulation.     

SPH方法在剪切式碰撞能量吸收器中的应用

剪切式碰撞能量吸收器是提高汽车被动安全性的一种实用新型设计。建立了剪切式碰撞能量吸收器数值分析模型,采用SPH方法数值模拟了碰撞吸能器的碰撞过程,研究了其碰撞吸能特性,并通过不同碰撞速度的台车碰撞试验进行了试验验证研究。数值模拟结果同试验结果相符,表明SPH方法在碰撞吸能器性能研究中是行之有效的数值计算方法。同时分析了在碰撞吸能器的复杂碰撞情况下,SPH方法相比有限元法的优势。     

螺旋送粉器内流场的三维数值模拟

为了给螺旋送粉器内流道优化提供理论依据,根据气-固两相流输送的相关理论,运用Eulerian-Lagrange数学方式,采用FLUENT软件,通过离散相模型对螺旋送粉器内固体颗粒的轨迹进行数值模拟,主要研究螺旋送粉器的转速对其输送性能的影响,及不同直径粉体颗粒轨迹。结果表明:气体速度和粉体颗粒大小是影响螺旋送粉器送粉性能的关键因素。     

大尺寸RPC管-HSC组合柱轴压性能试验研究

在活性粉末混凝土(RPC)预制管内部浇筑高强混凝土,形成RPC管-高强混凝土(HSC)组合柱,扩展了RPC管-混凝土组合柱(CFRT)这一新型组合结构的范围。开展了RPC管离心法成型的配合比与工艺研究,成功试制了RPC管。对4组大尺寸CFRT和1组箍筋约束高强混凝土柱试件进行轴压试验,试验参数为内部混凝土和纵筋配置率。结果表明:内部填充HSC的CFRT柱在整个受力过程中,预制管身基本完整,抗压性能显著优于用于对比的箍筋约束高强混凝土柱;CFRT柱的抗压强度随内部混凝土强度的提高而提高,但约束效应逐步降低,对CFRT柱的内部混凝土的强度应该有所限制;在内部混凝土中配置纵向钢筋,对CFRT的延性影响不大;基于Mander(1988)约束模型,提出了内部填充高强混凝土的CFRT柱轴向承载力计算方法,模型的预测结果与试验结果吻合较好。