Page 275 - Computational Fluid Dynamics for Engineers
P. 275

9.2  Basic  Concepts  in  Grid  Generation  and  Mapping              265






                                         (U)                   (U)
              (U)
                                          D                  C


                                          A                  B

                               d,i)     (1,1)                 (1.1)

                               X     °    '     2      3       ^
         (a)                         (b)
         Fig.  9.1.  Mapping  an  irregular  simply  connected  region  into  the  computational  domain
         as  a  rectangle,  (a)  Physical  plane,  (b)  computational  plane.


         placed  along  the  bottom  boundary  segment  AB  of  the  physical  domain,  I-grid
         points  must  be  placed  along  the  opposing  top  boundary  segment  DC  of  the
         physical  domain.  Similarly,  if J-grid  points  are  placed  along  the  right  boundary
        segment  BC  of  the  physical  domain,  J-grid  points  must  also  be  placed  along
        the  left  segment  AD. The  actual  values  of  £ and  rj in the  computational  domain
         are  immaterial  because  they  do  not  appear  in  the  final  expressions.  Thus,  the
        coordinates  of the grid  A in the computational  domain  can be taken  as £  =  1 and
        7]  =  1 and  the  mesh  size  as  A£  =  1 and  Ar\  =  1. In  this  way,  a  square  mesh  can
         be  constructed  over  the  rectangular  transformed  region  in  the  computational
        plane.
           In  the  above  example  of  mapping,  an  irregular  region  in  the  physical  plane
        is mapped  as  a  rectangular  region  into  the  computational  plane.  Depending  on
        the  choice  of  the  values  of  £,  77 along  the  boundary  segments  of  the  physical
        region,  other  configurations  can  be  generated  in  the  computational  plane.  To
        illustrate  this,  let  us  consider  an  L-shaped  irregular  region  ABCDEFA  in  the
        physical  plane,  as  shown  in  Fig.  9.2a.  One  possibility  is to  map  the  region  into
         an  L-shaped  regular  region,  as  shown  in  Fig.  9.2b.  This  is  accomplished  by
        choosing the  values  of  77  along  various  boundary  segments  in the  physical  plane
        as  follows:

         1.  With  77  =  constant,  vary  £ monotonically  along the  boundary  segments  AB,
            DC  and  FE  of the  physical  region,  and
         2.  With  £ =  constant,  vary  77 monotonically  along  the  boundary  segments  FA,
            ED  and  CB  of the  physical  region.
        With  these  requirements,  the  L-shaped  irregular  region  can  be mapped  into  the
        computational  plane  as an  L-shaped  regular  region. The  coordinates  of the  grid
        point  A  in  the  computational  plane  can  again  be  chosen  as  £ =  77  =  1 and  the
        mesh  size taken  as  A£  =  Arj  =  1.
   270   271   272   273   274   275   276   277   278   279   280