GROUP 1. Run title and other preliminaries
TEXT(1D Y; Test of u1*r = constant
TITLE
DISPLAY
The exact solution for this case is:
u1 * r = 0.01 for all y
regardless of the value of RINNER.
The numerical solution shows a maximum deviation of
6.5 % close to y=0.0, for rinner = 0.0 , but less than
1.0 % for all y, when rinner = 0.1 * yvlast.
ENDDIS
GROUP 3. X-direction grid specification
CARTES=F;XULAST=0.001
GROUP 4. Y-direction grid specification
GRDPWR(Y,20,1.,1.0)
RINNER=0.1*YVLAST
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1)
GROUP 9. Properties of the medium (or media)
ENUL=100.0
GROUP 11. Initialization of variable or porosity fields
FIINIT(U1)=0.01;IURINI=1
GROUP 13. Boundary conditions and special sources
PATCH(INNER,CELL,1,NX,1,1,1,1,1,1)
IURVAL=1 ! i.e interpret the 'val' of coval as u1*r, not u1
COVAL(INNER,U1,FIXVAL,0.01)
PATCH(OUTER,NORTH,1,NX,NY,NY,1,1,1,1)
COVAL(OUTER,P1,FIXP,0.0);COVAL(OUTER,U1,FIXVAL,0.01)
GROUP 15. Termination of sweeps
LSWEEP=100
RESREF(P1)=1.E-7;RESREF(U1)=1.E-7;RESREF(V1)=1.E-7
GROUP 17. Under-relaxation devices
LSG57=T ! i.e use PHOENICS-recommended settings
GROUP 22. Spot-value print-out
IXMON=-1;IYMON=6; TSTSWP=-1
GROUP 23. Field print-out and plot control
IURPRN=1 ! i.e. print u1*r in place of u1
ITABL=2
NYPRIN=1
PATCH(LONGPLOT,PROFIL,1,1,1,NY,1,1,1,1)
PLOT(LONGPLOT,P1,0.0,0.0);PLOT(LONGPLOT,U1,0.0,0.0)
GROUP 24. Dumps for restarts
isg40=1