界面追踪的Level Set和Particle Level Set方法  

The Level Set and Particle Level Set method for tracing interface

作　　者：谷汉斌[1] 李炎保[1] 李绍武[1] 张庆河[1] 

机构地区：[1]天津大学建筑工程学院

出　　处：《水动力学研究与进展（A辑）》2005年第2期152-160,共9页Chinese Journal of Hydrodynamics

基　　金：国家自然科学基金项目(50279028)

摘　　要：Level Set方法是20世纪90年代很流行的界面追踪技术,国外已有个别实例将此方法应用于不可压缩流体的计算。但其对界面具有抹圆的作用,造成梯度很大的界面失真。Particle Level Set方法可以弥补这一缺陷。本文首先介绍了这两种界面追踪的方法。然后采用五阶WENO格式离散空间导数,并采用三阶 TVD Runge Kutta方法离散时间导数,用计算机Fortran语言建立了界面追踪的模型。最后通过平移流场中的同心方框、旋转流场中的缺口圆和剪切流场中圆的三维界面的追踪算例,对两种方法进行了对比验证。同时,用剪切流场中圆的算例对比了 Particle Level Set和单纯拉哥朗日方法。结果表明,Particle Level Set对界面的描述更为精细,可以用于大梯度的界面追踪;Particle Level Set和单纯拉哥朗日方法的结果很接近。Level Set method is the popular technique of interface tracing in 1990s, and it has been used in a few samples of incompressible flow. However, it smears out the sharp interface so that it would distort the large gradient. The Particle Level Set method can remedy this default. First, two methods are introduced in brief. Then a 5th-order WENO scheme for discrete spacial derivative, a third TVD-Runge-Kutta scheme for discrete temporal derivative, and the interface tracing model are built up in Fortran. At last two methods are compared by 3D interface tracing samples which are homocentric square frame in parallel shift flow, gap circle in whirly flow and torus in shearing flow. In the mean time, the Particle Level Set method is compared with Lagragian schemes for torus in shearing flow. It is shown that Particle Level Set method is suited to large gradient interface deformation and compares favorably with purely Lagragian schemes for interface resolution.

关 键 词：LevelSet ParticleLevelSet界面追踪 

分 类 号：TV131.4[水利工程—水力学及河流动力学]