PHOTON USE p parphi 1 20 1 up z vec x 1 sh pause con wcrt x 1 fi;0.001 pause con h1 x 1 fi;0.001 ENDUSE GROUP 1. Run title and other preliminaries TEXT(Free Convec From Hot Vertical Wall TITLE DISPLAY This run is similar to the run above,except that the flow is driven by a buoyancy force rather than by a prescribed mass inflow. As a consequence,the grid expands at a different rate. | ^ ^ |/ ^ ^ ^ | | ^ |/ Pressure, | ^ ^ ^ ^ | | | | | | ^ |/ velocity |/ and | |/ enthalpy |/ fixed to | | |/ zero |gravity |/ | v |/ Wall; |/ heated with H1=1.0 | |/ |/ ^ | |/ | z |/ <--| | ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ |/ y ENDDIS Interesting variants include changes to: the magnitude of the added force; and to the wall boundary values. GXBUOY is used for this case. GROUP 4. Y-direction grid specification GRDPWR(Y,20,4.0E-3,1.0) *** Grid grows as (axial distance)**0.25 AZYV=0.25;ZWADD=1.0E-03 GROUP 5. Z-direction grid specification *** The flow is parabolic PARAB=T *** Axial grid distance varies as IZ**4.0 GRDPWR(Z,100,0.1,4.0) AZDZ=GRND (ZGRID OF DZ IS 0.5*YVLAST) GROUP 7. Variables stored, solved & named SOLVE(P1,V1,W1,H1) GROUP 8. Terms (in differential equations) & devices *** Deactivate the built-in source term in H1 equation TERMS(H1,N,Y,Y,Y,Y,Y) GROUP 9. Properties of the medium (or media) ENUL=1.E-5; PRNDTL(H1)=0.7 CP1 = 1.005000E+03 GROUP 11. Initialization of variable or porosity fields FIINIT(V1)=0.0; FIINIT(W1)=0.0; FIINIT(H1)=0.0 GROUP 13. Boundary conditions and special sources 1. Wall boundary: constant temperature WALL (CONST,SOUTH,1,1,1,1,1,NZ,1,1) COVAL(CONST,W1,1.0,0.0) COVAL(CONST,H1,1.0,1.0) 2. Free stream boundary PATCH(HIGHY,NORTH,1,1,NY,NY,1,NZ,1,1) COVAL(HIGHY,P1,FIXP,0.0);COVAL(HIGHY,W1,ONLYMS,0.0) COVAL(HIGHY,H1,ONLYMS,0.0) 3. Buoyancy source *** The added source in the Z-direction momentum equation is RHO*Volume*grav.*dvo1dt/cp*H1 (See GREX3, Group 13,Sec.15, and GXBUOY); the following data need to be set for this purpose: BUOYE=H1ref; BUOYC=grav. BUOYE=0.0; BUOYC=-9.81 DVO1DT=1.005E+02 PATCH(BUOYANCY,PHASEM,1,1,1,NY,1,NZ,1,1) COVAL(BUOYANCY,W1,FIXFLU,BOUSS) GROUP 14. Downstream pressure for PARAB=T *** The downstream pressure is set by the parameter PBAR with its default value of 0.0. IPARAB=1 GROUP 16. Termination of iterations LITHYD=10 GROUP 17. Under-relaxation devices RELAX(W1,FALSDT,0.1); RELAX(V1,FALSDT,0.1) GROUP 22. Spot-value print-out NPRMON=NZ/10; IYMON=3; NPLT=1; TSTSWP=2 GROUP 23. Field print-out and plot control NZPRIN=NZ/4 PATCH(FIXEDZ,PROFIL,1,1,1,NY,1,NZ,1,1) PLOT (FIXEDZ,W1,0.0,0.0); PLOT(FIXEDZ,H1,0.0,0.0) IDISPA=1;IDISPB=10;IDISPC=NZ #CONPROM #MAXABS #ENDPAUSE TSTSWP=-1