Fortran语言——涡量流函数法中心差分格式的二维方腔顶盖驱动计算要点详解.docxVIP

 PAGE \* MERGEFORMAT 23 涡量流函数法中心差分格式的二维方腔顶盖驱动计算 本题是关于粘性流体方腔顶盖驱动的问题。采用涡量流函数法，在均分网格下，用中心差分格式进行计算，结果与文献中所采用的其他方法和格式进行比较，认为中心差分格式符合计算的精度，但是其明显缺点是计算过程不稳定。 1、参数的无量纲化 令 X=xH; Y=yH; U=uut ; V=vut ; Ω=ww0 ; Ψ=φφ0 其中，由涡量的定义式： QUOTE X=xH; Y=yH; U=uut ; V=vut ; Ω=wwt;  将速度的导数项无量纲化，并将上述定义式代入，即 ω=?u?y-?v?x=ut?UH?Y-ut?VH?X=utH(?U?Y-?V?X) 令utH=w0 则涡量的无量纲形式可以写为： Ω=?U?Y-?V?X=ωω0 再由流函数的定义式：u=?ψ/?y 同理，令 ψ0=ut*H 可以得到流函数的无量纲表达式为：Ψ=ψψ0 流函数边界条件 X,Y∈AB Ψ=0 X,Y∈BC Ψ=0 X,Y∈CD Ψ=0 X,Y∈AD Ψ=0 涡量的边界条件（Thom公式） X,Y∈AB Ω1，j=2(Ψ2,j-Ψ1,j)δX2 X,Y∈BC Ωi，1=2(Ψi,2-Ψi,1)δY2 X,Y∈CD ΩL1,j=2(ΨL2,j-ΨL1,j)δX2 X,Y∈AD Ωi，M1=2(Ψi,M2-Ψi,M1)δY2 + 2 δY 2、方程的无量纲处理过程 2.1、非守恒型方程的处理 （1）将以上假设的各式代入到非守恒型方程组： 得到以下无量纲形式的方程： ①涡量控制方程 ?2Ω?X2+?2Ω?Y2=ReU?Ω?X+V?Ω?Y ②流函数控制方程 ?2Ψ?X2+?2Ψ?Y2-Ω=0 （2）采用有限差分法离散非守恒型涡量的无量纲控制方程： ReUPΩE-ΩW2ΔX+VPΩN-ΩS2ΔY=ΩW-2ΩP+ΩE?X2+ΩS-2ΩP+ΩN?Y2 整理得到：2?X2+2?Y2ΩP=1?X2+ReUP2ΔXΩW+1?X2-ReUP2ΔXΩE +1?Y2+ReVP2ΔYΩS+1?Y2-ReVP2ΔYΩN 简化为：4ΩP=1+ReUPΔX2ΩW+1-ReUPΔX2ΩE +1+ReVPΔY2ΩS+1-ReVPΔY2ΩN 写成一般形式是: apΩp=aWΩW+aEΩE+aSΩS+aNΩN 其中 ap=4 ；aW=1+ReUPΔX2 ；aE=1-ReUPΔX2 aS=1+ReVPΔY2 ； aN=1-ReVPΔY2 2.2、守恒型方程的处理 （1）将假设的各个无量纲量代入守恒型涡量方程 得到其无量纲形式为： Re[?UΩ?X+?VΩ?Y]=?2Ω?X2+?2Ω?Y2 （2）采用有限容积法离散守恒型涡量的无量纲控制方程： 对于扩散项，有 snwe[ ??X?Ω?X+??Y?Ω?Y]dXdY =?Ω?X|we?Y+?Ω?Y|sn?X =[ΩE-ΩPδXe-ΩP-Ωw(δX)w]?Y+ΩN-ΩPδYn-ΩP-ΩS(δY)s?X 对于对流项，有 snwe?UΩ?X+?VΩ?YdXdY=[(UΩ)e-(UΩ)w]?Y+[(VΩ)n-(VΩ)s] ?X 对于界面上的取值，采用如下形式： Ue≥0 时，Ωe=ΩP+(fe+-ΩP)； Ue≤0 时，Ωe=ΩE+(fe--ΩE) Uw≥0时，Ωw=ΩW+(fe+-ΩW)； Uw≤0 时，Ωw=ΩP+(fe--ΩP) Vn≥0时，Ωn=ΩP+(fe+-ΩP)； Vn≤0 时，Ωn=ΩN+(fe--ΩN) Vs≥0时，Ωs=ΩS+(fe+-ΩS)； Vs≤0 时，Ωs=ΩP+(fe--ΩP) 其中，界面上涡量的插值采用中心差分格式，令 fe+=fe-=ΩP+ΩE2 ；fw+=fw-=ΩP+Ωw2； fn+=fn-=ΩP+ΩN2 ；fs+=fs-=ΩP+ΩS2 所以，Ωe=ΩP+ΩE2 Ωw=ΩP+ΩW2 Ωn=ΩP+ΩN2 Ωs=ΩP+ΩS2 关于界面上速度的插值，可以用上一层次的流函数来表示，形式如下： Ue=Ψ