BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Francis Filbet

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

ResearcherUnclaimed source record
math.NAmath.APNumerical Analysismath-phmath.MPComputational Engineering, Finance, and Sciencephysics.flu-dyn

Catalog footprint

What is connected

27works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

27 published item(s)

preprint2022arXiv

Concentration phenomena in Fitzhugh-Nagumo's equations: A mesoscopic approach

We consider a spatially extended mesoscopic FitzHugh-Nagumo model with strong local interactions and prove that its asymptotic limit converges towards the classical nonlocal reaction-diffusion FitzHugh-Nagumo system. As the local interactions strongly dominate, the weak solution to the mesoscopic equation under consideration converges to the local equilibrium, which has the form of Dirac distribution concentrated to an averaged membrane potential. Our approach is based on techniques widely developed in kinetic theory (Wasserstein distance, relative entropy method), where macroscopic quantities of the mesoscopic model are compared with the solution to the nonlocal reaction-diffusion system. This approach allows to make the rigorous link between microscopic and reaction-diffusion models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the stability of conservative discontinuous Galerkin/Hermite spectral methods for the Vlasov-Poisson system

We study a class of spatial discretizations for the Vlasov-Poisson system written as an hyperbolic system using Hermite polynomials. In particular, we focus on spectral methods and discontinuous Galerkin approximations. To obtain L 2 stability properties, we introduce a new L 2 weighted space, with a time dependent weight. For the Hermite spectral form of the Vlasov-Poisson system, we prove conservation of mass, momentum and total energy, as well as global stability for the weighted L 2 norm. These properties are then discussed for several spatial discretizations. Finally, numerical simulations are performed with the proposed DG/Hermite spectral method to highlight its stability and conservation features.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A spectral collocation method for the Landau equation in plasma physics

In this paper we present a spectral collocation method for the fast evaluation of the Landau collision operator for plasma physics, which allows us to obtain spectrally accurate numerical solutions. The method is inspired by the seminal work [36], but it is specifically designed for Coulombian interactions, taking into account the particular structure of the operator. It allows us to reduce the number of discrete convolutions to provide an approximation of the Landau operator. Then, we show that the method preserves the total mass whereas momentum and energy are approximated with spectral accuracy. Numerical results for the Landau equation in three dimensions in velocity space are presented to illustrate the efficiency of the present approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotics of the three dimensional Vlasov equation in the large magnetic field limit

We study the asymptotic behavior of solutions to the Vlasov equation in the presence of a strong external magnetic field. In particular we provide a mathematically rigorous derivation of the guiding-center approximation in the general three dimensional setting under the action of large inhomogeneous magnetic fields. First order corrections are computed and justified as well, including electric cross field, magnetic gradient and magnetic curvature drifts. We also treat long time behaviors on two specific examples, the two dimensional case in carte-sian coordinates and a poloidal axi-symmetric geometry, the former for expository purposes. Algebraic manipulations that underlie concrete computations make the most of the linearity of the stiffest part of the system of characteristics instead of relying on any particular variational structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Conservative discontinuous Galerkin/Hermite Spectral Method for the Vlasov-Poisson System

We propose a class of conservative discontinuous Galerkin methods for the Vlasov-Poisson system written as a hyperbolic system using Hermite polynomials in the velocity variable. These schemes are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in general are hard to achieve within other numerical method frameworks for simulating the Vlasov-Poisson system. The proposed scheme employs discontinuous Galerkin discretization for both the Vlasov and the Poisson equations, resulting in a consistent description of the distribution function and electric field. Numerical simulations are performed to verify the order of accuracy and conservation properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence Analysis of Asymptotic Preserving Schemes for Strongly Magnetized plasmas

