LEVEL SET输运方程的求解方法及其对气-液两相流运动界面数值模拟的影响

用数值模拟方法来研究气-液两相流动与传热现象是当今多相流领域的一个热门课题.由于两相流固有的复杂性,气-液两相流界面迁移现象的数值模拟一直是两相流研究中的一大难点.本文介绍了捕捉气-液两相流相界面运动的水平集方法(Level Set)及其研究进展,介绍了求解Level Set输运方程的3种方法,即一般差分格式、Superbee-TVD格式和Runge-Kutta法-5阶WENO组合格式.结合主流场的求解,分别用这3种方法对4种典型相界面在5种流场中的迁移特性进行了模拟计算,并对计算结果进行了比较和分析.结果表明,Runge-Kutta法-5阶WENO组合格式求解Level Set输运方程的效果最好,在以后的计算中将主要采用这种组合格式来进行气-液相界面输运方程的求解.

气泡/液滴运动的Level Set方法模拟研究

利用Level Set方法，结合投影法求解了描述气泡/液滴运动的Navier-Stokes方程。对地面常重力场中不同大小的空气泡在高黏度糖浆溶液中的自由上升运动现象，数值模拟结果与实验观测结果符合甚好，表明该方法能够计算大密度比和黏度比$(>1000:1)$情况下的气液两相流动。而对等密度液滴的热毛细迁移现象的数值模拟结果同样能够与实验结果相一致，表明该方法同样适于研究具有Marangoni效应的两相流动现象，特别是在空间微重力环境中的气液两相传热现象中的局部流动与传热问题。

A Unified Tensor Level Set for Image Segmentation

This paper presents a new region-based unified tensor level set model for image segmentation. This model introduces a three-order tensor to comprehensively depict features of pixels, e.g., gray value and the local geometrical features, such as orientation and gradient, and then, by defining a weighted distance, we generalized the representative region-based level set method from scalar to tensor. The proposed model has four main advantages compared with the traditional representative method as follows. First, involving the Gaussian filter bank, the model is robust against noise, particularly the salt-and pepper-type noise. Second, considering the local geometrical features, e. g., orientation and gradient, the model pays more attention to boundaries and makes the evolving curve stop more easily at the boundary location. Third, due to the unified tensor pixel representation representing the pixels, the model segments images more accurately and naturally. Fourth, based on a weighted distance definition, the model possesses the capacity to cope with data varying from scalar to vector, then to high-order tensor. We apply the proposed method to synthetic, medical, and natural images, and the result suggests that the proposed method is superior to the available representative region-based level set method.

A Relay Level Set Method for Automatic Image Segmentation

This paper presents a new image segmentation method that applies an edge-based level set method in a relay fashion. The proposed method segments an image in a series of nested subregions that are automatically created by shrinking the stabilized curves in their previous subregions. The final result is obtained by combining all boundaries detected in these subregions. The proposed method has the following three advantages: 1) It can be automatically executed without human-computer interactions; 2) it applies the edge-based level set method with relay fashion to detect all boundaries; and 3) it automatically obtains a full segmentation without specifying the number of relays in advance. The comparison experiments illustrate that the proposed method performs better than the representative level set methods, and it can obtain similar or better results compared with other popular segmentation algorithms.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

T形微通道中互不相溶两相流数值模拟

采用摄动有限体积(PFV)算法和水平集(level set)技术对T形微通道内互不相溶两相流动进行了数值模拟研究。考察了两相界面张力和微通道壁面润湿性对流动的影响,精确地捕捉到了油水两相流动的界面。对一些典型的T形微通道油水两相流动进行了数值计算,模拟结果和实验结果吻合较好。分析总结出了微通道内两相流动过程中的一些基本规律,为微通道内的液液两相流动实验设计和工业应用提供了新的数值预测手段.

蒸汽-冷流体接触冷凝流动的数值模拟

