I n d on e s i an   Jo u r n al   o El e c t r i c al   En gi n e e r i n g   an d   C o m p u te r   S c i e n c e   V o l .   16 ,   N o .   2 N o v e m b e r   201 9 ,   pp.   101 8 ~ 10 25   IS S N :   2502 - 4752 ,   D O I :   10. 1 1591 / i j e e c s . v 1 6 .i 2 . pp 101 8 - 1025             1018       Jou r n al   h o m e pa ge ht t p: / / i ae s c or e . c om / j our na l s / i nde x . php/ i j e e c s   C o m p u t a t i o n a l   st u d y   o f   f l o w   a r o u n d   a   N A C A   0 0 1 2   b y   u si n g     R o e   FV M   sch e m e   a n d   d a v i s - y e e   T V D   schem e       F ati m ah   Y u s o p 1 ,   Zam r i   O m ar 2 ,   B am b an B as u n o 3 ,   N i k   N o r mu n i r M at   H as s an 4   1, 4 F a c ul t y   o f   E ng i ne e r i ng   T e c hno l o gy ,   U ni v e r s i t i   T un   H us s e i n   O nn   M a l a y s i a ,   M a l a y s i a   2 ,3 F a c ul t y   o f   M e c ha ni c a l   a nd   M a nuf a c t u r i ng   E ng i ne e r i ng ,   U ni v e r s i t i   T un  H us s e i O n M a l a y s i a ,   M a l a y s i a       A r ti c l e   I n fo     A B S TR A C T   Ar t i c l e   h i s t or y :   R e c e i v e J a n   12 ,   2 019   R e v i s e   M a r   10 ,   2019   A c c e pt e A p r   15 ,   201 9       C ur r e nt l y   C F D   ha b e e n   c o ns i de r e d   a s   a n   i m po r t a n t   t o o l   f or   s o l v i ng   e ng i ne e r i ng   pr o bl e m s .   T he   a p pl i c a t i o o f   C F D   ha be e us e i n t e ns i v e l y   i a i r c r a f t   i ndu s t r i e s   i n   de s i g a   ne w   a i r c r a f t   o r   i t he   e f f o r t   o f   i m pr o v e m e nt   o t he   e x i t i ng   a i r c r a f t .   I t e r m   of   C F D   c o m put e r   c o de ,   t he   C F D   c ode   di f f e r s   w i t a ny   o t he r s   m a y   due   t o   t he   di f f e r e nc e   i t h e   num e r i c a l   s c he m e   ha v e   be e us e d .   T he r e f o r e ,   t he   p r e s e n t   w o r pr e s e nt s   t h e   c o m pa r i s o r e s u l t   be t w e e n   t w o   de v e l o pe c o m put e r   c o de s   w i t h   A N S Y S - F L U E N T   s o f t w a r e   t o   t he   c a s e   o f   t r a n s o ni c   s t e a dy   f l o w   pa s t   t hr o ug a i r f o i l   N A C A   0012.   T he   f i r s t   c o m put e r   c o de   us e a   f i n i t e   d i f f e r e nc e   m e t ho w i t num e r i c a l   s c he m e   a c c o r di ng   t o   D a v i s - Y e e   T V D   s c he m e .   M e a nw h i l e ,   t h e   s e c o nd  c o m put e r   c o de   us e a   R o e s   c e l l   c e nt r e   f i n i t e   v o l um e   s c he m e .   T h e   f l o w   a na l y s i s   i s   c a r r i e o ut   a t   t w o   M a c num be r ,   M   ( 0. 6 &   0. 8) .   E a c M a c num b e r   a p pl i e t o   t w di f f e r e n t   a ng l e s   o f   a t t a c ks   (   &   ) .   T h e   f l o w   do m a i d i s c r e t i z e by   us e   of   C - t o po l ogy   w i t 193 x63  g r i po i n t s .   T he   c o m pa r i s o i t e r m   o f   t he   pr e s s ur e   c o e f f i c i e nt ,   a l o ng   t he   a i r f o i l   s ur f a c e   a r e   p r e s e n t e d .   F r o m   t h e   r e s u l t ,   i nd i c a t e d   t ha t   d e v e l o pe c o m put e r   c o de   i s   a b l e   t o   c a pt u r e   t he   pr e s e nc e   o f   s ho c w a v e   i t h e   f l o w   f i e l d .   Ke y w or d s :   Co m put a t i o na l   f l ui dy n a m i c s   upw i n d   s c h e m e   C opy r i gh t   ©   201 9   I n s t i t ut e   o f   A dv anc e E ng i ne e r i ng   and   S c i e nc e .     A l l   r i gh t s   r e s e r v e d .   Cor r e s pon di n g   Au t h or :   F a t i m a Y us o p,     F a c ul t y   of   E n g i n e e ri n T e c hn o l o gy ,   U n i v e r s i t i   T u n   H us s e i O nn  M a l a y s i a ,     P a go h   H i g h e r   E duc a t i o H ub ,   84600   P a go h ,   M u a r ,   J o h o r.   E m a i l :   f a t i m a h y @ ut h m . e du. m y       1.   I N TR O D U C TI O N     Co m put a t i o na l   f l ui dy n a m i c s   (CF D i s   ve r y   w i de s pr e a us e   i n   e n gi n e e r i n a c t i v i t i e s   w h i c i n v o l ve s   i n   s o l v i n f l ui f l ow   pr o b l e m s .   T h e   go ve rn i n e q ua t i o n s   o f l ui m o t i o n   c a l l e t h e   N a v i e r - S t o ke s   e qua t i o n s   a r e   w e l l   e s t a b l i s h e n e a rl y   o n e   a n ha l f   c e n t u r i e s   a go   [1].   T i n o c o   e s t i m a t e t ha t   di r e c t   s o l ut i o n   t o   t h e   N a v i e r - S t o ke s   i n   l i n e   w i t h   p r o g r e s s   i n   t h e   de v e l o pm e n t   of   c o m put i n g   po w e r   a nd  t h e   a l go r i t h m   i n   s o l v i n t h e   n o n - l i n e a r   p a r t i a l   di f f e r e n t i a l   e qua t i o n   c a n   b e   m a de   i t h e   n e xt   60  y e a r s   f r o m   n o w   [2].   H ow e ve r   f o r   pa r t i c ul a r   f l o w   pr o b l e m   e s pe c i a l l y   fo r   t h e   c a s e   of   h i g h   R e y n o l ds   n um b e r   f l o w   pa s t   t hr o ug h   a   s t r e a m l i n e   b o d y   a t   r e l a t i v e l y   l ow   a n g l e   o f   a t t a c k,   t h e   go v e r ni n e qu a t i o n   f o r   s uc h   k i n f l o w   c a n   b e   w e l l   p r e s e n t e by   t h e   E ul e r   e qu a t i o n s   [ 3,   4] .     F o r   t h e   c a s e   of   t r a n s o ni c   s t e a dy   f l ow   pr o b l e m s ,   t h e   E u l e r   e qua t i o n   c a n   b e h a v e   a s   E l l i pt i c   pa rt i a l   di f fe r e nt i a l   e qua t i o n   o r   h y pe r - bo l i c   pa rt i a l   di f f e r e n t i a l   e qua t i o n   de pe ndi n o n   t h e   l o c a l   M a c h   num b e r   [5] .     In   t h e   p r e s e n c e   s uc h   m i xe t y p e s   be h a v i o r ,   m a ke s   t h e   E ul e r   e qua t i o n   b e c o m e s   di ff i c ul t i e s   t b e   s o l ve d.   S uc h   di f f i c ul t y   c a n   b e   r e m ov e by   t r e a t i n t h e   go v e r ni n e qua - t i o n   o f l ui m o t i o n   i n   t h e   fo r m   t i m e   de pe n de n t   e qua t i o n .   I n   t h i s   r e s pe c t   t h e   E u l e e qua t i o b e h a v e s   a s   h y pe r b o l i c   e qua t i o w i t h   r e s pe c t   t o   t i m e   o v e r   t h e   w h o l e   f l ow   d o m a i n.   T h e r e   a r e   v a r i o us   m e t h o ds   h a b e e n   de v e l o p e fo r   s o l v i n t h e   E u l e r   e qua t i o n   i c o r r e s po n di ng  t o   t h e   t r a n s o ni c   s t e a dy   f l ow   pr o b l e m .     Evaluation Warning : The document was created with Spire.PDF for Python.
In do n e s i a J   E l e c   E ng  &   Co m S c i     IS S N :   2502 - 4752       Com put a t i on al   s t ud y   o f   f l ow   ar ound  N A C A   0012   by   us i n R oe   F V s c h e m e   and   …  ( F at i m ah   Y us op )   1019   T h e   p r e s e n t   w o r de a l s   w i t h   t h e   f l o w   pr o b l e m   r e l a t e t o   s t re a m - l i n e   b o d y   of   a   r e l a t i v e l y   l ow   a n g l e   of   a t t a c i n   w hi c h   t h e   v i s c o u s   e ff e c t   c a n   b e   n e gl e c t e d.   I n   s ol v i n t h i s   t y pe   o f   f l ow ,   t h e   t w o   c o m put e r   c o de s   a r e   de ve l o pe a n de s i g n e d.   T h e   f i r s t   c o m put e r   c o de   a r e   de ve l o pe b a s e o n   f i n i t e   di f f e r e n c e   m e t h o d   a c c o r di n t o   D a v i s - Y e e   T V D   s c h e m e   [6 ,   7] ,   w hi l e   t h e   s e c on d   c o m put e c o de   us e s   a   f i n i t e   v o l um e   a pp r o a c by   a do pt i n g   t h e   s c h e m e   i n t r o duc e by   Roe   [8,   9].       2.   R ES EA R C H   M ET H O D   T h e   m e t h o do l o g y   of   t h i s   s t udy   i s   b a s e o n   CF D   c o de   de ve l o pm e n t .   B a s i c a l l y ,   CF D   c o de   i n v o l v e s   t hr e e   pr o c e s s ,   c a l l e pr e - p r o c e s s i n g,   s o l v e r   a n po s t   pr o c e s s i n w i t h   s e que n c e   a s   de pi c t e i n   F i gu r e   1.     P r e - p r o c e s s i n i n v o l ve s   t h e   w a y   t h e   f l ow   d o m a i n   o f   t h e   fl o w   pr o b l e m   un de r   i n v e s t i ga t i o n   i s   di s c r e t i z e d.   S o l ve r   r e l a t e i n   m a nn e o f   t h e   go v e r n i ng  e qua t i o o f   f l ui m o t i o n   i s   s o l v e d.   T h e   po s t   p r o c e s s i n r e l a t e i t h e   r e s ul t   o f   f l ow   pr o pe r t i e s   pr e s e nt e [10] .   I n   t e xt   o f   s o l ve r ,   w hi c i n v o l v e s   h o w   t o   s o l ve   t h e   go v e r n i ng  e qua t i o n   o f   f l ui m o t i o n,   i t   c a us e   a   n u m e ri c a l   s c h e m e   d e ve l o p e a c c o r di n t o   f i n i t e   di f f e r e n c e   m e t h o (F D M o r   f i n i t e   v o l um e   m e t h o (F V M ).   T h e   p r e s e n t   w o r i n v o l v e s   t h e   s o l v e r   de v e l o pe b a s e o n   F D M   a s   w e l l   a s   F V M .   M e a n w hi l e ,   F i gu r e   s h o w s   t h e   ge n e ra l   f l ow   i n   s o l v e r   c o de .   T h e   pr o c e s s e s   s t a r t   w i t h   de f i ni n t h e   s i m ul a t i o n   pa ra m e t e r   u n t i l   ge n e ra t e s   t h e   f i na l   r e s ul t .             F i gu r e   1 .   G e n e ra l   f l o w   of   CF D   c o de   de v e l o pm e nt           F i gu r e   2 .   F l o w   of   CF D   s o l v e r       Evaluation Warning : The document was created with Spire.PDF for Python.
                               IS S N :   2502 - 4752   In do n e s i a J   E l e c   E ng  &   Co m S c i ,   V o l .   16 ,   N o .   2 N o v e m be r   2 019   :     10 18   -   10 25   1020   2. 1 .   S tr u c tu r e d   G r i d   G e n e r at i o n     T h i s   s t udy   w a s   i m pl e m e nt e C - t o po l o g y   i n   de f i n i ng  t h e   f l ow   d o m a i n   pa s t   t hr o ug h   a i r f o i l   a n h a v e   b e e n   us e by   t h e   E l l i p t i c   g r i ge n e r a t o r   a nd  a l ge b r a i c   ge n e ra t o r   a s   i t s   g r i i n i t i a l i z a t i o n   w i t h   s p a t i a l   di s c r e t i z a t i o n   i n   ξ   a nd  η  di r e c t i o n   w i t h   193  po i n t s   a nd  63  po i n t s   r e s pe c t i v e l y   gi v e s   t h e   r e s ul t   a s   s h o w n   i n   F i gu r e   3 .           F i gu r e   3 .   C - g r i d   o f   N A CA   0012       2. 2 .     G o v e r n i n Eq u ati o n     F o r   t h e   c o m pr e s s i b l e   f l ow   pa s t   t hr o ug h   a   s t r e a m l i n e   b o d y   a t   r e l a t i v e l y   l ow   a n gl e   o f   a t t a c a nd  hi g R e y n o l num b e r   c a b e   a pp r o xi m a t e a s   o n e   s o l v e s   a i n v i s c i f l ow   pr o b l e m .   H e n c e   s uc h   f l ow   ph e n o m e na   c a n   b e   w e l l   pr e s e n t e by   t h e   E ul e r   e qua t i o n s .   T h e   u n s t e a d y ,   t w di m e n s i o n a l   c o m pr e s s i b l e   i n v i s c i f l ow s   m a ke   t h e   E ul e r   e qua t i o n s   i n   di f f e r e n t i a l   a nd  i n   c o n s e r v a t i v e   fo r m   a n v e c t o r   n o t a t i o n   c a n   b e   w r i t t e   a s   [11 ,   12]:                                 (1)     In   a b ov e   e qua t i o n   Q   i s   a   s t a t e   v e c t o r   of   d e pe n de nt   v a ri a b l e s   w h i l e   E   a n F   a r e   t h e   f l ux  v e c t o r   i n   x   a n y   di r e c t i o r e s pe c t i v e l y .   T h e s e   t hr e e   v e c t o r s   a r e   de f i n e d   a s :         [               ]       [                          (           ) ]               [                               (           ) ]           (2)     a n d                             (3)     T h e   v a r i a b l e s   a r e   a ppe a ri n i n   ( 2 )   ha v e   de f i n i t i o n   a s   c om m o n l y   us e i n   f l ui dy n a m i c s   s t udy .     T h e   v a ri a b l e   ρ  r e p r e s e nt s   t h e   a i r   de n s i t y ,   a n d   v   a r e   t h e   c o m po n e n t   v e l oc i t i e s   i n   x   a n d   y   di r e c t i o n.   e t   i s   t h e   t o t a l   i nt e rn a l   e n e r gy   pe r   uni t   o f   m a s s   a n i s   t h e   p r e s s ur e .   Co n s i de r i n ( 2 ) ,   i t   i s   c l e a rl y   i n di c a t i n t h a t   ( 1 )   c o n s i s t   f o ur   e qua t i o n s   w i t h   f i v e   de pe n de n t   v a r i a b l e s .   O n e   a ddi t i o n a l   r e l a t i o n s h i i s   r e qui r e t o   m a ke   ( 1 )   a s   s o l v a b l e   e qua t i o n .   By   c o n s i de ri n t h e   f l ow   b e h a v e s   a s   a   pe r f e c t   ga s ,   t h e   a ddi t i o n a l   r e l a t i o n s h i c a n   b e   de f i n e f r o m   t h e   r e l a t i o n s h i p   b e t w e e n   t h e   i n t e rna l   e n e r gy   pe u n i t   o f   m a s s   e   a n d   o t h e f l o w   v a r i a b l e s   a s :             (       )   (4)     In  i nt e g r a l   f o r m ,   t h e   E ul e r   e qu a t i o n   c a b e   p r e s e n t e a s   [ 13 ]:     Evaluation Warning : The document was created with Spire.PDF for Python.
In do n e s i a J   E l e c   E ng  &   Co m S c i     IS S N :   2502 - 4752       Com put a t i on al   s t ud y   o f   f l ow   ar ound  N A C A   0012   by   us i n R oe   F V s c h e m e   and   …  ( F at i m ah   Y us op )   1021                        (                   )              (5)   In  a b o ve   e qua t i o n ,     i s   a   c e l l   v o l um e ,   n x   a nd  n y   a r e   t h e   c o m po n e nt   o f   t h e   n o rm a l   u n i t y   v e c t o r   o t h e   c e l l   f a c e s .   S   i s   t h e   c e l l   s u r f a c e   a r e a .   W hi l e   v a r i a b l e   Q ,   E   a n F   ha v e   t h e   s a m e   de f i n i t i o a s   i t h e   E ul e r   e qua t i o n   i n   i t s   di f f e r e n t i a l   f o r m .     2. 3 .     TV D   S c h e m e   T V D   s c h e m e   i s   a   h i g h - r e s o l ut i o n   s c h e m e   t ha t   w a s   de v e l o pe t o   s o l v e   t h e   p r o b l e m   r e l a t e t o   n o n l i n e a r   h y pe r bo l i c   c o n s e r v a t i o n   l a w s .   T h e   t e r m   T V D   w a s   i n t r o duc e by   H a r t e n   i n   19 83  [14] .   O v e r a l l ,   T V D   s c h e m e   i s   n o t   a s   s pe c i f i c   a s   a   s c h e m e   b ut   m o r e   a s   pr o pe r t i e s .   T h e r e f o r e ,   T V D   s c h e m e   c o n t a i n s   a   v a r i e t y   of   s c h e m e .   T ha t   m e a n s ,   T V D   s c h e m e   c a n   b e   s um m a ri z e d   t o   h a v e   t h e   fo l l ow i n p h y s i c a l   c h a ra c t e ri s t i c s     [15 - 17] .   a)   E n t r o py   c o n di t i o   A   de c r e a s e   o f   e n t r o py   a s s oc i a t e w i t h   e xpa n s i o n   s h o c ks   m us t   n o t   b e   a d m i t t e d.   b)   M o n o t o n i c   c o n di t i o   T hi s   c o n d i t i o n   m us t   b e   e n f o r c e t o   p re ve n t   o s c i l l a t o r y   b e h a v i o r   i t h e   n u m e r i c a l   s c h e m e .     c)   T o t a l   V a ri a t i o D i m i ni s h i ng  ( T V D   T h e   t o t a l   v a ri a t i o o f   a n y   ph y s i c a l l y   a dm i s s i b l e   s o l ut i o m us t   n o t   b e   a l l o w e t o   i n c r e a s e   t i m e .   B a s i c a l l y ,   c a n   b e   de f i n e a s   f o l l ow :          |       |          T hr o ug h   a b o ve   c r i t e r i a ,   i t   c a de v e l o a   m a nn e r   o n   h o w   t o   di s c r e t i z e   t h e   go v e r n i ng  e qua t i o o f   f l ui m o t i o n.   A m o n t h e m   a r e   T V D   s c h e m e   a c c o r di n t o   [18 - 20].   T h e   CF D   c o de   w h i c h   de v e l o pe b a s e o n   T V D   c r i t e ri a   ha b e e n   w i de l y   us e a s   CF D   c o de   fo r   a i r c r a f t   a e ro d y n a m i c s   a na l y s i s   [21].   In   t e rm   o f   f l ui f l o w   a n a l y s i s ,   T V D   s c h e m e   ha b e e n   us e e xt e n s i v e l y   e s pe c i a l l y   f o r   i n t e rna l   f l ow   pr o b l e m   s uc h   i n   [2 2,   2 3]  a nd  [24].   T h e r e f o r e ,   t h i s   s t u dy   de ve l o pe t h e   T V D   s c h e m e   b a s e o n   D a v i s - Y e e   s c h e m e   f o r   s o l ve   e xt e r na l   f l o w   t h a t   a b l e   t o   s o l v e   t h e   a e r o d y n a m i c s   a n a l y s i s   a s   s h o w n   i n   ( 6 ) .       ̅                 ̅                       [ (     )               (     )                 (     )               (     )             ]             [ (     )               (     )                 (     )               (     )               ]   (6)     In  a b o ve   e qua t i o n   Ф   i s   a   f l ux  l i m i t e v e c t o r   a n d   G   i s   l i m i t e r   f un c t i o n.     2. 3 .     R O S c h e m e   T h e   s e c o n c o de   i s   de ve l o pe by   R O E   s c h e m e .   T hi s   s c h e m e   i s   b a s e o n   a   c ha ra c t e ri s t i c   de c o m po s i t i o n   o f   t h e   f l ux  di f fe r e n c e s   [2 5].   T h e   R oe   s c h e m e   us e f l ux  fo r m u l a   a t   t h e   i n t e r f a c e   of   a   c o n t r o l   vo l um e   t h a t   i s   e qua l   t o   t h e   a v e r a ge   f l uxe s   of   l e f t   a n ri g h t   s t a t e s   m i n us   a   d i f f e r e n c i ng  t e r m   w hi c h   s p l i t s   t h e   di f fe r e n c e   o f   t h e   f l uxe s   o b o t s i de s   o f   t h e   c o nt r o l   v o l um e .   T h e   b a s i c   e qua t i o i s   e xp r e s s e a s   b e l ow   [8].     (     )     (     )     (   ̅      )         (           )   (7)     w h e r e   A R o e   i s   Ro e   m a t ri a n d   L   a nd  R   t h e   l e f t   a nd  ri g ht   s t a t e   r e s pe c t i v e l y .       3.   R ES U LTS   A ND  D IS C U S S I O N   T h e   s i m ul a t i o n   w a s   c o n duc t e i n   f o ur   di f f e r e n t   f l ow   c o n di t i o n s   w i t h   v a r y i n t h e   M a c h   n um b e   (0. a n d   0 . 65)  a n d   a ngl e   o f   a t t a c (0°   a nd  5 ° ).   Co m pa ri s o n   r e s ul t s   a r e   p r e s e n t e b e t w e e n   t h e   t w o   de v e l o pe d   c o d e s   (R O E   F V M   a n d   T V D   F D M a n d   w i t h   c o m m e r c i a l   c o de   A N S Y S - F L U E N T .     3. 1 .     A i r   F l o w   o n   N A C A   001 A i r fo i l   fo r   M =0 . 65   an d   α =0 °   T h e   c o m pa r i s o r e s ul t   i t e r m s   o f   pr e s s u r e   c o e ff i c i e n t   d i s t ri b ut i o n   a l o n t h e   a i r f o i l   s u r f a c e   i s   a s   s h o w n   i n   F i gu r e   4 .   B o t h   de v e l o p e c o m put e r   c o de s   a r e   i n   g oo a gr e e m e n t   w i t h   e a c h   o t h e r   c o m pa r e t o   t h e   ANSYS - F L U E N T   s of t w a r e   r e s ul t s   di s c r e p a n c i e s   a ppe a ri n a t   x   <   0. 4   c .   T h e   r e s ul t   a l s o   i n d i c a t e t h e   v a l ue   fo r   uppe r   s u r f a c e   a n l o w e r   s u r f a c e   a t   t h e   s a m e   po i n t   due   t o   t h e   s y m m e t r i c a l   a i r f o i l   w i t h     a n g l e   o f   a t t a c k.   In   t e rm s   of   pr e s s ur e   c o e ff i c i e n t   v a l ue   a t   s t a g n a t i o n   po i n t ,   de v e l o pe c o de   i s   n e a r l y   1. c o m pa r e t o   A N S Y S - F L U E N T   s of t w a r e   ob t a i n e a r o u n 0 . 1 .   T h e   c o m pa ri s o n   r e s ul t s   i n   v i e w   of   M a c h   n u m b e r   i s   s h o w n   i n     F i gu r e   5 .       Evaluation Warning : The document was created with Spire.PDF for Python.
                               IS S N :   2502 - 4752   In do n e s i a J   E l e c   E ng  &   Co m S c i ,   V o l .   16 ,   N o .   2 N o v e m be r   2 019   :     10 18   -   10 25   1022       F i gu r e   4 .   P r e s s ur e   c o e ff i c i e n t   d i s t r i b ut i o a t   α =     a n d   M   =   0 . 65         (a )     (b )     F i gu r e   5 .   T h e   M a c n u m b e r   c o n t o u r   n e a r   t h e   N A CA   0012 ,   a t   α=     a nd  M   =   0 . 65   t hr o ug (a )   T V D   F D M   a n d   (b Ro e   F V M       3. 2 .     A i r   F l o w   o n   N A C A   001 A i r fo i l   fo r   M =0 . 65   an d   α =5 °   A   c o m pa r i s o n   r e s ul t   b e t w e e n   t h e   de v e l o p e c o m put e r   c o de s   a n A N S Y S - F L U E N T   s of t w a r e   fo r   pr e s s u r e   c o e ff i c i e n t   a l o ng  t h e   a i r f o i l   s ur f a c e   a t   di ff e r e n t   a n g l e s   of   a t t a c i s   s h o w n   i n   F i gu r e   6.   A t   t hi s   a n gl e   of  a t t a c k,   b o t h   de v e l o pe c o m put e r   c o de s   (F D M   Co de   a n d   F V M   Co de pr e di c t   t h e   p r e s e n c e   o f   t h e   s h o c w a ve ,   w h i l e   A N S Y S - F L U E N T   s of t w a r e   do e s   n o t   i n d i c a t e   s uc h   p h e n o m e na .   A s   a   r e s ul t ,   t h e   p r e s s u r e   di s t r i b ut i o di s t r i b ut i o n   a l o n t h e   a i r f o i l   s ur f a c e   di ff e r s .   H ow e ve r ,   t h e   di f fe r e n c e   be t w e e n   t h e   T V D   F D M   c o de   a n t h e   Ro e   F V M   c o de   a r e   a s   n o t   m uc a s   t h e   d i f fe r e n c e   w i t h   A N S Y S - F L U E N T   s o f t w a r e .   T h e   di f f e r e n c e   o f   T V D   F D M   c o d e   a n d   R o e   F V M   c o d e   a r e   s i m pl y   due   t o   t h e i r   di f f e r e n c e   i p r e di c t i ng  t h e   l o c a t i o n   o f   s h o c w a ve .   H ow e ve r ,   b o t h   s c h e m e s   a r e   a b l e   t o   pr o duc e   a   s h a r p r e s s u re   gr a d i e nt   w i t h o ut   w i ggl e s   o r   o s c i l l a t i o n   a t   t h e i r   s h o c po i nt s .   T a b l e   1   s h o w s   t h e   s h o c po s i t i o o n   t h e   uppe s ur f a c e   o f   t h e   a i r f o i l   f o r   t h e s e   t hr e e   CF D   c o de s .   T h e   f l ow   pa t t e rn  o f   M a c h   n um b e r   di s t ri b ut i o n   o v e r   t h e   f l ow   f i e l fo r   t h e   T V D   F D M   a nd  R o e   F V M     i n   F i gu r e   7 .       T a b l e 1.   Co m pa ri s o n   r e s ul t   s h o c   po s i t i o n   o f   N A CA   0012   a t   α=       a n M   =   0 . 65   S c h e m e   S h o c k   L o c a t i o n   o U p p e r   S u r f a c e   T V D   F D M   0 . 2 2   c   Ro e   F V M   0 . 3   c   A N S Y S - F L U E N T   s o f t w a r e   N o   s h o c k   w a v e         F i gu r e   6 .   P r e s s ur e   c o e ff i c i e n t   d i s t r i b ut i o a t   α =     a n d   M   =   0 . 65   Evaluation Warning : The document was created with Spire.PDF for Python.
In do n e s i a J   E l e c   E ng  &   Co m S c i     IS S N :   2502 - 4752       Com put a t i on al   s t ud y   o f   f l ow   ar ound  N A C A   0012   by   us i n R oe   F V s c h e m e   and   …  ( F at i m ah   Y us op )   1023     (a )     (b )     F i gu r e   7 .   T h e   M a c n u m b e r   c o n t o u r   n e a r   t h e   N A CA   0012 ,   a t   α=     a nd  M   =   0 . 65   t hr o ug   (a T V D   F D M   a n d   (b R o e   F V M       3. 3 .     A i r   F l o w   o n   N A C A   001 A i r fo i l   fo r   M =0 . 8   an d   α =0 °   W h e n   t h e   M a c h   n u m b e r   w a s   i n c r e a s i n t o   0. a n d   t h e   a n gl e   o f   a t t a c w a s   ke pt   a t   α  =   0 ° ,     t hr e e   c o m put e c o de s   pr o duc e   t h e   pr e s e n c e   of   t h e   s h o c w a v e   a s   i n di c a t e i n   F i gu r e   8.   T h e   R oe   F V M   c o d e   pr e di c t s   t h e   po s i t i o n   o f   s h o c w a v e   n e a rl y   a t   t h e   s a m e   l o c a t i o a s   o b t a i n e by   A N S Y S - F L U E N T   s of t w a r e .   W h i l e   t h a t ,   T V D   F D M   c o de   i s   s l i gh t l y   f ur t h e r   do w n s t re a m .   H ow e ve r ,   i n   t e r m s   o n   h o w   t h e   p r e s s u r e   c oe ff i c i e n t   w a s   di s t ri b ut e a l o n t h e   a i r f o i l   s ur f a c e ,   t h e   t w de ve l o p e c o m put e r   c o de s   pr o duc e n e a rl y   t h e   s a m e   r e s ul t   b ut   di f f e r   s i g ni f i c a nt l y   w i t h   A N S Y S - F L U E N T   s o f t w a r e .   T hi s   i s   due   t o   t h e i di f f e r e n c e   i pr e di c t i n t h e   s h o c l o c a t i o n   a s   s h o w n   i n   T a b l e   2.   I t e r m s   o f   c a pa b i l i t y   i n   p r o duc i n s ha r p r e s s u r e   g r a d i e n t   a t   t h e   s h o c k,   t h e   de ve l o p e c o m put e r   c o de s   a r e   b e t t e r   t ha n   t h e   r e s ul t   pr o duc e by   A N S Y S - F L U E N T   s of t w a r e .   I n   t e rm s   o f   t h e   M a c h   n u m b e r ,   c o m pa ri s o n   r e s ul t s   a r e   s h o w n   i n   F i gu r e   9 .   T h e s e   t w o   f i gur e s   c l e a r l y   i n di c a t e   t ha t   t h e   T V D   F D M   c o de   a n R o e   F V M   c o de   gi v e   n e a rl y   s i m i l a r   f l o w   pa t t e rn.   H ow e v e r ,   t h e   s h o c ks   pr o duc e by   T V D   F D M   c o de   w e r e   b i gge t ha n   R o e   F V M   c ode .       T a b l e   2 .   Co m p a r i s o r e s ul t   s h o c po s i t i o o f   N A CA   0012  a t   α=   0 °   a n d   M   =   0 . 8   S c h e m e   S h o c k   l o c a t i o n   T V D   F D M   U p p e r   s u rf a c e     0 . 6 8   c   L o w e s u r f a c e   0 . 6 6   c   Ro e   F V M   U p p e r   s u rf a c e   0 . 5 2   c   L o w e s u r f a c e   0 . 5 0   c   A N S Y S Y - F L U E N T   s o ft w a re   U p p e r   s u rf a c e   0 . 5 1   c   L o w e s u r   a c e   0 . 5 1   c     - 1 . 2 - 1 . 0 - 0 . 8 - 0 . 6 - 0 . 4 - 0 . 2 0 . 0 0 . 2 0 . 4 0 . 6 0 . 8 1 . 0 1 . 2 0 0 . 2 0 . 4 0 . 6 0 . 8 1 Cp x /c R O E   FV M T V D F DM FL U E N T   2 n d   O r d e r   U p w in d     F i gu r e   8 .   P r e s s ur e   c o e ff i c i e n t   d i s t r i b ut i o f o r   N A CA 0012,   a t   α =     a nd  M   =   0 . 8         (a )     (b )     F i gu r e   9 .   T h e   M a c n u m b e r   c o n t o u r   n e a r   t h e   N A CA   0012 ,   a t   α=     a nd  M   =   0 . 8   t hr o ug   (a T V D   F D M   a n d   (b R o e   F V M   Evaluation Warning : The document was created with Spire.PDF for Python.
                               IS S N :   2502 - 4752   In do n e s i a J   E l e c   E ng  &   Co m S c i ,   V o l .   16 ,   N o .   2 N o v e m be r   2 019   :     10 18   -   10 25   1024   3. 4 .     A i r   fl o w   o n   N A C A   0012   ai r fo i l   fo r   M =0 . an d   α =5°   D i f fe r e n t   a ngl e s   o f   a t t a c ks   b ut   s a m e   M a c h   n u m b e r   a s   p r e v i o us l y   a r e   a pp l i e f o r   t h e   c a s e   o f   f l ow   pa s t   t hr o ugh   N A CA   001 2.   I n c r e a s i ng  a n g l e   of   a t t a c f r o m   α  =   0 t o   be c o m e   α  =   ,   h a b e e n   fo un t o   c a us e   t h e   s h o c w a v e   t o   oc c ur   a t   t h e   po s i t i o n   f u r t h e r   do w n s t r e a m   c o m pa r e t o   po s i t i o n   of   s h o c w a v e   a t   α  =   00.     In   t e rm s   o f   pr e s s ur e   c oe ff i c i e n t   a l o n t h e   a i r f o i l   s ur f a c e ,   t h e   c o m pa r i s o n   r e s ul t s   b e t w e e n   t hi s   t hr e e   s of t w a r e   a r e   s h o w n   i n   F i gu r e   10.   It   i s   n e c e s s a r y   t o   b e   n o t e t h a t   t h e   A N S Y S - F L U E N T   s of t w a r e   w h i c h   i ni t i a l l y   r u n   b y   us e   of   s e c o n o r de r   u pw i n s c h e m e   i s   f o un t o   n o t   w o r k.   T he   A N S Y S - F L U E N T   s of t w a r e   n e e ds   t o   us e   a   f i r s t   o r de upw i nd  s c h e m e .   I t hi s   r e s pe c t   o n e   c a i de nt i fy   t ha t   A N S Y S - F L U E N T   s o f t w a r e   i s   l a c ki ng  c o n s i s t e n c y   i n   t h e i r   num e r i c a l   s c h e m e   i n   us e       T a b l e   3 .   Co m p a r i s o r e s ul t   s h o c po s i t i o o f   N A CA   0012  a t   α=   5 °   a n d   M   =   0 . 8   S c h e m e     S h o c k   l o c a t i o n     T V D   F D M     U p p e r   s u rf a c e     0 . 8 0   c   L o w e s u r f a c e   n o   s h o c k   Ro e   F V M   U p p e r   s u rf a c e   0 . 8 3   c   L o w e s u r f a c e   n o   s h o c k   A N S Y S - F L U E N T   s o f t w a r e     U p p e r   s u rf a c e   0 . 6 3   c   L o w e s u r f a c e   n o   s h o c k       - 1 . 6 - 1 . 4 - 1 . 2 - 1 . 0 - 0 . 8 - 0 . 6 - 0 . 4 - 0 . 2 0 . 0 0 . 2 0 . 4 0 . 6 0 . 8 1 . 0 1 . 2 1 . 4 0 0 . 2 0 . 4 0 . 6 0 . 8 1 Cp x /c R O E   F V M T V F DM FL U E N T   1 s t   O r d e r   U p w i n d     F i gu r e   10 .   P r e s s u r e   c o e ff i c i e n t   d i s t ri b ut i o n,   Cp,   a t   α=     a n d   M   =   0 . 8       Co n s i de r i n F i gu r e   10  a b ov e   c o n c l ude s   t ha t   t h e   p r e s e n c e   of   s h o c k   w a v e   a t   t h e   uppe r   s u r f a c e   w i t t h e i r   s h o c l o c a t i o n   i s   a s   gi v e n   i n   T a b l e   3.   W h i l e   t ha t ,   a t   l o w e r   s ur f a c e ,   t h e r e   i s   n o   s h o c w a ve .   T h e   s h o c l o c a t i o n   pl a y s   a n   i m po rt a nt   r o l e ,   s i n c e   i t   w i l l   de t e rm i n e   t he   o ve r a l l   p r e s s ur e   c oe f f i c i e n t   di s t r i b ut i o n   l o o k - a l i ke s .   Co n s i de r i ng  t h e   r e s ul t   o f   s h o c l oc a t i o n,   t h e   t w o   d e ve l o pe c o m put e r s   a r e   a b l e   t o   pr o duc e   t h e   r e s ul t   w h i c h   i s   c l o s e   t o   e a c h   o t h e r   b ut   s l i g h t l y   di ff e r e n t   w i t h   t h e   r e s ul t   p r o duc e by   A N S Y S - F L U E N T   s o f t w a r e F i gu r e   11  s h o w   t h e   c o m pa r i s o n   o f   M a c h   n u m b e r   n e a r l y   t h e   s a m e   r e s ul t s   a t   t h e   gi v e n   f l o w   c o n di t i o n   α  =   5 a n M   =   0 . 8 .   T h e s e   t w o   f i gur e s   i n di c a t e   t ha t   t w o   de v e l o pe c o de s   a r e   i go o a g r e e m e nt   w i t e a c h   o t h e r .         (a )     (b )     F i gu r e   11 .   T h e   M a c n u m b e r   c o n t o u n e a r   t h e   N A CA   0012 ,   a t   α=     a nd  M   =   0 . f r o m     (a T V D   F D M   a n d   (b R o e   F V M       O v e r a l l ,   c o m pa r i s o n   w i t h   A N S Y S - F L U E N T   r e s ul t   i n di c a t e t h e   r e s ul t s   a r e   qui t e   di f f e r i n e s pe c i a l l y   t o   pr e di c t   t h e   s h o c k.   T h e   r e a s o n   i s   due   t h e   n u m e r i c a l   s c h e m e   t ha t   a ppl i e i n   t h e   s o f t w a r e .   E a c h   s c h e m e   ha s   t h e i r   o w n   c a pa b i l i t i e s   t o   c a pt u r e   a n p r e d i c t   t h e   di s c o n t i nui t i e s .   H i gh  r e s o l ut i o n   s uc h   T V D   s c h e m e   gi v e s   go o r e s ul t s   b ut   n e e s pe n d   a   l o t   o f   t i m e .       4.   C O N C LU S I O N   T h i s   p a pe r   p r e s e nt s   t h e   c o m put a t i o na l   s t udy   of   f l ow   a r o un d   a   N A CA   001 by   us i n g   R o e   F V M   S c h e m e   a nd  D a v i s - Y e e   T V D   S c h e m e .   T hr o ug h   c o m pa r i s o n   r e s ul t s   t o   t h e   c a s e   of   f l ow   pa s t   t hr o ugh   a   Evaluation Warning : The document was created with Spire.PDF for Python.
In do n e s i a J   E l e c   E ng  &   Co m S c i     IS S N :   2502 - 4752       Com put a t i on al   s t ud y   o f   f l ow   ar ound  N A C A   0012   by   us i n R oe   F V s c h e m e   and   …  ( F at i m ah   Y us op )   1025   s y m m e t ri c a l   a i r f o i l   N A CA   0012 ,   i t   c a n   b e   c o n c l ude t ha t   t he   de ve l o pe c o m put e r   c o de   i s   a b l e   t o   c a p t u r e   t h e   pr e s e n c e   o f   s h o c w a ve   i n   t h e   f l ow   f i e l d.   It   s e e m s   t h a t   t h e   de v e l o pe c o m put e r   c o de   m uc b e t t e r   i c a pt u r i n g   di s c o n t i n u i t y   f l ow   ph e n o m e na   c o m pa r e   t o   t h e   A N S Y S - F L U E N T   s of t w a r e   by   pr o duc i n a   s ha r p r e s s u r e   gra di e nt   a t   t h e   s h o c po i n t .   H ow e v e r ,   c o m pa r i s o n   w i t h   e xpe ri m e n t   i s   n e c e s s a r y   t o   e n s ur e   t ha t   t h e   de v e l o pe d   c o m put e r   c o de   pr e di c t   c o r r e c t l y   t h e   s h o c l o c a t i o n .   S i n c e   s h o c po s i t i o n   i s   v e r y   i m po r t a nt   f o r   ob t a i n i ng  a   c o r r e c t   v a l ue   o f   t h e   o ve r a l l   l i f t ,   d ra a n pi t c h i ng  m o m e n t   c o e ff i c i e n t   f o r   t h e   b o d y   i m m e r s e i n   t h e   c o m pr e s s i b l e   f l ow   f i e l d.   T hi s   r e p r e s e nt s   t h e   s ugge s t e f o r   t h e   f ut ur e   w o r k.       A C K N O WL ED G E M EN TS     T h e   a ut h o r   w o ul l i ke   t o   U n i v e r s i t i   T u n   H us s e i n   O nn   M a l a y s i a   (U T H M )   fo r   t h e   c o n t i n uo us   s uppo r t   i n   t h i s   p r o j e c t       R EF ER EN C ES   [ 1]   E .   N .   T i no c o ,   " T he   I m pa c t   o f   H i g P e r f o r m a nc e   C o m put i ng   a nd  C o m put a t i o na l   F l u i D y na m i c s   o A i r c r a f t   D e v e l o pm e nt   "   i 20t A nni v e r s ar y   D i nne r   &   Sy m pos i um   C A SC   W as hi nt on,   D C :   C o al a t i on  f o r   A c ade m i c   Sc i e n t i f i c   C om put at i on ,   2009 .   [ 2]   F .   T .   J o hns o n ,   E .   N .   T i no c o ,   a nd  N .   J .   Y u ,   " T h i r t y   y e a r s   o f   de v e l o pm e nt   a nd  a ppl i c a t i o o f   C F D   a t   B o e i ng   C o m m e r c i a l   A i r pl a ne s ,   S e a t t l e , "   C o m put e r s   &   F l u i ds ,   v o l .   3 4,   pp .   1 115 - 1151 ,   200 5.   [ 3]   A .   J a m e s o n,   W .   S c hm i dt ,   a n E .   T u r k e l ,   " N um e r i c a l   s o l ut i o ns   o f   t he   E u l e r   e q ua t i o ns   by   f i ni t e   v o l um e   m e t ho ds   us i ng   R ung e - K ut t a   t i m e - s t e pp i ng   s c he m e s , "   A I A A   P a pe r ,   v o l .   12 59 ,   p.   19 81,   1 981 .   [ 4]   M .   D r e l a ,   " N e w t o s o l ut i o o f   c o upl e v i s c o us / i nv i s c i m ul t i   e l e m e nt   a i r f o i l   f l o w s , "   i 21 s t   F l ui D y na m i c s ,   P l a s m a   D y na m i c s   a nd  L a s e r s   C o nf e r e nc e :   A m e r i c a I ns t i t ut e   o f   A e r o na ut i c s   a nd   A s t r o na u t i c s ,   1990 .   [ 5]   H .   G hua ng hui ,   " N um e r i c a l   s i m u l a t i o ns   o f   t he   s t e a dy   E ul e r   e qu a t i o o uns t r uc t u r e g r i d . "   v o l .   P h. D :   H o n g   H o n B a pt i s t   U ni v e r s i t y ,   2009 ,   p.   1 48 .   [ 6]   M .   H .   C he n,   C .   C .   H s u ,   a nd   W .   S hy y ,   " A s s e s s m e nt   o f   T V D   s c he m e s   f o r   i nv i s c i d   a n t ur bul e n t   f l o w   c o m put a t i o n, "   I nt e r n at i on al   J o ur n al   f or   N um e r i c al   M e t hod s   i n   F l ui ds ,   v o l .   12 ,   pp .   161 - 17 7,   19 91 .   [ 7]   K .   R .   Q ur e s h i   a nd  C .   H .   L e e ,   " B e h a v i o r   of   T V D   L i m i t e r s   o t he   S o l ut i o o f   N o n - L i ne a r   H y pe r bo l i c   E qua t i o n, "   M o de r n   P hy s i c s   L e t t e r s   B ,   v o l .   1 9,   pp .   150 7 - 1510 ,   2005 .   [ 8]   R .   A bg r a l l R ,   A e x t e n s i o o f   R o e ' s   upw i nd  s c he m e   t o   a l g e b r a i c   e qui l i b r i um   r e a l   g a s   m o de l s .   C o m pu t e r s   &   F l u i ds ,   19 ( 2) ,   171 - 18 2,   19 91.   [ 9]   P .   G l a i s t e r ,   " A   w e a f o r m ul a t i o o f   R o e ' s   s c he m e   f o r   t w o - di m e ns i o na l ,   uns t e a dy ,   c o m pr e s s i bl e   f l o w s   a nd  s t e a dy ,   s upe r s o ni c   f l o w s , "   C o m put e r s   &   M a t h e m a t i c s   w i t A ppl i c a t i o ns ,   v o l .   3 0,   pp .   85 - 93 ,   1995 .   [ 10]   C .   H i r s h ,   N um e r i c a l   C o m put a t i o o f   I nt e r n a l   a n E xt e r n a l   F l o w s   ( 2 nd  E d . ) .   G r e a t   B r i t i a n:   E l s e v i e r ,   200 7.   [ 11]   K .   A .   H o f f m a &   S .   T .   C h i a ng ,   C o m put a t i o na l   F l ui d   D y na m i c s   ( 4 t e d .   V o l .   I I ) .   U S A :   w w w . E E S bo o ks . c o m .   2 000 .   [ 12]   J .   C .   T a nne hi l l ,   D .   A .   A nde r s o n,   a nd  R .   H .   P l e t c he r ,   C o m put a t i o na l   F l u i M e c ha n i c s   a n H e a t   T r a n s f e r ,   2n e d.   U ni t e d   S t a t e   o f   A m e r i c a :   T a y l o r   &   F r a nc i s ,   19 97.   [ 13]   J .   B l a z e k,   C o m put a t i o na l   F l ui d   D y na m i c s :   P r i nc i pl e s   a n A ppl i c a t i o ns :   E l s e v i e r ,   200 5   [ 14]   A .   H a r t e n,   " H i g r e s o l u t i o s c he m e s   f o r   hy pe r bo l i c   c o ns e r v a t i o l a w s , "   J o ur n al   o f   C om put at i on al   P hy s i c s ,   v o l .   4 9,   pp.   35 7 - 393,   1 983 .   [ 15]   T .   J .   C h ung ,   C o m put a t i o na l   f l u i d   dy na m i c s :   C a m br i dg e   uni v e r s i t y   pr e s s ,   2 010 .   [ 16]   Y .   L .   C ha ng ,   " D e v e l o pm e nt   o f   a   C F D   c o de   us i ng   T V D   s c h e m e s   a n a dv a nc e t u r bu l e nc e   m o de l s   f o r   i nc o m pr e s s i b l e   f l o w   s i m ul a t i o ns . "   P h . D :   M i c hi g a n   T e c hno l o g i c a l   U ni v e r s i t y ,   1997 ,   p.   10 7.   [ 17]   S .   F .   D a v i s ,   A   S i m pl i f i e d   T V D   F i n i t e   D i f f e r e nc e   S c he m e   v i a   A r t i f i c i a l   V i s c o s i t y .   SI A M   J o ur n al   on   Sc i e nt i f i c   and   St a t i s t i c a l   C om pu t i ng ,   8( 1 ) ,   1 - 18 ,   198 7.     [ 18]   H .   C .   Y e e   a nd  A .   H a r t e n,   " I m pl i c i t   T V D   s c he m e s   f o r   h y pe r bo l i c   c o ns e r v a t i o l a w s   i c ur v i l i ne a r   c o o r di na t e s , "   A I A A   J o ur na l ,   v o l .   2 5,   pp .   266 - 27 4,   20 15 / 07 / 29   19 85.   [ 19]   S .   C ha kr a v a r t hy   a nd  S .   O s h e r ,   " A   ne w   c l a s s   o f   hi g a c c ur a c y   T V D   s c he m e s   f o r   hy pe r bo l i c   c o ns e r v a t i o l a w s , "     i 23r A e r o s pa c e   S c i e nc e s   M e e t i ng :   A m e r i c a I ns t i t u t e   o f   A e r o na ut i c s   a n A s t r o na ut i c s ,   1 985 .   [ 20]   K .   A nde r s o n,   J .   L .   T ho m a s ,   a nd  B .   V a L e e r ,   " C o m pa r i s o o f   f i ni t e   v o l um e   f l ux   v e c t o r   s pl i t t i ng   f o r   t he   E ul e r   e qua t i o ns , "   A I A A   J o ur na l ,   v o l .   24 ,   p p.   14 53 - 1460 ,   198 6.   [ 21]   A .   J a m e s o n,   R e qu i r e m e n t s   a n t r e n ds   o f   c o m put a t i o na l   f l u i dy na m i c s   a s   a   t o o l   f o r   a i r c r a f t   de s i g n .   P r oc e e di ngs   of   t he   12t N A L   s y m po s i um   o a i r c r af t   c om pu t a t i ona l   ae r ody n am i c s .   T o ky o ,   J a pa n.   19 94 .   [ 22]   H .   S h i g e ki ,   A .   J u s t i n ,   &   H .   K l a us ,   D e v e l o pm e n t   o f   a   m o di f i e R ung e - K ut t a   s c he m e   w i t h   T V D   l i m i t e r s   f o r   t he   i de a l   t w o - di m e ns i o na l   M H D   e qu a t i o ns .   3 6t h   A I A A   A e r o s pa c e   S c i e nc e s   M e e t i ng   a nd  E x hi bi t .   A m e r i c a I ns t i t u t e   o f   A e r o na ut i c s   a n A s t r o na ut i c s .   1 998 .     [ 23]   P .   G l a i s t e r ,   F l ux   di f f e r e nc e   s p l i t t i ng   f o r   t he   E u l e r   e qu a t i o ns   i n   o ne   s pa t i a l   c o - o r di na t e   w i t a r e a   v a r i a t i o n.   I nt e r n at i on al   J o ur n al   f or   N um e r i c al   M e t hod s   i n   F l ui ds .   1988 .   8:   97 - 119.     [ 24]   S .   I .   S o hn,   A   ne w   T V D - M U S C L   s c he m e   f o r   h y pe r bo l i c   c o ns e r v a t i o l a w s .   C o m put e r s   &   M a t he m a t i c s   w i t A ppl i c a t i o ns ,   50 ( 2) ,   231 - 24 8,   20 05 .     [ 25]   L .   M o t t ur a ,   L .   V i g e v a no ,   &   M .   Z a c c a nt i ,   A E v a l ua t i o o f   R o e ' s   S c he m e   G e ne r a l i z a t i o ns   f o r   E qu i l i b r i um   R e a l   G a s   F l o w s .   J ou r na l   o f   C om pu t a t i ona l   P hy s i c s ,   1 38( 2) ,   354 - 399 M   1991 .   Evaluation Warning : The document was created with Spire.PDF for Python.