1,721,026 research outputs found
An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity
Full text linkWe 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
Derivation and Stability Analysis of a Macroscopic Multilane Model for Traffic Flow
No full textThe mathematical modeling and the stability analysis of multilane traffic in the macroscopic scale is considered. We propose a new first order model derived from microscopic dynamics with lane changing, leading to a coupled system of hyperbolic balance laws. The macroscopic limit is derived without assuming ad hoc space and time scalings. The analysis of the stability of the equilibria of the model is discussed. The proposed numerical tests confirm the theoretical findings between the macroscopic and microscopic modeling, and the results of the stability analysis
Arbitrary order finite volume well-balanced schemes for the Euler equations with gravity
Full text linkThis work presents arbitrary high order well balanced finite volume schemes for the Euler equations with a prescribed gravitational field. It is assumed that the desired equilibrium solution is known, and we construct a scheme which is exactly well balanced for that particular equilibrium. The scheme is based on high order reconstructions of the fluctuations from equilibrium of density, velocity, and pressure, and on a well-balanced integration of the source terms, while no assumptions are needed on the numerical flux, beside consistency. This technique also allows one to construct well-balanced methods for a class of moving equilibria. Several numerical tests demonstrate the performance of the scheme on different scenarios, from equilibrium solutions to nonsteady problems involving shocks. The numerical tests are carried out with methods up to fifth order in one dimension, and third order accuracy in two dimensions
A consistent kinetic model for a two-component mixture of polyatomic molecules
Full text linkWe consider a multi component gas mixture with translational and internal energy degrees of freedom assuming that the number of particles of each species remains constant. We will illustrate the derived model in the case of two species, but the model can be easily generalized to multiple species. The two species are allowed to have different degrees of freedom in internal energy and are modelled by a system of kinetic ES-BGK equations featuring two interaction terms to account for momentum and energy transfer between the species. We prove consistency of our model: conservation properties, positivity of the temperature, H-theorem and convergence to a global equilibrium in the form of a global Maxwell distribution. Thus, we are able to derive the usual macroscopic conservation laws. For numerical purposes we apply the Chu reduction to the developed model for polyatomic gases and give an application for a gas consisting of a monoatomic and a diatomic species
Quinpi: Integrating Stiff Hyperbolic Systems with Implicit High Order Finite Volume Schemes
No full textMany interesting physical problems described by systems of hyperbolic conservation laws are stiff, and thus impose a very small time-step because of the restrictive CFL stability condition. In this case, one can exploit the superior stability properties of implicit time integration which allows to choose the time-step only from accuracy requirements, and thus avoid the use of small time-steps. We discuss an efficient framework to devise high order implicit schemes for stiff hyperbolic systems without tailoring it to a specific problem. The nonlinearity of high order schemes, due to space- and time-limiting procedures which control nonphysical oscillations, makes the implicit time integration difficult, e.g. because the discrete system is nonlinear also on linear problems. This nonlinearity of the scheme is circumvented as proposed in (Puppo et al., Comm. Appl. Math. & Comput., 2023) for scalar conservation laws, where a first order implicit predictor is computed to freeze the nonlinear coefficients of the essentially non-oscillatory space reconstruction, and also to assist limiting in time. In addition, we propose a novel conservative flux-centered a-posteriori time-limiting procedure using numerical entropy indicators to detect troubled cells. The numerical tests involve classical and artificially devised stiff problems using the Euler's system of gas-dynamics
Adaptivity for systems of conservation laws driven by an error indicator based on entropy production
No full textin: La matematica e le sue applicazioni, 13-2009, Dip. di Matematica, Politecnico di Torino. ISSN 1974-305
- …
