PHOTON USE
  p
  phi
 
 
 
  msg grid
  gr x 1
  msg Press  to continue
  pause
  cl
  msg 1st-phase (light fluid) velocity vectors
  set vec comp
  - vcrt wcrt
  vec x 1 sh
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg 2nd-phase (dense fluid) velocity vectors
  set vec comp
  - v2cr w2cr
  vec x 1 sh
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of pressure
  con p1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 1st-phase (light fluid) volume fraction
  con r1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 2nd-phase (dense fluid) volume fraction
  con r2 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 2nd-phase (dense fluid) viscosity
  con enul x 1 fi;0.002
  gr ou x 1
  msg Press e to END
  enduse
TEXT(2PHS NON NEWTONIAN FLOW IN DUCT :   W569
TITLE
  DISPLAY
    This case concerns plane, two-dimensional, incompressible
    two phase Non Newtonian flow through a 180 degree turnaround
    duct. The case is based on Library case b527
  ENDDIS
    GROUP 1. Run title
REAL(REYNO,WIN,DIAM);REYNO=3.E4;DIAM=30.;WIN=1.0
    GROUP 6. Body-fitted coordinates or grid distortion
BFC=T;NONORT=T
GSET(D,1,10,20,1.0,DIAM,100.0)
GSET(C,K21,F,K21,+,0.0,110.0,0.0)
GSET(C,K17,F,K21,+,0.0,0.0,-100.0,INC,-1.5)
GSET(C,K5,F,K17,RX,-3.14159,95.0,0.0,INC,1.0)
GSET(C,K1,F,K5,+,0.0,0.0,100.0,INC,1.5)
    Set wup=t to account better for the high curvature of
    the w resolute...
WUP=T
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);SOLUTN(P1,Y,Y,Y,N,N,N)
    Activate second phase variables
ONEPHS=F;SOLVE(V2,W2,R1,R2)
    Reynolds number and Coefficient of Drag for the spheres may
    be STOREd using STORE(REYN) and STORE(CD) respectively,
    activate the following line to do so.
    STORE(VREL,APRJ,CFIPS)
    GROUP 9. Properties of the medium (or media)
ENUL=WIN*DIAM/REYNO
    Switch to power-law viscosity fluid
ENULA=ENUL;ENULB=2.0;ENUL=STRAIN;STORE(ENUL,GEN1);ENUT=0.0
    Set phase densities
RHO1=1;RHO2=1000.0
    GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=1.E-10;FIINIT(W1)=WIN
FIINIT(W2)=WIN;FIINIT(R1)=0.5;FIINIT(R2)=0.5
    Activate drag law for spheres - diameter is CFIPB
CFIPS=GRND7;CFIPA=0.001;CFIPB=5.0E-3
    GROUP 13. Boundary conditions and special sources
    Inlet - 50/50 by volume of each phase
INLET(INLET,LOW,1,1,1,NY,1,1,1,1)
VALUE(INLET,P1,0.5*RHO1*WIN);VALUE(INLET,W1,WIN)
VALUE(INLET,P2,0.5*RHO2*WIN);VALUE(INLET,W2,WIN)
    Outlet - both phases allowed to leave
PATCH(OUTLET,HIGH,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,1.E-3,0.0);COVAL(OUTLET,P2,1.E-3*RHO2,0.0)
COVAL(OUTLET,V1,ONLYMS,0.0);COVAL(OUTLET,W1,ONLYMS,0.0)
COVAL(OUTLET,V2,ONLYMS,0.0);COVAL(OUTLET,W2,ONLYMS,0.0)
    Friction on internal surfaces - Phase 1 only by default
WALL(NORTH,NORTH,1,NX,NY,NY,1,NZ,1,1)
WALL(SOUTH,SOUTH,1,NX,1,1,1,NZ,1,1)
    GROUP 15. Termination of sweeps
LSWEEP=500
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.5)
RELAX(V1,FALSDT,1.0);RELAX(W1,FALSDT,1.0)
RELAX(V2,FALSDT,1.0);RELAX(W2,FALSDT,1.0)
RELAX(R1,LINRLX,.5);RELAX(R2,LINRLX,.5)
RELAX(ENUL,LINRLX,.1)
    GROUP 22. Spot-value print-out
IYMON=7;IZMON=15;TSTSWP=-1