/BATCH,LIST /title, initial development of flow over a bump. /com /com ------------------------------------------------------------------ /com ** Initial development of flow over a semicircular bump /com ** Reference: "Critical Free-Surface Flow Over a Semi-Circular /com ** Obstruction" By L.K. Forbes, J. Eng. Math., /com ** vol. 22, pp. 3-13, 1998. /com ** "Computation of free-surface flows using the /com ** finite-volume method and moving grids", /com ** By S. Muzaferija and M. Peric, Numerical Heat /com ** Transfer, Part B, 32:369-384, 1997. /com ------------------------------------------------------------------ ! ! ! Property and parameter input rho = 1.0e0 ! Density grav = 9.81 ! Gravity mu = 1.0e-6 ! Viscosity ! R = 0.03 ! Bump radius L = 2*R H = 0.075 ! Inlet height H0 = R ! Height above the Inlet outlen = 5.0*R ! Outlet length inlen = 5.0*R ! Inlet length ny1 = 8 ! Number of elements in line 1 ny2 = 2 ! Number of elements in line 2 ny3 = 4 ! Number of elements in line 3 ry2 = 1.0 ! Ratio of element size in line 2 ry3 = 1.0 ! Ratio of element size in line 3 nx1 = 15 ! Number of elements in line 13 nx2 = 15 ! Number of elements in line 21 nx3 = 20 ! Number of elements in line 17 rx1 = 0.5 ! Ratio of element size in line 13 rx2 = 1.0 ! Ratio of element size in line 25 rx3 = 1.5 ! Ratio of element size in line 17 nr1 = 10 ! Number of elements in radial direction near bump !! /prep7 et,1,141 k,1,-inlen,0.0 ! KP 1 - 4 are at the inlet k,2,-inlen,L k,3,-inlen,H k,4,-inlen,H+H0 ! l,1,2 ! Lines 1-3 are the inlet line l,2,3 l,3,4 ! k,5,-L,0.0 ! KP 5 - 8 are before the bump k,6,-L,L k,7,-L,H k,8,-L,H+H0 ! l,5,6 ! Lines 4-6 before the bump l,6,7 l,7,8 ! k,9,L,0.0 ! KP 9 - 12 are after the bump k,10,L,L k,11,L,H k,12,L,H+H0 ! l,9,10 ! Lines 7-9 after the bump l,10,11 l,11,12 ! k,13,outlen,0.0 ! KP 13 - 16 are at the outlet k,14,outlen,L k,15,outlen,H k,16,outlen,H+H0 ! l,13,14 ! Lines 10-12 at the outlet l,14,15 l,15,16 ! l,1,5 ! Lines 13-16 are horizontal lines near inlet l,2,6 l,3,7 l,4,8 ! l,9,13 ! Lines 17-20 are horizontal lines near outlet l,10,14 l,11,15 l,12,16 ! l,6,10 ! Lines 21-23 are horizontal lines near the bump l,7,11 l,8,12 ! k,17,R,0.0 ! KP 17 to 23 are for the semicircular bump k,18,0.0,0.0 circle,18,R,,17,180,4 lcomb,25,26,0 ! l,19,9 ! Lines 26, 28-30 in radial direction near bump l,20,10 l,22,6 l,23,5 ! ! Select lines and establish line divisions ! lsel,s,,,1,4,3 lesize,all,,,ny1 ! lsel,s,,,2,11,3 lesize,all,,,ny2,ry2 ! lsel,s,,,3,12,3 lesize,all,,,ny3,ry3 ! lsel,s,,,13,16 lesize,all,,,nx1,rx1 ! lsel,s,,,21,23 lesize,all,,,nx2,1.0/rx2 lsel,s,,,25 lesize,all,,,nx2,rx2 ! lsel,s,,,17,20 lesize,all,,,nx3,rx3 ! lsel,s,,,27 lesize,all,,,ny1 ! lsel,s,,,26 lsel,a,,,28,30 lesize,all,,,nr1 ! ! define the y-direction mesh lsel,s,,,24 lsel,a,,,7 lsel,a,,,10 lesize,all,,,nx2,1.5 ! alls ! Create areas.... a,1,5,6,2 ! Inlet a,2,6,7,3 a,3,7,8,4 a,9,13,14,10 ! Outlet a,10,14,15,11 a,11,15,16,12 a,6,10,11,7 ! Bump a,7,11,12,8 a,19,9,10,20 a,20,10,6,22 a,5,23,22,6 alls amesh,all ! ! Inlet boundary condition lsel,s,,,1,2 nsll,,1 d,all,vx,0.32 d,all,vy,0.0 ! nsel,s,loc,x,-inlen nsel,r,loc,y,0.0 d,all,vx,0.001 ! ! boundary volume fraction (SFL and SFE commands) sfl,1,vfrc,1,,1 sfl,2,vfrc,1,,1 !! !!!! The walls lsel,s,,,13 lsel,a,,,17 lsel,a,,,24,25 lsel,a,,,27 lsel,a,,,26 lsel,a,,,30 nsll,,1 d,all,vx,0 d,all,vy,0 alls /solu !!!! FLOTRAN Input acel,0.0,grav,0.0 FLDA,SOLU,VOF,ON FLDA,SOLU,FLOW,ON FLDA,SOLU,TRAN,ON FLDA,SOLU,TURB,ON FLDA,TIME,GLOB,3 FLDA,TIME,NUMB,50 FLDA,TIME,STEP,1e-2 FLDA,TIME,APPE,0.1 FLDA,TIME,PRES,1.0E-6 FLDA,ADVM,MOME,SUPG flda,relx,vx,0.5 flda,relx,vy,0.5 flda,relx,pres,1.0 flda,prot,dens,constant ! Constant density flda,prot,visc,constant ! Constant viscosity flda,nomi,dens,rho ! Density set earlier flda,nomi,visc,mu ! Viscosity set earlier ! solve