The present paper is devoted to the convergence analysis of a class of asymptotic preserving particle schemes [Filbet \& Rodrigues, SIAM J. Numer. Anal., 54 (2) (2016)] for the Vlasov equation with a strong external magnetic field. In this regime, classical Particle-In-Cell (PIC) methods are subject to quite restrictive stability constraints on the time and space steps, due to the small Larmor radius and plasma frequency. The asymptotic preserving discretization that we are going to study removes such a constraint while capturing the large-scale dynamics, even when the discretization (in time and space) is too coarse to capture fastest scales. Our error bounds are explicit regarding the discretization, stiffness parameter, initial data and time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence of a Finite Volume Scheme for a System of Interacting Species with Cross-Diffusion

In this work we present the convergence of a positivity preserving semi-discrete finite volume scheme for a coupled system of two non-local partial differential equations with cross-diffusion. The key to proving the convergence result is to establish positivity in order to obtain a discrete energy estimate to obtain compactness. We numerically observe the convergence to reference solutions with a first order accuracy in space. Moreover we recover segregated stationary states in spite of the regularising effect of the self-diffusion. However, if the self-diffusion or the cross-diffusion is strong enough, mixing occurs while both densities remain continuous.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A p-adaptive local discontinuous galerkin level set method for Willmore flow

The level set method is often used to capture interface behavior in two or three dimensions. In this paper, we present a combination of local discontinuous Galerkin (LDG) method and level set method for simulating Willmore flow. The LDG scheme is energy stable and mass conservative, which are good properties comparing with other numerical methods. In addition, to enhance the efficiency of the proposed LDG scheme and level set method, we employ a p-adaptive local discontinuous Galerkin technique, which applies high order polynomial approximations around the zero level set and low order ones away from the zero level set. A major advantage of the level set method is that the topological changes are well defined and easily performed. In particular, given the stiffness of Willmore flow, a high order semi-implicit Runge-Kutta method is employed for time discretization, which allows larger time step. The equations at the implicit time level are linear, we demonstrate an efficient and practical multi-grid solver to solve the equations. Numerical examples are given to illustrate the combination of the LDG scheme and level set method provides an efficient and practical approach when simulating the Willmore flow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the asymptotic limit of the three dimensional Vlasov-Poisson system for large magnetic field : formal derivation

This paper establishes the long time asymptotic limit of the three dimensional Vlasov-Poisson equation with strong external magnetic field. The guiding center approximation is investigated in the three dimensional case with a non-constant magnetic field. In the long time asymptotic limit, the motion can be split in two parts : one stationary flow along the lines of the magnetic field and the guiding center motion in the orthogonal plane of the magnetic field where classical drift velocities are recovered. We discuss in particular the effect of nonconstant external magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymptotically stable particle-in-cell methods for the Vlasov-Poisson system with a strong external magnetic field

This paper deals with the numerical resolution of the Vlasov-Poissonsystem with a strong external magnetic field by Particle-In-Cell(PIC) methods. In this regime, classical PIC methods are subject tostability constraints on the time and space steps related to the smallLarmor radius and plasma frequency. Here, we propose anasymptotic-preserving PIC scheme which is not subjected to theselimitations. Our approach is based on first and higher order semi-implicit numericalschemes already validated on dissipative systems. Additionally, when the magnitude of the external magneticfield becomes large, this method provides a consistent PICdiscretization of the guiding-center equation, that is, incompressibleEuler equation in vorticity form. We propose several numerical experiments which provide a solid validation of the method and its underlying concepts.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Hierarchy of Hybrid Numerical Methods for Multi-Scale Kinetic Equations

In this paper, we construct a hierarchy of hybrid numerical methods for multi-scale kinetic equations based on moment realizability matrices, a concept introduced by Levermore, Morokoff and Nadiga. Following such a criterion, one can consider hybrid scheme where the hydrodynamic part is given either by the compressible Euler or Navier-Stokes equations, or even with more general models, such as the Burnett or super-Burnett systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mixed semi-Lagrangian/finite difference methods for plasma simulations

In this paper, we present an efficient algorithm for the long time behavior of plasma simulations. We will focus on 4D drift-kinetic model, where the plasma's motion occurs in the plane perpendicular to the magnetic field and can be governed by the 2D guiding-center model. Hermite WENO reconstructions, already proposed in \cite{YF15}, are applied for solving the Vlasov equation. Here we consider an arbitrary computational domain with an appropriate numerical method for the treatment of boundary conditions. Then we apply this algorithm for plasma turbulence simulations. We first solve the 2D guiding-center model in a D-shape domain and investigate the numerical stability of the steady state. Then, the 4D drift-kinetic model is studied with a mixed method, i.e. the semi-Lagrangian method in linear phase and finite difference method during the nonlinear phase. Numerical results show that the mixed method is efficient and accurate in linear phase and it is much stable during the nonlinear phase. Moreover, in practice it has better conservation properties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On discrete functional inequalities for some finite volume schemes

We prove several discrete Gagliardo-Nirenberg-Sobolev and Poincaré-Sobolev inequalities for some approximations with arbitrary boundary values on finite volume meshes. The keypoint of our approach is to use the continuous embedding of the space $BV(Ω)$ into $L^{N/(N-1)}(Ω)$ for a Lipschitz domain $ Ω\subset \mathbb{R}^{N}$, with $N \geq 2$. Finally, we give several applications to discrete duality finite volume (DDFV) schemes which are used for the approximation of nonlinear and non isotropic elliptic and parabolic problems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Rescaling Velocity Method for Dissipative Kinetic Equations - Applications to Granular Media

We present a new numerical algorithm based on a relative energy scaling for collisional kinetic equations allowing to study numerically their long time behavior, without the usual problems related to the change of scales in velocity variables. It is based on the knowledge of the hydrodynamic limit of the model considered, but is able to compute solutions for either dilute or dense regimes. Several applications are presented for Boltzmann like equations. This method is particularly efficient for numerical simulations of the granular gases equation with dissipative energy: it allows to study accurately the long time behavior of this equation and is very well suited for the study of clustering phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conservative and non-conservative methods based on hermite weighted essentially-non-oscillatory reconstruction for Vlasov equations

We introduce a WENO reconstruction based on Hermite interpolation both for semi-Lagrangian and finite difference methods. This WENO reconstruction technique allows to control spurious oscillations. We develop third and fifth order methods and apply them to non-conservative semi-Lagrangian schemes and conservative finite difference methods. Our numerical results will be compared to the usual semi-Lagrangian method with cubic spline reconstruction and the classical fifth order WENO finite difference scheme. These reconstructions are observed to be less dissipative than the usual weighted essentially non- oscillatory procedure. We apply these methods to transport equations in the context of plasma physics and the numerical simulation of turbulence phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical Simulations of Kinetic Models for chemotaxis

We present a new algorithm based on a Cartesian mesh for the numerical approximation of kinetic models for chemosensitive movements set in an arbitrary geometry. We investigate the influence of the geometry on the collective behavior of bacteria described by a kinetic equation interacting with nutrients and chemoattractants. Numerical simulations are performed to verify accuracy and stability of the scheme and its ability to exhibit aggregation of cells and wave propagations. Finally some comparisons with experiments show the robustness and accuracy of such kinetic models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Inverse Lax-Wendroff method for boundary conditions of Boltzmann type models

In this paper we present a new algorithm based on a Cartesian mesh for the numerical approximation of kinetic models on complex geometry boundary. Due to the high dimensional property, numerical algorithms based on unstructured meshes for a complex geometry are not appropriate. Here we propose to develop an inverse Lax-Wendroff pro- cedure, which was recently introduced for conservation laws [S. Tan & C.W. Shu, JCP (2010)], to the kinetic equations. Applications in 1Dx3D and 2Dx3D of this algorithm for Boltzmann type operators (BGK, ES-BGK models) are then presented and numerical results illustrate the accuracy properties of this algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of an Asymptotic Preserving Scheme for Relaxation Systems

We study the convergence of a class of asymptotic preserving numerical schemes initially proposed by F. Filbet & S. Jin \cite{filb1} and G. Dimarco & L. Pareschi \cite{DimarcoP} in the context of nonlinear and stiff kinetic equations. Here, our analysis is devoted to the approximation of a system of transport equations with a nonlinear source term, for which the asymptotic limit is given by a conservation laws. We investigate the convergence of the approximate solution $(\ueps_h,\veps_h)$ to a nonlinear relaxation system, where $\eps>0$ is a physical parameter and $h$ represents the discretization parameter. Uniform convergence with respect to $\eps$ and $h$ is proven and error estimates are also obtained. Finally, several numerical tests are performed to illustrate the accuracy and efficiency of such a scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical study of a nonlinear heat equation for plasma physics

This paper is devoted to the numerical approximation of a nonlinear temperature balance equation, which describes the heat evolution of a magnetically confined plasma in the edge region of a tokamak. The nonlinearity implies some numerical difficulties, in particular long time behavior, when solved with standard methods. An efficient numerical scheme is presented in this paper, based on a combination of a directional splitting scheme and the IMEX scheme introduced in [Filbet and Jin]

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On deterministic approximation of the Boltzmann equation in a bounded domain

In this paper we present a fully deterministic method for the numerical solution to the Boltzmann equation of rarefied gas dynamics in a bounded domain for multi-scale problems. Periodic, specular reflection and diffusive boundary conditions are discussed and investigated numerically. The collision operator is treated by a Fourier approximation of the collision integral, which guarantees spectral accuracy in velocity with a computational cost of $N\,\log(N)$, where $N$ is the number of degree of freedom in velocity space. This algorithm is coupled with a second order finite volume scheme in space and a time discretization allowing to deal for rarefied regimes as well as their hydrodynamic limit. Finally, several numerical tests illustrate the efficiency and accuracy of the method for unsteady flows (Poiseuille flows, ghost effects, trend to equilibrium).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Numerical Scheme for the Quantum Boltzmann Equation Efficient in the Fluid Regime

Numerically solving the Boltzmann kinetic equations with the small Knudsen number is challenging due to the stiff nonlinear collision term. A class of asymptotic preserving schemes was introduced in [6] to handle this kind of problems. The idea is to penalize the stiff collision term by a BGK type operator. This method, however, encounters its own difficulty when applied to the quantum Boltzmann equation. To define the quantum Maxwellian (Bose-Einstein or Fermi- Dirac distribution) at each time step and every mesh point, one has to invert a nonlinear equation that connects the macroscopic quantity fugacity with density and internal energy. Setting a good initial guess for the iterative method is troublesome in most cases because of the complexity of the quantum functions (Bose-Einstein or Fermi-Dirac function). In this paper, we propose to penalize the quantum collision term by a 'classical' BGK operator instead of the quantum one. This is based on the observation that the classical Maxwellian, with the temperature replaced by the internal energy, has the same first five moments as the quantum Maxwellian. The scheme so designed avoids the aforementioned difficulty, and one can show that the density distribution is still driven toward the quantum equilibrium. Numerical results are present to illustrate the efficiency of the new scheme in both the hydrodynamic and kinetic regimes. We also develop a spectral method for the quantum collision operator.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An Asymptotic Preserving Scheme for the ES-BGK model

In this paper, we study a time discrete scheme for the initial value problem of the ES-BGK kinetic equation. Numerically solving these equations are challenging due to the nonlinear stiff collision (source) terms induced by small mean free or relaxation time. We study an implicit-explicit (IMEX) time discretization in which the convection is explicit while the relaxation term is implicit to overcome the stiffness. We first show how the implicit relaxation can be solved explicitly, and then prove asymptotically that this time discretization drives the density distribution toward the local Maxwellian when the mean free time goes to zero while the numerical time step is held fixed. This naturally imposes an asymptotic-preserving scheme in the Euler limit. The scheme so designed does not need any nonlinear iterative solver for the implicit relaxation term. Moreover, it can capture the macroscopic fluid dynamic (Euler) limit even if the small scale determined by the Knudsen number is not numerically resolved. We also show that it is consistent to the compressible Navier-Stokes equations if the viscosity and heat conductivity are numerically resolved. Several numerical examples, in both one and two space dimensions, are used to demonstrate the desired behavior of this scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources

