DISPLAY
 
 
  SOLID-PROPELLANT ROCKET
 
  3-dimensional, Polar, steady, elliptic simulation
 
 
  Demonstration of flow in an idealised 3D solid-propellant rocket
  with a (roughly) star-shaped charge.
 
 
    ^y         __--__                   --------------------  ^y
    |         /   |  \--solid              solid              |
    +--->x    \   \  /  propellant         propellant         +--->z
               \   \/                   --------------------
                 \/                     --------------------
 
 
   Time required for complete run appx. 1.5 min on PC/486-33
 
 
  enddis
 
  PHOTON USE
  p
 
 
 
  view z
  gr ou z 1
  msg Gas-velocity vectors
  vec z 5 sh
  msg            -
  msg Press  to continue
  pause
  vec off;red
  msg W1 contours at iz = nz/2
  gr ou x 1;con w1 z 5 fi;0.001
  msg            -
  msg Press  to continue
  pause
  con off
  view x
  msg pressure contours at ix=1
  con p1 x 1 fi;0.001
  msg            -
  msg Press  to continue
  pause
  con off;red
  msg velocity vectors at ix = 1
  vec x 1 sh
  msg            -
  msg Press  to continue
  pause
  vec off;gr off;red
  gr ou x m;view 1 -10 0
  vec x m sh
  msg velocity vectors at ix=nx
  msg            -
  msg Press e to END
  enduse
 
real(burnrate)
    GROUP 1. Run title and other preliminaries
TEXT(3D Solid-Propellant Rocket        
TITLE
libref=895
burnrate=0.1
mesg(Surface burning rate is specified as :burnrate:
mesg(                                      kg/m**2 s
mesg(If another rate is desired, enter it now please.
mesg(Otherwise, press RETURN
readvdu(burnrate,real,burnrate)
    GROUP 2. Transience; time-step specification
    GROUP 3. X-direction grid specification
  ** Activate polar coordinate system
CARTES=F
grdpwr(x,10,1.0,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,10,0.01,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,10,0.10,1.0)
    GROUP 6. Body-fitted coordinates or grid distortion
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1)
SOLUTN(P1,P,P,Y,P,P,P)
    GROUP 8. Terms (in differential equations) & devices
    GROUP 9. Properties of the medium (or media)
ENUL=1.E-5
    GROUP 10. Inter-phase-transfer processes and properties
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.1
conpor(prop,0.0,volume,-(nx/2+1),-nx,-(ny/2+1),-ny,-1,-nz)
    GROUP 13. Boundary conditions and special sources
WALL (WALL,NORTH,1,2,NY,NY,1,NZ,1,1)
COVAL(WALL,W1,1.0,0.0)
coval(wall,u1,1.0,0.0)
INLET(IN,north,1,nx/2,ny,NY,1,nz,1,1)
VALUE(IN,P1,burnrate);VALUE(IN,W1,0.0)
VALUE(IN,V1,0.0)
patch(in2,west,nx/2,nx/2,ny/2+1,ny,1,nz,1,1)
coval(in2,p1,fixflu,burnrate)
patch(in3,north,nx/2+1,nx,ny/2,ny/2,1,nz,1,1)
coval(in3,p1,fixflu,burnrate)
PATCH(OUTLET,HIGH,1,nx,1,NY/2,NZ,NZ,1,1)
COVAL(OUTLET,P1,FIXVAL,0.0)
COVAL(OUTLET,V1,ONLYMS,0.0);COVAL(OUTLET,W1,ONLYMS,0.0)
    GROUP 14. Downstream pressure for PARAB=.TRUE.
    GROUP 15. Termination of sweeps
LSWEEP=300
    GROUP 16. Termination of iterations
    GROUP 17. Under-relaxation devices
RELAX(W1,FALSDT,0.1);RELAX(V1,FALSDT,0.1)
RELAX(U1,FALSDT,0.1)
    GROUP 18. Limits on variables or increments to them
    GROUP 19. Data communicated by satellite to GROUND
    GROUP 20. Preliminary print-out
    GROUP 21. Print-out of variables
    GROUP 22. Spot-value print-out
NPRMON=20;IYMON=5
    GROUP 22. Spot-value print-out
IYMON=ny/2;IZMON=nz/2;IPLTL=LSWEEP
    GROUP 23. Field print-out and plot control
NYPRIN=2;NZPRIN=5;NPLT=1
PATCH(EXIT,PROFIL,1,1,1,NY,NZ-1,NZ-1,1,1)
PLOT(EXIT,W1,0.0,0.0)
PATCH(PIPE,CONTUR,1,1,1,NY,1,NZ,1,1)
PLOT(PIPE,W1,0.0,10.0);PLOT(PIPE,H1,0.0,10.0)
 
selref=t;resfac=1.e-2
TSTSWP=-1