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. create a square, and
3. to match them to the two sides of a constant-I plane,
(namely I2 and I1 respectively).
The relevant PIL commands will be displayed, then processed.
The first sets the grid type to BFC
Please enter the circle radius (default 0.25m)
Please enter the side-length of the square (default 0.5m)
The input data to be used in this example are..
The centre of the circle will be placed at:
x = 0.0
y = - square root of 2 times the radius
z = + square root of 2 times the radius
Now define points A,B,C and D, positioned as follows:
A at (0,-root,root)
B at (0,root,root)
C at (0,root,-root)
D at (0,-root,-root)
by the commands:
for A
for B
for C
and for D
Now draw four arcs (AB,BC,CD,DA) from A to D
to create a circle.
Arcs AB and CD must have equal numbers of intervals;
so must arcs BC and DA.
The intervals will be supposed uniform on each arc.
Please enter nab, the number of cells for AB and CD (default 5)
The command to create arc AB, with nab uniform intervals,
which also goes through the point (0,0,radius),
is:
Please enter nbc, the number of cells for BC and DA (default 5)
The command to create arc BC, with nbc cells
which also goes through the point (0,0,radius),
is:
The corresponding command for CD is:
and for DA:
Now define a frame ABCD which has corner points
A,B,C & D by the command:
Now set the BFC grid dimension as 1 x nab x nbc, with
reference length 1m x side x side, by the command:
Now match grid plane I2 on frame ABCD using the
trans-finite method, by the command:
Now create the square on plane I = 1
Define points E,F,G and H, with positions:
E at (-side,-side/2,side/2)
F at (-side,side/2,side/2)
G at (-side,side/2,-side/2)
H at (-side,-side/2,-side/2)
by the commands:
for E:
for F:
for G:
and for H:
Now draw four lines (EF, FG, GH, HE) from E to H
to create a square.
Lines EF & GH each have nab uniformly-distributed intervals,
and lines FG & HE each have nbc uniformly-distributed
intervals, so as to match the numbers of intervals on the
circle.
The commands are:
then
then
and finally
Now define a frame EFGH (square) which has corner points
E, F, G and H, by the command:
Now match grid plane I1 on frame EFGH using the trans-finite
by the command:
This completes the exercise.
However, you may now wish to type VIEW, so as to be able to
inspect, or indeed alter, the grid by means of the graphical
grid generator.
enddis
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(GSET COMMAND TUTORIAL 1 :B538)
************************************************************
************************************************************
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 = 5
YVLAST =0.5
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 5
ZWLAST =0.5
************************************************************
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, 6, 1, 6)
* 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)
LIJ = F ;LJK = T ;LIK = F
IMON = 1 ;JMON = 3 ;KMON = 3
************************************************************
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