Stable and efficient Petrov-Galerkin methods for a kinetic Fokker-Planck equation

by   Julia Brunken, et al.

We propose a stable Petrov-Galerkin discretization of a kinetic Fokker-Planck equation constructed in such a way that uniform inf-sup stability can be inferred directly from the variational formulation. Inspired by well-posedness results for parabolic equations, we derive a lower bound for the dual inf-sup constant of the Fokker-Planck bilinear form by means of stable pairs of trial and test functions. The trial function of such a pair is constructed by applying the kinetic transport operator and the inverse velocity Laplace-Beltrami operator to a given test function. For the Petrov-Galerkin projection we choose an arbitrary discrete test space and then define the discrete trial space using the same application of transport and inverse Laplace-Beltrami operator. As a result, the spaces replicate the stable pairs of the continuous level and we obtain a well-posed numerical method with a discrete inf-sup constant identical to the inf-sup constant of the continuous problem independently of the mesh size. We show how the specific basis functions can be efficiently computed by low-dimensional elliptic problems, and confirm the practicability and performance of the method for a numerical example.



There are no comments yet.


page 1

page 2

page 3

page 4


A Stable FE Method For the Space-Time Solution of the Cahn-Hilliard Equation

In its application to the modeling of a mineral separation process, we p...

A spatial discontinuous Galerkin method with rescaled velocities for the Boltzmann equation

In this paper we present a numerical method for the Boltzmann equation. ...

A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form

We present a stable discontinuous Galerkin (DG) method with a perfectly ...

Estimation and numerical validation of inf-sup constant for bilinear form (p, div u)

We give a derivation for the value of inf-sup constant for the bilinear ...

Fast auxiliary space preconditioners on surfaces

This work presents new uniform preconditioners for the discrete Laplace-...

A Unified Divergence Analysis on Projection Method for Symm's Integral Equation of the First Kind

A large amount of literatures analyze (SIE) which concerns the construct...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.