TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2014, Date 23/07/14 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY A box is filled with two immiscible fluids. The bottom fluid is denser than the upper. At time t=0, the interface between the two liquids is uniformly inclined, and is then let free to oscillate under the infulence of gravity. The vessel is 13m long, and 4m wide. The lower liquid is 0.2m deep, the upper liquid is 0.2m deep. The original inclination of the interface is +/0.035m with respect to the average interface height. The top liquid contains two anode blocks. The upper fluid depth under the anodes is 0.04m ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 498 ************************************************************ Group 1. Run Title TEXT(ESTER demo, Transient one-phase - Case 3) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,400.,100) * Modify regions RSET(T,1,100,1.) ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,5,9,5) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(15)=HI ;NAME(147)=PRPS NAME(148)=NPOR ;NAME(149)=EPOR NAME(150) =CPOR * Solved variables list SOLVE(P1,U1,V1,W1) * Stored variables list STORE(CPOR,EPOR,NPOR,PRPS,HI,SURF) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices U1AD =0. ;V1AD =0. W1AD = GRND ************************************************************ Group 9. Properties * Domain material index is 51 signifying: * etal SETPRPS(1, 51,ESTR) ENUT =0. DRH1DP =1.0E-10 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(HI)=1.001E-10 ;FIINIT(PRPS)=-1. FIINIT(NPOR)=1. ;FIINIT(EPOR)=1. FIINIT(CPOR)=1. No PATCHes used for this Group INIADD = F ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources This patch is attached to object PREF1 PATCH(PREF1, CELL, -1, 1, 0, 0, 0, 0, 1, 100) COVAL(PREF1, P1, 1., 0.) COVAL(PREF1, U1, 0., SAME) COVAL(PREF1, V1, 0., SAME) COVAL(PREF1, W1, 0., SAME) This patch is attached to object PREF2 PATCH(PREF2, CELL, -1, 1, 0, 0, 0, 0, 1, 100) COVAL(PREF2, P1, 1., 0.) COVAL(PREF2, U1, 0., SAME) COVAL(PREF2, V1, 0., SAME) COVAL(PREF2, W1, 0., SAME) EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 60 SARAH =1000. RESREF(P1)=1.0E-10 ;RESREF(U1)=1.0E-10 RESREF(V1)=1.0E-10 ;RESREF(W1)=1.0E-10 RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.9 ) RELAX(U1 ,FALSDT,8. ) RELAX(V1 ,FALSDT,8. ) RELAX(W1 ,FALSDT,8. ) ************************************************************ Group 18. Limits ************************************************************ Group 19. EARTH Calls To GROUND Station NAMGRD =ESTR PARSOL = F ISG62 = 1 CSG8 ='ESTR' CSG10 ='ESTR' SPEDAT(SET,GXMONI,TRANSIENT,L,F) SPEDAT(SET,MATERIAL,51,L,T) SPEDAT(SET,MATERIAL,52,L,T) SPEDAT(SET,MATERIAL,100,L,T) SPEDAT(SET,MATERIAL,151,L,T) SPEDAT(SET,MATERIAL,198,L,T) SPEDAT(SET,MATERIAL,199,L,T) LG( 5) = T IG( 1) = 2 IG( 2) = 4 IG( 3) = 1 IG( 4) = 1 IG( 6) = 5 RG( 1) =2300. RG( 2) =2100. RG( 3) =1. RG( 4) =4.1E+06 RG( 5) =200. RG( 6) =2.33E+04 RG( 7) =-9.81 RG( 11) =1. RG( 12) =100. RG( 13) =1. RG( 14) =1. RG( 16) =-0.035 RG( 17) =0.035 RG( 18) =2.608696E-07 RG( 19) =1.714286E-06 RG( 31) =9.09E+04 RG( 32) =999.869995 RG( 33) =7800. RG( 34) =9.08E+05 RG( 35) =2700. RG( 36) =1.98E+07 RG( 37) =2700. RG( 38) =2.86E+07 RG( 71) =7800. RG( 72) =2.17E+06 RG( 73) =8950. RG( 74) =1.34E+07 RG( 75) =7800. RG( 76) =6.45E+05 RG( 77) =7800. RG( 78) =8.7E+05 ************************************************************ Group 20. Preliminary Printout ************************************************************ Group 21. Print-out of Variables OUTPUT(HI ,Y,N,Y,N,Y,Y) OUTPUT(PRPS,N,N,N,N,N,N) OUTPUT(NPOR,Y,N,N,N,N,N) OUTPUT(EPOR,Y,N,N,N,N,N) OUTPUT(CPOR,Y,N,N,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 2 ;IYMON = 2 ;IZMON = 1 NPRMON = 100000 NPRMNT = 1 TSTSWP = -5 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NTPRIN = 10 ;ISTPRF = 1 ;ISTPRL = 100000 NPLT = 1 ISWPRF = 1 ;ISWPRL = 100000 IPROF = 3 PATCH(PLOT1, PROFIL, 1, 1, 1, 1, 2, 2, 1, 100) PLOT(PLOT1, HI, 0., 0.) PATCH(PLOT2, PROFIL, 1, 1, 4, 4, 2, 2, 1, 100) PLOT(PLOT2, HI, 0., 0.) PATCH(PLOT3, PROFIL, 1, 1, 9, 9, 2, 2, 1, 100) PLOT(PLOT3, HI, 0., 0.) ************************************************************ Group 24. Dumps For Restarts IDISPA = 10 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='c' GVIEW(P,0.57735,-0.57735,0.57735) GVIEW(UP,-0.408248,0.408248,0.816497) GVIEW(DEPTH,3674.429932) GVIEW(NEARPLANE,0.026) GVIEW(VDIS,6.531529) GVIEW(CENTRE,2.,6.5,0.2) > DOM, SIZE, 4.000000E+00, 1.300000E+01, 4.000000E-01 > DOM, MONIT, 1.012500E+00, 1.057142E+00, 5.000001E-02 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02 > GRID, RSET_X_1, 1, 1.000000E+00 > GRID, RSET_X_2, 1, 1.000000E+00 > GRID, RSET_X_3, 1, 1.000000E+00 > GRID, RSET_X_4, 1, 1.000000E+00 > GRID, RSET_X_5, 1, 1.000000E+00 > GRID, RSET_Y_1, 1, 1.000000E+00 > GRID, RSET_Y_2, 7, 1.000000E+00 > GRID, RSET_Y_3, 1, 1.000000E+00 > GRID, RSET_Z_1, 2, 1.000000E+00 > GRID, RSET_Z_2, 2, 1.000000E+00 > GRID, RSET_Z_3, 1, 1.000000E+00 > DOM, T_AMBIENT, 0.000000E+00 > DOM, CARBON-BLOCKS,NO > DOM, ANODES, YES > DOM, SIZE-X, 1.725000E+00, 2.500000E-01, 1.500001E-01 > DOM, SIZE-Y, 1.270000E+01, 0.000000E+00, 1.499991E-01 > DOM, SIZE-Z-1, 2.000000E-01, 1.600000E-01 > DOM, SIZE-Z-2, 2.000000E-01, 4.000001E-02 > DOM, NUM-X, 1, 1, 1 > DOM, NUM-Y, 7, 0, 1 > DOM, NUM-Z, 2, 2, 1, 1 > DOM, NUM-ANODES, 2, 1 > DOM, TAPPING-GAP, NO > DOM, LORENTZ, OFF > DOM, AUTOMESH, ON > DOM, INTERFACE, SLOPE_IN_Y > OBJ, NAME, INIELEC > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E-01 > OBJ, SIZE, 4.000000E+00, 1.300000E+01, 2.000000E-01 > OBJ, GEOMETRY, default > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, USER > OBJ, OPAQUE, 50 > OBJ, MATERIAL, 52,Bath > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, PREF1 > OBJ, POSITION, 1.500001E-01, 3.778571E+00, 0.000000E+00 > OBJ, SIZE, 1.725000E+00, 1.814286E+00, 1.000000E-01 > OBJ, GEOMETRY, default > OBJ, TYPE, CELLTYPE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 49 > OBJ, OPAQUE, 50 > OBJ, PATCHES, PREF1 > OBJ, NAME, PREF2 > OBJ, POSITION, 1.500001E-01, 3.778571E+00, 2.000000E-01 > OBJ, SIZE, 1.725000E+00, 1.814286E+00, 2.000003E-02 > OBJ, GEOMETRY, default > OBJ, TYPE, CELLTYPE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 49 > OBJ, OPAQUE, 50 > OBJ, PATCHES, PREF2 > OBJ, NAME, ANOD1 > OBJ, POSITION, 1.500001E-01, 1.499991E-01, 2.400000E-01 > OBJ, SIZE, 1.725000E+00, 1.270000E+01, 1.600000E-01 > OBJ, GEOMETRY, default > OBJ, TYPE, ANODE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 238 > OBJ, OPAQUE, 50 > OBJ, NAME, ANOD2 > OBJ, POSITION, 2.125000E+00, 1.499991E-01, 2.400000E-01 > OBJ, SIZE, 1.725000E+00, 1.270000E+01, 1.600000E-01 > OBJ, GEOMETRY, default > OBJ, TYPE, ANODE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 238 > OBJ, OPAQUE, 50 > OBJ, NAME, INTERFAC > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E-01 > OBJ, SIZE, 4.000000E+00, 1.300000E+01, 0.000000E+00 > OBJ, GEOMETRY, cube > OBJ, TYPE, NULL > OBJ, COLOR-MODE, DEFAULT2 > OBJ, OPAQUE, 50 > OBJ, NAME, FREESURF > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 4.000000E-01 > OBJ, SIZE, 4.000000E+00, 1.300000E+01, 0.000000E+00 > OBJ, GEOMETRY, cube > OBJ, GRID, N,N,N > OBJ, TYPE, NULL > OBJ, COLOR-MODE, DEFAULT2 > OBJ, OPAQUE, 50 ************************************************************ MENSAV(S,MENU,ESTR,ESTR) STOP