Conservative polynomial approximations and applications to Fokker-Planck equations

Abstract

We address the problem of constructing approximations based on orthogonal polynomials that preserve an arbitrary set of moments of a given function without loosing the spectral convergence property. To this aim, we compute the constrained polynomial of best approximation for a generic basis of orthogonal polynomials. The construction is entirely general and allows us to derive structure preserving numerical methods for partial differential equations that require the conservation of some moments of the solution, typically representing relevant physical quantities of the problem. These properties are essential to capture with high accuracy the long-time behavior of the solution. We illustrate with the aid of several numerical applications to Fokker-Planck equations the generality and the performances of the present approach.

26 pages, 17 figures, 7 tables

References in corpus (6)

• Fast algorithms for computing the Boltzmann collision operator
• Structure preserving schemes for nonlinear Fokker-Planck equations and applications
• Metriplectic Integrators for the Landau Collision Operator
• Fully Discrete Positivity-Preserving and Energy-Dissipating Schemes for Aggregation-Diffusion Equations with a Gradient Flow Structure
• Convolutive decomposition and fast summation methods for discrete-velocity approximations of the Boltzmann equation
• Convergence of spectral discretizations of the Vlasov-Poisson system