Source author record

Gabriella Puppo

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

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

6works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

An all Mach number finite volume method for isentropic two-phase flow

We present an implicit-explicit finite volume scheme for isentropic two phase flow in all Mach number regimes. The underlying model belongs to the class of symmetric hyperbolic thermodynamically compatible models. The key element of the scheme consists of a linearisation of pressure and enthalpy terms at a reference state. The resulting stiff linear parts are integrated implicitly, whereas the non-linear higher order and transport terms are treated explicitly. Due to the flux splitting, the scheme is stable under a CFL condition which determined by the resolution of the slow material waves and allows large time steps even in the presence of fast acoustic waves. Further the singular Mach number limits of the model are studied and the asymptotic preserving property of the scheme is proven. In numerical simulations the consistency with single phase flow, accuracy and the approximation of material waves in different Mach number regimes are assessed.

preprint2019arXiv

An all speed second order IMEX relaxation scheme for the Euler equations

We present an implicit-explicit finite volume scheme for the Euler equations. We start from the non-dimensionalised Euler equations where we split the pressure in a slow and a fast acoustic part. We use a Suliciu type relaxation model which we split in an explicit part, solved using a Godunov-type scheme based on an approximate Riemann solver, and an implicit part where we solve an elliptic equation for the fast pressure. The relaxation source terms are treated projecting the solution on the equilibrium manifold. The proposed scheme is positivity preserving with respect to the density and internal energy and asymptotic preserving towards the incompressible Euler equations. For this first order scheme we give a second order extension which maintains the positivity property. We perform numerical experiments in 1D and 2D to show the applicability of the proposed splitting and give convergence results for the second order extension.

preprint2019arXiv

An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity

We present an implicit-explicit well-balanced finite volume scheme for the Euler equations with a gravitational source term which is able to deal also with low Mach flows. To visualize the different scales we use the non-dimensionalized equations on which we apply a pressure splitting and a Suliciu relaxation. On the resulting model, we apply a splitting of the flux into a linear implicit and an non-linear explicit part that leads to a scale independent time-step. The explicit step consists of a Godunov type method based on an approximative Riemann solver where the source term is included in the flux formulation. We develop the method for a first order scheme and give an extension to second order. Both schemes are designed to be well-balanced, preserve the positivity of density and internal energy and have a scale independent diffusion. We give the low Mach limit equations for well-prepared data and show that the scheme is asymptotic preserving. These properties are numerically validated by various test cases.

preprint2016arXiv

Kinetic models for traffic flow resulting in a reduced space of microscopic velocities

The purpose of this paper is to study the properties of kinetic models for traffic flow described by a Boltzmann-type approach and based on a continuous space of microscopic velocities. In our models, the particular structure of the collision kernel allows one to find the analytical expression of a class of steady-state distributions, which are characterized by being supported on a quantized space of microscopic speeds. The number of these velocities is determined by a physical parameter describing the typical acceleration of a vehicle and the uniqueness of this class of solutions is supported by numerical investigations. This shows that it is possible to have the full richness of a kinetic approach with the simplicity of a space of microscopic velocities characterized by a small number of modes. Moreover, the explicit expression of the asymptotic distribution paves the way to deriving new macroscopic equations using the closure provided by the kinetic model.

preprint2015arXiv

Fundamental diagrams in traffic flow: the case of heterogeneous kinetic models

Experimental studies on vehicular traffic provide data on quantities like density, flux, and mean speed of the vehicles. However, the diagrams relating these variables (the fundamental and speed diagrams) show some peculiarities not yet fully reproduced nor explained by mathematical models. In this paper, resting on the methods of kinetic theory, we introduce a new traffic model which takes into account the heterogeneous nature of the flow of vehicles along a road. In more detail, the model considers traffic as a mixture of two or more populations of vehicles (e.g., cars and trucks) with different microscopic characteristics, in particular different lengths and/or maximum speeds. With this approach we gain some insights into the scattering of the data in the regime of congested traffic clearly shown by actual measurements.

preprint2014arXiv

Well-balanced high order schemes on non-uniform grids and entropy residuals

This paper is concerned with the construction of high order schemes on irregular grids for balance laws, including a discussion of an a-posteriori error indicator based on the numerical entropy production. We also impose well-balancing on non uniform grids for the shallow water equations, which can be extended similarly to other cases, obtaining schemes up to fourth order of accuracy with very weak assumptions on the regularity of the grid. Our results show the expected convergence rates, the correct propagation of shocks across grid discontinuities and demonstrate the improved resolution achieved with a locally refined non-uniform grid. The error indicator based on the numerical entropy production, previously introduced for the case of systems of conservation laws, is extended to balance laws. Its decay rate and its ability to identify discontinuities is illustrated on several tests. The schemes proposed in this work naturally can also be applied to systems of conservation laws.