懂peric非正交程序的帮个忙

wildgrass

在peric非正交网格流动求救程序中，采用高斯定律求节点梯度时， 在GRADCO
子程序中这个是什么意思，感觉是网格方向梯度与法线梯度之差，推半天也推不出来？
  FIE=FI(IJN)*FAC+FI(IJP)*FACP+DFXI*(XF-XI)+DFYI*(YF-YI)

harry

[这个贴子最后由harry在 2004/11/07 11:20am 第 1 次编辑]

哈哈，我也遇到了这个问题，斑竹CFDNHT给了我回答，下边是他的意见
是这样的，某个界面的值＝线形平均插值＋梯度修正值×距离 。不论是u，v，T 还是P，rho，都是这样处理的。你仔细看看peric的书《computational methods for flow dynamics》2002. 或其他peric的文章，这就是peric所提出的特殊的对非正交曲线界面值的求解公式。
--------------------------------------------------------------------------
下面是我的一点点评注
事实上在《computational methods for flow dynamics》，1996版P217提到了这种线性插值和中点估计准则，只不过这个方法在P299多重网格中才给予详细的解释，这下这个程序看懂一大部分了，基础比较差，费劲呀

harry

接着这个帖，我再问一个问题，在peric的k－e湍流计算中，他只保留了法向的切应力，其它的给忽略了，层流这样处理还可以，湍流好像不行吧，这样对计算的误差有多大呀，详细地说
对于速度U在peric程序中，应力txx留下用于计算法向的扩散流量，而txy却被忽略了，（常规的是把txy放在源项中，陶文铨的k－e程序中（subroutine 10）就有单独处理这个源项的子程序）
高人们提提意见呀

liuhuafeifei

Peric并没有忽略txy

harry

非常感谢指点，俺基础比较差，到处都是问题呀
确实txy没忽略，我回过头把peric的书的第七章又详细地看了一下，并对比了一下non－orthogonal和Cartesian的程序，终于明白了
在non－orthogonal程序的SUBROUTINE FLUXUV(IJP,IJN,IJ1,IJ2,FM,CAP,CAN,FAC,G)中，下边这个语句就是算源项的txy的
FDUE=VISI*(2.*DUXI*SX+(DUYI+DVXI)*SY)
我原来还以为这个语句是算网格线和法线方向的扩散对流的差，就是peric书上提到Muzaferija所提出的公式，所以怎么想也没想通，呵呵太笨了。

  

harry

[这个贴子最后由harry在 2004/11/09 08:45am 第 2 次编辑]

我想问的是，上边这个公式（即Muzaferijia提出的）在程序中用了吗？
另外一个问题是，对于peric的non-orthogonal的单层网格程序，明明说了不是多重网格，怎么还有粗网格和细网格之间转换呀！是不是这个在计算中没有用呀？

CFDNHT

本人曾在去年看过peric的程序和书，总结就是：书中几乎所有公式都在不同文件程序体现出来了，而且如果再多看看peric、demirdzic、muzaferija等几位大牛的论文（1986～1996），就会发觉程序丝路很容易懂。当然基本功就是要熟悉张量知识，能善于推导公式！－－－－我个人没看懂很多，更没用过。

凤翼天翔

这个贴子好。从中学到了不少东西。
我也在看Peric的程序，用它来做高雷诺数下的圆柱绕流。只是效果太不理想，是不是对流项格式的精度低了？
不知道user.f中求解总压力和总剪切力的语句有没有问题，可以用吗？
请指教，谢谢。

liuhuafeifei

下面引用由凤翼天翔在 2004/11/29 10:19am 发表的内容：
不知道user.f中求解总压力和总剪切力的语句有没有问题，可以用吗？
请 ...

没问题