TALK=T;RUN(1,1)
DISPLAY
This case is like 369, but use the y-direction instead of x.
An InForm-made box object moves uniformly along a duct, the
cross-section of which it wholly fills.
The duct contains a porous medium causing resistance to flow
The ends of the duct are partially open.
The fluid is compressible.
Parameters are provided which enable the effects on convergence
and accuracy to be noted of:
* box size divided by cell size;
* time to traverse a cell divided by time step;
* the compressibility factor;
* the resistance coefficient;
* pressure-drop coefficients at the ends of the duct.
PHOTON USE commands are included to assist display of results.
ENDDIS
photon use
p
parphi
con v1 x 1 z 2 m fi;0.01
pause
con floy x 1 fi;0.01
pause
con obid x 1 fi;0.01
pause
con p1 x 1 z 2 m fi;0.01
pause
con anp1 x 1 fi ;0.01
pause
con apv1 x 1 fi;0.01
pause
enduse
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(MOFOR - In-Form: BOX linear motion in Y )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = F
* Set overall time and no. of steps
TFIRST =0. ;TLAST =1.
FSTEP = 1 ;LSTEP = 50
TFRAC(1)=0.02 ;TFRAC(2)=0.04
TFRAC(3)=0.06 ;TFRAC(4)=0.08
TFRAC(5)=0.1 ;TFRAC(6)=0.12
TFRAC(7)=0.14 ;TFRAC(8)=0.16
TFRAC(9)=0.18 ;TFRAC(10)=0.2
TFRAC(11)=0.22 ;TFRAC(12)=0.24
TFRAC(13)=0.26 ;TFRAC(14)=0.28
TFRAC(15)=0.3 ;TFRAC(16)=0.32
TFRAC(17)=0.34 ;TFRAC(18)=0.36
TFRAC(19)=0.38 ;TFRAC(20)=0.4
TFRAC(21)=0.42 ;TFRAC(22)=0.44
TFRAC(23)=0.46 ;TFRAC(24)=0.48
TFRAC(25)=0.5 ;TFRAC(26)=0.52
TFRAC(27)=0.54 ;TFRAC(28)=0.56
TFRAC(29)=0.58 ;TFRAC(30)=0.6
TFRAC(31)=0.62 ;TFRAC(32)=0.64
TFRAC(33)=0.66 ;TFRAC(34)=0.68
TFRAC(35)=0.7 ;TFRAC(36)=0.72
TFRAC(37)=0.74 ;TFRAC(38)=0.76
TFRAC(39)=0.78 ;TFRAC(40)=0.8
TFRAC(41)=0.82 ;TFRAC(42)=0.84
TFRAC(43)=0.86 ;TFRAC(44)=0.88
TFRAC(45)=0.9 ;TFRAC(46)=0.92
TFRAC(47)=0.94 ;TFRAC(48)=0.96
TFRAC(49)=0.98 ;TFRAC(50)=1.
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 1
XULAST =2.
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 100
YVLAST =1.
YFRAC(1)=4.999995E-03 ;YFRAC(2)=0.015051
YFRAC(3)=0.025101 ;YFRAC(4)=0.035152
YFRAC(5)=0.045202 ;YFRAC(6)=0.055253
YFRAC(7)=0.065303 ;YFRAC(8)=0.075354
YFRAC(9)=0.085404 ;YFRAC(10)=0.095455
YFRAC(11)=0.105505 ;YFRAC(12)=0.115556
YFRAC(13)=0.125606 ;YFRAC(14)=0.135657
YFRAC(15)=0.145707 ;YFRAC(16)=0.155758
YFRAC(17)=0.165808 ;YFRAC(18)=0.175859
YFRAC(19)=0.185909 ;YFRAC(20)=0.19596
YFRAC(21)=0.20601 ;YFRAC(22)=0.216061
YFRAC(23)=0.226111 ;YFRAC(24)=0.236162
YFRAC(25)=0.246212 ;YFRAC(26)=0.256263
YFRAC(27)=0.266313 ;YFRAC(28)=0.276364
YFRAC(29)=0.286414 ;YFRAC(30)=0.296465
YFRAC(31)=0.306515 ;YFRAC(32)=0.316566
YFRAC(33)=0.326616 ;YFRAC(34)=0.336667
YFRAC(35)=0.346717 ;YFRAC(36)=0.356768
YFRAC(37)=0.366818 ;YFRAC(38)=0.376869
YFRAC(39)=0.386919 ;YFRAC(40)=0.39697
YFRAC(41)=0.40702 ;YFRAC(42)=0.417071
YFRAC(43)=0.427121 ;YFRAC(44)=0.437172
YFRAC(45)=0.447222 ;YFRAC(46)=0.457273
YFRAC(47)=0.467323 ;YFRAC(48)=0.477374
YFRAC(49)=0.487424 ;YFRAC(50)=0.497475
YFRAC(51)=0.507525 ;YFRAC(52)=0.517576
YFRAC(53)=0.527626 ;YFRAC(54)=0.537677
YFRAC(55)=0.547727 ;YFRAC(56)=0.557778
YFRAC(57)=0.567828 ;YFRAC(58)=0.577879
YFRAC(59)=0.587929 ;YFRAC(60)=0.59798
YFRAC(61)=0.60803 ;YFRAC(62)=0.618081
YFRAC(63)=0.628131 ;YFRAC(64)=0.638182
YFRAC(65)=0.648232 ;YFRAC(66)=0.658283
YFRAC(67)=0.668333 ;YFRAC(68)=0.678384
YFRAC(69)=0.688434 ;YFRAC(70)=0.698485
YFRAC(71)=0.708535 ;YFRAC(72)=0.718586
YFRAC(73)=0.728636 ;YFRAC(74)=0.738687
YFRAC(75)=0.748737 ;YFRAC(76)=0.758788
YFRAC(77)=0.768838 ;YFRAC(78)=0.778889
YFRAC(79)=0.788939 ;YFRAC(80)=0.79899
YFRAC(81)=0.80904 ;YFRAC(82)=0.819091
YFRAC(83)=0.829141 ;YFRAC(84)=0.839192
YFRAC(85)=0.849242 ;YFRAC(86)=0.859293
YFRAC(87)=0.869343 ;YFRAC(88)=0.879394
YFRAC(89)=0.889444 ;YFRAC(90)=0.899495
YFRAC(91)=0.909545 ;YFRAC(92)=0.919596
YFRAC(93)=0.929646 ;YFRAC(94)=0.939697
YFRAC(95)=0.949748 ;YFRAC(96)=0.959798
YFRAC(97)=0.969849 ;YFRAC(98)=0.979899
YFRAC(99)=0.989949 ;YFRAC(100)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =0.1
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(143)=FLOY ;NAME(144)=RSV1
NAME(145)=ANP1 ;NAME(146)=APV1
NAME(147)=SOV1 ;NAME(148)=RHO1
NAME(149)=PRPS ;NAME(150)=OBID
* 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,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(FLOY,Y,N,N,N,N,Y)
SOLUTN(RSV1,Y,N,N,N,N,Y)
SOLUTN(ANP1,Y,N,N,N,N,Y)
SOLUTN(APV1,Y,N,N,N,N,Y)
SOLUTN(SOV1,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(OBID,Y,N,N,N,N,Y)
DEN1 = 148
PRPS = 149
************************************************************
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(V1,Y,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
DENPCO = T
************************************************************
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(V1)=1.
PRT(V1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
FIINIT(FLOY)=1.0E-10 ;FIINIT(RSV1)=1.0E-10
FIINIT(ANP1)=1.0E-10 ;FIINIT(APV1)=1.0E-10
FIINIT(SOV1)=1.0E-10 ;FIINIT(RHO1)=1.0E-10
FIINIT(PRPS)=-1. ;FIINIT(OBID)=1.0E-10
PATCH(^OB3 ,INIVAL, 3, 0, 0, 0, 0, 0, 1, 1)
INIT(^OB3 ,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(RESIST ,VOLUME, 1, 1, 1, 100, 1, 1, 1, 50)
COVAL(RESIST ,V1 ,10. ,0. )
PATCH(INLET ,CELL , 1, 1, 1, 1, 1, 1, 1, 50)
COVAL(INLET ,P1 ,1. ,0. )
PATCH(OUTLET ,CELL , 1, 1, 100, 100, 1, 1, 1, 50)
COVAL(OUTLET ,P1 ,1. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 60 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,1.)
RELAX(FLOY,LINRLX,1.)
RELAX(RSV1,LINRLX,1.)
RELAX(ANP1,LINRLX,1.)
RELAX(APV1,LINRLX,1.)
RELAX(SOV1,LINRLX,1.)
RELAX(RHO1,LINRLX,0.1)
RELAX(PRPS,LINRLX,1.)
RELAX(OBID,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(FLOY)=1.0E+10 ;VARMIN(FLOY)=-1.0E+10
VARMAX(RSV1)=1.0E+10 ;VARMIN(RSV1)=-1.0E+10
VARMAX(ANP1)=1.0E+10 ;VARMIN(ANP1)=-1.0E+10
VARMAX(APV1)=1.0E+10 ;VARMIN(APV1)=-1.0E+10
VARMAX(SOV1)=1.0E+10 ;VARMIN(SOV1)=-1.0E+10
VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(OBID)=1.0E+10 ;VARMIN(OBID)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
CONWIZ = T
ISG52 = 2
ISG62 = 0
SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SPEDAT(SET,PROPERTY,RHO1,C,=1.0+1.0E-04*P1)
SPEDAT(SET,PROPERTY,DRH1DP,C,=1.0E-04)
SPEDAT(SET,MOVOB,BOX,C,=POS(0.0&0.15*YVLAST+0.666667*TIM&0.0&0&0&$)
SPEDAT(SET,MOVOB,BOX,C,0&0&0&0&0))
SPEDAT(SET,STORED,SOV1,C,=SORC(V1)!SOURCS)
SPEDAT(SET,STORED,APV1,C,=APCO(V1)!EQCOEF)
SPEDAT(SET,STORED,ANP1,C,=ANCO(P1)!EQCOEF)
SPEDAT(SET,STORED,RSV1,C,=RESI(V1)!RESIDU)
SPEDAT(SET,STORED,FLOY,C,=RHO1*V1!ZSLFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,^OB3,C,BOX)
SPEDAT(SET,OBJTYP,^OB3,C,BLOCKAGE)
SPEDAT(SET,BOX,MATERIAL,R,198.)
SPEDAT(SET,FACETDAT,NUMOBJ,I,3)
SPEDAT(SET,MATERIAL,198,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=2.807 ;EX(V1)=0.6666
EX(FLOY)=0.6598 ;EX(RSV1)=7.599E-06
EX(ANP1)=0.1535 ;EX(APV1)=1000.
EX(SOV1)=9.851001E-04 ;EX(RHO1)=0.9999
EX(PRPS)=0. ;EX(OBID)=0.
************************************************************
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(V1,Y,N,Y,Y,Y,Y)
OUTPUT(FLOY,Y,N,Y,N,N,N)
OUTPUT(RSV1,Y,N,Y,N,N,N)
OUTPUT(ANP1,Y,N,Y,N,N,N)
OUTPUT(APV1,Y,N,Y,N,N,N)
OUTPUT(SOV1,Y,N,Y,N,N,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(OBID,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 26 ;IZMON = 1
NPRMON = 60 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NTPRIN = 1 ;ISTPRF = 1 ;ISTPRL = 100000
NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000
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 = 0 ;IDISPC = 0
STOP