1D/3D Problem BURN3D1D test problem
TSA6 jws
GASFLOW V2.0
NOTES:
 
$Source: /home/cvs/cvsroot/input/burn3d1d.ingf,v $
$Date: 1996/12/18 15:19:54 $
$Locker:  $
$Revision: 1.1 $
$State: Exp $
$Author: jspore $
 
$Log: burn3d1d.ingf,v $
Revision 1.1  1996/12/18 15:19:54  jspore
Initial revision

Revision 1.6  1995/10/18 14:43:46  jspore
Ran all .ingf files through edgf to make for sure all namelist
ends started in column 2.  Fred also removes trailing blanks.

Revision 1.5  1995/05/15 16:27:32  jspore
Added extra null mesh input.  This allows SETLSIZE logic for reading
the MESHGN namelist to work correctly on INDY machines.
 
Revision 1.4  1995/05/08 21:29:02  jspore
Added header information.
 
 
 * 12/21/94
 * burn3d1d.ar  on  ODIN in directory /home/jspore/gasflow/gf2/test_problems
 * 1 3D block and 1 1D ducts.
 * 101-SY burp gas composition
 * No flow in or out to determine peak pressure and temperature.
 * BLOCK #1
      400 x 100 x 100 cm
      4 x 1 x 1 =  4 real  cells    m222 = 26   mrrr = 29
 * Duct # 1
      4 100 cm long cells with flow area of 100 X 100 cm
      4 real cells
 
 
   J=2 +-------+-------+-------+-------+
       |    26 |    27 |    28 |    29 |
       |       |       |       |       |
       |       |       |       |       |
   J=1 +-------+-------+-------+-------+
       1       2       3       4       5
                   3D BLOCK 1
 
   node 1              node 2
 
   connects             PBC
    3D block
  +...+---+---+---+---+...+
  .  1|  2|  3|  4|  5|  6.
  .   |   |   |   |   |  B.
  +...+---+---+---+---+...+
            D1
 
 
 -----------------------------------------------------------------------
     N E T W O R K    I N P U T
 -----------------------------------------------------------------------
 $innet
 
       netopt    = 2,
       ductdef(1,1) =   0.0,   0.0, 100.0, 100.0, 100.0, 100.0, 1, 2, 4,
                        1.0,     0,     0,
       cpnt(1,1)    =   0.0,   0.0,   0.0,
       cpnt(1,2)    = 400.0,   0.0,   0.0,
       nwcx =  5 , 5 , 1 , 2 , 1, 2 , 1 , 1 ,
       iwshear = 0,
 
 $end
 
 -----------------------------------------------------------------------
          B A S I C    I N P U T
 -----------------------------------------------------------------------
 $xput
 
      autot    = 1.0,  cyl = 0.0,   omg =   34.0 ,
      delt0    = 00.00100  ,
      deltmin  = 1.000e-08,
      deltmax  = 1.000e-03,
      epsi0    = 1.000e-08,
      epsimax  = 1.000e-08,
      epsimin  = 1.000e-08,
      gz       = -000.0,
      iobpl    = 0,
      itdowndt =  25,
      itupdt   = 100,
      itmax    = 100,
      lpr      = 1,
      maxcyc   = 04000,
      ittyfreq =  20,
      nu       = 0.15  ,
      pltdt    = 100.00,
      prtdt    =   3.000,
      twfin    =  003.00,
      tddt     = 100.0000,
      velmx    = 1.5,
      ibw      = 1,
      ibe      = 1,
      ibs      = 1,
      ibn      = 1,
      ibb      = 1,
      ibt      = 1,
 
      burndef(1,1) = 1, 2, 1, 2, 1, 2, 1, 1000.0, 1.0, 9999.0,
      iburn = 2,
 
      mat =  'n2' , 'o2' , 'h2',
             'n2o', 'h2o', 'co',
             'co2', 'nh3', 'ch4',
 
 
 ; 3D BLOCK  # 1----------------------------------------------
      gasdef(1,1) = 1 ,  'im1',  1 ,  'jm1' ,  1 , 'km1',   1 ,   ; I.C.
                    1.00000000e6, 300.0, 1, 0., 0.,
            'n2', 0.70895, 'h2', 0.00687, 'n2o', 0.01286,
            'h2o', 0.00474, 'co', 0.00154, 'nh3', 0.02780,
            'ch4', 0.00094, 'o2', 0.2363,
 ;          'n2', 0.72952, 'h2', 0.00344, 'n2o', 0.00643,
 ;          'h2o', 0.00237, 'co', 0.00077, 'nh3', 0.01390,
 ;          'ch4', 0.00047, 'o2', 0.2431,
 
     gasdef(1,2) =  0,  6,  1, -1,  -1, -1,   0 ,
                    1.0000000000e6, 300.00, 1, 0., 0.,
            'n2', 0.70895, 'h2', 0.00687, 'n2o', 0.01286,
            'h2o', 0.00474, 'co', 0.00154, 'nh3', 0.02780,
            'ch4', 0.00094, 'o2', 0.2363,
 ;          'n2', 0.72952, 'h2', 0.00344, 'n2o', 0.00643,
 ;          'h2o', 0.00237, 'co', 0.00077, 'nh3', 0.01390,
 ;          'ch4', 0.00047, 'o2', 0.2431,
 
      gasdef(1,3) =  5,  6,  1, -1,  -1, -1,   0 ,
                    1.0000000000e6, 300.00, 1, 0., 99999.,
                     'n2' , 1.00 ,
      vvalue = 0.0,
 
      vbc(1,1) =  2,  0, 1, 0,  0,  0, 0, 1, 0.0, 99999.,
 
 $end
------------------------------------------------------------------------
          M E S H
------------------------------------------------------------------------
 $meshgn
     iblock = 1,
     xgrid =    0.0, 100.0, 200.0, 300.0, 400.0,
     ygrid =    0.0, 100.0,
     zgrid =    0.0, 100.0,
 $end
 $meshgn
 $end
 
------------------------------------------------------------------------
          G R A P H I C S
------------------------------------------------------------------------
 $grafic
 
  thdt = 0.01,
  thp(1,1)  =  2,  2,  2, 1,  'tk'  ,  0,
               2,  2,  2, 1,  'pn'  ,  0,
               5,  2,  2, 1,  'tk'  ,  0,
               5,  2,  2, 1,  'pn'  ,  0,
               2,  1,  0, 0,  'tk'  ,  0,
               2,  1,  0, 0,  'pn'  ,  0,
               5,  1,  0, 0,  'tk'  ,  0,
               5,  1,  0, 0,  'pn'  ,  0,
 
 $end
------------------------------------------------------------------------
          P A R T I C L E S
------------------------------------------------------------------------
 $parts
 $end
 $rheat
 $end
 $special
 $end
 $specialp
 $end