Abstract In this paper we propose a new reformulation of the first order hyperbolic model for unsteady turbulent shallow water flows recently proposed in Gavrilyuk et al. (J Comput Phys 366:252–280, 2018). The novelty formulation forwarded here is use evolution variable that guarantees trace discrete Reynolds stress tensor to be always non-negative. mathematical particularly challenging because...

After introducing the Discontinuous Galerkin (DG) method a detailed misfit analysis on its numerical approximation is performed. We investigate the accuracy of the scheme, the element type (tetrahedrons and hexahedrons), the spatial sampling of the computational domain and the number of propagated wavelengths. As the error norm we chose a time-frequency representation, which illustrates the tim...

We further investigate the high order positivity-preserving discontinuous Galerkin (DG) methods for linear hyperbolic and radiative transfer equations developed in [14]. The DG methods in [14] can maintain positivity and high order accuracy, but they rely both on the scaling limiter in [15] and a rotational limiter, the latter may alter cell averages of the unmodulated DG scheme, thereby affect...

Many applications from geosciences require simulations of seismic waves in porous media. Biot's theory poroelasticity describes the coupling between solid and fluid phases introduces a stiff source term, thereby increasing computational cost motivating efficient methods utilising High-Performance Computing. We present novel realisation discontinuous Galerkin scheme with Arbitrary DERivative tim...

We develop a high-order positivity-preserving discontinuous Galerkin (DG) scheme for linear Vlasov-Boltzmann transport equations (Vlasov-BTE) under the action of quadratically confined electrostatic potentials. The solutions of such BTEs are positive probability distribution functions and it is very challenging to have a mass-conservative, high-order accurate scheme that preserves positivity of...

In this paper we carry out the extension of the ADER approach to multidimensional non-linear systems of conservation laws. We implement non-linear schemes of up to fourth order of accuracy in both time and space. Numerical results for the compressible Euler equations illustrate the very high order of accuracy and non-oscillatory properties of the new schemes. Compared to the state-of-art finite...

Abstract. We present a numerical scheme for the solution of a class of atmospheric models where high horizontal resolution is required while a coarser vertical structure is allowed. The proposed scheme considers a layering procedure for the original set of equations, and the use of high-order ADER finite volume schemes for the solution of the system of balance laws arising from the dimensional ...