PHOTON USE p n2 up z gr x 1 msg contours of den1 at time = 0.05s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n4 up z gr x 1 msg contours of den1 at time = 0.10s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n6 up z gr x 1 msg contours of den1 at time = 0.15s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n8 up z gr x 1 msg contours of den1 at time = 0.20s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n10 up z gr x 1 msg contours of den1 at time = 0.25s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n12 up z gr x 1 msg contours of den1 at time = 0.30s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n14 up z gr x 1 msg contours of den1 at time = 0.35s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n16 up z gr x 1 msg contours of den1 at time = 0.40s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n18 up z gr x 1 msg contours of den1 at time = 0.45s con den1 x 1 fi;.01 msg velocity vectors vec x 1 pause p n20 up z gr x 1 msg contours of den1 at time = 0.50s con den1 x 1 fi;.01 msg velocity vectors vec x 1 msg - msg Press e to END enduse GROUP 1. Run title and other preliminaries TEXT(SLUMPING OF A LIQUID COLUMN BY SEM:P106 TITLE DISPLAY FREE SURFACE EXAMPLES - Slumping of a column of liquid Scalar Equation Method (SEM) 2-dimensional (y-z), Cartesian, transient, elliptic simulation This case simulates the slumping of a liquid column that may arise,for example, when container walls break down. In addition an obstacle can be errected in the path of the liquid. liquid | column |------ z ^ | - - | | | - | +-----> | -- | | wall y | -- | | |_____|______|___ ///////////////// enddis GROUP 2. Transience; time-step specification STEADY=F;LSTEP=30;LSTEP=20;TFRAC(1)=-LSTEP;TFRAC(2)= 0.025 GROUP 4. Y-direction grid specification GRDPWR(Y,20,2.0,1.0);NY=10 GROUP 5. Z-direction grid specification GRDPWR(Z,16,0.5,1.0);NZ=10 GROUP 7. Variables stored, solved & named STORE(DEN1,PRPS);SOLVE(VFOL,SURN) SOLUTN(P1,Y,Y,Y,N,N,N);SOLUTN(V1,Y,Y,N,N,N,N) SOLUTN(W1,Y,Y,N,N,N,N) GROUP 8. Terms (in differential equations) & devices ** activate the "gas-and-liquid algorithm", ie volumetric continuity equation, and allow convection fluxes to be modified in GROUND GALA=T;TERMS(VFOL,N,N,N,N,P,P);TERMS(SURN,N,N,N,N,P,P) GROUP 9. Properties of the medium (or media) ** signal that density is to be computed by the HOL method and set the densities of the liquid and gas respectively GROUP 11. Initialization of variable or porosity fields FIINIT(P1)=0.0;FIINIT(V1)=0.0;FIINIT(W1)=0.0;FIINIT(SURN)=0.0 FIINIT(DEN1)=1.189;INIADD=F ** place the initial "pile" of liquid of which the slumping is to be simulated PATCH(LIQUID,INIVAL,1,NX,1,4,1,NZ-4,1,1) redefine patch for reduced grid PATCH(LIQUID,INIVAL,1,NX,1,NY/2,1,NZ/2,1,1) COVAL(LIQUID,SURN,ZERO,1.0);COVAL(LIQUID,DEN1,ZERO,998) ** Activate the next two lines to insert a containment wall PATCH(WALL,NORTH,1,1,3*NY/4,3*NY/4,1,NZ/4,1,LSTEP) COVAL(WALL,V1,FIXVAL,0.0) GROUP 13. Boundary conditions and special sources ** the pressure is held to zero along the open top boundary PATCH(REFP,CELL,1,NX,1,NY,NZ,NZ,1,LSTEP) COVAL(REFP,P1,FIXP,ZERO) ** provide for the gravity-force source of w1 PATCH(GRAV,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP) COVAL(GRAV,W1,FIXFLU,-9.81) ** time patch for transient term PATCH(TCON,CELL,1,NX,1,NY,1,NZ,1,LSTEP) COVAL(TCON,SURN,GRND,GRND) GROUP 15. Termination of sweeps LSWEEP=10;LITER(SURN)=1 GROUP 16. Termination of iterations RESREF(P1)=1.E-6;RESREF(V1)=1.E-6;RESREF(W1)=1.E-6 GROUP 17. Under-relaxation devices RELAX(V1,FALSDT,0.01);RELAX(W1,FALSDT,0.01) GROUP 19. Data communicated by satellite to GROUND ** provide for the dumping of field data at each time step, for us by PHOTON IPRPSA=67;IPRPSB=0;IDISPA=2;CSG1=N;SURF=T;RLOLIM=0.4;RUPLIM=0.6 VARMIN(SURN)=0.0; VARMAX(SURN)=1.0 SPEDAT(SET,GXMONI,TRANSIENT,L,F) GROUP 22. Spot-value print-out TSTSWP=-1;IYMON=NY/2;IZMON=NZ/2 GROUP 23. Field print-out and plot control NTPRIN=LSTEP/2;OUTPUT(P1,Y,Y,Y,Y,Y,Y);OUTPUT(V1,Y,Y,Y,Y,Y,Y) OUTPUT(W1,Y,Y,Y,Y,Y,Y);OUTPUT(SURN,Y,N,y,N,N,N) OUTPUT(DEN1,Y,N,y,N,N,N) ** provide contour-plot output on the line printer PATCH(MAP,CONTUR,1,1,1,NY,1,NZ,1,LSTEP) PLOT(MAP,SURN,ZERO,1)