VOF simulations of the contact angle dynamics during the drop spreading: Standard models and a new wetting force model

Introduction

In this study,a novel numerical implementation for the adhesion of liquid droplets impacting normally on solid dry surfaces is presented. The advantage of this new approach, compared to the majority of existing models, is that the dynamic contact angle forming during the surface wetting process is not inserted as a boundary condition, but is derived implicitly by the induced fluid flow characteristics (interface shape) and the adhesion physics of the gas–liquid-surface interface (triple line), starting only from the advancing and receding equilibrium contact angles. These angles are required in order to define the wetting properties of liquid phases when interacting with a solid surface.

Methodology

The physical model is implemented as a source term in the momentum equation of a Navier-Stokes CFD flow solver as an “adhesion-like” force which acts at the triple-phase contact line as a result of capillary interactions between the liquid drop and the solid substrate. The numerical simulations capture the liquid–air interface movement by considering the volume of fluid (VOF) method and utilizing an automatic local grid refinement technique in order to increase the accuracy of the predictions at the area of interest, and simultaneously minimize numerical diffusion of the interface.

Results

The proposed model is validated against previously reported experimental data of normal impingement of water droplets on dry surfaces at room temperature. A wide range of impact velocities, i.e. Weber numbers from as low as 0.2 up to 117, both for hydrophilic (θ_adv = 10° – 70°) and hydrophobic (θ_adv = 105° – 120°) surfaces, has been examined. Predictions include in addition to droplet spreading dynamics, the estimation of the dynamic contact angle; the latter is found in reasonable agreement against available experimental measurements.

Conclusion

It is thus concluded that theimplementation of this model is an effective approach for overcoming the need of a pre-defined dynamic contact angle law, frequently adopted as an approximate boundary condition for such simulations. Clearly, this model is mostly influential during the spreading phase for the cases of low We number impacts (We < ˜80) since for high impact velocities, inertia dominates significantly over capillary forces in the initial phase of spreading.     

Separation efficiency of liquid–solid undergoing vibration based on breakage of liquid bridge

Vibrating separation is a significant method for liquid–solid separation. A typical example is the vibrating screen to dewater wet granular matter. The properties of granular matter and the vibrating parameters significantly affect the separation efficiency. This study investigates the effect of vibration parameters in separation based on the breakage of large-scale liquid bridge numerically by using a calibrated simulation model. Through analysing the simulation results, the liquid bridge shape and the volume between two sphere particles for various particle sizes and particle distances were studied in the static condition under the effect of gravity. The results show a general reducing trend of liquid bridge volume when the radius ratio of two particles increases, particularly when the ratio increases to 5. Additionally, a set of vibrating motion was applied to the liquid bridge in the simulation model. A group of experiments were also performed to validate the simulation model with vibration. Then, the effect of vibrating peak acceleration, distance between spheres and radius on the separation efficiency which was reflected by the residual water were investigated. It is found that separation efficiency increased obviously with the peak acceleration and the increase slowed down after the peak acceleration over 1 m/s².     

On the thermodynamic behaviors and interactions between bubble pairs: A numerical approach

Thermodynamic behaviors and interactions between bubble pairs are important to better understand the cavitation phenomena. In this study, a compressible two-phase model, accounting for thermal effects to investigate the thermodynamic behaviors and interactions between bubble pairs, is developed in OpenFOAM. The volume of fluid (VOF) method is adopted to capture the interface. Validations are performed by comparing the simulation results of a single bubble and bubble pairs with corresponding experimental data. The dynamical behaviors of bubble pairs and their thermodynamic effect at different relative distances γ are investigated and discussed, which help reveal the bubble cloud dynamics. The quantitative analysis of γ effects on the maximum temperature during bubble collapse is performed with three distinct stages identified. For a single bubble collapsing near the rigid surface, the thermodynamic characteristics at different relative distances are similar to that of the bubble pairs, but the maximum temperature is higher since the single bubble can collapse to a smaller volume.     

Development of an efficient statistical volumes of fluid–Lagrangian particle tracking coupling method

The breakup of a liquid jet into irregular liquid structures and droplets leading to the formation of a dilute spray has been simulated numerically. To overcome the shortcomings of certain numerical methods in specific flow regimes, a combined approach has been chosen. The intact liquid core, its primary breakup and the dense spray regime are simulated using the volumes of fluid (VOF) method in combination with LES, whereas the Lagrangian particle tracking (LPT) approach in the LES context is applied to the dilute spray regime and the secondary breakup of droplets. A method has been developed to couple both simulation types on a statistical basis. This statistical coupling approach (SCA) reflects the dominating physical mechanisms of the two‐phase flow in each regime to a high degree. The main benefit of the SCA is computational efficiency as compared with the more straightforward approach where one follows each structure, denoted here as the direct coupling approach. The computational benefits stem from the reduction of computational time since the VOF simulation is run only until statistical convergence and not during the whole spray development. A second benefit using the SCA is the possibility to use the stochastic parcel method in the LPT simulation whereby a large number of droplets may be handled. The coupling approach is applied to the atomization of a fuel jet in a high pressure chamber, demonstrating the gain of efficiency of the SCA as compared with direct coupling approach. Copyright © 2014 The Authors. International Journal for Numerical Methods in Fluids published by John Wiley & Sons Ltd.     

