Modeling surface tension of a two‐dimensional droplet using smoothed particle hydrodynamics

Surface tension plays a significant role at the dynamic interface of free‐surface flows especially at the microscale in capillary‐dominated flows. A model for accurately predicting the formation of two‐dimensional viscous droplets in vacuum or gas of negligible density and viscosity resulting from axisymmetric oscillation due to surface tension is solved using smoothed particle hydrodynamics composed of the Navier‐Stokes system and appropriate interfacial conditions for the free‐surface boundaries. The evolution of the droplet and its free‐surface interface is tracked over time to investigate the effects of surface tension forces implemented using a modified continuous surface force method and is compared with those performed using interparticle interaction force. The dynamic viscous fluid and surface tension interactions are investigated via a controlled curvature model and test cases of nonsteady oscillating droplets; attention is focused here on droplet oscillation that is released from an initial static deformation. Accuracy of the results is attested by demonstrating that (i) the curvature of the droplet that is controlled; (ii) uniform distribution of fluid particles; (iii) clean asymmetric forces acting on the free surface; and (iv) nonsteady oscillating droplets compare well with analytical and published experiment findings. The advantage of the proposed continuous surface force method only requires the use of physical properties of the fluid, whereas the interparticle interaction force method is restricted by the requirement of tuning parameters.     

绕水翼超空化流动形态与速度分布

为揭示超空化流场结构特性,利用高速全流场显示技术,观察了绕hydronautics水翼的超空化流动形态,并利用数字粒子图像测速仪(DPIV)测量了其速度分布. 在测量空穴内部流速分布时,采用空化流场中的空化泡作为示踪粒子来显示流动结构. 结果表明:随着空化数的降低,超空化流动显现出了明显的阶段特征,其中水汽混合相和汽相的分布决定了空化区域的形态与流速分布;空化区和主流区的汽液交界面处存在着较大的速度梯度;低速分布区域随着空化数的降低由水翼吸力面中后部向水翼下游移动;在空化区域内部,水汽混合区的速度相对较低,而汽相区则与主流区有着相近的速度分布.关键词超空化水翼、DPIV、高速摄像、空化形态、流速分布      

Oldroyd-B黏弹性液滴弹跳行为的改进SPH模拟

基于光滑粒子流体动力学SPH(Smoothed Particle Hydrodynamics)方法对Oldroyd-B黏弹性液滴撞击固壁面产生的弹跳行为进行了模拟与分析。首先,为了解决SPH模拟黏弹性自由表面流出现的张力不稳定性问题,联合粒子迁移技术提出了一种改进SPH方法。然后,对Oldroyd-B黏弹性液滴撞击固壁面产生的铺展行为进行了改进SPH模拟,与文献结果的比较验证了方法的有效性。最后,通过降低Reynolds数捕捉到了液滴的弹跳行为;并在此基础上,分析了液滴黏度比、Weissenberg数和Reynolds数对液滴弹跳行为的影响。结果表明,改进SPH方法可有效地模拟黏弹性自由表面流问题;液滴黏度比、Weissenberg数和Reynolds数对液滴最大回弹高度均有显著的影响。     

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

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

基于表面力模型的液滴冲击弹性基底SPH模拟

采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s～3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。     

SPH固壁边界处理方法的改进

提出一种适用于光滑质点水动力学（SPH）方法的改进的边界处理方法。在这种方法中,边界粒子的压力可通过其周围的流体粒子的压力插值得到,从而改进了耦合边界法在边界上压力不准的问题。运用这种改进的边界处理方法模拟了二维方形水箱中的非线性晃荡问题以及二维楔形体自由入水问题。模拟结果与实验结果吻合较好,证明了此改进的边界处理方法是有效的。     

Velocity field measurements of cavitating flows

A particle Image Velocimetry (PIV) system has been developed to study the microfluid mechanics of cavitating flows. Planar PIV was used to examine the non-cavitating flow in the thin boundary layer near a hydrofoil surface for the cases of a naturally developing boundary layer and a boundary layer stimulated to turbulence by roughness near the foil leading edge. PIV was also used to examine the flow near the surface of individual cavitation bubbles and incipient attached cavitation. A system was devised to create a single nucleus in the flow upstream of a hydrofoil, and planar PIV was used to study the flow around the resulting traveling cavitation bubble. Velocity vectors were determined close to the solid surfaces and the gas/liquid interfaces of the bubbles. Seeding of the flow with particles did not result in the addition of active cavitation nuclei.     

Large‐scale SPH simulations of droplet impact onto a liquid surface up to the consequent formation of Worthington jet

