Accurate spectral numerical schemes for kinetic equations with energy diffusion
From MaRDI portal
(Redirected from Publication:349963)
continuous spectrumcontinuum kinetic calculationsFokker-Plank collisionsorthogonal polynomialsSturm-Liouville theory
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) Stochastic methods (Fokker-Planck, Langevin, etc.) applied to problems in time-dependent statistical mechanics (82C31) Statistical mechanics of plasmas (82D10)
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
- A direct method for the Boltzmann equation based on a pseudo-spectral velocity space discretization
- A Gaussian quadrature procedure for use in the solution of the Boltzmann equation and related problems
- 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 algorithm for computing Jordan chains and inverting analytic matrix functions
- An Eulerian gyrokinetic-Maxwell solver.
- 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
- 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?)
- 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
(10)- Krylov iterative methods and synthetic acceleration for transport in binary statistical media
- An efficient numerical method for solving the Boltzmann equation in multidimensions
- Kinetically consistent high-accuracy schemes.
- Accurate spectral numerical schemes for kinetic equations with energy diffusion
- The unified transform for mixed boundary condition problems in unbounded domains
- Computing the Dirichlet-Neumann operator on a cylinder
- A spectral transform method for singular Sturm-Liouville problems with applications to energy diffusion in plasma physics
- Computing generalized eigenfunctions in rigged Hilbert spaces
- Nonequilibrium effects for reactions with activation energy: convergence of the expansions of solutions of the Boltzmann and Lorentz Fokker Planck equations with Sonine and Maxwell polynomials as basis functions
- Relic neutrino freeze-out: dependence on natural constants
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)
