[Federico]

Hello everyone,

I am a graduate student working on annular two-phase flow. I am trying to implement my own Low Reynolds k-epsilon code in Fortran 90 for both, single-phase and two-phase flow, but I cannot get it right. I have already linearized the equations in order to have always positive variables along iterations, the problem seems to be in the wall boundary condition. I know it is a difficult task but as it is such a wide spread model I assume it cannot be that hard. I would really appreciate some advice. Thank you very much in advance,

Federico