Newton's Method for Steady and Unsteady Reacting ...

90 downloads 0 Views 2MB Size Report
ABSTRACT. This paper describes the application of Newton's method to low Mach number steady and unsteady laminar diffusion flames, which are ...
Newton’s Method for Steady and Unsteady Reacting Flows 

Wensheng Shen





Jun Zhang

Fuqian Yang

Department of Computer Science University of Kentucky Lexington, KY, 40506, USA

Department of Computer Science University of Kentucky Lexington, KY, 40408, USA

Department of Chemical and Materials Engineering University of Kentucky Lexington, KY, 40508, USA

[email protected]

[email protected]

[email protected]

ABSTRACT

                                                          

        !      "      "    "         " # $  %                                     !         " #$  %          &          

  "   $ '              "   ! (                 

 

          "       "        

   ! )      % "                

    "           !                                          !      *       +                  "              , " ( - #. ,   

-  " #   

   !

Categories and Subject Descriptors

- !/ 01 2 3 4 5 6 2 3 78 9 : ; < : 6 = > 3 7? @ A B C   D    ' %  



    E      ! F F G"H I  D     (     #    J     K  $ L +    K M  NH H I !        $             # O     ( ( E "H HP N IQ N  . ( E "H NH NP QF   ( ( O "H I N RP S R      D T '     D ' "O - H N"H N' E F IP S /      

K  $ #    '   

  O     !  -   #  !  C   C   !  .   C   !

Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. To copy otherwise, to republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. ACMSE ’06 Melbourne, Florida USA Copyright 200X ACM X-XXXXX-XX-X/XX/XX ... U 5.00.

General Terms      

Keywords

                  

1. INTRODUCTION

D        %         

           V                0/N  /R  /G  / /A ! .             W                   ! #     

                             W          

     ! )             

                        "       +                         "                                  !        %                                        0FA ! C "    +                    

+    W                   

%              0/N  /R  /G  / N  IA !                                          !                            0GA              " # $  %      "                 ,     #   $ 0NA       "             " #$    

    " &    !               "    "          " # $  %                        0NA        

          "               &    !   W              " 

          "    " # $    "              "      " # $ 

  !                 !  

                    #

"   N             #

  Q ! O "       +          #

  F !                #

  I             #

  S !

2.

LAMINAR DIFFUSION FLAME MODEL

                     B             !                 " # $  %                    "          

    0/N  /R  /G  FA !      

         !

2.1 Vorticity-Velocity Formulation

.                       O  $                 " #$  %    

                         !                "             " #$                           %                   " #$  %            & "         !        %            "                "      

+            0I  /HA B /  

 

 

  

  

/  

 

  





N





 



/





 





N



 



/ 



 

  



   







  



@

   



>

 

 

N

H

          

 

 





 

  

N



H

      



   



  





/



         

      







      





 

            



 





              

>



 

Q

 



    

         

 

    %   

   %      

  +              

    +          &              >  



  



 

 

 > 0 A     @  

      

   

 @ 0  A  !     

        

   +   

       

> 







      

)      

                      !     

         ,  $  #     0QA            0P  /N  /RA !     

            $                     "

                                         !                  

          +                    0/NA ! J     

                        "                     &            

+                  +   "             ! #    &     $                     +            +       + "        ! D    +                &         +     

        !                     +       

            0/GA ./ 0 /



2.2 Flame Sheet Model

 



N



   

/RA                +    ! "#$ %    ' ( )*          +  &     

          ,     

     -           !    

&               %           $        

     "               %        0IA       L   %      $     "         %   ! 

F

.1 01



.2 0 2 3 .4 0 4



I

.2 0 2  

  0 /  0 1  0 4   0 2            + "              

   . /  . 1  . 4   . 2          W      !                    "            

&    5          !       O  $       L         



       &    6  6  6  

7 8 9 :; 7 8 < :; 7 8 = :; 6     ?              +   7 8 > :; @ /  @ 1  @ 4   @ 2       W    

   +               %     !     %         %     

  

              &    0/RA ! .                                 

 

        !           % "     

  %             

                  "     

  "        %   0/N  /GA



        @                           !        % "                " # $   

%            +   !         %             &                            !               

              !                      +        0I  /N 



A 

 

A 



 

A 



/



   

 A   A @ @    

  

S     "

 

 "    %    +        

 

                               + "           !

2.3 Boundary Conditions

      +                

  " # $  %          %    +     +                 "

   

   +   

    

           

   

       !               " # $  % "          %           0F  SA B   +         HB H H     H         , B H      

H H       H  B   H     A            A         A B   ' + 

      



H  H  H ! H            





3.

