DISPLAY This case is created for demonstration using of (modsor In-Form statement. The case considers the flow near the thin pipe grid which is placed horizontally at north faces of cells for iy=iysh. The diameter of pipes is very small and equal diam=0.01 m. Thus the pipe grid is a subgrid object. The pressure drop at pipe grid can be calculated by formulas dP = f*2*rho*vel^2, [N/m^2] f=0.33*Re^(-1/5) and Re=abs(vel)*diam/enul The resistance at the pipe grid is simulated by updating the built-in source term of equation for v1 variable at iy=iysh cells. The inform variable is introduced for comparison two way of updating by in-form and by group 12 feature. ENDDIS PHOTON USE p gr ou x 1 gr x 1 y 10 11 z 1 20 col 4 VEC X 1 SH msg WATER RECIRCULATION NEAR THIN PIPE GRID msg Velocity vectors pause vec off;red con vabs x 1 fil;0.0001 gr x 1 y 10 11 z 1 20 col 0 msg Absolute velocity contours pause con off;red con p1 x 1 fil;0.0001 gr x 1 y 10 11 z 1 20 col 0 msg Pressure contours msg Press e to END enduse GROUP 1. Run title and other preliminaries TEXT(Modify built-in source for V1 TITLE ** define uniform grid initially GROUP 4. Y-direction grid specification NY=20;NZ=20 GRDPWR(Y,NY,.5,1.0) INTEGER(IYSH);IYSH=NY/2 ! y-location pipe grid GROUP 5. Z-direction grid specification GRDPWR(Z,NZ,1.,1.0) GROUP 7. Variables stored, solved & named SOLVE(P1,V1,W1);SOLUTN(P1,Y,Y,Y,P,P,P) STORE(VABS,VLSQ,PRPS) GROUP 9. Properties of the medium (or media) SETPRPS(1,67) ! water ENUT=.001 GROUP 13. Boundary conditions and special sources REAL(VELIN);VELIN=10. ! inlet velocity ** inflow boundary PATCH(INL,LOW,1,NX,1,IYSH,1,1,1,1) COVAL(INL,P1,FIXFLU,RHO1*VELIN) COVAL(INL,W1,ONLYMS,VELIN) ** outflow boundary PATCH(OUTL,HIGH,1,NX,IYSH+1,NY,NZ,NZ,1,1) COVAL(OUTL,P1,1.E3,0.0) REAL(DIAM); DIAM=0.01 ! pipe diameter ** the Reynolds number calculations (STORE var REYN is VABS*DIAM/ENUL) ! calculate Reynolds number ** The resistance at the pipe grid ! Modify build-in source of v1 variable BOOLEAN(INFORM); INFORM=T ! by means In-Form INFORM=F ! by means froup 12 feature INFORM=T IF(INFORM) THEN PATCH(MSOR1,CELL,1,NX,IYSH,IYSH,1,NZ,1,1) (MODSOR of V1 at MSOR1 is ANORTH*.33*REYN^(-0.2)*2*RHO1*VLSQ) ELSE PATCH(GP12SOR2,CELL,1,NX,IYSH,IYSH,1,NZ,1,1) COVAL(GP12SOR2,V1,3.0E+2,0.0) ENDIF GROUP 22. Monitor print-out IZMON=NZ/2;IYMON=IYSH;YPLS=T GROUP 24. DUMPS FOR RESTARTS LSWEEP=350;TSTSWP=-1;NPRINT=LSWEEP;NYPRIN=1;NZPRIN=1 LSG57=T ! It is important for inform=t DISTIL=T EX(P1 )=9.922E+03; EX(V1 )=1.723E+00; EX(W1 )=5.832E+00 EX(REYN)=6.303E+04; EX(VLSQ)=4.886E+01; EX(VABS)=6.341E+00 lsg57 STOP