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