TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2015, Date 20/04/16
CPVNAM=VDI; SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
Click here for results
ENDDIS
VRV USE
PROBE 5.000000E-01 5.000000E-02 8.500000E-01; PROBE ON
SLICE X
VECTOR ON
VECTOR REFERENCE 2.500000E+00
VARIABLE Velocity; VARIABLE RANGE 0.000000E+00 5.500000E+00
CONTOUR ON
CONTOUR FILL ON
MSG Press OK to animate the solution
PAUSE
ANIMATE
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 112
************************************************************
Group 1. Run Title
TEXT(Falling box with broken tip )
************************************************************
Echo save-block settings for Group 1
save1begin
BOOLEAN(MOVOB)
MOVOB=T
save1end
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.,1.,20)
* Modify regions
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,50,50)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(149)=OBID ;NAME(150)=PRPS
* Solved variables list
SOLVE(P1,V1,W1)
* Stored variables list
STORE(PRPS,OBID)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
ENUT =0.
DRH1DP =1.0E-05
DVO1DT =3.41E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS)=0.
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
No PATCHes used for this Group
EGWF = T
************************************************************
Echo save-block settings for Group 13
SAVE13BEGIN
IF(MOVOB) THEN
define moving objects by In-Form
(MOVOB of BLOCK is OFFSET(0&0&2.1))
(MOVOB of BLOCK is POS(0&0&0&60*tim&0&0))
(MOVOB of TIP is OFFSET(0&2.5&-0.5) with PARENT=BLOCK)
(MOVOB of TIP is POS(0&0&0&-180*tim&0&0))
ENDIF
SAVE13END
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 20
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,0.5 )
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
ISG62 = 0
SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET)
SPEDAT(SET,MOFOR,VELCON,R,100.)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
************************************************************
Echo save-block settings for Group 19
save19begin
IF(MOVOB) THEN
SPEDAT(SET,MOFOR,mofFILE,C,NOTSET)
ELSE
SPEDAT(SET,MOFOR,mofFILE,C,BROKNTIP.mof)
ENDIF
save19end
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
OUTPUT(V1 ,Y,Y,Y,N,Y,Y)
OUTPUT(W1 ,Y,Y,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = 9
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 10000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='M'
GVIEW(P,-1.,0.,0.)
GVIEW(UP,0.,1.,0.)
GVIEW(VDIS,3.427886)
GVIEW(CENTRE,0.5,2.5,2.5)
> DOM, SIZE, 1.000000E+00, 5.000000E+00, 5.000000E+00
> DOM, MONIT, 5.000000E-01, 5.000000E-02, 8.500000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_Y_1, 50, 1.000000E+00
> GRID, RSET_Z_1, 50, 1.000000E+00
> OBJ, NAME, BLOCK
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.600000E+00
> OBJ, SIZE, 1.000000E+00, 2.500000E+00, 5.000000E-01
> OBJ, GEOMETRY, cube
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, DOMAIN
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, INI_PRESS, 0.
> OBJ, NAME, B2
> OBJ, POSITION, 0.000000E+00, 4.800000E+00, 4.800000E+00
> OBJ, SIZE, 1.000000E+00, 1.000004E-01, 1.000004E-01
> OBJ, GEOMETRY, cubet
> OBJ, TYPE, PRESSURE_RELIEF
> OBJ, PRES_RELIEF, 1000. ,0.
> OBJ, NAME, TIP
> OBJ, POSITION, 0.000000E+00, 2.500000E+00, 1.600000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 5.000000E-01
> OBJ, GEOMETRY, cube
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, DOMAIN
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, INI_PRESS, 0.
STOP
DISTIL=T
EX(P1 )= 3.907E+01
EX(V1 )= 1.132E+00
EX(W1 )= 1.255E+00
EX(OBID)= 2.000E-01