The main task of numerical relativity is to solve Einstein equations with the aid of supercomputer. There are two main schemes to write Einstein equations explicitly as differential equations. One is based on 3 + 1 decomposition and reduces the Einstein equations to a Cauchy problem. The another takes the advantage of the characteristic property of Einstein equations and reduces them to a set of ordinary differential equations. The latter scheme is called characteristic formalism which is free of constraint equations in contrast to the corresponding Cauchy problem. Till now there is only one well developed code (PITT code) for characteristic formalism. In PITT code, special finite difference algorithm is adopted for the numerical calculation. And it is this special difference algorithm that restricts the numerical accuracy order to second-order. In addition, this special difference algorithm makes the popular Runge–Kutta method used in Cauchy problem not available. In this paper, we modify the equations for characteristic formalism. Based on our new set of equations, we can use usual finite difference method as done in usual Cauchy evolution. And Runge–Kutta method can also be adopted naturally. We develop a set of code in the framework of AMSS-NCKU code based on our new method and some numerical tests are done.
Numerical Treatments for Nonlinear Integro-Fractional Differential Equations of Volterra-Hammerstein Type using Runge-Kutta Method with the aid of Finite Difference Approximation