EULER方程求解运动翼型的问题

dyb1984

各位达人，帮帮忙吧，我要毕不了业了快。我用ROE,FVS,AUSM写的求解二维翼型的euler程序，定常计算时不同格式算出来的结果吻合的还行。当我改成非定常程序时并计算NACA0012翼型的CT5算例时，发现计算的不是很好。我现在有几个疑问希望大家能帮我看看（我上传了我的程序，希望大家帮我看看那里不对吧）：
1.用MUSCL方法插值单元界面左右两侧的变量时，使用单元中心点上的速度时是否需要考虑物面运动速度。即使用u还是u-uw。
2.用MUSCL方法插值时，如果遇到物面边界条件，则使用的是物面上的变量值，还是虚拟网格上的变量值。比如如下图所示：

当计算红色界面左侧的变量时，使用MUSCL插值应到底是使用（界面，网格2，网格3）还是（网格1，网格2， 网格3）的变量。
3.物面边界条件的虚拟网格1我存放的是物面上的值，虚拟网格0的值等于虚拟网格1处的值，即：

！PW由网格点2,3,4的压力及其距离物面的距离的插值得到
PW    = PN2*CS2 + PN3*CS3 + PN4*CS4  
RHOW  = ( PW / (PN2/RHON2**GAMMA) )**(1.D0/GAMMA)

！ S(IM,JM,KM,1:3,IDIR)分别是物面单位法向量的三个分量
！ Q(IN,JN,KN,:)是网格点2上的主变量
UBAR  = Q(IN,JN,KN,2) * S(IM,JM,KM,1,IDIR)                &
           + Q(IN,JN,KN,3) * S(IM,JM,KM,2,IDIR)                          &
           + Q(IN,JN,KN,4) * S(IM,JM,KM,3,IDIR)                &
           - VS(IM,JM,KM,IDIR) / S(IM,JM,KM,4,IDIR)
！
Q(IBC ,JBC ,KBC ,1) = RHOW
Q(IBC ,JBC ,KBC ,2) = Q(IN,JN,KN,2)-UBAR*S(IM,JM,KM,1,IDIR)
Q(IBC ,JBC ,KBC ,3) = Q(IN,JN,KN,3)-UBAR*S(IM,JM,KM,2,IDIR)
Q(IBC ,JBC ,KBC ,4) = Q(IN,JN,KN,4)-UBAR*S(IM,JM,KM,3,IDIR)
Q(IBC ,JBC ,KBC ,5) = PW
！
Q(IBC2,JBC2,KBC2,1) = Q(IBC ,JBC ,KBC ,1)
Q(IBC2,JBC2,KBC2,2) = Q(IBC ,JBC ,KBC ,2)
Q(IBC2,JBC2,KBC2,3) = Q(IBC ,JBC ,KBC ,3)
Q(IBC2,JBC2,KBC2,4) = Q(IBC ,JBC ,KBC ,4)
Q(IBC2,JBC2,KBC2,5) = Q(IBC ,JBC ,KBC ,5)
4.ROE格式中是如何考虑动网格影响的，我在单元左右界面通量的对流项和ROE平均变量
构成的对流项中考虑物面运动速度，即：
！UBARW为物面运动速度
UBARW  = VS(I)
UBARP  = S(I,1)*UP  + S(I,2)*VP  + S(I,3)*WP - UBARW
UBARM  = S(I,1)*UM + S(I,2)*VM + S(I,3)*WM - UBARW
UBARR  =  S(I,1)*UR + S(I,2)*VR + S(I,3)*WR - UBARW
下面是我计算的结果和别人计算的结果以及实验结果


[ 本帖最后由 dyb1984 于 2010-11-11 00:32 编辑 ]

dyb1984

跪求各位啦

dyb1984

顶

通流

算的挺好，怎么就毕不了业了？

onesupeng

用ISAAC？有说明书的吧？

onesupeng

原帖由 通流 于 2010-11-11 17:09 发表

                               
登录/注册后可看大图

算的挺好，怎么就毕不了业了？

算一算就可以毕业了？

dyb1984

我是在ISAAC的基础上该的，只有一个很简单的帮助文档那里有什么说明书啊。然后那个黑色的结果使用TUMFLMCODE算的，但是我没有源程序那个。

dyb1984

我现在换成VAN ALBADA限制器，感觉稍微好点。哎，很想找本有源程序的书看看大牛是怎么处理这种简单问题的，我也好有个对照。要不然按照自己的理解和计算结果验证感觉不是很严谨啊

dyb1984

我想问一下这里谁有这种非定长的euler程序，就是算物面运动的动网格程序，可以给我一份交流分享一下吗？还有那位达人有computational fluid dynamics: principles and applications 这本书的配套程序，小弟原用wilcox的turbulence model for CFD的源程序换