1D/3D Problem #1B * 1B_1D_a * PBC=Phi * 1D/3D * Non-Uniform P TSA6 JWS GASFLOW V2.0 c $Source: /home/cvs/cvsroot/input/1B_1D_a.ingf,v $ c $Date: 1998/09/29 16:10:58 $ c $Locker: $ c $Revision: 1.2 $ c $State: Exp $ c $Author: necker $ c c $Log: 1B_1D_a.ingf,v $ c Revision 1.2 1998/09/29 16:10:58 necker c Changes from Jay Spore/LANL: c Ductdef overspecified duct cells. c Changed itupdt to 24. c c Revision 1.1 1996/12/18 09:47:15 jspore c Initial revision c c Revision 1.5 1995/10/18 14:43:46 jspore c Ran all .ingf files through edgf to make for sure all namelist c ends started in column 2. Fred also removes trailing blanks. c c Revision 1.4 1995/05/15 16:27:32 jspore c Added extra null mesh input. This allows SETLSIZE logic for reading c the MESHGN namelist to work correctly on INDY machines. c c Revision 1.3 1995/05/01 16:21:29 jspore c Modified input to be consistent with the new CNVTINP routine. c This means no hollerith input and no embedded blanks in character c input. c Revision 1.2 95/04/09 23:56:55 jspore Removed the old particle input. With Billy's new aerosol model the old particle input will result in an input error. Revision 1.1 95/03/15 17:33:00 jspore Initial revision NOTES: * 7 OCT 94 * 1B_1D_a on ODIN in directory /home/jspore/gasflow/gf2/test_problems * 1 3D block, 1 duct * air * BLOCK #1 300 x 300 x 100 cm 3 x 3 x 1 = 9 real cells m222 = 32 5 x 5 x 3 = 75 total cells mrrr = 44 * DUCT #1 400 x (50 x 100) cm 4 cells W end: connected to 3D block E end: PBC J=4 +-------+-------+-------+ | 42 | 43 | 44 | | | | | | | | | J=3 +-------+-------+-------+ | 37 | 38 | 39 | | | | | | | | | J=2 +-------+-------+-------+-------+-------+-------+-------+ . . . . | 32 | 33 | 34 | 77 | 78 | 79 | 80 | 81 . | | | 1 | 2 | 3 | 4 | 5 | 6 . | | | | | | | | P=Phi . J=1 +-------+-------+-------+-------+-------+-------+-------+ . . . . 1 2 3 4 1 2 3 4 5 NODE 1 NODE 2 3D BLOCK NETWORK DUCT ----------------------------------------------------------------------- N E T W O R K I N P U T ----------------------------------------------------------------------- $innet cpnt = 300.0, 50.0, 50.0, 700.0, 50.0, 50.0, ductdef = 0.0, 0.0, 100.0, 100.0, 100.0, 100.0, 1, 2, 4, 1.0, 0.0, 0.0, netopt = 2, kopt1d3d = 1, nwcx = 4, 4, 1, 2, 1, 2, 1, 1, iwshear = 1, $end ----------------------------------------------------------------------- B A S I C I N P U T ----------------------------------------------------------------------- $xput pbc(1,1) = 2, 0, 1, 0, 0 , 0 , 0 , 0.00, 1.e99, pbc(1,2) = 1, 4, 1, 4, 2 , 2 , 1 , 0.00, 1.e99, vvalue = 1000.0,-500.0, 0.0, ; vbc(1,1) = 2, 0, 1, 0, 0 , 0 , 1 , 3, 0.00, 1.e99, autot = 1.0, cyl = 0.0, omg = 34.0 , delt0 = 00.01000 , deltmin = 1.000e-08, deltmax = 1.000e-00, epsi0 = 1.000e-08, epsimax = 1.000e-08, epsimin = 1.000e-08, gz = -000.0, iobpl = 0, itdowndt = 25, itupdt = 24, itmax = 100, lpr = 1, maxcyc = 05000, ittyfreq = 50, nu = 0.15 , pltdt = 025.00, prtdt = 025.000, twfin = 025.00, tddt = 100.0000, velmx = 1.5, ibw = 1, ibe = 1, ibs = 1, ibn = 1, ibb = 1, ibt = 5, pamb0 = 1.0e+06, mat = 'air', ; 3D BLOCK ------------------------------------------------- gasdef(1,1) = 1 , 'im1', 1 , 'jm1' , 1 , 'km1', 1 , ; I.C. 1.00100000e6, 300.0, 1, 0., 0., 'air', 1.00000, gasdef(1,2) = 0 , 1, 1 , 4 , 01 , 02, 1 , ; west B.C. 1.00100000e6, 300.0, 1, 0.0, 99999.0, 'air', 1.00000, gasdef(1,3) = 'im1','imax', 1, 4, 1 , 02, 1 , ; east B.C. 1.00100000e6, 300.0, 1, 0.0, 0.0, 'air', 1.00000, gasdef(1,4) = 1, 4, 0, 1, 1 , 02, 1 , ; south B.C 1.00100000e6, 300.0, 1, 0.0, 99999.0, 'air', 1.00000, gasdef(1,5) = 1, 4, 4, 5, 1 , 02, 1 , ; north B.C 1.00100000e6, 300.0, 1, 0.0, 99999.0, 'air', 1.00000, gasdef(1,6) = 1, 4, 1, 4, 0 , 1, 1 , ; bttm B.C. 1.00100000e6, 300.0, 1, 0.0, 99999.0, 'air', 1.00000, gasdef(1,7) = 1, 4, 1, 4, 2 , 3, 1 , ; top B.C. 1.00100000e6, 300.0, 1, 0.0, 99999.0, 'air', 1.00000, ; 1D NETWORK ---------------------------------------------- gasdef(1, 8) = 0 , 6, 1 , 0 , 0 , 0, 0 , ; D#1 I.C. 1.00000000e6, 300.0, 1, 0., 0., 'air', 1.00000, gasdef(1, 9) = 5 , 6, 1 , 0 , 0 , 0, 0 , ; D#1 E BC 1.00000000e6, 300.0, 1, 0.0, 99999.0, 'air', 1.00000, $end ------------------------------------------------------------------------ M E S H ------------------------------------------------------------------------ $meshgn iblock = 1, xgrid = 0.0, 100.0, 200.0, 300.0, ygrid = 0.0, 100.0, 200.0, 300.0, zgrid = 0.0, 100.0, $end $meshgn $end ------------------------------------------------------------------------ G R A P H I C S ------------------------------------------------------------------------ $grafic pnt(1, 1) = 1, 1, 2, 1, 'im1', 'jm1', 2, 1, 1, 2, 2, 1, 'im1', 2, 2, 1, 1, 1, 0, 0, ; duct #1 5, 1, 0, 0, ; duct #1 p1d = 3, 4, 'un' , 0, 3, 4, 'pn' , 0, 3, 4, 'diffp', 0, 5, 6, 'un' , 0, ; duct #1 5, 6, 'pn' , 0, ; duct #1 5, 6, 'diffp', 0, ; duct #1 v2d = 1, 2, 1, thdt = 1.000, thp(1,1) = 4, 2, 2, 1, 'pn', 0, 4, 2, 2, 1, 'un', 0, 4, 2, 2, 1, 'rn', 0, 4, 2, 2, 1, 'tk', 0, 2, 1, 0, 0, 'pn', 0, ; duct #1 5, 1, 0, 0, 'pn', 0, ; duct #1 2, 1, 0, 0, 'un', 0, ; duct #1 5, 1, 0, 0, 'un', 0, ; duct #1 $end ------------------------------------------------------------------------ P A R T I C L E S ------------------------------------------------------------------------ $parts $end $rheat $end $special $end $specialp $end