介绍了关于蒸汽-冷流体直接接触冷凝流动与传热的数值计算模型与部分研究结果。用Level Set方法确定蒸汽-冷流体接触界面的位置和形状,建立了对蒸汽和冷流体普遍适用的动量、能量和质量守恒方程,在能量和质量寺恒方程中增加了部分项用于计算蒸汽冷凝所产生的影响。用有限差分法在交错网格上离散控制方程,用Runge-Kutta法-五阶WENO组合格式求解Level Set输运方程,用压力修正的迭代Projection方法求解动量方程,而用SIMPLE方法求解温度控制方程。对算例的计算结果表明,本文所建立的数值计算模型能反映物理现象的宏观特性。根据计算结果,分析了本文模型的优缺点,并指出了今后改进的方向。

气-液两相流界面迁移现象的数值模拟研究

借助于Level Set函数,建立了气-液两相的统一控制方程组,并在交错网格中进行离散.用两种格式,即Superbee-TVD格式和5阶WENO格式求解Level Set函数的输运方程,用SIMPLER算法的思想对主流场控制方程的求解方法进行改进.数值实验结果表明,在求解Level Set的控制方程时,5阶WENO方法比Superbee-TVD格式的结果更准确;用改进的数值算法可成功实现对密度比大于1 000/1的气-液两相流界面迁移问题的数值模拟.对几种典型大密度比气-液两相流问题的计算结果与实际问题的物理规律完全一致,验证了该方法的有效性和可靠性.

利用Ghost Fluid方法模拟激波与柱形界面相互作用

利用Ghost Fluid方法(后面简称Ghost方法)和γ-model方法，在同样的时空离散精度条件下，对激波与柱形界面相互作用的二维可压缩流场进行了直接模拟，并与实验结果相比较．从模拟结果看，在短时间内， Ghost方法和γ-model方法模拟的结果与实验结果基本相同，两种方法均正确地模拟出界面的位置、激波的强度和速度．但随着时间的发展，具有较大数值耗散的γ-model方法的计算结果与实验差别越来越大；而数值耗散较小的Ghost方法能较为正确地模拟界面的运动．

不同重力条件下不相溶流体层间界面波演化规律的数值模拟研究

报道了关于不相溶流体层间界面波演化规律的数值模拟研究及结果，重点考察了重力条件对界面波演化特性的影响。考虑在深度方向无限扩展的互不相容的两个流体层，上层流体比下层的轻，但比下层的运动速度快；两层流体间的界面上存在正弦波形的初始扰动，并随流体流动而不断变化。本文采用Level Set方法来实现对运动的相界面的追踪，用有限差分法来离散控制方程组。为了提高数值算法的稳定性，采用三阶的Runge-Kutta法来离散时间导数，而采用五阶的WENO（Weighted Essentially Non-oscillatory）格式来离散一阶对流输运项，并用压力修正投影法（Pressure Correction Projection Method）来实现离散控制方程组的求解。为了提高对复杂非稳态过程的解的准确度，采用了嵌套的三层迭代循环。本文对一系列工况条件下的界面波演化过程进行了计算；除了研究重力的作用之外，还考察了流体密度、粘性、表面张力、初始界面波频率、振幅及波数对界面波演化特性的影响。其中，上下流体层的最大密度比和粘性比可达3000/1，而重力加速度在0~5g0（g0=9.8m/s^2）之间变化，上下流体层间的最大速度差为8m/s。研究结果表明，随着重力、流体密度比、流体粘性比及表面张力的增加，界面波的演化受到不同程度的抑制，而界面波的传播速度也与重力及流体的密度、粘性和表面张力等因素相关。

Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

An Improved Ce/Se Scheme For Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

摄动有限体积算法及其在两相流中应用

采用一种全新的摄动有限体积(PFV)算法和水平集(Level Set)技术对液液两相系统中液滴坠落进行数值模拟,数值结果表明,PFV新算法具有节点少、精度高,效率高,编程方便等优点,能成功模拟液液两相流动,为两相流动数值模拟提供了一种新的途径.