Convergence of the Two Point Flux Approximation method and the fitted Two Point Flux Approximation method for options pricing with local volatility function

In this paper, we deal with numerical approximations for solving the Black–Scholes Partial Differential Equation (PDE) for European and American options pricing with local volatility. This PDE is well-known to be degenerated. Local volatility model is a model where the volatility depends locally of both stock price and time. In contrast to constant volatility or time-dependent volatility models for which analytical representations of the exact solution is known for European Call options, there is no analytical solution for local volatility. 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 (FTPFA). The Fitted Two-Point Flux Approximation (FTPFA) 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 with the goal to handle the degeneracy of the PDE while the TPFA method is used on the rest of space domain. This combination yields our 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 theoretical results are provided.