PHOTON USE p p1 msg Drilling 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 Drilling 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 Drilling 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 Drilling 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 >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>> DISPLAY This input demonstrates the "drilling-milling" technique available in PLANT to generate marker distribution to fit the complex shapes on polar, CARTES=F, grids. The method rests on PLANT ability to handle virtually "live" mathematical expressions as the arguments of its functions. All cases of the file are wholly focused on calculation of marker disributions. Therefore, no other actions are supported by input data. Four different cases, in which the variants of "drilling" are used to get different shapes, are arranged as unsteady problem: at each time moment the new marker distribution is PLANTed calculated and dumped into specified file to be viewed by PHOTON commands supplied. ENDDIS 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 1. Run title and other preliminaries TEXT( DRILling CYLINDER: CARTES=F 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) GROUP 11. Initialization of variable or porosity fields iniadd=f;FIINIT(MARK)=0.0 NAMSAT=MOSG PLANTBEGIN 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 <<<<<<<<<<<<<<<<<<<<< PLANTEND GROUP 15. Termination of sweeps LSWEEP=1 SPEDAT(SET,GXMONI,TRANSIENT,L,F) GROUP 23. Print-out & plot control idispa=1; csg1=p; selref=t; resfac=1.e-2 dmpstk=t DISTIL=T EX(MARK)=1.812E-01 LIBREF=603 STOP