第9届国际两相流计算讨论会综述

国际两相流计算讨论会专门交流气固、气液、液固两相流动的数学模型,数值模拟方法,检验数模的测量技术,以及应用方面的最新成果。会议中用某些事先不通知的标准实验数据检验不同研究者的数学模型,公布其对照结果,加以评价。第9届会议于1999年4月13～16日在德国Merseburg城的Martin-Luther大学召开,与会代表有47人。会议发表论文42篇,中国清华大学周力行在会上发表论文2篇。

数值模拟两相流求解器interFoam的应用

两相流的数值模拟是近年来的计算流体力学(Computational Fluid Dynamics,简称CFD)领域研究的热点和前沿课题之一,inter Foam求解器是目前流行计算流体力学软件Open FOAM中能够模拟两相流问题,计算水气相互作用的一个两相流求解器,继承了C++语言面向对象的技术特点,得到了广泛的应用。

水煤浆气化文丘里管压降的计算分析

本文介绍了水煤浆气化文丘里管压降的组成部分和计算方法,并结合工程实例计算分析了文丘里管各部分压降的比重,得出了主要压降为液滴加速压降、喉管两相流动摩擦压降和收缩管扩张管的局部压降。

底吹脱氧、脱硫工艺的数值计算及试验分析

介绍了钢包底吹系统基本原理，建立钢包三维数学模型，进行Euler两相流计算，理论分析了底吹脱氧、脱硫的动力学条件。结合实验室底吹脱氧、脱硫热态试验，分析了钢包底吹脱氧及深脱硫的可行性及高效性。结果表明该脱氧、脱硫方法不仅能充分利用钢包底吹系统、节约生产成本，也为生产优质钢提供了新的工艺路线。

化工装置中再沸器安全阀的设计及选用

安全阀作为一种重要的安全装置,在各种工业领域中都有广泛的应用。它们在保护设备和人员安全方面发挥着关键作用。文章详细介绍了安全阀在化工装置中的应用,并探讨了在具体工作条件、介质特性的要求下正确计算和选用安全阀的方法。在对安全阀的类型、计算方法和选用标准进行深入分析的基础上提出了实用的建议,以确保装置的安全运行。文章同时还探讨了两相流介质中安全阀的计算,通过对单独计算和欧米伽法计算出来的气液两相的泄放面积进行对比,得出相对可靠的选用依据,可为今后相似工艺装置中的安全阀设计提供参考。

Lubrication properties of Nitinol 60 alloy used as high-speed rolling bearing and numerical simulation of flow pattern of oil-air lubrication

The friction and wear tests were performed using Nitinol 60 alloy pin sliding over GCr15 steel disk in a pin-on-disk tribometer system under PAO oil lubrication conditions. It was found that Nitinol 60 alloy can be lubricated well and has shown remarkable tribological performance. Average coefficient of friction （COF） of Nitinol 60 is 0.6 under dry friction; however, average COF decreases to 0.1 under PAO oil lubrication. SEM image of the worn surface shows that Nitinol 60 exhibits excellent wear resistance and the wear mechanism is mainly adhesive wear. Flow pattern of oil-air flow in oil pipe was simulated by FLUENT software with VOF model for acquiring working performance of oil-air lubrication. The optimum velocity of oil and air at the inlet was achieved, which provides the great proposal for the design of experiment of oil-air lubrication of Nitinol 60 alloy. The simulation results showed that the optimum annular flow of flow pattern was obtained when air velocity is 10 m/s and oil velocity is 0.05 m/s. The formation mechanism of annular flow was also discussed in the present study.

Numerical simulation of two-phase flow field in underwater sealing device based on dynamic mesh

In order to speed underwater launch of minor-caliber weapons,a sealing device can be set in front of underwater muzzle to separate water,preventing the muzzle from water immersion.By establishing and simplifying the model of underwater weapon sealing device and unstructured mesh computing domain model based on computational fluid dynamics（CFD）,dynamic mesh and user defined function（UDF）,the N-S equation is solved and the numerical analysis and calculation of the complex two-phase flow inside the sealing device are carried out.The results show that the gas discharged from the sealing device is conducive to the formation of the projectile supercavity.When the projectile is launched at 5munder water,the shock wave before and after the projectile has impact on the box body up to 100 MPa,therefore the sealing device must be strong enough.The research results have the vital significance to the design of underwater weapon sealing device and the formation of the projectile supercavitation.

Three-dimensional Computational Fluid Dynamics Modeling of Two-phase Flow in a Structured Packing Column

Characterizing the complex two-phase hydrodynamics in structured packed columns requires a power- ful modeling tool. The traditional two-dimensional model exhibits limitations when one attempts to model the de- tailed two-phase flow inside the columns. The present paper presents a three-dimensional computational fluid dy- namics （CFD） model to simulate the two-phase flow in a representative unit of the column. The unit consists of an CFD calculations on column packed with Flexipak 1Y were implemented within the volume of fluid （VOF） mathe- matical framework. The CFD model was validated by comparing the calculated thickness of liquid film with the available experimental data. Special attention was given to quantitative analysis of the effects of gravity on the hy- drodynamics. Fluctuations in the liquid mass flow rate and the calculated pressure drop loss were found to be quali- tatively in agreement with the experimental observations.

Computational fluid dynamics simulation of gas-liquid two phases flow in 320 m^3 air-blowing mechanical flotation cell using different turbulence models

According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.

Analysis of Flow Structure and Calculation of Drag Coefficient for Concurrent-up Gas-Solid Flow

This study investigates the heterogeneous structure and its influence on drag coefficient for concurrent-up gas-solid flow. The energy-minimization multi-scale (EMMS) model is modified to simulate the variation of structure parameters with solids concentration, showing the tendency for particles to aggregate to form clusters and for fluid to pass around clusters. The global drag coefficient is resolved into that for the dense phase, for the dilute phase and for the so-called inter-phase, all of which can be obtained from their respective phase-specific structure parameters. The computational results show that the drag coefficients of the different phases are quite different, and the global drag coefficient calculated from the EMMS approach is much lower than that from the correlation of Wen and Yu. The simulation results demonstrate that the EMMS approach can well describe the heterogeneous flow structure, and is very promising for incorporation into the two-fluid model or the discrete particle model as the closure law for drag coefficient.

Experimental Studies of the Effect of Wall Roughness on Particle Behavior in Gas-Particle Flows

The effect of wall roughness on particle behavior in two-phase flows in a horizontal backward-facing step is studied using a phase-Doppler particle anemometer. The results show that the wall roughness widens the particle velocity probability density distribution, enhances the redistribution of particle velocity into different directions, reduces the particle longitudinal mean velocity and increases the longitudinal and transverse fluctuation velocities and Reynolds shear stress. The effect of roughness on particle motion in the recirculation zone is weaker than that in the fully developed flow region. The effect of roughness for small particles is restricted only in the near-wall region, while that for large particles diffuses to the whole flow field.

Reconstructing bubble profiles from gas-liquid two-phase flow data using agglomerative hierarchical clustering method

The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved computational fluid dynamics (CFD) simulations. To obtain this information, an efficient bubble profile reconstruction method based on an improved agglomerative hierarchical clustering (AHC) algorithm is proposed in this paper. The reconstruction method is featured by the implementations of a binary space division preprocessing, which aims to reduce the computational complexity, an adaptive linkage criterion, which guarantees the applicability of the AHC algorithm when dealing with datasets involving either non-uniform or distorted grids, and a stepwise execution strategy, which enables the separation of attached bubbles. To illustrate and verify this method, it was applied to dealing with 3 datasets, 2 of them with pre-specified spherical bubbles and the other obtained by a surface-resolved CFD simulation. Application results indicate that the proposed method is effective even when the data include some non-uniform and distortion.

A General Thermal Equilibrium Discharge Flow Model

Based on isentropic flow and thermal equilibrium assumptions, a model was derived to calculate discharge flow rate, which unified the rules of room temperature water discharge, high temperature and high pressure water discharge, two-phase critical flow, saturated steam and superheated steam critical flow, and gave a method to calculate critical condition. Because of the influence of friction, the entropy is increased in the actual discharge process, and the discharge flow rate in thermal equilibrium condition can be obtained by the original model multiplied by an appropriate correction coefficient. The model calculated results agreed well with the experiment data of long nozzle critical flow.

CFD-DEM simulation of fluid-solid flow of a tapered column separation bed

Research on recycling waste Printed Circuit Boards(PCB) is at the forefront of preventing environmental pollution and finding ways to recycle resources.The Tapered Column Separation Bed(TCSB) is invented aiming at disposing the problem that fine particles of waste printed circuit boards cannot be separated efficiently so as to obtain further insight about the underlying mechanisms and demonstrate the separation feasibility in the tapered column separation bed.In this work,a Computational Fluid Dynamics(CFD) coupled with Discrete Element Method(DEM) model for two-phase flow has been extended to simulate the fluid-solid flow in the tapered column separation bed.Its validity is demonstrated by its successful capturing the key features of particles' flow pattern,velocity,the pressure distribution,the axial position with time and axial force for particles with different densities.Simulation results show that the plastic particles and resin particles become overflow,while copper particles,iron particles and aluminum particles successively become underflow,with a discharge water flow rate of 1 m^3/h,an obliquity of 30°.The simulated results agree reasonably well with the experimental observation.Using this equipment to separate waste PCBs is feasible,theoretically.

Cavitation Rates in Water with Dissolved Gas and Other Impurities

Our objective is a better understanding of the role of physical properties of real fluids in the thermodynamics of cavitation in impure water. An extension to the classical homogenous nucleation theory suitable for mixtures is presented in attempt to address the discrepancy between the theoretical predictions and practical observations of cavitation rates in water at normal temperatures. The extension takes into account the non-equilibrium (dissipative) effects involved in nuclei formation through a substance dependent correction coefficient to be determined experimentally. The theory of thermodynamic fluctuations is applied to derive the work of formation of a bubble nucleus. The value of the correction coefficient is estimated using preliminary experimental data from a convergent-divergent nozzle. An application of the results to the numerical prediction of the cavitation zones in a radial-flow water pump is shown.