In this study, the whole process of liquid droplet impact onto a liquid surface up to the consequent formation of the central column was simulated using the smoothed particle hydrodynamics method (SPH), and compared with an experiment using a high‐speed video camera. The surface tension tensor for the particle‐based expression was adequately included as the gradient of the surface tension and that enabled the simulation leading to the formations of crater and crown as well as the consequent central column. The simulated time series of the crater depth and diameter and crown height corresponded quantitatively well with the experimental result up to the rebound motion while discrepancies remained as a lower central column height in the simulation, and this seemed to be ascribed to the difficulty in realizing the complex surface structure that inevitably appeared in the fast rebound motion. Copyright © 2009 John Wiley & Sons, Ltd.     

空化水动力学非定常特性研究进展及展望

空化作为一种重要的复杂水动力学现象,具有明显的三维流动特征与剧烈的非定常特性,在水力机械、船舶推进器、水利工程中广泛存在,且通常会带来不利的影响,长期以来一直是水动力学领域研究的重点与难点课题之一.本文首先从实验测量和数值模拟两个角度,综述了空化水动力学非定常特性研究的发展概况, 分析了当前存在的问题.在空化实验研究中,主要介绍了空化水洞、空化流场测量以及多物理场同步测量等方面所取得的进展.在数值模拟方法中, 对目前的空化模型和湍流模型进行了分类介绍,并重点讨论了大涡模拟、验证和确认等在空化流模拟中的应用.之后以附着型空化为主, 同时兼顾云状空泡、空蚀、涡空化等,梳理了其研究中存在的几个关键科学问题,包括空化演变、空化流动的三维结构、失稳机制、空化不稳定性及其与低频压力脉动的联系、空化与旋涡的相互作用、空化与弹性水翼的流固耦合、空化对尾流场影响等.最后展望了空化水动力学的研究方向和未来发展趋势.      

Research progresses and prospects of unsteady hydrodynamics characteristics for cavitation

空化作为一种重要的复杂水动力学现象,具有明显的三维流动特征与剧烈的非定常特性,在水力机械、船舶推进器、水利工程中广泛存在,且通常会带来不利的影响,长期以来一直是水动力学领域研究的重点与难点课题之一.本文首先从实验测量和数值模拟两个角度,综述了空化水动力学非定常特性研究的发展概况, 分析了当前存在的问题.在空化实验研究中,主要介绍了空化水洞、空化流场测量以及多物理场同步测量等方面所取得的进展.在数值模拟方法中, 对目前的空化模型和湍流模型进行了分类介绍,并重点讨论了大涡模拟、验证和确认等在空化流模拟中的应用.之后以附着型空化为主, 同时兼顾云状空泡、空蚀、涡空化等,梳理了其研究中存在的几个关键科学问题,包括空化演变、空化流动的三维结构、失稳机制、空化不稳定性及其与低频压力脉动的联系、空化与旋涡的相互作用、空化与弹性水翼的流固耦合、空化对尾流场影响等.最后展望了空化水动力学的研究方向和未来发展趋势.     

Numerical simulation of droplet coalescence based on the SPH method

In this paper, the smoothed particle hydrodynamics (SPH) method is employed in modeling and numerical simulation of droplet coalescence. Considering the effect of tangential force on boundary material, besides normal force, tangential force is also introduced in the continuum surface force (CSF) model. The formation of droplet, the coalescence processes of two droplets and three droplets are simulated by the modified CSF model. The validity of the modified model is verified from the aspects of the morphological change of the droplet, the smoothness of free surface and the conservation of the centroid of the system. Compared with finite element method, the results of the modified CSF model show that tangential force plays a crucial role in the CSF model when dealing with model boundary with curves and sharp angles.     

绕水翼空化流动及振动特性的实验研究

空化是发生在水力机械内部的一种水动力现象,其发展具有显著的非定常特性.空化流动中空穴的脱落以及溃灭会诱发结构振动,对水力机械的效率、噪声、安全性等造成影响. 研究空化流动中结构的振动特性具有重要的工程意义. 采用实验的方法研究了绕NACA66 水翼空化流动的空穴形态和水翼振动特性. 实验在一闭式空化水洞中进行. 采用高速摄像技术观测不同空化阶段的空穴形态,应用多普勒激光测振仪测量水翼的振动速度,并通过一套同步系统实现了高速相机和多普勒激光测振仪的同步触发和测量. 采用小波分析方法对不同空化阶段下的空穴形态和水翼振动数据在时域和频域中的特性进行了分析.对云状空化阶段的同步测量结果进行了研究,分析了振动与空穴发展过程的联系. 结果表明,随着空化数的降低,流场经历了无空化、初生空化、片状空化和云状空化4个阶段,水翼的振动强度呈逐渐增大趋势. 在片状空化和云状空化阶段,空穴脱落导致水翼振动,诱发的振动频率与空穴脱落频率相同. 对于云状空化,在附着型空穴生长阶段水翼发生高频小幅度振动,在空穴脉动和断裂脱落期间水翼表现为低频大幅振动.      

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.     

Rheological model of volume stretching of newtonian liquids

