TALK=F;RUN( 1, 1);VDU=VGAMOUSE

TEXT( Library case Y613: Drilling CYLINDER-CARTES=F

  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    PLANT information :
     * Data input groups used: 11, 19
     * Ground groups planted : 11, 19-6
     * Headings used  : INIT??, SC06??
     * Functions used : SPHERE, XYCIRC
     * Commands used  : IF, REGION
  <<<<<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
    GROUP 2. Transience; time-step specification
STEADY=F
GRDPWR(T,4,4.,1.0)
    GROUP 3. X-direction grid specification
CARTES=F;GRDPWR(X,40,2.*3.14,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,40,10.,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,40,20.,1.0)
    GROUP 7. Variables stored, solved & named
STORE(MARK);SOLVE(C1)
    GROUP 11. Initialization of variable or porosity fields
iniadd=f;FIINIT(MARK)=0.0

NAMSAT=MOSG
       1. Drilling the sphere
PATCH(INI1,INIVAL,1,NX,1,NY,1,NZ,1,1)
   VAL=SPHERE(1.,10.,10.,10.,5.)
INIT (INI1,MARK,0.,GRND)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    Ih the above settings the sphere is "drilled-out" in cylinder.
    Crucial faeture of this example is the specification of center
    coordinates: unlike for cartesian cases, the line of X=0 and Y=0 in
    cylindrical polar coordinate system has got the coordinates Xr=RV2D
    and Yr=RV2D in reference frame. It is the latter in which the
    arguments of geometrical functions must be specified.
    <<<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
       2. Drilling the cylinder
   MARK=0.
  IF(ISTEP.EQ.2.AND.LSWEEP.EQ.1)
   MARK=SPHERE(1.,10.,10.,ZGNZ,5.)
  IF(ISTEP.EQ.2.AND.LSWEEP.EQ.1)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    This example seems to be even simpler than foregoing one. It deals
    with the "drilling-out" the cylinder cavity of 5 m diameter in
    coaxial cylindrical domain by spherical "drill bit". The "drilling "
    effect is activated by using the distances from Z=0 plane as
    coordinates of sphere centre.
  <<<<<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
       3. Drilling the bullet shape
   MARK=0.
  IF(ISTEP.EQ.3.AND.LSWEEP.EQ.1)
   MARK=SPHERE(1.,10.,10.,10.,5.)
  IF(ISTEP.EQ.3.AND.LSWEEP.EQ.1)
   MARK=XYCIRC(1.,10.,10.,5.)
  REGION(1,NX,1,NY,5,20,3,3) /LSWEEP.EQ.1
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    Another variant of "drilling" technique is used for specification of
    bullet-shape cavity. The shape is combination of stationary sphere
    and 2D circle "drill-bit". The "drilling" is activated by extents of
    REGION commands over which XYCIRC function is applied. The command
    PLACE may also be used with advantage to employ physical rather than
    cell numbers extents.
  <<<<<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
       4. Drilling the annulus
   MARK=0.
  IF(ISTEP.EQ.4.AND.LSWEEP.EQ.1)
   MARK=SPHERE(1.,15.,10.,ZGNZ,4.)
  IF(ISTEP.EQ.4.AND.LSWEEP.EQ.1)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    This examples differs from "drilling" the cylinder case mainly in
    shifting the axis of "drill bit" of 5 meter away from domain centre
    line.
  <<<<<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
    GROUP 15. Termination of sweeps
LSWEEP=1
    GROUP 23. Print-out & plot control
idispa=1; csg1=p; selref=t; resfac=1.e-2
  PHOTON USE
  p
  p1


  msg  Drill the sphere
  GR Z 1
  con mark x 1 fil;.01
  con mark x 20 fil;.01
  con mark z 20 fil;.01
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit Enter to continue
  pause
  p
  p2


  msg  Drill the cylinder
  gr z 1 y 1 19
  gr z m y 1 19
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit Enter to continue
  pause
  p
  p3


  msg  Drill the bullet
  gr z 5 y 1 19
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit Enter to continue
  pause
  p
  p4


  msg  Drill the annulus
  GR Z 1
  con mark Z 1 X 1 20 Y 1 M fil;.01
  con mark Z 10 X 1 20 Y 1 M fil;.01
  con mark Z m X 1 20 Y 1 M fil;.01
  GR OU Y M;GR OU Z M;GR OU Z M
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit E to finish
  enduse
STOP