In this paper, we propose a general framework to design asymptotic preserving schemes for the Boltzmann kinetic kinetic and related equations. Numerically solving these equations are challenging due to the nonlinear stiff collision (source) terms induced by small mean free or relaxation time. We propose to penalize the nonlinear collision term by a BGK-type relaxation term, which can be solved explicitly even if discretized implicitly in time. Moreover, the BGK-type relaxation operator helps to drive the density distribution toward the local Maxwellian, thus natually imposes an asymptotic-preserving scheme in the Euler limit. The scheme so designed does not need any nonlinear iterative solver or the use of Wild Sum. It is uniformly stable in terms of the (possibly small) Knudsen number, and can capture the macroscopic fluid dynamic (Euler) limit even if the small scale determined by the Knudsen number is not numerically resolved. It is also consistent to the compressible Navier-Stokes equations if the viscosity and heat conductivity are numerically resolved. The method is applicable to many other related problems, such as hyperbolic systems with stiff relaxation, and high order parabilic equations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Analysis of spectral methods for the homogeneous Boltzmann equation

The development of accurate and fast algorithms for the Boltzmann collision integral and their analysis represent a challenging problem in scientific computing and numerical analysis. Recently, several works were devoted to the derivation of spectrally accurate schemes for the Boltzmann equation, but very few of them were concerned with the stability analysis of the method. In particular there was no result of stability except when the method is modified in order to enforce the posivity preservation, which destroys the spectral accuracy. In this paper we propose a new method to study the stability of homogeneous Boltzmann equations perturbed by smoothed balanced operators which do not preserve positivity of the distribution. This method takes advantage of the "spreading" property of the collision, together with estimates on regularity and entropy production. As an application we prove stability and convergence of spectral methods for the Boltzmann equation, when the discretization parameter is large enough (with explicit bound).

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

High order resolution of the Maxwell-Fokker-Planck-Landau model intended for ICF applications

A high order, deterministic direct numerical method is proposed for the nonrelativistic $2D_{\bf x} \times 3D_{\bf v}$ Vlasov-Maxwell system, coupled with Fokker-Planck-Landau type operators. Such a system is devoted to the modelling of electronic transport and energy deposition in the general frame of Inertial Confinement Fusion applications. It describes the kinetics of plasma physics in the nonlocal thermodynamic equilibrium regime. Strong numerical constraints lead us to develop specific methods and approaches for validation, that might be used in other fields where couplings between equations, multiscale physics, and high dimensionality are involved. Parallelisation (MPI communication standard) and fast algorithms such as the multigrid method are employed, that make this direct approach be computationally affordable for simulations of hundreds of picoseconds, when dealing with configurations that present five dimensions in phase space.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Solving the Boltzmann equation in N log N

In [C. Mouhot and L. Pareschi, "Fast algorithms for computing the Boltzmann collision operator," Math. Comp., to appear; C. Mouhot and L. Pareschi, C. R. Math. Acad. Sci. Paris, 339 (2004), pp. 71-76], fast deterministic algorithms based on spectral methods were derived for the Boltzmann collision operator for a class of interactions including the hard spheres model in dimension three. These algorithms are implemented for the solution of the Boltzmann equation in two and three dimension, first for homogeneous solutions, then for general non homogeneous solutions. The results are compared to explicit solutions, when available, and to Monte-Carlo methods. In particular, the computational cost and accuracy are compared to those of Monte-Carlo methods as well as to those of previous spectral methods. Finally, for inhomogeneous solutions, we take advantage of the great computational efficiency of the method to show an oscillation phenomenon of the entropy functional in the trend to equilibrium, which was suggested in the work [L. Desvillettes and C. Villani, Invent. Math., 159 (2005), pp. 245-316].

0 citations0 reviewsNo code linkedOpen access