Page 391 - Computational Fluid Dynamics for Engineers
P. 391

384                                                            Appendix  B


                                    t
         10.14  Model  Problem  for he  Implicit  Method:
         Quasi-ID  Nozzle   (nozzle_imp.c)
         Generate  executable  nozzle.imp  on  unix  by  using
         cc  -o  nozzle_imp  nozzle_imp.c  b l o c k t r i . c  -lm
         Test  Cases  with  PC  executable  nozzleImpPC.exe  and  NMAX  =  100
         1. CFL  =  0.5,  e e  =  1, Output:  CFL0p5results.txt  and  CFL0p5residuals.txt
         2.  CFL  =  1.0,  e e  =  1, Output:  CFLlresults.txt  and  CFLlresiduals.txt
         3.  CFL  =  3.0,  e e  =  1, Output:  CFL3results.txt  and  CFL3residuals.txt
         4.  CFL  =  1.0,  s e  =  10, Output:  CFLlEelOresults.txt  and
                                      CFLlEelOresiduals.txt



         Chapter   11
         Incompressible    Navier-Stokes     Equations


         11.2.  INS2D  Method   for  Incompressible  Navier—Stokes  Equations

         ins2d. f  contains  main  program  which  reads  input  data  and  controls  the  itera-
         tions  and  computations.
        Residual.f  contains  subroutines  flux.de, phLih,  residualx,  flux-df,  phLjh  and
         residualy.
         Subroutines  fluxAE  and  fluxAF  are  for  numerical  flux  computations,  see  Eq.
         (11.2.15).
         Subroutines  phLih  and  phLjh  are  defined  according  to  Eq.  (11.2.17).
         Subroutines  residualx  and  residualy  calculates  righthandside  according  to  Eq.
         (11.2.22).
         solver.  f  contains  subroutines  solver,  invmatrixS,  tribandsyssolver  and  thomas
        which  use  the  block  elimination  method  to  solve  block-tridiagonal  system.
         Test  cases  for  sudden  expansion  laminar  duct  flow with  various  values  of  /?,  grid
         size  and  Reynolds  number.
         Test  Cases  for  Sudden  Expansion  Laminar  Duct  Flow:
         1.  (3 — 1,  grid  size  nx  =  200  and  ny  =  100,
           Input:  blg200X100Inp.txt
            Output:  blg200X1000ut_g.txt;  blg200X1000ut_Mass.txt;
           blg200X1000ut_p.txt;  blg200X1000ut_u.txt;   blg200X1000ut_v.txt
         2.  /3 =  50,  grid  size  nx  =  200  and  ny  =  100,
            Input:  b50g200X100Inp.txt
            Output:  b50g200X1000ut_g.txt;  b50g200X1000ut_Mass.txt;
                                                     t
           b50g200X1000ut_p.txt;   b50g200X1000ut_u. x t ;  b50g200X1000ut_v.txt
         3.  (5 =  100, grid  size  nx  =  200  and  ny  =  100,
           Input:  bl00g200X100Inp.txt
   386   387   388   389   390   391   392   393   394   395   396