10.1.1.35.2402 (1).ps
TRANSCRIPT
-
8/14/2019 10.1.1.35.2402 (1).ps
1/16
N u m e r i c a l S i m u l a t i o n o f 2 D S q u a r e D r i v e n C a v i t y U s i n g
F o u r t h O r d e r C o m p a c t F i n i t e D i e r e n c e S c h e m e s
J u n Z h a n g
y
D e p a r t m e n t o f C o m p u t e r S c i e n c e , U n i v e r s i t y o f K e n t u c k y ,
7 7 3 A n d e r s o n H a l l , L e x i n g t o n , K Y 4 0 5 0 6 - 0 0 4 6 , U S A
A p r i l 6 , 2 0 0 0
A b s t r a c t
F o u r t h o r d e r c o m p a c t n i t e d i e r e n c e s c h e m e s a r e e m p l o y e d w i t h m u l t i g r i d t e c h -
n i q u e s t o s i m u l a t e t h e t w o d i m e n s i o n a l s q u a r e d r i v e n c a v i t y o w w i t h s m a l l t o l a r g e
R e y n o l d s n u m b e r s . T h e g o v e r n i n g N a v i e r - S t o k e s e q u a t i o n i s l i n e a r i z e d i n s t r e a m f u n c -
t i o n a n d v o r t i c i t y f o r m u l a t i o n . T h e f o u r t h o r d e r c o m p a c t a p p r o x i m a t i o n s c h e m e s a r e
c o u p l e d w i t h f o u r t h o r d e r a p p r o x i m a t i o n s f o r v e l o c i t i e s a n d v o r t i c i t y b o u n d a r i e s . N u -
m e r i c a l s o l u t i o n s a r e o b t a i n e d f o r s q u a r e d r i v e n c a v i t y o w a t h i g h R e y n o l d s n u m b e r s
a n d a r e c o m p a r e d w i t h s o l u t i o n s o b t a i n e d b y o t h e r r e s e a r c h e r s u s i n g o t h e r a p p r o x i -
m a t i o n m e t h o d s .
K e y w o r d s : I n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s , s q u a r e d r i v e n c a v i t y p r o b l e m , c o n -
v e c t i o n d i u s i o n e q u a t i o n , m u l t i g r i d m e t h o d , f o u r t h o r d e r c o m p a c t d i s c r e t i z a t i o n .
M a t h e m a t i c s S u b j e c t C l a s s i c a t i o n : 6 5 N 0 6 , 6 5 N 2 2 , 6 5 F 1 0 .
1 I n t r o d u c t i o n
T h e N a v i e r - S t o k e s e q u a t i o n s h a v e b e e n u s e d t o m o d e l u i d d y n a m i c s p h e n o m e n a d e s c r i b -
i n g o w s o f a n i n c o m p r e s s i b l e v i s c o u s u i d . T h e s e e q u a t i o n s a r e h i g h l y n o n l i n e a r a n d a r e
v e r y d i c u l t t o s o l v e , e s p e c i a l l y w h e n t h e a p p r o x i m a t e s o l u t i o n s a r e r e q u i r e d t o h a v e a
h i g h a c c u r a c y . A N a v i e r - S t o k e s e q u a t i o n m a y b e l i n e a r i z e d i n s t r e a m f u n c t i o n a n d v o r t i c -
i t y f o r m u l a t i o n 6 , 1 8 ] . C o m p u t i n g a c c u r a t e s o l u t i o n o f a c o n v e c t i o n d i u s i o n e q u a t i o n
t h u s b e c o m e s a n i s s u e . W e b e g i n o u r d i s c u s s i o n w i t h s o m e s t a b l e a n d a c c u r a t e n u m e r i c a l
m e t h o d s f o r s o l v i n g t h e 2 D c o n v e c t i o n d i u s i o n e q u a t i o n w i t h h i g h R e y n o l d s n u m b e r s .
T h e g e n e r a l c o n v e c t i o n d i u s i o n e q u a t i o n i s o f t h e f o r m
@
2
U ( x ; y )
@ x
2
+
@
2
U ( x ; y )
@ y
2
+ p ( x ; y )
@ U ( x ; y )
@ x
+ q ( x ; y )
@ U ( x ; y )
@ y
= f ( x ; y ) ; ( x ; y ) 2 ; ( 1 )
T h i s p a p e r h a s b e e n a c c e p t e d f o r p u b l i c a t i o n i n C o m p u t e r a n d M a t h e m a t i c s w i t h A p p l i c a t i o n s
y
E - m a i l : j z h a n g @ c s . u k y . e d u . U R L : h t t p : / / w w w . c s . u k y . e d u / j z h a n g
1
-
8/14/2019 10.1.1.35.2402 (1).ps
2/16
w h e r e p ( x ; y ) a n d q ( x ; y ) a r e c o n t i n u o u s f u n c t i o n s o f v a r i a b l e s x a n d y . T h e m a g n i t u d e o f
p a n d q s i m u l a t e s t h e R e y n o l d s n u m b e r s ( d e n o t e d b y R e ) w h i c h d e t e r m i n e s t h e c o n v e c t i o n
s t r e n g t h o f a o w . i s a c o n v e x 2 D d o m a i n a n d @ i s t h e b o u n d a r y o f . T h i s e q u a t i o n
a p p e a r s i n m a n y t r a n s p o r t p r o b l e m s . I n t h e c u r r e n t s i t u a t i o n i t r e s u l t s f r o m n u m e r i c a l
s o l u t i o n o f t h e s t e a d y s t a t e i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s 6 ] .
S u p p o s e E q . ( 1 ) i s d i s c r e t i z e d b y s o m e n i t e d i e r e n c e s c h e m e a n d a l i n e a r s y s t e m
A
h
u
h
= f
h
( 2 )
i s o b t a i n e d , w h e r e h i s t h e u n i f o r m g r i d s p a c i n g o f t h e d i s c r e t i z e d s p a c e
h
. T h e l i n e a r
s y s t e m ( 2 ) i s u s u a l l y o f v e r y l a r g e d i m e n s i o n . F o r s u c h l a r g e s p a r s e l i n e a r s y s t e m s , t h e
u s e o f d i r e c t s o l u t i o n m e t h o d s b a s e d o n G a u s s i a n e l i m i n a t i o n i n c u r s p r o h i b i t i v e m e m o r y
a n d C P U c o s t a n d i t e r a t i v e s o l u t i o n m e t h o d s a r e u s u a l l y c o n s i d e r e d . T h e p e r f o r m a n c e o f
c l a s s i c a l i t e r a t i v e m e t h o d s s u c h a s J a c o b i a n d S O R i s s e n s i t i v e t o t h e m e s h s i z e h , t h e t y p e
o f b o u n d a r y c o n d i t i o n s a n d o t h e r f a c t o r s . T h e c o e c i e n t m a t r i x A
h
i s n o n s y m m e t r i c a n d
f a r f r o m d i a g o n a l l y d o m i n a n t i f R e i s l a r g e 2 2 ] . T h i s p r o p e r t y a d d s f u r t h e r d i c u l t y t o
c l a s s i c a l i t e r a t i v e m e t h o d s .
E q . ( 1 ) m a y b e d i s c r e t i z e d b y t h e c e n t r a l d i e r e n c e s c h e m e ( C D S ) , t h e r e s u l t i n g
l i n e a r s y s t e m ( 2 ) i s a 5 p o i n t f o r m u l a w i t h a t r u n c a t i o n e r r o r o f o r d e r h
2
. I n t h e c a s e o f
C D S , c l a s s i c a l i t e r a t i v e m e t h o d s f o r s o l v i n g t h e r e s u l t i n g l i n e a r s y s t e m m a y n o t c o n v e r g e
w h e n t h e c o n v e c t i v e t e r m s d o m i n a t e a n d t h e c e l l R e y n o l d s n u m b e r ( R e ) i s g r e a t e r t h a n a
c e r t a i n c o n s t a n t . F o r t h i s r e a s o n , t h e u p w i n d d i e r e n c e s c h e m e ( U P S ) h a s b e e n u s e d f o r
m a n y y e a r s d e s p i t e i t i s o n l y r s t o r d e r a c c u r a t e .
D u e t o i t s i m p o r t a n c e i n p r a c t i c a l a p p l i c a t i o n s , v a r i o u s a t t e m p t s h a v e b e e n m a d e
t o s o l v e t h e c o n v e c t i o n d i u s i o n e q u a t i o n a n d t h e i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s
w i t h i t e r a t i v e ( e s p e c i a l l y m u l t i g r i d ) m e t h o d s 3 , 1 6 , 1 7 , 1 8 ] . F o r e x a m p l e , d e Z e e u w 2 0 ]
d e v e l o p e d a b l a c k - b o x m u l t i g r i d s o l v e r w i t h s o m e m a t r i x - d e p e n d e n t p r o l o n g a t i o n s a n d
r e s t r i c t i o n s f o r s o l v i n g c o n v e c t i o n d o m i n a t e d p r o b l e m s . T h i s a n d o t h e r m e t h o d s h a v e
d e m o n s t r a t e d t h e e c i e n c y o f m u l t i g r i d t e c h n i q u e s i n s o l v i n g c o n v e c t i o n d o m i n a t e d p r o b -
l e m s d i s c r e t i z e d b y U P S .
R e c e n t l y , t h e r e h a s b e e n i n c r e a s i n g i n t e r e s t i n d e v e l o p i n g f o u r t h o r d e r c o m p a c t
s c h e m e s f o r s o l v i n g E q . ( 1 ) a n d t h e i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s w i t h l a r g e
R e y n o l d s n u m b e r s , s e e , e . g . , 2 , 4 , 6 , 9 , 1 0 , 1 4 ] . T h e s e s c h e m e s a r e s o m e w h a t s i m i l a r a n d
n u m e r i c a l r e s u l t s r e p o r t e d b y t h e s e r e s e a r c h e r s a r e n o t v e r y d i e r e n t . T h e r e a r e a t l e a s t
t h r e e a d v a n t a g e s s h a r e d b y t h e s e s c h e m e s :
1 . U n c o n d i t i o n a l s t a b i l i t y : a l t h o u g h t h e c o e c i e n t m a t r i c e s a r e n o l o n g e r d i a g o n a l l y
d o m i n a n t f o r l a r g e R e , t h e s c h e m e s h a v e b e e n s h o w n n u m e r i c a l l y s t a b l e f o r l a r g e
R e y n o l d s n u m b e r s 4 , 8 , 2 3 ] ;
2 . H i g h a c c u r a c y : i t h a s b e e n s h o w n t h a t t h e s e s c h e m e s d o p r o d u c e n u m e r i c a l s o l u t i o n
o f f o u r t h o r d e r a c c u r a c y f o r t h e c o n v e c t i o n d i u s i o n e q u a t i o n 4 , 8 , 2 1 , 2 3 ] a n d o f
h i g h a c c u r a c y f o r t h e N a v i e r - S t o k e s e q u a t i o n s w i t h s m a l l t o m e d i u m R e 1 , 6 , 1 0 ] ;
2
-
8/14/2019 10.1.1.35.2402 (1).ps
3/16
3 . E a s y b o u n d a r y t r e a t m e n t : s i n c e t h e c o m p u t a t i o n a l s t e n c i l i n v o l v e s o n l y t h e n e a r e s t
9 g r i d p o i n t s , t h e s c h e m e s a r e o f c o m p a c t t y p e a n d n o s p e c i a l f o r m u l a i s n e e d e d f o r
a p p r o x i m a t i n g g r i d p o i n t s n e a r t h e b o u n d a r i e s o f a r e g u l a r l y s t r u c t u r e d d o m a i n .
H o w e v e r , u n t i l r e c e n t l y , c o m p u t a t i o n a l a d v a n t a g e s o f t h e s e c o m p a c t s c h e m e s h a v e
n o t b e e n f u l l y i n v e s t i g a t e d . F o r e x a m p l e , i t i s n o t c l e a r i f t h e s e s c h e m e s c a n b e u s e d t o
s o l v e i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s o f v e r y l a r g e R e b e c a u s e o f t h e l i m i t a t i o n s o f
a v a i l a b l e c o m p u t e r p o w e r a n d o f t h e d i c u l t y o f s o l v i n g t h e r e s u l t i n g l i n e a r s y s t e m ( l a c k
o f d i a g o n a l d o m i n a n c e ) w i t h t r a d i t i o n a l S O R t y p e i t e r a t i v e m e t h o d s .
T o f u l l y i n v e s t i g a t e t h e p o t e n t i a l o f u s i n g t h e f o u r t h o r d e r c o m p a c t s c h e m e s f o r
s o l v i n g N a v i e r - S t o k e s e q u a t i o n s w i t h l a r g e R e , n o n t r a d i t i o n a l i t e r a t i v e m e t h o d s s e e m n e c -
e s s a r y . O n e c l a s s o f p r o m i s i n g m e t h o d s a r e m u l t i g r i d t e c h n i q u e s w h i c h h a v e b e e n s u c c e s s -
f u l l y u s e d w i t h t h e r s t a n d s e c o n d o r d e r d i s c r e t i z a t i o n s c h e m e s f o r s o l v i n g p r o b l e m s i n
c o m p u t a t i o n a l u i d d y n a m i c s ( i n c l u d i n g t h e d r i v e n c a v i t y p r o b l e m ) ( s e e , e . g . , 3 , 1 7 , 1 9 ] ) .
A p r e l i m i n a r y i n v e s t i g a t i o n o n c o m b i n i n g t h e f o u r t h o r d e r c o m p a c t s c h e m e s w i t h m u l t i -
g r i d t e c h n i q u e s w a s m a d e b y A l t a s a n d B u r r a g e 1 ] f o r d i u s i o n d o m i n a t e d o w p r o b l e m s .
M u l t i g r i d s o l u t i o n a n d a c c e l e r a t e d m u l t i g r i d s o l u t i o n m e t h o d s w i t h t h e f o u r t h o r d e r c o m -
p a c t s c h e m e s f o r s o l v i n g c o n v e c t i o n d o m i n a t e d p r o b l e m s a r e r e l a t i v e l y n e w 8 , 2 3 , 2 4 ]
T h e c u r r e n t p a p e r i s t o u s e m u l t i g r i d t e c h n i q u e s t o e v a l u a t e t h e f o u r t h o r d e r c o m -
p a c t s c h e m e s i n s o l v i n g s t e a d y s t a t e i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s f o r l a r g e R e ,
e x e m p l i e d b y s o l v i n g t h e s q u a r e d r i v e n c a v i t y p r o b l e m . T h i s f o l l o w s f r o m G u p t a ' s w o r k
6 ] u s i n g S O R i t e r a t i v e m e t h o d w i t h t h e s e s c h e m e s t o s o l v e d r i v e n c a v i t y p r o b l e m f o r
R e 2 ; 0 0 0 a n d Z h a n g ' s w o r k 2 1 ] e m p l o y i n g t h e f o u r t h o r d e r c o m p a c t s c h e m e a n d m u l t i -
g r i d t e c h n i q u e s t o c o m p u t e h i g h a c c u r a c y s o l u t i o n o f t h e c o n v e c t i o n d i u s i o n e q u a t i o n
w i t h v e r y l a r g e R e y n o l d s n u m b e r s .
W e p r e s e n t t h e 9 p o i n t c o m p a c t n i t e d i e r e n c e d i s c r e t i z a t i o n s c h e m e f o r E q . ( 1 )
a n d d i s c r e t i z a t i o n s c h e m e s f o r t h e s t r e a m f u n c t i o n a n d v o r t i c i t y f o r m u l a t i o n o f t h e i n c o m -
p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s i n S e c t i o n 2 . T h e m u l t i g r i d s o l v e r s f o r t h e c o n v e c t i o n
d i u s i o n ( a n d t h e P o i s s o n ) e q u a t i o n a n d f o r t h e i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s
a r e o u t l i n e d i n S e c t i o n 3 . I n S e c t i o n 4 , w e s o l v e t h e d r i v e n c a v i t y p r o b l e m a n d c o m -
p a r e o u r n u m e r i c a l r e s u l t s w i t h t h o s e o b t a i n e d b y o t h e r r e s e a r c h e r s u s i n g o t h e r m e t h o d s .
C o n c l u d i n g r e m a r k s a r e g i v e n i n S e c t i o n 5 .
2 F o u r t h O r d e r F i n i t e D i e r e n c e S c h e m e s
T h e a p p r o x i m a t e v a l u e o f a f u n c t i o n U ( x ; y ) a t a m e s h p o i n t ( x ; y ) i s d e n o t e d b y U
0
.
T h e a p p r o x i m a t e v a l u e s a t i t s e i g h t i m m e d i a t e n e i g h b o r i n g p o i n t s a r e d e n o t e d b y U
i
; i =
1 ; 2 ; : : : ; 8 . T h e c o m p u t a t i o n a l s t e n c i l o f t h e 9 c o m p a c t g r i d p o i n t s i s s h o w n a s
0
B
@
u
6
u
2
u
5
u
3
u
0
u
1
u
7
u
4
u
8
1
C
A
:
3
-
8/14/2019 10.1.1.35.2402 (1).ps
4/16
T h e d i s c r e t i z e d v a l u e s o f p
i
; q
i
a n d f
i
; i = 0 ; 1 ; : : : ; 4 ; h a v e t h e i r o b v i o u s m e a n i n g s . T h e
c o m p a c t n i t e d i e r e n c e f o r m u l a f o r a r e f e r e n c e p o i n t ( x ; y ) i n v o l v e s t h e n e a r e s t e i g h t
n e i g h b o r i n g m e s h p o i n t s w i t h t h e m e s h s p a c i n g h a n d i s g i v e n b y ( f o r d e t a i l e d p r o c e d u r e
t o d e r i v e t h e n i t e d i e r e n c e s c h e m e , s e e 4 ] ) :
8
X
i = 0
i
U
i
=
h
2
2
8 f
0
+ f
1
+ f
2
+ f
3
+ f
4
] +
h
3
4
p
0
( f
1
f
3
) + q
0
( f
2
f
4
) ] ; ( 3 )
w h e r e t h e c o e c i e n t s
i
; i = 0 ; 1 : : : ; 8 ; a r e g i v e n a s
1
= 4 +
h
4
4 p
0
+ 3 p
1
p
3
+ p
2
+ p
4
] +
h
2
8
4 p
2
0
+ p
0
( p
1
p
3
) + q
0
( p
2
p
4
) ] ;
2
= 4 +
h
4
4 q
0
+ 3 q
2
q
4
+ q
1
+ q
3
] +
h
2
8
4 q
2
0
+ p
0
( q
1
q
3
) + q
0
( q
2
q
4
) ] ;
3
= 4
h
4
4 p
0
p
1
+ 3 p
3
+ p
2
+ p
4
] +
h
2
8
4 p
2
0
p
0
( p
1
p
3
) q
0
( p
2
p
4
) ] ;
4
= 4
h
4
4 q
0
q
2
+ 3 q
4
+ q
1
+ q
3
] +
h
2
8
4 q
2
0
p
0
( q
1
q
3
) q
0
( q
2
q
4
) ] ;
5
= 1 +
h
2
( p
0
+ q
0
) +
h
8
( q
1
q
3
+ p
2
p
4
) +
h
2
4
p
0
q
0
;
6
= 1
h
2
( p
0
q
0
)
h
8
( q
1
q
3
+ p
2
p
4
)
h
2
4
p
0
q
0
;
7
= 1
h
2
( p
0
+ q
0
) +
h
8
( q
1
q
3
+ p
2
p
4
) +
h
2
4
p
0
q
0
;
8
= 1 +
h
2
( p
0
q
0
)
h
8
( q
1
q
3
+ p
2
p
4
)
h
2
4
p
0
q
0
;
0
= 2 0 + h
2
( p
2
0
+ q
2
0
) + h ( p
1
p
3
) + h ( q
2
q
4
) ] :
R e s u l t s o f n u m e r i c a l e x p e r i m e n t s i n 4 ] s h o w t h a t S O R i t e r a t i v e m e t h o d w i t h t h i s s c h e m e
c o n v e r g e s f o r a n y v a l u e s o f p ( x ; y ) a n d q ( x ; y ) t e s t e d .
W h e n R e 0 , E q . ( 1 ) r e d u c e s t o t h e P o i s s o n e q u a t i o n , a n d E q . ( 3 ) r e d u c e s t o t h e
w e l l k n o w n ( s i m p l e r ) M e h r s t e l l e n f o r m u l a 5 ] :
4 U
1
+ U
2
+ U
3
+ U
4
] + U
5
+ U
6
+ U
7
+ U
8
2 0 U
0
=
1
2
h
2
8 f
0
+ f
1
+ f
2
+ f
3
+ f
4
] : ( 4 )
M u l t i g r i d a p p l i c a t i o n s o f t h e M e h r s t e l l e n f o r m u l a h a v e b e e n i n v e s t i g a t e d b y S c h a e r 1 3 ] ,
G u p t a , K o u a t c h o u a n d Z h a n g 7 ] .
T h e N a v i e r - S t o k e s e q u a t i o n s r e p r e s e n t i n g t h e t w o d i m e n s i o n a l s t e a d y o w o f a n i n -
c o m p r e s s i b l e v i s c o u s u i d a r e g i v e n i n s t r e a m f u n c t i o n a n d v o r t i c i t y f o r m u l a t i o n a s f o l l o w s
3 , 6 , 1 8 ] :
@
2
( x ; y )
@ x
2
+
@
2
( x ; y )
@ y
2
= ( x ; y ) ; ( 5 )
@
2
( x ; y )
@ x
2
+
@
2
( x ; y )
@ y
2
R e
h
u ( x ; y )
@ ( x ; y )
@ x
+ v ( x ; y )
@ ( x ; y )
@ y
i
= 0 ; ( 6 )
u ( x ; y ) =
@ ( x ; y )
@ y
; v ( x ; y ) =
@ ( x ; y )
@ x
: ( 7 )
4
-
8/14/2019 10.1.1.35.2402 (1).ps
5/16
H e r e i s t h e s t r e a m f u n c t i o n , t h e v o r t i c i t y ; u a n d v a r e t h e v e l o c i t i e s ; R e i s t h e n o n d i -
m e n s i o n a l R e y n o l d s n u m b e r 6 ] .
T h e s t r e a m f u n c t i o n ( 5 ) i s a P o i s s o n e q u a t i o n a n d t h e f o u r t h o r d e r a p p r o x i m a t i o n i s
g i v e n b y E q . ( 4 ) w i t h U = a n d f = . T h e v o r t i c i t y e q u a t i o n ( 6 ) i s a s p e c i a l c a s e
o f t h e c o n v e c t i o n d i u s i o n e q u a t i o n ( 1 ) a n d t h e f o u r t h o r d e r a p p r o x i m a t i o n i n t h i s c a s e
m a y b e o b t a i n e d w i t h U = , f = 0 a n d p ( x ; y ) = R e u ( x ; y ) , q ( x ; y ) = R e v ( x ; y ) i n
E q . ( 3 ) .
T h e v e l o c i t i e s u , v a t a g r i d p o i n t ( x ; y ) a r e c a l c u l a t e d f r o m t h e d i s c r e t e a p p r o x i -
m a t i o n o f E q . ( 7 ) . I t h a s b e e n s h o w n 6 ] t h a t i t i s b e n e c i a l f o r b o t h c o n v e r g e n c e a n d
a c c u r a c y t o u s e t h e f o u r t h o r d e r a p p r o x i m a t i o n s f o r t h e v e l o c i t i e s . I n p a r t i c u l a r , G u p t a
5 ] d e r i v e d s o m e h i g h a c c u r a c y c o m p a c t a p p r o x i m a t i o n s f o r t h e g r a d i e n t s o f t h e s o l u t i o n
o f t h e P o i s s o n e q u a t i o n . S i n c e t h e s t r e a m f u n c t i o n e q u a t i o n ( 5 ) i s a P o i s s o n e q u a t i o n i n
, h i g h a c c u r a c y a p p r o x i m a t i o n s f o r t h e g r a d i e n t
x
a n d
y
c a n b e o b t a i n e d f r o m 5 ] ,
a n d t h e c o r r e s p o n d i n g f o u r t h o r d e r c o m p a c t a p p r o x i m a t i o n s f o r t h e v e l o c i t i e s a r e g i v e n
a s ( a l s o s e e 6 ] ) :
u
0
= (
2
4
) = 3 h + (
5
6
7
8
) = 1 2 h + h (
2
4
) = 1 2 ;
v
0
= (
3
1
) = 3 h (
5
6
7
+
8
) = 1 2 h + h (
3
1
) = 1 2 :
)
( 8 )
3 S o l u t i o n M e t h o d
I n t h i s p a p e r , w e a r e m a i n l y c o n c e r n e d w i t h t h e a c c u r a c y o f t h e c o m p u t e d s o l u t i o n , n o t
t h e c o m p u t a t i o n a l e c i e n c y o f t h e s o l u t i o n m e t h o d t h a t w e u s e d . W e j u s t o u t l i n e o u r
s o l u t i o n m e t h o d i n t h i s s e c t i o n a n d p o i n t o u t r e l e v a n t r e f e r e n c e s i n t h i s d i r e c t i o n 8 , 1 8 ,
1 9 , 2 0 , 2 1 , 2 3 , 2 7 ] .
O u r m u l t i g r i d s o l u t i o n m e t h o d f o r t h e c o n v e c t i o n d i u s i o n e q u a t i o n ( 1 ) ( t h e P o i s s o n
e q u a t i o n ( 4 ) i s a s p e c i a l c a s e ) i s a s f o l l o w s :
1 . S t a r t o n t h e n e g r i d w i t h s o m e i n i t i a l g u e s s a n d p e r f o r m
1
r e d - b l a c k G a u s s - S e i d e l
r e l a x a t i o n s w e e p s .
2 . C o m p u t e r e s i d u a l s o n t h e n e g r i d p o i n t s w h i c h a r e c o m m o n t o t h e c o a r s e g r i d a n d
m u l t i p l y t h e r e s i d u a l s b y a s c a l i n g f a c t o r a n d i n j e c t t h e s c a l e d r e s i d u a l s t o t h e
c o a r s e g r i d 2 3 , 2 5 , 2 6 ] .
3 . P e r f o r m m u l t i g r i d c y c l e s o n t h e c o a r s e g r i d .
4 . I n t e r p o l a t e t h e c o a r s e g r i d c o r r e c t i o n t o t h e n e g r i d b y b i l i n e a r i n t e r p o l a t i o n a n d
a d d t h e c o r r e c t i o n t o t h e n e g r i d s o l u t i o n .
5 . P e r f o r m
2
r e d - b l a c k G a u s s - S e i d e l r e l a x a t i o n s w e e p s o n t h e n e g r i d .
I f = 1 , t h e m u l t i g r i d a l g o r i t h m i s c a l l e d a V c y c l e . I t i s a W c y c l e f o r = 2 .
1
a n d
2
a r e t h e n u m b e r s o f p r e s m o o t h i n g a n d p o s t s m o o t h i n g s w e e p s .
5
-
8/14/2019 10.1.1.35.2402 (1).ps
6/16
F o r t h e s t r e a m f u n c t i o n e q u a t i o n ( 5 ) , a V ( 1 , 1 ) c y c l e a l g o r i t h m i s s u c i e n t t o o b t a i n
a c c u r a t e s o l u t i o n w i t h a c c e p t a b l e c o n v e r g e n c e a n d i s m o r e c o s t e e c t i v e t h a n a W ( 1 , 1 )
c y c l e a l g o r i t h m . T h e s m o o t h e r f o r t h e s t r e a m f u n c t i o n i s t h e r e d - b l a c k G a u s s - S e i d e l w i t h -
o u t a p a r a m e t e r . F o r t h e v o r t i c i t y e q u a t i o n ( 6 ) , w e u s e a W ( 1 , 1 ) c y c l e a l g o r i t h m . T h e
s m o o t h e r i s t h e r e d - b l a c k G a u s s - S e i d e l w i t h a d a m p i n g p a r a m e t e r !
i
2 ( 0 ; 1 ) ( a l l d a m p i n g
f a c t o r s w e r e f o u n d b y n u m e r i c a l e x p e r i m e n t s ) .
T h e N a v i e r - S t o k e s e q u a t i o n s ( 5 ) { ( 7 ) m a y b e s o l v e d b y a n e s t e d i n n e r - o u t e r i t e r a -
t i o n p r o c e d u r e ( s e e 6 ] ) w i t h d i e r e n t m u l t i g r i d c y c l i n g a l g o r i t h m s b e i n g a p p l i e d t o t h e
s t r e a m f u n c t i o n e q u a t i o n ( 5 ) a n d t h e v o r t i c i t y e q u a t i o n ( 6 ) . W h i l e i t i s s o m e t i m e s a d v a n -
t a g e o u s t o s o l v e E q s . ( 5 ) { ( 7 ) s i m u l t a n e o u s l y t o m a i n t a i n t h e p h y s i c a l c o u p l i n g b e t w e e n
t h e m 3 , 1 8 ] , t h e r e l a x a t i o n m e t h o d s ( s u c h a s t h e c o u p l e d s t r o n g l y i m p l i c i t l y ( C S I ) p r o -
c e d u r e 3 ] a n d i n c o m p l e t e L U d e c o m p o s i t i o n 1 8 ] ) u s e d t o a c c o m p l i s h t h i s c o u p l i n g a r e
e x p e n s i v e a n d d i c u l t t o i m p l e m e n t , c o m p a r i n g w i t h t h e p o i n t G a u s s - S e i d e l r e l a x a t i o n .
S i n c e o u r g o a l i s t o e v a l u a t e t h e u s e f u l n e s s o f t h e f o u r t h o r d e r c o m p a c t s c h e m e s f o r s o l v i n g
N a v i e r - S t o k e s e q u a t i o n s w i t h l a r g e R e , w e u s e d t h e s i m p l e s t p o i n t G a u s s - S e i d e l r e l a x a t i o n .
4 S o l u t i o n o f S q u a r e D r i v e n C a v i t y
T h e s t e a d y o w o f a n i n c o m p r e s s i b l e v i s c o u s u i d i n a s q u a r e c a v i t y = ( 0 ; 1 ) ( 0 ; 1 ) h a s
b e e n u s e d f o r m a n y y e a r s a s t h e m o d e l p r o b l e m b y r e s e a r c h e r s t o t e s t t h e i r n e w n u m e r i c a l
s c h e m e s a n d s o l u t i o n m e t h o d s 1 , 3 , 6 , 1 7 , 1 8 ] , a l t h o u g h t h e r e i s n o a n a l y t i c a l s o l u t i o n
a v a i l a b l e a n d t h e r e a r e s i n g u l a r i t i e s a t t w o o f i t s c o r n e r s . H i g h l y a c c u r a t e b e n c h m a r k
s o l u t i o n s a r e a v a i l a b l e i n t h e l i t e r a t u r e . I n p a r t i c u l a r , G h i a e t a l . 3 ] u s e d m u l t i g r i d
t e c h n i q u e s a n d g r i d p o i n t s o f 2 5 6 2 5 6 t o c o m p u t e n u m e r i c a l s o l u t i o n s f o r 1 0 0 R e
1 0 ; 0 0 0 . T h e i r s o l u t i o n s h a v e b e e n c o n s i d e r e d v e r y a c c u r a t e n o t o n l y b e c a u s e o f t h e s m a l l
g r i d s p a c i n g e m p l o y e d , b u t a l s o a s a r e s u l t o f t h e f a c t t h a t a n d w e r e c o m p u t e d d i r e c t l y
f r o m t h e s t r e a m f u n c t i o n a n d v o r t i c i t y f o r m u l a t i o n ( 5 ) { ( 7 ) .
T h e o w i s i n d u c e d b y t h e s l i d i n g m o t i o n o f t h e t o p w a l l ( y = 1 ) f r o m r i g h t t o l e f t
a n d i s d e s c r i b e d b y t h e N a v i e r - S t o k e s e q u a t i o n s ( 5 ) { ( 7 ) . T h e b o u n d a r y c o n d i t i o n s a r e
t h o s e o f n o s l i p : o n t h e s t a t i o n a r y w a l l s u = 0 a n d v = 0 ; o n t h e s l i d i n g w a l l u = 1 a n d
v = 0 ( s e e F i g . 1 ) .
I n o r d e r t o s o l v e t h e d r i v e n c a v i t y p r o b l e m , w e r e p l a c e t h e N a v i e r - S t o k e s e q u a t i o n s
( 5 ) { ( 6 ) b y t h e n i t e d i e r e n c e a p p r o x i m a t i o n s g i v e n i n E q s . ( 3 ) a n d ( 4 ) r e s p e c t i v e l y . T h e
v e l o c i t i e s , d e n e d i n E q . ( 7 ) , a r e c a l c u l a t e d b y u s i n g t h e f o u r t h o r d e r a p p r o x i m a t i o n s
( 8 ) . T h e u n i t s q u a r e i s c o v e r e d b y a g r i d o f u n i f o r m m e s h s i z e h = 1 = N . T h e d i s c r e t e
a p p r o x i m a t i o n s ( 3 ) a n d ( 4 ) a r e w r i t t e n a t e a c h o f t h e ( N 1 )
2
i n t e r i o r g r i d p o i n t s . Z e r o
v a l u e s a r e p r e s c r i b e d f o r o n t h e b o u n d a r y . T h e u s u a l a p p r o x i m a t i o n s f o r v o r t i c i t y o n
t h e b o u n d a r y a r e t h e J e n s e n ' s f o r m u l a s ( s e e 6 , 1 2 ] ) , w h i c h h a v e a l o c a l t r u n c a t i o n e r r o r
o f s e c o n d o r d e r . F o u r t h o r d e r a p p r o x i m a t i o n s c o u l d a l s o b e d e n e d t o o b t a i n b o u n d a r y
v a l u e s o f . I n p a r t i c u l a r , s o m e f o u r t h o r d e r a p p r o x i m a t i o n s a n a l o g o u s t o t h e J e n s e n ' s
6
-
8/14/2019 10.1.1.35.2402 (1).ps
7/16
-
8/14/2019 10.1.1.35.2402 (1).ps
8/16
-
8/14/2019 10.1.1.35.2402 (1).ps
9/16
h L o c a t i o n D i . i n D i . i n
1 / 1 6 0 . 1 0 2 9 0 1 3 . 3 0 7 2 7 2 ( 0 . 3 7 5 0 0 0 , 0 . 7 5 0 0 0 0 ) 0 . 5 9 % 4 . 3 7 %
1 / 3 2 0 . 1 0 3 3 5 4 3 . 2 8 2 3 9 5 ( 0 . 3 7 5 0 0 0 , 0 . 7 5 0 0 0 0 ) 0 . 1 5 % 3 . 5 9 %
1 / 6 4 0 . 1 0 3 5 0 1 3 . 1 0 8 6 4 3 ( 0 . 3 9 0 6 2 5 , 0 . 7 3 4 3 7 5 ) 0 . 0 1 % 1 . 9 0 %
1 / 1 2 8 0 . 1 0 3 5 1 1 3 . 1 6 8 7 4 5 ( 0 . 3 8 2 8 1 3 , 0 . 7 3 4 3 7 5 ) 0 . 0 0 % 0 . 0 0 %
T a b l e 1 : V a l u e s a n d l o c a t i o n s o f t h e p r i m a r y v o r t e x f o r R e = 1 0 0 u s i n g d i e r e n t d i s -
c r e t i z a t i o n s .
4 . 1 C o m p a r i s o n w i t h b e n c h m a r k s o l u t i o n s
U n l e s s o t h e r w i s e i n d i c a t e d e x p l i c i t l y , w e c o m p a r e o u r r e s u l t s w i t h t h o s e o b t a i n e d b y G h i a
e t a l . 3 ] a s t h e b e n c h m a r k s o l u t i o n s . O u r p r o b l e m w a s s e t u p s l i g h t l y d i e r e n t f r o m t h a t
o f G h i a e t a l . 3 ] a n d t h e u v e l o c i t y a t t h e t o p w a l l ( y = 1 ) i s d i e r e n t d u e t o t h e f a c t
t h a t t h e o w i n o u r p r o b l e m i s i n d u c e d b y t h e s l i d i n g m o t i o n o f t h e t o p w a l l f r o m r i g h t
t o l e f t 6 ] . I t w a s i n d u c e d f r o m l e f t t o r i g h t i n t h e m o d e l p r o b l e m s o l v e d b y G h i a e t a l .
i n 3 ] . T h e c o m p u t e d v a l u e a t a g r i d p o i n t ( x ; y ) l i s t e d i n o u r t a b l e s s h o u l d b e c o m p a r e d
w i t h t h a t a t t h e p o i n t ( 1 x ; y ) i n t a b l e s o f 3 ] .
W e r s t s o l v e d t h e m o d e l p r o b l e m f o r R e = 1 0 0 w i t h h = 1 = 1 6 ; 1 = 3 2 ; 1 = 6 4 ; 1 = 1 2 8 a n d
c o m p a r e o u r r e s u l t s w i t h t h o s e o f G h i a e t a l . w i t h h = 1 = 1 2 8 . T h e n u m e r i c a l v a l u e s a r e
g i v e n i n T a b l e 1 . I t c a n b e s e e n t h a t , i f a 5 % d e p a r t u r e f r o m t h e b e n c h m a r k s o l u t i o n i s
a c c e p t a b l e a s t h e e n g i n e e r i n g a c c u r a c y , o u r m e t h o d w i t h h = 1 = 1 6 c a n p r o d u c e a c c e p t a b l e
r e s u l t s , w h i l e G h i a e t a l . u s e d m u c h n e r d i s c r e t i z a t i o n ( h = 1 = 1 2 8 ) . N o t e t h a t w h e n t h e
m e s h i s r e n e d , t h e a c c u r a c y o f o u r s o l u t i o n i s i n c r e a s e d r a p i d l y .
W h e n R e i n c r e a s e s t o 1 ; 0 0 0 , r e s u l t s i n T a b l e 2 s h o w t h a t n e r m e s h i s n e e d e d t o
p r o d u c e a c c u r a t e s o l u t i o n . H o w e v e r , w e c a n s e e t h a t w i t h o n l y a q u a r t e r o f t h e n u m b e r
o f e q u a t i o n s u s e d b y G h i a e t a l . , o u r m e t h o d p r o d u c e d h i g h a c c u r a c y s o l u t i o n s .
W h e n R e 2 ; 0 0 0 , G u p t a 6 ] r e p o r t e d s l o w c o n v e r g e n c e w h e n t h e f o u r t h o r d e r c o m -
p a c t s c h e m e s w e r e u s e d w i t h S O R i t e r a t i o n . I t i s s h o w n i n T a b l e 3 t h a t o u r m u l t i g r i d
a l g o r i t h m c a n c o m p u t e v e r y a c c u r a t e s o l u t i o n f o r R e = 3 ; 2 0 0 w i t h h = 1 = 1 2 8 .
W i t h R e = 5 ; 0 0 0 , G h i a e t a l . h a d t o u s e h = 1 = 2 5 6 t o c o m p u t e a c c u r a t e s o l u t i o n ,
w e f o u n d ( s e e T a b l e 4 ) t h a t o u r f o u r t h o r d e r m u l t i g r i d m e t h o d c a n c o m p u t e s o l u t i o n t o
c o m p a r a b l e a c c u r a c y u s i n g h = 1 = 1 2 8 .
F o r R e = 7 ; 5 0 0 , w i t h h = 1 = 1 2 8 , o u r m e t h o d w a s s t i l l a b l e t o y i e l d s o l u t i o n a c c u r a t e
e n o u g h ( 5 % e r r o r ) t o c o m p a r e w i t h t h e b e n c h m a r k s o l u t i o n o f G h i a e t a l . 3 ] u s i n g
h = 1 = 2 5 6 .
4 . 2 C o m p a r i s o n o f h i g h a c c u r a c y s o l u t i o n s
I n S e c t i o n 4 . 1 , t h e d i e r e n c e s f o r t h e v o r t i c i t y v a l u e s a r e u s u a l l y l a r g e r t h a n t h a t o f t h e
s t r e a m f u n c t i o n v a l u e s f o r l a r g e R e , a s i n d i c a t e d b y d a t a i n T a b l e 3 f o r R e = 3 ; 2 0 0 a n d i n
9
-
8/14/2019 10.1.1.35.2402 (1).ps
10/16
h L o c a t i o n D i . i n D i . i n
1 / 1 6 0 . 0 9 1 2 0 8 3 . 0 7 3 3 3 2 ( 0 . 4 3 7 5 0 0 , 0 . 7 5 0 0 0 0 ) 2 3 . 3 6 % 4 8 . 5 7 %
1 / 3 2 0 . 1 0 6 7 2 3 2 . 0 2 6 7 2 4 ( 0 . 4 6 8 7 5 0 , 0 . 5 9 3 7 5 0 ) 1 3 . 3 2 % 2 . 0 2 %
1 / 6 4 0 . 1 1 7 3 8 6 2 . 0 5 6 6 0 8 ( 0 . 4 6 8 7 5 0 , 0 . 5 6 2 5 0 0 ) 1 . 3 6 % 0 . 5 8 %
1 / 1 2 8 0 . 1 1 8 8 0 6 2 . 0 6 6 7 7 7 ( 0 . 4 6 8 7 5 0 , 0 . 5 6 2 5 0 0 ) 0 . 1 7 % 0 . 0 9 %
T a b l e 2 : V a l u e s a n d l o c a t i o n s o f t h e p r i m a r y v o r t e x f o r R e = 1 ; 0 0 0 u s i n g d i e r e n t d i s -
c r e t i z a t i o n s .
h L o c a t i o n D i . i n D i . i n
1 / 6 4 0 . 1 0 9 0 9 4 1 . 8 9 3 5 3 4 ( 0 . 4 8 4 3 7 5 , 0 . 5 7 8 1 2 5 ) 9 . 9 5 % 2 . 9 3 %
1 / 1 2 8 0 . 1 2 0 1 5 7 1 . 9 4 8 9 3 4 ( 0 . 4 8 4 3 7 5 , 0 . 5 3 9 0 6 3 ) 0 . 1 8 % 4 . 0 1 %
T a b l e 3 : V a l u e s a n d l o c a t i o n s o f t h e p r i m a r y v o r t e x f o r R e = 3 ; 2 0 0 u s i n g d i e r e n t d i s -
c r e t i z a t i o n s .
h L o c a t i o n D i . i n D i . i n
1 / 1 2 8 0 . 1 1 8 1 2 1 1 . 9 0 6 2 1 4 ( 0 . 4 8 4 3 7 5 , 0 . 5 3 9 0 6 3 ) 0 . 7 1 % 2 . 4 8 %
T a b l e 4 : V a l u e s a n d l o c a t i o n s o f t h e p r i m a r y v o r t e x f o r R e = 5 ; 0 0 0 u s i n g h = 1 2 8 .
h L o c a t i o n D i . i n D i . i n
1 / 1 2 8 0 . 1 1 4 3 3 8 1 . 8 7 5 8 9 7 ( 0 . 4 9 2 1 8 8 , 0 . 5 5 4 6 8 8 ) 4 . 2 8 % 0 . 2 1 %
T a b l e 5 : V a l u e s a n d l o c a t i o n s o f t h e p r i m a r y v o r t e x f o r R e = 7 ; 5 0 0 u s i n g h = 1 = 1 2 8 .
1 0
-
8/14/2019 10.1.1.35.2402 (1).ps
11/16
A u t h o r s M e s h s i z e M e t h o d O r d e r
C u r r e n t a u t h o r 0 : 1 2 0 1 5 7 1 : 9 4 8 9 3 4 1 2 8 1 2 8 4 t h - o r d e r
C u r r e n t a u t h o r 0 : 1 2 1 6 5 9 1 : 9 5 9 4 8 7 2 5 6 2 5 6 4 t h - o r d e r
L i , T a n g & F o r n b e r g 1 0 ] 0 : 1 2 0 5 2 9 1 : 9 4 2 8 6 1 2 8 1 2 8 4 t h - o r d e r
N i s h i d a & S a t o f u k a 1 1 ] 0 : 1 2 1 1 5 4 1 : 9 5 0 7 8 1 2 8 1 2 8 6 t h - o r d e r
G h i a , G h i a & S h i n 3 ] 0 : 1 2 0 3 7 7 1 : 9 8 8 6 0 1 2 8 1 2 8 2 n d - o r d e r
T a b l e 6 : C o m p a r i s o n o f h i g h a c c u r a c y s o l u t i o n f o r R e = 3 ; 2 0 0 .
T a b l e 4 f o r R e = 5 ; 0 0 0 . T h i s m a y b e c a u s e d b y t h e i n a c c u r a c y o f t h e b e n c h m a r k s o l u t i o n
f o r l a r g e R e . A s i n d i c a t e d i n 1 1 ] a h i g h e r o r d e r o f a c c u r a c y i n s p a c e i s n e c e s s a r y f o r ( t i m e
d e p e n d e n t ) h i g h R e y n o l d s n u m b e r s i m u l a t i o n s .
I n T a b l e 6 , w e g i v e s o m e r e c e n t l y a v a i l a b l e h i g h a c c u r a c y r e s u l t s f r o m s e v e r a l h i g h e r
o r d e r m e t h o d s f o r R e = 3 ; 2 0 0 a n d c o m p a r e t h e r e s u l t s w i t h t h o s e o b t a i n e d b y o u r s c h e m e
w i t h h = 1 = 1 2 8 a n d h = 1 = 2 5 6 .
I t c a n b e s e e n f r o m T a b l e 6 t h a t G h i a e t a l . ' s s o l u t i o n i s r e l a t i v e l y i n a c c u r a t e c o m -
p a r e d w i t h t h e h i g h o r d e r , h i g h a c c u r a c y s o l u t i o n r e c e n t l y a v a i l a b l e i n t h e l i t e r a t u r e . I n
p a r t i c u l a r , t h e v o r t i c i t y v a l u e i s s h o w n t o h a v e a r e l a t i v e e r r o r o f 1 : 4 9 % w i t h r e s p e c t t o
o u r v e r y a c c u r a t e s o l u t i o n c o m p u t e d b y u s i n g o u r f o u r t h o r d e r a l g o r i t h m w i t h h = 1 = 2 5 6 .
I n c o n t r a s t , o u r s o l u t i o n w i t h h = 1 = 1 2 8 h a s o n l y a r e l a t i v e e r r o r o f 0 : 5 4 % i n v o r t i c i t y
v a l u e w i t h r e s p e c t t o o u r v e r y a c c u r a t e s o l u t i o n , b u t h a s a r e l a t i v e d i e r e n c e o f 4 : 0 1 % w i t h
r e s p e c t t o G h i a e t a l . ' s s o l u t i o n . I t c a n b e o b s e r v e d t h a t t h e s o l u t i o n f r o m o u r f o u r t h o r d e r
c o m p a c t a l g o r i t h m i s c o m p a r a b l e w i t h t h a t o f L i e t a l . ' s f o u r t h o r d e r c o m p a c t s c h e m e 1 0 ] ,
t h e r e l a t i v e e r r o r s i n b o t h s t r e a m f u n c t i o n a n d v o r t i c i t y v a l u e s a r e l e s s t h a n 0 : 3 2 % . F u r -
t h e r m o r e , t h e s o l u t i o n f r o m o u r f o u r t h o r d e r m e t h o d u s i n g h = 1 = 2 5 6 i s c o m p a r a b l e w i t h
N i s h i d a a n d S a t o f u k a ' s s i x t h - o r d e r s p e c t r a l m e t h o d 1 1 ] u s i n g h = 1 = 1 2 8 w i t h a d i e r e n c e
o f 0 : 4 5 % . T h i s c o m p a r i s o n a g a i n s u p p o r t s t h e c l a i m m a d e b y N i s h i d a a n d S a t o f u k a 1 1 ]
t h a t h i g h e r o r d e r a l g o r i t h m i s n e c e s s a r y ( a t l e a s t b e n e c i a l ) f o r h i g h R e y n o l d s n u m b e r
c o m p u t a t i o n s .
4 . 3 S o l u t i o n c o n t o u r s
T h e s t r e a m l i n e s a n d t h e v o r t i c i t y c o n t o u r s w i t h R e = 5 ; 0 0 0 a n d 7 ; 5 0 0 f o r h = 1 = 1 2 8 a r e
p r e s e n t e d i n F i g s . 2 a n d 3 . T h e s e g u r e s a r e g r a p h i c a l l y c o m p a r a b l e w i t h t h e w e l l k n o w n
g u r e s o b t a i n e d b y G h i a e t a l . 3 ] , t a k i n g i n t o a c c o u n t t h e d i e r e n c e i n t h e d e n i t i o n o f
t h e v e l o c i t y ( s e e S e c t i o n 4 . 1 ) .
5 C o n c l u d i n g R e m a r k s
F o u r t h o r d e r c o m p a c t d i s c r e t i z a t i o n f o r m u l a s h a v e b e e n u s e d i n c o n j u n c t i o n w i t h m u l t i g r i d
t e c h n i q u e s t o s i m u l a t e t h e t w o d i m e n s i o n a l s q u a r e d r i v e n c a v i t y o w a n d m a y b e u s e d t o
1 1
-
8/14/2019 10.1.1.35.2402 (1).ps
12/16
s o l v e m o r e g e n e r a l s t e a d y s t a t e i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n . F o r t h e d r i v e n
c a v i t y p r o b l e m , o u r c o m p u t e d s o l u t i o n s u p t o R e = 7 ; 5 0 0 w i t h c o a r s e r d i s c r e t i z a t i o n s
a r e c o m p a r e d w e l l w i t h t h e b e n c h m a r k s o l u t i o n s o b t a i n e d b y o t h e r r e s e a r c h e r s u s i n g n e r
d i s c r e t i z a t i o n s .
R e f e r e n c e s
1 ] I . A l t a s a n d K . B u r r a g e , A h i g h a c c u r a c y d e f e c t - c o r r e c t i o n m u l t i g r i d m e t h o d f o r
t h e s t e a d y i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s , J . C o m p u t . P h y s . 1 1 4 , 2 2 7 { 2 3 3
( 1 9 9 4 ) .
2 ] S . C . R . D e n n i s a n d J . D . H u d s o n , C o m p a c t h
4
n i t e - d i e r e n c e a p p r o x i m a t i o n s t o
o p e r a t o r s o f N a v i e r - S t o k e s t y p e , J . C o m p u t . P h y s . 8 5 , 3 9 0 { 4 1 6 ( 1 9 8 9 ) .
3 ] U . G h i a , K . N . G h i a a n d C . T . S h i n , H i g h - R e s o l u t i o n s f o r i n c o m p r e s s i b l e o w u s i n g
t h e N a v i e r - S t o k e s e q u a t i o n s a n d a m u l t i g r i d m e t h o d , J . C o m p u t . P h y s . 4 8 , 3 8 7 { 4 1 1
( 1 9 8 2 ) .
4 ] M . M . G u p t a , R . P . M a n o h a r , a n d J . W . S t e p h e n s o n , A s i n g l e c e l l h i g h o r d e r s c h e m e f o r
t h e c o n v e c t i o n - d i u s i o n e q u a t i o n w i t h v a r i a b l e c o e c i e n t s , I n t . J . N u m e r . M e t h o d s
F l u i d s 4 , 6 4 1 { 6 5 1 ( 1 9 8 4 ) .
5 ] M . M . G u p t a , A f o u r t h - o r d e r P o i s s o n s o l v e r , J . C o m p u t . P h y s . 5 5 , 1 6 6 { 1 7 2 ( 1 9 8 4 ) .
6 ] M . M . G u p t a , H i g h a c c u r a c y s o l u t i o n s o f i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s , J .
C o m p u t . P h y s . 9 3 , 3 4 3 { 3 5 9 ( 1 9 9 1 ) .
7 ] M . M . G u p t a , J . K o u a t c h o u a n d J . Z h a n g , C o m p a r i s o n o f s e c o n d a n d f o u r t h o r d e r
d i s c r e t i z a t i o n s f o r t h e m u l t i g r i d P o i s s o n s o l v e r , J . C o m p u t . P h y s . 1 3 2 ( 2 ) , 2 2 6 { 2 3 2
( 1 9 9 7 ) .
8 ] M . M . G u p t a , J . K o u a t c h o u a n d J . Z h a n g , A c o m p a c t m u l t i g r i d s o l v e r f o r t h e
c o n v e c t i o n - d i u s i o n e q u a t i o n , J . C o m p u t . P h y s . 1 3 2 ( 1 ) , 1 2 3 { 1 2 9 ( 1 9 9 7 ) .
9 ] T . Y . H o u a n d B . T . R . W e t t o n , S t a b l e f o u r t h o r d e r s t r e a m - f u n c t i o n m e t h o d s f o r
i n c o m p r e s s i b l e o w s w i t h b o u n d a r i e s , ( p r e p r i n t , 1 9 9 3 ) .
1 0 ] M . L i , T . T a n g a n d B . F o r n b e r g , A c o m p a c t f o u r t h - o r d e r n i t e d i e r e n c e s c h e m e f o r
t h e s t e a d y i n c o m p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s , I n t . J . N u m e r . M e t h o d s F l u i d s
2 0 , 1 1 3 7 { 1 1 5 1 ( 1 9 9 5 ) .
1 1 ] H . N i s h i d a a n d N . S a t o f u k a , H i g h e r - o r d e r s o l u t i o n s o f s q u a r e d r i v e n c a v i t y o w u s i n g
a v a r i a b l e - o r d e r m u l t i - g r i d m e t h o d , I n t . J . N u m e r . M e t h o d s E n g . 3 4 , 6 3 7 { 6 5 3 ( 1 9 9 2 ) .
1 2 ] P . J . R o a c h e , C o m p u t a t i o n a l F l u i d D y n a m i c s , H e r m o s a , A l b u q u e r q u e , N M , 1 9 7 3 .
1 3 ] S . S c h a e r , H i g h o r d e r m u l t i - g r i d m e t h o d s , M a t h . C o m p . 4 3 , 8 9 { 1 1 5 ( 1 9 8 4 ) .
1 2
-
8/14/2019 10.1.1.35.2402 (1).ps
13/16
1 4 ] W . F . S p o t z a n d G . F . C a r e y , H i g h - o r d e r c o m p a c t s c h e m e f o r t h e s t r e a m - f u n c t i o n
v o r t i c i t y e q u a t i o n s , I n t . J . N u m e r . M e t h o d s E n g . 3 8 , 3 4 9 7 { 3 5 1 2 ( 1 9 9 5 ) .
1 5 ] W . F . S p o t z , A c c u r a c y a n d p e r f o r m a n c e o f n u m e r i c a l w a l l b o u n d a r y c o n d i t i o n s f o r
s t e a d y , 2 D i n c o m p r e s s i b l e s t r e a m - f u n c t i o n v o r t i c i t y , I n t . J . N u m e r . M e t h o d s F l u i d s ,
2 8 ( 4 ) , 7 3 7 { 7 5 7 ( 1 9 9 8 ) .
1 6 ] C . P . T z a n o s , H i g h e r - o r d e r d i e r e n c e m e t h o d w i t h a m u l t i g r i d a p p r o a c h f o r t h e s o -
l u t i o n o f t h e i n c o m p r e s s i b l e o w e q u a t i o n s a t h i g h R e y n o l d s n u m b e r s , N u m e r . H e a t
T r a n s f e r : P a r t B , 2 2 , 1 7 9 { 1 9 8 ( 1 9 9 2 ) .
1 7 ] M . C . T h o m p s o n a n d J . H . F e r z i g e r , A n a d a p t i v e m u l t i g r i d t e c h n i q u e f o r t h e i n c o m -
p r e s s i b l e N a v i e r - S t o k e s e q u a t i o n s , J . C o m p u t . P h y s . 8 2 , 9 4 { 1 2 1 ( 1 9 8 9 ) .
1 8 ] P . W e s s e l i n g , M u l t i g r i d s o l u t i o n o f t h e N a v i e r - S t o k e s e q u a t i o n s i n t h e v o r t i c i t y -
s t r e a m f u n c t i o n f o r m u l a t i o n , i n E c i e n t S o l u t i o n s o f E l l i p t i c S y s t e m s , N o t e s o n N u -
m e r i c a l F l u i d M e c h a n i c s , V o l . 1 0 , W . H a c k b u s c h , E d . , V i e w e g , B r a u n s c h w e i g , 1 9 8 4 ,
p . 1 4 5 { 1 5 4 .
1 9 ] N . G . W r i g h t a n d P . H . G a s k e l l , A n e c i e n t m u l t i g r i d a p p r o a c h t o s o l v i n g h i g h l y
r e c i r c u l a t i n g o w s , C o m p u t e r s & F l u i d s , 2 4 , 6 3 { 7 9 ( 1 9 9 5 ) .
2 0 ] P . M . d e Z e e u w , M a t r i x - d e p e n d e n t p r o l o n g a t i o n s a n d r e s t r i c t i o n s i n a b l a c k b o x m u l t i -
g r i d s o l v e r , J . C o m p u t . A p p l . M a t h . 3 3 , 1 { 2 7 ( 1 9 9 0 ) .
2 1 ] J . Z h a n g , A c c e l e r a t e d m u l t i g r i d h i g h a c c u r a c y s o l u t i o n o f t h e c o n v e c t i o n - d i u s i o n
e q u a t i o n w i t h h i g h R e y n o l d s n u m b e r , N u m e r . M e t h o d s P a r t i a l D i e r e n t i a l E q u a t i o n s ,
1 3 , 7 7 - 9 2 ( 1 9 9 7 ) .
2 2 ] J . Z h a n g , O n c o n v e r g e n c e o f i t e r a t i v e m e t h o d s f o r a f o u r t h - o r d e r d i s c r e t i z a t i o n
s c h e m e , A p p l . M a t h . L e t t . 1 0 ( 2 ) , 4 9 { 5 5 ( 1 9 9 7 ) .
2 3 ] J . Z h a n g , M u l t i g r i d A c c e l e r a t i o n T e c h n i q u e s a n d A p p l i c a t i o n s t o t h e N u m e r i c a l S o -
l u t i o n o f P a r t i a l D i e r e n t i a l E q u a t i o n s , P h . D . T h e s i s , T h e G e o r g e W a s h i n g t o n U n i -
v e r s i t y , 1 9 9 7 .
2 4 ] J . Z h a n g , O n c o n v e r g e n c e a n d p e r f o r m a n c e o f i t e r a t i v e m e t h o d s w i t h f o u r t h - o r d e r
c o m p a c t s c h e m e s , N u m e r . M e t h o d s P a r t i a l D i e r e n t i a l E q u a t i o n s , 1 4 , 2 6 3 { 2 8 0
( 1 9 9 8 ) .
2 5 ] J . Z h a n g , R e s i d u a l s m o o t h i n g t e c h n i q u e s i n m u l t i g r i d , I : e q u i v a l e n c e p r o o f , A p p l .
M a t h . C o m p u t . 8 6 ( 2 - 3 ) , 2 8 3 { 3 0 3 ( 1 9 9 7 ) .
2 6 ] J . Z h a n g , R e s i d u a l s c a l i n g t e c h n i q u e s i n m u l t i g r i d , I I : p r a c t i c a l a p p l i c a t i o n s , A p p l .
M a t h . C o m p u t . 9 0 ( 2 - 3 ) , 2 2 9 { 2 9 2 ( 1 9 9 8 ) .
1 3
-
8/14/2019 10.1.1.35.2402 (1).ps
14/16
2 7 ] J . Z h a n g , A n o t e o n a n a c c e l e r a t e d h i g h a c c u r a c y m u l t i g r i d s o l u t i o n o f t h e c o n v e c t i o n -
d i u s i o n e q u a t i o n w i t h h i g h R e y n o l d s n u m b e r , N u m e r . M e t h o d s P a r t i a l D i e r e n t i a l
E q u a t i o n s , 1 6 ( 1 ) , 1 { 1 0 ( 2 0 0 0 ) .
1 4
-
8/14/2019 10.1.1.35.2402 (1).ps
15/16
0 0.2 0.4 0.6 0.8 10
0.2
0.4
0.6
0.8
1
0 0.2 0.4 0.6 0.8 10
0.2
0.4
0.6
0.8
1
F i g u r e 2 : S t r e a m l i n e a n d e q u i v o r t i c i t y c u r v e s f o r R e = 5 ; 0 0 0 w i t h h = 1 = 1 2 8 .
1 5
-
8/14/2019 10.1.1.35.2402 (1).ps
16/16
0 0.2 0.4 0.6 0.8 10
0.2
0.4
0.6
0.8
1
0 0.2 0.4 0.6 0.8 10
0.2
0.4
0.6
0.8
1
F i g u r e 3 : S t r e a m l i n e a n d e q u i v o r t i c i t y c u r v e s f o r R e = 7 ; 5 0 0 w i t h h = 1 = 1 2 8 .
1 6