Some Numerical Methods and Comparisons for Solving Mathematical Model of Surface Decontamination by Disinfectant Solution

A mathematical model is considered to determine the effectiveness of disinfectant 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 accuracy 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