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