运行VOF法模拟LNG液舱的晃荡特性

液舱晃荡引起的载荷己成为航行中载液船舶安全性评估的重要内容之一、以及船舶结构动力学的一个热点.文中首先运用VOF法建立LNG晃荡液舱模型,给出对LNG液舱晃荡特性数值模拟方法,运用模拟结果讨论不同工况(如装载率、频率、横摇角度幅值)下晃荡动压特性,在此基础上给出了LNG液舱的晃荡载荷特性与抑制晃荡的措施.     

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

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

Finite element analysis of flow and heat transfer with moving free surface using fixed grid system

A numerical study was performed on flow and heat transfer involving moving free surfaces that occurs in mold filling processes such as casting and injection molding. In these problems, the calculation domain changes continuously and the numerical treatment of the moving interface tends to cause artificial diffusion. Among the solution algorithms based on the Eulerian method, the volume-of-fluid (VOF) method was used because the method is simple and efficient in handling the complex flow patterns inside the cavity. To solve the transport equation of free surface without artificial smearing of the interface the baby-cell method was employed in the geometric reconstruction of the free surface. Furthermore, a predictor–corrector method was adopted in the time integration of volume-of-fluid (VOF) transport equation to increase the accuracy. The proposed scheme was verified through several benchmark problems. In order to show the capability of the proposed method, several three-dimensional mold filling processes were solved. The current algorithm was applied to the floating body problem. Three-dimensional floating body problems were tested.     

Film flow around a fast rotating roller

In this study, the film thickness around the roller is numerically estimated using the volume of fluid (VOF) method to clarify the film-formation process around the rotating roller. Parametric studies were performed to compare the effects of ink properties (viscosity, surface tension) and operational conditions (roller rotation speed, initial immersed angle) on film thickness. The viscosity of the ink and the speed of rotation of the roller were found to be the dominant factors that determine the ink film thickness. In addition, a correlation equation is proposed to predict the thickness of the ink film around a printing roller rotating at a speed of 20–30 rad/s, as a function of angular position, angular velocity, and viscosity.     

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.     

Modelling and simulation of wave transformation in porous structures using VOF based two-phase flow model

This paper represents the results of wave transformation in porous structures and hydraulic performance of a vertical porous seawall. The study was carried out using a VOF based two-phase numerical hydrodynamic model. The model was developed by coupling an ordinary porous flow model based on extended Navier–Stokes equations for porous media, and a two-phase flow model. A unique solution domain was established with proper treatment of the interface boundary between water, air and the structure. The VOF method with an improved fluid advection algorithm was used to trace the interface between water and air. The resistance to flow caused by the presence of structural material was modeled in terms of drag and inertia forces. The parameters that govern resistance to flow in a porous media were calibrated for a typical structural setup and then the computational efficacy of the model was evaluated for several wave and structural conditions other than the calibrated setup. A set of comparisons of wave properties in and around the structure showed that the model reproduced reasonably good agreement between computed results and measured data. The model was then applied to investigate wave transformation in a vertical porous structure. The role of porosity and width of a structure in reducing wave reflection and increasing energy dissipation was investigated. It is confirmed that there exists an optimum value of structure width and porosity that can maximize hydraulic performances of a porous seawall.     

振荡液滴内部流态

在外界来流作用下,液滴在固体表面上呈现周期性振荡特性.利用数值方法模拟平板上二维液滴在气流剪切作用下的界面及内部流动特性,重构二维液滴内部流场,着重认识液滴内部速度分布和压力分布.     

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.