Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

A Structure-Preserving Penalization Method for the Single-species Rosenbluth-Fokker-Planck Equation

Hamad El Kahza, Luis Chacón, William Taitano, Jingmei Qiu, Jingwei Hu

TL;DR

The paper develops a structure-preserving explicit-implicit-null (EIN) penalization framework for the stiff, nonlinear single-species Rosenbluth-Fokker-Planck equation. It combines a generalized Chang-Cooper spatial discretization that preserves the Maxwellian equilibrium with a variable-coefficient, conservatively reformulated penalization operator to stabilize time integration, and an adaptive positivity-preserving timestep. The approach enforces strict conservation of mass, momentum, and energy, preserves the equilibrium null space, and guarantees a discrete maximum principle, yielding unconditionally stable, positivity-preserving simulations across linear and nonlinear anisotropic diffusion tests and Rosenbluth-Fokker-Planck relaxation scenarios. Numerical results demonstrate robust performance, accurate relaxation rates, and sustained positivity even with large timesteps, highlighting potential for multi-species extensions and Vlasov coupling. The work advances reliable, structure-aware numerical solvers for Coulomb collisional dynamics in plasmas.

Abstract

The Rosenbluth-Fokker-Planck (RFP) equation describes Coulomb collisional dynamics within and across species in plasmas. It belongs to the broader class of anisotropic-diffusion-advection equations, whose numerical approximation is highly-nontrivial due to its nonlinearity, stiffness, and structural properties such as conservation and entropy dissipation (hence with the Maxwellian distribution as the equilibrium state). In this paper, we propose a structure-preserving penalization scheme for the stiff, single-species RFP equation. The scheme features three novel components: 1) a novel generalization of the well-known Chang-Cooper discretization for the RFP equation that is equilibrium-preserving and enables positivity while preserving mass, momentum, and energy; 2) an easy-to-invert isotropic variable-coefficient penalization operator to deal with the temporal stiffness without resorting to a fully implicit scheme, borrowing ideas from explicit-implicit-null (EIN) methods, and 3) an adaptive timestepping strategy that preserves the positivity of the full penalized scheme. The resulting scheme conserves mass, momentum, and energy strictly, is unconditionally stable, and robustly positivity preserving. The scheme is demonstrated with linear and nonlinear anisotropic diffusion examples of increasing complexity, including several single-species RFP examples.

A Structure-Preserving Penalization Method for the Single-species Rosenbluth-Fokker-Planck Equation

TL;DR

The paper develops a structure-preserving explicit-implicit-null (EIN) penalization framework for the stiff, nonlinear single-species Rosenbluth-Fokker-Planck equation. It combines a generalized Chang-Cooper spatial discretization that preserves the Maxwellian equilibrium with a variable-coefficient, conservatively reformulated penalization operator to stabilize time integration, and an adaptive positivity-preserving timestep. The approach enforces strict conservation of mass, momentum, and energy, preserves the equilibrium null space, and guarantees a discrete maximum principle, yielding unconditionally stable, positivity-preserving simulations across linear and nonlinear anisotropic diffusion tests and Rosenbluth-Fokker-Planck relaxation scenarios. Numerical results demonstrate robust performance, accurate relaxation rates, and sustained positivity even with large timesteps, highlighting potential for multi-species extensions and Vlasov coupling. The work advances reliable, structure-aware numerical solvers for Coulomb collisional dynamics in plasmas.

Abstract

The Rosenbluth-Fokker-Planck (RFP) equation describes Coulomb collisional dynamics within and across species in plasmas. It belongs to the broader class of anisotropic-diffusion-advection equations, whose numerical approximation is highly-nontrivial due to its nonlinearity, stiffness, and structural properties such as conservation and entropy dissipation (hence with the Maxwellian distribution as the equilibrium state). In this paper, we propose a structure-preserving penalization scheme for the stiff, single-species RFP equation. The scheme features three novel components: 1) a novel generalization of the well-known Chang-Cooper discretization for the RFP equation that is equilibrium-preserving and enables positivity while preserving mass, momentum, and energy; 2) an easy-to-invert isotropic variable-coefficient penalization operator to deal with the temporal stiffness without resorting to a fully implicit scheme, borrowing ideas from explicit-implicit-null (EIN) methods, and 3) an adaptive timestepping strategy that preserves the positivity of the full penalized scheme. The resulting scheme conserves mass, momentum, and energy strictly, is unconditionally stable, and robustly positivity preserving. The scheme is demonstrated with linear and nonlinear anisotropic diffusion examples of increasing complexity, including several single-species RFP examples.
Paper Structure (30 sections, 145 equations, 17 figures, 2 algorithms)

This paper contains 30 sections, 145 equations, 17 figures, 2 algorithms.

Table of Contents

  1. Introduction
  2. The single-species Rosenbluth-Fokker-Planck equation
  3. EIN penalization method for anisotropic diffusion
  4. Stability analysis of the semi-discrete scheme
  5. Positivity analysis of the semi-discrete scheme
  6. Adaptive positivity-preserving time stepping
  7. Maximum-principle-preserving spatial discretization
  8. Structure-preserving EIN penalization scheme for the single-species RFP operator
  9. Conservative penalization approach
  10. Null-space preservation property of conservative penalization operators $\mathcal{L}_\beta$
  11. Structure-preserving spatial discretization of the EIN-penalized RFP equation
  12. Structure-preserving discretization of the penalization operator $\mathcal{L}_\beta$
  13. Structure-preserving spatial discretization of the RFP operator
  14. Maximum-principle- and null-space-preserving RFP spatial discretization
  15. Strict discrete conservation strategy in RFP
  16. ...and 15 more sections

Figures (17)

  • Figure 1: Evolution of the solution for the anisotropic ring test at four representative times.
  • Figure 2: (a) Minimum value of the solution in time for a fixed time-step ratio $\Delta t = 100\,\Delta t_{\mathrm{CFL}}$. (b) Minimum value of the solution in time for the adaptive time-stepping strategy. (c) Evolution of the adaptive time-step ratio $N_{CFL}=\Delta t_n / \Delta t_{\mathrm{CFL}}$ as a function of time.
  • Figure 3: Log-log temporal and spatial convergence of the penalized implicit scheme for the variable-coefficient anisotropic diffusion test. (a) The $L_2$ error at the final time is plotted as a function of the time-step size $\Delta t$, together with a reference line of slope 1. (b) Spatial convergence study: the $L_2$ error at the final time is plotted as a function of the mesh size, together with the corresponding reference line of slope 2.
  • Figure 4: Evolution of the solution for the anisotropic ring test at four representative times.
  • Figure 5: (a) Minimum value of the solution in time for a fixed time-step ratio with $N_{CFL}=100$. (b) Minimum value of the solution in time for the adaptive time-stepping strategy. (c) Evolution of the adaptive time-step ratio $\Delta t_n / \Delta t_{\mathrm{CFL}}$ as a function of time.
  • ...and 12 more figures