H  

  ?     "          $    H  ?   / !   *       +  ' %  ! G   //                    0/RA ! O 

  

         "             " # $  %       $   

              ! ! > 0   A         A ! )             

          0      A  ! .                      %   $        F  F    > 

NEWTON’S METHOD

                      "                     &  

     0/F  /SA !       #  $  0/N  /G  I  N  /HA     

         "        !                              %           "          !            

    %         W                            %                     %   $ !                     

*      +    +         W    &          !    

    0/N  /RA                  *       +      V            !                 

                         ! O               "              !               %     

      "       +    ! L   >  H        %                            %        

+       "     >            >   0RA   >  



 >



 >  > 



>





> 





>  





R         

         >  

    $           ! ,           >  H   ! ! ' % ! R      

  

 

  "                

      





>   >  

 > 



>



>







>



 >     



H  /





>     >     



G

 

H  /





P

 

      *       +        

 "        !            

   %              B >   >  

 > 



>



? >





>

? 

>  



?



H  /



>     >     



/H 

 



H  /





// 



   

/   /    /    /    

/  /     /     /      

/  /     /     /      

/  /     /     /      

! ! ! "

/N  

'         F  F   $     

&       +    0/RA  

#





 #







# 

#



# 

 

/Q 

 

  #            "        #    $   

 >            #   $ %  

      $     %    $ %           "  !                 N"           $   

  

     

"       &    "         "

 !   

&            0IA                    

              %   "        N"      

 > &  

    0IA

