流函数-涡量法的二维方腔流数值模拟(共13页).doc

精选优质文档-倾情为你奉上流函数- 涡量法的二维方腔流数值模拟基本方程：1 在直角坐标系下，不可压非定常流体所满足的流函数涡量形式的N-S方程为其中为雷诺数差分格式：采用FTCS格式有：对于本问题，将方腔四边同时分为等分，则有故2 在直角坐标系下，不可压定常流体所满足的流函数涡量形式的N-S方程为其中为雷诺数差分格式：采用FTCS格式有：对于本问题，将方腔四边同时分为等分，则有，则有即边界条件：在腔体的两侧和顶边，（第二式由泰勒级数展开得到）在底边（第二式由泰勒级数展开得到）其中代表边界，代表与边界相邻的节点。而即Matlab程序为：1 不可压非定常流体clear;%参数设置Re=10; %雷诺数取10,100，500，1000L=1;