A Fractional-Differential Approach to Numerical Simulation of Electron-Induced Charging of Ferroelectrics

The paper proposes a fractional-differential modification of the mathematical model of the process of nonstationary charging of polar dielectric materials under conditions of irradiation with medium-energy electron beams. The mathematical formalization is based on a spherically symmetric diffusion–drift equation with a fractional time derivative. An implicit finite-difference scheme is constructed using the Caputo derivative approximation. An application program has been developed in Matlab software that implements the designed computational algorithm. Verification of an approximate solution of the problem is demonstrated using a test example. The results of computational experiments to evaluate the characteristics of field effects of injected charges in ferroelectrics when varying the order of fractional differentiation in subdiffusion regimes are presented.