论文部分内容阅读
Level Set方法是处理封闭运动界面随时间演化过程几何拓扑变化的有效工具,已经成为处理广义界面问题最重要的方法之一,广泛应用于各种与曲线演化相关的科学研究和工程计算领域,具有重要的理论研究意义和应用背景。随着科学研究的不断深入,以及工程计算问题提出的更高需求,有效地捕捉和追踪运动界面变得更为重要。多介质界面和图像轮廓界面是广义界面问题中最普遍的问题,本文采用理论分析和数值模拟相结合的方法,探究Level Set方法在重新初始化方面的改进及其在欧拉多介质界面计算问题中的应用;探究Level Set方法在变分理论方面的改进及其在图像轮廓界面捕捉问题中的应用。主要开展工作如下:(1)详细阐述Level Set方法的基本思想和基本理论,分别介绍Level Set方法在多介质界面和图像轮廓界面处理领域中的国内外研究现状及其发展趋势,归纳Level Set方法求解问题的一般步骤,分析Level Set方法的优势与不足。(2)针对重新初始化造成的界面位置偏离问题,提出一种固定界面Level Set方法,主要改进在于将传统Level Set方法中固定的光滑参数修正为变化的光滑参数,通过固定界面有效地改善了Level Set方法的质量守恒性,具有较高的分辨率。固定界面Level Set方法的基本思想是在重新初始化相邻两次迭代过程中,以固定界面为基本原则,构建保持界面不动满足的等式条件,从而推导光滑参数满足的计算公式。讨论和分析了重新初始化、常数流场、旋转流场和剪切流场数值算例,结果表明,与传统Level Set方法相比,固定界面Level Set方法改善了流体体积的守恒性,能够捕捉到尖锐、狭窄、细长的界面,且界面的几何变化没有影响到周围流场的速度变化。(3)讨论固定界面Level Set方法在计算欧拉多介质流体计算问题中的应用,有效地抑制了流体界面附近产生的数值振荡和耗散,提高了间断处的分辨率,改善了气泡的质量守恒性,具有较强地捕捉流体界面能力。简述相关基本理论、控制方程及其数值方法,归纳固定界面Level Set方法求解欧拉多介质流体问题的详细步骤。讨论和分析了气-气激波管、气-液激波管和激波作用气泡问题的数值模拟,结果表明,与传统Level Set方法相比,固定界面Level Set方法有效地消除了激波处的数值耗散,密度、等压和速度的数值解更加吻合精确解,能够更准确地捕捉激波、膨胀波和接触间断的位置,能够更精细地描述气泡在激波作用下的变形过程,并降低气泡质量的损失。(4)讨论固定界面Level Set方法在爆轰波阵面传播计算问题中的应用,有效地避免了由于误差累积导致的爆轰波阵面偏离现象,具有较高的计算精度和分辨率,能够较为准确地描述爆轰波阵面的传播情况。简述相关基本理论、控制方程及其数值方法,归纳固定界面Level Set方法求解爆轰波阵面传播问题的详细步骤。讨论和分析了中心一点起爆、一点起爆、两点起爆和圆弧炸药平面起爆问题的数值模拟,结果表明,与传统Level Set方法相比,固定界面Level Set方法能够更为准确地捕捉到合理地爆轰波阵面,数值计算结果精度更高。(5)针对传统变分Level Set方法处理图像轮廓界面出现的边界粘连问题,提出一种自适应变分Level Set方法,主要改进在于将外部能量泛函中的常数修正为依赖于图像信息的函数,自适应变分Level Set方法的演化速度快,数值计算稳定,有效性和适用性较强。新参数函数借助于图像的二阶导数值,自动确定零等值线演化的方向,自动调整零等值线的演化速度,能够确保演化曲线的长度和面积同时取得全局最小值。讨论和分析了人造图像、医学图像、多目标图像轮廓界面等数值算例,结果表明,与传统变分Level Set方法相比,自适应变分Level Set方法对演化曲线的初始位置不敏感,可以实现自适应运动,能够同时捕捉图像的内外轮廓界面、狭窄深度凹陷区域边界以及多目标物体的边界,特别适用于小距离间隔多目标物体边界轮廓及其弱边界的图像轮廓界面捕捉。