TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.5, Date 10/02/03
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
TEXT(MOFOR by In-Form: 2D motion of 2 objects
libref=360
TITLE
DISPLAY
This input file shows how MOFOR, the moving-frame-of-reference
feature, can be used without a .MOF file, when suitable In-Form
statements are used to define its motion, as described
here.
Two cylinders move through a fluid. The first follows a parabolic
trajectory from left to right, the second a horizontal trajectory
from right to left.
'(MOVOB' statements dictate these motions of by use of the POS()
function, of which the arguments set:
X, Y and Z coordinates and rotation angles about X, Y and Z axes,
as functions of the time variable, TIM.
The variable OBID (= object identifier) marks the positions, for
ease of display by VIEWER, PHOTON or WinPHOTON.
ENDDIS
VRV USE
* Domain scaling factors
SCALE 1.000000E+00 1.000000E+00 1.000000E+00
* Settings for current slice
PROBE 1.070000E+01 3.500000E+00 5.000000E-02; PROBE ON
SLICE Z
* View and up directions
VIEW 0.000000E+00 0.000000E+00 1.000000E+00
UP 0.000000E+00 1.000000E+00 0.000000E+00
* View centre
VIEW CENTRE 1.100001E+01 3.500000E+00 5.000000E-02
* View size
VIEW SIZE 1.038916E+01
* View perspective
VIEW DEPTH 3.000000E+00;VIEW TILT 0.8
DOMAIN ON
VARIABLE Velocity
VARIABLE RANGE 0.000000E+00 1.400000E+01
VECTOR ON; CONTOUR OFF; SURFACE OFF
VECTOR REFERENCE 5.000000E+00
ANIMATE START 1 END 20 INTERVAL 1
ENDUSE
photon use
p
m1;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m2;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m3;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m4;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m5;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m6;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m7;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m8;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m9;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m10;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m11;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m12;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m13;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m14;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m15;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m16;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m17;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m18;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m19;
gr ou z 1
ve z 1 sh
surf obid z 1
upause 2
p
m20;
gr ou z 1
ve z 1 sh
surf obid z 1
enduse
************************************************************
IRUNN = 1 ;LIBREF = 360
************************************************************
Group 1. Run Title
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.000000E+00,2.000000E+00,20)
* Modify regions
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,110,35,1,.01)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(149) =OBID ; NAME(150) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 )
* Stored variables list
STORE(PRPS,OBID)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
************************************************************
Echo InForm settings for Group 7
INFORM7BEGIN
** Definition of the VR moving objects by In-Form
char(vel,gravt,times)
vel=10.; gravt=9.81; times=tim
(MOVOB of CYL1 is POS(:times:*:vel:,:times:*:vel:-0.5*:gravt:*:t$
imes:^2,0,0,0,0))
(MOVOB of CYL2 is POS(-:times:*:vel:,0,0,0,0,0))
store(obol)
(stored var obol is old(obid))
INFORM7END
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
SETPRPS(1, 0)
ENUT = 0.000000E+00
DRH1DP = 1.000000E-05
DVO1DT = 3.410000E-03
EMISS = 0.000000E+00 ;SCATT = 0.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = 0.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
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 20
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 5.000000E-01)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
CSG1 ='M'
ISG62 = 0
SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
EX (P1 ) = 2.451000E+01 ;EX (U1 ) = 3.334000E+00
EX (V1 ) = 1.756000E+00 ;EX (OBID) = 3.896000E-02
EX (PRPS) = 0.000000E+00 ;EX (OBOL) = 3.948000E-02
************************************************************
Group 21. Print-out of Variables
OUTPUT(V1 ,Y,Y,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 54 ;IYMON = 18 ;IZMON = 1
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
NOWIPE = T
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='M'
************************************************************
Echo InForm settings for Group 25
INFORM25BEGIN
DEBUG=f
MOFOR=t
formula=f
INFORM25END
GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 2.200000E+01, 7.000000E+00, 1.000000E-01
> DOM, MONIT, 1.070000E+01, 3.500000E+00, 5.000000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, RELIEF
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.000000E-01, 2.000000E-01, 1.000000E-01
> OBJ, CLIPART, cubet
> OBJ, ROTATION24, 1
> OBJ, TYPE, PRESSURE_RELIEF
> OBJ, PRES_RELIEF, 1.000000E+03, 0.000000E+00
> OBJ, NAME, CYL1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ, CLIPART, cylinder
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, -1
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, INI_PRESS, 0.000000E+00
> OBJ, SCAL_FIXF, 0.000000E+00
> OBJ, NAME, CYL2
> OBJ, POSITION, 2.100000E+01, 3.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ, CLIPART, cylinder
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, -1
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, INI_PRESS, 0.000000E+00
> OBJ, SCAL_FIXF, 0.000000E+00
> OBJ, NAME, WALL
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.200000E+01, 0.000000E+00, 1.000000E-01
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
libref=360
STOP