Numerical modeling of 3-D flow on porous broad crested weirs

Gabion weirs with optional design as a broad crested weirs are suitable structures to reduce flash flood with a minimal negative impact on the water environment. In the present study, the 3-D flow was simulated around gabion weirs with respect to free-surface water. The Reynolds-averaged Navier–Stokes equations are solved to predict water surface over the gabion weir. The VOF method with the geometric reconstruction scheme was applied to treat the complex free-surface flow. Simulations were performed using three variants of the k–ε and the RSM models to find the water level and velocity distribution profile and results are compared with several experimental data available in the literature. The structured mesh was used for all domains with high dense mesh near the solid region. A comparison between experimental data and simulations indicates that the k–ε model can be used to predict the complex flow and water level with high accuracy.     

A mesh-dependent model for applying dynamic contact angles to VOF simulations

Typical VOF algorithms rely on an implicit slip that scales with mesh refinement, to allow contact lines to move along no-slip boundaries. As a result, solutions of contact line phenomena vary continuously with mesh spacing; this paper presents examples of that variation. A mesh-dependent dynamic contact angle model is then presented, that is based on fundamental hydrodynamics and serves as a more appropriate boundary condition at a moving contact line. This new boundary condition eliminates the stress singularity at the contact line; the resulting problem is thus well-posed and yields solutions that converge with mesh refinement. Numerical results are presented of a solid plate withdrawing from a fluid pool, and of spontaneous droplet spread at small capillary and Reynolds numbers.     

倾斜上升管内气液两相弹状流壁面传质特性研究

采用VOF模型对倾角为45°、80°、85°三种情况下倾斜上升管内弹状流的壁面传质特性进行了研究.传质特性通过其与壁面切应力的类比关系来体现.数值模拟结果表明,低混合物流速时,上管壁面切应力在液膜区有明显波动,而下管壁面切应力分布则比较光滑.随着混合流速的增大,液膜区上下壁面切应力分布趋于一致.管子下壁面切应力平均值大于管子上壁面,在Taylor气泡运动速度较低时较为突出.随着Taylor气泡速度的增大,管子上下壁面的切应力平均值趋于相同.相同的混合流速下倾斜角度越大,上下管壁的切应力分布越趋于近似.下降液膜区的壁面切应力平均值大于Taylor气泡尾迹区域.根据Chilton-Colburn的类比关系,壁面切应力的规律完全适用于壁面传质系数.     

气相中双液滴正碰过程数值模拟研究

液滴碰撞现象普遍存在于动力装置燃烧室喷嘴的下游区域,影响燃料的雾化性能。为了揭示相同直径的双液滴中心碰撞机理,求解了轴对称坐标系下的N-S方程,采用VOF(Volume of Fluid)方法捕捉液滴碰撞过程中气液自由表面的演化规律。利用Qian等提供的实验结果对计算模型进行数值校验,验证了模型的准确性。在此基础上,研究了环境压强对液滴碰撞反弹后不同结果(分离和融合)的影响,分析了环境压强和Weber数对液滴碰撞分离的影响。结果表明,液滴在碰撞反弹后的状态(分离或融合)是由液滴间气膜压强与环境气动阻力共同作用的结果,环境压强对液滴碰撞分离过程基本没有影响;Weber数越大,碰撞过程中变形的幅度越大。     

流体控制方程的超声空化泡动力学模拟*

本文基于流体动力学控制方程和VOF模型,在FLUENT 14.5软件环境下对超声空化泡进行数值模拟。首先研究了超声空化泡一个周期内的形态变化,并且利用空化泡形态变化的最大面积、最小面积、膨胀时间、收缩时间等数值结果分析超声参数对空化效果的影响。同时探究了双频超声作用下空化泡运动的变化,计算结果表明:在其他条件相同的情况下,在1~5MPa范围内,超声声压幅值为3MPa时空化效果最好;当超声频率大于20kHz时,空化效果随着超声频率的增大而降低。对于频率相同的双频超声,较声压幅值为其两倍的单频超声有更好的空化效果;对于频率不同的双频超声,空化效果受到频率差的影响。     

A stencil-like volume of fluid （VOF） method for tracking free interface

A stencil-like volume of fluid （VOF） method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions are possible in 2D cases, and the interface can be reconstructed by only requiring the known local volume fraction information. On the other hand, the fluid-occupying-length is defined on each side of the stencil, through which a unified fluid-occupying volume model and a unified algorithm can be obtained to solve the interface advection equation. The method is suitable for the arbitrary geometry of the grid cell, and is extendible to 3D cases. Typical numerical examples show that the current method can give ＂sharp＂ results for tracking free interface.