A computational fluid dynamics code was employed to obtain detailed characteristics of horizontal wavy stratified two-phase flow for two geometries (pipe and channel). By treating the gas and the liquid phases separately, the calculations are performed utilizing information concerning the pressure drop and the liquid film thickness, which can either be obtained experimentally or calculated using one of the available correlations. The basic convergence criterion is the continuity of the time-averaged interfacial velocity and shear stress between the two phases. A commercial CFD code (CFX(R)) was used for the calculation of the velocity profiles of both phases and the distribution of the shear stresses, i.e. liquid-to-all, interfacial. The calculated values are found to be in good agreement with those obtained experimentally.