PHOTON USE p gr ou Z 1 SET VEC REF 0.005;ve Z 1 sh msg Velocity distribution: msg Pressto continue pause;clear gr ou z 1 msg Temperature distribution : con h1 z 1 fil;0.0004 msg Press e to END ENDUSE TEXT(NATURAL CONVECTION IN POROUS MEDIA:131 DISPLAY This case shows how the PLANT option should be used to introduce the non-linear momentum sinks of porous resistive media governed by Ergun's drag law. The problem is similar to one for empty partitioned cavern published by S.V.Zhubrin and A.P.Khrupov in PHOENICS Journal of CFD and its applications, vol.7, N1, pp. 93-106, 1994 ENDDIS Special data ============ REAL(BETA,TREF,RHOREF) RHOREF= 0.7;BETA=1./(273.+230.);TREF=503. * 'Phasem' force = (rho - rhoref) * grav BUOYD=RHOREF BUOYB=-9.81 REAL(POROS,DIAM,LENGTH);POROS=0.5;DIAM=0.02 LENGTH=DIAM*POROS/(1.-POROS) * GROUP 3. X-direction grid specification. NX=30 GRDPWR(X,NX,0.1,1.) INTEGER(NXF01,NXL01); NXF01=1; NXL01=8 INTEGER(NXF02,NXL02); NXF02=9; NXL02=13 INTEGER(NXF03,NXL03); NXF03=14; NXL03=17 INTEGER(NXF04,NXL04); NXF04=18; NXL04=22 INTEGER(NXF05,NXL05); NXF05=23; NXL05=30 * GROUP 4. Y-direction grid specification. NY=30 GRDPWR(Y,NY,0.1,1.) INTEGER(NYF01,NYL01); NYF01=1; NYL01=8 INTEGER(NYF02,NYL02); NYF02=9; NYL02=22 INTEGER(NYF03,NYL03); NYF03=23; NYL03=30 * GROUP 7. Variables stored, solved & named. SOLVE(P1);SOLUTN(P1,Y,Y,Y,N,N,N) SOLVE(U1,V1,H1) STORE(RHO1,EPOR,VPOR,NPOR,PRPS) * GROUP 8. Terms (in differential equations) & devices. TERMS(H1,N,Y,Y,N,Y,N) * GROUP 9. Properties of the medium (or media). ENUL= 3.72E-05;PRNDTL(H1)=0.678 * 'Phasem' force = (rho - rhoref) * grav RHO1=RECSCAL;RHO1B=2.78E-3 * GROUP 11. Initialization of variable or porosity fields. FIINIT(H1)=TREF;FIINIT(P1)=1.e-5;FIINIT(RHO1)=RHOREF FIINIT(EPOR)=POROS;FIINIT(NPOR)=POROS;FIINIT(VPOR)=POROS CONPOR(0.0,CELL,-NXF03,-NXL03,1,-NYL01,1,1) CONPOR(0.0,CELL,-NXF03,-NXL03,-NYF03,NY,1,1) * GROUP 13. Boundary conditions and special sources. PATCH(REFP,CELL,NX,NX,NY,NY,1,1,1,1) COVAL(REFP,P1,FIXP,0.0) COVAL(REFP,U1,ONLYMS,0.0) COVAL(REFP,V1,ONLYMS,0.0) COVAL(REFP,H1,ONLYMS,SAME) * WALL Boundary Condition, Named NORTH WALL(NOR,NORTH,1,NXL02,NY,NY,1,1,1,1) WALL(NOR1,NORTH,NXF04,NX,NY,NY,1,1,1,1) * WALL Boundary Condition, Named SOUTH WALL(SOU,SOUTH,1,NXL02,1,1,1,1,1,1) WALL(SOU1,SOUTH,NXF04,NX,1,1,1,1,1,1) * WALL Boundary Condition, Named WEST WALL(WES,WEST,1,1,1,NY,1,1,1,1) COVAL(WES,H1,1./PRNDTL(H1),508.0) * WALL Boundary Condition, Named EAST WALL(EAS,EAST,NX,NX,1,NY,1,1,1,1) COVAL(EAS,H1,1./PRNDTL(H1),498.0) * GRAVITY FORCE Boundary Condition PATCH(BUOY,PHASEM,1,NX,1,NY,1,1,1,1) COVAL(BUOY,V1,FIXFLU,DENSDIFF) * Ergun's resitance law PATCH(ERG1,PHASEM,1,NX,1,NY,1,NZ,1,1000) COVAL(ERG1,U1,150.*ENUL/LENGTH**2,0.0) PLANTBEGIN PATCH(ERG2,PHASEM,1,NX,1,NY,1,NZ,1,1000) CO=RG(2)*ABS(U1) COVAL(ERG2,U1,GRND,0.0) PATCH(ERG4,PHASEM,1,NX,1,NY,1,NZ,1,1000) CO=RG(2)*ABS(V1) COVAL(ERG4,V1,GRND,0.0) PLANTEND PATCH(ERG3,PHASEM,1,NX,1,NY,1,NZ,1,1000) COVAL(ERG3,V1,150.*ENUL/LENGTH**2,0.0) * GROUP 15. Termination of sweeps. LSWEEP=100 * GROUP 16. Termination of iterations. LITER(P1)=30;ENDIT(P1)=5.0E-6 LITER(H1)=1;ENDIT(H1)=5.E-6 ISOLX=0;ISOLY=0 * GROUP 17. Under-relaxation devices. RELAX(U1,FALSDT,1.0E-2) RELAX(V1,FALSDT,1.0E-2) * GROUP 18. Limits on variables or increments to them. VARMIN(H1)=490.0;VARMAX(H1)=510.0 VARMIN(U1)=-100.0;VARMAX(U1)=100.0 VARMIN(V1)=-100.0;VARMAX(V1)=100.0 * GROUP 19. Data communicated by satellite to GROUND. NAMSAT=MOSG RG(1)=150.*ENUL/LENGTH**2 RG(2)=1.75/LENGTH * GROUP 20. Preliminary print-out. * GROUP 21. Print-out of variables. OUTPUT(P1,Y,N,N,y,Y,Y);OUTPUT(U1,Y,N,N,y,Y,Y) OUTPUT(V1,Y,N,N,y,Y,Y);OUTPUT(H1,Y,N,N,y,Y,Y) * GROUP 22. Spot-value print-out. IXMON=9;IYMON=20 * GROUP 23. Field print-out and plot control. ITABL=3 tstswp=-1 dmpstk=t DISTIL=T EX(P1)=8.714E-03; EX(U1)=7.371E-04; EX(V1)=7.447E-04 EX(H1)=5.021E+02; EX(PRPS)=1.501E+01; EX(NPOR)=5.000E-01 EX(VPOR)=4.644E-01; EX(EPOR)=5.000E-01; EX(RHO1)=6.643E-01 LIBREF=131