BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Maurizio Tavelli

Maurizio Tavelli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisMathematical Software

Catalog footprint

What is connected

5works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

An all Froude high order IMEX scheme for the shallow water equations on unstructured Voronoi meshes

We propose a novel numerical method for the solution of the shallow water equations in different regimes of the Froude number making use of general polygonal meshes. The fluxes of the governing equations are split such that advection and acoustic-gravity sub-systems are derived, hence separating slow and fast phenomena. This splitting allows the nonlinear convective fluxes to be discretized explicitly in time, while retaining an implicit time marching for the acoustic-gravity terms. Consequently, the novel schemes are particularly well suited in the low Froude limit of the model, since no numerical viscosity is added in the implicit solver. Besides, stability follows from a milder CFL condition which is based only on the advection speed and not on the celerity. High order time accuracy is achieved using the family of semi-implicit IMEX Runge-Kutta schemes, while high order in space is granted relying on two discretizations: (i) a cell-centered finite volume (FV) scheme for the nonlinear convective contribution on the polygonal cells; (ii) a staggered discontinuous Galerkin (DG) scheme for the solution of the linear system associated to the implicit discretization of the pressure sub-system. Therefore, three different meshes are used, namely a polygonal Voronoi mesh, a triangular subgrid and a staggered quadrilateral subgrid. The novel schemes are proved to be Asymptotic Preserving (AP), hence a consistent discretization of the limit model is retrieved for vanishing Froude numbers, which is the given by the so-called "lake at rest" equations. Furthermore, the novel methods are well-balanced by construction, and this property is also demonstrated. Accuracy and robustness are then validated against a set of benchmark test cases with Froude numbers ranging in the interval $\Fr \approx [10^{-6};5]$, hence showing that multiple time scales can be handled by the novel methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems

ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for solving systems of first-order hyperbolic partial differential equations (PDEs). Hyperbolic PDEs are typically derived from the conservation laws of physics and are useful in a wide range of application areas. Applications powered by ExaHyPE can be run on a student's laptop, but are also able to exploit thousands of processor cores on state-of-the-art supercomputers. The engine is able to dynamically increase the accuracy of the simulation using adaptive mesh refinement where required. Due to the robustness and shock capturing abilities of ExaHyPE's numerical methods, users of the engine can simulate linear and non-linear hyperbolic PDEs with very high accuracy. Users can tailor the engine to their particular PDE by specifying evolved quantities, fluxes, and source terms. A complete simulation code for a new hyperbolic PDE can often be realised within a few hours - a task that, traditionally, can take weeks, months, often years for researchers starting from scratch. In this paper, we showcase ExaHyPE's workflow and capabilities through real-world scenarios from our two main application areas: seismology and astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A staggered space-time discontinuous Galerkin method for the three-dimensional incompressible Navier-Stokes equations on unstructured tetrahedral meshes

In this paper we propose a novel arbitrary high order accurate semi-implicit space-time DG method for the solution of the three-dimensional incompressible Navier-Stokes equations on staggered unstructured curved tetrahedral meshes. As typical for space-time DG schemes, the discrete solution is represented in terms of space-time basis functions. This allows to achieve very high order of accuracy also in time, which is not easy to obtain for the incompressible Navier-Stokes equations. Similar to staggered finite difference schemes, in our approach the discrete pressure is defined on the primary tetrahedral grid, while the discrete velocity is defined on a face-based staggered dual grid. A very simple and efficient Picard iteration is used in order to derive a space-time pressure correction algorithm that achieves also high order of accuracy in time and that avoids the direct solution of global nonlinear systems. Formal substitution of the discrete momentum equation on the dual grid into the discrete continuity equation on the primary grid yields a very sparse five-point block system for the scalar pressure, which is conveniently solved with a matrix-free GMRES algorithm. From numerical experiments we find that the linear system seems to be reasonably well conditioned, since all simulations shown in this paper could be run without the use of any preconditioner. For a piecewise constant polynomial approximation in time and proper boundary conditions, the resulting system is symmetric and positive definite. This allows us to use even faster iterative solvers, like the conjugate gradient method. The proposed method is verified for approximation polynomials of degree up to four in space and time by solving a series of typical 3D test problems and by comparing the obtained numerical results with available exact analytical solutions, or with other numerical or experimental reference data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A staggered semi-implicit discontinuous Galerkin method for the two dimensional incompressible Navier-Stokes equations

In this paper we propose a new spatially high order accurate semi-implicit discontinuous Galerkin (DG) method for the solution of the two dimensional incompressible Navier-Stokes equations on staggered unstructured curved meshes. While the discrete pressure is defined on the primal grid, the discrete velocity vector field is defined on an edge-based dual grid. The flexibility of high order DG methods on curved unstructured meshes allows to discretize even complex physical domains on rather coarse grids. Formal substitution of the discrete momentum equation into the discrete continuity equation yields one sparse block four-diagonal linear equation system for only one scalar unknown, namely the pressure. The method is computationally efficient, since the resulting system is not only very sparse but also symmetric and positive definite for appropriate boundary conditions. Furthermore, all the volume and surface integrals needed by the scheme presented in this paper depend only on the geometry and the polynomial degree of the basis and test functions and can therefore be precomputed and stored in a preprocessor stage, which leads to savings in terms of computational effort for the time evolution part. In this way also the extension to a fully curved isoparametric approach becomes natural and affects only the preprocessing step. The method is validated for polynomial degrees up to $p=3$ by solving some typical numerical test problems and comparing the numerical results with available analytical solutions or other numerical and experimental reference data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A staggered space-time discontinuous Galerkin method for the incompressible Navier-Stokes equations on two-dimensional triangular meshes

In this paper we propose a novel arbitrary high order accurate semi-implicit space-time discontinuous Galerkin method for the solution of the two dimensional incompressible Navier-Stokes equations on staggered unstructured triangular meshes. Isoparametric finite elements are used to take into account curved domain boundaries. The discrete pressure is defined on the primal triangular grid and the discrete velocity field is defined on an edge-based staggered dual grid. While staggered meshes are state of the art in classical finite difference approximations of the incompressible Navier-Stokes equations, their use in the context of high order DG schemes is novel and still quite rare. Formal substitution of the discrete momentum equation into the discrete continuity equation yields a sparse four-point block system for the scalar pressure, which is conveniently solved with a matrix-free GMRES algorithm. A very simple and efficient Picard iteration is then used in order to achieve high order of accuracy also in time, which is in general a non-trivial task in the context of high order discretizations for the incompressible Navier-Stokes equations. The flexibility and accuracy of high order space-time DG methods on curved unstructured meshes allows to discretize even complex physical domains with very coarse grids in both, space and time. The use of a staggered grid allows to avoid the use of Riemann solvers in several terms of the discrete equations and significantly reduces the total stencil size of the linear system that needs to be solved for the pressure. The proposed method is validated for approximation polynomials of degree up to $p=4$ in space and time by solving a series of typical numerical test problems and by comparing the obtained numerical results with available exact analytical solutions or other numerical reference data.

0 citations0 reviewsNo code linkedOpen access