Page 84 - Basic Structured Grid Generation
P. 84

Differential geometry of surfaces in E 3  73

                        Exercise 17. Making use of eqn (3.26), show that if A is a surface vector with covari-
                                                                 ˜
                        ant components given by the 2×1 column vector A (with elements A α ) and background
                        cartesian components given by the 3 × 1 column vector A (with elements A i ), then
                                                          √
                                                       ˜
                                                      A =   aC −1 A                       (3.150)
                        or, equivalently,
                                                                  α
                                                     A α = A i ∂y i /∂u .                 (3.151)
                        Exercise 18. Show by direct matrix multiplication using eqns (3.143) and (3.147) that

                                                     1   T            T

                                                 C  √ J     = I 3 − NN ,                  (3.152)
                                                      a
                        where here N stands for the column vector of cartesian components of the surface
                                                                               1   T
                        normal vector and I 3 is the unit 3 × 3 matrix, and hence that √ J  is not a right
                                                                                a
                        inverse for C.
                                                           1
                                                              2
                          Now consider a surface vector field V(u ,u ), defined at all points of the surface and
                        having the property that it is everywhere tangential to the surface. Then the divergence
                        of V is given by an expression analogous to eqn (1.134):
                                                   ∂V        ∂V       ∂V
                                                α         1         2
                                        ∇· V = a ·     = a ·     + a ·
                                                   ∂u α      ∂u 1     ∂u 2
                                                1            ∂V             ∂V
                                             = √    (a 2 × N) ·  + (N × a 1 ) ·           (3.153)
                                                 a           ∂u 1           ∂u 2
                                                1     ∂V i
                                             = √ C iα   α                                 (3.154)
                                                 a    ∂u
                        in terms of the background cartesian components V i of V. (Here the index i is summed
                        from 1 to 3, while α is summed from 1 to 2.) These expressions are non-conservative.
                        A conservative form follows by using eqn (3.144):

                                        1     ∂                ∂
                                ∇· V = √       ((a 2 × N) · V) +  ((N × a 1 ) · V) − 2κ m N · V,
                                         a ∂u 1               ∂u 2
                        but now the last term vanishes because V is a tangential vector. Hence we have
                                             1     ∂                ∂
                                     ∇· V = √       ((a 2 × N) · V) +  ((N × a 1 ) · V)   (3.155)
                                              a ∂u 1               ∂u 2
                                             1   ∂
                                          = √      (C iα V i ).                           (3.156)
                                              a ∂u α
                                                      1
                                                         2
                          For any surface scalar field ϕ(u ,u ), the gradient ∇ϕ must be a surface vector
                        field, according to eqn (3.139). Thus we can combine the results of eqns (3.141) and
                                                               2
                        (3.156) to give a formula for the Laplacian ∇ ϕ =∇ · ∇ϕ:
                                                     1  ∂     1       ∂ϕ
                                               2
                                             ∇ ϕ = √        √ C iα C iβ    .              (3.157)
                                                     a ∂u α   a       ∂u β
   79   80   81   82   83   84   85   86   87   88   89