The ADER approach to solve hyperbolic equations to very high order of accuracy has seen explosive developments in the last few years,including both methodological aspects as well as very ambitious applications.In spit...The ADER approach to solve hyperbolic equations to very high order of accuracy has seen explosive developments in the last few years,including both methodological aspects as well as very ambitious applications.In spite of methodological progress,the issues of efficiency and ease of implementation of the solution of the associated generalized Riemann problem(GRP)remain the centre of attention in the ADER approach.In the original formulation of ADER schemes,the proposed solution procedure for the GRP was based on(i)Taylor series expansion of the solution in time right at the element interface,(ii)subsequent application of the Cauchy-Kowalewskaya procedure to convert time derivatives to functionals of space derivatives,and(iii)solution of classical Riemann problems for high-order spatial derivatives to complete the Taylor series expansion.For realistic problems the Cauchy-Kowalewskaya procedure requires the use of symbolic manipulators and being rather cumbersome its replacement or simplification is highly desirable.In this paper we propose a new class of solvers for the GRP that avoid the Cauchy-Kowalewskaya procedure and result in simpler ADER schemes.This is achieved by exploiting the history of the numerical solution that makes it possible to devise a time-reconstruction procedure at the element interface.Still relying on a time Taylor series expansion of the solution at the interface,the time derivatives are then easily calculated from the time-reconstruction polynomial.The resulting schemes are called ADER-TR.A thorough study of the linear stability properties of the linear version of the schemes is carried out using the von Neumann method,thus deducing linear stability regions.Also,via careful numerical experiments,we deduce stability regions for the corresponding non-linear schemes.Numerical examples using the present simplified schemes of fifth and seventh order of accuracy in space and time show that these compare favourably with conventional ADER methods.This paper is restricted to the one-dimensional scalar case with source term,but preliminary results for the one-dimensional Euler equations indicate that the time-reconstruction approach offers significant advantages not only in terms of ease of implementation but also in terms of efficiency for the high-order range schemes.展开更多
The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is base...The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is based on the approach of Vila[25]and another is proposed in this article which mimics the classical 1D LaxWendroff scheme.In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods,demonstrated via the von Neumann stability analysis.Moreover,the issue of the consistency for the equations of gas dynamics is analyzed.An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates.This not only provides an improvement in accuracy of the numerical solutions,but also assures that the consistency condition on the gradient of the kernel function is satisfied using an equidistant distribution of particles in Lagrangian mass coordinates.Three different Riemann solvers are implemented for the first-order Godunov type SPH schemes in Lagrangian coordinates,namely the Godunov flux based on the exact Riemann solver,the Rusanov flux and a new modified Roe flux,following the work of Munz[17].Some well-known numerical 1D shock tube test cases[22]are solved,comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the first-order Godunov finite volume method in Eulerian coordinates and the standard SPH scheme with Monaghan’s viscosity term.展开更多
基金G.I.Montecinos thanks the National Chilean Fund for Scientific and Technological Development,FONDECYT(Fondo Nacional de Desarrollo Científico y Tecnológico),in the frame of the project for Initiation in Research 11180926
文摘The ADER approach to solve hyperbolic equations to very high order of accuracy has seen explosive developments in the last few years,including both methodological aspects as well as very ambitious applications.In spite of methodological progress,the issues of efficiency and ease of implementation of the solution of the associated generalized Riemann problem(GRP)remain the centre of attention in the ADER approach.In the original formulation of ADER schemes,the proposed solution procedure for the GRP was based on(i)Taylor series expansion of the solution in time right at the element interface,(ii)subsequent application of the Cauchy-Kowalewskaya procedure to convert time derivatives to functionals of space derivatives,and(iii)solution of classical Riemann problems for high-order spatial derivatives to complete the Taylor series expansion.For realistic problems the Cauchy-Kowalewskaya procedure requires the use of symbolic manipulators and being rather cumbersome its replacement or simplification is highly desirable.In this paper we propose a new class of solvers for the GRP that avoid the Cauchy-Kowalewskaya procedure and result in simpler ADER schemes.This is achieved by exploiting the history of the numerical solution that makes it possible to devise a time-reconstruction procedure at the element interface.Still relying on a time Taylor series expansion of the solution at the interface,the time derivatives are then easily calculated from the time-reconstruction polynomial.The resulting schemes are called ADER-TR.A thorough study of the linear stability properties of the linear version of the schemes is carried out using the von Neumann method,thus deducing linear stability regions.Also,via careful numerical experiments,we deduce stability regions for the corresponding non-linear schemes.Numerical examples using the present simplified schemes of fifth and seventh order of accuracy in space and time show that these compare favourably with conventional ADER methods.This paper is restricted to the one-dimensional scalar case with source term,but preliminary results for the one-dimensional Euler equations indicate that the time-reconstruction approach offers significant advantages not only in terms of ease of implementation but also in terms of efficiency for the high-order range schemes.
文摘The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is based on the approach of Vila[25]and another is proposed in this article which mimics the classical 1D LaxWendroff scheme.In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods,demonstrated via the von Neumann stability analysis.Moreover,the issue of the consistency for the equations of gas dynamics is analyzed.An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates.This not only provides an improvement in accuracy of the numerical solutions,but also assures that the consistency condition on the gradient of the kernel function is satisfied using an equidistant distribution of particles in Lagrangian mass coordinates.Three different Riemann solvers are implemented for the first-order Godunov type SPH schemes in Lagrangian coordinates,namely the Godunov flux based on the exact Riemann solver,the Rusanov flux and a new modified Roe flux,following the work of Munz[17].Some well-known numerical 1D shock tube test cases[22]are solved,comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the first-order Godunov finite volume method in Eulerian coordinates and the standard SPH scheme with Monaghan’s viscosity term.
