Researcher profile

Francis Filbet

Francis Filbet contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 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.

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.

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.

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.

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.

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.

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.

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.

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.

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]

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).

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.

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.

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).