Numerical solutions of a coupled system of nonlinear partial differential equationsmodelling the ... more Numerical solutions of a coupled system of nonlinear partial differential equationsmodelling 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 timederivative are compared to solutions obtained with aRoe-Sweby scheme coupled to a minmod limiter and a TVD approximation to the time derivative. Results fromboth 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.TheRoe-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.
Numerical solutions of a coupled system of nonlinear partial differential equationsmodelling the ... more Numerical solutions of a coupled system of nonlinear partial differential equationsmodelling 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 timederivative are compared to solutions obtained with aRoe-Sweby scheme coupled to a minmod limiter and a TVD approximation to the time derivative. Results fromboth 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.TheRoe-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.
Uploads