关于非正交网格修正的问题

dragonercfd

关于非正交网格修正的一些问题
1、请问在peric的fluxuv子程序中，fdue=visi*(2.*duxi*sx+(duyi+dvxi)*sy)这个扩散项的显式通量表达式是如何推导出来的啊？
我是根据附件中公式来推导的：

扩散项的显式通量表达式就是old中的第一项，怎么推导也应该是fdue=visi*(duxi*sx+duyi*sy)啊，为什么会多出dvxi来呢？
而在peric的fluxsc子程序中，扩散项的显式通量表达式却为：
fdfie=visi*rpr(ifi)*(dfxi*sx+dfyi*sy)，rpr(ifi)是什么东东啊？

harry

对于层流，你说得对，可以把剪切应力张量忽略掉，由于是非正交，peric的程序在处理时没有忽略，他把源项放到这里了而已
对于第二个问题，我在程序中没找到呀，你看peric的程序要用叉乘或点乘想想，一般都想得通
我现在没看这个程序了，搞自己课题了，原来就看了也不很明白，有机会交流

doowh

rpr(ifi)应该是Prandtl数的倒数，＝1/Pr

dragonercfd

楼上的兄弟，非常谢谢，第二个问题我现在明白了，这是标量方程中的广义扩散系数的普朗特数的倒数；第一个问题，是不是因为考虑到轴对称的效果呢？因为在柱坐标系下的速度梯度的计算会多出一些项来，不知道这种理解是否正确？