TALK=T;RUN(1,1)
DISPLAY
This case exemplifies some of the basic techniques
available in In-Form for field initialisations.
Three examples are provided, each for a different xy-plane
in a single 3D box.
Example 1: Initialization of stagnation-point flow on IZ=1
===============================================
PATCH(INI1,INIVAL,1,NX,1,NY,1,1,1,1)
(INITIAL of U1 at INI1 is XU-10.)
(INITIAL of V1 at INI1 is -(YV-10.))
Example 2: Solid-body rotation flow on IZ=2
================================
PATCH(INI2,INIVAL,1,NX,1,NY,2,2,1,1)
(INITIAL of U1 at INI2 is YG-10.)
(INITIAL of V1 at INI2 is -(XG-10.))
Example 3: The values on IZ=1 and IZ=2 are added together
and placed on IZ=3
==============================================
PATCH(INI3,INIVAL,1,NX,1,NY,3,3,1,1)
(INITIAL of U1 at INI3 is -U1[&&-1]+U1[&&1])
(INITIAL of V1 at INI3 is -V1[&&2]+V1[&&-2])
The Q1 contains PHOTON USE commands
ENDDIS
PHOTON USE
p
msg
msg
msg All examples which follow are assembled as different slabs
msg and number of slabs of one 3D domain called
msg
msg INITBOX.
msg
msg The series of PHOTON plots will follow.
msg
msg They illustrate some of the existing initialisation
msg techniques available via In-Form.
msg
msg The presentation covers:
msg
msg * Flow field initialistions by formulae
msg * Manipulating with initial fields
msg * MARKing sub-domains to create arbitrary initial fields
msg * Geometry initialisations
msg
msg Hit enter to proceed to 1st slab.
msg
pause
clear
vi z
gr ou z 1
set vec ref 25
vec z 1 sh
msg
msg
msg The 1st slab of the box is initialised by stagnation
msg point flow
msg Hit enter to proceed to 2nd slab.
pause
vec off;red
vec z 2 sh
msg
msg
msg The 2nd slab is initialised by solid body rotation flow
msg
msg Hit enter to proceed to 3rd slab.
pause
vec off;red
vec z 3 sh
msg
msg The flow in the 3rd slab is constructed by adding the
msg velocity distributions of 1st and 2nd slabs.
msg
msg Hit E to finish.
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Three examples of initialization :758 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 758
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 40
XULAST =20.
XFRAC(1)=0.025 ;XFRAC(2)=0.05
XFRAC(3)=0.075 ;XFRAC(4)=0.1
XFRAC(5)=0.125 ;XFRAC(6)=0.15
XFRAC(7)=0.175 ;XFRAC(8)=0.2
XFRAC(9)=0.225 ;XFRAC(10)=0.25
XFRAC(11)=0.275 ;XFRAC(12)=0.3
XFRAC(13)=0.325 ;XFRAC(14)=0.35
XFRAC(15)=0.375 ;XFRAC(16)=0.4
XFRAC(17)=0.425 ;XFRAC(18)=0.45
XFRAC(19)=0.475 ;XFRAC(20)=0.5
XFRAC(21)=0.525 ;XFRAC(22)=0.55
XFRAC(23)=0.575 ;XFRAC(24)=0.6
XFRAC(25)=0.625 ;XFRAC(26)=0.65
XFRAC(27)=0.675 ;XFRAC(28)=0.7
XFRAC(29)=0.725 ;XFRAC(30)=0.75
XFRAC(31)=0.775 ;XFRAC(32)=0.8
XFRAC(33)=0.825 ;XFRAC(34)=0.85
XFRAC(35)=0.875 ;XFRAC(36)=0.9
XFRAC(37)=0.925 ;XFRAC(38)=0.95
XFRAC(39)=0.975 ;XFRAC(40)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 40
YVLAST =20.
YFRAC(1)=0.025 ;YFRAC(2)=0.05
YFRAC(3)=0.075 ;YFRAC(4)=0.1
YFRAC(5)=0.125 ;YFRAC(6)=0.15
YFRAC(7)=0.175 ;YFRAC(8)=0.2
YFRAC(9)=0.225 ;YFRAC(10)=0.25
YFRAC(11)=0.275 ;YFRAC(12)=0.3
YFRAC(13)=0.325 ;YFRAC(14)=0.35
YFRAC(15)=0.375 ;YFRAC(16)=0.4
YFRAC(17)=0.425 ;YFRAC(18)=0.45
YFRAC(19)=0.475 ;YFRAC(20)=0.5
YFRAC(21)=0.525 ;YFRAC(22)=0.55
YFRAC(23)=0.575 ;YFRAC(24)=0.6
YFRAC(25)=0.625 ;YFRAC(26)=0.65
YFRAC(27)=0.675 ;YFRAC(28)=0.7
YFRAC(29)=0.725 ;YFRAC(30)=0.75
YFRAC(31)=0.775 ;YFRAC(32)=0.8
YFRAC(33)=0.825 ;YFRAC(34)=0.85
YFRAC(35)=0.875 ;YFRAC(36)=0.9
YFRAC(37)=0.925 ;YFRAC(38)=0.95
YFRAC(39)=0.975 ;YFRAC(40)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 3
ZWLAST =0.5
ZFRAC(1)=0.333333 ;ZFRAC(2)=0.666667
ZFRAC(3)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(3)=U1 ;NAME(5)=V1
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(U1,Y,N,N,N,N,Y)
SOLUTN(V1,Y,N,N,N,N,Y)
************************************************************
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)
FIINIT(U1)=1.0E-10 ;FIINIT(V1)=1.0E-10
PATCH(INI1 ,INIVAL, 1, 40, 1, 40, 1, 1, 1, 1)
PATCH(INI2 ,INIVAL, 1, 40, 1, 40, 2, 2, 1, 1)
PATCH(INI3 ,INIVAL, 1, 40, 1, 40, 3, 3, 1, 1)
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
RELAX(U1,LINRLX,1.)
RELAX(V1,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,INITIAL,U1!INI1,C,=XU-10.)
SPEDAT(SET,INITIAL,V1!INI1,C,=-(YV-10.))
SPEDAT(SET,INITIAL,U1!INI2,C,=YG-10.)
SPEDAT(SET,INITIAL,V1!INI2,C,=-(XG-10.))
SPEDAT(SET,INITIAL,U1!INI3,C,=-U1[&&-1]+U1[&&1])
SPEDAT(SET,INITIAL,V1!INI3,C,=-V1[&&2]+V1[&&-2])
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(U1)=5.485 ;EX(V1)=5.485
************************************************************
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
OUTPUT(U1,Y,N,Y,N,N,N)
OUTPUT(V1,Y,N,Y,N,N,N)
************************************************************
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
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
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