Hi Andy & Agg,

Thank you for your help. As per the first method suggested by Andy, where in we impose a constant mean shear stress and therby a constant Pr gradient,(correct me) so in that case for a channel flow, we calculate the shear stess based on the frictional velocity(deduced from frictional Re) and then use the realtion Tau = hydraulic perimeter * dP/dX

to calcualte the dP/dX and add them as source term into the Streamwise momentum equation.

The second method where in mass flow is constant and wall stress and/or pressure drop is adjusted iteratively to maintain an instantaeous force balance, is something new to me, can you kindly shed some more light on this issue.

Agg, thanks a lot for the Reference, I had searched a lot for the Literature u have quoted, I was unable to find it, if possible can you mail me a copy of this paper or can you quote some other references in this regard.

Thanks Lots,

Glen

[siyamak]

Hello,

Help please!

I wrote a 2D finite volume code for a rectangular using Neumann boundary in bottom, Dirichlet boundary at top and periodic boundary in x direction. I get some results and everything seems to be fine other than by results tilted toward right boundary and I'm assuming that is because of NOT applying correct periodic boundary condition. Could you please let me know how to build a coefficient matrix which have periodic boundary in x direction?

I appreciate if anybody can help me with that.