PHOTON USE
p
msg pressure contours and velocity vectors
msg (for case 276, pressure represents velocity potential)
con p1 x 1 fi;0.05; con p1 y 1 fi;0.05
set ref vec 40.0; vec x 1; vec y 1
msg
msg Press e to END
ENDUSE
****** TO LOAD CASE: TYPE L(274) ******
GROUP 1. Run title and other preliminaries
TEXT(3D Nav.Sto.Flow Around A Moving Van
TITLE
DISPLAY
Note: January 2009, In order that this could be used with
the grid-refinement module, SPINTO, the earler use of CONPOR
with its indicial arguments has been replaced by use of initial-
value patches such as:
** Vehicle body
patch(BODY1,inival,object,body1, 1,1)
coval(BODY1,prps,0.0,198)
wherein the arguments 'object' and 'body1' place the initial
values in the space indicated by the geometric arguments of
RSET(B,BODY1,,,,,)
This is a 3-dimensional Cartesian example in which a solid
body immersed in an unconfined flow is represented by means of
'porosities'. The turbulence of the flow is crudely represented
by a prescribed constant turbulent viscosity within the fluid,
together with a logarithmic law of the wall for skin friction.
The lower (south) boundary of the flow domain, representing
the road surface, is modelled as a moving wall.
An impression of the geometry can be gained by examining
the following sketch. |
______________ |______
|**************| |******|
^ __|**************| ^ |******|
| |*****************| | |******|
y| |*****************| y| |******|
| ~**~~~~~~~~~~~~~**~ | |~~~~**
** ** | **
-------------------------------------------
z----> x--->
ENDDIS
GROUP 2. Transience; time-step specification
GROUP 3. X-direction grid specification
nx=8; ny=12; nz=20
incl(spinxyz1) ! move this line to the left when using SPINTO
* Write objects
RSET(D,VAN,5.0,6.0,20.,1,0)
RSET(M,nx,ny,nz)
RSET(B,FWHL, 1., 0, 5., 0.5, 0.5, 1., 12, 0)
RSET(B,RWHL, 1., 0, 9., 0.5, 0.5, 1., 12, 0)
RSET(B,BODY1, 0, 0.5, 5., 1.5, 1., 5., 12, 0)
RSET(B,BODY2, 0, 1.5, 6., 1.5, 1.5, 4., 12, 0)
RSET(B,ROAD, 0, 0, 0 , 5., 0, 20., 10, 0)
RSET(B,UPSTR, 0, 0, 0, 5., 6., 0, 7, 0)
RSET(B,DWSTR, 0, 0, 20., 5., 6., 0, 5, 0)
GROUP 6. Body-fitted coordinates or grid distortion
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1)
**Select whole-field solution procedure for pressure,
and point-by-point for velocities
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)
STORE(PRPS,IMB1)
GROUP 8. Terms (in differential equations) & devices
GROUP 9. Properties of the medium (or media)
ENUL=1.E-5; ENUT=1.0E-3
GROUP 10. Inter-phase-transfer processes and properties
GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=0.0; FIINIT(W1)=14.0; fiinit(prps)=-1
** Vehicle body
PATCH(BODY1,INIVAL,OBJECT,BODY1, 1,1)
COVAL(BODY1,PRPS,0.0,198)
PATCH(BODY2,INIVAL,OBJECT,BODY2, 1,1)
COVAL(BODY2,PRPS,0.0,198)
** Front wheel
PATCH(FWHL,INIVAL,OBJECT,FWHL, 1,1)
COVAL(FWHL,PRPS,0.0,198)
** Rear wheel
PATCH(RWHL,INIVAL,OBJECT,RWHL, 1,1)
COVAL(RWHL,PRPS,0.0,198)
GROUP 13. Boundary conditions and special sources
** Upstream boundary
PATCH(UPSTR,LOW,OBJECT,UPSTR,1,1)
COVAL(UPSTR,P1,FIXFLU,14.0); COVAL(UPSTR,W1,ONLYMS,14.)
** Downstream boundary
PATCH(DWSTR,HIGH,object,dwstr,1,1)
COVAL(DWSTR,P1,0.1,0.)
** Road surface
PATCH(ROAD,SWALL,OBJECT,ROAD,1,1)
COVAL(ROAD,W1,LOGLAW,14.0)
GROUP 15. Termination of sweeps
LSWEEP=200
GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,1.0E-1);RELAX(V1,FALSDT,1.0E-1)
RELAX(W1,FALSDT,1.0E-1)
VARMIN(U1)=-1.E11;VARMAX(U1)=1.0
VARMIN(V1)=-1.E11;VARMAX(V1)=1.0
VARMIN(W1)=-1.E11;VARMAX(W1)=1.0
GROUP 21. Print-out of variables
YZPR=T
GROUP 22. Spot-value print-out
IXMON=5;IYMON=5;IZMON=5
GROUP 22. Spot-value print-out
TSTSWP=-1
GROUP 23. Field print-out and plot control
NXPRIN=NX/2;NYPRIN=NY/5;NZPRIN=NZ/10
PATCH(PROFILE,PROFIL,NX-1,NX-1,NY/2,NY/2,1,NZ,1,1)
PLOT(PROFILE,W1,0.0,0.0);PLOT(PROFILE,P1,0.0,0.0)
PLOT(PROFILE,U1,0.0,0.0);PLOT(PROFILE,V1,0.0,0.0)
PATCH(XEQ1,CONTUR,1,1,1,NY,1,NZ,1,1);PLOT(XEQ1,P1,0.0,20.0)
PLOT(XEQ1,U1,0.0,20.0);PLOT(XEQ1,V1,0.0,20.0)
PLOT(XEQ1,W1,0.0,20.0)
PATCH(XEQ5,CONTUR,5,5,1,NY,1,NZ,1,1);PLOT(XEQ5,P1,0.0,20.0)
PLOT(XEQ5,W1,0.0,20.0)
PATCH(WAKE,CONTUR,1,NX,1,NY,NZ-5,NZ-5,1,1)
PLOT(WAKE,W1,0.0,20.0)
#endpause
#maxmin