TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.5, Date 16/12/02
CPVNAM=VDI;SPPNAM=Core
************************************************************
PHOTON USE
p; ; ; ; ;
msg Contours of temperature Tem1
do jj = 1, 6
con tem1 j jj fi;0.005
enddo
ENDUSE
Echo DISPLAY / USE settings
DISPLAY
This case uses the geometry and finer grid volumes of F304.
Temperatures at opposite boundaries along x-direction is
fixed as boundary conditions. The solution is one dimensional
with linear profile
T(x) = (T2-T1)(x/L)+T1
at x=0 T=T1 ; at x=L T=T2
In-Form facility is used for printing coefficients, residuals
and corrections and extinction coeffieint.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 321
************************************************************
Group 1. Run Title
TEXT(Simple f304 (tem1, 3D) )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,8,8,8)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(148) =TEM1 ; NAME(149) =PRPS
NAME(150) =ASAP
* Solved variables list
SOLVE(TEM1)
* Stored variables list
STORE(ASAP,PRPS)
STORE(APT1,AET1,AWT1,ANT1,AST1)
STORE(ALT1,AHT1,GAT1,TRES,TCOR)
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Echo InForm settings for Group 7
inform7begin
** temperature equation coefficients
(STORED of APT1 is APCO(TEM1))
(STORED of AET1 is AECO(TEM1))
(STORED of AWT1 is AWCO(TEM1))
(STORED of ANT1 is ANCO(TEM1))
(STORED of AST1 is ASCO(TEM1))
(STORED of ALT1 is ALCO(TEM1))
(STORED of AHT1 is AHCO(TEM1))
** Exchange coefficients of TEM1, temperature
(STORED of GAt1 is GAMM(TEM1))
** Print-out residuals of temperature
(STORED of TRES is RESI(TEM1))
** Print-out of temperature corrections
(STORED of TCOR is CORR(TEM1))
inform7end
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05
RHO1 = 1.000000E+00
ENUL = 1.000000E+00
CP1 = 1.000000E+00
ENUT = 0.000000E+00
PRNDTL(TEM1) = 1.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(TEM1) = 1.000000E+03 ;FIINIT(PRPS) = -1.000000E+00
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
No PATCHes used for this Group
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 30
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL=F
EX(GAT1)= 1.000E+00 ; EX(AHT1)= 1.250E+18
EX(ALT1)= 8.750E-01 ; EX(AST1)= 8.750E-01
EX(ANT1)= 8.750E-01 ; EX(AWT1)= 2.188E-01
EX(AET1)= 2.188E-01 ; EX(APT1)= 1.250E-01
EX(TEM1)= 1.500E+03 ; EX(PRPS)= 1.000E+00
EX(TCOR)= 2.862E-05 ;EX(TRES)= 2.094E-01
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 2 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
NYPRIN = 1
NZPRIN = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.000000E+00,-9.887711E-01,1.494381E-01)
GVIEW(UP,0.000000E+00,1.494381E-01,9.887711E-01)
> DOM, SIZE, 8.000000E+00, 4.000000E+00, 4.000000E+00
> DOM, MONIT, 1.500000E+00, 2.500000E-01, 2.500000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, B1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 4.000000E+00, 4.000000E+00
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, NAME, B2
> OBJ, POSITION, 8.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 4.000000E+00, 4.000000E+00
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, -1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, NAME, B3
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 8.000000E+00, 0.000000E+00, 4.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, -1
> OBJ, TYPE, PLATE
> OBJ, NAME, B4
> OBJ, POSITION, 0.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, SIZE, 8.000000E+00, 0.000000E+00, 4.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, -1
> OBJ, TYPE, PLATE
> OBJ, NAME, B5
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 8.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, B6
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 4.000000E+00
> OBJ, SIZE, 8.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, B7
> OBJ, POSITION, 2.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, SIZE, 4.000000E+00, 2.000000E+00, 2.000000E+00
> OBJ, CLIPART, fine
> OBJ, ROTATION24, 1
> OBJ, TYPE, FINE_GRID_VOL
> OBJ, FINEFACTS, 2.000000E+00, 2.000000E+00, 2.000000E+00
> OBJ, NAME, PLT2
> OBJ, POSITION, 8.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 4.000000E+00, 4.000000E+00
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00, 1.000000E+03
> OBJ, NAME, PLT1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 4.000000E+00, 4.000000E+00
> OBJ, CLIPART, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_TEMP, 0.000000E+00, 2.000000E+03
SINGLE SWEEP SOLUTION
integer(icase)
icase=2
if(icase.eq.2) then
LSWEEP = 1
RELAX(TEM1,LINRLX, 1.0000)
LITER (TEM1) = 1
spedat(set,fgmb_itr,tem1,i,10000)
endif
STOP