Library case H121:
THE EFFECT OF ENCLOSURE ORIENTATION
talk=t;run(1,1) PHOTON USE p use patgeo up .54 .84 0 red scal msg msg Temperature contours msg cont tem1 z 1 fil;.01 msg Pressto continue pause msg Velocity field vec z 1 sh msg MSG Type E to end ENDUSE DISPLAY Library case H121: THE EFFECT OF ENCLOSURE ORIENTATION This case models simple two-dimensional flow through a tilted box (angle is 1 rad). Air has two apertures for flow and the buoyancy force is generated from the temperature difference; single solid in the middle of the domain generates heat. ENDDIS ************************************************************ CPVNAM=VDI;SPPNAM=HOTBOX ************************************************************ IRUNN = 1 ;LIBREF = 121 ************************************************************ Group 1. Run Title TEXT(The Effect Of Enclosure Orientation ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,10,10,1) * Set overall domain extent: * xulast yvlast zwlast name XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00 RSET(D,BOX ) * Set objects: x0 y0 z0 * dx dy dz name XPO= 4.000000E-01; YPO= 4.000000E-01; ZPO= 0.000000E+00 XSI= 2.000000E-01; YSI= 2.000000E-01; ZSI= 1.000000E+00 RSET(B,BLK ) XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00 XSI= 0.000000E+00; YSI= 2.000000E-01; ZSI= 1.000000E+00 RSET(B,APT1 ) XPO= 8.000000E-01; YPO= 1.000000E+00; ZPO= 0.000000E+00 XSI= 2.000000E-01; YSI= 0.000000E+00; ZSI= 1.000000E+00 RSET(B,APT2 ) * Modify default grid RSET(X,1,4,1.000000E+00,G) RSET(X,2,2,1.000000E+00) RSET(X,3,2,1.000000E+00) RSET(X,4,2,1.000000E+00) RSET(Y,1,2,1.000000E+00,G) RSET(Y,2,2,1.000000E+00) RSET(Y,3,2,1.000000E+00) RSET(Y,4,4,1.000000E+00) ************************************************************ Group 6. Body-Fitted coordinates * X-cyclic boundaries switched ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME( 46) =BLOK ; NAME( 47) =PRPS NAME( 49) =DEN1 ; NAME( 50) =TEM1 * Solved variables list SOLVE(P1 ,U1 ,V1 ,TEM1) * Stored variables list STORE(DEN1,PRPS,BLOK) * Additional solver options SOLUTN(TEM1,Y,Y,Y,N,N,Y) IVARBK = -1 ;ISOLBK = 1 ************************************************************ Group 8. Terms & Devices NEWRH1 = T NEWENL = T ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0 ************************************************************ Group 9. Properties SETPRPS(1, 2) RHO1 = GRND5 RHO1A = 0.000000E+00 ;RHO1B = 3.484300E-03 RHO1C = 7.142900E-01 PRESS0 = 1.000000E+05 TEMP0 = 2.730000E+02 CP1 = 1.004000E+03 ENUL = GRND2 ;ENUT = 0.000000E+00 ENULA =-4.946790E-06 ;ENULB = 4.583940E-08 ENULC = 8.097400E-11 DRH1DP = GRND5 DVO1DT = 3.330000E-03 PRNDTL(TEM1) = -2.630000E-02 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(BLOK) = 1.000000E+00 ;FIINIT(PRPS) = 2.000000E+00 FIINIT(TEM1) = 1.000000E+01 No PATCHes used for this Group INIADD = F ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1) COVAL (BUOYANCY,U1 , FIXFLU , GRND2 ) COVAL (BUOYANCY,V1 , FIXFLU , GRND2 ) BUOYA =-5.300000E+00 ; BUOYB =-8.255000E+00 BUOYC = 0.000000E+00 BUOYD = 1.231200E+00 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 100 RESFAC = 1.000000E-02 ************************************************************ Group 16. Terminate Iterations LITER (TEM1) = 15 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) RELAX(U1 ,FALSDT, 5.000000E-03) RELAX(V1 ,FALSDT, 5.000000E-03) RELAX(TEM1,FALSDT, 1.000000E+04) ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 ************************************************************ Group 19. EARTH Calls To GROUND Station NAMGRD =HTBX ASAP = T ************************************************************ Group 20. Preliminary Printout ECHO = T ************************************************************ Group 21. Print-out of Variables OUTPUT(BLOK,N,N,N,N,N,N) OUTPUT(PRPS,N,N,N,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 4 ;IYMON = 5 ;IZMON = 1 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts NOWIPE = T GVIEW(P,-2.999549E-02,-1.079130E-02,9.994916E-01) GVIEW(UP,4.879574E-01,8.725354E-01,2.406454E-02) > DOM, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, MONIT, 3.500000E-01, 4.500000E-01, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > DOM, RELAX, 1.000000E+00 > OBJ1, NAME, BLK > OBJ1, POSITION, 4.000000E-01, 4.000000E-01, 0.000000E+00 > OBJ1, SIZE, 2.000000E-01, 2.000000E-01, 1.000000E+00 > OBJ1, CLIPART, cube4 > OBJ1, ROTATION24, 1 > OBJ1, TYPE, BLOCKAGE > OBJ1, MATERIAL, 112 > OBJ1, HEAT_FLVX, 0.000000E+00, 3.000000E+02 > OBJ2, NAME, APT1 > OBJ2, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ2, SIZE, 0.000000E+00, 2.000000E-01, 1.000000E+00 > OBJ2, CLIPART, cube12 > OBJ2, ROTATION24, 1 > OBJ2, TYPE, APERTURE > OBJ2, PRESSURE, 0.000000E+00 > OBJ2, TEMPERATURE, 1.000000E+01 > OBJ2, COEFFICIENT, 1.000000E+00 > OBJ3, NAME, APT2 > OBJ3, POSITION, 8.000000E-01, 1.000000E+00, 0.000000E+00 > OBJ3, SIZE, 2.000000E-01, 0.000000E+00, 1.000000E+00 > OBJ3, CLIPART, cube12 > OBJ3, ROTATION24, 1 > OBJ3, TYPE, APERTURE > OBJ3, PRESSURE, 0.000000E+00 > OBJ3, TEMPERATURE, 1.000000E+01 > OBJ3, COEFFICIENT, 1.000000E+00 M stop