' > & '



)/ &  & &  ,  ) , , *(+ # - -  -



/ $

>&

 

/F 

 

   +           "   "  . . . $         $   

       

. . .          " 

           " 

               $   

 > 

 ! O             

       ' > & '  /H  !         "

             K    

    

   !

4. FINITE DIFFERENCE APPROXIMATION AND LINEAR SOLVER                   &        " %    !     " #$  %        "          

    %    & "     "      

  "      

 "

    ! D           "         

     ! D    +    

 "        "  &               "      +               ! O   "                 &     

  %     

      Q              

            0NA !

,   & "       

  "               "

 ! J     "        +          &                  









  

 #





#  # 



 #



 #  #



 #   #

N

 

 #  #

/I 





#  # 



 #   #



  



  7

   

 



 

 

  

 7  

#







  7 





 

 

#



#

 

 



#



    7   



 

#



 

 # 





 





#

# ' 





#

N

#



'





 





/S 



 #

#





#



  #  H '* ( 



N



#



# #

 H '* ( 

#

/R 

#



   # 

 

 #  #







#







 #



 #  # 

/G  







  '   '* (        $     +        !   

  "                        "      







  # #               7  



  

#    #  #   #    #   #               7    #

# #   #     



 

/      









6

1

2

3



92 6 =  5



= : 7?3 @ 7  

             "       "           V         

      

  !     &       V               0G  SA !     V   

    +   V  "          V       ,  H N       +     V       , 1 N I  !         "              +                   H !I N !   +                 H NQ N ! ,       +             Q I    !          1     +      H   R I     

     

   

H 

Q H    +   

                                      +   V   , 1   

+                     

   0GA  !                 "            /NP  /S /! J   

    Q  Q     +        O   /     P               *       +      ' % ! N /  !

/P  

      7     7   7       7  

   



5. RESULTS AND DISCUSSION



O                     "           

    $ '                               ' % ! NH  ! O           H  

  

      

       ,  

   

   $       !       7     7    7       7 

5



  & "              

"                    "  +   0/N  /QA  

4

 7@ > 5



 

9







8

7

NH   .           %        

      %        +           %        + %     "           $                         !       "    %                   

          !       *       +  "                 " V   

      , " ( - # . ,  0/IA           $ -  " #    

   0IA      

    !   , " ( - # . ,       "       +       %       + "

   "           !





           





































































































! ! ! ! ! ! ! ! ! ! ! "

N / 

.      

  Q       "         A             +

   F F     $       ' % ! /N  ! L    .: .                        .      

          +   

"   

                   $           +             .: . . 

       P      +       .:         *      +  P     ! O .: . .      " /NP  /S /               $      GQ  H RS               N  P P H  RQS ! O      

       "    

                     0/  /GA    +   +  NN .1     3 . 1   . 1                   W   .    . 1  .  1       .  1   $   /  N  N   / 

 !  )    

             

%   V  

  ! J         &      "

. 



/

%

< 


5    4 5 35 6 = 5 2 3 > 5 = :   5 :  32 7? 5  ; : 6 3 4 5 935 2    932 35 5  > 2 3 7: ? 9 ;: 34 5  7  > 97: ?  2 6 5 

                 " 

     "                           A   NP G K  - !             

      %                       N"                  &    ' % ! /F        / H  /H    , " ( - #. ,   "            N"          

       %             0IA !     "            

"         0F  I  SA    &      /NP  /S /        Q Q  F /!     "       

  "          &      

                     

      "                    !     "               O  N        +    +    "       N !SQ     NH IH K 

      ( C J   #  ",  " /HH        IH H    # C. E (     N - ,      / /!G     !             IIP  /H P R   /IS G    /  N   Q 

      /   &        Q        !    "                   "       /NP  /S / !            "             O  Q                      &          ! .  H /N I          "     &             "                 NH IH K     +    "   N !S /    !       IH     ( C J 

 

     "            "   %               "  %          "     !     

    

     "  %       &            *       +           !

6. CONCLUSIONS

)                    "                ! O              "           " #$  %                  "  &  %    

      &    

10

9

9

9

9

9

8

8

8

8

8

7

7

7

7

7

6

6

6

6

6

5

5

5

4

4

4

4

3

3

3

3

3

2

2

2

2

2

1

1

1

1

-1

0 r(cm) 

1

0 -2

2



 7@ > 5   4 5  5 = 5 ?  5 ? 8  : ;  H / 9  2 ?  5   H /N I 9  



-1

0 r(cm) 

1

 26 5

2

0 -2

-1

1

0 -2

2

1 -1

0

r(cm)

1

0 -2

2

-1







H H N I 9    

1

2







0

r(cm)





H H I 9  8  

 



H H R I 9    



REFERENCES

0/A , ! ,      ! #  $ ! L  

    &            +            !         N Q  BN N / NIG  /PP G !  0NA , ! ,      ! #  $ ! J    +                B         

    &    !                              Q Q N  Q Q I ! (          NHH /! 0QA # ! ,  $   ! #     ! D       !               NH /H  BPP G  /HHF  /P NG ! 0FA ( ! D      . ! '     ! #  $ !                          ! 

   !      !!         F B/F P  /IF  /PP G ! 0IA . ! '  ! "   # "          $   % & ! C  D     M  J      

0 r(cm)



35 6 = 5 2 3 > 5 : ? 3 76 5  2 

     !        

       

          "     

%    ! '            %     "                     "        !                    W                                       "          !     "                "    " # $  %           "         , " ( - # . ,     

      - #         "    " #$  %           

 

  

             !            W         "       

  "    " # $  %                    " #$  %             "     !

7.

5

4

0 -2



5

z(cm)

407.50 626.50 845.50 1064.50 1283.49 1502.49 1721.49 1940.49

10

z(cm)

407.50 626.50 845.50 1064.49 1283.49 1502.49 1721.49 1940.49

10

z(cm)

407.50 626.50 845.50 1064.50 1283.50 1502.50 1721.50 1940.50

10

z(cm)

407.49 626.47 845.44 1064.42 1283.40 1502.38 1721.35 1940.33

10

z(cm)

407.50 626.50 845.50 1064.49 1283.49 1502.49 1721.49 1940.49



  (   

   /P PF !

0SA # ! K  * ! '      ( ! D      ! C

               +              !    !     Q G BNS P  NRP  NHH Q ! 0RA D ! K    D ! K   ! *   "

   "$  

    B               ! ( ! )    /PQ N  BQ I R  Q P R  *   NH HF ! 0GA E !        !      ! #   $    . ! #   ! (        +           "     +                !  &  # *   !    +    ,    -             S P Q  RH N !   (          /PP G ! 0PA * !      !             (  . /  %   ! C  D     J     )     #   )     NH H I !

 !  !       *         #      #. !   0      % & &   !!    1    . /  %       ! C  D     M 

0/HA  !   $ !

J     



  (   

    NH HQ !

0/ /A C ! C    L !     * ! ,   C ! (     ( ! (  &   ) ! O     * ! * 

! .  

 V

          "    

     !       /F H B/NQ  /S G  /P PG ! 0/NA  ! #   $  E !       D ! K   !            "      +                !       S R BG I  /N N  /P GP ! 0/QA - ! # ! 2        %   .    ! (     J    C     M $  /PG I ! 0/FA ( !  D    !      * ! .   ! (                              

  !  (     NG I  BPQ R  P QP  /PP H ! 0/IA  !     ! , "   B        

        "                    !   (      !  /Q N  BS Q / SF F  /P P N ! 0/SA  !   $  $    !                       !  (     N R R  BG G I  GP / /P GP ! 0/RA M ! 3  ! 2         0     

1    . /  %   &   .     $   $        ! C  D     M  J   

     (   

   /P P / ! 0/GA M ! 3    ! #   $ ! C                         !       P H BNG P  Q /Q  /PP Q !