Efficient Finite Difference WENO Scheme for Hyperbolic Systems withNon-conservativeProducts  

作　　者：Dinshaw S.Balsara Deepak Bhoriya Chi-Wang Shu Harish Kumar 

机构地区：[1]Physics Department,University of Notre Dame,Notre Dame,USA [2]ACMS Department,University of Notre Dame,Notre Dame,USA [3]Division of Applied Mathematics,Brown University,Providence,USA [4]Department of Mathematics,Indian Institute of Technology,Delhi,India

出　　处：《Communications on Applied Mathematics and Computation》2024年第2期907-962,共56页应用数学与计算数学学报（英文）

基　　金：support via NSF grants NSF-19-04774,NSF-AST-2009776,NASA-2020-1241;NASA grant 80NSSC22K0628.DSB;HK acknowledge support from a Vajra award,VJR/2018/00129;a travel grant from Notre Dame International;support via AFOSR grant FA9550-20-1-0055;NSF grant DMS-2010107.

摘　　要：Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we hav

关 键 词：Hyperbolic PDEs Numerical schemes Non-conservative products Stiff source terms Finite difference WENO 

分 类 号：O24[理学—计算数学] O29[理学—数学]