TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 01/04/04
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
This case is like 375, but use the xz-plane instead of xy.
A wide thin box moves diagonally across a two-dimensional space
at uniform velocity. A box is an InForm-made object.
The domain contains a porous medium causing resistance to flow.
The boundaries of the domain are open.
The fluid is slightly 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 domain.
PHOTON USE commands are included to assist display of results.
ENDDIS
photon use
p
p1;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p2;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p3;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p4;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p5;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p6;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p7;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p8;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p9;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p10;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p11;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p12;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p13;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p14;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p15;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p16;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p17;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p18;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p19;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
upause 1
p
p20;
gr ou y 1
set prop off
ve y 1 sh
surf obid y 0.1
enduse
************************************************************
IRUNN = 1 ;LIBREF = 377
************************************************************
GROUP 1. Run title and other preliminaries
TEXT(MOFOR/In-Form: diagonal motion in XZ
save1begin
Initial data of problem
-----------------------
delaration of parameters
REAL(BZSZFAC,BXSZFAC,ANGL,TIMEFAC,COMPFAC,PCO,RESCO) ! delaration
REAL(DMSIZE,PI) ! of parameters
! settings
BZSZFAC=9.0 ! box Z size divided by cell size
BXSZFAC=1.0 ! box X size divided by cell size
ANGL=-45. ! BOX rotation angle about Y axis
TIMEFAC=0.75 ! time to traverse a cell divided by time step
COMPFAC=0.01 ! compressibility factor
PCO=1.0 ! pressure coefficient
RESCO=1.0 ! resistance to flow coefficient
DMSIZE=1.0 ! domain size
PI=3.14159 ! Pi number
save1end
TLAST=10.0 ! PIL variables
ZWLAST=DMSIZE
XULAST=DMSIZE
LSTEP=20
NZ=40
NX=40
GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,LSTEP,TLAST,1.0)
GROUP 3. X-direction grid specification
GRDPWR(X,NX,XULAST,1.0)
GROUP 4. Y-direction grid specification
GRDPWR(Y,1,0.1,1.0)
GROUP 5. Z-direction grid specification
GRDPWR(Z,NZ,ZWLAST,1.0)
GROUP 7. Variables stored, solved & named
SOLVE(P1,W1,U1)
SOLUTN(P1,Y,Y,Y,N,N,Y)
STORE(OBID,PRPS)
GROUP 13. Boundary conditions and special sources
** Resistance
PATCH(RESIST,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(RESIST,W1,RESCO,0.0)
COVAL(RESIST,U1,RESCO,0.0)
** open boundaries
PATCH(low z,CELL,1,NX,1,1,1,1,1,LSTEP)
COVAL(low z,P1,PCO,0.0)
PATCH(high z,CELL,1,NX,1,1,NZ,NZ,1,LSTEP)
COVAL(high z,P1,PCO,0.0)
PATCH(low x,CELL,1,1,1,1,1,NZ,1,LSTEP)
COVAL(low x,P1,PCO,0.0)
PATCH(high x,CELL,NX,NX,1,1,1,NZ,1,LSTEP)
COVAL(high x,P1,PCO,0.0)
** The space within which box may move
PATCH(IMOFOR,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
INFORM13BEGIN
store(rho1)
(property rho1 is 1.0+ :compfac:*p1) ! density depends on pressure
(property drh1dp is :compfac:)
denpco=t
** Moving box
REAL(ZCELL,XCELL,TSTEP,TCELL,REALNZ,REALNX,REALNT,WVEL,UVEL) ! declaration
! of reals
realnz=nz; realnx=nx
realnt=lstep
tstep=tlast/realnt ! size of time step
zcell=zwlast/realnz ! Z size of cell
xcell=xulast/realnx ! X size of cell
tcell=tstep*timefac ! time to traverse cell
wvel=zcell/tcell ! W1 velocity of box motion
uvel=xcell/tcell ! U1 velocity of box motion
REAL(ZSIZE,XSIZE,RTMP,SINA,COSA)
CHAR(XPS,YPS,ZPS,VEL)
real(ang); ang=pi*angl/180.
sina=sin(ang)
cosa=cos(ang)
zsize=zcell*bzszfac
xsize=xcell*bxszfac
vel=:wvel:
rtmp=.1*zwlast-0.5*zsize*cosa+0.5*xsize*sina
zps=:vel:*tim+(:rtmp:)
vel=:uvel:
rtmp=.1*xulast-0.5*zsize*sina-0.5*xsize*cosa
xps=:vel:*tim+(:rtmp:)
yps=.0
vel
zps
xps
(MOVOB of BOX is POS(:XPS:&:YPS:&:ZPS:&0&:ANGL:&0))
INFORM13END
SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET)
ISG62 = 0
GROUP 15. Termination of sweeps
LSWEEP=100
GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.1)
RELAX(U1,FALSDT,0.1*TCELL)
RELAX(W1,FALSDT,0.1*TCELL)
RELAX(RHO1,LINRLX,0.1)
varmax(p1)=1.0
varmin(p1)=-1.0
varmin(rho1)=0.01
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
GROUP 22. Spot-value print-out
NPRMON=LSWEEP;IXMON=NX/4+1;IYMON=1
IZMON=-1 ! use negative value to cause graphics monitor
! to display maximum and minumum values
GROUP 23. Print-out & plot control
TSTSWP=-1
IDISPA=1; CSG1=P
LIBREF=377
SELREF=T; RESFAC=1.E-2
NTPRIN=1
DISTIL=T
EX(P1)=1.790E-03; EX(U1)=8.849E-03; EX(W1)=8.848E-03
EX(RHO1)=1.000E+00
************************************************************
GVIEW(P,0.000000E+00,1.000000E+00,0.000000E+00)
GVIEW(UP,1.000000E+00,0.000000E+00,0.000000E+00)
> DOM, SIZE, 1.000000E+00, 1.000000E-01, 1.000000E+00
> DOM, MONIT, 2.625000E-01, 5.000000E-02, 1.250000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, LOWZ
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E-01, 2.500000E-02
> OBJ, GEOMETRY, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, HIGHZ
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 9.750000E-01
> OBJ, SIZE, 1.000000E+00, 1.000000E-01, 2.499998E-02
> OBJ, GEOMETRY, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, LOWX
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.500000E-02, 1.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, HIGHX
> OBJ, POSITION, 9.750000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.499998E-02, 1.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, CELLTYPE
> OBJ, NAME, BOX
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.500000E-02, 1.000000E-01, 2.250000E-01
> OBJ, GEOMETRY, cube14
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
lsg57=t
#maxabs
STOP