burn1d
Home Page >>
Verification >>
burn1d
2 Ducts / 1 Junction Test
TSA-6 JWS
GF2
$Source: /home/cvs/cvsroot/input/burn1d.ingf,v $
$Date: 1996/12/18 10:51:53 $
$Locker: $
$Revision: 1.1 $
$State: Exp $
$Author: jspore $
$Log: burn1d.ingf,v $
Revision 1.1 1996/12/18 10:51:53 jspore
Initial revision
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
-----------
* BURN1D.AR test case
* 20 DEC 94
* 2 ducts, 1 jct
* file: INGF stored in a library on
* /089534/gasflow/gf2/test_problems/burn1d.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: flammable mixture
* A flow loss of one at cell 3 duct 1.
* WEST B.C.: specified zero velocity
* EAST B.C.: specified zero velocity
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
netopt = 2,
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,
iwshear = 0,
$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.00300,
deltmin = 1.000e-08,
deltmax = 1.000e-03,
nrsdump = 0,
epsi0 = 1.000e-08,
epsimax = 1.000e-08,
epsimin = 1.000e-08,
itdowndt = 024,
itupdt = 024,
itmax = 100,
maxcyc = 4000,
twfin = 3.000,
lpr = 1,
iobpl = 1,
ittyfreq = 20,
pltdt = 020.000,
prtdt = 3.0 ,
tddt = 5000.00,
velmx = 1.5,
burndef(1,1) = 1, 2, 1, 0, 0, 0, 0, 1000.0, 1.0, 5.0,
iburn = 2,
mat = 'n2' , 'o2' , 'h2',
'n2o', 'h2o', 'co',
'co2', 'nh3', 'ch4',
gasdef(1,1) = 0, 6, 1, -1, -1, -1, 1 ,
1.0000000000e6, 300.00, 1, 0., 0.,
'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, 5, 2, -2, -2, -2, 1 ,
1.0000000000e6, 300.00, 1, 0., 0.,
'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) = 0, 1, 1, -1, -1, -1, 1 ,
1.0000000000e6, 300.00, 1, 0., 99999.,
'n2' , 1.00 ,
gasdef(1,4) = 4, 5, 2, -1, -1, -1, 1 ,
1.0000000000e6, 300.00, 1, 0., 99999.,
'n2' , 1.00 ,
vvalue = 0.0,
vbc(1,1) = 1, 0, 1, 0, 0, 0, 0, 1, 0.0, 99999.,
vbc(1,2) = 2, 0, 2, 0, 0, 0, 0, 1, 0.0, 99999.,
$end
------------------------------------------------------------------------
M E S H
------------------------------------------------------------------------
$meshgn
$end
------------------------------------------------------------------------
G R A P H I C S
------------------------------------------------------------------------
$grafic
thdt = 0.0100,
thp(1,1) = 1, 1,-1,-1, 'pn' , 0,
thp(1,2) = 5, 1,-1,-1, 'pn' , 0,
thp(1,3) = 1, 2,-2,-2, 'tk' , 0,
thp(1,4) = 4, 2,-2,-2, 'tk' , 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,
You would like more information? Goto the
contact page.
