Page 98 - Computational Fluid Dynamics for Engineers
P. 98

84                                                    3.  Turbulence  Models



                                    (£m)o  =  aUe<5*7tr7                    (3.2.3)

         Here  <5* denotes  the  displacement  thickness
                                                 u
                                             1       dy                    (3.2.4a)
                                      -j
                                        Jo
         and  7  accounts  for  the  intermittency  of the  outer  region  and  is represented  by
                                                [y-Y]
                                  7  =   1 - e r f                        (3.2.4b)
                                                 y/2a
         where  Y  and  a  are  general  intermittency  parameters,  with  Y  denoting  the
         value  of  y  for  which  7  =  0.5,  and  a  the  standard  deviation.  The  dimensionless
         intermittency  parameters  Y/6*  and  a/6*  are  expressed  as  functions  of  H  as
         shown  in  Fig.  3.1a.  The  variation  of  the  ratio  of  boundary-layer  thickness  6  to
         <5* with  H  is shown  in  Fig.  3.1b.  The  parameter  a  is  calculated  from

                                            °-° 168                        (3.2.5)
                                a  =                   1.5
                                     1
                                       -'(£'£
        Here  subscript  m  denotes  the  location  where  turbulent  shear  is maximum.  The
        parameter  f3 is given  by
                                          6
                                                     Rt  <  1.0
                                   l  +   2Rt(2-Rt)
                            0=<                                            (3.2.6)
                                   l  +  Rt
                                                     Rt  >  1-0
                                     Rt
        Here  Rt  denotes the ratio  of wall shear to maximum  Reynolds shear  stress  [2,3],

                                    Rt  =       = —                        (3.2.7)
                                          (-Q(u'v'))m
           In Eqs.  (3.2.1)  7t r  is an intermittency  factor  which  represents the  streamwise
        region from the onset  of transition to turbulent  flow. It  is defined  by the  following
        expression
                                                       rx  dx
                            7tr   1  — exp   -G(x-x tr)  I                 (3.2.*
                                                          u e
                                                      JXu
        where  xt Y  is the  location  of the  start  of the  transition  and  the  factor  G  is  given
        empirically  by
                                     G        R M                           3 2 9a
                                       = &$ £                              ( - - )
        with  R XtT  denoting  Reynolds  number,  R Xtr  =  (u ex/v) tr,  and  C  a  constant  with
         a recommended  value  of 60 for  high Reynolds  flows.  For  lower Reynolds  number,
        C  is given  by  [2, 3]
                                C 2  -  213(log Rx tr  -  4.7323)         (3.2.9b)
   93   94   95   96   97   98   99   100   101   102   103