TALK=T;RUN(1,1) 2ph428

  PHOTON USE
  p
  p1
 
 
 
 
  view z
  norm
  msg                                 Oil volume fractions and
  msg                                 velocity vectors
  msg                                         Time step 1
  msg
  msg                                         Press RETURN after
  msg                                         each pause
  gr z m
  con r2 z m fi;10
  pause
  set vec comp;au2 bv2 cw2
  vec z m
  pause
  p
  p2
 
 
  view z
  norm
  msg                                    Time step 2
  gr z m
  con r2 z m fi;10
  pause
  set vec comp;au2 bv2 cw2
  vec z m
  pause
  p
  p3
 
 
  view z
  norm
  msg                                    Time step 3
  gr z m
  con r2 z m fi;10
  pause
  set vec comp;au2 bv2 cw2
  vec z m
  pause
  p
  p4
 
 
  view z
  norm
  msg                                    Time step 4
  gr z m
  con r2 z m fi;10
  pause
  set vec comp;au2 bv2 cw2
  vec z m
  pause
  p
  p5
 
 
  view z
  norm
  msg                                    Time step 5
  gr z m
  con r2 z m fi;10
  pause
  set vec comp;au2 bv2 cw2
  vec z m
  pause
  p
  p6
 
 
 
  view z
  norm
  msg                                    Time step 6
  gr z m
  con r2 z m fi;10
  pause
  set vec comp;au2 bv2 cw2
  vec z m
  pause
  p
  p7
 
 
  view z
  norm
  msg                                    Time step 7
  gr z m
  con r2 z m fi;10
  pause
  set vec comp;au2 bv2 cw2
  vec z m
  msg
  msg Press e to END
  enduse
 
    Group 1. Run Title and Number
