2 Ducts / 1 Junction Test TSA-6 jws GF2 $Source: /home/cvs/cvsroot/input/blower.ingf,v $ $Date: 1996/12/18 10:51:53 $ $Locker: $ $Revision: 1.1 $ $State: Exp $ $Author: jspore $ $Log: blower.ingf,v $ Revision 1.1 1996/12/18 10:51:53 jspore Initial revision Revision 1.6 1996/07/11 15:06:37 jspore BLOWDEF had one too many inputs. Revision 1.5 1995/10/30 16:18:37 jspore Changed autot from 0.0 to 1.0. Revision 1.4 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.3 1995/05/08 21:29:02 jspore Added header information. ------------------------------------------------------------------------ N O T E S ----------- * BLOWER.AR test case * 12 AUG 93 -- date of latest change * 30 AUG 94 -- input modified for GF2 * 2 ducts, 1 jct * file: INGF stored in an ar library file on * /089534/gasflow/gasnet/test_problems/blower.ar * duct 1: 4 x 100 = 400 cm, duct 2: 3 x 100 = 300 cm, jct 1: 100 cm/side * Number of real cells = 4 + 3 + 1 = 8 * Cartesian coordinates, uniform grid * Working fluid: air * Blower at cell 4 duct 1 * WEST B.C.: specified pressure * EAST B.C.: specified pressure node 1 node 2 node 3 PBC PBC +...+---+---+---+---+---+---+---+---+...+ . 1| 2| 3| 4| 5| 12| 8| 9| 10| 11. . B | | | | |J1 | | | | B . +...+---+---+---+---+---+---+---+---+...+ D1 D2 ------------------------------------------------------------------------ N E T W O R K P A R A M E T E R S ------------------------------------------------------------------------ $innet ductdef(1,1) = 0.0, 0.0, 100.0, 100.0, 100.0, 100.0, 1, 2, 4, 1.0, 0, 0, ductdef(1,2) = 0.0, 0.0, 100.0, 100.0, 100.0, 100.0, 2, 3, 3, 1.0, 0, 0, defjct(1,1) = 100.0, 100.0, 100.0, 2, cpnt(1,1) = 0.0, 0.0, 0.0, cpnt(1,2) = 450.0, 0.0, 0.0, cpnt(1,3) = 800.0, 0.0, 0.0, blowdef(1,1) = 3 , 4 , 1 , 1 , 1, 1, 40.0, 0.001, 100000.0, 50.0, blwtb(1,1,1) = 0.0, 3000.0, 1.0e05, 200.0, 2.0e05, 100.0, 3.0e05, 50.0, 4.0e05, 25.0, 4.5e05, 10.0, 5.0e05, 0.0, 5.5e05, -10.0, spdtb(1,1,1) = 0.0, 0.0, 10.0, 10.0, 20.0, 30.0, 30.0, 40.0, 50.0, 40.0, 100.0, 40.0, iwshear = 0, netopt = 2, $end ------------------------------------------------------------------------ B A S I C I N P U T ------------------------------------------------------------------------ $xput muoption = 1, nu = 0.150 , gz = -000.0, autot = 1.0, delt0 = 0.1000, deltmin = 1.000e-04, deltmax = 1.000e-01, nrsdump = 0, epsi0 = 1.000e-08, epsimax = 1.000e-08, epsimin = 1.000e-08, itdowndt = 024, itupdt = 024, itmax = 025, maxcyc = 2000, twfin = 100.000, lpr = 1, iobpl = 1, ittyfreq = 20, pltdt = 020.000, prtdt = 10.0 , tddt = 5000.00, velmx = 1.5, mat = 'air' , gasdef(1,1) = 0, 6, 1, -1, -1, -1, 1 , 1.0000000000e6, 300.00, 1, 0., 0., 'air' , 1.00000, gasdef(1,2) = 0, 5, 2, -2, -2, -2, 1 , 1.0000000000e6, 300.00, 1, 0., 0., 'air' , 1.00000, gasdef(1,3) = 0, 1, 1, -1, -1, -1, 1 , 0.9999000000e6, 300.00, 1, 0., 9999999., 'air' , 1.00000, gasdef(1,4) = 4, 5, 2, -2, -2, -2, 1 , 1.0000000000e6, 300.00, 1, 0., 9999999., 'air' , 1.00000, vvalue = 123.45, pbc(1,1) = 1, 0, 1, 0, 0, 0, 0, 0.0, 99999., pbc(1,2) = 2, 0, 2, 0, 0, 0, 0, 0.0, 99999., $end ------------------------------------------------------------------------ M E S H ------------------------------------------------------------------------ $meshgn $end ------------------------------------------------------------------------ G R A P H I C S ------------------------------------------------------------------------ $grafic pnt(1,1) = 1, 1, -1, -1, pnt(1,2) = 5, 1, -1, -1, thdt = 0.0100, thp(1,1) = 1, 1,-1,-1, 'un' , 0, thp(1,2) = 5, 1,-1,-1, 'un' , 0, thp(1,3) = 1, 2,-2,-2, 'un' , 0, thp(1,4) = 4, 2,-2,-2, 'un' , 0, $end ------------------------------------------------------------------------ P A R T I C L E S ------------------------------------------------------------------------ $parts $end ------------------------------------------------------------------------ H E A T T R A N S F E R & C O N D E N S A T I O N ------------------------------------------------------------------------ $rheat $end ------------------------------------------------------------------------ S P E C I A L ------------------------------------------------------------------------ $special $end c----------------------------------------------------------------------- c I N F O R M A T I O N c----------------------------------------------------------------------- c***** ENERGY c 1 erg = 1.e-7 J = 1 gm-cm**2/sec**2 c 1 J/gm-K = 1.e7 erg/gm-K = 1.e7 cm**2/sec**2-K c 1 J/gm = 1.e7 erg/gm = 1.e7 cm**2/sec**2 c***** PRESSURE c 1 dyne = 1 gm-cm/sec**2 c 1 dyne/cm**2 = 1 gm/(cm-sec**2) = 0.1 Pa c 1 dyne/cm**2 = 1.e-7 J/cm**3 c 1 bar = 1.e5 Pa = 1.e6 dynes/cm**2 = 0.98692 atm c 1 atm = 1.013253 bars c***** TEMPERATURE c TK = (TF + 459.67)/1.8 c TF = 1.8*TK - 459.67 c ugasc = 8.31441e7 gm-cm**2/(s**2-gm-mole-K) universal gas constant. c c Gas proprties at T = 300 K. c c ---------------------------------------------------------------------- c Component Specie M R Cv gamma Hf0 c number (J/gm-K) (J/gm-K) Cp/Cv (J/gm) c ---------------------------------------------------------------------- c 1 C 12.011 0.0000 0.709932 0.0000 0.0 c 2 CO 28.011 0.2968 0.74318 1.399 -3946.353 c 3 CO2 44.011 0.1889 0.6544 1.288 -8941.445 c 4 H2 2.016 4.12418 10.0849 1.409 0.0 c 5 H2O 18.016 0.46152 1.4108 1.327 -13422.87 c 6 N2 28.013 0.2968 0.7448 1.4 0.0 c 7 N2O 44.013 0.18892 0.6912 1.273 1864.182 c 8 O2 32.000 0.25983 0.6618 1.393 0.0 c 9 air 28.97 0.287 0.716802 1.4 c 10 AR 39.948 0.20813 0.310685 1.6699 0.0 c 11 HE 4.003 2.07705 3.114010 1.667 0.0 c 12 NH3 17.0305 0.48821 1.534400 1.3182 -2695.075 c 13 CH4 16.0430 0.51826 1.704795 1.304 -4667.0 c 14 HO 17.008 0.48885 1.274243 1.384 2320.29 c 15 H 1.008 8.24842 12.37272 1.667 216256.3 c 16 HO2 33.008 0.25189 0.805001 1.3129 633.7857 c 17 NO 30.0065 0.277087 0.717513 1.3862 3009.039 c 18 O 16.000 0.519651 0.849825 1.6115 15574.68 c 19 NH 15.0145 0.553759 1.387414 1.3991 22571.78 c 20 HNO 31.0145 0.268081 0.848794 1.3158 3210.73 For the Nth defined fluid volume: gasdef( 1, N) = i1 Starting mesh index in the I-direction. gasdef( 2, N) = i2 Ending mesh index in the I-direction. gasdef( 3, N) = j1 Starting mesh index in the J-direction. gasdef( 4, N) = j2 Ending mesh index in the J-direction. gasdef( 5, N) = k1 Starting mesh index in the K-direction. gasdef( 6, N) = k2 Ending mesh index in the K-direction. gasdef( 7, N) = iblk block number (1 for now) gasdef( 8, N) = P Pressure (dynes/cm**2). gasdef( 9, N) = T Temperature (K). gasdef(10, N) = iopt EOS option flag (1 for mass fracs, 2 for mole fracs) gasdef(11, N) = tstart Start time for this specification (sec). gasdef(12, N) = tend End time for this specification (sec). gasdef(13, N) = mn1 First gas specie (hollerith) gasdef(14, N) = xfrac mass (or volume) fraction of material MN1. gasdef(15, N) = mn2 Second gas specie (hollerith) gasdef(16, N) = xfrac mass (or volume) fraction of material MN2. etc. NOTE: One pair of variables (MN, XFRAC) is defined for each gas specie. If iopt = 1 , XFRAC is interpreted as a mass fraction; if iopt = 2, XFRAC is interpreted as a volume fraction. The sum of the XFRACs must be 1.0000 or an error will occur. At least one material must be defined. veldef(1,1) = 01, 01, 1,14, 09, 17, 1, 1, 0.0, 500.00,