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 !