Novel Numerical Investigation of Reaction Diffusion Equation Arising in Oil Price Modeling

Abstract

Consideration is given to a reaction–diffusion free boundary value problem with one or two turning points arising in oil price modeling. First, an exact (analytical) solution to the reduced problem (i.e., no diffusion term) was obtained for some given parameters. The space–time Chebyshev pseudospectral and superconsistent Chebyshev collocation method is proposed for both reaction diffusion (RDFBP) and reduced free boundary value problem. Error bounds on the discrete L2–norm and Sobolev norm (Hp) are presented. Adaptively graded intervals were introduced and used according to the value of turning points to avoid the twin boundary layers phenomena. Excellent convergent (spectrally) and stable results for some special turning points were obtained for both reduced and RDFBP equations on an adaptively graded interval and this has been documented for the first time.