operator splitting for convection-dominated nonlinear

0 downloads 0 Views 731KB Size Report
influx of modern software tools such as object-orientation into scientific computing. ... of small parts that can be fitted together independently (like Lego pieces), ... pecially to capture sharp gradient variations in the solutions of (1). The.
In "Godunov Methods: Theory & Applications", Edited Review, E. F. Toro (editor), Kluwer Academic/Plenum Publishers, 2000, pp. 469-475

         

         

 

                          ! " # $  

    

       % &  %$ ' ())* %   ! " + +# $ ,      

 -      .  / . % 0 )*) %  )12 .    ! " + # 3 $ 

      

       .  / . % 0 )(1 %  )12 .    ! "   # $ 

   

     Æ                            

 





 

 







  



  



                     ¼

       ! "                   #  $                      %              ! &                       

               '  





 



Æ  Æ 

 







 ' 

(       %               Æ                )           ! "             %            ! *                    ! *             Æ                 ! + %                #      )                               ! *    %        ,       

   - )     #    ! & - )      .       -               ! "   -      %    #    

      ! &   %       - )

       ! / %     

 

         #        . + %               .   0 "                 )    1               

       #              /   !% 222! "        Æ          1

   Æ      3             !               )                     ! "               ) )   )  . 

   

4

   /   /  % 2251 6  % 2241 7 % 22   

)                 ! "                                Æ ! *        #         ,       %    %       %  %  Æ , %  )   !

  ! "#  $% &             Æ            . %      # ! "  8 )  . 9      

             .     .             ! *          :  ! $       

 !&#

    



   

  

 

; .                 %  $            6  ! "                   6         % !!%             )            %   .      )  ! <          %  #   6            %   ! (       %      6    )      .          ,  )  ! &     %      .       )        ! "   .            )     #                     6  %  /   /  % 225! (       .           )      ,                           ! "   .                        Æ              !!%  )          # +% 222! "         .              $    /   6  % 224! =    .         %    3                ! &   %               1               '

'!(" 

   

> Initial Data

0.8

Flux Function

0.6

1

0.5 0.4

0.2

0.4

0.5

0

−0.2

0.3

−0.4

0

−0.6

0.2 −0.5

−0.8 −2

−1.5

−1

−0.5

0

0.5

1

1.5

2

−1.5

−1

−0.5

0

0.5

1

1.5

2

−1.5

−1

−0.5

0

0.5

1

1.5

2

0.8

0.1

0.6

−1 −2

0.4

0 −1.5

−1

−0.5

0

0.5

1

1.5

2

−1

−0.5

0

0.5

0.2

1

0

Solution

Fronts in (x,t)−plane

−0.2

1

1

−0.4

−0.6

0.8

−0.8 −2

0.5 0.8

0.6 0.6

0

0.4

0.4

0.2

0

−0.5 0.2

−0.2

−0.4

−1 −2

−0.6

−1.5

−1

−0.5

0

0.5

1

1.5

2

0 −2

−1.5

−1

−0.5

0

0.5

1

1.5

2

−0.8 −2

     ! "#$    "# %& '(#) * + , -. /#"  ''0 )    1  , 23 /#"     !+ 4)  5$0 ''0 ) 6  , %-7 237 1 -. ) '  ()  , 38.9

          -        ' %  ( ! "             +% / %  ?  % 225         %  +% 222! (      %         +% / %  ?  % 2251 +% 222     )                   -  ! (               $(+ ) ! ( 4       79           #           ! @        

 %          

! ! "     6                          .       ! " #  . ,      ! @     ,   .              )     %    .    ! "                        .!

'!(" 

   

A

1

1

1

0

0

0

−1

−1

−1

0

1

−1

−1

0

1

−1

0

1

  # 1 (1:   $1 +  !   ;   , 8   :# 1  1 $" 8. 1  878- #$ "   )  ,  3 "# 17  1"+17 1 3 11 )'$1 * + %- #7 237 1 -. "'#  ! "'"  !+ 1

Suggest Documents