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
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. The example is
an up-to-date version of that reported in the report PDR/CHAM UK/5
An impression of the geometry can be gained by examining
the following sketch. |
______________ |______
|**************| |******|
^ __|**************| ^ |******|
| |*****************| | |******|
y| |*****************| y| |******|
| ~**~~~~~~~~~~~~~**~ | |~~~~**
** ** | **
-------------------------------------------
z----> x--->
ENDDIS
#pause
GROUP 2. Transience; time-step specification
GROUP 3. X-direction grid specification
nx=8; ny=12; nz=20
incl(spinxyz1)
* 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)
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.)
COVAL(DWSTR,U1,ONLYMS,0.0);COVAL(DWSTR,V1,ONLYMS,0.0)
COVAL(DWSTR,W1,ONLYMS,0.0)
** Road surface
patch(ROAD,swall,object,road,1,1)
COVAL(ROAD,W1,LOGLAW,14.0)
** Following Patches are only necessary if BODY definition
above is set with positive limits **
GROUP 15. Termination of sweeps
LSWEEP=200
LSWEEP=20
GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,1.0e-4);RELAX(V1,FALSDT,1.0e-4)
RELAX(W1,FALSDT,1.0e-4)
SARAH=0.005
varmin(u1)=-1.e11;varmax(u1)=1.0
varmin(v1)=-1.e11;varmax(v1)=1.0
varmin(w1)=-1.e11;varmax(w1)=2.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
#conprom
lsweep=10
lsweep=50
lsweep=5
incl(spinxyz3)