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