分享一个程序：方柱绕流(fortran)

sjtugxl

有一段我还没看懂，想请教一下大家，可不可以给讲讲是什么意思？
公式的出处在哪里？需要阅读什么论文和资料？谢谢。

程序段在下面，为什么求P(x,y)是这样的？
          IF (wall(x,y)) THEN
            U(x,y) = 0.d0
            V(x,y) = 0.d0
            P(x,y) = density * cs
          ELSE
            Dxy=f(0,x,y)+f(1,x,y)+f(2,x,y)+f(3,x,y)+f(4,x,y)
     &                  +f(5,x,y)+f(6,x,y)+f(7,x,y)+f(8,x,y)
            U(x,y) = (f(1,x,y) + f(5,x,y) + f(8,x,y)
     &              -(f(3,x,y) + f(6,x,y) + f(7,x,y))) / Dxy
            V(x,y) = (f(2,x,y) + f(5,x,y) + f(6,x,y)
     &              -(f(4,x,y) + f(7,x,y) + f(8,x,y))) / Dxy
            P(x,y) = 2.0d0*(Dxy/Dxy0-1.0d0) * cs/Uw/Uw
          END IF

[ 本帖最后由 sjtugxl 于 2011-5-4 12:33 编辑 ]

feixiang9

这里应该是计算了一个使用动压头即无量纲化之后的压力差吧！
计算公式是：
p=cs*(rho-rho_0)/(1/2*(rho_0)*(u_0)^2)

sjtugxl

请问为什么要做这样的处理呢？

rho/3 不可以吗？

tsfed

rho/3只是方柱表面的压力（常数），而流场其他地方格点处的rho不是常数。

onesupeng

2楼的理解正确

tsfed

涡函数是 w=偏v/偏x-偏u/偏y  ，为什么还要除以x方向的长度xmax？

onesupeng

原帖由 tsfed 于 2011-5-5 02:57 发表

                               
登录/注册后可看大图

涡函数是 w=偏v/偏x-偏u/偏y  ，为什么还要除以x方向的长度xmax？

这是没有意义的处理，本来应该是想无量纲化

sjtugxl

rho不是随着迁移变化的吗？rho(x,y)不是常数吧？

为什么方柱表面的压力就是常数？

原帖由 tsfed 于 2011-5-5 09:57 发表

                               
登录/注册后可看大图

rho/3只是方柱表面的压力（常数），而流场其他地方格点处的rho不是常数。

sjtugxl

请教，为什么要求压力差呢？

原帖由 onesupeng 于 2011-5-5 10:52 发表

                               
登录/注册后可看大图

2楼的理解正确

tsfed

谢谢

onesupeng

原帖由 sjtugxl 于 2011-5-5 04:17 发表

                               
登录/注册后可看大图

请教，为什么要求压力差呢？

不可压缩流动里面的压力是相对的

sjtugxl

定义方柱的时候，程序中是这样的。是不是应该改一下？

1.         a=21
   
2.         b=21       
   
3.         x1=50
   
4.         x2=x1+a
   
5.         y1=(Ymax+1-b)/2
   
6.         y2=(Ymax+1+b)/2
   
7.       DO x = x1, x2
   
8.            wall(x,y1) = .true.
   
9.          wall(x,y2) = .true.
   
10.       END DO
    
11.       DO y = y1, y2
    
12.            wall(x1,y) = .true.
    
13.            wall(x2,y) = .true.
    
14.         END DO

复制代码

1.         a=21
   
2.         b=21       
   
3.         x1=50
   
4.         x2=x1+a
   
5.         y1=(Ymax+1-b)/2
   
6.         y2=(Ymax+1+b)/2
   
7.         DO x=x1,x2
   
8.                 DO y=y1,y2
   
9.                         wall(x,y)=.true.
   
10.                 END DO
    
11.         END DO

复制代码

onesupeng

应该两种都可以

前者是一个框，你是一个实体

sjtugxl

为什么方柱表面压力为常数？
那求解阻力的时候怎么求？

原帖由 tsfed 于 2011-5-5 09:57 发表

                               
登录/注册后可看大图

rho/3只是方柱表面的压力（常数），而流场其他地方格点处的rho不是常数。

onesupeng

表面压力那个是没有道理的，忽略它