In this paper, a new numerical technique is proposed for the simulations of advection-diffusion-reaction type elliptic and parabolic interface models. The proposed technique comprises of the Haar wavelet collocation method and the finite difference method. In this technique, the spatial derivative is approximated by truncated Haar wavelet series, while for temporal derivative, the finite difference formula is used. The diffusion coefficients, advection coefficients, and reaction coefficients are considered discontinuously across the fixed interface. The newly established numerical technique is applied to both linear and nonlinear benchmark interface models. In the case of linear interface models, Gauss elimination method is used, whereas for nonlinear interface models, the nonlinearity is removed by using the quasi-Newton linearization technique. The
L
∞
errors are calculated for different number of collocation points. The obtained numerical results are compared with the immersed interface method. The stability and convergence of the method are also discussed. On the whole, the numerical results show more efficiency, better accuracy, and simpler applicability of the newly developed numerical technique compared to the existing methods in literature.