声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2821|回复: 14

[分享]几个流体力学的例子

 关闭 [复制链接]
发表于 2005-12-16 20:04 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
/title, Compressible Flow in a Converging Diverging Nozzle<BR>/prep7<BR>! Assign Element Type<BR>!<BR>et,1,141<BR>!<BR>! Create the Model Geometry<BR>!<BR>/pnum,line,1<BR>!k,1,0,2.432<BR>k,2,0,0<BR>k,3,1,2.432<BR>k,4,1,0<BR>k,5,2,2.232<BR>k,6,2,0<BR>k,7,5,0<BR>k,8,5,.7<BR>k,9,6,0<BR>k,10,6,.5<BR>k,11,14,0<BR>k,12,14,.8<BR>l,2,4<BR>l,3,4<BR>l,1,3<BR>l,1,2<BR>l,5,8<BR>l2tan,3,5<BR>l,10,12<BR>l2tan,5,7<BR>l,4,6<BR>l,6,7<BR>l,7,9<BR>l,9,11<BR>a,2,4,3,1<BR>a,4,6,5,3<BR>a,6,7,8,5<BR>a,7,9,10,8<BR>a,9,11,12,10<BR>!<BR>! !Mesh the Model and Compress Numbers<BR>!<BR>ntran=24<BR>rtran=10<BR>na=10<BR>nb=10<BR>nc=25<BR>nd=10<BR>ne=100<BR>re=-1.5<BR>lesize,4,,,ntran,rtran<BR>lesize,2,,,ntran,rtran<BR>lesize,13,,,ntran,1/rtran<BR>lesize,14,,,ntran,1/rtran<BR>lesize,15,,,ntran,1/rtran<BR>lesize,16,,,ntran,1/rtran<BR>lesize,1,,,na<BR>lesize,3,,,na<BR>lesize,6,,,nb<BR>lesize,9,,,nb<BR>lesize,5,,,nc<BR>lesize,10,,,nc<BR>lesize,8,,,nd<BR>lesize,11,,,nd<BR>lesize,7,,,ne,re<BR>lesize,12,,,ne,re<BR>amesh,all<BR>numcmp,node<BR>numcmp,elem<BR>!<BR>! Apply Boundary Conditions<BR>!<BR>lsel,s,,,1<BR>lsel,a,,,9,12,1<BR>nsll,,1<BR>d,all,vy,0<BR>lsel,all<BR>lsel,s,,,4<BR>nsll,,1<BR>d,all,pres,6.13769e+6<BR>lsel,all<BR>lsel,s,,,3<BR>lsel,a,,,5,8,1<BR>nsll,,1<BR>d,all,vx,0<BR>d,all,vy,0<BR>lsel,all<BR>lsel,s,,,16<BR>nsll,,1<BR>d,all,pres,0<BR>!<BR>! Set Fluid Properties<BR>!<BR>alls<BR>/solu<BR>flda,prot,dens,air-cm<BR>flda,prot,visc,air-cm<BR>flda,prot,cond,air-cm<BR>flda,prot,spht,air-cm<BR>flda,vary,dens,t<BR>!<BR>! Set Flow Environment Parameters<BR>!<BR>flda,pres,refe,1.01325e+6<BR>flda,temp,ttot,550<BR>flda,temp,nomi,550<BR>!<BR>! Set Solution Options and CFD Solver Controls<BR>!<BR>flda,solu,turb,t<BR>flda,solu,comp,t<BR>flda,meth,pres,3<BR>flda,conv,pres,1e-12<BR>!<BR>! Solve<BR>!<BR>flda,relx,dens,1<BR>flda,stab,mome,1<BR>flda,stab,pres,1e-4<BR>flda,stab,visc,10<BR>flda,iter,exec,20<BR>solve<BR>flda,stab,visc,1<BR>solve<BR>flda,stab,visc,.1<BR>solve<BR>flda,stab,visc,.01<BR>solve<BR>flda,stab,visc,.001<BR>solve<BR>flda,stab,visc,.0001<BR>solve<BR>flda,stab,visc,.00001<BR>solve<BR>flda,stab,visc,0<BR>flda,iter,exec,160<BR>solve<BR>flda,stab,mome,10<BR>flda,stab,pres,1e-2<BR>flda,iter,exec,100<BR>solve<BR>flda,stab,mome,1e+15<BR>flda,stab,pres,1e+15<BR>flda,iter,exec,400<BR>solve<BR>
回复
分享到:

