2plates
Home Page >>
Verification >>
2plates
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
You would like more information? Goto the
contact page.
