DISPLAY

  Library case 653:

  Tutorial:

  Water flow into a bund (SEM)  

  ENDDIS
  
 TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2018, Date 17/01/18
 CPVNAM=VDI; SPPNAM=Core
 ************************************************************
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 1. Run Title
 TEXT(Water pouring into a bund)
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.,0.5,500)
    * Modify regions
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,28,12,17)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(144)=PRPS ;NAME(145)=VFOL
 NAME(146)=SURN ;NAME(147)=EPKE
 NAME(148)=DEN1 ;NAME(149)=EL1
 NAME(150) =ENUT
    * Solved variables list
 SOLVE(P1,U1,V1,W1,VFOL,SURN)
    * Stored variables list
 STORE(ENUT,EL1,DEN1,EPKE,PRPS)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 TURMOD(KECHEN)
 
 ************************************************************
  Group 8. Terms & Devices
 TERMS (VFOL,N,N,N,N,Y,N)
 TERMS (SURN,N,N,N,N,Y,N)
 GALA = T
 RUPLIM =0.999 ;RLOLIM =1.0E-03
 ************************************************************
  Group 9. Properties
 PRESS0 =1.01325E+05 ;TEMP0 =273.
 RHO1 =998.22998
 ENUL =1.006E-06
 CP1 =4181.799805
 DVO1DT =1.18E-04
 PRT(KE)=0.75 ;PRT(EP)=1.15
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1)=0. ;FIINIT(PRPS)=-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
 
 PATCH(BUOYANCY, PHASEM, 0, 0, 0, 0, 0, 0, 1, 500)
 COVAL(BUOYANCY, U1, FIXFLU, GRND1)
 COVAL(BUOYANCY, V1, FIXFLU, GRND1)
 COVAL(BUOYANCY, W1, FIXFLU, GRND1)
 
 BUOYA =0. ; BUOYB =0.
 BUOYC =-9.81
 EGWF = T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 20
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=200 ;LITER(SURN)=1
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,1. )
 RELAX(KE ,LINRLX,0.5 )
 RELAX(EP ,LINRLX,0.5 )
 KELIN = 3
 ************************************************************
  Group 18. Limits
 VARMAX(SURN)=1. ;VARMIN(SURN)=0.
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK = T
 SURF = T
 PARSOL = F
 CONWIZ = T
 IENUTA = 2
 IPRPSA = 67
 ISURFA = 1
 ISG62 = 1
 SPEDAT(SET,OUTPUT,NOFIELD,L,T)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(VFOL,Y,Y,Y,N,Y,Y)
 OUTPUT(SURN,Y,Y,Y,N,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 22 ;IYMON = 6 ;IZMON = 2
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 ISWPRF = 1 ;ISWPRL = 100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 IDISPA = 10 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='N'
 
 GVIEW(P,-1.060937E-03,-0.999779,0.020971)
 GVIEW(UP,0.023874,0.02094,0.999495)
 GVIEW(VDIS,1.146957)
 GVIEW(CENTRE,0.5,0.5,0.15)
 
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 3.000000E-01
> DOM,    MONIT,       7.555072E-01, 4.642857E-01, 3.190332E-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,    -14, 1.100000E+00,G
> GRID,   RSET_X_2,     -3, 1.100000E+00,G
> GRID,   RSET_X_3,    -11, 1.100000E+00,G
> GRID,   RSET_Y_1,      2,-1.100000E+00,G
> GRID,   RSET_Y_2,     -4, 1.100000E+00,G
> GRID,   RSET_Y_3,     -4, 1.100000E+00,G
> GRID,   RSET_Y_4,      2, 1.100000E+00,G
> GRID,   RSET_Z_1,      3,-1.100000E+00,G
> GRID,   RSET_Z_2,     -2, 1.100000E+00,G
> GRID,   RSET_Z_3,    -12, 1.100000E+00,G
> DOM,    INI_AMB,    YES
> DOM,    INI_BUOY,   YES
 
> OBJ,    NAME,        BEAM
> OBJ,    POSITION,    5.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E-01, TO_END,       6.000000E-02
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
 
> OBJ,    NAME,        WATER-IN
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 5.000000E-01, 9.000000E-02
> OBJ,    GEOMETRY,    cube3t
> OBJ,    TYPE,        INLET
> OBJ,    FLUID,       HEAVY
> OBJ,    VELOCITY,    2. ,0. ,0.
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    TURB-INTENS, 5.
 
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    AT_END,       2.000000E-01, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 6.000000E-01, 9.000000E-02
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube12t
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,     P_AMBIENT
> OBJ,    COEFFICIENT, 1000.
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    TURBULENCE,  SAME , SAME
 
> OBJ,    NAME,        FREE-TOP
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, AT_END
> OBJ,    SIZE,        TO_END,       TO_END,       0.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube12t
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,     P_AMBIENT
> OBJ,    COEFFICIENT, 1000.
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    TURBULENCE,  SAME , SAME
STOP