An equation of state for a volumetrically stretched cavitating liquid medium that holds in the entire span of volume concentrations of bubbles ranging from cavitation nuclei to the stage of formation of a cellular foam structure is obtained based on a proposed macrorheological model. The dependence of the modulus of volume elasticity of a liquid on the volume concentration of bubbles is plotted, and a method for estimating the relaxation time for tensile stresses in cavitating liquid media is proposed. Lavrent'ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 1, pp. 35–46, January–February, 1998.     

侧斜与负载对螺旋桨无空化和空化水动力性能的影响

为了定量分析空化初始航速的影响因素,首先分析侧斜和负载对螺旋桨无空化和有空化时性能的影响。以NSRDC4381无侧斜桨和4383100％侧斜桨为对象,采用改进Sauer空化模型和修正SST湍流模型,对空化崩溃性能、空化初生性能和无空化时正车、倒车以及紧急倒车敞水性能进行了计算与比较。结果表明,预报两个桨的敞水性能曲线和多空化数下的空化崩溃性能曲线均与实验值吻合较好。在中度负载区间（J=0．5-0．9）内,侧斜对正车和紧急倒车时敞水性能以及空化时推力和力矩崩溃性能均无明显影响,但会使倒车敞水性能显著下降。在重载和轻载条件下,侧斜均能明显改善空化初生性能。侧斜一定时,负载会直接影响尾流湍流速度脉动量和涡核集中区分布,影响轴面速度流管收缩程度,进而影响无空化和有空化条件下的推进性能。     

Numerical simulation of multiphase cavitating flows around an underwater projectile

The present simulation investigates the multiphase cavitating flow around an underwater projectile. Based on the Homogeneous Equilibrium Flow assumption, a mixture model is applied to simulate the multiphase cavitating flow including ventilated cavitation caused by air injection as well as natural cavitation that forms in a region where the pressure of liquid falls below its vapor pressure. The transport equation cavitating model is applied. The calculations are executed based on a suite of CFD code. The hydrodynamics characteristics of flow field under the interaction of natural cavitation and ventilated cavitation is analyzed. The results indicate that the ventilated cavitation number is under a combined effect of the natural cavitation number and gas flow rate in the multiphase cavitating flows.     

Numerical study of a single drop impact onto a liquid film up to the consequent formation of a crown

In this paper, the whole dynamic process of a single drop impact onto a thin liquid surface up to the consequent formation of a thin crown is numerically studied using the smoothed particle hydrodynamics (SPH) method. Especially, the gravity, artificial viscosity, and surface tension are introduced into the model. The obtained SPH numerical results are compared with experimental results. The numerical model of the SPH method is valid for simulating the dynamic process of a single drop impact onto a liquid surface. Meanwhile, it is found that the whole dynamic process mainly depends on the depth of the liquid pool and the initial velocity of the droplet.     

Effect of fumed silica nanoparticles on the morphology and rheology of immiscible polymer blends

We examined the effect of interfacially active particles on the morphology and rheology of droplet/matrix blends of two immiscible homopolymers. Experiments were conducted on polybutadiene/polydimethylsiloxane (10/90) blend and the inverse system. The effects of fumed silica nanoparticles, at low particle loadings (0.1–2.0 wt%), were examined by direct flow visualization and by rheology. Fumed silica nanoparticles were found to significantly affect the morphology of polymer blends, inducing droplet cluster structure and decreasing the droplet size, regardless of which phase wets the particles preferentially. This is surprising in light of much past research that shows that particles are capable of bridging and thus induce droplet cluster structure in droplet/matrix systems only when they are preferentially wetted by the continuous phase. Therefore, there should exist other possible mechanisms responsible for these droplet cluster structures except for the bridging mechanism. We proposed a particle-flocculating mechanism based on the fact that fumed silica particles readily flocculate due to their high aspect ratio, fractal-like shape, or interparticle attractions. Optical microscopy also reveals that the clustering structure becomes more extensive, and the droplet sizes in the clusters become smaller when the particle loading is increased. Rheologically, the chief effect of particles is to change the flow behavior from a liquid-like rheology to gel-like behavior. This gel-like behavior can be attributed to droplet clustering. Moreover, it should be emphasized that such gel-like behavior can be seen in the blends regardless of which phase wets the particles preferentially, suggesting that, once again, bridging is not the only cause of droplet clustering.     

DROPLET COLLISION AND COALESCENCE MODEL

A new droplet collision and coalescence model was presented,a quick-sort method for locating collision partners was also devised and based on theoretical and experimental results,further advancement was made to the droplet collision outcome. The advantages of the two implementations of smoothed particle hydrodynamics (SPH) method were used to limit the collision of droplets to a given number of nearest droplets and define the probability of coalescence,numerical simulations were carried out for model validation.Results show that the model presented is mesh-independent and less time consuming,it can not only maintains the system momentum conservation perfectly,but not susceptible to initial droplet size distribution as well.