Pierre Degond

Pierre Degond, Alexei Lozinski, Jacek Narski et al.

The concern of the present work is the introduction of a very efficient Asymptotic Preserving scheme for the resolution of highly anisotropic diffusion equations. The characteristic features of this scheme are the uniform convergence with respect to the anisotropy parameter $0<\eps <<1$, the applicability (on cartesian grids) to cases of non-uniform and non-aligned anisotropy fields $b$ and the simple extension to the case of a non-constant anisotropy intensity $1/\eps$. The mathematical approach and the numerical scheme are different from those presented in the previous work [Degond et al. (2010), arXiv:1008.3405v1] and its considerable advantages are pointed out.

Pierre Degond, Jiale Hua, Laurent Navoret

In this paper, we study the numerical simulations for Euler system with maximal density constraint. This model is developed in [1, 3] with the constraint introduced into the system by a singular pressure law, which causes the transition of different asymptotic dynamics between different regions. To overcome these difficulties, we adapt and implement two asymptotic preserving (AP) schemes originally designed for low Mach number limit [2,4] to our model. These schemes work for the different dynamics and capture the transitions well. Several numerical tests both in one dimensional and two dimensional cases are carried out for our schemes.

Pierre Degond, Fabrice Deluzet, Alexei Lozinski et al.

The present paper introduces an efficient and accurate numerical scheme for the solution of a highly anisotropic elliptic equation, the anisotropy direction being given by a variable vector field. This scheme is based on an asymptotic preserving reformulation of the original system, permitting an accurate resolution independently of the anisotropy strength and without the need of a mesh adapted to this anisotropy. The counterpart of this original procedure is the larger system size, enlarged by adding auxiliary variables and Lagrange multipliers. This Asymptotic-Preserving method generalizes the method investigated in a previous paper [arXiv:0903.4984v2] to the case of an arbitrary anisotropy direction field.

Pierre Degond, Giacomo Dimarco

In the present work, we present a novel numerical algorithm to couple the Direct Simulation Monte Carlo method (DSMC) for the solution of the Boltzmann equation with a finite volume like method for the solution of the Euler equations. Recently we presented in [14],[16],[17] different methodologies which permit to solve fluid dynamics problems with localized regions of departure from thermodynamical equilibrium. The methods rely on the introduction of buffer zones which realize a smooth transition between the kinetic and the fluid regions. In this paper we extend the idea of buffer zones and dynamic coupling to the case of the Monte Carlo methods. To facilitate the coupling and avoid the onset of spurious oscillations in the fluid regions which are consequences of the coupling with a stochastic numerical scheme, we use a new technique which permits to reduce the variance of the particle methods [11]. In addition, the use of this method permits to obtain estimations of the breakdowns of the fluid models less affected by fluctuations and consequently to reduce the kinetic regions and optimize the coupling. In the last part of the paper several numerical examples are presented to validate the method and measure its computational performances.

Pierre Degond, Piotr Minakowski, Laurent Navoret et al.

We are interested in the numerical simulations of the Euler system with variable congestion encoded by a singular pressure. This model describes for instance the macroscopic motion of a crowd with individual congestion preferences. We propose an asymptotic preserving (AP) scheme based on a conservative formulation of the system in terms of density, momentum and density fraction. A second order accuracy version of the scheme is also presented. We validate the scheme on one-dimensional test-cases and extended here to higher order accuracy. We finally carry out two dimensional numerical simulations and show that the model exhibit typical crowd dynamics.

Alina Chertock, Pierre Degond, Jochen Neusser

The Navier-Stokes-Korteweg (NSK) equations are a classical diffuse-interface model for compressible two-phase flow. As direct numerical simulations based on the NSK system are quite expensive and in some cases even impossible, we consider a relaxation of the NSK system, for which robust numerical methods can be designed. However, time steps for explicit numerical schemes depend on the relaxation parameter and therefore numerical simulations in the relaxation limit are very inefficient. To overcome this restriction, we propose an implicit-explicit asymptotic-preserving finite volume method. We prove that the new scheme provides a consistent discretization of the NSK system in the relaxation limit and demonstrate that it is capable of accurately and efficiently computing numerical solutions of problems with realistic density ratios and small interfacial widths.

Stéphane Brull, Pierre Degond, Fabrice Deluzet

This paper is devoted to the numerical approximation of a degenerate anisotropic elliptic problem. The numerical method is designed for arbitrary space-dependent anisotropy directions and does not require any specially adapted coordinate system. It is also designed to be equally accurate in the strongly and the mildly anisotropic cases. The method is applied to the Euler-Lorentz system, in the drift-fluid limit. This system provides a model for magnetized plasmas.

Pierre Degond, Marina A. Ferreira, Sébastien Motsch

We consider algorithms that, from an arbitrarily sampling of $N$ spheres (possibly overlapping), find a close packed configuration without overlapping. These problems can be formulated as minimization problems with non-convex constraints. For such packing problems, we observe that the classical iterative Arrow-Hurwicz algorithm does not converge. We derive a novel algorithm from a multi-step variant of the Arrow-Hurwicz scheme with damping. We compare this algorithm with classical algorithms belonging to the class of linearly constrained Lagrangian methods and show that it performs better. We provide an analysis of the convergence of these algorithms in the simple case of two spheres in one spatial dimension. Finally, we investigate the behaviour of our algorithm when the number of spheres is large.

Pierre Degond, Silke Henkes, Hui Yu

Motivated by recent experimental and computational results that show a motility-induced clustering transition in self-propelled particle systems, we study an individual model and its corresponding Self-Organized Hydrodynamic model for collective behaviour that incorporates a density-dependent velocity, as well as inter-particle alignment. The modal analysis of the hydrodynamic model elucidates the relationship between the stability of the equilibria and the changing velocity, and the formation of clusters. We find, in agreement with earlier results for non-aligning particles, that the key criterion for stability is $(ρv(ρ))'> 0$, i.e. a non-rapid decrease of velocity with density. Numerical simulation for both the individual and hydrodynamic models with a velocity function inspired by experiment demonstrates the validity of the theoretical results.

Pierre Degond, Fabrice Deluzet, David Doyen

In this article, we design Asymptotic-Preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit, this limit being characterized by a Debye length negligible compared to the space scale of the problem. These methods are consistent discretizations of the Vlasov-Maxwell system which, in the quasi-neutral limit, remain stable and are consistent with a quasi-neutral model (in this quasi-neutral model, the electric field is computed by means of a generalized Ohm law). The derivation of Asymptotic-Preserving methods is not straightforward since the quasi-neutral model is a singular limit of the Vlasov-Maxwell model. The key step is a reformulation of the Vlasov-Maxwell system which unifies the two models in a single set of equations with a smooth transition from one to another. As demonstrated in various and demanding numerical simulations, the Asymptotic-Preserving methods are able to treat efficiently both quasi-neutral plasmas and non-neutral plasmas, making them particularly well suited for complex problems involving dense plasmas with localized non-neutral regions.

Pierre Degond, Fabrice Deluzet, Claudia Negulescu

In this article we introduce an asymptotic preserving scheme designed to compute the solution of a two dimensional elliptic equation presenting large anisotropies. We focus on an anisotropy aligned with one direction, the dominant part of the elliptic operator being supplemented with Neumann boundary conditions. A new scheme is introduced which allows an accurate resolution of this elliptic equation for an arbitrary anisotropy ratio.

Pierre Degond, Giacomo Dimarco, Luc Mieussens

This paper collects the efforts done in our previous works [P. Degond, S. Jin, L. Mieussens, A Smooth Transition Between Kinetic and Hydrodynamic Equations, J. Comp. Phys., 209 (2005) 665--694.],[P.Degond, G. Dimarco, L. Mieussens, A Moving Interface Method for Dynamic Kinetic-fluid Coupling, J. Comp. Phys., Vol. 227, pp. 1176-1208, (2007).],[P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)] to build a robust multiscale kinetic-fluid solver. Our scope is to efficiently solve fluid dynamic problems which present non equilibrium localized regions that can move, merge, appear or disappear in time. The main ingredients of the present work are the followings ones: a fluid model is solved in the whole domain together with a localized kinetic upscaling term that corrects the fluid model wherever it is necessary; this multiscale description of the flow is obtained by using a micro-macro decomposition of the distribution function [P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)]; the dynamic transition between fluid and kinetic descriptions is obtained by using a time and space dependent transition function; to efficiently define the breakdown conditions of fluid models we propose a new criterion based on the distribution function itself. Several numerical examples are presented to validate the method and measure its computational efficiency.

Pierre Degond, Giacomo Dimarco, Lorenzo Pareschi

In this work we propose a new approach for the numerical simulation of kinetic equations through Monte Carlo schemes. We introduce a new technique which permits to reduce the variance of particle methods through a matching with a set of suitable macroscopic moment equations. In order to guarantee that the moment equations provide the correct solutions, they are coupled to the kinetic equation through a non equilibrium term. The basic idea, on which the method relies, consists in guiding the particle positions and velocities through moment equations so that the concurrent solution of the moment and kinetic models furnishes the same macroscopic quantities.