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.04m deep. The original
inclination of the interface is +/0.035m with respect to the
average interface height.
The expected (and calculated) oscillation period is 148.5s.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 498
************************************************************
Group 1. Run Title
TEXT(ESTER demo, Transient one-phase - Case 1)
************************************************************
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,1,9,4)
************************************************************
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,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, 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, 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(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(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,198,L,T)
SPEDAT(SET,MATERIAL,199,L,T)
SPEDAT(SET,MATERIAL,151,L,T)
LG( 5) = T
IG( 1) = 2
IG( 2) = 4
IG( 3) = 1
IG( 4) = 1
IG( 6) = 4
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(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 = 1 ;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 ='b'
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.529723)
GVIEW(CENTRE,2.,6.5,0.12)
> DOM, SIZE, 4.000000E+00, 1.300000E+01, 2.400000E-01
> DOM, MONIT, 2.000000E+00, 2.166667E+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_Y_1, 9, 1.000000E+00
> GRID, RSET_Z_1, 2, 1.000000E+00
> GRID, RSET_Z_2, 2, 1.000000E+00
> DOM, T_AMBIENT, 0.000000E+00
> DOM, CARBON-BLOCKS,NO
> DOM, ANODES, NO
> DOM, SIZE-X, 7.700000E-01, 4.000000E-02, 1.600000E-01
> DOM, SIZE-Y, 1.400000E+00, 1.800000E-01, 1.000000E-01
> DOM, SIZE-Z-1, 2.000000E-01, 0.000000E+00
> DOM, SIZE-Z-2, 3.999998E-02, 3.999998E-02
> DOM, NUM-X, 1, 1, 1
> DOM, NUM-Y, 1, 1, 1
> DOM, NUM-Z, 2, 2, 1, 1
> DOM, NUM-ANODES, 0, 0
> 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, 3.999998E-02
> 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, 0.000000E+00, 4.333333E+00, 0.000000E+00
> OBJ, SIZE, 4.000000E+00, 1.444444E+00, 1.000000E-01
> OBJ, GEOMETRY, default
> OBJ, TYPE, CELLTYPE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 49
> OBJ, OPAQUE, 50
> OBJ, PATCHES, SKIP, PREF1
> OBJ, NAME, PREF2
> OBJ, POSITION, 0.000000E+00, 4.333333E+00, 2.000000E-01
> OBJ, SIZE, 4.000000E+00, 1.444444E+00, 2.000000E-02
> OBJ, GEOMETRY, default
> OBJ, TYPE, CELLTYPE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 49
> OBJ, OPAQUE, 50
> OBJ, PATCHES, PREF2
> 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, 2.400000E-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