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 quadratic velocity
U(t)= t - t**2.0
or with varing accelaration
a = 1-2.0*t
The method is implemented by In-Form formulea written in
the Q1 file.
The above velocity is used to specify the mass flow and
velocity boundary conditions at INLET; and the above
accelation is used to provide 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 Pressure; CON ON
VECTOR ON
ANIMATE 1 20 1
PAUSE
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 805
************************************************************
Group 1. Run Title
TEXT(Moving sphere with varing accelaration )
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.000000E+00,1.000000E+00,20)
* 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
SAVE13BEGIN
patch(in,low,1,nx,1,ny,1,1,1,lstep)
(source of p1 at in is (tim-tim^2.0)*rho1)
(source of w1 at in is (tim-tim^2.0) with onlyms)
patch(acel,phasem,1,nx,1,ny,1,nz,1,lstep)
(source of w1 at acel is 1.-2*(tim-0.5*dt))
SAVE13END
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 20
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
************************************************************
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'
EX(U1)=2.633E-04; EX(V1)=2.633E-04; EX(W1)=4.720E-04
GVIEW(P,-9.984754E-01,5.497228E-02,4.992418E-03)
GVIEW(UP,5.497159E-02,9.984879E-01,-2.748602E-04)
> 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