TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2007, Date 21/12/07
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
This parameterised Q1 permits study of the angled-plate
feature described in Encyclopaedia Article 'angled plates'.
Cases 1 to 4 concern variously-placed single plates.
Case 5 concerns two plates, positioned so as to represent a
tent or a pitched roof.
The SAVEnBEGIN, SAVnEND, insertions ensure that , if the file
is read by the VR-Editor and then re-written by it, the
parameter-using features are preserved.
ENDDIS
VRV USE
SLICE Y
VIEW 0.000000E+00 -1.000000E+00 0.000000E+00
UP 0.000000E+00 0.000000E+00 1.000000E+00
VARIABLE Pressure
CONTOUR ON
PAUSE
* Start of frame
VARIABLE Velocity
VECTOR ON; CONTOUR OFF
PAUSE
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(angled plate(s) in xz plane caseno=1 )
************************************************************
Echo InForm settings for Group 1
SAVE1BEGIN
real(u1in,w1in,rhoin)
boolean(porous,lowin,hout,block1,block2)
integer(caseno)
u1in=1.0
w1in=0.0
porous=f
xpos=0.3
zpos=0.3
matnum=0
mesg(caseno=:caseno: OK? If not choose other
readvdu(caseno,int,0)
caseno
case caseno of
when 1
+ block1=t
+ block2=f
+ xsiz=0.4
+ zsiz=0.4
+ w1in=u1in
+ lowin=t
+ hout=t
when 2
+ block1=t
+ block2=f
+ xsiz=0.4
+ zsiz=0.4
+ w1in=0.0
+ lowin=f
+ hout=f
when 3
+ block1=t
+ block2=f
+ xsiz=0.4
+ zsiz=0.2
+ w1in=0.0
+ lowin=f
+ hout=f
when 4
+ block1=f
+ block2=t
+ xsiz=0.4
+ zsiz=0.2
+ w1in=0.0
+ lowin=f
+ hout=f
when 5
+ block1=t
+ block2=t
+ xsiz=0.2
+ zsiz=0.2
+ w1in=0.0
+ lowin=f
+ hout=f
orelse
+ block1=t
+ block2=t
+ xsiz=0.2
+ zsiz=0.1
+ lowin=f
+ hout=f
endcase
if(matnum.eq.0) then
rhoin=1.189
else
matnum=67
rhoin=1000.0
endif
masflo=rhoin*u1in
TEXT(angled plate(s) in xz plane caseno=:caseno:
SAVE1END
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,51,1,51)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(149) =RHO1 ; NAME(150) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,W1 )
* Stored variables list
STORE(PRPS,RHO1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(U1 ,Y,Y,Y,N,N,Y)
SOLUTN(W1 ,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
ENUT = 0.000000E+00
DVO1DT = 3.410000E-03
************************************************************
Echo InForm settings for Group 9
SAVE9BEGIN
SETPRPS(1, MATNUM)
SAVE9END
************************************************************
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
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 300
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
PARSOL = T
CONWIZ = T
ISG50 = 1
ISG52 = 2
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 3 ;IYMON = 1 ;IZMON = 3
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='P'
CSG2 ='X'
************************************************************
Echo InForm settings for Group 25
SAVE25BEGIN
GVIEW(P,0.000000E+00,-1.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)
> DOM, SIZE, 1.000000E+00, 2.000000E-01, 1.000000E+00
> DOM, MONIT, 5.000000E-02, 1.000000E-01, 5.000000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
if(block1) then
> OBJ, NAME, BLOCK
> OBJ, POSITION, xpos, 0.000000E+00, zpos
> OBJ, SIZE, xsiz, 2.000000E-01, zsiz
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, wedge2
> OBJ, ROTATION24, 2
> OBJ, TYPE, BLOCKAGE
> OBJ, WIREFRAME, YES
> OBJ, MATERIAL, 299
endif
if(block2) then
> OBJ, NAME, Block2
> OBJ, POSITION, 1.0-xpos-xsiz, 0.000000E+00, xpos
> OBJ, SIZE, xsiz, 2.000000E-01, zsiz
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, wedge2
> OBJ, ROTATION24, 4
> OBJ, TYPE, BLOCKAGE
> OBJ, WIREFRAME, YES
> OBJ, MATERIAL, 299
endif
if(lowin) then
> OBJ, NAME, LINLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 2.000000E-01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, u1in, 0.000000E+00, w1in
> OBJ, voluflow, w1in
> OBJ, TEMPERATURE, 0.000000E+00
endif
> OBJ, NAME, winLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 2.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, u1in, 0.000000E+00, w1in
> OBJ, MASSFLOW, masflo
> OBJ, voluflow, u1in
> OBJ, TEMPERATURE, 0.000000E+00
if(hout) then
> OBJ, NAME, HOUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ, SIZE, 1.000000E+00, 2.000000E-01, 0.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, COEFFICIENT, 1.000000E-03*rhoin
endif
> OBJ, NAME, EOUTLET
> OBJ, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 2.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, COEFFICIENT, 1.000000E-03*rhoin
if(porous) then
PATCH(POROS,CELL,0,0,0,0,0,0,1,1)
COVAL(POROS,U1,1.E3,0.)
COVAL(POROS,W1,1.E3,0.)
endif
NXPRIN=1
NZPRIN=1
(INTREAD CUTCDB is -1)
INTEGER(IPBPRF,IPBPRL,NPBPRI)
IPBPRF=1;IPBPRL=5;NPBPRI=2
(INTREAD CUTCPR is :IPBPRF: :IPBPRL: :NPBPRI:)
SAVE25END
STOP