迎风紧致格式求解Hamilton—Jacobi方程.pdf

维普资讯 计 算 物 理 第 22卷 第 2期 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Vo1．22．No．2 2005年 3月 Mal．，2005 [文章编号] 1001246X(2005)02．O117．06 迎风紧致格式求解 Hamilton—Jacobi方程 田保林， 傅德薰， 马延文， 李新亮 (中国科学院力学研究所 LNM 国家重点实验室，北京 100080) [摘 要] 基于Hamilton．Jacobi(H．J)方程和双曲型守恒律之间的关系，将三阶和五阶迎风紧致格式推广应用 于求解H．J方程，建立了高精度的H．J方程求解方法．给出了一维和二维典型数值算例的计算结果，其中包括一 个平面激波作用下的Richtmyer-Meshkov界面不稳定性 问题 ．数值试验表明，在解的光滑区域该方法具有高精度 ， 而在导数不连续的不光滑区域也获得了比较好的分辨效果 ．相 比于同阶精度的WENO格式 ，本方法具有更小的 数值耗散，从而有利于多尺度复杂流动的模拟 中H．J方程的求解 ． [关键词] Hamilton．Jacobi方程；迎风紧致格式 [中图分类号] O35 [文献标识码] A O 引言 Hamilton．Jacobi(H．J)方程在许多实际问题 中都有应用，如常用的界面追踪技术——kve1Set方法 中所用 的控制方程就属于H．J方程n，在优化控制理论、几何光学和图像处理等方面也有广泛应用 ． n维空间下一般形式的H．J方程为 f≯。+ ( ，… ， ，t，≯ ， … ， ≯ )=0， ． {【 (1) ≯(，0)= ≯()． 方程 (1)的解通常是连续的，但是即使给定充分光滑的初值 ，它的解也可能是不光滑的，而双 曲型方程 的 解则可能是不连续的．在满足熵条件的情况下，方程 (1)具有唯一解 ．关于方程 (1)解的存在性和唯一性 问题 ， Crandall_2等做了专 门的论述 ． H．J方程和双曲型守恒律有着紧密的联系，在一维情形下H．J方程可以由相应的双 曲型守恒律方程积分 获得 ．对于一维 H．J方程 声。+H(声)=0， (2) 如果令 “= ≯，那么方程 (2)就可以转换成双曲型守恒律 “。+H(“) =0． (3) 对多维情形虽没有相应的转化关系，但也可以从形式上认为H．J方程是双曲型守恒律 的积分形式 ． 对于双曲型守恒律方程的求解 ，已经发展 了很多成熟的迎风型通量分裂方法，诸如全局和局部的 lax． Friedrichs格式 (LF，LLF)、Godunov型格式、具有熵修正的 Roe格式 (RF)等 ．基于 H．J方程和双 曲型方程的联 系，可 以将这些方法应用到 H．J方程的求解 中．文 [3，4]以这些单调格式为基架，以高精度的 ENO格式和 WENO格式求解 H．J方程 ．文 [7～10]从粘性解 的角度分析了双 曲型守恒律与H．J方程的关系，并将其发展的 求解双曲型守恒律方程的二阶精度 的松弛格式 (relaxationschem