GASFLOW - burn3d1d

burn3d1d

Home Page >> Verification >> burn3d1d

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