Abstract:We introduce a family of various finite volume discretization schemes for the Fokker--Planck operator, which are characterized by different weight functions on the edges. This family particularly includes the well-established Scharfetter--Gummel discretization as well as the recently developed square-root approximation (SQRA) scheme. We motivate this family of discretizations both from the numerical and the modeling point of view and provide a uniform consistency and error analysis. Our main results state that the convergence order primarily depends on the quality of the mesh and in second place on the quality of the weights. We show by numerical experiments that for small gradients the choice of the optimal representative of the discretization family is highly non-trivial while for large gradients the Scharfetter--Gummel scheme stands out compared to the others.

What problem does this paper attempt to address?

### What problems does this paper attempt to solve? This paper aims to study and analyze a series of finite - volume methods for discretizing the Fokker - Planck operator, with particular attention to the influence of different edge - weight functions on this discretization family. Specifically, the author introduces a class of finite - volume discretization schemes with different weight functions based on the Stolarsky mean and conducts a unified analysis of the consistency and convergence of these schemes. #### Main problems include: 1. **Discretization consistency**: - Investigate whether this family of discretization schemes can consistently approximate the continuous Fokker - Planck equation as the grid is refined. - Analyze the contribution of different weight functions in the discretization scheme to consistency. 2. **Convergence order**: - Explore that the convergence order of the discretization scheme mainly depends on the grid quality and the weight quality. - Provide numerical experiments to verify the selection of the optimal discretization scheme in different gradient situations. 3. **Comparison between the Scharfetter - Gummel (SG) and Square - Root Approximation (SQRA) schemes**: - The SG scheme performs excellently in the strong - gradient region, while the SQRA scheme may be better in the weak - gradient region. - Study the respective advantages and disadvantages of these two schemes as special cases of the Stolarsky - mean discretization family. 4. **Generalized Stolarsky - mean discretization scheme**: - Introduce and analyze the generalized discretization scheme based on the Stolarsky mean. - Propose a new discretization framework that includes not only the existing SG and SQRA schemes but also other potential discretization methods. 5. **Error estimation and numerical simulation**: - Give the error estimate of the discretization scheme and verify the theoretical results through numerical experiments. - Show the influence of different parameter selections (such as α and β) on the discretization effect. #### Specific mathematical formulas The Fokker - Planck equation can be expressed as: \[ -\nabla\cdot(\kappa\nabla u)-\nabla\cdot(\kappa u\nabla V) = f \] where \( \kappa>0 \) is the diffusion coefficient, \( V:\Omega\rightarrow\mathbb{R} \) is the given potential function, and \( f \) describes the possible source or sink terms. The discretized flux \( J(u, V) \) can be written as: \[ J(u, V)=-\kappa(\nabla u + u\nabla V) \] For the SG and SQRA schemes, the discretized forms are as follows: \[ B(V_i, V_j)= \begin{cases} \frac{V_i - V_j}{e^{V_i - V_j}-1}&\text{(SG)}\\ e^{-\frac{1}{2}(V_i - V_j)}&\text{(SQRA)} \end{cases} \] These formulas are used to describe the flux calculation methods under different discretization schemes. #### Summary The main objective of this paper is to reveal the consistency and convergence properties of a series of finite - volume discretization schemes based on the Stolarsky mean through systematic analysis, thereby providing theoretical basis and practical guidance for selecting the optimal discretization scheme.

Consistency and convergence for a family of finite volume discretizations of the Fokker--Planck operator

Discretizing the Fokker-Planck equation with second-order accuracy: a dissipation driven approach

A New High-Order Compact Scheme Of Unstructured Finite Volume Method

Maximum-Principle-Satisfying Third Order Discontinuous Galerkin Schemes for Fokker-Planck Equations

The Entropy Satisfying Discontinuous Galerkin Method for Fokker–Planck Equations

Convergence and Quasi-Optimality of Adaptive FEM with Inhomogeneous Dirichlet Data

On the order of accuracy of finite-volume schemes on unstructured meshes

The Finite Volume Scheme Preserving Maximum Principle for Two-Dimensional Time-Fractional Fokker–Planck Equations on Distorted Meshes

Adaptive vertex-centered finite volume methods with convergence rates

Convergent discretizations for the Nernst–Planck–Poisson system

Strong Convergence of the JKO Scheme for the Fokker–Planck Equation

Nonstandard finite difference schemes for linear and non-linear Fokker–Planck equations

Spectral Analysis of High Order Continuous FEM for Hyperbolic PDEs on Triangular Meshes: Influence of Approximation, Stabilization, and Time-Stepping

The design of conservative finite element discretisations for the vectorial modified KdV equation

Finite element discretizations for variable-order fractional diffusion problems

A centered limited finite volume approximation of the momentum convection operator for low‐order nonconforming face‐centered discretizations

Superconvergent Non-Polynomial Approximations

A Mixed FEM for a Time-Fractional Fokker–Planck Model

Long-time behavior for discretization schemes of Fokker-Planck equations via couplings

On full linear convergence and optimal complexity of adaptive FEM with inexact solver

Convergence rates for a finite volume scheme of the stochastic heat equation