求解水平集延拓方程的半拉格朗日方法.pdfVIP

Candidate Liu Yonghui Supervisor Associate Professor Xu Jianjun College Mathematics and Computational Science Program Computational Mathematics Specialization Computational Fluid Dynamics Degree Master of Science University Xiangtan University Date April 15, 2014 万方数据 湘潭大学 学位论文原创性声明 本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研 究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文 不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研 究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完 全意识到本声明的法律后果由本人承担。 作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，同 意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允 许论文被查阅和借阅。本人授权湘潭大学可以将本学位论文的全部或 部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复 制手段保存和汇编本学位论文。 涉密论文按学校规定处理。 作者签名： 日期： 年 月 日 导师签名： 日期： 年 月 日 万方数据 摘要 本文中，我们用半拉格朗日（S-L）方法求解一个可将界面上的量延拓到包含 界面的一个邻域中的Hamilton-Jacobian方程。水平集函数被用来表示界面。其中 的一个S-L方法是一阶方法，另一个是形式上的二阶方法。一阶方法采用一阶欧 拉方法来寻找出发点，双线性插值被用来近似出发点处的延拓量；形式上的二阶 方法采用二阶Runge-Kutta格式来寻找出发点，三阶ENO格式被用来进行插值。 最后我们给出了数值算例，很明显一阶S-L方法比二阶S-L方法效果要好，这是由 于穿过界面法线速度的间断不连续性所造成的。作为比较，文中也给出了经典的 高阶欧拉方法的数值计算。 关键词：水平集函数; 半拉格朗日方法; 插值; 高阶欧氏方法 I 万方数据 In this paper, we study semi-Lagrangian(S-L) methods for solving a Hamilton- Jacobian equation, which extends an interface quantity into a neighborhood of the interface. A l