使用道具 举报

 楼主| 发表于 2005-12-16 20:05 | 显示全部楼层

回复:(meil32)[分享]几个流体力学的例子

<P>/BATCH,LIST<BR>/title,broken dam problem by the VOF solution.<BR>/com<BR>/com ------------------------------------------------------------------<BR>/com **  Reference: "Part IV: An Experimental Study of the Collapse of<BR>/com **              Liquid Columns on a Rigid Horizontal Plane" <BR>/com **              By Martin and Moyce, Phil. Trans. Roy. Soc. <BR>/com **              (London,A244, pp. 312-324, 1952.<BR>/com **              "Volume of Fluid (VOF) Method for the Dynamics of<BR>/com **              Free Boundaries" By C.W. Hirt and B.D. Nichols,<BR>/com **              J. Comput. Physics, Vol. 39, pp. 201-225, 1981.<BR>/com **<BR>/prep7  <BR>!!!!!   Define some dimensions and gravity              <BR>LX = 5.0                   !   x length<BR>LY = 2.5                   !   y length<BR>DELX = 0.2                 !   x spacing<BR>DELY = 0.05                !   y spacing<BR>grav = 9.81                !   gravity<BR>et,1,141<BR>!!!!  meshing<BR>rect,0,LX,0,LY<BR>lsel,s,,,2,4,2<BR>lesize,all,DELY<BR>lsel,s,,,1,3,2<BR>lesize,all,DELX<BR>amesh,all</P>
<P>nsel,s,loc,y,0.0                    ! wall boundary conditions     <BR>nsel,a,loc,x,0.0<BR>d,all,vx,0.0<BR>d,all,vy,0.0</P>
<P>nsel,s,loc,x,0.0,1.0                ! setup initial VFRC<BR>nsel,r,loc,y,0.0,2.0<BR>esln,s,1<BR>ice,all,vfrc,1.0</P>
<P>alls                                ! initial conditions for vx &amp; vy<BR>ic,all,vx,0<BR>ic,all,vy,0<BR>save<BR>fini</P>
<P>/solu<BR>!!!!  FLOTRAN input<BR>acel,0.0,grav,0<BR>FLDA,NOMI,DENS,1.0e3<BR>FLDA,NOMI,VISC,1.0d-3<BR>FLDA,SOLU,FLOW,T<BR>FLDA,SOLU,TRAN,T<BR>FLDA,SOLU,turb,F<BR>flda,solu,vof,t<BR>flda,time,glob,5<BR>flda,time,numb,70<BR>flda,time,step,0.01<BR>flda,time,appe,0.1<BR>flda,advm,mome,supg<BR>flda,relx,vx,0.5<BR>flda,relx,vy,0.5<BR>flda,relx,pres,1.0<BR>solve   <BR>fini<BR></P>
 楼主| 发表于 2005-12-16 20:05 | 显示全部楼层

回复:(meil32)[分享]几个流体力学的例子

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

回复:(meil32)[分享]几个流体力学的例子

<P>finish<BR>/clear</P>
<P>/filname,trans_test,1<BR>/prep7</P>
<P>et,1,141<BR>mp,ex,1,1<BR>mp,kxx,2,190<BR>mp,dens,2,2700<BR>mp,c,2,900</P>
<P>x0=0<BR>x1=15<BR>x21=x1+0.2<BR>x22=x21+0.3<BR>x23=x22+1.5<BR>x3=x23+20</P>
<P>thick=0.05<BR>y0=0<BR>y1=0.2<BR>y2=y1+0.15<BR>y3=10<BR>dt=0.0001</P>
<P>rect,x0,x1,y0,y3<BR>rect,x23,x3,y0,y2<BR>rect,x23,x3,y2,y2+thick<BR>rect,x23,x3,y2+thick,y3<BR>rect,x22,x23,y2,y2+thick<BR>rect,x22,x23,y2+thick,y3<BR>aglue,all</P>
<P>k,,x1+thick<BR>k,,x21,y1<BR>k,,x21,y1+thick</P>
<P>bsplin,2,13,20<BR>bsplin,17,10,9</P>
<P>a,2,20,17,9<BR>a,2,20,24,3</P>
<P>lsel,s,length,,0.05<BR>lesi,all,,,4,,,,,0</P>
<P>lsel,s,loc,x,x0+dt,x1-dt<BR>lesi,all,,,50,10</P>
<P>lsel,s,loc,x,x1+dt,x22-dt<BR>lesi,all,,,10,1</P>
<P>lsel,s,loc,x,x22+dt,x23-dt<BR>lesi,all,,,15,-4</P>
<P>lsel,s,loc,x,x23+dt,x3-dt<BR>lesi,all,,,60,10</P>
<P>lsel,s,loc,y,y2+dt,y3-dt<BR>lesi,all,,,40,10</P>
<P>lsel,s,loc,x,x23<BR>lsel,a,loc,x,x3<BR>lsel,r,loc,y,y0+dt,y1-dt<BR>lesi,all,,,8,-3</P>
<P>FLST,5,5,4,ORDE,5   <BR>FITEM,5,1   <BR>FITEM,5,4   <BR>FITEM,5,7   <BR>FITEM,5,11  <BR>FITEM,5,15  <BR>CM,_Y1,LINE <BR>LSEL, , , ,P51X <BR>*GET,_z1,LINE,,COUNT<BR>*SET,_z2,0  <BR>*DO,_z5,1,_z1   <BR>*SET,_z2,LSNEXT(_z2)<BR>*GET,_z3,LINE,_z2,ATTR,NDNX <BR>*GET,_z4,LINE,_z2,ATTR,SPNX <BR>*get,_z6,line,_z2,attr,kynd <BR>*IF,_z3,GT,0,THEN   <BR>*IF,_z4,NE,0,THEN   <BR>LESIZE,_z2,,,_z3,1/_z4,,,,_z6   <BR>*ENDIF  <BR>*ENDIF  <BR>*ENDDO  <BR>CMSEL,S,_Y1 <BR>CMDELE,_Y1  </P>
<P>alls<BR>MSHAPE,0,2D <BR>MSHKEY,1</P>
<P>asel,s,loc,x,x1,x23<BR>asel,r,loc,y,y0,y2+thick<BR>mat,2<BR>amesh,all<BR>asel,inve<BR>mat,1<BR>amesh,all</P>
<P>nsel,s,loc,y,0<BR>d,all,vy,0</P>
<P>nsel,s,loc,x,x23<BR>nsel,r,loc,y,y0,y2<BR>d,all,vx,0<BR>d,all,vy,0</P>
<P>nsel,s,loc,x,x0<BR>d,all,vx,0.8*20.05*(293**0.5)<BR>d,all,vy,0<BR>d,all,temp,1.16*293<BR>d,all,pres,0</P>
<P>nsel,s,loc,y,y3<BR>d,all,vy</P>
<P>nsel,s,loc,x,x3<BR>d,all,pres,0</P>
<P>lsel,s,,,8,10,2<BR>lsel,a,,,29<BR>dl,all,,temp,300,1<BR>alls<BR>tunif,293<BR>/solu</P>
<P>FLDATA1,SOLU,TRAN,1 <BR>FLDATA1,SOLU,FLOW,1 <BR>FLDATA1,SOLU,TEMP,1 <BR>FLDATA1,SOLU,TURB,1 <BR>FLDATA1,SOLU,COMP,1 <BR>FLDATA1,SOLU,VOF,0  <BR>FLDATA1,SOLU,SFTS,0 <BR>FLDATA1,SOLU,IVSH,0 <BR>FLDATA1,SOLU,SWRL,0 <BR>FLDATA1,SOLU,SPEC,0 <BR>FLDATA1,SOLU,ALE,0  <BR>FLDATA1,SOLU,RDSF,1 </P>
<P>/COM,,Transient Analysis,1  <BR>*SET,_z4,-999   <BR>*SET,_z14,1 <BR>!*  <BR>FLDATA4,TIME,STEP,0.01, <BR>FLDATA4,TIME,ISTEP,0<BR>FLDATA4,TIME,NUMB,1000, <BR>FLDATA4,TIME,TEND,1.0e06<BR>FLDATA4,TIME,GLOB,8,<BR>FLDATA4,TIME,VX,0.01,   <BR>FLDATA4,TIME,VY,0.01,   <BR>FLDATA4,TIME,VZ,0.01,   <BR>FLDATA4,TIME,PRES,1e-006,   <BR>FLDATA4,TIME,TEMP,1e-006,   <BR>FLDATA4,TIME,ENKE,0.01, <BR>FLDATA4,TIME,ENDS,0.01, <BR>FLDATA4A,STEP,OVER,10,  <BR>FLDATA4,TIME,OVER,0 <BR>FLDATA4A,STEP,APPE,10,  <BR>FLDATA4,TIME,APPE,1.0e6 <BR>FLDATA4A,STEP,SUMF,10,  <BR>FLDATA4,TIME,SUMF,1.0e6 <BR>FLDATA4,TIME,BC,0   <BR>FLDATA4,TIME,TEND,1000000,  </P>
<P>/COM,,Transient Analysis,1  <BR>FLDATA12,PROP,DENS,4<BR>FLDATA13,VARY,DENS,1<BR>FLDATA12,PROP,VISC,4<BR>FLDATA13,VARY,VISC,0<BR>FLDATA12,PROP,COND,4<BR>FLDATA13,VARY,COND,0<BR>FLDATA12,PROP,SPHT,4<BR>FLDATA13,VARY,SPHT,0<BR>!*  <BR>FLDATA7,PROT,DENS,AIR-SI<BR>FLDATA8,NOMI,DENS,-1<BR>FLDATA9,COF1,DENS,0 <BR>FLDATA10,COF2,DENS,0<BR>FLDATA11,COF3,DENS,0<BR>FLDATA7,PROT,VISC,AIR-SI<BR>FLDATA8,NOMI,VISC,-1<BR>FLDATA9,COF1,VISC,0 <BR>FLDATA10,COF2,VISC,0<BR>FLDATA11,COF3,VISC,0<BR>FLDATA12,PROP,IVIS  <BR>FLDATA7,PROT,COND,AIR-SI<BR>FLDATA8,NOMI,COND,-1<BR>FLDATA9,COF1,COND,0 <BR>FLDATA10,COF2,COND,0<BR>FLDATA11,COF3,COND,0<BR>FLDATA7,PROT,SPHT,AIR-SI<BR>FLDATA8,NOMI,SPHT,-1<BR>FLDATA9,COF1,SPHT,0 <BR>FLDATA10,COF2,SPHT,0<BR>FLDATA11,COF3,SPHT,0</P>
<P>!FLDA,STAB,VISC,10       <BR>FLDA,STAB,MOME,1.0      <BR>FLDA,STAB,TURB,1      <BR>FLDA,STAB,PRES,1.e-3    <BR>FLDA,METH,PRES,3     <BR>FLDA,SRCH,PRES,10     <BR>FLDA,MAXI,PRES,1000     <BR>flda,capp,pres,t<BR>flda,capp,pmin,-9.5E+4<BR>solve<BR></P>
 楼主| 发表于 2005-12-16 20:06 | 显示全部楼层

回复:(meil32)[分享]几个流体力学的例子

<P><BR>/batch,list<BR>/filename,spec<BR>/com ** Analysis Types Laminar and Turbulent, Incompressible<BR>/com       2D - quadrilateral and triangular elements<BR>/com ** Features    Thermal, Multiple Species Transport<BR>/com ** Options   Dilute mixtures, composite property types (CGAS, CMIX)<BR>/com ** Variable laminar properties<BR>/com ** Turbulence<BR>/com ** Construction of geometry<BR>/prep7<BR>!!!!!!!  Define some dimensions - SI unite are used<BR>lenin=.3<BR>half=lenin/2.<BR>width=.1<BR>hfwid=width/2.<BR>outlen=1.2<BR>!!!!!!!  Define inlet and outlet rectangles<BR>rect,-lenin,0,-hfwid,hfwid<BR>rect,.2,.3,.25,.25+lenin<BR>rect,.2,.3,-.25,-.25-lenin<BR>rect,.4,.4+outlen,-.15,.15<BR>!!!!!!!  Lines that border the mixing area<BR>I2tan,-3,-8<BR>I2tan,-6,-15<BR>I2tan,-13,-10<BR>I2tan,-12,-1<BR>!!!!!!!  Parameters defined for meshing<BR>nlcurv=16<BR>rlcurv=1<BR>nscurv=9<BR>nispan=8<BR>rispan=-1.5<BR>nospan=13<BR>rospan=.2<BR>nilen=10<BR>rilen=-2<BR>nolen=24<BR>rolen=4<BR>!!!!!!!<BR>!flst,2,2,4,orde,2<BR>!fitem,2,17<BR>!fitem,2,20<BR>lsel,s,,,17,20,3<BR>lesize,all,,,nlcurv,rlcurv<BR>lsel,s,,,18,19<BR>!fitem,2,18<BR>!fitem,2,-19<BR>lesize,all,,,nscurv,rscurv<BR>!flst,2,3,4,orde,3<BR>!fitem,2,2<BR>!fitem,2,5<BR>!fitem,2,11<BR>lsel,s,,,2,5,3<BR>lsel,a,,,11<BR>lesize,all,,,nispan,1<BR>!flst,2,3,4,orde,3<BR>!fitem,2,4<BR>!fitem,2,7<BR>!fitem,2,9<BR>lsel,s,,,4,7,3<BR>lsel,a,,,9<BR>lesize,all,,,nispan,rispan<BR>!flst,2,6,4,orde,6<BR>!fitem,2,1<BR>!fitem,2,3<BR>!fitem,2,6<BR>!fitem,2,8<BR>!fitem,2,10<BR>!fitem,2,12<BR>lsel,s,,,1,3,2<BR>lsel,a,,,6,12,2<BR>lesize,all,,,nilen,rilen<BR>lsel,s,,,13,16<BR>lesize,13,,,nolen,rolen<BR>lesize,15,,,nolen,1./rolen<BR>lesize,16,,,nospan,rospan<BR>lesize,14,,,nospan,rospan<BR>alls<BR>!!!!!!!  Define the mixing area...<BR>!flst,2,8,3<BR>!fitem,2,2<BR>!fitem,2,12<BR>!fitem,2,11<BR>!fitem,2,13<BR>!fitem,2,16<BR>!fitem,2,6<BR>!fitem,2,5<BR>!fitem,22,3<BR>a,2,12,11,13,16,6,5,3<BR>alls<BR>/com ** CONSTRUCTION OF THE MESH<BR>!!!!!!! Put triangles in the mixing area (5)<BR>asel,s,,,5<BR>et,1,141<BR>mshape,1,2d<BR>amesh,5<BR>!!!!!!! Put a mapped mesh in the rectangles (quads)<BR>asel,s,,,1,4<BR>mshape,0,2d<BR>mshkey,1<BR>amesh,all<BR>/com ** BOUNDARY CONDITIONS<BR>!!!!!!! Wall boundary conditions<BR>lesel,s,,,,1,3,2<BR>lsel,a,,,6,12,2<BR>lsel,a,,,13,17,2<BR>lsel,a,,,18,20<BR>nsll,,1<BR>d,all,vx,0<BR>d,all,vy,0<BR>!!!!!!! Define velocities, temperature at top, bottom, left<BR>vtop=.1<BR>vbot=.1<BR>vlef=.1<BR>lsel,s,,,7<BR>nsll<BR>d,all,vx,0<BR>d,all,vy,-vtop<BR>nsll,,1<BR>d,all,temp,300<BR>lsel,s,,,4<BR>nsll<BR>d,all,vx,vlef<BR>d,all,vy,0<BR>nsll,,1<BR>d,all,temp,400<BR>lsel,s,,,9<BR>nsll<BR>d,all,vx,0<BR>d,all,vy,vbot<BR>nsll,,1<BR>d,all,temp,300<BR>lsel,s,,,14<BR>nsll,,1<BR>d,all,pres,0<BR>flda,conv,iter,5<BR>save<BR>finish<BR>alls<BR>/solu<BR>/com **  FLOTRAN INPUT<BR>flda,iter,exec,40<BR>flda,temp,nomi,300 ! Initial temperature field is 300K<BR>flda,prot,dens,air-si ! Initial solution will be for AIR<BR>flda,vary,dens,true<BR>flda,prot,visc,air-si<BR>flda,vary,visc,true<BR>flda,prot,cond,air-si<BR>flda,vary,cond,true<BR>flda,prot,spht,air-si<BR>save<BR>solve ! Solve 40 iterations with air at 300K<BR>finish ! Prepare for species transport...<BR>/prep7<BR>flda,prot,dens,cmix ! Composite mixture for density<BR>flda,prot,visc,cmix ! Composite mixture for viscosity<BR>keyopt,1,1,3 ! Specify that there are 3 species</P>
<P>flda,solu,spec,t ! Turn on solution for multiple species transport<BR>msdata,2 ! Specify species 2 (N2) as the algebraic species<BR>!<BR>!!!!!!!  Species Property Input<BR>msspec,1,o2,31.999<BR>msprop,1,DENS,GAS,1.2998,300,1.01325E+5<BR>msvary,1,dens,t<BR>msnomf,1,.3<BR>msprop,1,VISC,CONSTANT,1.2067E-5<BR>msprop,1,mdif,CONSTANT,2.149E-5<BR>msprop,1,cond,CONSTANT,.02674<BR>mscap,1,1<BR>!<BR>msspec,2,n2,28.018<BR>msprop,2,DENS,GAS,1.1138,300,1.01325E+5<BR>msvary,2,dens,t<BR>msnomf,2,.3<BR>msprop,2,VISC,CONSTANT,1,786E-5<BR>msvary,3,dens,t<BR>msnomf,3,.4<BR>msprop,3,VISC,CONSTANT,8.94E-6<BR>msprop,3,mdif,CONSTANT,4.964E-5<BR>msprop,3,cond,CONSTANT,.1815<BR>mscap,3,1<BR>!<BR>msrelx,1,1.0<BR>msrelx,3,1.0<BR>msmeth,1,3<BR>mssolu,1,,,2,1.e-8<BR>alls<BR>!!!!!!!  Species boundary conditions<BR>lsel,s,,,7<BR>nsll,,1<BR>d,all,o2,1.0<BR>d,all,n2,0.0<BR>d,all,h2,0.0<BR>lsel,s,,,4<BR>nsll,,1<BR>d,all,o2,0.0<BR>d,all,n2,1.0<BR>d,all,h2,0.0<BR>lsel,s,,,9<BR>nsll,,1<BR>d,all,o2,0.0<BR>d,all,n2,0.0<BR>d,all,h2,1.0<BR>alls<BR>!!!!!!!  End of species boundary conditions<BR>flda,iter,exec,20 ! Ask for 20 global iterations<BR>flda,conv,outp,land ! Adjust convergence monitor output style<BR>save<BR>finish<BR>/solu<BR>solve ! 20 iterations with species activated<BR>!<BR>! Prepare for energy solution<BR>flda,solu,temp,t   ! Achieve a constant flow temperature solution<BR>flda,solu,flow,f   ! Freeze the flowfield<BR>flda,meth,temp3   ! Activate PCCR solver<BR>flda,conv,temp,1.e-10   ! Convergence criterion for PCCR<BR>flda,iter,exe,5   ! Need only a few iterations for solution...<BR>flda,relx,temp,1.0   ! no relaxation on temperature<BR>save<BR>solve<BR>! ! Prepare for coupled solution</P>
<P>flda,iter,exec,35   ! Achieve a flow and thermal solution<BR>flda,solu,flow,t<BR>solve<BR>flda,iter,exec,50   ! Run 50 more iterations to refine solution<BR>save<BR>solve<BR>finish<BR></P>
发表于 2005-12-17 01:26 | 显示全部楼层
搂住可否提供例题说明
发表于 2005-12-28 19:36 | 显示全部楼层

谢谢

谢谢!收下了[em02]
发表于 2006-1-26 10:42 | 显示全部楼层
谢谢楼主
发表于 2006-2-23 19:26 | 显示全部楼层

搂住可否提供例题说明

搂住可否提供例题说明
发表于 2006-3-2 13:53 | 显示全部楼层
全是命令啊,怎么说也得差不多给个说明啊!
发表于 2006-4-1 10:11 | 显示全部楼层
<P>求示楼主加说明附件</P>
发表于 2006-5-31 12:48 | 显示全部楼层
<P>看不懂啊</P>
发表于 2006-5-31 14:08 | 显示全部楼层
谢谢
发表于 2006-6-3 20:26 | 显示全部楼层
<P>真的是好东西啊<BR>向共享广义者下跪!!!!</P>
发表于 2006-9-1 19:43 | 显示全部楼层
为什么不打包放在一块呢?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 21:26 , Processed in 0.068988 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表