TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2009, Date 21/07/09
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
autoplot use
AUTOPLOT
fi
monplt.csv 2
da 1 isweep wave
isweep wbar
da 1 isweep w1-1
pl 1
colf 2
enduse
VRV USE
* Start of frame
PROBE 2.500000E-01 2.500000E-01 4.500000E+00
SLICE X
SLICE SAVE
SLICE Y
SLICE SAVE
VARIABLE Pressure
VECTOR ON
VECTOR COLOUR MULTI
CONTOUR ON
PAUSE
enduse
This case is based on Library case 274, which deals with the
flow past a moving van.
An impression of the geometry can be gained by examining
the following sketch. |
______________ |______
|**************| |******|
^ __|**************| ^ |******|
| |*****************| | |******|
y| |*****************| y| |******|
| ~**~~~~~~~~~~~~~**~ | |~~~~**
** ** | **
-------------------------------------------
z----> x--->
Now, In-Form is used to
# calculate the average Z-velocity 15m downstream of the
inlet plane,
# print this value to inforout at the end of the run
# produce a file monplt.csv which contains the variation
of average velocity, and Z-velocity at a location
with sweep.
To see the related tutorial, click
here.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 274
************************************************************
Group 1. Run Title
TEXT(3D Nav.Sto.Flow Around A Moving Van )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,8,12,21)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(149)=IMB1 ;NAME(150)=PRPS
* Solved variables list
SOLVE(P1,U1,V1,W1)
* Stored variables list
STORE(PRPS,IMB1)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(U1,Y,Y,N,Y,N,Y)
SOLUTN(V1,Y,Y,N,Y,N,Y)
SOLUTN(W1,Y,Y,N,Y,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
RHO1 =1.
ENUL =1.0E-05
CP1 =1.
ENUT =1.0E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(W1)=14.
FIINIT(IMB1)=1.001E-10 ;FIINIT(PRPS)=-1.
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
************************************************************
Echo save-block settings for Group 13
save13begin
(make totar is 0.0)
(make w1tot is 0.0)
(store1 of totar at avew1 is sum(ahigh))
(store1 of w1tot at avew1 is sum(w1*ahigh))
(print of w1bar is w1tot/totar)
(table in monplt.csv is get(w1tot/totar, w1{1,1,15}) with head(wbar$
, w1-1)!sweep)
save13end
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(U1 ,FALSDT,0.1 )
RELAX(V1 ,FALSDT,0.1 )
RELAX(W1 ,FALSDT,0.1 )
************************************************************
Group 18. Limits
VARMAX(U1)=1. ;VARMIN(U1)=-1.0E+11
VARMAX(V1)=1. ;VARMIN(V1)=-1.0E+11
VARMAX(W1)=1. ;VARMIN(W1)=-1.0E+11
************************************************************
Group 19. EARTH Calls To GROUND Station
ISG50 = 1
ISG52 = 1
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 5 ;IZMON = 5
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 4
NYPRIN = 2
NZPRIN = 2
YZPR = T
ISWPRF = 1 ;ISWPRL = 100000
PATCH(PROFILE, PROFIL, 7, 7, 6, 6, 1, 20, 1, 1)
PLOT(PROFILE, P1, 0., 0.)
PLOT(PROFILE, U1, 0., 0.)
PLOT(PROFILE, V1, 0., 0.)
PLOT(PROFILE, W1, 0., 0.)
PATCH(XEQ1, CONTUR, 1, 1, 1, 12, 1, 20, 1, 1)
PLOT(XEQ1, P1, 0., 20.)
PLOT(XEQ1, U1, 0., 20.)
PLOT(XEQ1, V1, 0., 20.)
PLOT(XEQ1, W1, 0., 20.)
PATCH(XEQ5, CONTUR, 5, 5, 1, 12, 1, 20, 1, 1)
PLOT(XEQ5, P1, 0., 20.)
PLOT(XEQ5, W1, 0., 20.)
PATCH(WAKE, CONTUR, 1, 8, 1, 12, 15, 15, 1, 1)
PLOT(WAKE, W1, 0., 20.)
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.6239,0.47984,0.616849)
GVIEW(UP,-0.339038,0.877352,-0.339568)
> DOM, SIZE, 5.000000E+00, 6.000000E+00, 2.000000E+01
> DOM, MONIT, 2.550000E+00, 2.250000E+00, 4.500000E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, RSET_X_1, 2, 1.000000E+00
> GRID, RSET_X_2, 1, 1.000000E+00
> GRID, RSET_X_3, 5, 1.000000E+00
> GRID, RSET_Y_1, 1, 1.000000E+00
> GRID, RSET_Y_2, 2, 1.000000E+00
> GRID, RSET_Y_3, 3, 1.000000E+00
> GRID, RSET_Y_4, 6, 1.000000E+00
> GRID, RSET_Z_1, 5, 1.000000E+00
> GRID, RSET_Z_2, 1, 1.000000E+00
> GRID, RSET_Z_3, 3, 1.000000E+00
> GRID, RSET_Z_4, 1, 1.000000E+00
> GRID, RSET_Z_5, 10, 1.000000E+00
> GRID, RSET_Z_6, 1, 1.000000E+00
> DOM, T_AMBIENT, 0.000000E+00
> OBJ, NAME, BODY1
> OBJ, POSITION, 0.000000E+00, 5.000000E-01, 5.000000E+00
> OBJ, SIZE, 1.500000E+00, 1.000000E+00, 5.000000E+00
> OBJ, GEOMETRY, cube14
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, BODY2
> OBJ, POSITION, 0.000000E+00, 1.500000E+00, 6.000000E+00
> OBJ, SIZE, 1.500000E+00, 1.500000E+00, 4.000000E+00
> OBJ, GEOMETRY, cube14
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, FWHL
> OBJ, POSITION, 1.000000E+00, 0.000000E+00, 5.000000E+00
> OBJ, SIZE, 5.000000E-01, 5.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube14
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, RWHL
> OBJ, POSITION, 1.000000E+00, 0.000000E+00, 9.000000E+00
> OBJ, SIZE, 5.000000E-01, 5.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube14
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, UPSTR
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.000000E+00, 6.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 0. ,0. ,14.
> OBJ, NAME, DWSTR
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+01
> OBJ, SIZE, 5.000000E+00, 6.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 0.1
> OBJ, NAME, ROAD
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.000000E+00, 0.000000E+00, 2.000000E+01
> OBJ, GEOMETRY, cube11
> OBJ, TYPE, PLATE
> OBJ, VELOCITY, 0. ,0. ,14.
> OBJ, NAME, AVEW1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.500000E+01
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, default
> OBJ, TYPE, USER_DEFINED
> OBJ, GRID, NO
STOP