```

GROUP 1. Run title and other preliminaries

TEXT(Powre-Law_2D Fluid Ellip. Lam Pipe
TITLE

The problem concerns the steady laminar flow of a power-law
pseudo-plastic non-Newtonian fluid in a pipe. For fully-developed
flow the analytical solution for the pressure drop is given by:

dp/dz = 4.*rho*win**2*[(2+6n)/n]**n/Re/D

where D is the pipe diameter, win the mean velocity and Re the
power-law Reynolds number, defined by:

Re = D**n*win**(2-n)*rho/K

The problem is solved by use of the 2d y-z elliptic solver for a
Reynolds number of 10 and a power-law index of 0.5. For these
conditions in fully developed flow, the expected pressure drop
and centre-line axial velocity are about 6.33 N/m**3 and 1.66m/s,
respectively. The PHOENICS predictions show fairly good
agreement with these results.

REAL(RIN,DIN,WIN,AIN,DPDZ);RIN=0.1;DIN=2.*RIN;WIN=1.0
GROUP 2. Transience; time-step specification
CARTES=F;XULAST=0.1;AIN=XULAST*RIN*RIN/2.
GROUP 4. Y-direction grid specification
NY=12;GRDPWR(Y,NY,RIN,-1.5)
GROUP 5. Z-direction grid specification
NZ=15;GRDPWR(Z,NZ,10.*RIN,1.9)
GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);STORE(VISL,GEN1)
SOLUTN(P1,Y,Y,Y,P,P,P)
GROUP 9. Properties of the medium (or media)
RHO1=1.0;ENUT=0.
REAL(REY,POWER);REY=10.;POWER=0.5
ENULA=WIN**(2.0-POWER)*DIN**POWER/REY
ENULB=POWER
** enul = enula*(dwdy)*[0.5*(enulb-1)]
ENUL=STRAIN;GENK=T
DPDZ=4.*RHO1*WIN**2/DIN/REY*((2.+6.*POWER)/POWER)**POWER
DPDZ
REY
POWER
GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=WIN
GROUP 13. Boundary conditions and special sources

PATCH(INLET,LOW,1,NX,1,NY,1,1,1,LSTEP)
COVAL(INLET,P1,FIXFLU,WIN*RHO1);COVAL(INLET,W1,ONLYMS,WIN)

PATCH(OUTL,HIGH,1,NX,1,NY,NZ,NZ,1,LSTEP)
COVAL(OUTL,P1,1.E2,0.0);COVAL(OUTL,W1,ONLYMS,SAME)

PATCH(WALL,NWALL,1,NX,NY,NY,1,NZ,1,LSTEP)
COVAL(WALL,W1,1.0,0.0)

GROUP 15. Termination of sweeps
LSWEEP=50
GROUP 16. Termination of iterations
RESREF(P1)=1.E-12*WIN*AIN
RESREF(W1)=RESREF(P1)*RHO1*WIN;RESREF(V1)=RESREF(W1)
GROUP 17. Under-relaxation devices
REAL(DTF);DTF=ZWLAST/WIN/NZ
RELAX(P1,LINRLX,1.0)
RELAX(W1,FALSDT,DTF);RELAX(V1,FALSDT,DTF)
RELAX(VISL,LINRLX,1.0)
GROUP 22. Spot-value print-out
IYMON=NY/2;IZMON=NZ/2;TSTSWP=-1
GROUP 23. Field print-out and plot control
NPLT=5;NYPRIN=1;NZPRIN=1
GROUP 24. Dumps for restarts
ENULA
ENULB
```