TALK=T;RUN(1,1)
   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
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Moving sphere with constant accelaration)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =10.
 FSTEP = 1 ;LSTEP = 10
 TFRAC(1)=0.1 ;TFRAC(3)=0.3
 TFRAC(5)=0.5 ;TFRAC(7)=0.7
 TFRAC(9)=0.9
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 18
 XULAST =4.
 XFRAC(1)=0.05 ;XFRAC(2)=0.1
 XFRAC(3)=0.15 ;XFRAC(4)=0.2
 XFRAC(5)=0.25 ;XFRAC(6)=0.305556
 XFRAC(7)=0.361111 ;XFRAC(8)=0.416667
 XFRAC(9)=0.472222 ;XFRAC(10)=0.527778
 XFRAC(11)=0.583333 ;XFRAC(12)=0.638889
 XFRAC(13)=0.694444 ;XFRAC(14)=0.75
 XFRAC(15)=0.8125 ;XFRAC(16)=0.875
 XFRAC(17)=0.9375 ;XFRAC(18)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 18
 YVLAST =4.
 YFRAC(1)=0.05 ;YFRAC(2)=0.1
 YFRAC(3)=0.15 ;YFRAC(4)=0.2
 YFRAC(5)=0.25 ;YFRAC(6)=0.305556
 YFRAC(7)=0.361111 ;YFRAC(8)=0.416667
 YFRAC(9)=0.472222 ;YFRAC(10)=0.527778
 YFRAC(11)=0.583333 ;YFRAC(12)=0.638889
 YFRAC(13)=0.694444 ;YFRAC(14)=0.75
 YFRAC(15)=0.8125 ;YFRAC(16)=0.875
 YFRAC(17)=0.9375 ;YFRAC(18)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 31
 ZWLAST =10.
 ZFRAC(1)=0.028571 ;ZFRAC(2)=0.057143
 ZFRAC(3)=0.085714 ;ZFRAC(4)=0.114286
 ZFRAC(5)=0.142857 ;ZFRAC(6)=0.171429
 ZFRAC(7)=0.2 ;ZFRAC(8)=0.233333
 ZFRAC(9)=0.266667 ;ZFRAC(10)=0.3
 ZFRAC(11)=0.333333 ;ZFRAC(12)=0.366667
 ZFRAC(13)=0.4 ;ZFRAC(14)=0.433333
 ZFRAC(15)=0.466667 ;ZFRAC(16)=0.5
 ZFRAC(17)=0.533333 ;ZFRAC(18)=0.566667
 ZFRAC(19)=0.6 ;ZFRAC(20)=0.633333
 ZFRAC(21)=0.666667 ;ZFRAC(22)=0.7
 ZFRAC(23)=0.733333 ;ZFRAC(24)=0.766667
 ZFRAC(25)=0.8 ;ZFRAC(26)=0.833333
 ZFRAC(27)=0.866667 ;ZFRAC(28)=0.9
 ZFRAC(29)=0.933333 ;ZFRAC(30)=0.966667
 ZFRAC(31)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(149)=VLSQ ;NAME(150)=PRPS
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,Y,N,N,Y)
 SOLUTN(W1,Y,Y,Y,N,N,Y)
 SOLUTN(VLSQ,Y,N,N,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 PRPS = 150
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(P1,Y,Y,Y,N,Y,Y)
 TERMS(U1,Y,Y,Y,Y,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1. ;TMP1 =0. ;EL1 =0.
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.0E-05 ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1.
 CP1 =1005. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
 FIINIT(VLSQ)=1.0E-10 ;FIINIT(PRPS)=-1.
 
 PATCH(^OB2 ,INIVAL, 2, 0, 0, 0, 0, 0, 1, 1)
 INIT(^OB2 ,PRPS,0. ,198. )
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(IN ,LOW , 1, 18, 1, 18, 1, 1, 1, 10)
 COVAL(IN ,P1 ,In-Form:source - see Grp 19)
 COVAL(IN ,W1 ,In-Form:source - see Grp 19)
 
 PATCH(ACEL ,PHASEM, 1, 18, 1, 18, 1, 31, 1, 10)
 COVAL(ACEL ,W1 ,In-Form:source - see Grp 19)
 
 PATCH(OB1 ,HIGH , 1, 18, 1, 18, 31, 31, 1, 10)
 COVAL(OB1 ,P1 ,1000. ,0. )
 COVAL(OB1 ,U1 ,0. ,0. )
 COVAL(OB1 ,V1 ,0. ,0. )
 COVAL(OB1 ,W1 ,0. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 50 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=250 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(VLSQ,LINRLX,1.)
 RELAX(PRPS,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(VLSQ)=1.0E+10 ;VARMIN(VLSQ)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 CALFOR = T
 ISG62 = 1
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 SPEDAT(SET,SOURCE,R1!IN,C,=TIM*RHO1)
 SPEDAT(SET,SOURCE,W1!IN,C,=TIM!ONLYMS)
 SPEDAT(SET,SOURCE,W1!ACEL,C,=1.0)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,DRAG,MOMCENX,R,0.)
 SPEDAT(SET,DRAG,MOMCENY,R,0.)
 SPEDAT(SET,DRAG,MOMCENZ,R,0.)
 SPEDAT(SET,OBJNAM,!OB1,C,B2)
 SPEDAT(SET,OBJTYP,!OB1,C,OUTLET)
 SPEDAT(SET,ARATIO,!OB1,R,1.)
 SPEDAT(SET,OBJNAM,^OB2,C,B3)
 SPEDAT(SET,OBJTYP,^OB2,C,BLOCKAGE)
 SPEDAT(SET,B3,MATERIAL,R,198.)
 SPEDAT(SET,FACETDAT,NUMOBJ,I,2)
 SPEDAT(SET,MATERIAL,198,L,T)
 RG( 1) =2.
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=9.743 ;EX(U1)=0.2959
 EX(V1)=0.2959 ;EX(W1)=10.07
 EX(VLSQ)=0. ;EX(PRPS)=3.354
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 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
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 2 ;ISTPRF = 1 ;ISTPRL = 100000
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = T
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
 IDISPA = 1 ;IDISPB = 1 ;IDISPC = 100
 CSG1    ='M'
STOP