We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniq...We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that,taking into account the boundary conditions,fills the ghost cells with appropriate values,so that a standard reconstruction can be applied also in the boundary cells.In Naumann et al.(Appl.Math.Comput.325:252–270.https://doi.org/10.1016/j.amc.2017.12.041,2018),motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network,a different technique was explored that avoids the use of ghost cells,but instead employs for the boundary cells a different stencil,biased towards the interior of the domain.In this paper,extending that approach,which does not make use of ghost cells,we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids.In several numerical tests,we compare the novel reconstruction with the standard approach using ghost cells.展开更多
We present a new conservative semi-Lagrangian finite difference weighted essentially non-oscillatory scheme with adaptive order.This is an extension of the conservative semi-Lagrangian(SL)finite difference WENO scheme...We present a new conservative semi-Lagrangian finite difference weighted essentially non-oscillatory scheme with adaptive order.This is an extension of the conservative semi-Lagrangian(SL)finite difference WENO scheme in[Qiu and Shu,JCP,230(4)(2011),pp.863-889],in which linear weights in SL WENO framework were shown to not exist for variable coefficient problems.Hence,the order of accuracy is not optimal from reconstruction stencils.In this paper,we incorporate a recent WENO adaptive order(AO)technique[Balsara et al.,JCP,326(2016),pp.780-804]to the SL WENO framework.The new scheme can achieve an optimal high order of accuracy,while maintaining the properties of mass conservation and non-oscillatory capture of solutions from the original SL WENO.The positivity-preserving limiter is further applied to ensure the positivity of solutions.Finally,the scheme is applied to high dimensional problems by a fourth-order dimensional splitting.We demonstrate the effectiveness of the new scheme by extensive numerical tests on linear advection equations,the Vlasov-Poisson system,the guiding center Vlasov model as well as the incompressible Euler equations.展开更多
We propose a high order finite difference linear scheme combined with ahigh order bound preserving maximum-principle-preserving (MPP) flux limiter tosolve the incompressible flow system. For such problem with highly o...We propose a high order finite difference linear scheme combined with ahigh order bound preserving maximum-principle-preserving (MPP) flux limiter tosolve the incompressible flow system. For such problem with highly oscillatory structure but not strong shocks, our approach seems to be less dissipative and much lesscostly than a WENO type scheme, and has high resolution due to a Hermite reconstruction. Spurious numerical oscillations can be controlled by the weak MPP fluxlimiter. Numerical tests are performed for the Vlasov-Poisson system, the 2D guidingcenter model and the incompressible Euler system. The comparison between the linearand WENO type schemes, with and without the MPP flux limiter, will demonstrate thegood performance of our proposed approach.展开更多
In this article we present a new class of high order accurate Arbitrary-Eulerian-Lagrangian(ALE)one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimension...In this article we present a new class of high order accurate Arbitrary-Eulerian-Lagrangian(ALE)one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes.A WENO reconstruction algorithm is used to achieve high order accuracy in space and a high order one-step time discretization is achieved by using the local space-time Galerkin predictor proposed in[25].For that purpose,a new element-local weak formulation of the governing PDE is adopted on moving space-time elements.The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes.Moreover,a polynomial mapping defined by the same local space-time basis functions as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element.To maintain algorithmic simplicity,the final ALE one-step finite volume scheme uses moving triangular meshes with straight edges.This is possible in the ALE framework,which allows a local mesh velocity that is different from the local fluid velocity.We present numerical convergence rates for the schemes presented in this paper up to sixth order of accuracy in space and time and show some classical numerical test problems for the two-dimensional Euler equations of compressible gas dynamics.展开更多
基金MIUR-PRIN project 2017KKJP4X“Innovative numerical methods for evolutionary partial differential equations and applications”.Gabriella Puppo acknowledges also the support of 2019 Ateneo Sapienza research project no.RM11916B51CD40E1.
文摘We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that,taking into account the boundary conditions,fills the ghost cells with appropriate values,so that a standard reconstruction can be applied also in the boundary cells.In Naumann et al.(Appl.Math.Comput.325:252–270.https://doi.org/10.1016/j.amc.2017.12.041,2018),motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network,a different technique was explored that avoids the use of ghost cells,but instead employs for the boundary cells a different stencil,biased towards the interior of the domain.In this paper,extending that approach,which does not make use of ghost cells,we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids.In several numerical tests,we compare the novel reconstruction with the standard approach using ghost cells.
文摘We present a new conservative semi-Lagrangian finite difference weighted essentially non-oscillatory scheme with adaptive order.This is an extension of the conservative semi-Lagrangian(SL)finite difference WENO scheme in[Qiu and Shu,JCP,230(4)(2011),pp.863-889],in which linear weights in SL WENO framework were shown to not exist for variable coefficient problems.Hence,the order of accuracy is not optimal from reconstruction stencils.In this paper,we incorporate a recent WENO adaptive order(AO)technique[Balsara et al.,JCP,326(2016),pp.780-804]to the SL WENO framework.The new scheme can achieve an optimal high order of accuracy,while maintaining the properties of mass conservation and non-oscillatory capture of solutions from the original SL WENO.The positivity-preserving limiter is further applied to ensure the positivity of solutions.Finally,the scheme is applied to high dimensional problems by a fourth-order dimensional splitting.We demonstrate the effectiveness of the new scheme by extensive numerical tests on linear advection equations,the Vlasov-Poisson system,the guiding center Vlasov model as well as the incompressible Euler equations.
基金T.Xiong acknowledges support by NSFC grant No.11971025NSF grant of Fujian Province No.2019J06002the Strategic Priority Research Program of Chinese Academy of Sciences Grant No.XDA25010401.
文摘We propose a high order finite difference linear scheme combined with ahigh order bound preserving maximum-principle-preserving (MPP) flux limiter tosolve the incompressible flow system. For such problem with highly oscillatory structure but not strong shocks, our approach seems to be less dissipative and much lesscostly than a WENO type scheme, and has high resolution due to a Hermite reconstruction. Spurious numerical oscillations can be controlled by the weak MPP fluxlimiter. Numerical tests are performed for the Vlasov-Poisson system, the 2D guidingcenter model and the incompressible Euler system. The comparison between the linearand WENO type schemes, with and without the MPP flux limiter, will demonstrate thegood performance of our proposed approach.
基金the European Research Council(ERC)under the European Union’s Seventh Framework Programme(FP7/2007-2013)with the research project STiMulUs,ERC Grant agreement no.278267.
文摘In this article we present a new class of high order accurate Arbitrary-Eulerian-Lagrangian(ALE)one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes.A WENO reconstruction algorithm is used to achieve high order accuracy in space and a high order one-step time discretization is achieved by using the local space-time Galerkin predictor proposed in[25].For that purpose,a new element-local weak formulation of the governing PDE is adopted on moving space-time elements.The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes.Moreover,a polynomial mapping defined by the same local space-time basis functions as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element.To maintain algorithmic simplicity,the final ALE one-step finite volume scheme uses moving triangular meshes with straight edges.This is possible in the ALE framework,which allows a local mesh velocity that is different from the local fluid velocity.We present numerical convergence rates for the schemes presented in this paper up to sixth order of accuracy in space and time and show some classical numerical test problems for the two-dimensional Euler equations of compressible gas dynamics.