2D turbulent flow between two parallel plates TSA10 JWS GASFLOW V2.1 NOTES: $Source: /home/cvs/cvsroot/input/2plates.ingf,v $ $Date: 1998/09/29 14:43:32 $ $Locker: $ $Revision: 1.2 $ $State: Exp $ $Author: necker $ $Log: 2plates.ingf,v $ Revision 1.2 1998/09/29 14:43:32 necker Changes from Jay Spore/LANL: Changed the input to Re = 10000 and increased length of the mesh to 50*H or 2500 cm to ensure that fully developed profile is obtained. Empirical solution for wall shear coefficient for two plates is according to Wilcox, "Turbulence Modeling for CFD" is: Cf = 0.0706(ReH)**(-0.25) = tauw/(0.5*rho*Uavg**2) Cf = 0.0706(10000)**(-0.25) = 0.00706 GASFLOW solution at wall at k=96, j=2, i=2 Cf = nu*(U(2)-U(1))/(0.5*dely(2)*Uavg**2) = 0.153*1.1262/(0.5*0.05*30.6**2) = 0.00736 Turned van Leer on and changed wallset to 6. Revision 1.1 1998/07/30 14:03:46 necker Added test problem from Jay Spore/LANL: Turbulent flow between two flat plates. * 1 3D block -- 2D flow in the y-z direction * air * BLOCK #1 2 x 50 x 2500 cm 1 x 50 x 100 = 5000 real cells * Top boundary is fixed pressure. * Bottom boundary is fixed inlet velocity, * North and South boundaries are rigid no-slip. * West and East boundaries are rigid free-slip. * Re = 30.6*50/.153 = 10000.0 * Rex = 30.6*x/.153 = 200*x * L = 50*H = 2500, required to get fully developed velocity profile. ----------------------------------------------------------------------- N E T W O R K I N P U T ----------------------------------------------------------------------- $innet $end ----------------------------------------------------------------------- B A S I C I N P U T ----------------------------------------------------------------------- $xput pamb0 = 1.0e+06, idiffmom = 1, impdiff = 1, tmodel = 'ko', ifvl = 1, autot = 1.0, cyl = 0.0, delt0 = 0.01000 , deltmin = 1.000e-08, deltmax = 1.000e-01, epsi0 = 1.000e-08, epsimax = 1.000e-08, epsimin = 1.000e-08, tkeamb0 = 0.01, sclamb0 = 1.0, gz = 0.0, iobpl = 0, itdowndt = 100, itupdt = 95, itmax = 150, lpr = 1, maxcyc = 100000, ittyfreq = 5, nu = 0.153, pltdt = 10.00, prtdt = 50.000, twfin = 50.00, tddt = 1000.0000, velmx = 1.5, ibw = 1, ibe = 1, ibs = 2, ibn = 2, ibb = 1, ibt = 5, wi = 15.3, iwallfunc = 1, wallset = 6, mat = 'air', vvalue = 30.6, vbc(1,1) = 1, 'im1', 1, 'jm1', 1, 1, 1, 1, 0.0, 1.0e+8, gasdef(1,1) = 1 ,'im1', 1 ,'jm1', 1 ,'km1', 1 , ; I.C. 1.0e6, 300.0, 1, 0., 0., 'air', 1.00000, gasdef(1,2) = 0 , 1, 1 ,'jm1', 01 ,'km1', 1 , ; west B.C. 1.0e6, 300.0, 1, 0.,1.0e+99, 'air', 1.00000, gasdef(1,3) = 'im1','imax', 1,'jm1', 1 ,'km1', 1 , ; east B.C. 1.0e6, 300.0, 1, 0.,1.0e+99, 'air', 1.00000, gasdef(1,4) = 1,'im1', 0, 1, 1 ,'km1', 1 , ; south B.C 1.0e6, 300.0, 1, 0.,1.0e+99, 'air', 1.00000, gasdef(1,5) = 1,'im1','jm1','jmax', 1 ,'km1', 1 , ; north B.C 1.0e6, 300.0, 1, 0.,1.0e+99, 'air', 1.00000, gasdef(1,6) = 1,'im1', 1,'jm1', 0 , 1, 1 , ; bttm B.C. 1.0e6, 300.0, 1, 0.,1.0e+99, 'air', 1.00000, gasdef(1,7) = 1,'im1', 1,'jm1','km1','kmax', 1 , ; top B.C. 1.0e6, 300.0, 1, 0.,1.0e+99, 'air', 1.00000, $end ------------------------------------------------------------------------ M E S H ------------------------------------------------------------------------ $meshgn iblock = 1, nkx=1, nky = 2, nkz = 1, xl(1) = 0.0, xc(1) = 0.0, nxl(1)= 0, nxr(1)= 1, dxmn(1)=1.e9, xl(2)= 2.0, yl(1) = 0.0, yc(1) = 0.0, nyl(1)= 0, nyr(1)= 25, dymn(1)=0.05, yl(2)= 25.0, yc(2) = 50.0, nyl(2)=25, nyr(2)= 0, dymn(2)=0.05, yl(3)= 50.0, zl(1) = 0.0, zc(1) = 0.0, nzl(1)= 0, nzr(1)=100, dzmn(1)=1.e9, zl(2)=2500.0, $end $meshgn $end ------------------------------------------------------------------------ G R A P H I C S ------------------------------------------------------------------------ $grafic pnt(1, 1) = 2, 1, 40, 1, 2, 'jmax', 40, 1, 2, 1, 2, 1, 2, 'jmax', 2, 1, 2, 1, 1, 1, 2, 'jm1', 'km1', 1, 2, 1, 96, 1, 2, 'jmax', 96, 1, 2, 1, 96, 1, 2, 'jm1', 'km1', 1, p1d = 1, 2, 'wn', 0, 3, 4, 'wn', 0, 7, 8, 'wn', 0, 1, 2, 'tke', 0, 3, 4, 'tke', 0, 7, 8, 'tke', 0, 1, 2, 'eps', 0, 3, 4, 'eps', 0, 7, 8, 'eps', 0, 1, 2, 'mu', 0, 3, 4, 'mu', 0, 7, 8, 'mu', 0, v2d = 5, 6, 1, 9, 10, 1, thdt = 0.5, thp(1,1) = 2, 16, 2, 1, 'wn', 0, 2, 16, 40, 1, 'wn', 0, 2, 16, 98, 1, 'wn', 0, 2, 17, 98, 1, 'wn', 0, 2, 17, 40, 1, 'wn', 0, 2, 17, 2, 1, 'wn', 0, 2, 18, 98, 1, 'wn', 0, 2, 18, 40, 1, 'wn', 0, 2, 18, 2, 1, 'wn', 0, 2, 2, 2, 1, 'tke', 0, 2, 51, 2, 1, 'tke', 0, 2, 2, 40, 1, 'tke', 0, 2, 51, 40, 1, 'tke', 0, 2, 2, 98, 1, 'tke', 0, 2, 51, 98, 1, 'tke', 0, 2, 2, 2, 1, 'eps', 0, 2, 51, 2, 1, 'eps', 0, 2, 2, 40, 1, 'eps', 0, 2, 51, 40, 1, 'eps', 0, 2, 2, 98, 1, 'eps', 0, 2, 51, 98, 1, 'eps', 0, 2, 2, 2, 1, 'mu', 0, 2, 51, 2, 1, 'mu', 0, 2, 2, 40, 1, 'mu', 0, 2, 51, 40, 1, 'mu', 0, 2, 2, 98, 1, 'mu', 0, 2, 51, 98, 1, 'mu', 0, 2, 2, 2, 1, 'nu', 0, 2, 51, 2, 1, 'nu', 0, 2, 2, 40, 1, 'nu', 0, 2, 51, 40, 1, 'nu', 0, 2, 2, 98, 1, 'nu', 0, 2, 51, 98, 1, 'nu', 0, $end ------------------------------------------------------------------------ P A R T I C L E S ------------------------------------------------------------------------ $parts $end $rheat $end $special $end $specialp $end