Page 59 - Basics of Fluid Mechanics and Introduction to Computational Fluid Dynamics
P. 59

44       BASICS  OF  FLUID  MECHANICS  AND  INTRODUCTION  TO  CFD


              a  nonlinear  equation  written  in  conservative  form,  Le.,  in  a  domain  D

             of  the  plane  (z,t),  namely


                                                  ur  +  (f  (u)),  =0


             Or


                                                       divF  = 0

             where  F  =  (f  (uw)  ,u)  and  “div”  is  the  space-time  divergence  operator.
             If  ®  is  a  smooth  function  with  compact  support  in  the  plane  (z,t),  then
             the  above  differential  equation  leads  to  the  fulfilment,  for  any  ®  ,  of  the
              “orthogonality”  relation  {  @divFdzdt  =  0  which  comes,  by  integrating
                                              D
             by  parts,  to  [  grad®  -Fdzdt  =0.
                              D
                 If  wu  is  a  smooth  function  the  last  relation  is  eguivalent  with  the  given

             differential  equation;  but  if  it  is  not  smooth  enough,  the  last  equality
             keeps  its  sense  while  the  differential  equation  does  not.
                 We  will  say  that  u  is  a  weak  solution  of  the  differential  equation  if  it
             satisfies  {  grad®  -  Fdzdt  =  0  ,  for  any  smooth  function  ®  with  compact
                        D
              support.  Obviously,  if we  want  to  join  also  the  initial  conditions  u(z,0)  =
              uo,  then,  integrating  on  D,  (a  subdomain  of  D  from  the  half-plane  t¢  >  0)
              we  get

                               | sraas  -  Fdzdt  +  /  ®  (x, 0)  up  (x)  dz  =  0

                              Dt                        DNOz

              and  if  ®  has  its  support  far  from  the  real  axis,  the  last  term  would
             disappear  again.
                 So  we  have  both  a  differential        and  a  weak  form  for  the  considered
             equation.  We  will  also  have  an  integral  form  if  we  integrate  the  initial
                                                                                                  b
             equation  along  an  interval  [a,b]  of  the  real  axis,  precisely  a  fudz  =
                                                                                                 a
              f  (u)|®.

                 Of  course  we  should  ask  if  a  weak  solution  satisfies  necessarily  the  in-
             tegral  form  of  the  equation  ?  Provided  that  the  same  quantities,  which
             showing  up  in  the  conservative  form  of  the  equation  are  kept  for  the  in-
             tegral  form  too,  the  answer  is  affirmative.  That  is  why  the  weak  solution
             will  be  basically  the  target  of  our  searches.
                 Let  us  now  investigate  the  properties  of  the  weak  solutions  of  the
              conservation  law  uz  +  (f  (u)),,  =  0  in  the  neighborhood  of  a  jump  dis-

              continuity  (i.e.,  of  first  order,  the  only  ones  with  physical  sense).  Let  u
              be  a  weak  solution  along  the  smooth  curve  ©  in  the  plane  (z,t).  Let  ®
   54   55   56   57   58   59   60   61   62   63   64