Page 199 - Computational Fluid Dynamics for Engineers
P. 199

6.3  Finite-Difference  Method                                        187



         and  with  fij  =  0

                                                     0
                                     F 2

                               F  =                                      (6.3.16b)
                                     F   >  Fj  =    0
                                     J
                                    Fj               0

         the  Aj,  Bj  and  Cj  matrices  of the  coefficient  matrix  A  in  Eq.  (4.5 •7)  are  given
         by
                                     a\  «2
                                        a\
                                     a 3    a 2
                                        as  a\  a?
                               Ai  =                                     (6.3.17a)

                                                  «3  a*  ct2



                                     Oil  CX2
                                     c*3  a\  a.2
                                        az  a\  0:2
                                  =                                      (6.3.17b)
                               A 3
                                                 « 3  d\  OL2


                                    1  al  a 2
                                     az  a\  a.2
                                        as  a\  a2
                               Aj  =                                     (6.3.17c)

                                                  «3  a\  a 2
                                                         at,
                                                     a* 2
                              Bj  =  Cj  =  Ii      2<j<J- 1             (6.3.17d)
                                         -H
                                      Ci          •II                    (6.3.17e)

         where
                                   1                4      +        4
              =  a\  +  - ,
                         a^  =  as  -  -OL2,  al  =  aj  +  - a 2 ,  a^  =  ai  +  - a 2  (6.3.18)
           a 1
           Table  6.1  and  Fig.  6.6  show  the  results  obtained  for  a  grid  composed  of  91
        points  in the  ^-direction  and  46 points  in the  //-direction  with  770 — 0.1.  Several
   194   195   196   197   198   199   200   201   202   203   204