一点疑惑关于Blazek的程序

qiongquangu

1. void BndConds::ZeroResiduals( const Geometry &geometry, Equations equsType, CONSVARS rhs[] )
   
2. {
   
3.   int i, ib, ibn, ibegn, iendn;
   
4. 
   
5.   ibegn = 0;
   
6. 
   
7.   for (ib=0; ib<geometry.nSegs; ib++)
   
8.   {
   
9.     iendn = geometry.ibound[ib].bnodeIndex;
   
10. 
    
11.     // symmetry boundary
    
12. 
    
13.     if (geometry.btype[ib]>=500 && geometry.btype[ib]<600)
    
14.     {
    
15.       if (geometry.btype[ib]-500 < 2) // x=const. line -> x-component
    
16.       {
    
17.         for (ibn=ibegn; ibn<=iendn; ibn++)
    
18.         {
    
19.           i = geometry.bnode[ibn].node;
    
20.           rhs[i].xmom = 0.0;
    
21.         }
    
22.       }
    
23.       else                            // y=const. line -> y-component
    
24.       {
    
25.         for (ibn=ibegn; ibn<=iendn; ibn++)
    
26.         {
    
27.           i = geometry.bnode[ibn].node;
    
28.           rhs[i].ymom = 0.0;
    
29.         }
    
30.       }
    
31. 
    
32.     // viscous (no-slip) wall
    
33.     }
    
34.     else if ((geometry.btype[ib]>=300 && geometry.btype[ib]<400) &&
    
35.              equsType==Equations::NavierStokes)
    
36.     {
    
37.       for (ibn=ibegn; ibn<=iendn; ibn++)
    
38.       {
    
39.         i = geometry.bnode[ibn].node;
    
40.         rhs[i].xmom = 0.0;       // velocity components = 0
    
41.         rhs[i].ymom = 0.0;
    
42.       }
    
43.     }
    
44.     ibegn = iendn + 1;
    
45.   } // boundary
    
46. }

复制代码

为什么要对粘性壁面的残差进行归零？ 粘性壁面也是有动量通量和粘性通量的，但是Blazek附带的C++程序，将其通量残差归零，不明白怎么回事，希望高手指点

qiongquangu

已经搞明白，因为要保证粘性壁面的速度为零，那么必须让壁面上的控制体的动量通量残差（包括对流通量和粘性通量）为0，否则，因为通量的存在，会使得粘性壁面上的速度不为0.不满足无滑移壁面条件