涡量流函数法处理台阶流

sanqianchi

我用Fortran程序不是很熟，在用涡量流函数方法模拟台阶流时，初始场的涡量和流函数应怎样给定？流函数边界又怎么更新呢？还有入口和出口速度又是如何更新的呢？
   我的初始涡量边界是根据涡量定义给定的，初始速度场为u(x,y)=1，v(x,y)=0.上下边界采用固壁边界条件。
   希望有高人能指点一二，急啊....