A mathematical model is considered to determine the effectiveness of disin-
fectant solution for surface decontamination. The decontamination process involved the
diffusion of bacteria into disinfectant solution and the reaction of the disinfectant killing
effect. The mathematical model is a reaction-diffusion type. Finite difference method and
method of lines with fourth-order Runge-Kutta method are utilized to solve the model
numerically. To obtain stable solutions, von Neumann stability analysis is employed to
evaluate the stability of finite difference method. For stiff problem, Dormand-Prince
method is applied as the estimated error of fourth-order Runge-Kutta method. MATLAB
programming is selected for the computation of numerical solutions. From the results
obtained, fourth-order Runge-Kutta method has a larger stability region and better ac-
curacy of solutions compared to finite difference method when solving the disinfectant
solution model. Moreover, a numerical simulation is carried out to investigate the effect
of different thickness of disinfectant solution on bacteria reduction. Results show that
thick disinfectant solution is able to reduce the dimensionless bacteria concentration more
effectively.