The aim of this study is to develop a combined approximation technique to nd a numerical solution to the foam drainage equation in various absorption and distillation processes. In this approach, rst, discretization of time is performed with the aid of the Taylor expansion series. Hence, a collocation method based on novel Bessel polynomials is utilized for the space variable. Thus, the solution is found by solving a linear system of algebraic equations in each time step. Numerical simulations are provided to check the accuracy and eciency of the presented algorithm. The numerical results are compared with exact solutions as well as the outcomes of other existing numerical methods.