Numerical simulation of churn flow in a vertical pipe