LevelSet方法模拟界面运动

来源 :中国科学院工程热物理研究所 | 被引量 : 17次 | 上传用户:lconan
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
模拟界面运动的数值方法主要分为两种类型:界面追踪(Interface Tracking)和界面捕捉(Interface Capturing)。这两种方法的主要区别是,界面追踪为Lagrange型方法,界面捕捉为Euler型方法。VOF方法和LevelSet方法是界面捕捉中最经典的两种方法,本文主要研究LevelSet方法在界面捕捉中的应用。本文对VOF方法和LevelSet方法进行了详细的介绍,采用Taylor-Galerkin方法推导出LevelSet控制方程的有限元离散格式。使用C++语言编写了TG有限元方法和WENO有限差分方法求解LevelSet控制方程的代码。采用TG有限元LevelSet方法、WENO有限差分LevelSet方法和Youngs-VOF方法分别对常数、旋转、剪切流场中的界面运动算例进行数值模拟,经对比可知,TG-LevelSet方法和WENO-LevelSet方法各有优缺点,TG-LevelSet方法对剪切流场中界面模拟的还原性更好一些,WENO-LevelSet方法对旋转流场中缺口圆模拟的稳定性更好一些,这两种方法的求解结果均好于Youngs-VOF方法。许多物理现象(如火焰面、激波面等)在界面处会发生物理量间断,正确处理这些界面处物理量的间断,对数值模拟的结果存在至关重要的影响。直接求解界面附近物理量导数值,会造成计算结果发散或间断抹平、界面模糊化等虚假计算结果,鉴于以上原因,本文采用Ghost技术,使界面附近物理量导数的求解合理化,从而获得精确的间断面捕捉。结合Rankin-Hugoniot间断条件和火焰特性,给出适用于火焰界面的间断条件。SMAC方法是使用原始变量求解非稳态Navier-Stokes方程的MAC类方法中较成功的方法,通过计算经典方腔驱动算例,验证作者编写的SMAC代码是正确的。采用LevelSet方法模拟界面运动变化,SMAC方法求解流体力学控制方程,Ghost技术处理界面间断,对一维、二维火焰界面合并问题进行数值模拟,得到较好结果,从而验证使用LevelSet方法处理该类问题是适用的。论文较为详细的介绍了处理表面张力的常用模型(CSF模型)和以界面压力间断来处理表面张力的Ghost模型。Fluent商业软件采用结合CSF模型的VOF方法模拟考虑表面张力的界面运动,本文采用结合Ghost技术的LevelSet方法。通过使用上述两种方法对液滴上升合并问题的模拟,可以看出LwFlow程序和Fluent商业软件得到的计算结果非常相近,进而说明LwFlow程序中,采用LevelSet方法、SMAC方法和Ghost技术相结合的思路求解液滴合并问题是合适的,同时通过相对质量变化率的比较,可以看出在此问题中LevelSet方法的守恒性要好于VOF方法。最后对LwFlow程序中涉及的格式、算法等进行了简要介绍。
其他文献
矿物燃料燃烧排放的CO2引起的全球气候变暖问题在逐渐威胁着人类的生存环境,进行CO2的减排工作迫在眉睫。Oxy-steam燃烧作为一种新型的燃烧方式,能够克服现有的以O2/CO2为代表
日常生活或临床手术中,金属异物可能进入或滞留体内,对人体产生一定的危害,多数情况下必需取出。体内金属异物的准确定位是其快速取出的首要条件,现有的定位方法中金属探测仪是一种简单便携又无损害的定位方式。此外,金属探测仪可以检测隐藏的金属物体,广泛应用于国防、公安、海关的安全检查以及建筑施工中。但是目前大多数金属探测仪的抗干扰能力差,测量灵敏度较低,不能对探测目标定位。针对这些问题,本文利用各种微弱信号
X射线齿科全景数字化成像系统是当今放射医疗科学研究应用的热点。目前,市场上存在直接数字化全景成像仪器,但是主要依靠进口,因此价格昂贵,应用范围小。中、小型医院目前仍普遍
为了适应医学影像检查技术的飞速发展,开展以计算机辅助诊断或智能化诊断为目标的医学图像处理和分析研究,已经成为目前这个领域的一个研究热点和将来发展的主要趋势。在颅脑CT