Giuseppe Orlando

Luca Arpaia, Giuseppe Orlando, Christian Ferrarin et al.

We present a discontinuous finite element method for the shallow water equations which exploits high-resolution realistic bathymetry data without any regularity assumption, also in the case of high-order discretizations. We prove a number of mathematical properties specific to the proposed method that is well-balanced, mass-conserving and positivity-preserving under a mild CFL condition also in the presence of wet-dry fronts. The method includes a consistent conservative discretization for passive tracers. We use a high-order Discontinuous Galerkin (DG) method as implemented in the deal.II library. This environment provides efficient and native parallelization techniques and automatically handles non-conforming meshes to implement adaptive strategies which are tested in a coastal environment. Idealized test cases show the robustness in presence of irregular bathymetries also with under-resolved features at the grid scale. A benchmark with realistic bathymetry and a complex domain shows the potential of the proposed discretization for adaptive simulations of coastal flows.

Letizia Bottani, Tommaso Benacchio, Giuseppe Orlando et al.

We present a high-order implicit-explicit discontinuous Galerkin (IMEX-DG) solver for the compressible Euler equations to account for rotational effects within a fully compressible atmospheric framework. Time integration follows a second-order additive Runge-Kutta scheme, treating stiff acoustic modes implicitly and advective terms explicitly. The solver is built on the deal.II finite element library, combining matrix-free operator evaluation, adaptive non-conforming meshes capabilities, and distributed-memory parallelism. Two alternative treatments of the rotational and gravitational source terms within the solution strategy, based on nonlinear fixed-point iterations, are introduced and compared in terms of accuracy, robustness, and computational efficiency. A discrete analysis of the rotational operator is also carried out in order to derive a formulation suitable for efficient matrix-free implementation and to avoid inconsistent naive discretisations. The proposed formulation is validated through convergence studies on rotating inertia-gravity wave benchmarks and further assessed in fully three-dimensional simulations of stratified flow over orography on both uniform and adaptive meshes. The numerical results show that the rotating IMEX-DG framework has the expected accuracy and stability properties while correctly capturing the asymmetry and wave structures induced by rotation in large-scale atmospheric flows.