1D/3D A Problems
Home Page >>
Verification >>
1D/3D A Problems
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
You would like more information? Goto the
contact page.
