A SBP-SAT method for the numerical solution and its application in the wave motion
YANG Zailin1, 2,SUN Cheng1,JIANG Guanxixi1,YANG Yong1,2
1.College of Aerospace and Civil Engineering, Harbin Engineering University, Harbin 150080, China;
2.Key Laboratory of Advanced Material of Ship and Mechanics, Ministry of Indus-try and Information Technology, Harbin Engineering University, Harbin 150080, China
摘要基于分部求和(summation by parts)方法和同时逼近项(simultaneous approximation terms)技术建立的有限差分方法,具有更高的精度和稳定性。同时在介质几何不连续、参数突变条件具有较大的优势。国内对SBP-SAT方法的相关研究目前较少,论文对该方法的研究背景,方法发展过程进行了介绍并基于SBP-SAT方法和弹性波动理论,结合初边值条件,推导出曲线网格条件下的弹性波动SBP-SAT离散方程。最后,通过数值模拟实现地震波传播过程,介绍该方法在地震数值模拟领域中的应用价值和前景。
Abstract:Summation by parts with simultaneous approximation terms have a set of properties that leading to provable time stability and higher order spatial discretization of partial differential equations. The finite difference methods based on this theory can establish the scheme that have incomparable advantages for simulating geometric discontinuity or complex media. At present, there are few researches on SBP-SAT methods in China. Therefore, this paper provides the history of SBP-SAT methods and its application to the numerical solution of partial differential equations. Then, the discrete equations of elastic wave under the curvilinear coordinates are derived. At last, combining initial value conditions and boundary conditions, a simple simulation for imitating seismic wave propagation is functioning. The potential and the application value are present.
杨在林1,2,孙铖1,蒋关希曦1,杨勇1,2. SBP-SAT方法及其在波动领域的应用[J]. 振动与冲击, 2020, 39(12): 150-157.
YANG Zailin1, 2,SUN Cheng1,JIANG Guanxixi1,YANG Yong1,2. A SBP-SAT method for the numerical solution and its application in the wave motion. JOURNAL OF VIBRATION AND SHOCK, 2020, 39(12): 150-157.
[1] 牟永光, 裴正林. 三维复杂介质地震数值模拟[M]. 石油工业出版社. 2005.
MU Yongguang, PEI Zhenglin. Numerical simulation of seismic wave propagation[M]. Progress in Geophysics, 19(4):933-941. . 2004
[2] 廖振鹏. 工程波动理论导论[M]. 第2版. 北京:科学出版社,2002.136-285.
LIAO Zhenpeng. Introduction to Wave Motion Theories in Engineering(2nd edn) [M].Beijing Science Press. 2002
[3] Alterman Z, Karal F C Jr. Propagation of elastic waves in layered media by finite difference methods[J]. Bulletin of the Seismological Society of America, 1968,58(1):367-398.
[4] 王勖成, 邵敏编著. 有限元法基本原理和数值方法. 第二版[M]. 清华大学出版社. 1999.
WANG Xucheng, Shao min. Basic principle and numerical method of finite element method(2nd edn) [M]. Tsinghua University Press .1999
[5] Berg P, If F, Skovgaard O. A spectral method for seismic wave propagation in elastic media[J]. Wave Motion, 1990, 12(5):415-427.
[6] Michel Bouchon, Michel Campillo, Stephane Gaffet.. A boundary integral equation-discrete Wavenumber reprentation method to study wave propagation in multilayered media having irregular interfaces[J]. Geophysics. 1989, 54(9): 1134-1140.
[7] 杨在林, 孙铖, 王耀,等. 含方形凹陷半无限非均匀介质波动问题 FDM 模拟[J]. 地震工程学报, 2015, 37(2):553-558.
YANG Zailin, Sun Cheng, Wang Yao, et al. 2015. Simulation of wave motion in a semi-infinite inhomogeneous medium with a rectangular depression[J], China earthquake engineering journal, 37(2):553-558.
[8] Yang, Zailin, Wang, et al. Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method[J]. Earthquake Engineering and Engineering Vibration, 2013, 12(4):569-576.
[9] Kreiss H O, Scherer G. Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations[J]. Mathematical Aspects of Finite Elements in Partial Differential Equations, 1974:195-212.
[10] Boom P D, Zingg D W. A generalized framework for nodal first derivative summation-by-parts operators[J]. Journal of Computational Physics, 2014, 266(6):214-239.
[11] Hicken, J. E, Zingg, et al. Summation-by-parts operators and high-order quadrature[J]. Journal of Computational & Applied Mathematics, 2013, 237(1):111-125.
[12] Gustafsson B. The convergence rate for difference approximations to mixed initial boundary value problems[J]. Mathematics of Computation, 1975, 29(130):396-406.
[13] Diener P, Dorband E N, Schnetter E, et al. Optimized High-Order Derivative and Dissipation Operators Satisfying Summation by Parts, and Applications in Three-dimensional Multi-block Evolutions[J]. Journal of Scientific Computing, 2007, 32(1):109-145.
[14] Kreiss H O, Scherer G. On the Existence of Energy Estimates for Difference Approximations for Hyperbolic Systems[J]. Department of Information Technology Uppsala University, 1977.
[15] Olsson P. High-order difference methods and data-parallel implementation[D]. Uppsala University, 1992.
[16] Godela Scherer. On energy estimates for difference approximations to hyperbolic partial differential equations[D]. Uppsala University, 1977.
[17] Bo S. Summation by Parts for Finite Difference Approximations for d/dx[J]. Journal of Computational Physics, 1994, 110(1):47-67.
[18] Carpenter M H, Nordström J, Gottlieb D. A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy[M]. Academic Press Professional, Inc. 1999.
[19] Mattsson K. Summation by parts operators for finite difference approximations of second derivatives. Academic Press Professional, Inc. 2004.
[20] Svärd M, Nordström J. On the order of accuracy for difference approximations of initial-boundary value problems[J]. Journal of Computational Physics, 2006, 218(1):333-352.
[21] Mattsson K, Svärd M, Shoeybi M. Stable and accurate schemes for the compressible Navier–Stokes equations[J]. Journal of Computational Physics, 2008, 227(4):2293-2316.
[22] Kamakoti R, Pantano C. High-Order Narrow Stencil Finite-Difference Approximations of Second-Order Derivatives Involving Variable Coefficients[J]. 2009, 31(6):4222-4243.
[23] Mattsson K. Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients[J]. Journal of Scientific Computing, 51(3):650-682, 2012.
[24] Mattsson K, Svärd M, Nordström, et al. Stable and accurate artificial dissipation[J]. Journal of Scientific Computing, 2004, 21(1):57-79.
[25] Nordström, Jan M. Conservative Finite Difference Formulations, Variable Coefficients, Energy Estimates and Artificial Dissipation[J]. Journal of Scientific Computing, 29(3):375-404, 2006.
[26] Fisher T C, Carpenter M H. High-order entropy stable finite difference schemes for nonlinear conservation laws: Finite domains[J]. Journal of Computational Physics, 2013, 252(C):518-557.
[27] Fisher T C, Carpenter M H, Nordström J, et al. Discretely conservative finite-difference formulations for nonlinear conservation laws in split form: Theory and boundary conditions[J]. Journal of Computational Physics, 2013, 234(1):353-375.
[28] Olsson P, Oliger J. Energy and maximum norm estimates for nonlinear conservation laws[J]. 1994.
[29] Sandham N D, Li Q, Yee H C. Entropy Splitting for High-Order Numerical Simulation of Compressible Turbulence[J]. Journal of Computational Physics, 2002, 178(2):307-322.
[30] Yee H C, Sjögreen B. Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations[J]. Fluid Mechanics & Its Applications, 2004, 66:141-198.
[31] Yee H C, Vinokur M, Djomehri M J. Entropy splitting and numerical dissipation[J]. Journal of Computational Physics, 2000, 162(1):33-81.
[32] M. Svärd. On Coordinate Transformations for Summation-by-Parts Operators[J]. Journal of Scientific Computing, 2004, 20(1):29-42.
[33] Hicken J E, Zingg D W. Superconvergent Functional Estimates from Summation-By-Parts Finite-Difference Discretizations[J]. Society for Industrial and Applied Mathematics, 33(2):893-922, 2011.
[34] Hicken, J. E, Zingg, et al. Summation-by-parts operators and high-order quadrature[J]. Journal of Computational & Applied Mathematics, 2013, 237(1):111-125.
[35] Fernández D C D R, Hicken J E, Zingg D W. Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations[J]. Computers & Fluids, 2014, 95(95):171-196.
[36] Fernández D C D R, Boom P D, Zingg D W. A generalized framework for nodal first derivative summation-by-parts operators[J]. Journal of Computational Physics, 2014, 266(6):214-239.
[37] Ranocha H. Generalised summation-by-parts operators and variable coefficients[J]. Journal of Computational Physics, 2018, 362:20-48.
[38] Hicken J E, Fernández D C D R, Zingg D W. Multidimensional Summation-By-Parts Operators: General Theory and Application to Simplex Elements[J]. Mathematics, 2015, 38(4):A1935-A1958.
[39] Fernández D C D R, Hicken J E, Zingg D W. Simultaneous Approximation Terms for Multi-dimensional Summation-by-Parts Operators[J]. Journal of Scientific Computing, 2016:1-28.
[40] Funaro D. Domain decomposition methods for pseudo spectral approximations[J]. Numerische Mathematik, 1987, 52(3):329-344.
[41] Funaro D, Gottlieb D. A New Method of Imposing Boundary Conditions in Pseudospectral Approximations of Hyperbolic Equations[J]. Mathematics of Computation, 1988, 51(184):599-613.
[42] Carpenter M H, Gottlieb D, Abarbanel S. Time-Stable Boundary Conditions for Finite-Difference Schemes Solving Hyperbolic Systems: Methodology and Application to High-Order Compact Schemes[J]. Journal of Computational Physics, 1994, 111(2):220-236.
[43] Carpenter M H, Nordström J, Gottlieb D. Revisiting and Extending Interface Penalties for Multi-domain Summation-by-Parts Operators[J]. Journal of Scientific Computing, 2010, 45(1-3):118-150.
[44] Jan Nordström, Carpenterb M H. Boundary and Interface Conditions for High-Order Finite-Difference Methods Applied to the Euler and Navier-Stokes Equations[M]. Academic Press Professional, Inc. 1999.
[45] Jan Nordström, Carpenterb M H. High-Order Finite Difference Methods, Multidimensional Linear Problems, and Curvilinear Coordinates[J]. Journal of Computational Physics, 2001, 173(1):149-174.
[46] Mattsson K. Boundary Procedures for Summation-by-Parts Operators[J]. Journal of Scientific Computing, 2003, 18(1):133-153.
[47] Svärd M, Carpenter M H, Nordström J. A stable high-order finite difference scheme for the compressible Navier–Stokes equations, far-field boundary conditions[J]. Journal of Computational Physics, 2007, 227(10):4805-4824.
[48] Svärd M, Nordström J. A stable high-order finite difference scheme for the compressible Navier–Stokes equations : No-slip wall boundary conditions[J]. Journal of Computational Physics, 227(10):4805-4824, 2008.
[49] Nordström J, Gong J, Weide E V D. A stable and conservative high order multi-block method for the compressible Navier-Stokes equations[J]. Journal of Computational Physics, 2009, 228(24):9020-9035.
[50] Berg J, Nordström J. Stable Robin solid wall boundary conditions for the Navier-Stokes equations[J]. Journal of Computational Physics, 2011, 230(19):7519-7532.
[51] Mattsson K, Ham F, Iaccarino G. Stable Boundary Treatment for the Wave Equation on Second-Order Form[J]. Journal of Scientific Computing, 2009, 41(3):366-383.
[52] Duru K, Kozdon J E, Kreiss G. Boundary conditions and stability of a perfectly matched layer for the elastic wave equation in first order form[J]. Journal of Computational Physics, 2015, 303(C):372-395.
[53] Mattsson K, Carpenter M H. Stable and Accurate Interpolation Operators for High-Order Multi-Block Finite-Difference Methods[J]. Siam J.sci.comput, 2010, 32(4):2298-2320.
[54] Kozdon J E, Wilcox L C. Stable Coupling of Nonconforming, High-Order Finite Difference Methods[J]. Siam Journal on Scientific Computing, 2015, 38(2):A923-A952.
[55] Friedrich L, Fernández D C D R, Winters A R, et al. Conservative and Stable Degree Preserving SBP Operators for Non-conforming Meshes[J]. Journal of Scientific Computing, 2017:1-30.
[56] Wang S, Virta K, Kreiss G. High Order Finite Difference Methods for the Wave Equation with Non-conforming Grid Interfaces[J]. Journal of Scientific Computing, 2016.
[57] Wang S. An Improved High Order Finite Difference Method for Non-conforming Grid Interfaces for the Wave Equation[J]. Journal of Scientific Computing, 2017:1-18.
[58] Kreiss H O, Petersson N A. Difference Approximations for the Second Order Wave Equation[J]. Siam Journal on Numerical Analysis, 2003, 40(5):1940-1967.
[59] Kreiss H O, Petersson N A. A Second Order Accurate Embedded Boundary Method for the Wave Equation with Dirichlet Data.[J]. Siam Journal on Scientific Computing, 2006, 27(27):1141-1167.
[60] Appelo D, Petersson N A. A Stable Finite Difference Method for the Elastic Wave Equation on Complex Geometries with Free Surfaces[J]. Communications in Computational Physics, 2009, 5(1):84-107.
[61] Sjögreen B. A Fourth Order Accurate Finite Difference Scheme for the Elastic Wave Equation in Second Order Formulation[J]. Journal of Scientific Computing, 2012, 52(1):17-48.
[62] Sjögreen B, Petersson N A. Source Estimation by Full Wave Form Inversion[J]. Journal of Scientific Computing, 2014, 59(1):247-276.
[63] Anders P N, Sjögreen B. Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains[J]. Communications in Computational Physics, 2014, 16(4):913-955.
[64] Petersson N A, Sjögreen B. Wave propagation in anisotropic elastic materials and curvilinear coordinates using a summation-by-parts finite difference method[J]. Journal of Computational Physics, 2015, 299:820-841.
[65] Petersson N A, Sjögreen B. High Order Accurate Finite Difference Modeling of Seismo-Acoustic Wave Propagation in a Moving Atmosphere and a Heterogeneous Earth Model Coupled Across a Realistic Topography[J]. Journal of Scientific Computing, 2017(4):1-34.
[66] Dovgilovich L, Sofronov I. High-accuracy finite-difference schemes for solving elastodynamic problems in curvilinear coordinates within multiblock approach[J]. Applied Numerical Mathematics, 2015, 93(C):176-194.
[67] Kozdon J E, Nordström J. Interaction of Waves with Frictional Interfaces Using Summation-by-Parts Difference Operators: Weak Enforcement of Nonlinear Boundary Conditions[J]. Journal of Scientific Computing, 2012, 50(2):341-367.
[68] Kozdon J E, Nordström J. Simulation of Dynamic Earthquake Ruptures in Complex Geometries Using High-Order Finite Difference Methods[J]. Journal of Scientific Computing, 2013, 55(1):92-124.
[69] Duru K, Dunham E M. Dynamic earthquake rupture simulations on nonplanar faults embedded in 3D geometrically complex, heterogeneous elastic solids[J]. Journal of Computational Physics, 2016, 305:185-207.
[70] O'Reilly O, Nordström J, Kozdon J E, et al. Simulation of Earthquake Rupture Dynamics in Complex Geometries Using Coupled Finite Difference and Finite Volume Methods[J]. Communications in Computational Physics, 2015, 17(2):337-370.
[71] Erickson B A, Dunham E M, Khosravifar A. A finite difference method for off-fault plasticity throughout the earthquake cycle[J]. Journal of Mechanics Physics of Solids, 2017, 109.50-77.
[72] Torberntsson K, Stiernström V, Mattsson K, et al. A finite difference method for earthquake sequences in poroelastic solids[J]. Computational Geosciences, 2018, 22(5):1351-1370.
[73] Gao L, Fernandez D C D R, Carpenter M, et al. SBP-SAT finite difference discretization of acoustic wave equations on staggered block-wise uniform grids[J]. 2018.