TEXT( Oil slick in a gulf with wind: W428
TITLE
  DISPLAY
  2-dimensional (x-y), polar, transient, elliptic simulation
 
  This demonstration concerns the motion of a wind-driven oil slick.
  Geometry of the problem is simplified using polar-grid option.
 
  The two-phase option of PHOENICS is employed.
  ENDDIS
 
    Group 2. Transience
 STEADY=F;LSTEP=7;TLAST=1.0E+4
 TFRAC (  1) =  1.000E-01 ;TFRAC (  2) =  2.000E-01
 TFRAC (  3) =  3.000E-01 ;TFRAC (  4) =  4.000E-01
 TFRAC (  5) =  5.000E-01 ;TFRAC (  6) =  6.000E-01
 TFRAC (  7) =  7.000E-01
    Group 3. X-Direction Grid Spacing
 CARTES=F;NX=20;XULAST=1.58
 grdpwr(x,nx,xulast,1.0)
    Group 4. Y-Direction Grid Spacing
 NY=10;YVLAST=0.2E4;RINNER=1.0E4
 YFRAC(1)=-10.0;YFRAC(2)=0.1
 grdpwr(y,ny,yvlast,1.0)
    Group 5. Z-Direction Grid Spacing
 NZ=2;ZWLAST=1.0
 ZFRAC(1)=-2.0;ZFRAC(2)=0.5
    Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS=F;STORE(VPOR,EPOR,NPOR)
 SOLUTN(P1,Y,Y,Y,N,N,N)
 SOLVE(U1,U2,V1,V2,W1,W2,R1,R2)
 SOLUTN(u1,Y,Y,Y,N,N,N)
 SOLUTN(v1,Y,Y,Y,N,N,N)
 SOLUTN(w1,Y,Y,Y,N,N,N)
 SOLUTN(u2,Y,Y,Y,N,N,N)
 SOLUTN(v2,Y,Y,Y,N,N,N)
 SOLUTN(w2,Y,Y,Y,N,N,N)
    Group 9. Properties
 RHO1=1000.0;RHO2=800.0
    Group 10.Inter-Phase Transfer Processes
 CFIPS=1.0E+3;RLOLIM=0.001
    Group 11.Initialise Var/Porosity Fields
 INIADD=F;FIINIT(R1)=1.0;FIINIT(R2)=0.0
 PATCH(BED,LINVLY,1,NX,1,NY/2,1,NZ,1,lstep)
 INIT(BED,VPOR,0.0001,0.01)
 INIT(BED,EPOR,0.0001,0.01)
 INIT(BED,NPOR,0.0001,0.01)
 PATCH(BED2,LINVLY,1,NX,NY/2+1,NY,1,NZ,1,lstep)
 INIT(BED2,VPOR,-0.0001,0.1)
 INIT(BED2,EPOR,-0.0001,0.1)
 INIT(BED2,NPOR,-0.0001,0.1)
 PATCH(START,INIVAL,NX-1,NX,NY-1,NY,NZ,NZ,1,1)
 INIT(START,R1,0.0,0.99)
 INIT(START,R2,0.0,0.01)
    Group 13. Boundary & Special Sources
 PATCH(WIND,HIGH,1,NX,NY/2,NY,NZ,NZ,1,LSTEP)
 COVAL(WIND,U1,0.001,-5.0)
 PATCH(GRAVITY,PHASEM, 1, 1, 1, 1, 1, 1, 1,lstep)
 COVAL(GRAVITY,W1,FIXFLU,-9.81)
 COVAL(GRAVITY,W2,FIXFLU,-9.81)
 PATCH(REF ,CELL, nx,nx, 1, ny, nz,nz, 1,lstep)
 COVAL(REF ,p1, 1.e10, 0.0)
   COVAL(REF ,V1, 0.0, 0.0)
   COVAL(REF ,R1, 1.0, 0.0)
    Group 15. Terminate Sweeps
 LSWEEP=100
 RESREF(R2)=1.0E-10
    Group 17. Relaxation
 RELAX(U1,FALSDT,1.0E2);RELAX(U2,FALSDT,1.0E2)
 RELAX(V1,FALSDT,1.0E2);RELAX(V2,FALSDT,1.0E2)
 RELAX(W1,FALSDT,0.1);RELAX(W2,FALSDT,0.1)
 RELAX(R1,LINRLX,5.0E-1);RELAX(R2,LINRLX,5.0E-1)
    Group 18. Limits on variables values or increments to them
 VARMIN(R1)=1.E-10;VARMIN(R2)=1.E-10
    Group 19. EARTH Calls To GROUND Station
 PHS2PA=2.0E3
 CSG1=PHI;CSG2=XYZ;IDISPA=1
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    Group 21. Print-out of Variables
 OUTPUT(P1,Y,N,N,Y,Y,Y);OUTPUT(P2,Y,N,N,N,N,N)
 OUTPUT(W1,Y,Y,Y,Y,Y,Y);OUTPUT(W2,Y,Y,Y,Y,Y,Y)
 OUTPUT(R1,N,N,N,N,N,N);OUTPUT(R2,Y,N,Y,Y,Y,Y)
 OUTPUT(VPOR,N,N,N,N,N,N);OUTPUT(EPOR,N,N,N,N,N,N)
 OUTPUT(NPOR,N,N,N,N,N,N)
    Group 22. Monitor Print-Out
 IXMON=NX-1;IYMON=NY-1;IZMON=NZ
 UWATCH=T;HIGHLO=F
    Group 23.Field Print-Out & Plot Control
 NTPRIN=1;NXPRIN=1;NYPRIN=1;IXPRL=10;IYPRL=10
 IPLTF=2;IPLTL=LSWEEP;ITABL=1;NCOLPF=50
 PATCH(MAP,CONTUR,1,NX,1,NY,NZ,NZ,1,LSTEP)
 PLOT(MAP,R2,0.0, 1.0E+01)
    Group 24. Dumps For Restarts
SELREF=T; RESFAC=1.E-3
TSTSWP=-1
NAME(U2)=AU2;NAME(V2)=BV2;NAME(W2)=CW2