Accurate spectral numerical schemes for kinetic equations with energy diffusion
From MaRDI portal
Publication:349963
DOI10.1016/J.JCP.2015.03.039zbMATH Open1349.76652OpenAlexW1997773080MaRDI QIDQ349963FDOQ349963
orthogonal polynomialscontinuous spectrumcontinuum kinetic calculationsFokker-Plank collisionsSturm-Liouville theory
Stochastic methods (Fokker-Planck, Langevin, etc.) applied to problems in time-dependent statistical mechanics (82C31) Statistical mechanics of plasmas (82D10) Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs (65M70) Turbulent transport, mixing (76F25) Ionized gas flow in electromagnetic fields; plasmic flow (76X05)
Abstract: We examine the merits of using a family of polynomials that are orthogonal with respect to a non-classical weight function to discretize the speed variable in continuum kinetic calculations. We consider a model one-dimensional partial differential equation describing energy diffusion in velocity space due to Fokker-Planck collisions. This relatively simple case allows us to compare the results of the projected dynamics with an expensive but highly accurate spectral transform approach. It also allows us to integrate in time exactly, and to focus entirely on the effectiveness of the discretization of the speed variable. We show that for a fixed number of modes or grid points, the non-classical polynomials can be many orders of magnitude more accurate than classical Hermite polynomials or finite-difference solvers for kinetic equations in plasma physics. We provide a detailed analysis of the difference in behavior and accuracy of the two families of polynomials. For the non-classical polynomials, if the initial condition is not smooth at the origin when interpreted as a three-dimensional radial function, the exact solution leaves the polynomial subspace for a time, but returns (up to roundoff accuracy) to the same point evolved to by the projected dynamics in that time. By contrast, using classical polynomials, the exact solution differs significantly from the projected dynamics solution when it returns to the subspace. We also explore the connection between eigenfunctions of the projected evolution operator and (non-normalizable) eigenfunctions of the full evolution operator, as well as the effect of truncating the computational domain.
Recommendations
- New velocity-space discretization for continuum kinetic calculations and Fokker-Planck collisions
- Numerical Solution of the Boltzmann Equation I: Spectrally Accurate Approximation of the Collision Operator
- scientific article; zbMATH DE number 1494171
- scientific article; zbMATH DE number 1426678
- The spectral collocation method for the kinetic equation with the nonlinear two-dimensional Coulomb collisional operator
Cites work
- scientific article; zbMATH DE number 1183040 (Why is no real title available?)
- scientific article; zbMATH DE number 1049347 (Why is no real title available?)
- scientific article; zbMATH DE number 1131587 (Why is no real title available?)
- scientific article; zbMATH DE number 1987549 (Why is no real title available?)
- A Gaussian quadrature procedure for use in the solution of the Boltzmann equation and related problems
- A direct method for the Boltzmann equation based on a pseudo-spectral velocity space discretization
- A spectral transform method for singular Sturm-Liouville problems with applications to energy diffusion in plasma physics
- Accurate spectral numerical schemes for kinetic equations with energy diffusion
- Additive Runge-Kutta schemes for convection-diffusion-reaction equations
- An Eulerian gyrokinetic-Maxwell solver.
- An algorithm for computing Jordan chains and inverting analytic matrix functions
- AstroGK: astrophysical gyrokinetics code
- Comparison of initial value and eigenvalue codes for kinetic toroidal plasma instabilities
- Construction of Gauss-Christoffel Quadrature Formulas
- Fourth-Order Time-Stepping for Stiff PDEs
- Half-Range Generalized Hermite Polynomials and the Related Gaussian Quadratures
- New velocity-space discretization for continuum kinetic calculations and Fokker-Planck collisions
- On Generating Orthogonal Polynomials
- On the Nature of the Spectrum of Singular Second Order Linear Differential Equations
- Singular Second-Order Operators: The Maximal and Minimal Operators, and Selfadjoint Operators in Between
- Titchmarsh-Weyl m -functions for second-order Sturm-Liouville problems with two singular endpoints
Cited in
(8)- An efficient numerical method for solving the Boltzmann equation in multidimensions
- A spectral transform method for singular Sturm-Liouville problems with applications to energy diffusion in plasma physics
- Accurate spectral numerical schemes for kinetic equations with energy diffusion
- Relic neutrino freeze-out: dependence on natural constants
- Kinetically consistent high-accuracy schemes.
- Computing the Dirichlet-Neumann operator on a cylinder
- Krylov iterative methods and synthetic acceleration for transport in binary statistical media
- The unified transform for mixed boundary condition problems in unbounded domains
This page was built for publication: Accurate spectral numerical schemes for kinetic equations with energy diffusion
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q349963)
