資源描述：
《流函數(shù)-渦量法的二維方腔流數(shù)值模擬》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。
1、流函數(shù)-渦量法的二維方腔流數(shù)值模擬基本方程：t———時(shí)間步長,s;γ———流體的運(yùn)動(dòng)粘度,m2/s;———渦量,s-1;———流函數(shù),m2/s;v———速度矢量,m/s;差分格式：采用FTCS格式，對于有：采用FTCS格式，對于有：邊界條件：速度的邊界條件：固體壁面處設(shè)置無滑移邊界條件,即u=0,v=0;流函數(shù)的邊界條件：在固體壁面及平板驅(qū)動(dòng)處的流函數(shù)=0;渦量的邊界條件：平板驅(qū)動(dòng)處：左右壁面上:;網(wǎng)格劃分：采用等距結(jié)構(gòu)化網(wǎng)格劃分（40*40）編程計(jì)算：本算例采用MATLAB進(jìn)行編譯，其主要優(yōu)勢是語言簡單，可以方便地描繪出方腔環(huán)流的等值線圖等。主要語句：whilenorm(c1
2、)>1e-4
3、norm(c2)>1e-4n=n+1;O1=O;E1=E;t=t+dt;fori=2:Iforj=2:JO(i,j)=O(i,j)-dt*(u(i,j)*(O(i+1,j)-O(i-1,j))/(2*dx)+v(i,j)*(O(i,j+1)-O(i,j-1))/(2*dy))+dt/re*((O(i+1,j)-2*O(i,j)+O(i-1,j))/dx^2+(O(i,j+1)-2*O(i,j)+O(i,j-1))/dy^2);endendfori=2:Iforj=2:JE(i,j)=E(i,j)+0.25*(E(i+1,j)+E(i-1,j)+E(i,j+1)+E
4、(i,j-1)-4*E(i,j)+dx^2*O(i,j));endend%0為渦量，E流函數(shù)%計(jì)算效果圖：分別設(shè)置雷諾數(shù)為200，500,1000雷諾數(shù)Re=200時(shí)的流函數(shù)圖及速度矢量圖雷諾數(shù)Re=500時(shí)的流函數(shù)圖及速度矢量圖雷諾數(shù)Re=1000時(shí)的流函數(shù)圖及速度矢量圖
