TALK=T;RUN( 1, 1)
DISPLAY
Simulation of laminar steady flow around a sphere in a uniform
stream at a Reynolds number of 40 .
The grid is cartesian and rather coarse (18*18*31), but the flow
follows the curved surface smoothly because the PARSOL (i.e.
cut-cell) technique is used.
The fluid is atmospheric air, treated as incompressible.
The variables computed are:
* p1, the pressure,
* u1, the x-direction velocity,
* v1, the y-direction velocity, and
* w1, the z-direction velocity.
This file contains a macro of commands which cause the Viewer
(when the macro button is pressed) to display contours
automatically.
ENDDIS
VRV USE
* Start of frame
VARIABLE Pressure; CON ON
msg Pressure contours
PAUSE
* Start of frame
VARIABLE Velocity; VEC ON
msg Pressure contours superimposed by velocity vectors
PAUSE
* Start of frame
* Setting object visibility and painting status
OBJECT SHOW TYPE BLOCKAGE
OBJECT PAINT TYPE BLOCKAGE ON
VARIABLE Velocity; VEC OFF
CONTOUR OFF
msg Velocity contours on the surface of the sphere
PAUSE
* Start of frame
* Setting object visibility and painting status
OBJECT SHOW TYPE BLOCKAGE
OBJECT PAINT TYPE BLOCKAGE ON
VARIABLE Pressure; VEC OFF
msg Pressure contours on the surface of the sphere
PAUSE
ENDUSE
************************************************************
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 805
************************************************************
Group 1. Run Title
TEXT(Sphere In A Uniform Stream; Re=40
TITLE
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,18,18,31,0.001)
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 4.000000E+00; YSI= 4.000000E+00; ZSI= 1.000000E+01
RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(150) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 )
* Stored variables list
STORE(PRPS,VLSQ)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(U1 ,Y,Y,Y,N,N,Y)
SOLUTN(V1 ,Y,Y,Y,N,N,Y)
SOLUTN(W1 ,Y,Y,Y,N,N,Y)
inform7begin
(stored var STGP is P1+0.5*RHO1*VLSQ with SWPFIN)
inform7end
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
RHO1 = 1.000000E+00
CP1 = 1.005000E+03
ENUL = 5.000000E-02 ;ENUT = 0.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
No PATCHes used for this Group
FIINIT(W1) = 1.0
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
Imbalance patches enclosing the sphere, which enable the drag on
the sphere to be deduced.
All boxes should report the same momentum balance, if the solution
has converged
PATCH(IMBL3&2,VOLUME,3,NX-2,3,NY-2, 3,NZ-2,1,1) ! the largest box
COVAL(IMBL3&2,W1,0,0)
PATCH(IMBL4&3,VOLUME,3,NX-2,3,NY-2, 4,NZ-3,1,1)
COVAL(IMBL4&3,W1,0,0)
PATCH(IMBL5&4,VOLUME,3,NX-2,3,NY-2, 5,NZ-4,1,1)
COVAL(IMBL5&4,W1,0,0)
PATCH(IMBL810,VOLUME,3,NX-2,3,NY-2, 8,NZ-10,1,1)
COVAL(IMBL810,W1,0,0)
PATCH(IMBL815,VOLUME,5,NX-4,5,NY-4, 8,NZ-15,1,1) ! the smallest box
COVAL(IMBL815,W1,0,0)
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 30
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1) = 250
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-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
************************************************************
Group 20. Preliminary Printout
ECHO = T
distil=t
EX(P1 )= 1.298E-01
EX(U1 )= 3.240E-02
EX(V1 )= 3.335E-02
EX(W1 )= 1.012E+00
EX(PRPS)= 3.354E+00
************************************************************
Group 21. Print-out of Variables
OUTPUT(U1 ,Y,Y,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = nx/2 ;IYMON = ny/2 ;IZMON = - 3*nz/4
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
NOWIPE = T
GVIEW(P,-9.928086E-01,0.000000E+00,-1.197123E-01)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 4.000000E+00, 4.000000E+00, 1.000000E+01
> DOM, MONIT, 5.555560E-01, 5.555560E-01, 4.838710E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> DOM, RELAX, 5.000000E-01
> OBJ1, NAME, B1
> OBJ1, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ1, SIZE, 4.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ1, CLIPART, cube3
> OBJ1, ROTATION24, 1
> OBJ1, TYPE, INLET
> OBJ1, PRESSURE, 0.000000E+00
> OBJ1, VELOCITY, 0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ1, TEMPERATURE, 0.000000E+00
> OBJ2, NAME, B2
> OBJ2, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+01
> OBJ2, SIZE, 4.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ2, CLIPART, cube12
> OBJ2, ROTATION24, 1
> OBJ2, TYPE, OUTLET
> OBJ2, PRESSURE, 0.000000E+00
> OBJ2, TEMPERATURE, 0.000000E+00
> OBJ2, COEFFICIENT, 1.000000E+03
> OBJ3, NAME, B3
> OBJ3, POSITION, 1.000000E+00, 1.000000E+00, 2.200000E+00
> OBJ3, SIZE, 2.000000E+00, 2.000000E+00, 2.000000E+00
> OBJ3, ARBORIEN, 0.000000E+00, 0.000000E+00, 0.000000E+01
> OBJ3, CLIPART, smoosphe
> OBJ3, ROTATION24, 1
> OBJ3, TYPE, BLOCKAGE
> OBJ3, MATERIAL, 198
STOP