TALK=T;RUN(1,1)
display
This Q1 demonstrates the use of 5 operations used
for the creation of BFC grids, namely:
* defining points,
* defining lines,
* defining frames,
* setting grid dimensions, and
* matching a grid to a frame.
Those operations are used here in order to:
1. create a circle,
2. dimension a grid,
3. match the grid to the circle, and
4. use the Laplace method to improve orthogonality.
The relevant PIL commands will be displayed, then processed.
The first sets the grid type to BFC
The coordinates of the four points used for defining the
circle will be taken as:
A (0,1,1); B(0,2,1); C(0,2,2); D(0,1,2)
The command to define point A is:
The command to define point B is:
The command to define point C is:
The command to define point D is:
The command to draw the arc L1 from A to B which goes through
(0.5*(ya+yb) , is:
The command to draw the arc L2 from B to C
The command to draw the arc L3 from C to D
The command to draw the arc L4 from A to D
The command to define frame F1 which has corner points
A,B,C,D is:
The command to set the BFC grid dimension as 1x20x20 is:
The command to match grid plane I1 on frame F1 using the
trans-finite method is:
Now the PIL command to copy to the I2 plane from the I1
plane by shifting it through a distance (0.5,0,0) is:
By this time, the BFC grid has already set, however,
to make it as orthogonal as possible, a Laplace
solver is used in the following operations-
The PIL command to match grid plane I1 on frame F1 using
the Laplace solver with 20 iterations is:
Finally the PIL command to copy to the I2 plane from the I1
plane by shifting it through a distance (0.5,0,0) is:
enddis
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(GRID_ON_CIRCLE :B544)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 1
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 1
XULAST =1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 20
YVLAST =1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 20
ZWLAST =1.
************************************************************
Group 6. Body-Fitted Coordinates
BFC = T ;NONORT = F
NCRT = -1
RSTGEO = F ;SAVGEO = F
UUP = F ;VUP = F ;WUP = F
NGEOM =CHAM
NAMXYZ =CHAM
ANGMIN =20.
DOMAIN( 1, 1, 1, 21, 1, 21)
* Set fixed sub-domain
FIXDOM( 1, 0, 0, 0, 0, 0, 0)
FIXDOM( 2, 0, 0, 0, 0, 0, 0)
FIXDOM( 3, 0, 0, 0, 0, 0, 0)
FIXDOM( 4, 0, 0, 0, 0, 0, 0)
FIXDOM( 5, 0, 0, 0, 0, 0, 0)
FIXDOM( 6, 0, 0, 0, 0, 0, 0)
FIXDOM( 7, 0, 0, 0, 0, 0, 0)
FIXDOM( 8, 0, 0, 0, 0, 0, 0)
FIXDOM( 9, 0, 0, 0, 0, 0, 0)
FIXDOM( 10, 0, 0, 0, 0, 0, 0)
MSWP = 20
LIJ = F ;LJK = T ;LIK = F
IMON = 1 ;JMON = 11 ;KMON = 11
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.0E-05 ;ENUT =0.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
No PATCHes yet used for this Group
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
No PATCHes yet used for this Group
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP