一种保持界面位置不动的水平集函数隐式重新初始化方法.pdfVIP

第28卷 第5期 计 算 物 理 V01．28．NO．5 2011年 9月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Sep．，2011 文章编号：1001-246X(2011)05-0667—10 一 种保持界面位置不动的水平集 函数 隐式重新初始化方法 张 彬 ， 刘小民 (西安交通大学能源与动力工程学 院，陕西 西安 710049) 摘 要 ：以保持界面不动为主要 目标 ，推导水平集函数重新初始化方程 中光滑参数的取值公式 ，得到一种新型的 隐式重新初始化方法．证明应用该方法进行重新初始化时，保证界面附近水平集函数节点值符号不变时间步长 只需满足原始的CFL条件即可．最后 ，将该方法与原有初始化方法结合 ，得到一种准确且快速 的隐式重新初始 化方法．用数值算例验证该方法 的有效性 ． 关键词 ：水平集 ；重新初始化 ；光滑参数 ；CFL条件 中图分类号 ：0351 文献标识码 ：A 0 引言 水平集方法是一种十分优秀的运动界面追踪方法，它最早 由Osher和 Sethian…于 1988年提出．目前 ，水 平集方法已在物理学、固体力学、流体力学-z 、晶体生长、燃烧、流场模拟、多相流、图形 图像处理和半导体 加工等众多领域得到了广泛的应用 ．文 [1，6]都曾经提到水平集方法的数值实现的一个关键 问题是 ：水 平集函数经过演化后 ，如何在保持界面位置不动的情况下，使水平集 函数重新变回符号距离函数．此问题也 就是本文要讨论的水平集函数 的重新初始化问题．重新初始化不是为了一定要使水平集 函数变成符号距离 函数 ，而是因为当水平集函数的梯度过大或过小时，水平集函数的数值求解都会有很大麻烦 ；而将水平集 函数重新初始化为符号距离函数是解决这种数值问题的一种很好的方法． Chopp 最早提出了将水平集函数重新初始化为符号距离函数 的想法 ，以克服水平集 函数演化过程 中 由于水平集函数的梯度过大或过小而引起的数值问题．此后水平集函数重新初始化的算法大量涌现，到 目前 为止，应用比较广泛的重新初始化算法大体可以分为 求解偏微分方程的隐式方法和基于快速行进法的显 式方法两大类．Sussman等 ¨成功地构造了一个用于水平集函数重新初始化 的偏微分方程 ，通过迭代求解该 偏微分方程可以隐式地实现符号距离函数 的重新初始化 ，此方法就是求解偏微分方程式的隐式方法的最初 形式．Sussman等所构造偏微分方程给出了偏微分方程重新初始化的一个十分重要的框架 ，之后很多人都在 此基础上发展了这种重新初始化方法．Peng等 改进了该偏微分方程 ，使其在应用于远离符号距离函数的 水平集 函数的重新初始化时能够得到更好 的效果．另外 ，上面所提到的显式重新初始化方法是一种基于快速 行进法思想的重新初始化方法．快速行进法 (FastMarchingMethod)的思想最早 由Sethian_l提出的，是一种 显式追踪界面的方法，但是要实现这种方法必须事先显式地知道界面的具体位置． 除了上面所提到的两种主要方法之外 ，Chan121和 Li 等提 出了一种无需重新初始化 的变分水平集方 法．这种方法通过在原水平集演化方程的一侧加入一项能量泛函和一些正则项 ，纠正了水平集函数对符号距 离函数的偏离，从而使水平集函数在演化的同时也保证了对符号距离函数的近似．之后，Duan等 ¨ 发展 了该变分水平集方法 ，并将其应用于 Stokes和 Navier．Stokes方程的拓扑优化研究 中，得到了满意的结果．虽 然变分水平集方法可以省去重新初始化过程 ，但是它也存在着一些 的缺陷．一方面，由于所添加 的能量泛函 是水平集函数的二阶导数项 ，所以在曲率较大处求解水平集演化方程时能量泛函的数值计算误差会很大 ，并 收稿 日期 ：2010—11—11；修 回日期 ：2011—02—17 基金项 目：中央高校基本科研业务费专项资金项 目(xj和 中国博士后科学基金 (20100481329)资助项 目 作者简介 ：张彬 (1986一)男，山东蒙阴，硕士生，主要从事水平集