一种基于几何重构的Youngs-VOF耦合水平集追踪方法

针对界面追踪方法中拉格朗日方法和欧拉--拉格朗日方法计算效率低、不适用大变形、不能应用于三维数值计算模型等问题,研究了一种效率高、界面清晰、适用于三维模型的计算气液两相界面迁移特性的欧拉运动界面追踪方法,该方法将"米"状相邻单元Youngs方法用于运动界面重构,将Youngs-VOF和水平集通过几何方法耦合,提高运动界面精度,克服了VOF和水平集方法存在的缺陷,避免了利用高阶导数本身的稳定性去求解水平集对流方程和距离函数方程."米"状相邻单元Youngs方法避免了数值耗散、数值色散性以及非线性效应引起的捕捉界面模糊的情况.Youngs-VOF耦合水平集方法既保证了计算界面时的稳定性,与拉格朗日方法相比又提高了计算效率.利用Youngs-VOF耦合水平集方法与VOF方法对单个气泡在水中上升过程数值计算与实验对比并对经典剪切流场中圆形运动界面模型的数值计算,验证了Youngs-VOF耦合水平集方法的有效性并比VOF方法捕捉界面更清晰、锐利;通过对溃坝--自由表面流动过程数值计算并与实验进行对比,验证了Youngs-VOF耦合水平集方法的稳定性以及对三维数值模型的适用性.      

孤立波与多个潜水障碍物相互作用的数值研究

本文以完整二维Navier-Stokes方程为控制方程,采用VOF界面跟踪技术和差分方法,数值计算孤立波与多个淹没水下物体相互作用的问题。本文对潜水物体的高度接近水深时,孤立波通过水下孤立直立方柱、两个间距较大的水下直立方柱和两个间距较小的水下直立方柱等三种情况分别进行了计算,给出了波形随时间的演化图,可以看到反射波、前传波和跟随振荡型小波列的生成。对孤立波通过水下孤立直立方柱情形的计算结果,与实验结果和势流理论结果进行了比较。     

Modelling of steel scrap movement

Uneven movement of steel scrap during electric arc furnace melting causes operational problems, such as cave-ins. To understand the phenomenon, a small-scale apparatus has been built to simulate the flow of scrap that may be interlocked. It was found that flow starts at the point where the suspended length of the pile bottom reaches the average scrap size. The force balance in a suspended scrap pile has been analyzed to estimate the effective cohesion of interlocking. A numerical model has been developed to simulate the scrap flow.     

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.     

Numerical investigation of the deformation mechanism of a bubble or a drop rising or falling in another fluid

A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） method in the frame of two incompressible and immiscible viscous fluids under the action of gravity, taking into consideration of surface tension effects. A comparison of the numerical results by this method with those by other works indicates the validity of the method. In the frame of inviseid and incompressible fluids without taking into consideration of surface tension effects, the mechanisms of the generation of the liquid jet and the transition from spherical shape to toroidal shape during the bubble or drop deformation, the increase of the ring diameter of the toroidal bubble or drop and the decrease of its cross-section area during its motion, and the effects of the density ratio of the two fluids on the deformation of the bubble or drop are analysed both theoretically and numerically.     

Application of the high-resolution Godunov method to the multi-fluid flow calculations

In this paper, we have numerically solved the multi-fluid problems using an operator-split two-step high-resolution Godunov PPM (parabolic piecewise method) for the flow in complex geometries. By using the front capturing method,the PPM integrator captures the interface in the solution process. The basic multi-fluid integrator is coupled to a Cartesian grid algorithm where a VOF (volume of fluid) representation of the fluid interface is also used. As an application of this method, we test the 2D interracial advection example and simulate an experimental hypervelocity launcher model from Sandia National Laboratories. The computational design of the hypervelocity launcher is also given in the paper.     

Numerical simulations of gas–liquid–solid flows in a hydrocyclone separator

The flow behavior in hydrocyclones is quite complex. In this study, the computational fluid dynamics (CFD) method was used to simulate the flow fields inside a hydrocyclone in order to investigate its separation efficiency. In the computational fluid dynamics study of hydrocyclones, the air-core dimension is a key to predicting the mass split between the underflow and overflow. In turn, the mass split influences the prediction of the size classification curve. Three models, the model, the Reynolds stress model (RSM) without considering the air-core, and the Reynolds stress turbulence model with the volume of fluid (VOF) multiphase model for simulating the air-core, were compared in terms of their predictions of velocity, axial and tangential velocity distributions, and separation proportion. The RSM with air-core simulation model, since it reproduces some detailed features of the turbulence and multiphase, clearly predicted the experimental data more closely than did the other two models.     

A new volume‐preserving and continuous interface reconstruction method for 2D multi‐material flow

A new two‐dimensional interface reconstruction method that ensures continuity of the interface and preserves volume fractions is presented here. It is made of two steps, first, the minimization of a cost functional based on volume fractions least square errors by using dynamic programming, a fast and efficient scheme well known in image processing, and then a local correction phase. In each cell, the interface is made of two line segments joining two edges. This new interface reconstruction method, called Dynamic Programming Interface Reconstruction has been coupled with various advection schemes, among them the Lagrange + remap scheme. With a reasonable computational cost, it has been observed in various test cases that Dynamic Programming Interface Reconstruction is more accurate and less diffusive compared with other existing classical reconstruction methods. Copyright © 2017 John Wiley & Sons, Ltd.     

An analytical interface reconstruction algorithm in the PLIC‐VOF method for 2D polygonal unstructured meshes

Piecewise linear interface calculation (PLIC) schemes have been extensively employed in the volume‐of‐fluid (VOF) method for interface capturing in numerical simulations of multiphase flows. Polygonal unstructured meshes are often adopted because of their geometric flexibility and superiority in gradient calculation. An analytical interface reconstruction algorithm in the PLIC‐VOF method for arbitrary convex polygonal cells has been proposed in this study. The line interface at a given orientation within a polygonal cell is located by an analytical technique. It has been tested successfully for four different geometric shapes that are common in polygonal meshes. The computational efficiency of the present algorithm has been compared with several published schemes in the literature. The proposed algorithm has been shown to yield higher accuracy with reduction in computational complexity. A numerical simulation of a dam‐breaking problem has been performed using the proposed analytical PLIC technique on polygonal meshes. The results are in good agreement with experimental data available in the literature, which serves as a demonstration of its performance in a real multiphase flow.     

An evaluation of interface capturing methods in a VOF based model for multiphase flow of a non-Newtonian ceramic in tape casting

The aim of the present study is to evaluate the different interface capturing methods as well as to find the best approach for flow modeling of the ceramic slurry in the tape casting process. The conventional volume of fluid (VOF) method with three different interpolation methods for interface capturing, i.e. the Geometric Reconstruction Scheme (GRS), High Resolution Interface Capturing (HRIC) and Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM), are investigated for the advection of the VOF, both for Newtonian and non-Newtonian cases. The main purpose is to find the best method for the free surface capturing during the flow of a ceramic slurry described by a constitutive power law equation in the tape casting process. First the developed model is tested against well-documented and relevant solutions from literature involving free surface tracking and subsequently it is used to investigate the flow of a La_0.85Sr_0.15MnO₃ (LSM) ceramic slurry modeled with the Ostwald de Waele power law. Results of the modeling are compared with corresponding experimental data and good agreement is found.