TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2007, Date 23/11/07
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
This case describes how the 'moving grid' method (see MOFOR
entry in Encyclopedia, POLIS) is applied to the simulation
of a moving sphere with a constant accelaration.
The velocity used here is a linear function of time
U(t)= t
and the accelaration is
a = 1.0
The method is implemented by In-Form formulea written in
the Q1 file.
The above velocity is used to provide boundary conditions
for the mass
flow rate and the velocity at INLET; and the accelation
is used to set
the domain body forces to ensure that the domain moves
with the sphere.
This file contains a marco of commands which cause the
Viewer when the Macro button is pressed to display the
animation automatically.
ENDDIS
VRV USE
DOMAIN ON
* Setting object visibility and painting status
OBJECT SHOW TYPE BLOCKAGE
OBJECT PAINT TYPE BLOCKAGE OFF
OBJECT WIREFRAME TYPE BLOCKAGE OFF
OBJECT SHOW TYPE OUTLET
OBJECT PAINT TYPE OUTLET OFF
OBJECT WIREFRAME TYPE OUTLET OFF
VARIABLE Velocity; CON ON
VECTOR ON
ANIMATE 1 20 1
PAUSE
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 805
************************************************************
Group 1. Run Title
TEXT(Moving sphere with constant accelaration)
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.000000E+00,1.000000E+01,10)
* Modify regions
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,18,18,31)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(149) =VLSQ ; NAME(150) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 )
* Stored variables list
STORE(PRPS,VLSQ)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(V1 ,Y,Y,Y,N,N,Y)
SOLUTN(W1 ,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
RHO1 = 1.000000E+00
ENUL = 1.000000E-05
CP1 = 1.005000E+03
ENUT = 0.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = -1.000000E+00
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 InForm settings for Group 13
Inform13BEGIN
patch(in,low,1,nx,1,ny,1,1,1,lstep)
(source of p1 at in is tim*rho1)
(source of w1 at in is tim with onlyms)
patch(acel,phasem,1,nx,1,ny,1,nz,1,lstep)
(source of w1 at acel is 1.0)
Inform13END
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 50
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (P1 ) = 250
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
NAMSAT =MOSG
ASAP = T
PARSOL = T
CALFOR = T
IDISPB = 1 ;IDISPC = 100
RG( 1) = 2.000000E+00
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
EX (P1 ) = 1.298000E-01 ;EX (U1 ) = 0.000000E+00
EX (V1 ) = 3.335000E-02 ;EX (W1 ) = 1.012000E+00
EX (VLSQ) = 0.000000E+00 ;EX (PRPS) = 3.354000E+00
************************************************************
Group 21. Print-out of Variables
OUTPUT(VLSQ,N,N,Y,N,N,N)
OUTPUT(PRPS,N,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 10 ;IYMON = 10 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NTPRIN = 2 ;ISTPRF = 1 ;ISTPRL = 100000
YZPR = T
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 1 ;IDISPB = 1 ;IDISPC = 100
CSG1 ='M'
DISTIL=T
EX(P1 )= 9.743E+00;EX(U1 )= 2.959E-01
EX(V1 )= 2.959E-01;EX(W1 )= 1.007E+01
GVIEW(P,-1.000000E+00,0.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 4.000000E+00, 4.000000E+00, 1.000000E+01
> DOM, MONIT, 2.000000E+00, 2.000000E+00, 1.428571E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 5, 1.000000E+00
> GRID, RSET_X_2, 9, 1.000000E+00
> GRID, RSET_X_3, 4, 1.000000E+00
> GRID, RSET_Y_1, 5, 1.000000E+00
> GRID, RSET_Y_2, 9, 1.000000E+00
> GRID, RSET_Y_3, 4, 1.000000E+00
> GRID, RSET_Z_1, 7, 1.000000E+00
> GRID, RSET_Z_2, 6, 1.000000E+00
> GRID, RSET_Z_3, 18, 1.000000E+00
> DOM, MOMCEN, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, NAME, B2
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+01
> OBJ, SIZE, 4.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, ROT-MODE, OLD
> OBJ, GEOMETRY, cube12
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, B3
> OBJ, POSITION, 1.000000E+00, 1.000000E+00, 2.000000E+00
> OBJ, SIZE, 2.000000E+00, 2.000000E+00, 2.000000E+00
> OBJ, GEOMETRY, sphere
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
STOP