GASFLOW - aerosol

aerosol

Home Page >> Verification >> aerosol

Simple Aerosol Test 5x5x5
TSA10 JWS 11/97
GASFLOW

$Id: aerosol.ingf,v 1.1 1998/01/14 13:13:36 necker Exp $

$Log: aerosol.ingf,v $
Revision 1.1 1998/01/14 13:13:36 necker
Simple test problem for the aerosol model.

-----------------------------------------------------------------------
$innet
$end
$xput

      ihystat = 1,
      idiffmom = 1,
      tmodel   = 'alg',
      idiffme = 0,
      prandtl = 0.70,
      schmidt = 1.0,
      autot    = 1.0,
      cyl      = 0.0,
      delt0    = 0.1000 ,
      deltmin = 1.000e-05,
      deltmax = 1.000e-01,
      epsi0    = 1.000e-08,
      epsimax = 1.000e-08,
      epsimin = 1.000e-08,
      gz       = -980.0,
      iobpl    = 0,
      itdowndt = 90,
      itupdt   = 85,
      itmax    = 100,
      lpr      = 2,
      maxcyc   = 500000,
      ittyfreq =     20,
      nu       =   0.15,
      pltdt    = 20.0,
      prtdt    = 100.0,
      twfin    = 40.0,
      tddt     = 40.0,
      velmx    = 1.5,
      ibw      = 5,
      ibe      = 3,
      ibs      = 2,
      ibn      = 2,
      ibb      = 2,
      ibt      = 2,

; vbc(1,1) =   1, 2, 1, 2, 'km1', 'km1', 1, 1, 0.00, 1.e99,
; pbc(1,1) =   1, 1, 1, 2,      1,      1, 1,     0.00, 1.e99,
; vvalue = 0.0,

     mat =   'air',

      gasdef(1,1) =     0, 'imax',    0, 'jmax',     0, 'kmax',   1,   ;I.C.
                    1.013500e6, 295.00, 1, 0., 0.,
                    'air', 1.00000,
      gasdef(1,2) =     0, 'imax',    0, 'jmax',     0,       1,   1,   ;Bottom B.C.'s
                    1.013500e6, 295.00, 1, 0., 999999.0,
                     'air', 1.0,
      gasdef(1,3) =     0, 'imax',    0, 'jmax', 'km1', 'kmax',   1,   ;Top B.C.'s
                    1.013500e6, 295.00, 1, 0., 999999.0,
                     'air', 1.0,
      gasdef(1,4) =     0,      1,    0, 'jmax',     0, 'kmax',   1,   ;West B.C.'s
                    1.013500e6, 295.00, 1, 0., 999999.0,
                     'air', 1.0,
      gasdef(1,5) = 'im1', 'imax',    0, 'jmax',     0, 'kmax',   1,   ;East B.C.'s
                    1.013500e6, 295.00, 1, 0., 999999.0,
                     'air', 1.0,
      gasdef(1,6) =     0, 'imax',    0,      1,     0, 'kmax',   1,   ;Sout B.C.'s
                    1.013500e6, 295.00, 1, 0., 999999.0,
                     'air', 1.0,
      gasdef(1,7) =     0, 'imax', 'jm1', 'jmax',     0, 'kmax',   1,   ;North B.C.'s
                    1.013500e6, 295.00, 1, 0., 999999.0,
                     'air', 1.0,

$end
------------------------------------------------------------------------
          M E S H
------------------------------------------------------------------------
$meshgn
     iblock = 1,
     nkx = 1,
     xl(1) =   0.0, xc(1) = 0.0, nxl(1) = 0, nxr(1)= 5, dxmn(1) = 9999.,
     xl(2) = 250.0,

     nky = 1,
     yl(1) =   0.0, yc(1) = 0.0, nyl(1) = 0, nyr(1)= 5, dymn(1) = 9999.,
     yl(2) = 250.0,

     nkz = 1,
     zl(1) =   0.0, zc(1) = 0.0, nzl(1) = 0, nzr(1)= 5, dzmn(1) = 9999.,
     zl(2) = 250.0,
$end
$meshgn

$end
------------------------------------------------------------------------
          G R A P H I C S
------------------------------------------------------------------------
$grafic

     pnt(1, 1) = 1, 1, 1,   1,
     pnt(1, 2) = 6, 6, 1,   1,
     pnt(1, 3) = 1, 1, 2,   1,
     pnt(1, 4) = 6, 6, 2,   1,
     pnt(1, 5) = 1, 1, 3,   1,
     pnt(1, 6) = 6, 6, 3,   1,
     pnt(1, 7) = 1, 1, 4,   1,
     pnt(1, 8) = 6, 6, 4,   1,
     pnt(1, 9) = 1, 1, 5,   1,
     pnt(1,10) = 6, 6, 5,   1,
     pnt(1,11) = 6, 6, 1,   1,
     pnt(1,12) = 6, 6, 7,   1,
     pnt(1,13) = 1, 1, 1,   1,
     pnt(1,14) = 6, 6, 6,   1,
     pnt(1,15) = 1, 1, 1,   1,
     pnt(1,16) = 1, 6, 6,   1,
     pnt(1,17) = 2, 1, 1,   1,
     pnt(1,18) = 2, 6, 6,   1,
     pnt(1,19) = 3, 1, 1,   1,
     pnt(1,20) = 3, 6, 6,   1,
     pnt(1,21) = 4, 1, 1,   1,
     pnt(1,22) = 4, 6, 6,   1,
     pnt(1,23) = 5, 1, 1,   1,
     pnt(1,24) = 5, 6, 6,   1,
     pnt(1,25) = 6, 1, 1,   1,
     pnt(1,26) = 6, 6, 6,   1,
     pnt(1,27) = 1, 1, 6,   1,
     pnt(1,28) = 6, 6, 6,   1,
     pnt(1,29) = 3, 3, 1,   1,
     pnt(1,30) = 3, 3, 6,   1,
     pnt(1,31) = 1, 3, 1,   1,
     pnt(1,32) = 1, 3, 6,   1,

;   c2d       = 17, 18, 'pn', 0,
;               19, 20, 'pn', 0,
;               21, 22, 'pn', 0,
;               23, 24, 'pn', 0,
;               25, 26, 'pn', 0,

     v3d       = 13, 14, 1, 1,

     v2d       = 3, 4, 1,   ; k=2 2D velcotiy
                  5, 6, 1,   ; k=3 2D velcotiy
                  7, 8, 1,   ; k=4 2D velcotiy
                  9, 10, 1,   ; k=5 2D velcotiy
                 27, 28, 1,   ; k=6 2D velcotiy

     p1d       = 11, 12, 'pn', 0,
                 11, 12, 'wn', 0,
                 29, 30, 'diffp', 0,
                 31, 32, 'diffp', 0,

     thdt = 0.10,

      viewcrds(1,1) = 135.0, 134.0, 120.0, 135.0, -500.0, 300.0,
      viewcrds(1,2) = 135.0, 134.0, 120.0, 135.0, 134.0, 500.0,

     thp(1, 1) = 1, 1, 1, 1, 'wn', 0,
     thp(1, 2) = 6, 6, 6, 1, 'wn', 0,
     thp(1, 3) = 2, 2, 2, 1, 'wn', 0,
     thp(1, 4) = 2, 2, 2, 1, 'un', 0,
     thp(1, 5) = 2, 2, 2, 1, 'vn', 0,
     thp(1, 6) = 2, 2, 2, 1, 'pn', 0,

$end
------------------------------------------------------------------------
          P A R T I C L E S
------------------------------------------------------------------------
$parts
$end
$rheat
$end
$special
$end
$specialp
$end