build
Home Page >>
Verification >>
build
EPA Windtunnel Simulation - 3D Block Building
TSA-4 Brown
------------------------------------------------------------------------
N O T E S
-----------
* FLOW AROUND SINGLE CUBIC 3D BUILDING
*********************************************************************
* 3-D domain 720 m x 420 m x 240 m
* Number of real cells = 40 x 30 x 22 = 26,400
* Cartesian coordinates, rectilinear non-uniform grid
* blgs. = 60 m long x 60m wide x 60 m high
* upstream = 180 m, downstream = 480 m
*********************************************************************
* Working fluid: air
* k-epsilon turbulence model
* van Leer numerical scheme turned on
* no gravity vector
*
* WEST B.C. (-x): specified inflow velocity (700 cm/s)
* EAST B.C. (+x): zero-gradient (continuative)
* SOUTH B.C. (-y): free-slip wall (side wall)
* NORTH B.C. (+y): free-slip wall (side wall)
* BOTTOM B.C. (-z): no-slip wall (ground-level surface)
* TOP B.C. (+z): free-slip wall (top of domain)
------------------------------------------------------------------------
------------------------------------------------------------------------
$xput
ifvl = 1,
tmodel = 'ke',
cmuke = 0.09,
c1ke = 1.44,
c2ke = 1.92,
sigmak = 1.0,
sigmae = 1.3,
idiffme = 1,
idiffmom = 1,
muoption = 1,
nu = 0.150 ,
; gz = -980.0,
maxcyc = 50000
twfin = 10.0,
autot = 1.0,
delt0 = 0.0300,
deltmin = 1.000e-05,
deltmax = 1.000e-00,
nrsdump = 0,
epsi0 = 1.000e-05,
epsimax = 1.000e-05,
epsimin = 1.000e-05,
itdowndt = 50,
itupdt = 100,
itmax = 500,
lpr = 1,
iobpl = 1,
ittyfreq = 50,
pltdt = 2.0,
prtdt = 10.0,
tddt = 5000.0,
velmx = 2.0,
mat = 'air' ,
gasdef(1,1) = 1, 41, 1, 31, 1, 23, 1 ,
1.0132500000e6, 300.00, 1, 0., 0.,
'air' , 1.00000,
gasdef(1,2) = 0, 1, 1, 31, 1, 23, 1 ,
1.0132500000e6, 300.00, 1, 0., 500.0,
'air' , 1.00000,
gasdef(1,3) = 41, 42, 1, 31, 1, 23, 1 ,
1.0132500000e6, 300.00, 1, 0., 500.0,
'air' , 1.00000,
mobs(1,1) = 11, 23, 10, 22, 1, 11, 1, 0,
ibw = 1,
ibe = 3,
ibn = 1,
ibs = 1,
ibb = 2,
ibt = 1,
ui=400.
vi=0.
wi=0.
nslipdef(1,1) = 11, 23, 10, 22, 1, 11, 1, 'both',
vvalue = 700.,
vbc(1,1) = 01, 01, 01,31, 01, 23, 1, 1, 0.0, 500.00,
$end
------------------------------------------------------------------------
M E S H
------------------------------------------------------------------------
$meshgn
iblock = 1,
nkx=3,
xl(1)= 0., xc(1)= 18000., nxl(1)= 10, nxr(1)= 0, dxmn(1)= 200.,
xl(2)= 18000., xc(2)= 18000., nxl(2)= 0, nxr(2)= 12, dxmn(2)= 200.,
xl(3)= 24000., xc(3)= 24000., nxl(3)= 0, nxr(3)= 18, dxmn(3)= 600.,
xl(4)= 72000.,
nky=4,
yl(1)=-21000., yc(1) =-2800., nyl(1)= 10, nyr(1)= 0, dymn(1)= 200.,
yl(2)= -2800., yc(2) =-2800., nyl(2)= 0, nyr(2)= 5, dymn(2)= 300.,
yl(3)= 0., yc(3) = 2800., nyl(3)= 5, nyr(3)= 0, dymn(3)= 300.,
yl(4)= 2800., yc(4) = 2800., nyl(4)= 0, nyr(4)= 10, dymn(4)= 200.,
yl(5)= 21000.,
nkz=2,
zl(1)= 0., zc(1)= 0., nzl(1)= 0, nzr(1)= 7, dzmn(1)= 200.,
zl(2)= 5000., zc(2)= 6000., nzl(2)= 3, nzr(2)= 12, dzmn(2)= 200.,
zl(3)= 24000.,
$end
------------------------------------------------------------------------
G R A P H I C S
------------------------------------------------------------------------
$grafic
pnt(1,1) = 2, 22, 01, 1,
pnt(1,2) = 2, 22, 23, 1,
pnt(1,3) = 9, 22, 01, 1,
pnt(1,4) = 9, 22, 23, 1,
pnt(1,5) = 17, 22, 01, 1,
pnt(1,6) = 17, 22, 23, 1,
pnt(1,7) = 26, 22, 01, 1,
pnt(1,8) = 26, 22, 23, 1,
pnt(1,9) = 38, 22, 01, 1,
pnt(1,10)= 38, 22, 23, 1,
p1d(1,1) = 1,2, 'un', 0
p1d(1,2) = 3,4, 'un', 0
p1d(1,3) = 5,6, 'un', 0
p1d(1,4) = 7,8, 'un', 0
p1d(1,5) = 9,10, 'un', 0
p1d(1,6) = 1,2, 'tke', 0
p1d(1,7) = 3,4, 'tke', 0
p1d(1,8) = 5,6, 'tke', 0
p1d(1,9) = 7,8, 'tke', 0
p1d(1,10)= 9,10, 'tke', 0
p1d(1,11)= 1,2, 'eps', 0
p1d(1,12)= 3,4, 'eps', 0
p1d(1,13)= 5,6, 'eps', 0
p1d(1,14)= 7,8, 'eps', 0
p1d(1,15)= 9,10, 'eps', 0
pnt(1,11) = 1, 22, 01, 1,
pnt(1,12) = 41, 22, 23, 1,
pnt(1,13) = 1, 1, 2, 1,
pnt(1,14) = 41, 31, 2, 1,
pnt(1,15) = 1, 1, 7, 1,
pnt(1,16) = 41, 31, 7, 1,
pnt(1,17) = 1, 1, 12, 1,
pnt(1,18) = 41, 31, 12, 1,
c2d(1,1) =11,12, 'un' , 0,
c2d(1,2) =11,12, 'vn' , 0,
c2d(1,3) =11,12, 'tke' , 0,
c2d(1,4) =11,12, 'eps' , 0,
c2d(1,5) =11,12, 'tk' , 0,
c2d(1,6) =11,12, 'pn' , 0,
c2d(1,7) =13,14, 'un' , 0,
c2d(1,8) =13,14, 'vn' , 0,
c2d(1,9) =13,14, 'tke' , 0,
c2d(1,10)=13,14, 'eps' , 0,
c2d(1,11)=13,14, 'tk' , 0,
c2d(1,12)=13,14, 'pn' , 0,
v2d(1,1) = 11, 12, 1,
v2d(1,2) = 13, 14, 1,
v2d(1,3) = 15, 16, 1,
v2d(1,4) = 17, 18, 1,
pnt(1,19) = 1, 1, 1, 1,
pnt(1,20) = 41, 31, 23, 1,
v3d(1,1) = 19, 20, 1, 1,
thdt = 0.5,
thp(1,1) = 3, 22, 2, 1, 'un' , 0,
thp(1,2) = 3, 22, 20, 1, 'un' , 0,
thp(1,3) = 40, 22, 2, 1, 'un' , 0,
thp(1,4) = 40, 22, 20, 1, 'un' , 0,
thp(1,5) = 3, 22, 2, 1, 'vn' , 0,
thp(1,6) = 3, 22, 20, 1, 'vn' , 0,
thp(1,7) = 40, 22, 2, 1, 'vn' , 0,
thp(1,8) = 40, 22, 20, 1, 'vn' , 0,
thp(1,9) = 3, 22, 2, 1, 'wn' , 0,
thp(1,10)= 3, 22, 20, 1, 'wn' , 0,
thp(1,11)= 40, 22, 2, 1, 'wn' , 0,
thp(1,12)= 40, 22, 20, 1, 'wn' , 0,
thp(1,13)= 3, 22, 2, 1, 'tke' , 0,
thp(1,14)= 40, 22, 2, 1, 'tke' , 0,
thp(1,15)= 3, 22, 2, 1, 'eps' , 0,
thp(1,16)= 40, 22, 2, 1, 'eps' , 0,
thp(1,17)= 3, 22, 2, 1, 'mu' , 0,
thp(1,18)= 40, 22, 2, 1, 'mu' , 0,
thp(1,19)= 3, 22, 20, 1, 'mu' , 0,
thp(1,20)= 40, 22, 20, 1, 'mu' , 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
$innet
$end
You would like more information? Goto the
contact page.
