Convergence of the Two Point Flux Approximation and a novel fitted Two-Point Flux Approximation method for pricing options

In this paper, we deal with numerical approximations for solving the Black-Scholes Partial Differential Equation (PDE). This PDE is well known to be degenerated. The space discretization is performed using the classical finite volume method with Two-Point Flux Approximation (TPFA) and a novel scheme called fitted Two-Point Flux Approximation (TPFA).The fitted Two-Point Flux Approximation (TPFA) combines the fitted finite volume method and the standard TPFA method. More precisely the fitted finite volume method is used when the stock price approaches zero in order to handle the degeneracy of the PDE and the TPFA method is used on the rest of space domain. This combination yields our novel fitted TPFA scheme. The Euler method is used for the time discretization. We provide the rigorous convergence proofs of the two fully discretized schemes. Numerical experiments to support our theoretical results are provided.