Synchronized flow and wide moving jams from balanced vehicular traffic

Recently we proposed an extension to the traffic model of Aw, Rascle and Greenberg. The extended traffic model can be written as a hyperbolic system of balance laws and numerically reproduces the reverse $\lambda$ shape of the fundamental diagram of traffic flow. In the current work we analyze the steady state solutions of the new model and their stability properties. In addition to the equilibrium flow curve the trivial steady state solutions form two additional branches in the flow-density diagram. We show that the characteristic structure excludes parts of these branches resulting in the reverse $\lambda$ shape of the flow-density relation. The upper branch is metastable against the formation of synchronized flow for intermediate densities and unstable for high densities, whereas the lower branch is unstable for intermediate densities and metastable for high densities. Moreover, the model can reproduce the typical speed of the downstream front of wide moving jams. It further reproduces a constant outflow from wide moving jams, which is far below the maximum free flow. Applying the model to simulate traffic flow at a bottleneck we observe a general pattern with wide moving jams traveling through the bottleneck.Comment: 10 pages, 12 figure

Uniqueness of the compactly supported weak solutions of the relativistic Vlasov-Darwin system

We use optimal transportation techniques to show uniqueness of the compactly supported weak solutions of the relativistic Vlasov-Darwin system. Our proof extends the method used by Loeper in J. Math. Pures Appl. 86, 68-79 (2006) to obtain uniqueness results for the Vlasov-Poisson system.Comment: AMS-LaTeX, 21 page

Numerical modeling of landquakes

International audienceThe Thurwieser landslide that occurred in Italy in 2004 and the seismic waves it generated have been simulated and compared to the seismic signal recorded a few tens of kilometers from the landslide source (i.e., landquake). The main features of the low frequency seismic signal are reproduced by the simulation. Topography effects on the flowing mass have a major impact on the generated seismic signal whereas they weakly affect low-frequency wave propagation. Simulation of the seismic signal makes it possible to discriminate between possible alternative scenarios for flow dynamics and to provide first estimates of the rheological parameters during the flow. As landquakes are continuously recorded by seismic networks, our results provide a new way to collect data on the dynamics and rheology of natural flows

Phase appearance or disappearance in two-phase flows

This paper is devoted to the treatment of specific numerical problems which appear when phase appearance or disappearance occurs in models of two-phase flows. Such models have crucial importance in many industrial areas such as nuclear power plant safety studies. In this paper, two outstanding problems are identified: first, the loss of hyperbolicity of the system when a phase appears or disappears and second, the lack of positivity of standard shock capturing schemes such as the Roe scheme. After an asymptotic study of the model, this paper proposes accurate and robust numerical methods adapted to the simulation of phase appearance or disappearance. Polynomial solvers are developed to avoid the use of eigenvectors which are needed in usual shock capturing schemes, and a method based on an adaptive numerical diffusion is designed to treat the positivity problems. An alternate method, based on the use of the hyperbolic tangent function instead of a polynomial, is also considered. Numerical results are presented which demonstrate the efficiency of the proposed solutions

Numerical approximation of the Euler-Poisson-Boltzmann model in the quasineutral limit

This paper analyzes various schemes for the Euler-Poisson-Boltzmann (EPB) model of plasma physics. This model consists of the pressureless gas dynamics equations coupled with the Poisson equation and where the Boltzmann relation relates the potential to the electron density. If the quasi-neutral assumption is made, the Poisson equation is replaced by the constraint of zero local charge and the model reduces to the Isothermal Compressible Euler (ICE) model. We compare a numerical strategy based on the EPB model to a strategy using a reformulation (called REPB formulation). The REPB scheme captures the quasi-neutral limit more accurately

On Nonlinear Stochastic Balance Laws

We are concerned with multidimensional stochastic balance laws. We identify a class of nonlinear balance laws for which uniform spatial $BV$ bounds for vanishing viscosity approximations can be achieved. Moreover, we establish temporal equicontinuity in $L^1$ of the approximations, uniformly in the viscosity coefficient. Using these estimates, we supply a multidimensional existence theory of stochastic entropy solutions. In addition, we establish an error estimate for the stochastic viscosity method, as well as an explicit estimate for the continuous dependence of stochastic entropy solutions on the flux and random source functions. Various further generalizations of the results are discussed

Probabilistic analysis of the upwind scheme for transport

We provide a probabilistic analysis of the upwind scheme for multi-dimensional transport equations. We associate a Markov chain with the numerical scheme and then obtain a backward representation formula of Kolmogorov type for the numerical solution. We then understand that the error induced by the scheme is governed by the fluctuations of the Markov chain around the characteristics of the flow. We show, in various situations, that the fluctuations are of diffusive type. As a by-product, we prove that the scheme is of order 1/2 for an initial datum in BV and of order 1/2-a, for all a>0, for a Lipschitz continuous initial datum. Our analysis provides a new interpretation of the numerical diffusion phenomenon

Simulation of a liquid-vapour compressible flow by a Lattice Boltzmann Method

International audienceThis work is devoted to the numerical resolution of a compressible three-phase flow with phase transition by a Lattice-Boltzmann Method (LBM). The flow presents complex features and large variations of physical quantities. The LBM is a robust numerical method that is entropy stable and that can be extended to second order accuracy without additional numerical cost. We present preliminary numerical results, which confirm its competitiveness compared to other Finite Volume methods

A splitting approach for the fully nonlinear and weakly dispersive Green-Naghdi model

The fully nonlinear and weakly dispersive Green-Naghdi model for shallow water waves of large amplitude is studied. The original model is first recast under a new formulation more suitable for numerical resolution. An hybrid finite volume and finite difference splitting approach is then proposed. The hyperbolic part of the equations is handled with a high-order finite volume scheme allowing for breaking waves and dry areas. The dispersive part is treated with a classical finite difference approach. Extensive numerical validations are then performed in one horizontal dimension, relying both on analytical solutions and experimental data. The results show that our approach gives a good account of all the processes of wave transformation in coastal areas: shoaling, wave breaking and run-up

A mathematical model for unsteady mixed flows in closed water pipes

We present the formal derivation of a new unidirectional model for unsteady mixed flows in non uniform closed water pipe. In the case of free surface incompressible flows, the \FS-model is formally obtained, using formal asymptotic analysis, which is an extension to more classical shallow water models. In the same way, when the pipe is full, we propose the \Pres-model, which describes the evolution of a compressible inviscid flow, close to gas dynamics equations in a nozzle. In order to cope the transition between a free surface state and a pressured (i.e. compressible) state, we propose a mixed model, the \PFS-model, taking into account changes of section and slope variation