An alternative scheme to calculate the strain rate tensor for the LES applications in the LBM.
Discrete-velocity approximations represent a popular way for computing the Boltzmann collision operator. The direct numerical evaluation of such methods involve a prohibitive cost, typically O(N2d + 1) where d is the dimension of the velocity space. In this paper, following the ideas introduced in [C. Mouhot and L. Pareschi, C. R. Acad. Sci. Paris Sér. I Math. 339 (2004) 71–76, C. Mouhot and L. Pareschi, Math. Comput. 75 (2006) 1833–1852], we derive fast summation techniques for the evaluation of...
In this work, we consider the computation of the boundary conditions for the linearized Euler–Poisson derived from the BGK kinetic model in the small mean free path regime. Boundary layers are generated from the fact that the incoming kinetic flux might be far from the thermodynamical equilibrium. In [2], the authors propose a method to compute numerically the boundary conditions in the hydrodynamic limit relying on an analysis of the boundary layers....
We propose here a model and a numerical scheme to compute the motion of rigid particles interacting through the lubrication force. In the case of a particle approaching a plane, we propose an algorithm and prove its convergence towards the solutions to the gluey particle model described in [B. Maury, ESAIM: Proceedings 18 (2007) 133–142]. We propose a multi-particle version of this gluey model which is based on the projection of the velocities onto a set of admissible velocities. Then, we describe...
We propose here a model and a numerical scheme to compute the motion of rigid particles interacting through the lubrication force. In the case of a particle approaching a plane, we propose an algorithm and prove its convergence towards the solutions to the gluey particle model described in [B. Maury, ESAIM: Proceedings18 (2007) 133–142]. We propose a multi-particle version of this gluey model which is based on the projection of the velocities onto a set of admissible velocities. Then, we describe...
We present a hybrid finite-volume-particle numerical method for computing the transport of a passive pollutant by a flow. The flow is modeled by the one- and two-dimensional Saint-Venant system of shallow water equations and the pollutant propagation is described by a transport equation. This paper is an extension of our previous work [Chertock, Kurganov and Petrova, J. Sci. Comput. (to appear)], where the one-dimensional finite-volume-particle method has been proposed. The core idea behind the...
We present a hybrid finite-volume-particle numerical method for computing the transport of a passive pollutant by a flow. The flow is modeled by the one- and two-dimensional Saint-Venant system of shallow water equations and the pollutant propagation is described by a transport equation. This paper is an extension of our previous work [Chertock, Kurganov and Petrova, J. Sci. Comput. (to appear)], where the one-dimensional finite-volume-particle method has been proposed. The core idea behind the...
A new numerical scheme called particle-in-wavelets is proposed for the Vlasov-Poisson equations, and tested in the simplest case of one spatial dimension. The plasma distribution function is discretized using tracer particles, and the charge distribution is reconstructed using wavelet-based density estimation. The latter consists in projecting the Delta distributions corresponding to the particles onto a finite dimensional linear space spanned by...
The aim of this paper is to present a finite volume kinetic method to compute the transport of a passive pollutant by a flow modeled by the shallow water equations using a new time discretization that allows large time steps for the pollutant computation. For the hydrodynamic part the kinetic solver ensures – even in the case of a non flat bottom – the preservation of the steady state of a lake at rest, the non-negativity of the water height and the existence of an entropy inequality. On an other...
The aim of this paper is to present a finite volume kinetic method to compute the transport of a passive pollutant by a flow modeled by the shallow water equations using a new time discretization that allows large time steps for the pollutant computation. For the hydrodynamic part the kinetic solver ensures – even in the case of a non flat bottom – the preservation of the steady state of a lake at rest, the non-negativity of the water height and the existence of an entropy inequality. On an other...
Both the porous medium equation and the system of isentropic Euler equations can be considered as steepest descents on suitable manifolds of probability measures in the framework of optimal transport theory. By discretizing these variational characterizations instead of the partial differential equations themselves, we obtain new schemes with remarkable stability properties. We show that they capture successfully the nonlinear features of the flows, such as shocks and rarefaction waves for...