Annular flow region boundary conditions

The CNMMR model with laminar flow liquid stream in the annular region consists of three ordinary differential equations for the gas in the tube core and two partial differential equations for the liquid in the annular region. These equations are coupled through the diffusion-reaction equations inside the membrane and boundary conditions. The model can be solved by first discretizing the liquid-phase mass balance equations in the radial direction by the orthogonal collocation technique. The resulting equations are then solved by a semi-implicit integration procedure [Harold etal., 1989]. 

Calculate the stream function for axisymmetric fully developed creeping viscous flow of an incompressible Newtonian fluid in the annular region between two concentric tubes. This problem is analogous to axial flow on the shell side of a double-pipe heat exchanger. It is not necessary to solve algebraically for all the integration constants. However, you must include all the boundary conditions that allow one to determine a unique solution for i/f. Express your answer for the stream function in terms of ... 

The well-posedness boundary (ZRC) (included in Figures 10, 11, 13) represents the limit of operational conditions (U, U, ) for which the governing set of continuity and momentum equations is still well-posed with respect to all wave modes. Hence, it is considered as an upper bound for the stratified-wavy flow pattern. Indeed, the data of stratified-wavy/annular transition follows the ZRC curve in the region of H < 0.5. 

---