Numerical investigation of thin film spreading driven by surfactant using upwind schemes

Numerical solutions of a coupled system of nonlinear partial differential equations modelling the effects of surfactant on the spreading of a thin film on a horizontal substrate are investigated. A CFL condition is obtained from a von Neumann stability analysis of a linearised system of equations. Numerical solutions obtained from a Roe upwind scheme with a third-order TVD Runge-Kutta approximation to the time derivative are compared to solutions obtained with a Roe-Sweby scheme coupled to a minmod limiter and a TVD approximation to the time derivative. Results from both of these schemes are compared to a Roe upwind scheme and a BDF approximation to the time derivative. In all three cases high-order approximations to the spatial derivatives are employed on the interior points of the spatial domain. The Roe-BDF scheme is shown to be an efficient numerical scheme for capturing sharp changes in gradient in the free surface profile and surfactant concentration. Numerical simulations of an initial exponential free surface profile coupled with initial surfactant concentrations for both exogenous and endogenous surfactants are considered.