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. Then, uniform u1 is imposed
in the direction of the temperature gradient. The soulution
is one dimensional with exponential profile
(exp((u1*rho1*cp1)/cond)*x) - 1
T(x) =(------------------------------------------ )*(T2-T1)+T1
( exp((u1*rho1*cp1)/cond)*L) -1
at x=0 T=T1 ; at x=L T=T2
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 322
************************************************************
Group 1. Run Title
TEXT(Simple f304 (tem1 & uniform u1, 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,1.000000E+02)
************************************************************
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,W1 ,V1 ,U1 )
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 ANT1 is ALCO(TEM1))
(STORED of AST1 is AHCO(TEM1))
** Exchange coefficients
of TEM1, temperature
(STORED of GAt1 is GAMM(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(U1 ) = 5.000000E-01 ;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
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL=F
EX(U1 )= 5.000E-01 ; EX(V1 )= 1.000E-10
EX(W1 )= 1.000E-10 ; EX(GAT1)= 1.000E+00
EX(AST1)= 1.250E+18 ; EX(ANT1)= 8.750E-01
EX(AWT1)= 2.734E-01 ; EX(AET1)= 1.641E-01
EX(APT1)= 1.250E-01 ; EX(TEM1)= 1.762E+03
EX(TCOR)= 1.000E-10 ; EX(TRES)= 1.000E-10
EX(AHT1)= 1.000E-10 ; EX(ALT1)= 1.000E-10
EX(PRPS)= 1.000E+00
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 2 ;IYMON = 8 ;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,-1.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)
> DOM, SIZE, 8.000000E+00, 4.000000E+00, 4.000000E+00
> DOM, MONIT, 1.500000E+00, 3.750000E+00, 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