HIGH ORDER WEIGHT COMPACR NONLINEAR SCHEME FOR TIME-DEPENDENT HAMILTON-JACOBI EQUATIONS
-
摘要:Hamilton-Jacobi (HJ) 方程是一类重要的非线性偏微分方程, 在物理学、流体力学、图像处理、微分几何、金融数学、最优化控制理论等方面有着广泛的应用. 由于HJ方程的弱解存在但不唯一, 且解的导数可能出现间断, 导致其数值求解具有一定的难度. 本文提出了非稳态HJ方程的7阶精度加权紧致非线性格式 (WCNS). 该格式结合了Hamilton函数的Lax-Friedrichs型通量分裂方法和一阶空间导数左、右极限值的高阶精度混合节点和半节点型中心差分格式. 基于7点全局模板和4个4点子模板推导了半节点函数值的高阶线性逼近和4个低阶线性逼近, 以及全局模板和子模板的光滑度量指标. 为避免间断附近数值解产生非物理振荡以及提高格式稳定性, 采用WENO型非线性插值方法计算半节点函数值. 时间离散采用3阶TVD型Runge-Kutta方法. 通过理论分析验证了WCNS格式对于光滑解具有最佳的7阶精度. 为方便比较, 经典的7阶WENO格式也被推广用于求解HJ方程. 数值结果表明, 本文提出的WCNS格式能够很好地模拟HJ方程的精确解, 且在光滑区域能够达到7阶精度; 与经典的同阶WENO格式相比, WCNS格式在精度、收敛性和分辨率方面更优, 计算效率略高.
-
关键词:
- Hamilton-Jacobi方程/
- 7阶WCNS/
- WENO插值/
- 高分辨率/
- 间断捕捉能力
Abstract:Hamilton-Jacobi (HJ) equations are an important class of nonlinear partial differential equations. They are often used in various applications, such as physics, fluid mechanics, image processing, differential geometry, financial mathematics, optimal control theory, and so on. Because the weak solutions of the HJ equations exist but are not unique, and the spatial derivatives of the solutions may be discontinuous, numerical difficulties arise in numerical solutions of these equations. This paper presents a seventh-order weighted compact nonlinear scheme (WCNS) for the time-dependent HJ equations. This scheme is composed of the monotone Lax-Friedrichs flux splitting method for the Hamilton functions and the high-order hybrid cell-node and cell-edge central differencing for the left and right limits of first-order spatial derivatives in the numerical Hamilton functions. A high-order linear approximation scheme and four low-order linear approximation schemes for the unknowns at half nodes are derived based on a seven-point global stencil and four four-point sub-stencils, respectively. The smoothness indicators of the global stencil and four sub-stencils are also derived. In order to avoid non-physical oscillations of numerical solutions near the discontinuities and improve the numerical stability of the designed scheme, the WENO-type nonlinear interpolation technique is adopted to compute the unknowns at half nodes. The third-order TVD Runge-Kutta method is used for time discretization. The presented WCNS scheme is verified to have the optimal seventh order of accuracy for smooth solutions by theoretical analysis. For the sake of comparison, the classical seventh-order WENO scheme for solving hyperbolic conservation laws is also extended to solve the HJ equations. Numerical results show that the presented WCNS scheme can well simulate the exact solutions and can achieve seventh-order accuracy in smooth regions. Compared with the classical WENO scheme of the same order, the presented WCNS scheme has better accuracy, convergence and resolution, and its computational efficiency is slightly higher. -
表 1一维情形时两种格式的数值误差、精度阶和CPU运行时间
Table 1.Numerical errors, convergence rates and CPU time obtained with two schemes for the 1D case
Method $N$ ${L_1}$error Order ${L_\infty }$error Order CPU time WCNS 60 2.00 × 10−8 — 3.72 × 10−8 — 0.078125 80 2.67 × 10−9 7.00 4.90 × 10−9 7.04 0.218750 100 5.60 × 10−10 7.00 1.02 × 10−9 7.05 0.453125 120 1.56 × 10−10 7.00 2.81 × 10−10 7.05 0.796875 140 5.33 × 10−11 6.98 9.52 × 10−11 7.03 1.312500 WENO 60 7.73 × 10−8 — 5.94 × 10−7 — 0.078125 80 1.33 × 10−8 6.13 1.31 × 10−7 5.24 0.203125 100 3.36 × 10−9 6.15 4.09 × 10−8 5.23 0.453125 120 1.14 × 10−9 5.94 1.57 × 10−8 5.24 0.828125 140 4.56 × 10−10 5.93 6.97× 10−9 5.27 1.093750 表 2二维情形时两种格式的数值误差、精度阶和CPU运行时间
Table 2.Numerical errors, convergence rates and CPU time obtained with two schemes for the 2D case
Method $M \times N$ ${L_1}$error Order ${L_\infty }$error Order CPU time WCNS 60 × 60 2.95× 10−8 — 4.80× 10−8 — 1.578125 80 × 80 3.95× 10−9 6.99 6.42× 10−9 6.99 4.843750 100 × 100 8.28 × 10−10 7.00 1.34× 10−9 7.01 12.65625 120 × 120 2.31 × 10−10 7.00 3.74 × 10−10 7.02 27.85937 140 × 140 7.86 × 10−11 7.00 1.27 × 10−10 7.01 54.92187 WENO 60 × 60 2.96× 10−8 — 1.88 × 10−7 — 1.140625 80 × 80 4.48× 10−9 6.57 3.84× 10−8 5.52 3.765625 100 × 100 1.07× 10−9 6.41 1.11× 10−8 5.55 9.984375 120 × 120 3.33 × 10−10 6.41 3.94× 10−9 5.70 22.09375 140 × 140 1.24 × 10−10 6.41 1.80× 10−9 5.09 43.14062 -
[1] 刘瑛, 杜光勋, 全权等. 基于Hamilton-Jacobi方程的飞行器机动动作可达集分析. 自动化学报, 2016, 42(3): 347-357 (Liu Ying, Du Guangxun, Quan Quan, et al. Reachability calculation for aircraft maneuver using hamilton-jacobi function.Acta Automatica Sinica, 2016, 42(3): 347-357 (in Chinese)doi:10.16383/j.aas.2016.c140888 [2] Chan CK, Lau KS, Zhang BL. Simulation of a premixed turbulent flame with the discrete vortex method.International Journal for Numerical Methods in Engineering, 2000, 48(4): 613-627doi:10.1002/(SICI)1097-0207(20000610)48:4<613::AID-NME900>3.0.CO;2-7 [3] Huang L, Wong SC, Zhang M, et al. Revisiting Hughes’ dynamic continuum model for pedestrian flow and the development of an efficient solution algorithm.Transportation Research Part B:Methodological, 2009, 43(1): 127-141doi:10.1016/j.trb.2008.06.003 [4] Lefèvre V, Garnica A, Lopez-Pamies O. A WENO finite-difference scheme for a new class of Hamilton-Jacobi equations in nonlinear solid mechanics.Computer Methods in Applied Mechanics and Engineering, 2019, 349(1): 17-44 [5] 陈苏婷, 李霞. 非紧空间上折现Hamilton-Jacobi方程的黏性解. 华东师范大学学报, 2022, 2022(2): 9-15 (Chen Suting, Li Xia. The viscosity solution of the discounted Hamilton-Jacobi equation in non-compact space.Proceedings of the National Academy of Sciences, 2022, 2022(2): 9-15 (in Chinese) [6] Crandall MG, Lions PL. Viscosity solutions of Hamilton-Jacobi equations.Transactions of the American Mathematical Society, 1983, 277: 1-42doi:10.1090/S0002-9947-1983-0690039-8 [7] Acker F, Borges RBR, Costa B. An improved WENO-Z scheme.Journal of Computational Physics, 2016, 313: 726-753doi:10.1016/j.jcp.2016.01.038 [8] 骆信, 吴颂平. 改进的五阶WENO-Z+ 格式. 力学学报, 2019, 51(6): 1927-1939 (Luo Xin, Wu Songping. An improved fifth-order WENO-Z + scheme.Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(6): 1927-1939 (in Chinese)doi:10.6052/0459-1879-19-249 [9] Balsara DS, Shu CW. Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy.Journal of Computational Physics, 2000, 160(2): 405-452doi:10.1006/jcph.2000.6443 [10] Osher S, Sethian JA. Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations.Journal of Computational Physics, 1988, 79(1): 12-49doi:10.1016/0021-9991(88)90002-2 [11] Jiang GS, Peng D. Weighted ENO schemes for Hamilton-Jacobi equations.SIAM Journal on Scientific Computing, 2000, 21(6): 2126-2143doi:10.1137/S106482759732455X [12] Bryson S, Levy D. High-order central WENO schemes for multidimensional Hamilton-Jacobi equations.SIAM Journal on Numerical Analysis, 2003, 41(4): 1339-1369doi:10.1137/S0036142902408404 [13] Qiu J, Shu CW. Hermite WENO schemes for Hamilton–Jacobi equations.Journal of Computational Physics, 2005, 204(1): 82-99doi:10.1016/j.jcp.2004.10.003 [14] Ha Y, Kim CH, Yang H, et al. A sixth-order weighted essentially non-oscillatory schemes based on exponential polynomials for Hamilton−Jacobi equations.Journal of Scientific Computing, 2018, 75(3): 1675-1700doi:10.1007/s10915-017-0603-8 [15] Cheng X, Feng J. A sixth-order finite difference WENO scheme for Hamilton−Jacobi equations.International Journal of Computer Mathematics, 2019, 96(3): 568-584doi:10.1080/00207160.2018.1447665 [16] Zhu J, Zheng F, Qiu J. New Finite Difference Hermite WENO Schemes for Hamilton-Jacobi Equations.Journal of Scientific Computing, 2020, 83(1): 1-21doi:10.1007/s10915-020-01189-x [17] Rathan S. L1-type smoothness indicators based weighted essentially nonoscillatory scheme for Hamilton-Jacobi equations.International Journal for Numerical Methods in Fluids, 2020, 92(12): 1927-1947doi:10.1002/fld.4855 [18] Abedian R. A symmetrical WENO-Z scheme for solving Hamilton−Jacobi equations.International Journal of Modern Physics C, 2020, 31(3): 2050039doi:10.1142/S0129183120500394 [19] Kim CH, Ha Y, Yang H, et al. A third-order WENO scheme based on exponential polynomials for Hamilton-Jacobi equations.Applied Numerical Mathematics, 2021, 165: 167-183doi:10.1016/j.apnum.2021.01.020 [20] Deng X, Zhang H. Developing high-order weighted compact nonlinear schemes.Journal of Computational Physics, 2000, 165(1): 22-44doi:10.1006/jcph.2000.6594 [21] Deng X, Maekawa H. Compact high-order accurate nonlinear schemes.Journal of Computational Physics, 1997, 130(1): 77-91doi:10.1006/jcph.1996.5553 [22] Nonomura T, Fujii K. Effects of difference scheme type in high-order weighted compact nonlinear schemes.Journal of Computational Physics, 2009, 228(10): 3533-3539doi:10.1016/j.jcp.2009.02.018 [23] Zhang S, Jiang S, Shu CW. Development of nonlinear weighted compact schemes with increasingly higher order accuracy.Journal of Computational Physics, 2008, 227(15): 7294-7321doi:10.1016/j.jcp.2008.04.012 [24] Deng X, Jiang Y, Mao M, et al. A family of hybrid cell-edge and cell-node dissipative compact schemes satisfying geometric conservation law.Computers&Fluids, 2015, 116: 29-45 [25] Nonomura T, Fujii K. Robust explicit formulation of weighted compact nonlinear scheme.Computers&Fluids, 2013, 85: 8-18 [26] Wong ML, Lele SK. High-order localized dissipation weighted compact nonlinear scheme for shock-and interface-capturing in compressible flows.Journal of Computational Physics, 2017, 339: 179-209doi:10.1016/j.jcp.2017.03.008 [27] Zhao G, Sun M, Xie S, et al. Numerical dissipation control in an adaptive WCNS with a new smoothness indicator.Applied Mathematics and Computation, 2018, 330: 239-253doi:10.1016/j.amc.2018.01.019 [28] 洪正, 叶正寅. 各向同性湍流通过正激波的演化特征研究. 力学学报, 2018, 50(6): 1356-1367 (Hong Zheng, Ye Zhengyin. Study on evolution characteristics of isotropic turbulence passing through a normal shock wave.Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(6): 1356-1367 (in Chinese)doi:10.6052/0459-1879-18-129 [29] Jiang YQ, Zhou SG, Zhang X, et al. High-order weighted compact nonlinear scheme for one-and two-dimensional Hamilton-Jacobi equations.Applied Numerical Mathematics, 2022, 171: 353-368doi:10.1016/j.apnum.2021.09.012 [30] Hiejima T. A high-order weighted compact nonlinear scheme for compressible flows.Computers&Fluids, 2022, 232: 105199 [31] Jiang YQ, Zhou SG, Zhang X, et al. High order all-speed semi-implicit weighted compact nonlinear scheme for the isentropic Navier–Stokes equations.Journal of Computational and Applied Mathematics, 2022, 411: 114272doi:10.1016/j.cam.2022.114272 [32] Ma Y, Yan ZG, Liu H, et al. Improved weighted compact nonlinear scheme for implicit large-eddy simulations.Computers&Fluids, 2022, 240: 105412 [33] Gottlieb S, Shu CW, Tadmor E. Strong stability-preserving high-order time discretization methods.SIAM Review, 2001, 43(1): 89-112doi:10.1137/S003614450036757X