Page 109 - Introduction to Computational Fluid Dynamics
P. 109

P1: IWV
                                                                                                11:7
                                        0 521 85326 5
                                                                                   May 25, 2005
            0521853265c04
                           CB908/Date
                     88
                                                                               2D BOUNDARY LAYERS
                            are many variants, all turbulence models of this type stem from a dimensionally
                            correct representation
                                                           µ t ∝ ρ l v ,                       (4.90)


                            where v is the representative velocity fluctuation scale in the transverse direction
                            y and l is a representative length scale. Two turbulence models used extensively for
                            boundary layer calculations are described in the following.
                            4.8.1 Mixing Length Model


                            Since v is responsible for transverse momentum transfer, it may be written in
                            dimensionally correct form as

                                                                   ∂u

                                                          v = l m                              (4.91)
                                                                    ∂y
                            so that

                                                                    ∂u
                                                         µ t = ρ l  2       ,                  (4.92)
                                                                   ∂y
                                                                m
                            where l m is called Prandtl’s mixing length. Now, because the velocity gradient can
                            be evaluated from the solution of the momentum equation, l m must be prescribed to
                            complete evaluation of µ t . Kays and Crawford [33], after extensive investigations
                            of a variety of wall-boundary-layer flows have prescribed the following formulas:

                                                                 y            y
                                               ⎧                  +
                                                 κ y 1 − exp −         , for    < 0.2,         (4.93)
                                               ⎪
                                               ⎪
                                                                 A +          δ
                                               ⎨
                                          l m =
                                                               y
                                               ⎪
                                               ⎩ 0.085δ   for    ≥ 0.2,                        (4.94)
                                               ⎪
                                                               δ
                            where y is the normal distance from the wall, δ is the velocity-boundary-layer
                            thickness and κ = 0.41. Further,

                                                    yu τ         τ w         ∂u
                                               +
                                              y =      , u τ =      ,τ w = µ    | w .          (4.95)
                                                     ν            ρ          ∂y
                               Finally, the value of A is sensitised to effects of suction or blowing and local
                                                   +
                            pressure gradient in a generalised manner as
                                                                  
                −1
                                            +
                                                         +
                                           A = 25 a v + bp / 1 + c v     +    + 1    ,         (4.96)
                                                                +
                                                         w               w
                            where
                                                            3
                                                +
                                                                          +
                                               p = µ  dp 
 τ ρ    −0.5  ,  v =  v w ,          (4.97)
                                                            w
                                                                          w
                                                       dx                     u τ
                            and a = 7.1, b = 4.25, and, c = 10.0. If p > 0 then b = 2.9 and c = 0.
                                                                 +
   104   105   106   107   108   109   110   111   112   113   114